Maximum likelihood analysis of mammalian p53 indicates the presence of positively selected sites and higher tumorigenic mutations in purifying sites

Large Scale Analyses and Visualization of Adaptive Amino Acid Changes Projects

When changes at few amino acid sites are the target of selection, adaptive amino acid changes in protein sequences can be identified using maximum-likelihood methods based on models of codon substitution (such as codeml). Although such methods have been employed numerous times using a variety of different organisms, the time needed to collect the data and prepare the input files means that tens or hundreds of coding regions are usually analyzed. Nevertheless, the recent availability of flexible and easy to use computer applications that collect relevant data (such as BDBM) and infer positively selected amino acid sites (such as ADOPS), means that the entire process is easier and quicker than before. However, the lack of a batch option in ADOPS, here reported, still precludes the analysis of hundreds or thousands of sequence files. Given the interest and possibility of running such large-scale projects, we have also developed a database where ADOPS projects can be stored. Therefore, this study also presents the B+ database, which is both a data repository and a convenient interface that looks at the information contained in ADOPS projects without the need to download and unzip the corresponding ADOPS project file. The ADOPS projects available at B+ can also be downloaded, unzipped, and opened using the ADOPS graphical interface. The availability of such a database ensures results repeatability, promotes data reuse with significant savings on the time needed for preparing datasets, and effortlessly allows further exploration of the data contained in ADOPS projects.

The P72R Polymorphism of p53 Predisposes to Obesity and Metabolic Dysfunction

p53 is well known for its tumor suppressor role, but this protein also has a poorly understood role in the regulation of metabolism. Human studies have implicated a common polymorphism at codon 72 of p53 in diabetic and pre-diabetic phenotypes. To understand this role, we utilized a humanized mouse model of the p53 codon 72 variants and monitored these mice following challenge with a high-fat diet (HFD). Mice with the arginine 72 (R72) variant of p53 developed more-severe obesity and glucose intolerance on a HFD, compared to mice with the proline 72 variant (P72). R72 mice developed insulin resistance, islet hypertrophy, increased infiltration of immune cells, and fatty liver disease. Gene expression analyses and studies with small-molecule inhibitors indicate that the p53 target genes Tnf and Npc1l1 underlie this phenotype. These results shed light on the role of p53 in obesity, metabolism, and inflammation.

p53 gene discriminates two ecologically divergent sister species of pine voles

Genes with relevant roles in the differentiation of closely-related species are likely to have diverged simultaneously with the species and more accurately reproduce the species tree. The Lusitanian (Microtus lusitanicus) and Mediterranean (M. duodecimcostatus) pine voles are two recently separated sister species with fossorial lifestyles whose different ecological, physiological and morphological phenotypes reflect the better adaptation of M. duodecimcostatus to the underground habitat. Here we asked whether the differentiation of M. lusitanicus and M. duodecimcostatus involved genetic variations within the tumour suppressor p53 gene, given its role in stress-associated responses. We performed a population-genetic analysis through sequencing of exons and introns of p53 in individuals from sympatric and allopatric populations of both the species in the Iberian Peninsula in which a unidirectional introgression of mitochondrial DNA was previously observed. We were able to discriminate the two species to a large extent. We show that M. duodecimcostatus is composed of one genetically unstructured group of populations sharing a P53 protein that carries a mutation in the DNA-binding region not observed in M. lusitanicus, raising the possibility that this mutation may have been central in the evolutionary history of M. duodecimcostatus. Our results provide suggestive evidence for the involvement of a master transcription factor in the separation of M. lusitanicus and M. duodecimcostatus during Microtus radiation in the Quaternary presumably via a differential adaptive role of the novel p53 in M. duodecimcostatus.

Divergent positive selection in rhodopsin from lake and riverine cichlid fishes

Studies of cichlid evolution have highlighted the importance of visual pigment genes in the spectacular radiation of the African rift lake cichlids. Recent work, however, has also provided strong evidence for adaptive diversification of riverine cichlids in the Neotropics, which inhabit environments of markedly different spectral properties from the African rift lakes. These ecological and/or biogeographic differences may have imposed divergent selective pressures on the evolution of the cichlid visual system. To test these hypotheses, we investigated the molecular evolution of the dim-light visual pigment, rhodopsin. We sequenced rhodopsin from Neotropical and African riverine cichlids and combined these data with published sequences from African cichlids. We found significant evidence for positive selection using random sites codon models in all cichlid groups, with the highest levels in African lake cichlids. Tests using branch-site and clade models that partitioned the data along ecological (lake, river) and/or biogeographic (African, Neotropical) boundaries found significant evidence of divergent selective pressures among cichlid groups. However, statistical comparisons among these models suggest that ecological, rather than biogeographic, factors may be responsible for divergent selective pressures that have shaped the evolution of the visual system in cichlids. We found that branch-site models did not perform as well as clade models for our data set, in which there was evidence for positive selection in the background. One of our most intriguing results is that the amino acid sites found to be under positive selection in Neotropical and African lake cichlids were largely nonoverlapping, despite falling into the same three functional categories: spectral tuning, retinal uptake/release, and rhodopsin dimerization. Taken together, these results would imply divergent selection across cichlid clades, but targeting similar functions. This study highlights the importance of molecular investigations of ecologically important groups and the flexibility of clade models in explicitly testing ecological hypotheses.

ADOPS - Automatic Detection Of Positively Selected Sites

Maximum-likelihood methods based on models of codon substitution have been widely used to infer positively selected amino acid sites that are responsible for adaptive changes. Nevertheless, in order to use such an approach, software applications are required to align protein and DNA sequences, infer a phylogenetic tree and run the maximum-likelihood models. Therefore, a significant effort is made in order to prepare input files for the different software applications and in the analysis of the output of every analysis. In this paper we present the ADOPS (Automatic Detection Of Positively Selected Sites) software. It was developed with the goal of providing an automatic and flexible tool for detecting positively selected sites given a set of unaligned nucleotide sequence data. An example of the usefulness of such a pipeline is given by showing, under different conditions, positively selected amino acid sites in a set of 54 Coffea putative S-RNase sequences. ADOPS software is freely available and can be downloaded from http://sing.ei.uvigo.es/ADOPS.

Codon substitution models based on residue similarity and their applications

Codon models are now widely used to draw evolutionary inferences from alignments of homologous sequence data. Incorporating physicochemical properties of amino acids into codon models, two novel codon substitution models describing the evolution of protein-coding DNA sequences are presented based on the similarity scores of amino acids. To describe substitutions between codons a continue-time Markov process is used. Transition/transversion rate bias and nonsynonymous codon usage bias are allowed in the models. In our implementation, the parameters are estimated by maximum-likelihood (ML) method as in previous studies. Furthermore, instantaneous mutations involving more than one nucleotide position of a codon are considered in the second model. Then the two suggested models are applied to five real data sets. The analytic results indicate that the new codon models considering physicochemical properties of amino acids can provide a better fit to the data comparing with existing codon models, and then produce more reliable estimates of certain biologically important measures than existing methods.