Looking for Darwin in all the wrong places: the misguided quest for positive selection at the nucleotide sequence level

Evolution of communication protocols using an artificial regulatory network

I describe the Utrecht Machine (UM), a discrete artificial regulatory network designed for studying how evolution discovers biochemical computation mechanisms. The corresponding binary genome format is compatible with gene deletion, duplication, and recombination. In the simulation presented here, an agent consisting of two UMs, a sender and a receiver, must encode, transmit, and decode a binary word over time using the narrow communication channel between them. This communication problem has chicken-and-egg structure in that a sending mechanism is useless without a corresponding receiving mechanism. An in-depth case study reveals that a coincidence creates a minimal partial solution, from which a sequence of partial sending and receiving mechanisms evolve. Gene duplications contribute by enlarging the regulatory network. Analysis of 60,000 sample runs under a variety of parameter settings confirms that crossover accelerates evolution, that stronger selection tends to find clumsier solutions and finds them more slowly, and that there is implicit selection for robust mechanisms and genomes at the codon level. Typical solutions associate each input bit with an activation speed and combine them almost additively. The parents of breakthrough organisms sometimes have lower fitness scores than others in the population, indicating that populations can cross valleys in the fitness landscape via outlying members. The simulation exhibits back mutations and population-level memory effects not accounted for in traditional population genetics models. All together, these phenomena suggest that new evolutionary models are needed that incorporate regulatory network structure.

Impact of recurrent gene duplication on adaptation of plant genomes

BACKGROUND

Recurrent gene duplication and retention played an important role in angiosperm genome evolution. It has been hypothesized that these processes contribute significantly to plant adaptation but so far this hypothesis has not been tested at the genome scale.

RESULTS

We studied available sequenced angiosperm genomes to assess the frequency of positive selection footprints in lineage specific expanded (LSE) gene families compared to single-copy genes using a dN/dS-based test in a phylogenetic framework. We found 5.38% of alignments in LSE genes with codons under positive selection. In contrast, we found no evidence for codons under positive selection in the single-copy reference set. An analysis at the branch level shows that purifying selection acted more strongly on single-copy genes than on LSE gene clusters. Moreover we detect significantly more branches indicating evolution under positive selection and/or relaxed constraint in LSE genes than in single-copy genes.

CONCLUSIONS

In this - to our knowledge -first genome-scale study we provide strong empirical support for the hypothesis that LSE genes fuel adaptation in angiosperms. Our conservative approach for detecting selection footprints as well as our results can be of interest for further studies on (plant) gene family evolution.

Distinguishing the effects of selection from demographic history in the genetic variation of two sister passerines based on mitochondrial-nuclear comparison

Determining the mechanisms responsible for the distribution of genetic diversity in natural populations has occupied a central role in molecular evolution. Our study was motivated by the unprecedented observation that a widespread Eurasian flycatcher, Ficedula albicilla, exhibited no variation at the mitochondrial DNA (mtDNA) ND2 gene in 75 individuals sampled over a 5000-km distance. In contrast, its sister species, F. parva, had low but considerably higher levels of mtDNA variation. We assessed whether natural selection or demographic factors could explain the absence of mtDNA variation in F. albicilla. Eighteen nuclear genes were sequenced to estimate the two species' phylogeographic histories, and for comparison to the mtDNA data. Multilocus coalescence analyses suggested that F. albicilla experienced a population expansion perhaps following a population bottleneck. Simulations based on this demographic history, however, did not replicate the extremely low level of mtDNA variation. Historical range changes based on ecological niche models also failed to explain the observed mtDNA patterns. Neutrality tests (DHEW and ML-HKA) suggested a non-neutral pattern in the mtDNA of F. albicilla. We found a transmembrane-skewed distribution of nonsynonymous substitutions between the two species, three of which caused functional change; the results implied that positive selection could have targeted mtDNA. Several lines of evidence support selection rather than demographic history as the main force influencing the patterns of mtDNA variation. Despite the influence of natural selection, many of the phylogeographic inferences derived from mtDNA were robust, including species limits and a high level of gene flow among populations within species.

Generation of divergent uroplakin tetraspanins and their partners during vertebrate evolution: identification of novel uroplakins

BACKGROUND

The recent availability of sequenced genomes from a broad array of chordates (cephalochordates, urochordates and vertebrates) has allowed us to systematically analyze the evolution of uroplakins: tetraspanins (UPK1a and UPK1b families) and their respective partner proteins (UPK2 and UPK3 families).

RESULTS

We report here: (1) the origin of uroplakins in the common ancestor of vertebrates, (2) the appearance of several residues that have statistically significantly positive dN/dS ratios in the duplicated paralogs of uroplakin genes, and (3) the existence of strong coevolutionary relationships between UPK1a/1b tetraspanins and their respective UPK2/UPK3-related partner proteins. Moreover, we report the existence of three new UPK2/3 family members we named UPK2b, 3c and 3d, which will help clarify the evolutionary relationships between fish, amphibian and mammalian uroplakins that may perform divergent functions specific to these different and physiologically distinct groups of vertebrates.

CONCLUSIONS

Since our analyses cover species of all major chordate groups this work provides an extremely clear overall picture of how the uroplakin families and their partner proteins have evolved in parallel. We also highlight several novel features of uroplakin evolution including the appearance of UPK2b and 3d in fish and UPK3c in the common ancestor of reptiles and mammals. Additional studies of these novel uroplakins should lead to new insights into uroplakin structure and function.

Speciation in fungal and oomycete plant pathogens

The process of speciation, by definition, involves evolution of one or more reproductive isolating mechanisms that split a single species into two that can no longer interbreed. Determination of which processes are responsible for speciation is important yet challenging. Several studies have proposed that speciation in pathogens is heavily influenced by host-pathogen dynamics and that traits that mediate such interactions (e.g., host mobility, reproductive mode of the pathogen, complexity of the life cycle, and host specificity) must lead to reproductive isolation and ultimately affect speciation rates. In this review, we summarize the main evolutionary processes that lead to speciation of fungal and oomycete plant pathogens and provide an outline of how speciation can be studied rigorously, including novel genetic/genomic developments.

The effects of linkage on comparative estimators of selection

Background

A major goal of molecular evolution is to determine how natural selection has shaped the evolution of a gene. One approach taken by methods such as K_A/K_S and the McDonald-Kreitman (MK) test is to compare the frequency of non-synonymous and synonymous changes. These methods, however, rely on the assumption that a change in frequency of one mutation will not affect changes in frequency of other mutations.

Results

We demonstrate that linkage between sites can bias measures of selection based on synonymous and non-synonymous changes. Using forward simulation of a Wright-Fisher process, we show that hitch-hiking of deleterious mutations with advantageous mutations can lead to overestimation of the number of adaptive substitutions, while background selection and clonal interference can distort the site frequency spectrum to obscure the signal for positive selection. We present three diagnostics for detecting these effects of linked selection and apply them to the human influenza (H3N2) hemagglutinin gene.

Conclusion

Various forms of linked selection have characteristic effects on MK-type statistics. The extent of background selection, hitch-hiking and clonal interference can be evaluated using the diagnostic statistics presented here. The diagnostics can also be used to determine how well we expect the MK statistics to perform and whether one form of the statistic may be preferable to another.

Molecular evolutionary analysis of vertebrate transducins: a role for amino acid variation in photoreceptor deactivation

Transducin is a heterotrimeric G protein that plays a critical role in phototransduction in the rod and cone photoreceptor cells of the vertebrate retina. Rods, highly sensitive cells that recover from photoactivation slowly, underlie dim-light vision, whereas cones are less sensitive, recover more quickly, and underlie bright-light vision. Transducin deactivation is a critical step in photoreceptor recovery and may underlie the functional distinction between rods and cones. Rods and cones possess distinct transducin α subunits, yet they share a common deactivation mechanism, the GTPase activating protein (GAP) complex. Here, we used codon models to examine patterns of sequence evolution in rod (GNAT1) and cone (GNAT2) α subunits. Our results indicate that purifying selection is the dominant force shaping GNAT1 and GNAT2 evolution, but that GNAT2 has additionally been subject to positive selection operating at multiple phylogenetic scales; phylogeny-wide analysis identified several sites in the GNAT2 helical domain as having substantially elevated dN/dS estimates, and branch-site analysis identified several nearby sites as targets of strong positive selection during early vertebrate history. Examination of aligned GNAT and GAP complex crystal structures revealed steric clashes between several positively selected sites and the deactivating GAP complex. This suggests that GNAT2 sequence variation could play an important role in adaptive evolution of the vertebrate visual system via effects on photoreceptor deactivation kinetics and provides an alternative perspective to previous work that focused instead on the effect of GAP complex concentration. Our findings thus further the understanding of the molecular biology, physiology, and evolution of vertebrate visual systems.

Evolutionary biochemistry: revealing the historical and physical causes of protein properties

The repertoire of proteins and nucleic acids in the living world is determined by evolution; their properties are determined by the laws of physics and chemistry. Explanations of these two kinds of causality - the purviews of evolutionary biology and biochemistry, respectively - are typically pursued in isolation, but many fundamental questions fall squarely at the interface of fields. Here we articulate the paradigm of evolutionary biochemistry, which aims to dissect the physical mechanisms and evolutionary processes by which biological molecules diversified and to reveal how their physical architecture facilitates and constrains their evolution. We show how an integration of evolution with biochemistry moves us towards a more complete understanding of why biological molecules have the properties that they do.

Directional Darwinian Selection in proteins

Background

Molecular evolution is a very active field of research, with several complementary approaches, including dN/dS, HON90, MM01, and others. Each has documented strengths and weaknesses, and no one approach provides a clear picture of how natural selection works at the molecular level. The purpose of this work is to present a simple new method that uses quantitative amino acid properties to identify and characterize directional selection in proteins.

Methods

Inferred amino acid replacements are viewed through the prism of a single physicochemical property to determine the amount and direction of change caused by each replacement. This allows the calculation of the probability that the mean change in the single property associated with the amino acid replacements is equal to zero (H₀: μ = 0; i.e., no net change) using a simple two-tailed t-test.

Results

Example data from calanoid and cyclopoid copepod cytochrome oxidase subunit I sequence pairs are presented to demonstrate how directional selection may be linked to major shifts in adaptive zones, and that convergent evolution at the whole organism level may be the result of convergent protein adaptations.

Conclusions

Rather than replace previous methods, this new method further complements existing methods to provide a holistic glimpse of how natural selection shapes protein structure and function over evolutionary time.

Signatures of rapid evolution in urban and rural transcriptomes of white-footed mice (Peromyscus leucopus) in the New York metropolitan area

Urbanization is a major cause of ecological degradation around the world, and human settlement in large cities is accelerating. New York City (NYC) is one of the oldest and most urbanized cities in North America, but still maintains 20% vegetation cover and substantial populations of some native wildlife. The white-footed mouse, Peromyscusleucopus, is a common resident of NYC's forest fragments and an emerging model system for examining the evolutionary consequences of urbanization. In this study, we developed transcriptomic resources for urban P. leucopus to examine evolutionary changes in protein-coding regions for an exemplar "urban adapter." We used Roche 454 GS FLX+ high throughput sequencing to derive transcriptomes from multiple tissues from individuals across both urban and rural populations. From these data, we identified 31,015 SNPs and several candidate genes potentially experiencing positive selection in urban populations of P. leucopus. These candidate genes are involved in xenobiotic metabolism, innate immune response, demethylation activity, and other important biological phenomena in novel urban environments. This study is one of the first to report candidate genes exhibiting signatures of directional selection in divergent urban ecosystems.

Whole-genome sequencing identifies a recurrent functional synonymous mutation in melanoma

Synonymous mutations, which do not alter the protein sequence, have been shown to affect protein function [Sauna ZE, Kimchi-Sarfaty C (2011) Nat Rev Genet 12(10):683-691]. However, synonymous mutations are rarely investigated in the cancer genomics field. We used whole-genome and -exome sequencing to identify somatic mutations in 29 melanoma samples. Validation of one synonymous somatic mutation in BCL2L12 in 285 samples identified 12 cases that harbored the recurrent F17F mutation. This mutation led to increased BCL2L12 mRNA and protein levels because of differential targeting of WT and mutant BCL2L12 by hsa-miR-671-5p. Protein made from mutant BCL2L12 transcript bound p53, inhibited UV-induced apoptosis more efficiently than WT BCL2L12, and reduced endogenous p53 target gene transcription. This report shows selection of a recurrent somatic synonymous mutation in cancer. Our data indicate that silent alterations have a role to play in human cancer, emphasizing the importance of their investigation in future cancer genome studies.

Genomic basis of ecological niche divergence among cryptic sister species of non-biting midges

BACKGROUND

There is a lack of understanding the evolutionary forces driving niche segregation of closely related organisms. In addition, pinpointing the genes driving ecological divergence is a key goal in molecular ecology. Here, larval transcriptome sequences obtained by next-generation-sequencing are used to address these issues in a morphologically cryptic sister species pair of non-biting midges (Chironomus riparius and C. piger).

RESULTS

More than eight thousand orthologous open reading frames were screened for interspecific divergence and intraspecific polymorphisms. Despite a small mean sequence divergence of 1.53% between the sister species, 25.1% of 18,115 observed amino acid substitutions were inferred by α statistics to be driven by positive selection. Applying McDonald-Kreitman tests to 715 alignments of gene orthologues identified eleven (1.5%) genes driven by positive selection.

CONCLUSIONS

Three candidate genes were identified as potentially responsible for the observed niche segregation concerning nitrite concentration, habitat temperature and water conductivity. Additionally, signs of positive selection in the hydrogen sulfide detoxification pathway were detected, providing a new plausible hypothesis for the species' ecological differentiation. Finally, a divergently selected, nuclear encoded mitochondrial ribosomal protein may contribute to reproductive isolation due to cytonuclear coevolution.

Phylogenomic analyses of nuclear genes reveal the evolutionary relationships within the BEP clade and the evidence of positive selection in Poaceae

BEP clade of the grass family (Poaceae) is composed of three subfamilies, i.e. Bambusoideae, Ehrhartoideae, and Pooideae. Controversies on the phylogenetic relationships among three subfamilies still persist in spite of great efforts. However, previous evidence was mainly provided from plastid genes with only a few nuclear genes utilized. Given different evolutionary histories recorded by plastid and nuclear genes, it is indispensable to uncover their relationships based on nuclear genes. Here, eleven species with whole-sequenced genome and six species with transcriptomic data were included in this study. A total of 121 one-to-one orthologous groups (OGs) were identified and phylogenetic trees were reconstructed by different tree-building methods. Genes which might have undergone positive selection and played important roles in adaptive evolution were also investigated from 314 and 173 one-to-one OGs in two bamboo species and 14 grass species, respectively. Our results support the ((B, P) E) topology with high supporting values. Besides, our findings also indicate that 24 and nine orthologs with statistically significant evidence of positive selection are mainly involved in abiotic and biotic stress response, reproduction and development, plant metabolism and enzyme etc. from two bamboo species and 14 grass species, respectively. In summary, this study demonstrates the power of phylogenomic approach to shed lights on the evolutionary relationships within the BEP clade, and offers valuable insights into adaptive evolution of the grass family.

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 ω. 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 ω 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 ∼80%) and low GC (∼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.

Reconstruction of ancestral metabolic enzymes reveals molecular mechanisms underlying evolutionary innovation through gene duplication

Gene duplications are believed to facilitate evolutionary innovation. However, the mechanisms shaping the fate of duplicated genes remain heavily debated because the molecular processes and evolutionary forces involved are difficult to reconstruct. Here, we study a large family of fungal glucosidase genes that underwent several duplication events. We reconstruct all key ancestral enzymes and show that the very first preduplication enzyme was primarily active on maltose-like substrates, with trace activity for isomaltose-like sugars. Structural analysis and activity measurements on resurrected and present-day enzymes suggest that both activities cannot be fully optimized in a single enzyme. However, gene duplications repeatedly spawned daughter genes in which mutations optimized either isomaltase or maltase activity. Interestingly, similar shifts in enzyme activity were reached multiple times via different evolutionary routes. Together, our results provide a detailed picture of the molecular mechanisms that drove divergence of these duplicated enzymes and show that whereas the classic models of dosage, sub-, and neofunctionalization are helpful to conceptualize the implications of gene duplication, the three mechanisms co-occur and intertwine.

Accumulation of slightly deleterious mutations in the mitochondrial genome: a hallmark of animal domestication

The hypothesis that domestication leads to a relaxation of purifying selection on mitochondrial (mt) genomes was tested by comparative analysis of mt genes from dog, pig, chicken, and silkworm. The three vertebrate species showed mt genome phylogenies in which domestic and wild isolates were intermingled, whereas the domestic silkworm (Bombyx mori) formed a distinct cluster nested within its closest wild relative (Bombyx mandarina). In spite of these differences in phylogenetic pattern, significantly greater proportions of nonsynonymous SNPs than of synonymous SNPs were unique to the domestic populations of all four species. Likewise, in all four species, significantly greater proportions of RNA-encoding SNPs than of synonymous SNPs were unique to the domestic populations. Thus, domestic populations were characterized by an excess of unique polymorphisms in two categories generally subject to purifying selection: nonsynonymous sites and RNA-encoding sites. Many of these unique polymorphisms thus seem likely to be slightly deleterious; the latter hypothesis was supported by the generally lower gene diversities of polymorphisms unique to domestic populations in comparison to those of polymorphisms shared by domestic and wild populations.

LEC-2, a highly variable lectin in the lichen Peltigera membranacea

Lectins are a diverse group of carbohydrate binding proteins often involved in cellular interactions. A lectin gene, lec-2, was identified in the mycobiont of the lichen Peltigera membranacea. Sequencing of lec-2 open reading frames from 21 individual samples showed an unexpectedly high level of polymorphism in the deduced protein (LEC-2), which was sorted into nine haplotypes based on amino acid sequence. Calculations showed that the rates of nonsynonymous versus synonymous nucleotide substitutions deviated significantly from the null hypothesis of neutrality, indicating strong positive selection. Molecular modeling revealed that most amino acid replacements were around the putative carbohydrate-binding pocket, indicating changes in ligand binding. Lectins have been thought to be involved in the recognition of photobiont partners in lichen symbioses, and the hypothesis that positive selection of LEC-2 is driven by variation in the Nostoc photobiont partner was tested by comparing mycobiont LEC-2 haplotypes and photobiont genotypes, as represented by the rbcLX region. It was not possible to pair up the two types of marker sequences without conflicts, suggesting that positive selection of LEC-2 was not due to variation in photobiont partners.

Molecular evolution of the genomic RNA of Apple stem grooving capillovirus

The complete genome of the German isolate AC of Apple stem grooving virus (ASGV) was sequenced. It encodes two overlapping open reading frames (ORFs), similarly to previously described ASGV isolates. Two regions of high variability were detected between the ASGV isolates, variable region 1 (V1, from amino acids (aa) 532 to 570), and variable region 2 (V2, from aa 1,583 to 1,868). The phylogenetic analysis of the V1 and V2 regions suggested that the ASGV diversity was structured by host plant species rather than geographical origin. The dN/dS ratio between nonsynonymous and synonymous nucleotide substitution rates varied greatly along the ASGV genome. Most of ORF1 showed predominant negative selection except for the two regions V1 and V2. V1 showed an elevated dN and an average dS when compared to the ORF1 background but no significant positive selection was detected. The V2 region of ORF1 showed an elevated dN and a low dS when compared to the ORF1 background with an average dN/dS ≈ 3.0 indicative of positive selection. However, the V2 area includes overlapping ORFs, making the dN/dS estimate biased. Joint estimates of the selection intensity in the different ORFs by a recent method indicated that this region of ORF1 was in fact evolving close to neutrality. This was convergent with previous results showing that introduction of stop codons in this region of ORF1 did not impair plant infection. These data suggest that the elimination of a stop codon caused the overprinting of a novel coding region over the ancestral ORF.

Synthesis of Experimental Molecular Biology and Evolutionary Biology: An Example from the World of Vision

Natural selection has played an important role in establishing various phenotypes, but the molecular mechanisms of phenotypic adaptation are not well understood. The slow progress is a consequence of mutagenesis experiments in which present-day molecules were used and of the limited scope of statistical methods used to detect adaptive evolution. To fully appreciate phenotypic adaptation, the precise roles of adaptive mutations during phenotypic evolution must be elucidated through the engineering and manipulation of ancestral phenotypes. Experimental and quantum chemical analyses of dim-light vision reveal some surprising results and provide a foundation for a productive study of the adaptive evolution of various phenotypes.

Balancing selection and genetic drift at major histocompatibility complex class II genes in isolated populations of golden snub-nosed monkey (Rhinopithecus roxellana)

BACKGROUND

Small, isolated populations often experience loss of genetic variation due to random genetic drift. Unlike neutral or nearly neutral markers (such as mitochondrial genes or microsatellites), major histocompatibility complex (MHC) genes in these populations may retain high levels of polymorphism due to balancing selection. The relative roles of balancing selection and genetic drift in either small isolated or bottlenecked populations remain controversial. In this study, we examined the mechanisms maintaining polymorphisms of MHC genes in small isolated populations of the endangered golden snub-nosed monkey (Rhinopithecus roxellana) by comparing genetic variation found in MHC and microsatellite loci. There are few studies of this kind conducted on highly endangered primate species.

RESULTS

Two MHC genes were sequenced and sixteen microsatellite loci were genotyped from samples representing three isolated populations. We isolated nine DQA1 alleles and sixteen DQB1 alleles and validated expression of the alleles. Lowest genetic variation for both MHC and microsatellites was found in the Shennongjia (SNJ) population. Historical balancing selection was revealed at both the DQA1 and DQB1 loci, as revealed by excess non-synonymous substitutions at antigen binding sites (ABS) and maximum-likelihood-based random-site models. Patterns of microsatellite variation revealed population structure. FST outlier analysis showed that population differentiation at the two MHC loci was similar to the microsatellite loci.

CONCLUSIONS

MHC genes and microsatellite loci showed the same allelic richness pattern with the lowest genetic variation occurring in SNJ, suggesting that genetic drift played a prominent role in these isolated populations. As MHC genes are subject to selective pressures, the maintenance of genetic variation is of particular interest in small, long-isolated populations. The results of this study may contribute to captive breeding and translocation programs for endangered species.

Molecular adaptation and resilience of the insect's nuclear receptor USP

Background

The maintenance of biological systems requires plasticity and robustness. The function of the ecdysone receptor, a heterodimer composed of the nuclear receptors ECR (NR1H1) and USP (NR2B4), was maintained in insects despite a dramatic divergence that occurred during the emergence of Mecopterida. This receptor is therefore a good model to study the evolution of plasticity. We tested the hypothesis that selection has shaped the Ligand-Binding Domain (LBD) of USP during evolution of Mecopterida.

Results

We isolated usp and cox1 in several species of Drosophilidae, Tenebrionidae and Blattaria and estimated non-synonymous/synonymous rate ratios using maximum-likelihood methods and codon-based substitution models. Although the usp sequences were mainly under negative selection, we detected relaxation at residues located on the surface of the LBD within Mecopterida families. Using branch-site models, we also detected changes in selective constraints along three successive branches of the Mecopterida evolution. Residues located at the bottom of the ligand-binding pocket (LBP) underwent strong positive selection during the emergence of Mecopterida. This change is correlated with the acquisition of a large LBP filled by phospholipids that probably allowed the stabilisation of the new Mecopterida structure. Later, when the two subgroups of Mecopterida (Amphiesmenoptera: Lepidoptera, Trichoptera; Antliophora: Diptera, Mecoptera, Siphonaptera) diverged, the same positions became under purifying selection. Similarly, several positions of the heterodimerisation interface experienced positive selection during the emergence of Mecopterida, rapidly followed by a phase of constrained evolution. An enlargement of the heterodimerisation surface is specific for Mecopterida and was associated with a reinforcement of the obligatory partnership between ECR and USP, at the expense of homodimerisation.

Conclusions

In order to explain the episodic mode of evolution of USP, we propose a model in which the molecular adaptation of this protein is seen as a process of resilience for the maintenance of the ecdysone receptor functionality.

Sequence conservation in the Ancylostoma secreted protein-2 of Necator americanus (Na-ASP-2) from hookworm infected individuals in Thailand

The Ancylostoma secreted protein-2 of Necator americanus (Na-ASP-2) was one of the promising vaccine candidates against the most prevalent human hookworm species as adverse vaccine reaction has compromised further human vaccine trials. To elucidate the gene structure and the extent of sequence diversity, we determined the complete nucleotide sequence of the Na-asp-2 gene of individual larvae from 32 infected subjects living in 3 different endemic areas of Thailand. Sequence analysis revealed that the gene encoding Na-ASP-2 comprised 8 exons. Of 3 nucleotide substitutions in these exons, only one causes an amino acid change from leucine to methionine. A consensus conserved GT and AG at the 5' and the 3' boundaries of each intron was observed akin to those found in other eukaryotic genes. Introns of Na-asp-2 contained 23 nucleotide substitutions and 0-18 indels. The mean number of nucleotide substitutions per site (d) in introns was not significantly different from the mean number of synonymous substitutions per synonymous site (d(S)) in exons whereas d in introns was significantly exceeded d(N) (the mean number of nonsynonymous substitutions per nonsynonymous site) in exons (p<0.05), suggesting that introns and synonymous sites in exons may evolve at a similar rate whereas functional constraints at the amino acid could limit amino acid substitutions in Na-ASP-2. A recombination site was identified in an intron near the 3' portion of the gene. The positions of introns and the intron phases in the Na-asp-2 gene comparing with those in other pathogenesis-related-1 proteins of Loa loa, Onchocerca volvulus, Heterodera glycines, Caenorhabditis elegans and human were relatively conserved, suggesting evolutionary conservation of these genes. Sequence conservation in Na-ASP-2 may not compromise further vaccine design if adverse vaccine effects could be resolved whereas microheterogeneity in introns of this locus may be useful for population genetics analysis of N. americanus.

Cancer in light of experimental evolution

Cancer initiation, progression, and the emergence of therapeutic resistance are evolutionary phenomena of clonal somatic cell populations. Studies in microbial experimental evolution and the theoretical work inspired by such studies are yielding deep insights into the evolutionary dynamics of clonal populations, yet there has been little explicit consideration of the relevance of this rapidly growing field to cancer biology. Here, we examine how the understanding of mutation, selection, and spatial structure in clonal populations that is emerging from experimental evolution may be applicable to cancer. Along the way, we discuss some significant ways in which cancer differs from the model systems used in experimental evolution. Despite these differences, we argue that enhanced prediction and control of cancer may be possible using ideas developed in the context of experimental evolution, and we point out some prospects for future research at the interface between these traditionally separate areas.

Transcriptome analysis of a North American songbird, Melospiza melodia

An effective way to understand the genomics of divergence in non-model organisms is to use the transcriptome to identify genes associated with divergence. We examine the transcriptome of the song sparrow (Melospiza melodia) and contrast it with the avian models zebra finch (Taeniopygia guttata) and chicken (Gallus gallus). We aimed to (i) obtain a functional annotation of a substantial portion of the song sparrow transcriptome; (ii) compare transcript divergence; (iii) efficiently characterize single nucleotide polymorphism/indel markers possibly fixed between song sparrow subspecies; and (iv) identify the most common set of transcripts in birds using the zebra finch as a reference. Using two individuals from each of three populations, whole-body mRNA was normalized and sequenced (110Mb total). The assembly yielded 38 539 contigs [N50 (the length–weighted median) = 482 bp]; 4574 were orthologous to both model genomes and 3680 are functionally annotated. This low-coverage scan of the song sparrow transcriptome revealed 29 982 SNPs/indels, 1402 fixed between populations and subspecies. Referencing zebra finch and chicken, we identified 43 and 5 fast-evolving genes, respectively. We also identified the most common set of transcripts present in birds with respect to zebra finch. This study provides new insight into songbird transcriptomes, and candidate markers identified here may help research in songbirds (oscine Passeriformes), a frequently studied group.

Conserved and quickly evolving immunome genes have different evolutionary paths

Genetic, transcript, and protein level variations have important functional and evolutionary consequences. We performed systematic data collection and analysis of copy-number variations, single-nucleotide polymorphisms, disease-causing variations, messenger RNA splicing variants, and protein posttranslational modifications for the genes and proteins essential for human immune system. Information about polymorphic and evolutionarily fixed genetic variations was used to group immunome genes to the most conserved and the most quickly changing ones under directed selection during the recent immunome evolution. Gene Ontology terms related to adaptive immunity are associated with gene groups subject to recent directing selection. In addition, several other characteristics of the immunome genes and proteins in these two categories have statistically significant differences. The presented findings question the usability of directed mouse genes as models for human diseases and conditions and shed light on the fine tuning of human immunity and its diverse functions.

Modeling HIV-1 drug resistance as episodic directional selection

The evolution of substitutions conferring drug resistance to HIV-1 is both episodic, occurring when patients are on antiretroviral therapy, and strongly directional, with site-specific resistant residues increasing in frequency over time. While methods exist to detect episodic diversifying selection and continuous directional selection, no evolutionary model combining these two properties has been proposed. We present two models of episodic directional selection (MEDS and EDEPS) which allow the a priori specification of lineages expected to have undergone directional selection. The models infer the sites and target residues that were likely subject to directional selection, using either codon or protein sequences. Compared to its null model of episodic diversifying selection, MEDS provides a superior fit to most sites known to be involved in drug resistance, and neither one test for episodic diversifying selection nor another for constant directional selection are able to detect as many true positives as MEDS and EDEPS while maintaining acceptable levels of false positives. This suggests that episodic directional selection is a better description of the process driving the evolution of drug resistance.

When exposed to treatment, HIV-1 and other rapidly evolving viruses have the capacity to acquire drug resistance mutations (DRAMs), which limit the efficacy of antivirals. There are a number of experimentally well characterized HIV-1 DRAMs, but many mutations whose roles are not fully understood have also been reported. In this manuscript we construct evolutionary models that identify the locations and targets of mutations conferring resistance to antiretrovirals from viral sequences sampled from treated and untreated individuals. While the evolution of drug resistance is a classic example of natural selection, existing analyses fail to detect the majority of DRAMs. We show that, in order to identify resistance mutations from sequence data, it is necessary to recognize that in this case natural selection is both episodic (it only operates when the virus is exposed to the drugs) and directional (only mutations to a particular amino-acid confer resistance while allowing the virus to continue replicating). The new class of models that allow for the episodic and directional nature of adaptive evolution performs very well at recovering known DRAMs, can be useful at identifying unknown resistance-associated mutations, and is generally applicable to a variety of biological scenarios where similar selective forces are at play.

Looking for Darwin in genomic sequences--validity and success of statistical methods

The use of codon substitution models to compare synonymous and nonsynonymous substitution rates is a widely used approach to detecting positive Darwinian selection affecting protein evolution. However, in several recent papers, Hughes and colleagues claim that codon-based likelihood-ratio tests (LRTs) are logically flawed as they lack prior hypotheses and fail to accommodate random fluctuations in synonymous and nonsynonymous substitutions Friedman and Hughes (2007) also used site-based LRTs to analyze 605 gene families consisting of human and mouse paralogues. They found that the outcome of the tests was largely determined by irrelevant factors such as the GC content at the third codon positions and the synonymous rate d(S), but not by the nonsynonymous rate d(N) or the d(N)/d(S) ratio, factors that should be related to selection. Here, we reanalyze those data. Contra Friedman and Hughes, we found that the test results are related to sequence length and the average d(N)/d(S) ratio. We examine the criticisms of Hughes and suggest that they are based on misunderstandings of the codon models and on statistical errors. Our analyses suggest that codon-based tests are useful tools for comparative analysis of genomic data sets.

Evolution of adaptive phenotypic traits without positive Darwinian selection

Recent evidence suggests the frequent occurrence of a simple non-Darwinian (but non-Lamarckian) model for the evolution of adaptive phenotypic traits, here entitled the plasticity-relaxation-mutation (PRM) mechanism. This mechanism involves ancestral phenotypic plasticity followed by specialization in one alternative environment and thus the permanent expression of one alternative phenotype. Once this specialization occurs, purifying selection on the molecular basis of other phenotypes is relaxed. Finally, mutations that permanently eliminate the pathways leading to alternative phenotypes can be fixed by genetic drift. Although the generality of the PRM mechanism is at present unknown, I discuss evidence for its widespread occurrence, including the prevalence of exaptations in evolution, evidence that phenotypic plasticity has preceded adaptation in a number of taxa and evidence that adaptive traits have resulted from loss of alternative developmental pathways. The PRM mechanism can easily explain cases of explosive adaptive radiation, as well as recently reported cases of apparent adaptive evolution over ecological time.

Independent evolution of functional MHC class II DRB genes in New World bat species

Genes of the major histocompatibility complex (MHC) play a pivotal role in the vertebrate immune system and are attractive markers for functional, fitness-related, genetic variation. Although bats (Chiroptera) represent the second largest mammalian order and are prone to various emerging infectious diseases, little is known about MHC evolution in bats. In the present study, we examined expressed MHC class II DRB sequences (exons 1 to 4) of New World bat species, Saccopteryx bilineata, Carollia perspicillata, Noctilio albiventris and Noctilio leporinus (only exon 2). We found a wide range of copy number variation of DRB loci with one locus detected in the genus Noctilio and up to ten functional loci observed in S. bilineata. Sequence variation between alleles of the same taxa was high with evidence for positive selection. We found statistical support for recombination or gene conversion events among sequences within the same but not between bat species. Phylogenetic relationships among DRB alleles provided strong evidence for independent evolution of the functional MHC class II DRB genes in the three investigated species, either by recent gene duplication, or homogenization of duplicated loci by frequent gene conversion events. Phylogenetic analysis of all available chiropteran DRB exon 2 sequences confirmed their monophyletic origin within families, but revealed a possible trans-species mode of evolution pattern in congeneric bat species, e.g. within the genera Noctilio and Myotis. This is the first study investigating phylogenetic relationships of MHC genes within bats and therefore contributes to a better understanding of MHC evolution in one of the most dominant mammalian order.

The nearly neutral and selection theories of molecular evolution under the fisher geometrical framework: substitution rate, population size, and complexity

The general theories of molecular evolution depend on relatively arbitrary assumptions about the relative distribution and rate of advantageous, deleterious, neutral, and nearly neutral mutations. The Fisher geometrical model (FGM) has been used to make distributions of mutations biologically interpretable. We explored an FGM-based molecular model to represent molecular evolutionary processes typically studied by nearly neutral and selection models, but in which distributions and relative rates of mutations with different selection coefficients are a consequence of biologically interpretable parameters, such as the average size of the phenotypic effect of mutations and the number of traits (complexity) of organisms. A variant of the FGM-based model that we called the static regime (SR) represents evolution as a nearly neutral process in which substitution rates are determined by a dynamic substitution process in which the population's phenotype remains around a suboptimum equilibrium fitness produced by a balance between slightly deleterious and slightly advantageous compensatory substitutions. As in previous nearly neutral models, the SR predicts a negative relationship between molecular evolutionary rate and population size; however, SR does not have the unrealistic properties of previous nearly neutral models such as the narrow window of selection strengths in which they work. In addition, the SR suggests that compensatory mutations cannot explain the high rate of fixations driven by positive selection currently found in DNA sequences, contrary to what has been previously suggested. We also developed a generalization of SR in which the optimum phenotype can change stochastically due to environmental or physiological shifts, which we called the variable regime (VR). VR models evolution as an interplay between adaptive processes and nearly neutral steady-state processes. When strong environmental fluctuations are incorporated, the process becomes a selection model in which evolutionary rate does not depend on population size, but is critically dependent on the complexity of organisms and mutation size. For SR as well as VR we found that key parameters of molecular evolution are linked by biological factors, and we showed that they cannot be fixed independently by arbitrary criteria, as has usually been assumed in previous molecular evolutionary models.

mtDNA haplotypes differ in their probability of being eliminated by a mass die-off in an abundant seabird

In this study, we take advantage of a natural experiment--a 2004 mass die-off of the Common Murre in Alaska to determine whether closely related mtDNA haplotypes differ in their probability of being eliminated during such a short term but a marked event removing hundreds of thousands of individuals. We sequenced complete mtDNA ND2 gene (1041 bp) for 168 Common Murres sampled from seven breeding colonies across Alaska before the 2004 die-off and 127 dead murres washed ashore during the die-off. We found little mtDNA variation and lack of geographical structuring among the seven Common Murre breeding colonies in Alaska. A comparison of the single-dominant mtDNA haplotype's frequency between live murres sampled on breeding colonies before the die-off (73.2%; 95% confidence interval 66.3-79.9%) and dead murres sampled during the die-off (59.1%; 95% confidence interval 50.4-67.4%; Fisher's exact P=0.01) showed that carriers of the dominant haplotype were significantly less likely to die than carriers of other haplotypes. At the same time, the ratio of non-synonymous to synonymous substitutions did not differ between live (10:35) and dead birds (18:34; Fisher's exact P=0.26), indicating that non-synonymous substitutions were as likely to be eliminated as synonymous substitutions. These results are consistent with the possibility of positive selection on the dominant mtDNA haplotype during the die-off.

Heterogeneity in genetic diversity among non-coding loci fails to fit neutral coalescent models of population history

Inferring aspects of the population histories of species using coalescent analyses of non-coding nuclear DNA has grown in popularity. These inferences, such as divergence, gene flow, and changes in population size, assume that genetic data reflect simple population histories and neutral evolutionary processes. However, violating model assumptions can result in a poor fit between empirical data and the models. We sampled 22 nuclear intron sequences from at least 19 different chromosomes (a genomic transect) to test for deviations from selective neutrality in the gadwall (Anas strepera), a Holarctic duck. Nucleotide diversity among these loci varied by nearly two orders of magnitude (from 0.0004 to 0.029), and this heterogeneity could not be explained by differences in substitution rates alone. Using two different coalescent methods to infer models of population history and then simulating neutral genetic diversity under these models, we found that the observed among-locus heterogeneity in nucleotide diversity was significantly higher than expected for these simple models. Defining more complex models of population history demonstrated that a pre-divergence bottleneck was also unlikely to explain this heterogeneity. However, both selection and interspecific hybridization could account for the heterogeneity observed among loci. Regardless of the cause of the deviation, our results illustrate that violating key assumptions of coalescent models can mislead inferences of population history.

Allelic imbalance in Drosophila hybrid heads: exons, isoforms, and evolution

Unraveling how regulatory divergence contributes to species differences and adaptation requires identifying functional variants from among millions of genetic differences. Analysis of allelic imbalance (AI) reveals functional genetic differences in cis regulation and has demonstrated differences in cis regulation within and between species. Regulatory mechanisms are often highly conserved, yet differences between species in gene expression are extensive. What evolutionary forces explain widespread divergence in cis regulation? AI was assessed in Drosophila melanogaster-Drosophila simulans hybrid female heads using RNA-seq technology. Mapping bias was virtually eliminated by using genotype-specific references. Allele representation in DNA sequencing was used as a prior in a novel Bayesian model for the estimation of AI in RNA. Cis regulatory divergence was common in the organs and tissues of the head with 41% of genes analyzed showing significant AI. Using existing population genomic data, the relationship between AI and patterns of sequence evolution was examined. Evidence of positive selection was found in 30% of cis regulatory divergent genes. Genes involved in defense, RNAi/RISC complex genes, and those that are sex regulated are enriched among adaptively evolving cis regulatory divergent genes. For genes in these groups, adaptive evolution may play a role in regulatory divergence between species. However, there is no evidence that adaptive evolution drives most of the cis regulatory divergence that is observed. The majority of genes showed patterns consistent with stabilizing selection and neutral evolutionary processes.

Evolution of viral genomes: interplay between selection, recombination, and other forces

RNA viruses evolve very rapidly, often recombine, and are subject to strong host (immune response) and anthropogenic (antiretroviral drugs) selective forces. Given their compact and extensively sequenced genomes, comparative analysis of RNA viral data can provide important insights into the molecular mechanisms of adaptation, pathogenicity, immune evasion, and drug resistance. In this chapter, we present an example-based overview of recent advances in evolutionary models and statistical approaches that enable screening viral alignments for evidence of adaptive change in the presence of recombination, detecting bursts of directional adaptive evolution associated with the phenotypic changes, and detecting of coevolving sites in viral genes.

Molecular spandrels: tests of adaptation at the genetic level

Although much progress has been made in identifying the genes (and, in rare cases, mutations) that contribute to phenotypic variation, less is known about the effects that these genes have on fitness. Nonetheless, genes are commonly labelled as 'adaptive' if an allele has been shown to affect a phenotype with known or suspected functional importance or if patterns of nucleotide variation at the locus are consistent with positive selection. In these cases, the 'adaptive' designation may be premature and may lead to incorrect conclusions about the relationships between gene function and fitness. Experiments to test targets and agents of natural selection within a genomic context are necessary for identifying the adaptive consequences of individual alleles.

CpG site degeneration triggered by the loss of functional constraint created a highly polymorphic macaque drug-metabolizing gene, CYP1A2

Background

Elucidating the pattern of evolutionary changes in drug-metabolizing genes is an important subject not only for evolutionary but for biomedical research. We investigated the pattern of divergence and polymorphisms of macaque CYP1A1 and CYP1A2 genes, which are major drug-metabolizing genes in humans. In humans, CYP1A2 is specifically expressed in livers while CYP1A1 has a wider gene expression pattern in extrahepatic tissues. In contrast, macaque CYP1A2 is expressed at a much lower level than CYP1A1 in livers. Interestingly, a previous study has shown that Macaca fascicularis CYP1A2 harbored unusually high genetic diversity within species. Genomic regions showing high genetic diversity within species is occasionally interpreted as a result of balancing selection, where natural selection maintains highly diverged alleles with different functions. Nevertheless many other forces could create such signatures.

Results

We found that the CYP1A1/2 gene copy number and orientation has been highly conserved among mammalian genomes. The signature of gene conversion between CYP1A1 and CYP1A2 was detected, but the last gene conversion event in the simian primate lineage occurred before the Catarrhini-Platyrrhini divergence. The high genetic diversity of macaque CYP1A2 therefore cannot be explained by gene conversion between CYP1A1 and CYP1A2. By surveying CYP1A2 polymorphisms in total 91 M. fascicularis and M. mulatta, we found several null alleles segregating in these species, indicating functional constraint on CYP1A2 in macaques may have weakened after the divergence between humans and macaques. We propose that the high genetic diversity in macaque CYP1A2 is partly due to the degeneration of CpG sites, which had been maintained at a high level by purifying selection, and the rapid degeneration process was initiated by the loss of functional constraint on macaque CYP1A2.

Conclusions

Our findings show that the highly polymorphic CYP1A2 gene in macaques has not been created by balancing selection but by the burst of CpG site degeneration after loss of functional constraint. Because the functional importance of CYP1A1/2 genes is different between humans and macaques, we have to be cautious in extrapolating a drug-testing data using substrates metabolized by CYP1A genes from macaques to humans, despite of their somewhat overlapping substrate specificity.

dN/dS-based methods detect positive selection linked to trade-offs between different fitness traits in the coat protein of potato virus Y

The dN/dS ratio between nonsynonymous and synonymous substitution rates has been used extensively to identify codon positions involved in adaptive processes. However, the accuracy of this approach has been questioned, and very few studies have attempted to validate experimentally its predictions. Using the coat protein (CP) of Potato virus Y (PVY; genus Potyvirus, family Potyviridae) as a case study, we identified several candidate positively selected codon positions that differed between clades. In the CP of the N clade of PVY, positive selection was detected at codon positions 25 and 68 by both the softwares PAML and HyPhy. We introduced nonsynonymous substitutions at these positions in an infectious cDNA clone of PVY and measured the effect of these mutations on virus accumulation in its two major cultivated hosts, tobacco and potato, and on its efficiency of transmission from plant to plant by aphid vectors. The mutation at codon position 25 significantly modified the virus accumulation in the two hosts, whereas the mutation at codon position 68 significantly modified the virus accumulation in one of its hosts and its transmissibility by aphids. Both mutations were involved in adaptive trade-offs. We suggest that our study was particularly favorable to the detection of adaptive mutations using dN/dS estimates because, as obligate parasites, viruses undergo a continuous and dynamic interaction with their hosts that favors the recurrent selection of adaptive mutations and because trade-offs between different fitness traits impede (or at least slow down) the fixation of these mutations and maintain polymorphism within populations.

Evolutionary genetics of an S-like polymorphism in Papaveraceae with putative function in self-incompatibility

BACKGROUND

Papaver rhoeas possesses a gametophytic self-incompatibility (SI) system not homologous to any other SI mechanism characterized at the molecular level. Four previously published full length stigmatic S-alleles from the genus Papaver exhibited remarkable sequence divergence, but these studies failed to amplify additional S-alleles despite crossing evidence for more than 60 S-alleles in Papaver rhoeas alone.

METHODOLOGY/PRINCIPAL FINDINGS

Using RT-PCR we identified 87 unique putative stigmatic S-allele sequences from the Papaveraceae Argemone munita, Papaver mcconnellii, P. nudicuale, Platystemon californicus and Romneya coulteri. Hand pollinations among two full-sib families of both A. munita and P. californicus indicate a strong correlation between the putative S-genotype and observed incompatibility phenotype. However, we also found more than two S-like sequences in some individuals of A. munita and P. californicus, with two products co-segregating in both full-sib families of P. californicus. Pairwise sequence divergence estimates within and among taxa show Papaver stigmatic S-alleles to be the most variable with lower divergence among putative S-alleles from other Papaveraceae. Genealogical analysis indicates little shared ancestral polymorphism among S-like sequences from different genera. Lack of shared ancestral polymorphism could be due to long divergence times among genera studied, reduced levels of balancing selection if some or all S-like sequences do not function in incompatibility, population bottlenecks, or different levels of recombination among taxa. Preliminary estimates of positive selection find many sites under selective constraint with a few undergoing positive selection, suggesting that self-recognition may depend on amino acid substitutions at only a few sites.

CONCLUSIONS/SIGNIFICANCE

Because of the strong correlation between genotype and SI phenotype, sequences reported here represent either functional stylar S-alleles, tightly linked paralogs of the S-locus or a combination of both. The considerable complexity revealed in this study shows we have much to learn about the evolutionary dynamics of self-incompatibility systems.

Association of Mc1r variants with ecologically relevant phenotypes in the European ocellated lizard, Lacerta lepida

A comprehensive knowledge on the genetic basis of coloration is crucial to understand how new colour phenotypes arise and how they contribute to the emergence of new species. Variation in melanocortin-1 receptor (Mc1r), a gene that has been reported as a target for repeated evolution in a wide range of vertebrate taxa, was assessed in European ocellated lizards (Lacerta lepida) to search for associations with melanin-based colour phenotypes. Lacerta lepida subspecies' distribution is associated with the three major bio-climatic regions in the Iberian Peninsula. A nonconserved and derived substitution (T162I) was associated with the L. l. nevadensis phenotype (prevalence of brown scales). Another substitution (S172C) was associated with the presence of black scales in both L. l. lepida and L. l. iberica, but no mutations were found to be associated with the higher proportion of black in L. l. iberica. Extensive genotyping of Mc1r along the contact zone between L. l. nevadensis and L. l. lepida revealed low gene flow (only two hybrids detected). The implications of these findings are discussed in the context of previous knowledge about the evolutionary history of ocellated lizards.

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.

Adaptive evolution of the venom-targeted vWF protein in opossums that eat pitvipers

The rapid evolution of venom toxin genes is often explained as the result of a biochemical arms race between venomous animals and their prey. However, it is not clear that an arms race analogy is appropriate in this context because there is no published evidence for rapid evolution in genes that might confer toxin resistance among routinely envenomed species. Here we report such evidence from an unusual predator-prey relationship between opossums (Marsupialia: Didelphidae) and pitvipers (Serpentes: Crotalinae). In particular, we found high ratios of replacement to silent substitutions in the gene encoding von Willebrand Factor (vWF), a venom-targeted hemostatic blood protein, in a clade of opossums known to eat pitvipers and to be resistant to their hemorrhagic venom. Observed amino-acid substitutions in venom-resistant opossums include changes in net charge and hydrophobicity that are hypothesized to weaken the bond between vWF and one of its toxic snake-venom ligands, the C-type lectin-like protein botrocetin. Our results provide the first example of rapid adaptive evolution in any venom-targeted molecule, and they support the notion that an evolutionary arms race might be driving the rapid evolution of snake venoms. However, in the arms race implied by our results, venomous snakes are prey, and their venom has a correspondingly defensive function in addition to its usual trophic role.

Identification of fast-evolving genes in the scleractinian coral Acropora using comparative EST analysis

To identify fast-evolving genes in reef-building corals, we performed direct comparative sequence analysis with expressed sequence tag (EST) datasets from two acroporid species: Acropora palmata from the Caribbean Sea and A. millepora from the Great Barrier Reef in Australia. Comparison of 589 independent sequences from 1,421 A. palmata contigs, with 10,247 A. millepora contigs resulted in the identification of 196 putative homologues. Most of the homologous pairs demonstrated high amino acid similarities (over 90%). Comparisons of putative homologues showing low amino acid similarities (under 90%) among the Acropora species to the near complete datasets from two other cnidarians (Hydra magnipapillata and Nematostella vectensis) implied that some were non-orthologous. Within 86 homologous pairs, 39 exhibited dN/dS ratios significantly less than 1, suggesting that these genes are under purifying selection associated with functional constraints. Eight independent genes showed dN/dS ratios exceeding 1, while three deviated significantly from 1, suggesting that these genes may play important roles in the adaptive evolution of Acropora. Our results also indicated that CEL-III lectin was under positive selection, consistent with a possible role in immunity or symbiont recognition. Further studies are needed to clarify the possible functions of the genes under positive selection to provide insight into the evolutionary process of corals.

A case-by-case evolutionary analysis of four imprinted retrogenes

Retroposition is a widespread phenomenon resulting in the generation of new genes that are initially related to a parent gene via very high coding sequence similarity. We examine the evolutionary fate of four retrogenes generated by such an event; mouse Inpp5f_v2, Mcts2, Nap1l5, and U2af1-rs1. These genes are all subject to the epigenetic phenomenon of parental imprinting. We first provide new data on the age of these retrogene insertions. Using codon-based models of sequence evolution, we show these retrogenes have diverse evolutionary trajectories, including divergence from the parent coding sequence under positive selection pressure, purifying selection pressure maintaining parent-retrogene similarity, and neutral evolution. Examination of the expression pattern of retrogenes shows an atypical, broad pattern across multiple tissues. Protein 3D structure modeling reveals that a positively selected residue in U2af1-rs1, not shared by its parent, may influence protein conformation. Our case-by-case analysis of the evolution of four imprinted retrogenes reveals that this interesting class of imprinted genes, while similar in regulation and sequence characteristics, follow very varied evolutionary paths.

Toll-like receptors (TLRs) and mannan-binding lectin (MBL): on constant alert in a hostile environment

In the beginning were neither B cells nor T cells nor antibodies, but innate immune defense alone. The primary functional theme of innate immunity is the distinction between self and non-self, which is maintained by a vast number of cellular and subcellular components. In this context, the immense importance of the Toll-like receptors (TLRs) is well established. Positive (Darwinian) selection seems to be acting on the ligand-binding domains of these molecules, suggesting a selection pattern similar to that previously observed in the MHC proteins. In sharp contrast to TLRs, the biological significance of mannan-binding lectin (MBL) is controversial, and, concerning humans, it has been suggested that low concentration of MBL in serum represents a selective advantage. In this mini-review, based on a doctoral thesis, evolutionary aspects of TLRs and MBL are discussed.