Relaxed Functional Constraints Research Articles

Peculiar patterns of amino acid substitution and conservation in the fast evolving tunicate Oikopleura dioica

We analyze the patterns and rates of amino acid evolution in tunicates with special interest on the extremely fast evolving Oikopleura dioica. We show that this species, on average, is twice as fast as the already fast evolving Ciona intestinalis. The acceleration in both species seems to be affected by similar evolutionary forces yet to different extent, since a substantial proportion of the most and less accelerated genes are orthologous between the two species.Among the possible causes that underlie the genome wide acceleration in Oikopleura, relaxation of functional constraints appears to be an important one, since all amino acids exhibit surprisingly homogenous levels of divergence. Such homogeneity, however, is not observed in Ciona.Apart from the genome wide acceleration, detailed analysis of functional groups of genes revealed that genes associated with regulatory functions (transcription regulators, chromatin remodeling proteins and metabolic regulators), have been subjected to an even more extreme process of acceleration, suggesting that adaptive evolution is the most probable cause of their unusual exacerbated rates.Another remarkable observation is that cysteine is among the less conserved amino acids, contrary to what is commonly observed in other species. The possible causes of this particular behavior are discussed.

Degeneration of an Intracellular Ion Channel in the Primate Lineage by Relaxation of Selective Constraints

Ion channel genes are highly conserved and are rarely degenerated in the primate lineage leading to humans. So far, the only well-characterized ion channel known to be degenerated in primates is the plasma membrane transient receptor potential channel TRPC2, possibly due to changes in the pheromone signaling. Here, by analyzing the sequence data from ten primate species, we have determined the degeneration process of the TPC3 gene that encodes a member of the two-pore channel (TPC) family recently implicated in Ca(2+) release by nicotinic acid adenine dinucleotide phosphate from intracellular acidic stores in animals. We show that degeneration of TPC3 likely began in the common ancestors of Apes and Old World monkeys through a conserved inactivating mutation, followed by additional deleterious mutations resulting in the generation of a TPC3 pseudogene in the descendant catarrhine lineage. Located at a chromosome recombination hot spot, catarrhine TPC3 pseudogenes underwent a series of lineage-specific rearrangements, including exon deletion and duplication. In contrast, we identify near full-length TPC3 sequences in New World monkeys and Prosimians and show that the gene is subjected to strong purifying selection and therefore likely functional. Our data provide the first evidence for relaxed functional constraints for an intracellular ion channel in primates and shed novel insights into the evolution and regulation of Ca(2+) signaling in the primate lineage.

Microevolutionary dynamics of a macroevolutionary key innovation in a Lepidopteran herbivore

BackgroundA molecular population genetics understanding is central to the study of ecological and evolutionary functional genomics. Population genetics identifies genetic variation and its distribution within and among populations, it reveals the demographic history of the populations studied, and can provide indirect insights into historical selection dynamics. Here we use this approach to examine the demographic and selective dynamics acting of a candidate gene involved in plant-insect interactions. Previous work documents the macroevolutionary and historical ecological importance of the nitrile-specifier protein (Nsp), which facilitated the host shift of Pieridae butterflies onto Brassicales host plants ~80 Myr ago.ResultsHere we assess the microevolutionary dynamics of the Nsp gene by studying the within and among-population variation at Nsp and reference genes in the butterfly Pieris rapae (Small Cabbage White). Nsp exhibits unexpectedly high amounts of amino acid polymorphism, unequally distributed across the gene. The vast majority of genetic variation exists within populations, with little to no genetic differentiation among four populations on two continents. A comparison of synonymous and nonsynonymous substitutions in 70 randomly chosen genes among P. rapae and its close relative Pieris brassicae (Large Cabbage White) finds Nsp to have a significantly relaxed functional constraint compared to housekeeping genes. We find strong evidence for a recent population expansion and no role for strong purifying or directional selection upon the Nsp gene.ConclusionsThe microevolutionary dynamics of the Nsp gene in P. rapae are dominated by recent population expansion and variation in functional constraint across the repeated domains of the Nsp gene. While the high amounts of amino acid diversity suggest there may be significant functional differences among allelic variants segregating within populations, indirect tests of selection could not conclusively identify a signature of historical selection. The importance of using this information for planning future studies of potential performance and fitness consequences of the observed variation is discussed.

Novel method for discerning the action of selection during evolution

A common problem in molecular comparative geno- mics is the identification of genes that are under positive, adaptive selection [1]. Such genes are likely to be crucial for speciation, species differentiation, and func- tional specialization. However, discerning the difference between positive selection and relaxation of func- tional constraint can be difficult using current methods. Both processes generally increase the rate of ami- no acid change relative to synonymous changes within coding regions, and unless the amino acid rate is over- whelmingly high across an entire gene, the signature of positive selection can be obscured [2]. Some methodologies do not explicitly determine the difference be- tween a relaxation of functional constraint and positive selection, leaving researchers to determine via other means whether the trajectory of a gene has been specialization or creation of a new function, or removal from the genome via a process of degeneration.

Extensive Functional Diversification of thePopulusGlutathioneS-Transferase Supergene Family

Identifying how genes and their functions evolve after duplication is central to understanding gene family radiation. In this study, we systematically examined the functional diversification of the glutathione S-transferase (GST) gene family in Populus trichocarpa by integrating phylogeny, expression, substrate specificity, and enzyme kinetic data. GSTs are ubiquitous proteins in plants that play important roles in stress tolerance and detoxification metabolism. Genome annotation identified 81 GST genes in Populus that were divided into eight classes with distinct divergence in their evolutionary rate, gene structure, expression responses to abiotic stressors, and enzymatic properties of encoded proteins. In addition, when all the functional parameters were examined, clear divergence was observed within tandem clusters and between paralogous gene pairs, suggesting that subfunctionalization has taken place among duplicate genes. The two domains of GST proteins appear to have evolved under differential selective pressures. The C-terminal domain seems to have been subject to more relaxed functional constraints or divergent directional selection, which may have allowed rapid changes in substrate specificity, affinity, and activity, while maintaining the primary function of the enzyme. Our findings shed light on mechanisms that facilitate the retention of duplicate genes, which can result in a large gene family with a broad substrate spectrum and a wide range of reactivity toward different substrates.

Loss of teeth and enamel in tetrapods: fossil record, genetic data and morphological adaptations.

Since their recruitment in the oral cavity, approximately 450 million years ago, teeth have been subjected to strong selective constraints due to the crucial role that they play in species survival. It is therefore quite surprising that the ability to develop functional teeth has subsequently been lost several times, independently, in various lineages. In this review, we concentrate our attention on tetrapods, the only vertebrate lineage in which several clades lack functional teeth from birth to adulthood. Indeed, in other lineages, teeth can be absent in adults but be functionally present in larvae and juveniles, can be absent in the oral cavity but exist in the pharyngeal region, or can develop on the upper jaw but be absent on the lower jaw. Here, we analyse the current data on toothless (edentate) tetrapod taxa, including information available on enamel-less species. Firstly, we provide an analysis of the dispersed and fragmentary morphological data published on the various living taxa concerned (and their extinct relatives) with the aim of tracing the origin of tooth or enamel loss, i.e. toads in Lissamphibia, turtles and birds in Sauropsida, and baleen whales, pangolins, anteaters, sloths, armadillos and aardvark in Mammalia. Secondly, we present current hypotheses on the genetic basis of tooth loss in the chicken and thirdly, we try to answer the question of how these taxa have survived tooth loss given the crucial importance of this tool. The loss of teeth (or only enamel) in all of these taxa was not lethal because it was always preceded in evolution by the pre-adaptation of a secondary tool (beak, baleens, elongated adhesive tongues or hypselodonty) useful for improving efficiency in food uptake. The positive selection of such secondary tools would have led to relaxed functional constraints on teeth and would have later compensated for the loss of teeth. These hypotheses raise numerous questions that will hopefully be answered in the near future.

Evolution of X-Degenerate Y Chromosome Genes in Greater Apes: Conservation of Gene Content in Human and Gorilla, But Not Chimpanzee

Compared with the X chromosome, the mammalian Y chromosome is considerably diminished in size and has lost most of its ancestral genes during evolution. Interestingly, for the X-degenerate region on the Y chromosome, human has retained all 16 genes, while chimpanzee has lost 4 of the 16 genes since the divergence of the two species. To uncover the evolutionary forces governing ape Y chromosome degeneration, we determined the complete sequences of the coding exons and splice sites for 16 gorilla Y chromosome genes of the X-degenerate region. We discovered that all studied reading frames and splice sites were intact, and thus, this genomic region experienced no gene loss in the gorilla lineage. Higher nucleotide divergence was observed in the chimpanzee than the human lineage, particularly for genes with disruptive mutations, suggesting a lack of functional constraints for these genes in chimpanzee. Surprisingly, our results indicate that the human and gorilla orthologues of the genes disrupted in chimpanzee evolve under relaxed functional constraints and might not be essential. Taking mating patterns and effective population sizes of ape species into account, we conclude that genetic hitchhiking associated with positive selection due to sperm competition might explain the rapid decline in the Y chromosome gene number in chimpanzee. As we found no evidence of positive selection acting on the X-degenerate genes, such selection likely targets other genes on the chimpanzee Y chromosome.

Genome evolution: Gene duplication and the resolution of adaptive conflict

Genome evolution: Gene duplication and the resolution of adaptive conflict

Comparative genomic analysis of teleost fish bmal genes

Bmal1 (Brain and muscle ARNT like 1) gene is a key circadian clock gene. Tetrapods also have the second Bmal gene, Bmal2. Fruit fly has only one bmal1/cycle gene. Interrogation of the five teleost fish genome sequences coupled with phylogenetic and splice site analyses found that zebrafish have two bmal1 genes, bmal1a and bmal1b, and bmal2a; Japanese pufferfish (fugu), green spotted pufferfish (tetraodon) and Japanese medaka fish each have two bmal2 genes, bmal2a and bmal2b, and bmal1a; and three-spine stickleback have bmal1a and bmal2b. Syntenic analysis further indicated that zebrafish bmal1a/bmal1b, and fugu, tetraodon and medaka bmal2a/bmal2b are ancient duplicates. Although the dN/dS ratios of these four fish bmal duplicates are all <1, implicating they have been under purifying selection, the Tajima relative rate test showed that fugu, tetraodon and medaka bmal2a/bmal2b have asymmetric evolutionary rates, suggesting that one of these duplicates have been subject to positive selection or relaxed functional constraint. These results support the notion that teleost fish bmal genes were derived from the fish-specific genome duplication (FSGD), divergent resolution following the duplication led to retaining different ancient bmal duplicates in different fishes, which could have shaped the evolution of the complex teleost fish timekeeping mechanisms.

Evolution of Soldier-Specific Venomous Protease in Social Aphids

In social aphids of the genus Tuberaphis a cysteine protease gene of the family cathepsin B exhibits soldier-specific expression and intestinal protease production. The product is orally excreted and injected by soldier nymphs into natural enemies, thereby exerting an insecticidal activity. In an attempt to gain insights into when and how the novel venomous protease for the altruistic caste has evolved, we investigated the soldier-specific type (S-type) and nonspecific type (N-type) cathepsin B genes from social and nonsocial aphids. All the social aphids examined, representing the genera Tuberaphis, Astegopteryx, and Cerataphis, possessed both the S-type and N-type genes. Phylogenetically distant nonsocial aphids also possessed cathepsin B genes allied to the S-type and the N-type, indicating the evolutionary origin of these genes in the common ancestor of extant aphids. In Tuberaphis species the S-type genes exhibited significant soldier-specific expression and accelerated molecular evolution whereas the N-type genes did not. In Astegopteryx and Cerataphis species, meanwhile, both the S-type and N-type genes exhibited neither remarkable soldier-specific expression nor accelerated molecular evolution. These results suggest that the S-type gene acquired the soldier-specific expression and the venom function after divergence of the genus Tuberaphis. On the structural model of the S-type protease of Tuberaphis styraci the accelerated molecular evolution was associated with the molecular surface rather than the catalytic cleft, suggesting that the venom activity was probably acquired by relatively minor modifications on the molecular surface rather than by generation of a novel active site. In Cerataphis jamuritsu the S-type gene was, although containing a stop codon, structurally almost intact and still transcribed, suggesting recent pseudogenization of the gene copy and possible relevance to relaxed functional constraint in the highly multiplied protease gene family. On the basis of these results we suggest that the massive amplification in aphid cathepsin B genes might have predisposed the evolution of venomous protease in the social aphid lineage and argue that gene duplication, accelerated molecular evolution, and acquisition of novel gene function must have played considerable roles in the evolution of complex biological systems including insect sociality.

Hominoid seminal protein evolution and ancestral mating behavior

Hominoid mating systems show extensive variation among species. The degree of sexual dimorphism in body size and canine size varies among primates in accordance with their mating system, as does the testes size and the consistency of ejaculated semen, in response to differing levels of sperm competition. To investigate patterns of evolution at hominoid seminal proteins and to make inferences regarding the mating systems of extinct taxa, we sequenced the entire coding region of the prostate-specific transglutaminase (TGM4) gene in human, chimpanzee, bonobo, western lowland gorilla, eastern lowland gorilla, orangutan, and siamang, including multiple humans, chimps, and gorillas. Partial DNA sequence of the coding regions was also obtained for one eastern lowland gorilla at the semenogelin genes (SEMG1 and SEMG2), which code for the predominant proteins in semen. Patterns of nucleotide variation and inferred protein sequence change were evaluated within and between species. Combining the present data with previous studies demonstrates a high rate of amino acid substitutions, and low intraspecific variation, at seminal proteins in Pan, presumably driven by strong sperm competition. Both gorilla species apparently possess nonfunctional TGM4, SEMG1, and SEMG2 genes, suggesting that gorillas have had low sperm competition, and therefore their current polygynous mating system, for a long time before their divergence. Similarly, orangutans show longstanding stasis at TGM4, which may be interpreted as evidence for an unchanging mating system for most of their evolution after their divergence from African apes. In contrast to the great apes, the data from humans could be interpreted as evidence of fluctuations between different mating systems or alternatively as a relaxed functional constraint in these proteins. It is our hope that this study is a first step toward developing a model to predict ancestral mating systems from extant molecular data to complement interpretations from the fossil record.

Comparative Analysis of Period Genes in Teleost Fish Genomes

Period (Per) is a canonical circadian clock gene. The fruit fly, an invertebrate, has one per gene, while the human, a tetrapod vertebrate, has three Per genes. Per1, Per2, and Per3 of the tetrapods were generated from two rounds of ancient genome duplications from the ancestral chordate Per gene. Searching for five teleost fish genomes in a combination of phylogenetic, splicing site, and syntenic analyses revealed that zebrafish have two per1 genes, per1a and per1b, one per2, and one per3; medaka, fugu, and tetraodon each have two per2 genes, per2a and per2b, one per1, and one per3; sticklebacks also have per2a, per2b, and one per1 but lack per3; and per1a/per1b in zebrafish and per2a/per2b in madaka, fugu, tetraodon, and stickleback are ancient duplicates. While the dN/dS ratios of the five fish per duplicates are all <1, suggesting that they likely have been subject to purifying selection, the Tajima relative rate test showed that zebrafish per1a/per1b and fugu and medaka per2a/per2b have asymmetric evolutionary rates, implicating that one of these duplicates might have been under positive selection or relaxed functional constraint. Further, in situ hybridization showed that zebrafish per1a and per1b clearly have distinct patterns of temporal and spatial expression. These results support the notion that extra copies of teleost per genes were generated from the fish-specific genome duplication, and divergent resolution after the duplication resulted in retention of different per duplicates in different fish, most of which have diverged significantly.

Comparative analysis of teleost fish genomes reveals preservation of different ancient clock duplicates in different fishes

Clock ( Circadian locomotor output cycle kaput) was the first vertebrate circadian clock gene identified in a mouse forward genetics mutagenesis screen. It encodes a bHLH-PAS protein that is highly conserved throughout evolution. Tetrapods also have the second Clock gene, Clock2 or Npas2 ( Neuronal PAS domain protein 2). Conversely, the fruit fly, an invertebrate, has only one clock gene. Interrogation of the five teleost fish genome databases revealed that the zebrafish and the Japanese pufferfish (fugu) each have three clock genes, whereas the green spotted pufferfish (tetraodon), the Japanese medaka fish and the three-spine stickleback each have two clock genes. Phylogenetic and splice site analyses indicated that zebrafish and fugu each have two clock1 genes, clock1a and clock1b and one clock2; tetraodon also have clock1a and clock1b but do not have clock2; and medaka and stickleback each have clock1b and one clock2. Genome neighborhood analysis further showed that clock1a/ clock1b in zebrafish, fugu and tetraodon is an ancient duplicate. While the dN/dS ratios of these three fish clock duplicates are all < 1, indicating that purifying selection has acted upon them; the Tajima relative rate test showed that all three fish clock duplicates have asymmetric evolutionary rates, implicating that one of these duplicates have been under positive selection or relaxed functional constraint. These results support the view that teleost fish clock genes were generated from an ancient genome-wide duplication, and differential gene loss after the duplication resulted in retention of different ancient duplicates in different teleost fishes, which could have contributed to the evolution of the distinct fish circadian clock mechanisms.

Molecular anatomy of the cytoplasmic domain of bovine growth hormone receptor, a quantitative trait locus

Quantitative trait loci (QTL) studies have indicated growth hormone receptor (GHR) as a candidate gene affecting cattle milk yield and composition. In order to characterize genetic variation at GHR in cattle, we studied European and East African breeds with different histories of selection, and Bos grunniens, Ovis aries, Sus scrofa, Bison bison and Rangifer tarandus as references. We sequenced most of the cytoplasmic domain (900 bp of exon 10), 89 bp of exon 8, including the putative causative mutation for the QTL effect, and 390 bp of intron 8 for comparison. In the cytoplasmic domain, seven synonymous and seven non-synonymous single nucleotide polymorphisms (SNP) were identified in cattle. Three non-synonymous SNPs were found in sheep and one synonymous SNP in yak, while other studied species were monomorphic. Three major haplotypes were observed, one unique to African breeds, one unique to European breeds and one shared. Bison and yak haplotypes are derivatives of the European haplotype lineage. Most of the exon 10 non-synonymous cattle SNPs appear at phylogenetically highly conserved sites. The polymorphisms in exon 10 cluster around a ruminant-specific tyrosine residue, suggesting that this site may act as an additional signalling domain of GHR in ruminants. Alternative explanations for the persistent polymorphism include balancing selection, hitch-hiking, pleiotropic or sexually antagonistic fitness effects or relaxed functional constraints.

Mitochondrial Genome Sequence Evolution in Chlamydomonas

The mitochondrial genomes of the Chlorophyta exhibit significant diversity with respect to gene content and genome compactness; however, quantitative data on the rates of nucleotide substitution in mitochondrial DNA, which might help explain the origin of this diversity, are lacking. To gain insight into the evolutionary forces responsible for mitochondrial genome diversification, we sequenced to near completion the mitochondrial genome of the chlorophyte Chlamydomonas incerta, estimated the evolutionary divergence between Chlamydomonas reinhardtii and C. incerta mitochondrial protein-coding genes and rRNA-coding regions, and compared the relative evolutionary rates in mitochondrial and nuclear genes. Synonymous and nonsynonymous substitution rates do not differ significantly between the mitochondrial and nuclear protein-coding genes. The mitochondrial rRNA-coding regions, however, are evolving much faster than their nuclear counterparts, and this difference might be explained by relaxed functional constraints on the mitochondrial translational apparatus due to the small number of proteins synthesized in Chlamydomonas mitochondria. Substitution rates at synonymous sites in a nonstandard mitochondrial gene (rtl) and at intronic and synonymous sites in nuclear genes expressed at low levels suggest that the mutation rate is similar in these two genetic compartments. Potential evolutionary forces shaping mitochondrial genome evolution in Chlamydomonas are discussed.

Mammalian specific mouse genes are evolving faster than mouse genes conserved across other eukaryotic lineages

Positive selection is usually considered in the context of a higher rate of substitutions in non-synonymous as compared to synonymous sites in complete coding sequences of genes or individual positions. We show that genes conserved in eukaryota, coelomata, and bilateria, that is, proteins that arose earlier in evolution as compared to mammalia specific genes evolve slowly and are subjected to negative selection. This finding supports the notion that evolutionary rates progressively diminish with the age of a gene. The data suggests that in both intron-containing and intronless genes synonymous sites may be subject to some degree of selection that is indicative of a relative acceleration of amino-acid substitution, which could be due to a relaxation of functional constraints and/or directional selection.

Different patterns of selection on the nuclear genes IRBP and DMP-1 affect the efficiency but not the outcome of phylogeny estimation for didelphid marsupials

Selection at the protein-level can influence nucleotide substitution patterns for protein-coding genes, which in turn can affect their performance as phylogenetic characters. In this study, we compare two protein-coding nuclear genes that appear to have evolved under markedly different selective constraints and evaluate how selection has shaped their phylogenetic signal. We sequenced 1100+ bp of exon 6 of the gene encoding dentin matrix protein 1 (DMP1) from most of the currently recognized genera of New World opossums (family: Didelphidae) and compared these data to an existing matrix of sequences from the interphotoreceptor retinoid-binding protein gene (IRBP) and morphological characters. In comparison to IRBP, DMP1 has far fewer sites under strong purifying selection and exhibits a number of sites under positive directional selection. Furthermore, selection on the DMP1 protein appears to conserve short, acidic, serine-rich domains rather than primary amino acid sequence; as a result, DMP1 has significantly different nucleotide substitution patterns from IRBP. Using Bayesian methods, we determined that DMP1 evolves almost 30% faster than IRBP, has 2.5 times more variable sites, has less among-site rate heterogeneity, is skewed toward A and away from CT (IRBP has relatively even base frequencies), and has a significantly lower rate of change between adenine and any other nucleotide. Despite these different nucleotide substitution patterns, estimates of didelphid relationships based on separate phylogenetic analyses of these genes are remarkably congruent whether patterns of nucleotide substitution are explicitly modeled or not. Nonetheless, DMP1 contains more phylogenetically informative characters per unit sequence and resolves more nodes with higher support than does IRBP. Thus, for these two genes, relaxed functional constraints and positive selection appear to improve the efficiency of phylogenetic estimation without compromising its accuracy.

Characterization of nuclear 18S rRNA gene sequence diversity and expression in an individual lake sturgeon (Acipenser fulvescens)

The unexpected discovery of extensive intraindividual variation of the 18S rRNA gene in nine North American species of sturgeon (and its absence in the North American paddlefish) was described in a previous paper by Krieger and Fuerst (2002b). As part of the study, genomic 18S rRNA genes from a lake sturgeon (Acipenser fulvescens) individual were amplified by polymerase chain reaction (PCR), cloned and sequenced, resulting in the identification of at least 17 different sequence variants. Because of the number of variants detected and the improbability that such a large number of different alleles are required for survival, these sequence variants were subjected to further analyses (sequence comparison, phylogenetic and relative rate) to examine the possibility that some variants may be pseudogenes. The cDNA produced by reverse transcriptase (RT)-PCR amplification of 18S rRNA isolated from lake sturgeon liver tissue was also sequenced. The topology of the phylogenetic tree produced suggests a split of the sequence variants into two groups: paddlefish-like (PL) alleles and non-paddlefish-like (NPL) alleles. Relative rate comparisons of these two groups indicate that the alleles within the NPL allele group are evolving more quickly than those in the paddlefish-like (PL) group, and so are likely under relaxed functional constraints. Less than one-third of nucleotide changes in the non-paddlefish-like (NPL) alleles occur at positions that are considered to be highly conserved in a universal eukaryotic model of rRNA structure. No substitutions at eukaryotic universally conserved sites occur in the PL alleles. Analysis of the sequence of the RT-PCR product shows that a single lake sturgeon 18S rDNA sequence variant is expressed in major quantities in the liver, and that this single product is identical to the allele most similar in sequence to the paddlefish. These results provide support for the idea that many variants are unimportant for proper cellular function or are pseudogenes. Although concerted evolution has apparently failed to homogenize the many rDNA sequence variants found within sturgeon, the detection of a single major expressed sequence variant in lake sturgeon indicates that the expression of the variants is in fact under the control of selective factors.

Diversification of furanocoumarin-metabolizing cytochrome P450 monooxygenases in two papilionids: Specificity and substrate encounter rate.

Diversification of cytochrome P450 monooxygenases (P450s) is thought to result from antagonistic interactions between plants and their herbivorous enemies. However, little direct evidence demonstrates the relationship between selection by plant toxins and adaptive changes in herbivore P450s. Here we show that the furanocoumarin-metabolic activity of CYP6B proteins in two species of swallowtail caterpillars is associated with the probability of encountering host plant furanocoumarins. Catalytic activity was compared in two closely related CYP6B4 and CYP6B17 groups in the polyphagous congeners Papilio glaucus and Papilio canadensis. Generally, P450s from P. glaucus, which feeds occasionally on furanocoumarin-containing host plants, display higher activities against furanocoumarins than those from P. canadensis, which normally does not encounter furanocoumarins. These P450s in turn catalyze a larger range of furanocoumarins at lower efficiency than CYP6B1, a P450 from Papilio polyxenes, which feeds exclusively on furanocoumarin-containing host plants. Reconstruction of the ancestral CYP6B sequences using maximum likelihood predictions and comparisons of the sequence and geometry of their active sites to those of contemporary CYP6B proteins indicate that host plant diversity is directly related to P450 activity and inversely related to substrate specificity. These predictions suggest that, along the lineage leading to Papilio P450s, the ancestral, highly versatile CYP6B protein presumed to exist in a polyphagous species evolved through time into a more efficient and specialized CYP6B1-like protein in Papilio species with continual exposure to furanocoumarins. Further diversification of Papilio CYP6Bs has likely involved interspersed events of positive selection in oligophagous species and relaxation of functional constraints in polyphagous species.

Detecting excess radical replacements in phylogenetic trees

There are a few instances in which positive Darwinian selection has been convincingly demonstrated at the molecular level. In this study, we present a novel test for detecting excess of radical amino-acid replacements. Such excess is usually indicative of positive Darwinian selection, but may also be due to relaxed functional constraints or model misspecification. In our test, each amino-acid replacement is characterized in terms of a physicochemical distance, i.e., the degree of dissimilarity between the exchanged amino-acid residues. By using phylogenetic trees based on protein sequences, our test identifies statistically significant deviations of the mean physicochemical distance from the random expectation, either along a taxonomic lineage or across a subtree. The mean inferred distance is calculated as the average physicochemical distance over all possible ancestral sequence reconstructions weighted by their likelihood. Our method substantially improves over previous approaches by taking into account the stochastic process, tree phylogeny, among-site rate variation, and alternative ancestral reconstructions. We provide a fast linear time algorithm for applying this test to all branches and all subtrees of a given phylogenetic tree. We validate this approach by applying it to two well-studied datasets: the MHC class I glycoproteins serving as a positive control, and the house-keeping gene carbonic anhydrase I serving as a negative control.