Pervasive positive selection on duplicated and nonduplicated vertebrate protein coding genes.

A stringent branch-site codon model was used to detect positive selection in vertebrate evolution. We show that the test is robust to the large evolutionary distances involved. Positive selection was detected in 77% of 884 genes studied. Most positive selection concerns a few sites on a single branch of the phylogenetic tree: Between 0.9% and 4.7% of sites are affected by positive selection depending on the branches. No functional category was overrepresented among genes under positive selection. Surprisingly, whole genome duplication had no effect on the prevalence of positive selection, whether the fish-specific genome duplication or the two rounds at the origin of vertebrates. Thus positive selection has not been limited to a few gene classes, or to specific evolutionary events such as duplication, but has been pervasive during vertebrate evolution.

A critical role for the T cell receptor alpha-chain connecting peptide domain in positive selection of CD1-independent NKT cells.

Natural killer T (NKT) cells are a subset of mature alpha beta TCR(+) cells that co-express NK lineage markers. Whereas most NKT cells express a canonical Valpha14/Vbeta8.2 TCR and are selected by CD1d, a minority of NKT cells express a diverse TCR repertoire and develop independently of CD1d. Little is known about the selection requirements of CD1d-independent NKT cells. We show here that NKT cells develop in RAG-deficient mice expressing an MHC class II-restricted transgenic TCR (Valpha2/Vbeta8.1) but only under conditions that lead to negative selection of conventional T cells. Moreover development of NKT cells in these mice is absolutely dependent upon an intact TCR alpha-chain connecting peptide domain, which is required for positive selection of conventional T cells via recruitment of the ERK signaling pathway. Collectively our data demonstrate that NKT cells can develop as a result of high avidity TCR/MHC class II interactions and suggest that common signaling pathways are involved in the positive selection of CD1d-independent NKT cells and conventional T cells.

Improving the performance of positive selection inference by filtering unreliable alignment regions.

Errors in the inferred multiple sequence alignment may lead to false prediction of positive selection. Recently, methods for detecting unreliable alignment regions were developed and were shown to accurately identify incorrectly aligned regions. While removing unreliable alignment regions is expected to increase the accuracy of positive selection inference, such filtering may also significantly decrease the power of the test, as positively selected regions are fast evolving, and those same regions are often those that are difficult to align. Here, we used realistic simulations that mimic sequence evolution of HIV-1 genes to test the hypothesis that the performance of positive selection inference using codon models can be improved by removing unreliable alignment regions. Our study shows that the benefit of removing unreliable regions exceeds the loss of power due to the removal of some of the true positively selected sites.

gcodeml: A Grid-enabled Tool for Detecting Positive Selection in Biological Evolution

One of the important questions in biological evolution is to know if certain changes along protein coding genes have contributed to the adaptation of species. This problem is known to be biologically complex and computationally very expensive. It, therefore, requires efficient Grid or cluster solutions to overcome the computational challenge. We have developed a Grid-enabled tool (gcodeml) that relies on the PAML (codeml) package to help analyse large phylogenetic datasets on both Grids and computational clusters. Although we report on results for gcodeml, our approach is applicable and customisable to related problems in biology or other scientific domains.

An adult thymic stromal-cell suspension model for in vitro positive selection.

Presented here is a cell-suspension model for positive selection using thymocytes from alphabeta-TCR (H-2Db-restricted) transgenic mice specific to the lymphocytic choriomeningitis virus (LCMV) on a nonselecting MHC background (H-2d or TAP-1 -/-), cocultured with freshly isolated adult thymus stromal cells of the selecting MHC type. The thymic stromal cells alone induced positive selection of functional CD4- CD8+ cells whose kinetics and efficiency were enhanced by nominal peptide. Fibroblasts expressing the selecting MHC alone did not induce positive selection; however, together with nonselecting stroma and nominal peptide, there was inefficient positive. These results suggest multiple signaling in positive selection with selection events able to occur on multiple-cell types. The ease with which this model can be manipulated should greatly facilitate the resolution of the mechanisms of positive selection in normal and pathological states.

Comparative genomics of chemosensory protein genes reveals rapid evolution and positive selection in ant-specific duplicates.

Gene duplications can have a major role in adaptation, and gene families underlying chemosensation are particularly interesting due to their essential role in chemical recognition of mates, predators and food resources. Social insects add yet another dimension to the study of chemosensory genomics, as the key components of their social life rely on chemical communication. Still, chemosensory gene families are little studied in social insects. Here we annotated chemosensory protein (CSP) genes from seven ant genomes and studied their evolution. The number of functional CSP genes ranges from 11 to 21 depending on species, and the estimated rates of gene birth and death indicate high turnover of genes. Ant CSP genes include seven conservative orthologous groups present in all the ants, and a group of genes that has expanded independently in different ant lineages. Interestingly, the expanded group of genes has a differing mode of evolution from the orthologous groups. The expanded group shows rapid evolution as indicated by a high dN/dS (nonsynonymous to synonymous changes) ratio, several sites under positive selection and many pseudogenes, whereas the genes in the seven orthologous groups evolve slowly under purifying selection and include only one pseudogene. These results show that adaptive changes have played a role in ant CSP evolution. The expanded group of ant-specific genes is phylogenetically close to a conservative orthologous group CSP7, which includes genes known to be involved in ant nestmate recognition, raising an interesting possibility that the expanded CSPs function in ant chemical communication.

TCR-alpha CDR3 loop audition regulates positive selection.

How positive selection molds the T cell repertoire has been difficult to examine. In this study, we use TCR-beta-transgenic mice in which MHC shapes TCR-alpha use. Differential AV segment use is directly related to the constraints placed on the composition of the CDR3 loops. Where these constraints are low, efficient selection of alphabeta pairs follows. This mode of selection preferentially uses favored AV-AJ rearrangements and promotes diversity. Increased constraint on the alpha CDR3 loops leads to inefficient selection associated with uncommon recombination events and limited diversity. Further, the two modes of selection favor alternate sets of AJ segments. We discuss the relevance of these findings to the imprint of self-MHC restriction and peripheral T cell activation.

Patterns of positive selection in seven ant genomes.

The evolution of ants is marked by remarkable adaptations that allowed the development of very complex social systems. To identify how ant-specific adaptations are associated with patterns of molecular evolution, we searched for signs of positive selection on amino-acid changes in proteins. We identified 24 functional categories of genes which were enriched for positively selected genes in the ant lineage. We also reanalyzed genome-wide data sets in bees and flies with the same methodology to check whether positive selection was specific to ants or also present in other insects. Notably, genes implicated in immunity were enriched for positively selected genes in the three lineages, ruling out the hypothesis that the evolution of hygienic behaviors in social insects caused a major relaxation of selective pressure on immune genes. Our scan also indicated that genes implicated in neurogenesis and olfaction started to undergo increased positive selection before the evolution of sociality in Hymenoptera. Finally, the comparison between these three lineages allowed us to pinpoint molecular evolution patterns that were specific to the ant lineage. In particular, there was ant-specific recurrent positive selection on genes with mitochondrial functions, suggesting that mitochondrial activity was improved during the evolution of this lineage. This might have been an important step toward the evolution of extreme lifespan that is a hallmark of ants.

Promoter IV of the class II transactivator gene is essential for positive selection of CD4+ T cells.

Major histocompatibility complex class II (MHCII) expression is regulated by the transcriptional coactivator CIITA. Positive selection of CD4(+) T cells is abrogated in mice lacking one of the promoters (pIV) of the Mhc2ta gene. This is entirely due to the absence of MHCII expression in thymic epithelia, as demonstrated by bone marrow transfer experiments between wild-type and pIV(-/-) mice. Medullary thymic epithelial cells (mTECs) are also MHCII(-) in pIV(-/-) mice. Bone marrow-derived, professional antigen-presenting cells (APCs) retain normal MHCII expression in pIV(-/-) mice, including those believed to mediate negative selection in the thymic medulla. Endogenous retroviruses thus retain their ability to sustain negative selection of the residual CD4(+) thymocytes in pIV(-/-) mice. Interestingly, the passive acquisition of MHCII molecules by thymocytes is abrogated in pIV(-/-) mice. This identifies thymic epithelial cells as the source of this passive transfer. In peripheral lymphoid organs, the CD4(+) T-cell population of pIV(-/-) mice is quantitatively and qualitatively comparable to that of MHCII-deficient mice. It comprises a high proportion of CD1-restricted natural killer T cells, which results in a bias of the V beta repertoire of the residual CD4(+) T-cell population. We have also addressed the identity of the signal that sustains pIV expression in cortical epithelia. We found that the Jak/STAT pathways activated by the common gamma chain (CD132) or common beta chain (CDw131) cytokine receptors are not required for MHCII expression in thymic cortical epithelia.

Selectome update: quality control and computational improvements to a database of positive selection.

Selectome (http://selectome.unil.ch/) is a database of positive selection, based on a branch-site likelihood test. This model estimates the number of nonsynonymous substitutions (dN) and synonymous substitutions (dS) to evaluate the variation in selective pressure (dN/dS ratio) over branches and over sites. Since the original release of Selectome, we have benchmarked and implemented a thorough quality control procedure on multiple sequence alignments, aiming to provide minimum false-positive results. We have also improved the computational efficiency of the branch-site test implementation, allowing larger data sets and more frequent updates. Release 6 of Selectome includes all gene trees from Ensembl for Primates and Glires, as well as a large set of vertebrate gene trees. A total of 6810 gene trees have some evidence of positive selection. Finally, the web interface has been improved to be more responsive and to facilitate searches and browsing.

The branch-site test of positive selection is surprisingly robust but lacks power under synonymous substitution saturation and variation in GC

Positive selection is widely estimated from protein coding sequence alignments by the nonsynonymous-to-synonymous ratio omega. Increasingly elaborate codon models are used in a likelihood framework for this estimation. Although there is widespread concern about the robustness of the estimation of the omega ratio, more efforts are needed to estimate this robustness, especially in the context of complex models. Here, we focused on the branch-site codon model. We investigated its robustness on a large set of simulated data. First, we investigated the impact of sequence divergence. We found evidence of underestimation of the synonymous substitution rate for values as small as 0.5, with a slight increase in false positives for the branch-site test. When dS increases further, underestimation of dS is worse, but false positives decrease. Interestingly, the detection of true positives follows a similar distribution, with a maximum for intermediary values of dS. Thus, high dS is more of a concern for a loss of power (false negatives) than for false positives of the test. Second, we investigated the impact of GC content. We showed that there is no significant difference of false positives between high GC (up to similar to 80%) and low GC (similar to 30%) genes. Moreover, neither shifts of GC content on a specific branch nor major shifts in GC along the gene sequence generate many false positives. Our results confirm that the branch-site is a very conservative test.

Optimization strategies for fast detection of positive selection on phylogenetic trees.

MOTIVATION: The detection of positive selection is widely used to study gene and genome evolution, but its application remains limited by the high computational cost of existing implementations. We present a series of computational optimizations for more efficient estimation of the likelihood function on large-scale phylogenetic problems. We illustrate our approach using the branch-site model of codon evolution. RESULTS: We introduce novel optimization techniques that substantially outperform both CodeML from the PAML package and our previously optimized sequential version SlimCodeML. These techniques can also be applied to other likelihood-based phylogeny software. Our implementation scales well for large numbers of codons and/or species. It can therefore analyse substantially larger datasets than CodeML. We evaluated FastCodeML on different platforms and measured average sequential speedups of FastCodeML (single-threaded) versus CodeML of up to 5.8, average speedups of FastCodeML (multi-threaded) versus CodeML on a single node (shared memory) of up to 36.9 for 12 CPU cores, and average speedups of the distributed FastCodeML versus CodeML of up to 170.9 on eight nodes (96 CPU cores in total).Availability and implementation: ftp://ftp.vital-it.ch/tools/FastCodeML/. CONTACT: selectome@unil.ch or nicolas.salamin@unil.ch.

Evolutionary trajectories of primate genes involved in HIV pathogenesis.

The current availability of five complete genomes of different primate species allows the analysis of genetic divergence over the last 40 million years of evolution. We hypothesized that the interspecies differences observed in susceptibility to HIV-1 would be influenced by the long-range selective pressures on host genes associated with HIV-1 pathogenesis. We established a list of human genes (n = 140) proposed to be involved in HIV-1 biology and pathogenesis and a control set of 100 random genes. We retrieved the orthologous genes from the genome of humans and of four nonhuman primates (Pan troglodytes, Pongo pygmaeus abeli, Macaca mulatta, and Callithrix jacchus) and analyzed the nucleotide substitution patterns of this data set using codon-based maximum likelihood procedures. In addition, we evaluated whether the candidate genes have been targets of recent positive selection in humans by analyzing HapMap Phase 2 single-nucleotide polymorphisms genotyped in a region centered on each candidate gene. A total of 1,064 sequences were used for the analyses. Similar median K(A)/K(S) values were estimated for the set of genes involved in HIV-1 pathogenesis and for control genes, 0.19 and 0.15, respectively. However, genes of the innate immunity had median values of 0.37 (P value = 0.0001, compared with control genes), and genes of intrinsic cellular defense had K(A)/K(S) values around or greater than 1.0 (P value = 0.0002). Detailed assessment allowed the identification of residues under positive selection in 13 proteins: AKT1, APOBEC3G, APOBEC3H, CD4, DEFB1, GML, IL4, IL8RA, L-SIGN/CLEC4M, PTPRC/CD45, Tetherin/BST2, TLR7, and TRIM5alpha. A number of those residues are relevant for HIV-1 biology. The set of 140 genes involved in HIV-1 pathogenesis did not show a significant enrichment in signals of recent positive selection in humans (intraspecies selection). However, we identified within or near these genes 24 polymorphisms showing strong signatures of recent positive selection. Interestingly, the DEFB1 gene presented signatures of both interspecies positive selection in primates and intraspecies recent positive selection in humans. The systematic assessment of long-acting selective pressures on primate genomes is a useful tool to extend our understanding of genetic variation influencing contemporary susceptibility to HIV-1.

Targeted expression of human CD1d in transgenic mice reveals independent roles for thymocytes and thymic APCs in positive and negative selection of Valpha14i NKT cells.

CD1d-dependent invariant Valpha14 (Valpha14i) NKT cells are innate T lymphocytes expressing a conserved semi-invariant TCR, consisting, in mice, of the invariant Valpha14-Jalpha18 TCR alpha-chain paired mostly with Vbeta8.2 and Vbeta7. The cellular requirements for thymic positive and negative selection of Valpha14i NKT cells are only partially understood. Therefore, we generated transgenic mice expressing human CD1d (hCD1d) either on thymocytes, mainly CD4+ CD8+ double positive, or on APCs, the cells implicated in the selection of Valpha14i NKT cells. In the absence of the endogenous mouse CD1d (mCD1d), the expression of hCD1d on thymocytes, but not on APCs, was sufficient to select Valpha14i NKT cells that proved functional when activated ex vivo with the Ag alpha-galactosyl ceramide. Valpha14i NKT cells selected by hCD1d on thymocytes, however, attained lower numbers than in control mice and expressed essentially Vbeta8.2. The low number of Vbeta8.2+ Valpha14i NKT cells selected by hCD1d on thymocytes was not reversed by the concomitant expression of mCD1d, which, instead, restored the development of Vbeta7+ Valpha14i NKT cells. Vbeta8.2+, but not Vbeta7+, NKT cell development was impaired in mice expressing both hCD1d on APCs and mCD1d. Taken together, our data reveal that selective CD1d expression by thymocytes is sufficient for positive selection of functional Valpha14i NKT cells and that both thymocytes and APCs may independently mediate negative selection.