Remarkable expansions of an X-linked reproductive homeobox gene cluster in rodent evolution

Rhox is a recently identified cluster of 12 X-linked homeobox genes in mice. The expression pattern of Rhox genes during postnatal testis development corresponds to their chromosomal position, much like the colinear gene regulation of the Hox gene clusters during animal embryonic development. We here report the identification of 18 additional Rhox genes and 3 pseudogenes in mice. Comparative analyses of the mouse, rat, human, dog, cow, opossum, and chicken genomes suggest that the Rhox cluster originated in the common ancestor of primates and rodents. It subsequently underwent two remarkable expansions, first in the common ancestor of mice and rats and then in mice. Positive selection promoting amino acid substitutions was detected in some young Rhox genes, suggesting adaptive functional diversification. The recent expansions of the Rhox cluster provide an opportunity to study the mechanism and origin of colinear gene regulation, but they may also undermine the utility of mouse models for understanding the development and physiology of the human reproductive system.     

Overestimation of nonsynonymous/synonymous rate ratio by reverse-translation of aligned amino acid sequences

In the analysis of protein-coding nucleotide sequences, the ratio of the number of nonsynonymous substitutions to that of synonymous substitutions (d(N)/d(S)) is used as an indicator for the direction and magnitude of natural selection operating at the amino acid sequence level. The d(S) and d(N) values are estimated based on the comparison of homologous codons, which are often identified by converting (reverse-translating) aligned amino acid sequences into codon sequences. In this method, however, homologous codons may be mis-identified when frame-shifts occurred or amino acid sequences were mis-aligned, which may lead to overestimation of the d(N)/d(S) ratio. Here the effect of reverse-translating aligned amino acid sequences on the estimation of d(N)/d(S) ratio was examined through a large-scale analysis of protein-coding nucleotide sequences from vertebrate species. Apparently, 1-9% of codon sites that were identified as homologous with reverse-translation contained non-homologous codons, where the d(N)/d(S) ratio was unduly high. By correcting the d(N)/d(S) ratio for these codon sites, it was inferred that the ratio was 5-43% overestimated with reverse-translation. These results suggest that caution should be exerted in the study of natural selection using the d(N)/d(S) ratio by reverse-translating aligned amino acid sequences.     

The Rates of Molecular Evolution in Rodent and Primate Mitochondrial DNA

A higher rate of molecular evolution in rodents than in primates at synonymous sites and, to a lesser extent, at amino acid replacement sites has been reported previously for most nuclear genes examined. Thus in these genes the average ratio of amino acid replacement to synonymous substitution rates in rodents is lower than in primates, an observation at odds with the neutral model of molecular evolution. Under Ohta's mildly deleterious model of molecular evolution, these observations are seen as the consequence of the combined effects of a shorter generation time (driving a higher mutation rate) and a larger effective population size (resulting in more effective selection against mildly deleterious mutations) in rodents. The present study reports the results of a maximum-likelihood analysis of the ratio of amino acid replacements to synonymous substitutions for genes encoded in mitochondrial DNA (mtDNA) in these two lineages. A similar pattern is observed: in rodents this ratio is significantly lower than in primates, again consistent only with the mildly deleterious model. Interestingly the lineage-specific difference is much more pronounced in mtDNA-encoded than in nuclear-encoded proteins, an observation which is shown to run counter to expectation under Ohta's model. Finally, accepting certain fossil divergence dates, the lineage-specific difference in amino acid replacement-to-synonymous substitution ratio in mtDNA can be partitioned and is found to be entirely the consequence of a higher mutation rate in rodents. This conclusion is consistent with a replication-dependent model of mutation in mtDNA. Received: 24 September 1999 / Accepted: 18 September 2000     

Adaptive evolution of the IgA hinge region in primates

IgA is a major component that prevents the penetration of pathogenic bacteria into mucosal surfaces. The IgA antibody is cleaved at the IgA hinge region with high specificity by IgA-specific proteases produced by several pathogenic bacteria. We conducted a genomic sequence analysis of the IgA genes of a wide spectrum of primates, including the first intron and second exon, which consist of the hinge region and the CH2 domain, to find evidence of positive selection. Because the hinge region is quite small, we combined the largest collection of sequences that could be clearly aligned and evaluated the total numbers of synonymous and nonsynonymous substitutions on the phylogenetic tree. The nonsynonymous to synonymous substitution ratio (d(N)/d(S) test) showed that hominoids, Old World monkeys, and New World monkeys have d(N)/d(S) ratios of 5.4, 6.3, and 4.2, respectively. Fisher's exact probability tests showed statistical significance for the Old World monkey. Because the substitution rates of the flanking sequences are more or less similar to the synonymous rates of the hinge region, these high values of d(N)/d(S) should be the result of positive selection at the hinge region. Combining the high sequence variability in each population and the highly accelerated nonsynonymous substitution rates in the hinge region, we conclude that this unusual IgA evolution is a molecular evidence of adaptive evolution possibly caused by the host-parasite relationship.     

The effect of tandem substitutions on the correlation between synonymous and nonsynonymous rates in rodents.

Nonsynonymous substitutions in DNA cause amino acid substitutions while synonymous substitutions in DNA leave amino acids unchanged. The cause of the correlation between the substitution rates at nonsynonymous (K(A)) and synonymous (K(S)) sites in mammals is a contentious issue, and one that impacts on many aspects of molecular evolution. Here we use a large set of orthologous mammalian genes to investigate the causes of the K(A)-K(S) correlation in rodents. The strength of the K(A)-K(S) correlation exceeds the neutral theory expectation when substitution rates are estimated using algorithmic methods, but not when substitution rates are estimated by maximum likelihood. Irrespective of this methodological uncertainty the strength of the K(A)-K(S) correlation appears mostly due to tandem substitutions, an excess of which is generated by substitutional nonindependence. Doublet mutations cannot explain the excess of tandem synonymous-nonsynonymous substitutions, and substitution patterns indicate that selection on silent sites is the likely cause. We find no evidence for selection on codon usage. The nature of the relationship between synonymous divergence and base composition is unclear because we find a significant correlation if we use maximum-likelihood methods but not if we use algorithmic methods. Finally, we find that K(S) is reduced at the start of genes, which suggests that selection for RNA structure may affect silent sites in mammalian protein-coding genes.     

Parallel acceleration of evolutionary rates in symbiont genes underlying host nutrition

The overproduction of essential amino acids by Buchnera aphidicola, the primary bacterial mutualist of aphids, is considered an adaptation for increased production of nutrients that are lacking in aphids' diet of plant sap. Given their shared role in host nutrition, amino acid biosynthetic genes of Buchnera are expected to experience parallel changes in selection that depend on host diet quality, growth rate, and population structure. This study evaluates the hypothesis of parallel selection across biosynthetic pathways by testing for correlated changes in evolutionary rates at biosynthetic genes of Buchnera. Previous studies show fast evolutionary rates at tryptophan biosynthetic genes among Buchnera associated with the aphid genus Uroleucon and suggest reduced purifying selection on symbiont nutritional functions in this aphid group. Here, we test for parallel rate acceleration at other amino acid biosynthetic genes of Buchnera-Uroleucon, including those for leucine (leuABC) and isoleucine/valine biosynthesis (ilvC). Ratios of nonsynonymous to synonymous substitutions (d(N)/d(S)) were estimated using codon-based maximum-likelihood methods that account for the extreme AT compositional bias of Buchnera sequences. A significant elevation in d(N)/d(S) at biosynthetic loci but not at two housekeeping genes sampled (dnaN and tuf) suggests reduced host-level selection on biosynthetic capabilities of Buchnera-Uroleucon. In addition, the discovery of trpEG pseudogenes in Buchnera-U. obscurum further supports reduced selection on amino acid biosynthesis.     

Sexual selection and the adaptive evolution of mammalian ejaculate proteins

An elevated rate of substitution characterizes the molecular evolution of reproductive proteins from a wide range of taxa. Although the selective pressures explaining this rapid evolution are yet to be resolved, recent evidence implicates sexual selection as a potentially important explanatory factor. To investigate this hypothesis, we sought evidence of a high rate of adaptive gene evolution linked to postcopulatory sexual selection in muroid rodents, a model vertebrate group displaying a broad range of mating systems. Specifically, we sequenced 7 genes from diverse rodents that are expressed in the testes, prostate, or seminal vesicles, products of which have the potential to act in sperm competition. We inferred positive Darwinian selection in these genes by estimation of the ratio of nonsynonymous (d(N), amino acid changing) to synonymous (d(S), amino acid retaining) substitution rates (omega = d(N)/d(S)). Next, we tested whether variation in this ratio among lineages could be attributed to interspecific variation in mating systems, as inferred from the variation in these rodents' relative testis sizes (RTS). Four of the 7 genes examined (Prm1, Sva, Acrv1, and Svs2, but not Svp2, Msmb, or Spink3) exhibit unambiguous evidence of positive selection. One of these, the seminal vesicle-derived protein Svs2, also shows some evidence for a concentration of positive selection in those lineages in which sperm competition is common. However, this was not a general trend among all the rodent genes we examined. Using the same methods, we then reanalyzed previously published data on 2 primate genes, SEMG1 and SEMG2. Although SEMG2 also shows evidence of positive selection concentrated in lineages subject to high levels of sperm competition, no such trend was found for SEMG1. Overall, despite a high rate of positive selection being a feature of many ejaculate proteins, these results indicate that the action of sexual selection potentially responsible for elevated evolutionary rates may be difficult to detect on a gene-by-gene basis. Although the extreme diversity of reproductive phenotypes exhibited in nature attests to the power of sexual selection, the extent to which this force predominates in driving the rapid molecular evolution of reproductive genes therefore remains to be determined.     

Natural selection on the influenza virus genome

Influenza viruses are the etiological agents of influenza. Although vaccines and drugs are available for the prophylaxis and treatment of influenza virus infections, the generation of escape mutants has been reported. To develop vaccines and drugs that are less susceptible to the generation of escape mutants, it is important to understand the evolutionary mechanisms of the viruses. Here natural selection operating on all the proteins encoded by the H3N2 human influenza A virus genome was inferred by comparing the numbers of synonymous (d(S) [D(S)]) and nonsynonymous (d(N) [D(N)]) substitutions per site. Natural selection was also inferred for the groups of functional amino acid sites involved in B-cell epitopes (BCEs), T-cell epitopes (TCEs), drug resistance, and growth in eggs. The entire region of PB1-F2 was positively selected, and positive selection also appeared to operate on BCEs, TCEs, and growth in eggs. The frequency of escape mutant generation appeared to be positively correlated with the d(N)/d(S) (D(N)/D(S)) values for the targets of vaccines and drugs, suggesting that the amino acid sites under strong functional constraint are suitable targets. In particular, TCEs may represent candidate targets because the d(N)/d(S) (D(N)/D(S)) values were small and negative selection was inferred for many of them.     

Codon-substitution models for detecting molecular adaptation at individual sites along specific lineages

The nonsynonymous (amino acid-altering) to synonymous (silent) substitution rate ratio (omega = d(N)/d(S)) provides a measure of natural selection at the protein level, with omega = 1, >1, and <1, indicating neutral evolution, purifying selection, and positive selection, respectively. Previous studies that used this measure to detect positive selection have often taken an approach of pairwise comparison, estimating substitution rates by averaging over all sites in the protein. As most amino acids in a functional protein are under structural and functional constraints and adaptive evolution probably affects only a few sites at a few time points, this approach of averaging rates over sites and over time has little power. Previously, we developed codon-based substitution models that allow the omega ratio to vary either among lineages or among sites. In this paper we extend previous models to allow the omega ratio to vary both among sites and among lineages and implement the new models in the likelihood framework. These models may be useful for identifying positive selection along prespecified lineages that affects only a few sites in the protein. We apply those branch-site models as well as previous branch- and site-specific models to three data sets: the lysozyme genes from primates, the tumor suppressor BRCA1 genes from primates, and the phytochrome (PHY) gene family in angiosperms. Positive selection is detected in the lysozyme and BRCA genes by both the new and the old models. However, only the new models detected positive selection acting on lineages after gene duplication in the PHY gene family. Additional tests on several data sets suggest that the new models may be useful in detecting positive selection after gene duplication in gene family evolution.     

Inferring natural selection operating on conservative and radical substitution at single amino acid sites

Natural selection operating on amino acid substitution at single amino acid sites can be detected by comparing the rates of synonymous (r(S)) and nonsynonymous (r(N)) nucleotide substitution at single codon sites. Amino acid substitutions can be classified as conservative or radical according to whether they retain the properties of the substituted amino acid. Here methods for comparing the rates of conservative (r(C)) and radical (r(R)) nonsynonymous substitution with r(S) at single codon sites were developed to detect natural selection operating on these substitutions at single amino acid sites. A method for comparing r(C) and r(R) at single codon sites was also developed to detect biases toward these substitutions at single amino acid sites. Charge was used as the property of the amino acids. In a computer simulation, false-positive rates of these methods were always < 5%, unless termination sites were included in the computation of the numbers of sites and estimates of transition/transversion rate ratio were highly biased. The frequency of detection of natural selection operating on conservative substitution was almost independent of the presence of natural selection operating on radical substitution, and vice versa. Natural selection operating specifically on conservative and radical substitution was detected more efficiently by comparing r(S) with r(C) and r(S) with r(R) than by comparing r(S) with r(N). These methods also appeared to be robust against the occurrence of recombination during evolution. In an analysis of class I human leukocyte antigen, negative selection operating on conservative substitution, but not positive selection operating on radical substitution, was observed at some of the codon sites with r(R) > r(C), suggesting that r(R) > r(C) may not necessarily be an indicator of positive selection operating on radical substitution.     

Maximum-likelihood analysis of molecular adaptation in abalone sperm lysin reveals variable selective pressures among lineages and sites

Maximum-likelihood models of codon substitution were used to analyze sperm lysin genes of 25 abalone (HALIOTIS:) species to identify lineages and amino acid sites under diversifying selection. The models used the nonsynonymous/synonymous rate ratio (omega = d(N)/d(S)) as an indicator of selective pressure and allowed the ratio to vary among lineages or sites. Likelihood ratio tests suggested significant variation in selective pressure among lineages. The variable selective pressure provided an explanation for the previous observation that the omega ratio is >1 in comparisons of closely related species and <1 in comparisons of distantly related species. Computer simulations demonstrated that saturation of nonsynonymous substitutions and constraint on lysin structure were unlikely to account for the observed pattern. Lineages linking closely related sympatric species appeared to be under diversifying selection, while lineages separating distantly related species from different geographic locations were associated with low evolutionary rates. The selective pressure indicated by the omega ratio was found to vary greatly among amino acid sites in lysin. Sites under potential diversifying selection were identified. Ancestral lysins were inferred to trace the route of evolution at individual sites and to provide lysin sequences for future laboratory studies.     

Diversification of olfactory receptor genes in the Japanese medaka fish, Oryzias latipes

Vertebrate olfactory receptors (OR) exists as the largest multigene family, scattered throughout the genome in clusters. Studies have shown that different animals possess remarkably diverse set of OR genes to recognize diverse odor molecules. In order to examine the evolutionary process of OR diversification, we examined three OR gene subfamilies from Japanese medaka fish (seven lines sampled from four populations). For each subfamily, the sequences of ancestral genes were inferred based on distance method. Examination of d(N)/d(S) ratios for each branch of phylogenetic trees suggested that purifying selection is the major force of evolution in medaka OR genes. However, for the mfOR1 and mfOR2 paralogous gene pairs, a nonrandom distribution of fixed amino acid changes and the d(N)>d(S) in a branch suggested that diversifying selection occurred after gene duplication. The fixed amino acid changes were observed in the third, fifth and sixth transmembrane domains, which has been predicted to serve as a ligand-binding pocket in a structural model. Compatibility test suggested that interlocus recombinations involving the fourth transmembrane domain occurred between the mfOR1 and mfOR2 gene pairs. The pattern of nucleotide substitutions in other OR genes agrees with the hypothesis that a limited number of amino acid residues are involved in odorant binding. Such comparative analyses of paralogous OR genes should provide bases for understanding the evolution, the structure, and the functional specificity of OR genes.     

The rate of adaptive evolution in enteric bacteria

Here we estimate the rate of adaptive substitution in a set of 410 genes that are present in 6 Escherichia coli and 6 Salmonella enterica genomes. We estimate that more than 50% of amino acid substitutions in this set of genes have been fixed by positive selection between the E. coli and S. enterica lineages. We also show that the proportion of adaptive substitutions is uncorrelated with the rate of amino acid substitution or gene function but that it may be correlated with levels of synonymous codon usage bias.     

Within-host evolution of CD8+-TL epitopes encoded by overlapping and non-overlapping reading frames of simian immunodeficiency virus

In order to understand the impact of overlapping reading frames on natural selection by host CD8+ T lymphocytes (CD8(+)-TL), we analyzed the pattern of nucleotide substitution in simian immunodeficiency virus (SIV) genomes sampled from populations at time of death in 35 rhesus monkeys. Both the mean number of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)) and the mean number of synonymous nucleotide substitutions per synonymous site (d(S)) were elevated in overlap regions in comparison to non-overlap regions. Mean d(N) exceeded mean d(S) in CD8(+)-TL epitopes restricted by the host's class I major histocompatibility complex molecules. This pattern, which is indicative of positive Darwinian selection favoring amino acid changes in these epitopes, was seen in both overlap and non-overlap regions; but mean d(N) was particularly elevated in restricted CD8(+)-TL epitopes encoded in overlap regions. Amino acid changes from the inoculum were defined as parallel if the same amino acid change occurred at the same site independently in two or more monkeys, and a surprisingly high proportion (71.9%) of observed amino acid changes throughout the SIV genome occurred in parallel in different monkeys. The proportion of parallel changes in restricted epitopes encoded by overlapping reading frames was still higher (80%), supporting the hypothesis that the interaction of positive selection and overlapping reading frames enhances the probability of convergent or parallel amino acid change.     

Genetic diversity of near genome-wide hepatitis C virus sequences during chronic infection: evidence for protein structural conservation over time

Infection with hepatitis C virus (HCV) is one of the leading causes of chronic hepatitis, liver cirrhosis and end-stage liver disease worldwide. The genetics of HCV infection in humans and the disease course of chronic hepatitis C are both remarkably variable. Although the response to interferon treatment is largely dependent on HCV genotypes, whether or not a relationship exists between HCV genome variability and clinical course of hepatitis C disease still remains unknown. To more thoroughly understand HCV genome evolution over time in association with disease course, near genome-wide HCV genomes present in 9 chronically infected participants over 83 total study years were sequenced. Overall, within HCV genomes, the number of synonymous substitutions per synonymous site (d(S)) significantly exceeded the number of non-synonymous substitutions per site (d(N)). Although both d(S) and d(N) significantly increased with duration of chronic infection, there was a highly significant decrease in d(N)/d(S) ratio in HCV genomes over time. These results indicate that purifying selection acted to conserve viral protein structure despite persistence of high level of nucleotide mutagenesis inherent to HCV replication. Based on liver biopsy fibrosis scores, HCV genomes from participants with advanced fibrosis had significantly greater d(S) values and lower d(N)/d(S) ratios compared to participants with mild liver disease. Over time, viral genomes from participants with mild disease had significantly greater annual changes in d(N), along with higher d(N)/d(S) ratios, compared to participants with advanced fibrosis. Yearly amino acid variations in the HCV p7, NS2, NS3 and NS5B genes were all significantly lower in participants with severe versus mild disease, suggesting possible pathogenic importance of protein structural conservation for these viral gene products.     

Molecular phylogenies and evolution of crt genes in algae

The carotenoids constitute the most widespread class of pigments in nature. Most previous work has concentrated on the identification and characterization of their chemical physical properties and bioavailability. In recent years, significant amounts of research have been conducted in an attempt to analyze the genes and the molecular regulation of the genes involved in the biosynthesis of carotenoids. However, it is important not to lose sight of the early evolution of carotenoid biosynthesis. One of the major obstacles in understanding the evolution of the respective enzymes and their patterns of selection is a lack of a well-supported phylogenic analysis. In the present research, a major long-term objective was to provide a clearer picture of the evolutionary history of genes, together with an evaluation of the patterns of selection in algae. These phylogenies will be important in studies characterizing the evolution of algae. The gene sequences of the enzymes involved in the major steps of the carotenoid biosynthetic pathway in algae (cyanobacteria, rhofophyta, chlorophyta) have been analyzed. Phylogenetic relationships among protein-coding DNA sequences were reconstructed by neighbor-joining (NJ) analysis for the respective carotenoid biosynthetic pathway genes (crt) in algae. The analysis also contains an estimation of the rate of nonsynonymous nucleotide substitutions per nonsynonymous site (d(N)), synonymous nucleotide substitution per synonymous site (d(S)), and the ratio of nonsynonmous (d(N)/d(S)) for the test of selection patterns. The phylogenetic trees show that the taxa of some genera have a closer evolutionary relationship with other genera in some gene sequences, which suggests a common ancient origin and that lateral gene transfer has occurred among unrelated genera. The d(N) values of crt genes in the early pathway are relatively low, while those of the following steps are slightly higher, while the d(N) values of crt genes in chlorophyta are higher than those in cyanobacteria. Most of the d(N)/d(S) values exceed 1. The phylogenetic analysis revealed that lateral gene transfer may have taken place across algal genomes and the d(N) values suggest that most of the early crt genes are well conserved compared to the later crt genes. Furthermore, d(N) values also revealed that the crt genes of chlorophyta are more evolutionary than cyanobacteria. The amino acids' changes are mostly adaptive evolution under the influence of positive diversity selection.