Structure,diversity, and evolution of the T-cell receptor VB gene repertoire in primates

AB T-cell receptors (TCR) that recognize major histocompatibility complex (MHC)/peptide antigen complexes regulate humoral and cellular arms of the adaptive immune response. Antigen binding sites of MHC and immunoglobulin heavy chain variable regions(Igh-V) are subject to diversity enhancing selection. We sought to establish whether positive Darwinian selection has driven diversity of TCRBV chains in the primate lineage by sequencing rearranged TCR from rhesus monkeys and chimpanzees and comparing them with those of humans. Rates of synonymous (silent) and nonsynonymous (replacement) substitutions indicate selection against amino acid replacements in TCRBV frameworks, and relaxation of these constraints in putative MHC/peptide contact sites. The lack of positive selection for variability in likely ligand contact sites suggests that mechanisms generating somatic diversity in TCR junctional regions have relaxed the pressure for selection of variability in the TCR V region encoded in the germline.     

Detecting the form of selection in the outer membrane protein C of Enterobacter aerogenes strains and Salmonella species

The types of selective pressure operating on the outer membrane protein C (ompC) of Enterobacter aerogenes strains, the causative agent for nosocomial infections, and Salmonella sp., the hazardous pathogen are investigated using the maximum likelihood-based codon substitution models. Although the rate of amino acid replacement to the silent substitution (omega) across the entire codon sites of ompC of E. aerogenes (omega=0.3194) and Salmonella sp. (omega=0.2047) indicate that the gene is subjected to purifying selection (i.e. omega<1), approximately 3.7% of ompC codon sites in E. aerogenes (omega=21.52) are under the influence of positive Darwinian selection (i.e. omega>1). Such contrast in the intensity of selective pressures in both pathogens could be associated with the differential response to the adverse environmental changes. In E. aerogenes, majority of the positively selected sites are located in the hypervariable cell-surface-exposed domains whereas the trans-membrane domains are functionally highly constrained.     

Codon-substitution models for heterogeneous selection pressure at amino acid sites

Comparison of relative fixation rates of synonymous (silent) and nonsynonymous (amino acid-altering) mutations provides a means for understanding the mechanisms of molecular sequence evolution. The nonsynonymous/synonymous rate ratio (omega = d(N)d(S)) is an important indicator of selective pressure at the protein level, with omega = 1 meaning neutral mutations, omega < 1 purifying selection, and omega > 1 diversifying positive selection. Amino acid sites in a protein are expected to be under different selective pressures and have different underlying omega ratios. We develop models that account for heterogeneous omega ratios among amino acid sites and apply them to phylogenetic analyses of protein-coding DNA sequences. These models are useful for testing for adaptive molecular evolution and identifying amino acid sites under diversifying selection. Ten data sets of genes from nuclear, mitochondrial, and viral genomes are analyzed to estimate the distributions of omega among sites. In all data sets analyzed, the selective pressure indicated by the omega ratio is found to be highly heterogeneous among sites. Previously unsuspected Darwinian selection is detected in several genes in which the average omega ratio across sites is <1, but in which some sites are clearly under diversifying selection with omega > 1. Genes undergoing positive selection include the beta-globin gene from vertebrates, mitochondrial protein-coding genes from hominoids, the hemagglutinin (HA) gene from human influenza virus A, and HIV-1 env, vif, and pol genes. Tests for the presence of positively selected sites and their subsequent identification appear quite robust to the specific distributional form assumed for omega and can be achieved using any of several models we implement. However, we encountered difficulties in estimating the precise distribution of omega among sites from real data sets.     

Evaluation of whether accelerated protein evolution in chordates has occurred before, after, or simultaneously with gene duplication

Gene duplication and loss are predicted to be at least of the order of the substitution rate and are key contributors to the development of novel gene function and overall genome evolution. Although it has been established that proteins evolve more rapidly after gene duplication, we were interested in testing to what extent this reflects causation or association. Therefore, we investigated the rate of evolution prior to gene duplication in chordates. Two patterns emerged; firstly, branches, which are both preceded by a duplication and followed by a duplication, display an elevated rate of amino acid replacement. This is reflected in the ratio of nonsynonymous to synonymous substitution (mean nonsynonymous to synonymous nucleotide substitution rate ratio [Ka:Ks]) of 0.44 compared with branches preceded by and followed by a speciation (mean Ka:Ks of 0.23). The observed patterns suggest that there can be simultaneous alteration in the selection pressures on both gene duplication and amino acid replacement, which may be consistent with co-occurring increases in positive selection, or alternatively with concurrent relaxation of purifying selection. The pattern is largely, but perhaps not completely, explained by the existence of certain families that have elevated rates of both gene duplication and amino acid replacement. Secondly, we observed accelerated amino acid replacement prior to duplication (mean Ka:Ks for postspeciation preduplication branches was 0.27). In some cases, this could reflect adaptive changes in protein function precipitating a gene duplication event. In conclusion, the circumstances surrounding the birth of new proteins may frequently involve a simultaneous change in selection pressures on both gene-copy number and amino acid replacement. More precise modeling of the relative importance of preduplication, postduplication, and simultaneous amino acid replacement will require larger and denser genomic data sets from multiple species, allowing simultaneous estimation of lineage-specific fluctuations in mutation rates and adaptive constraints.     

The ratio of replacement to silent divergence and tests of neutrality

Comparisons of replacement to silent divergence have been used in a variety of studies aimed at detecting selection. Here, such comparisons are shown to be very sensitive to the pattern of rate variation in replacement sites. Saturation may play an important role even at surprisingly low levels of divergence if the substitution rate varies across replacement sites. For example, saturation in replacement sites may be of importance in the evolution of the HIV-1 envelope gene. However, the pattern of saturation in replacement and silent sites may, in itself, provide valuable insight into the causes of DNA evolution. 210 DNA sequences from 15 different loci/systematic groups are analyzed, and evidence for positive selection is demonstrated in at least one of these data sets, through an analysis of the distribution of substitution rates along the sequence.     

Adaptive evolution after duplication of penaeidin antimicrobial peptides

Penaeidin antimicrobial peptides in penaeid shrimps are an important component of their innate immune system that provides immunity against infection caused by several gram-positive bacteria and filamentous fungal species. Despite the knowledge on the identification and characterization of these peptides in penaeid shrimps, little is known about the evolutionary pattern of these peptides and the underlying genetic mechanisms that maintain high sequence diversities in the penaeidin gene family. Based on the phylogenetic analyses and maximum likelihood-based codon substitution analyses, here we present the convincing evidence that multiple copies of penaeidins have evolved by gene duplication, and positive Darwinian selection (adaptive evolution) is the likely cause of accelerated rate of amino acid substitutions among these duplicated genes. While the average ratio of non-synonymous to synonymous substitutions (omega) for the entire coding region of both active domains is 0.9805, few codon sites showed significantly higher omega (3.73). The likelihood ratio tests that compare models incorporating positive selection (omega>1) at certain codon sites with models not incorporating positive selection (omega<1), failed to reject (p=0) the evidence of positive Darwinian selection. The rapid adaptive evolution of this gene family might be directed by the pathogens and the faster rate of amino acid substitutions in the N-terminal proline-rich and C-terminal cysteine-rich domains could be due to their direct involvement in the protection against pathogens. When the host expose to different habitats/environment an accelerated rate of amino acid substitutions in both the active domains may also be expected.     

Maximum likelihood methods for detecting adaptive evolution after gene duplication

The rapid accumulation of genomic sequences in public databases will finally allow large scale studies of gene family evolution, including evaluation of the role of positive Darwinian selection following a duplication event. This will be possible because recent statistical methods of comparing synonymous and nonsynonymous substitution rates permit reliable detection of positive selection at individual amino acid sites and along evolutionary lineages. Here, we summarize maximum-likelihood based methods, and present a framework for their application to analysis of gene families. Using these methods, we investigated the role of positive Darwinian selection in the ECP-EDN gene family of primates and the Troponin C gene family of vertebrates. We also comment on the limitations of these methods and discuss directions for further improvements.     

Positive Darwinian Selection Drives the Evolution of the Morphology-Related Gene, EPCAM, in Particularly Species-Rich Lineages of African Cichlid Fishes

The study of genetic evolution within the context of adaptive radiations offers insights to genes and selection pressures that result in rapid morphological change. Cichlid fishes are very species-rich and variable in coloration, behavior, and morphology, and so provide a classical model system for studying the genetics of adaptive radiation. In this study, we researched the evolution of the epithelial cell adhesion molecule (EPCAM), a candidate gene for the adaptive evolution of morphology broadly, and skin development specifically, in fishes. We compared EPCAM gene sequences from a rapidly speciating African cichlid lineage (the haplochromines), a species-poor African lineage (Nile tilapia Oreochromis niloticus), and a very young adaptive radiation in the Neotropics (sympatric crater lake Midas cichlids, Amphilophus sp.). Our results, based on a hierarchy of evolutionary analyses of nucleotide substitution, demonstrate that there are different selection pressures on the EPCAM gene among the cichlid lineages. Several waves of positive natural selection were identified not only on the terminal branches, but also on ancestral branches. Interestingly, significant positive or directional selection was found in the haplochromine cichlids only but not the comparatively species-poor tilapia lineage. We hypothesize that the strong signal of selection in the ancestral African cichlid lineage coincided with the transition from riverine to lacustrine habitat. The two neotropical species for which we collected new sequence data were invariant in the EPCAM locus. Our results suggest that functional changes promoted by positive Darwinian selection are widespread in the EPCAM gene during African cichlid evolution.     

Gene duplication,allelic diversity,selection processes and adaptive value of MHC class II <Emphasis Type="Italic">DRB</Emphasis> genes of the bank vole, <Emphasis Type="Italic">Clethrionomys glareolus</Emphasis>

The generation and maintenance of allelic polymorphism in genes of the major histocompatibility complex (MHC) is a central issue in evolutionary genetics. Recently, the focus has changed from ex situ to in situ populations to understand the mechanisms that determine adaptive MHC polymorphism under natural selection. Birth-and-death evolution and gene conversion events are considered to generate sequence diversity in MHC genes, which subsequently is maintained by balancing selection through parasites. The ongoing arms race between the host and parasites leads to an adaptive selection pressure upon the MHC, evident in high rates of non-synonymous vs synonymous substitution rates. We characterised the MHC class II DRB exon 2 of free living bank voles, Clethrionomys glareolus by single-strand conformation polymorphism and direct sequencing. Unlike other arvicolid species, the DRB locus of the bank vole is at least quadruplicated. No evidence for gene conversion events in the Clgl-DRB sequences was observed. We found not only high allelic polymorphism with 26 alleles in 36 individuals but also high rates of silent polymorphism. Exceptional for MHC class II genes is a purifying selection pressure upon the majority of MHC-DRB sequences. Further, we analysed the association between certain DRB alleles and the parasite burden with gastrointestinal trichostrongyle nematodes Heligmosomum mixtum and Heligmosomoides glareoli and found significant quality differences between specific alleles with respect to infection intensity. Our findings suggest a snapshot in an evolutionary process of ongoing birth-and-death evolution. One allele cluster has lost its function and is already silenced, another is loosing its adaptive value in terms of gastrointestinal nematode resistance, while a third group of alleles indicates all signs of classical functional MHC alleles.     

Gene expression intensity shapes evolutionary rates of the proteins encoded by the vertebrate genome

Natural selection leaves its footprints on protein-coding sequences by modulating their silent and replacement evolutionary rates. In highly expressed genes in invertebrates, these footprints are seen in the higher codon usage bias and lower synonymous divergence. In mammals, the highly expressed genes have a shorter gene length in the genome and the breadth of expression is known to constrain the rate of protein evolution. Here we have examined how the rates of evolution of proteins encoded by the vertebrate genomes are modulated by the amount (intensity) of gene expression. To understand how natural selection operates on proteins that appear to have arisen in earlier and later phases of animal evolution, we have contrasted patterns of mouse proteins that have homologs in invertebrate and protist genomes (Precambrian genes) with those that do not have such detectable homologs (vertebrate-specific genes). We find that the intensity of gene expression relates inversely to the rate of protein sequence evolution on a genomic scale. The most highly expressed genes actually show the lowest total number of substitutions per polypeptide, consistent with cumulative effects of purifying selection on individual amino acid replacements. Precambrian genes exhibit a more pronounced difference in protein evolutionary rates (up to three times) between the genes with high and low expression levels as compared to the vertebrate-specific genes, which appears to be due to the narrower breadth of expression of the vertebrate-specific genes. These results provide insights into the differential relationship and effect of the increasing complexity of animal body form on evolutionary rates of proteins.     

EP300 contributes to high-altitude adaptation in Tibetans by regulating nitric oxide production

The genetic adaptation of Tibetans to high altitude hypoxia likely involves a group of genes in the hypoxic pathway,as suggested by earlier studies.To test the adaptive role of the previously reported candidate gene EP300 (histone acetyltransferase p300),we conducted resequencing of a 108.9 kb gene region of EP300 in 80 unrelated Tibetans.The allele-frequency and haplotype-based neutrality tests detected signals of positive Darwinian selection on EP300 in Tibetans,with a group of variants showing allelic divergence between Tibetans and lowland reference populations,including Han Chinese,Europeans,and Africans.Functional prediction suggested the involvement of multiple EP300 variants in gene expression regulation.More importantly,genetic association tests in 226 Tibetans indicated significant correlation of the adaptive EP300 variants with blood nitric oxide (NO) concentration.Collectively,we propose that EP300 harbors adaptive variants in Tibetans,which might contribute to high-altitude adaptation through regulating NO production.     

Adaptive molecular evolution of a defence gene in sexual but not functionally asexual evening primroses

Theory predicts that sexual reproduction provides evolutionary advantages over asexual reproduction by reducing mutational load and increasing adaptive potential. Here, we test the latter prediction in the context of plant defences against pathogens because pathogens frequently reduce plant fitness and drive the evolution of plant defences. Specifically, we ask whether sexual evening primrose plant lineages (Onagraceae) have faster rates of adaptive molecular evolution and altered gene expression of a class I chitinase, a gene implicated in defence against pathogens, than functionally asexual evening primrose lineages. We found that the ratio of amino acid to silent substitutions (K(a) /K(s) = 0.19 vs. 0.11 for sexual and asexual lineages, respectively), the number of sites identified to be under positive selection (four vs. zero for sexual and asexual lineages, respectively) and the expression of chitinase were all higher in sexual than in asexual lineages. Our results are congruent with the conclusion that a loss of sexual recombination and segregation in the Onagraceae negatively affects adaptive structural and potentially regulatory evolution of a plant defence protein.     

Varying signals of the effects of natural selection during teleost growth hormone gene evolution.

The growth hormone (GH) gene of teleost fish exhibits a higher degree of variability compared with other vertebrate groups. However, the different selective constraints at the sequence level are not well understood. In this study, maximum-likelihood (ML) models of codon substitutions were used to investigate Darwinian adaptive evolution of the GH gene in teleost fishes. Complete GH gene sequences of 54 fish species were classified into 4 orders, and the variable nature of GH was examined by determining the dN and dS rate variation and the rates of molecular evolution for each teleost order. The results indicate that although the overall evolution rate for teleost GH is high ((1.15 +/- 0.01) x 10(-9) substitutions/(aa site x y)) compared with the "slow phases" in mammals ((0.21 to 0.28 +/- 0.05) x 10(-9)), the vital structure of this gene has been retained. While the majority of the amino acid changes appear to be due to relaxation of purifying selection, some positively selected sites were detected in regions with no specifically identified role in protein function. The positively selected regions observed in salmoniformes lineage suggests a possible role for positive selection driving functional divergence in paralogous forms of the GH gene after whole-genome duplication in this lineage.     

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.     

Molecular Adaptation of rbcL in the Heterophyllous Aquatic Plant Potamogeton

Background

Heterophyllous aquatic plants show marked phenotypic plasticity. They adapt to environmental changes by producing different leaf types: submerged, floating and terrestrial leaves. By contrast, homophyllous plants produce only submerged leaves and grow entirely underwater. Heterophylly and submerged homophylly evolved under selective pressure modifying the species-specific optima for photosynthesis, but little is known about the evolutionary outcome of habit. Recent evolutionary analyses suggested that rbcL, a chloroplast gene that encodes a catalytic subunit of RuBisCO, evolves under positive selection in most land plant lineages. To examine the adaptive evolutionary process linked to heterophylly or homophylly, we analyzed positive selection in the rbcL sequences of ecologically diverse aquatic plants, Japanese Potamogeton.

Principal Findings

Phylogenetic and maximum likelihood analyses of codon substitution models indicated that Potamogeton rbcL has evolved under positive Darwinian selection. The positive selection has operated specifically in heterophyllous lineages but not in homophyllous ones in the branch-site models. This suggests that the selective pressure on this chloroplast gene was higher for heterophyllous lineages than for homophyllous lineages. The replacement of 12 amino acids occurred at structurally important sites in the quaternary structure of RbcL, two of which (residue 225 and 281) were identified as potentially under positive selection.

Conclusions/Significance

Our analysis did not show an exact relationship between the amino acid replacements and heterophylly or homophylly but revealed that lineage-specific positive selection acted on the Potamogeton rbcL. The contrasting ecological conditions between heterophyllous and homophyllous plants have imposed different selective pressures on the photosynthetic system. The increased amino acid replacement in RbcL may reflect the continuous fine-tuning of RuBisCO under varying ecological conditions.     

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.     

Substitution rates in a new Silene latifolia sex-linked gene, SlssX/Y

Dioecious white campion Silene latifolia has sex chromosomal sex determination, with homogametic (XX) females and heterogametic (XY) males. This species has become popular in studies of sex chromosome evolution. However, the lack of genes isolated from the X and Y chromosomes of this species is a major obstacle for such studies. Here, I report the isolation of a new sex-linked gene, Slss, with strong homology to spermidine synthase genes of other species. The new gene has homologous intact copies on the X and Y chromosomes (SlssX and SlssY, respectively). Synonymous divergence between the SlssX and SlssY genes is 4.7%, and nonsynonymous divergence is 1.4%. Isolation of a homologous gene from nondioecious S. vulgaris provided a root to the gene tree and allowed the estimation of the silent and replacement substitution rates along the SlssX and SlssY lineages. Interestingly, the Y-linked gene has higher synonymous and nonsynonymous substitution rates. The elevated synonymous rate in the SlssY gene, compared with SlssX, confirms our previous suggestion that the S. latifolia Y chromosome has a higher mutation rate, compared with the X chromosome. When differences in silent substitution rate are taken into account, the Y-linked gene still demonstrates significantly faster accumulation of nonsynonymous substitutions, which is consistent with the theoretical prediction of relaxed purifying selection in Y-linked genes, leading to the accumulation of nonsynonymous substitutions and genetic degeneration of the Y-linked genes.     

Molecular Evolution in a Multisite Nearly Neutral Mutation Model

A simple nearly neutral mutation model of protein evolution was studied using computer simulation assuming a constant population size. In this model, a gene consists of a finite number of codons and there is no recombination within a gene. Each codon has two replacement and one silent sites. The fitness of a gene was determined multiplicatively by amino acids specified by codons (the independent multicodon model). Nucleotide diversity at replacement sites decreases as selection becomes stronger. A reduction of nucleotide diversity at silent sites also occurs as selection intensifies but the magnitude of the reduction is not a monotone function of the intensity of selection. The dispersion index is close to one. The average value of Tajima's and Fu and Li's statistics are negative and their absolute values increases as selection intensifies. However, their powers of detecting selection under the present model were not high unless the number of sites is large or mutation rate is high. The MK test was shown to detect intermediate selection fairly well. For comparison, the house-of-cards model was also investigated and its behavior was shown to be more sensitive to changes of population size than that of the independent multicodon model. The relevance of the present model for explaining protein evolution was discussed comparing its prediction and recent DNA data. Received: 24 May 1999 / Accepted: 17 August 1999     

The genomic rate of molecular adaptation of the human influenza A virus

Quantifying adaptive evolution at the genomic scale is an essential yet challenging aspect of evolutionary biology. Here, we develop a method that extends and generalizes previous approaches to estimate the rate of genomic adaptation in rapidly evolving populations and apply it to a large data set of complete human influenza A virus genome sequences. In accord with previous studies, we observe particularly high rates of adaptive evolution in domain 1 of the viral hemagglutinin (HA1). However, our novel approach also reveals previously unseen adaptation in other viral genes. Notably, we find that the rate of adaptation (per codon per year) is higher in surface residues of the viral neuraminidase than in HA1, indicating strong antibody-mediated selection on the former. We also observed high rates of adaptive evolution in several nonstructural proteins, which may relate to viral evasion of T-cell and innate immune responses. Furthermore, our analysis provides strong quantitative support for the hypothesis that human H1N1 influenza experiences weaker antigenic selection than H3N2. As well as shedding new light on the dynamics and determinants of positive Darwinian selection in influenza viruses, the approach introduced here is applicable to other pathogens for which densely sampled genome sequences are available, and hence is ideally suited to the interpretation of next-generation genome sequencing data.