Background: Natural selection is a potent evolutionary force that shapes phenotypic variation to match ecological conditions. However, we know little about the year-to-year consistency of selection, or how inter-annual variation in ecology shapes adaptive landscapes and ultimately adaptive radiations. Methods: Here we combine remote sensing data, field experiments, and a four-year study of natural selection to show that changes in vegetation structure associated with a severe drought altered both habitat use and natural selection in the brown anole, Anolis sagrei. Results: In natural populations, lizards increased their use of vegetation in wet years and this was correlated with selection on limb length but not body size. By contrast, a die-back of vegetation caused by drought was followed by reduced arboreality, selection on body size, and relaxed selection on limb length. With the return of the rains and recovery of vegetation, selection reverted back to pre-drought pattern of selection acting on limb length but not body size. To test for the impact of vegetation loss on natural selection during the drought, we experimentally removed vegetation on a separate study island in a naturally wet year. The experiment revealed similar inter-annual changes in selection on body size but not limb length. Conclusion: Our results illustrate the dynamic nature of ecology driving natural selection on Anolis morphology and emphasize the importance of inter-annual environmental variation in shaping adaptive variation. In addition, results illustrate the utility of using remote sensing data to examine ecology's role in driving natural selection.

Background: The multiplicity or loss of the vitellogenin (vtg) gene family in vertebrates has been argued to have broad implications for the mode of reproduction (placental or non-placental), cleavage pattern (meroblastic or holoblastic) and character of the egg (pelagic or benthic). Earlier proposals for the existence of three forms of vertebrate vtgs present conflicting models for their origin and subsequent duplication. Results: By integrating phylogenetics of novel vtg transcripts from old and modern teleosts with syntenic analyses of all available genomic variants of non-metatherian vertebrates we identify the gene orthologies between the Sarcopterygii (tetrapod branch) and Actinopterygii (fish branch). We argue that the vertebrate vtg gene cluster originated in proto-chromosome m, but that vtg genes have subsequently duplicated and rearranged following whole genome duplications. Sequencing of a novel fourth vtg transcript in labrid species, and the presence of duplicated paralogs in certain model organisms supports the notion that lineage-specific gene duplications frequently occur in teleosts. The data show that the vtg gene cluster is more conserved between acanthomorph teleosts and tetrapods, than in ostariophysan teleosts such as the zebrafish. The differential expression of the labrid vtg genes are further consistent with the notion that neofunctionalized Aa-type vtgs are important determinants of the pelagic or benthic character of the eggs in acanthomorph teleosts. Conclusion: The vertebrate vtg gene cluster existed prior to the separation of Sarcopterygii from Actinopterygii >450 million years ago, a period associated with the second round of whole genome duplication. The presence of higher copy numbers in a more highly expressed subcluster is particularly prevalent in teleosts. The differential expression and latent neofunctionalization of vtg genes in acanthomorph teleosts is an adaptive feature associated with oocyte hydration and spawning in the marine environment.

Background: Section Calochroi is one of the most species-rich lineages in the genus Cortinarius (Agaricales, Basidiomycota) and is widely distributed across boreo-nemoral areas, with some extensions into meridional zones. Previous phylogenetic studies of Calochroi (incl. section Fulvi) have been geographically restricted; therefore, phylogenetic and biogeographic relationships within this lineage at a global scale have been largely unknown. In this study, we obtained DNA sequences from a broad taxon sampling from Europe, Central America and North America. We inferred intra- and interspecific phylogenetic relationships as well as major morphological evolutionary trends within section Calochroi based on 576 ITS sequences, 230 ITS + 5.8S + D1/D2 sequences, and a combined dataset of ITS + 5.8S + D1/D2 and RPB1 sequences of a representative subsampling of 58 species. Results: More than 100 species were identified by integrating DNA sequences with morphological, macrochemical and ecological data. We show that the section Calochroi is consistently resolved with relatively high nodal support into at least seven major lineages: Calochroi, Caroviolacei, Dibaphi, Elegantiores, Napi, Pseudoglaucopodes and Splendentes; whereas Rufoolivacei and Sulfurini appeared polyphyletic. A close relationship between Dibaphi, Elegantiores, Napi and Splendentes was consistently supported. Combinations of specific morphological, pigmentation and molecular characters appear useful in circumscribing clades. Conclusions: Our analyses demonstrate that Calochroi is an exclusively northern hemispheric lineage, where species follow their host trees throughout their natural ranges within and across continents. Results of this study have contributed substantially to defining European species in this group and will help to either identify or to name new species occurring across the northern hemisphere. Major groupings are in partial agreement with earlier morphology-based and molecular phylogenetic hypotheses, but some relationships were unexpected, based on external morphology. In such cases, their true affinities appear to have been obscured by the repeated appearance of similar features among distantly related species. Therefore, further taxonomic studies are needed to evaluate the consistency of species concepts and interpretations of morphological features in a more global context. Ancestral analyses indicate that there are two major evolutionary trends within section Calochroi: (1) the development of bright pigments evolved independently multiple times and (2) the evolution of abruptly marginate to flattened stipe bulbs represents an autapomorphy of the Calochroi clade.

Background: The metazoan taxon Syndermata comprising Rotifera (in the classical sense of Monogononta+Bdelloidea+Seisonidea) and Acanthocephala has raised several theories connected to the phylogeny of these animal groups and the included subtaxa. While the monophyletic origin of Syndermata as well as of Acanthocephala is well established based on morphological and molecular data, the questions for a monophyletic origin of Monogononta, Bdelloidea, and Seisonidea and for the acanthocephalan sister group remain to be solved. Comparing the alternative hypotheses reveals that testing the phylogenetic validity of Eurotatoria (Monogononta+Bdelloidea) is the key to unravel the phylogenetic relations within Syndermata as a prerequisite for reconstructing the evolution of the acanthocephalan endoparasitism. Results: Here we present our results from a phylogenomic approach testing i) the monophyletic origin of monogononts and bdelloids and ii) their phylogenetic relations to acanthocephalans. For this analysis we have generated EST libraries of Pomphorhynchus laevis, Echinorhynchus truttae (both Palaeacanthocephala) and Brachionus plicatilis (Monogononta). By extending these data with database entries of Brachionus plicatilis (Monogononta), Philodina roseola (Bdelloidea) and 25 additional metazoan species, we conducted phylogenetic reconstructions based on 79 ribosomal proteins using maximum likelihood and bayesian approaches. Our findings strongly suggest that Eurotatoria are not monophyletic and that bdelloids are closer related to acanthocephalans than are monogononts. Conclusions: Mapping the morphological character evolution onto the molecular phylogeny suggests the emergence of a retractable anterior end (rostrum, proboscis) before the separation of Acanthocephala. This might have represented a key event leading to the later evolution of the acanthocephalan endoparasitism, given the enormous relevance of the proboscis for the anchoring of adult acanthocephalans to the definitive hosts' intestinal wall. If one further combines the present findings with previous support for Pararotatoria (Seisonidea+Acanthocephala) the rotatory organ might have undergone reduction before the separation of Acanthocephala, too.

Background: Differences in plant annual/perennial habit are hypothesized to cause a generation time effect on divergence rates. Previous studies that compared rates of divergence for internal transcribed spacer (ITS1 and ITS2) sequences of nuclear ribosomal DNA (nrDNA) in angiosperms have reached contradictory conclusions about whether differences in generation times (or other life history features) are associated with divergence rate heterogeneity. We compared annual/perennial ITS divergence rates using published sequence data, employing sampling criteria to control for possible artifacts that might obscure any actual rate variation caused by annual/perennial differences. Results: Relative rate tests employing ITS sequences from 16 phylogenetically-independent annual/perennial species pairs rejected rate homogeneity in only a few comparisons, with annuals more frequently exhibiting faster substitution rates. Treating branch length differences categorically (annual faster or perennial faster regardless of magnitude) with a sign test often indicated an excess of annuals with faster substitution rates. Annuals showed an approximately 1.6-fold rate acceleration in nucleotide substitution models for ITS. Relative rates of three nuclear loci and two chloroplast regions for the annual Arabidopsis thaliana compared with two closely related Arabidopsis perennials indicated that divergence was faster for the annual. In contrast, A. thaliana ITS divergence rates were sometimes faster and sometimes slower than the perennial. In simulations, divergence rate differences of at least 3.5-fold were required to reject rate constancy in >80 % of replicates using a nucleotide substitution model observed for the combination of ITS1 and ITS2. Simulations also showed that categorical treatment of branch length differences detected rate heterogeneity >80% of the time with a 1.5-fold or greater rate difference. Conclusions: Although rate homogeneity was not rejected in many comparisons, in cases of significant rate heterogeneity annuals frequently exhibited faster substitution rates. Our results suggest that annual taxa may exhibit a less than 2-fold rate acceleration at ITS. Since the rate difference is small and ITS lacks statistical power to reject rate homogeneity, further studies with greater power will be required to adequately test the hypothesis that annual and perennial plants have heterogeneous substitution rates. Arabidopsis sequence data suggest that relative rate tests based on multiple loci may be able to distinguish a weak acceleration in annual plants. The failure to detect rate heterogeneity with ITS in past studies may be largely a product of low statistical power.

Background: The filamentous fungus Ashbya gossypii grows into a multicellular mycelium that is distinct from the unicellular morphology of its closely related yeast species. It has been proposed that genes important for cell cycle regulation play central roles for such phenotypic differences. Because A. gossypii shares an almost identical set of cell cycle genes with the typical yeast Saccharomyces cerevisiae, the differences might occur at the level of orthologous gene regulation. Codon usage patterns were compared to identify orthologous genes with different gene regulation between A. gossypii and nine closely related yeast species. Results: Here we identified 3,151 orthologous genes between A. gossypii and nine yeast species. Two groups of genes with significant differences in codon usage (gene translation efficiency) were identified between A. gossypii and yeasts. 333 genes (Group I) and 552 genes (Group II) have significantly higher translation efficiency in A. gossypii and yeasts, respectively. Functional enrichment and pathway analysis show that Group I genes are significantly enriched with cell cycle functions whereas Group II genes are biased toward metabolic functions. Conclusions: Because translation efficiency of a gene is closely related to its functional importance, the observed functional distributions of orthologous genes with different translation efficiency might account for phenotypic differentiation between A. gossypii and yeast species. The results shed light on the mechanisms for pseudohyphal growth in pathogenic yeast species.

Background: The high polymorphism rate in the human ABO blood group gene seems to be related to susceptibility to different pathogens. It has been estimated that all genetic variation underlying the human ABO alleles appeared along the human lineage, after the divergence from the chimpanzee lineage. A paleogenetic analysis of the ABO blood group gene in Neandertals allows us to directly test for the presence of the ABO alleles in these extinct humans. Results: We have analysed two male Neandertals that were retrieved under controlled conditions at the El Sidron site in Asturias (Spain) and that appeared to be almost free of modern human DNA contamination. We find a human specific diagnostic deletion for blood group O (O01 haplotype) in both Neandertal individuals. Conclusions: These results suggest that the genetic change responsible for the O blood group in humans predates the human and Neandertal divergence. A potential selective event associated with the emergence of the O allele may have therefore occurred after humans separated from their common ancestor with chimpanzees and before the human-Neandertal population divergence.

Background: Genetic breaks separating regional lineages of marine organisms with potentially high broadcasting abilities are generally attributed either to dispersal barriers such as currents or upwelling, or to behavioural strategies promoting self-recruitment. We investigated whether such patterns could potentially also be explained by adaptations to different environmental conditions by studying two morphologically distinguishable genetic lineages of the estuarine mudprawn Upogebia africana across a biogeographic disjunction in south-eastern Africa. The study area encompasses a transition between temperate and subtropical biotas, where the warm, southward-flowing Agulhas Current is deflected away from the coast, and its inshore edge is characterised by intermittent upwelling. To determine how this phylogeographic break is maintained, we estimated gene flow among populations in the region, tested for isolation by distance as an indication of larval retention, and reared larvae of the temperate and subtropical lineages at a range of different temperatures. Results: Of four populations sampled, the two northernmost exclusively included the subtropical lineage, a central population had a mixture of both lineages, and the southernmost estuary had only haplotypes of the temperate lineage. No evidence was found for isolation by distance, and gene flow was bidirectional and of similar magnitude among adjacent populations. In both lineages, the optimum temperature for larval development was at about 23oC, but a clear difference was found at lower temperatures. While larvae of the temperate lineage could complete development at temperatures as low as 12oC, those of the subtropical lineage did not complete development below 17oC. Conclusions: The results indicate that both southward dispersal of the subtropical lineage inshore of the Agulhas Current, and its establishment in the temperate province, may be limited primarily by low water temperatures. There is no evidence that the larvae of the temperate lineage would survive less well in the subtropical province than in their native habitat, and their exclusion from this region may be due to a combination of upwelling, short larval duration with limited dispersal potential near the coast, plus transport away from the coast of larvae that become entrained in the Agulhas Current. This study shows how methods from different fields of research (genetics, physiology, oceanography and morphology) can be combined to study phylogeographic patterns.

Background: Mesoamerica represents a complex biogeographical region in the world, mostly due to its geological history. This complexity has lead to interesting biogeographical patterns that have shaped the actual diversity and distribution of the fauna in the region. The fish genus Astyanax represents an accurate model to access to these patterns due to its high phenotypic plasticity and because is the most diverse and is the widest distributed freshwater fish in America. In our study we used mitochondrial and nuclear DNA data to establish the phylogenetic relationships within the genus in Mesoamerica and develop historical biogeographical hypothesis to explain its current distribution. Results: Two independent analyses (Cyt b, 1140 bp and including 208 individuals from 147 localities and a combined matrix with Cyt b, 16S, COI and RAG1, for a subset of data) yielded similar topologies, recovering six major groups with phylogeographic structure. North America and Upper Central America were monophyletic, while Middle Central America represented a rapid radiation where relationships were not fully resolved. Lower Central America lineages showed a fragmented structure, with a very restricted distribution and high molecular divergence. All Bramocharax samples were included in different Astyanax lineages, with less than 1% of divergence with sympatric or allopatric Astyanax, giving additional evidence of the homoplasic nature of trophic specializations, which have evolved during independent evolutionary episodes and in different places. We observed a higher taxonomic diversity compared to previous phylogenetic studies. Mesoamerica colonization Astyanax was previous to the final closure of the Isthmus of Panama (6-8 Mya), leading to a fragmented configuration of Lower Central America. Additionally in upper part of Mesoamerica we observed two independent colonization tracks showing a turnover lineage process. Conclusions: Our results support a multiple and independent origin of morphological traits in this genus, where the Bramocharax genus morphotype represents a recurrent trophic adaptations of Astyanax. The presence of Astyanax in Mesoamerica was during the Miocene (~8 Mya) suggesting an incipient land bridge connecting South America and Central America, previous to the final closure of the Panama Isthmus (~3.3 Mya).

Background: Present day distributions of Palearctic taxa in northern latitudes mainly result from populations having survived in local patches during the Late Pleistocene and/or from recolonizing populations from southern temperate refugia. If well-studied Mediterranean and eastern European refugia are widely accepted, some recent biogeographical assumptions still remain unclear, such as the occurrence of multiple glacial refugia in Iberia and cryptic refugia in northern Europe during the last glaciations. The Lusitanian snail Elona quimperiana has a remarkably disjunct distribution, limited to northwestern France (Brittany), northwestern Spain and the Basque Country. By describing the phylogeographical structure of this species across its entire range, the present study attempts to identify refugia and subsequent recolonization routes. Results: Results based on 16S and COI gene sequences showed that the low genetic diversity observed in the Brittany populations should be associated with a recent demographic expansion. By contrast, populations from Spain exhibit several differentiated lineages and are characterized by demographic equilibrium, while the Basque populations are the only ones harboring typical distinct haplotypes. The center of the star-like networks of both gene sequences is occupied by a common ancestral-like haplotype found in Brittany and Spain, which might have originated from the middle of Northern Spain (i.e. Asturias, eastern Lugo and western Cantabria). Estimates of the divergence time between the Spain-Brittany and Basque lineages strongly suggest that E. quimperiana survived the Pleistocene glaciations in distinct refugia on the Iberian Peninsula, one of which is situated in Picos de Europa, and the other in the Basque Country. The occurrence of a northern refugium in France cannot be rejected as of yet. Conclusions: Present results confirm the Iberian origin of the land snail E. quimperiana and strongly support the emerging phylogeographic hypothesis of multiple refugia in Iberia during the last glaciations. The scenario of a spatial expansion of E. quimperiana from an Iberian refuge located in Asturias to northern areas provides the most probable explanation for the present distribution of this land snail. By harboring distinct haplotypes, the Basque Country populations appear to be of great importance in terms of potential adaptation, long term persistence and hence, the conservation of E. quimperiana.

Background: Animals can gain protection against predators and parasites by living in groups. The encounter-dilution effect provides protection when the probability of detection of a group does not increase in proportion to group size (i.e. encounter effect), so that predators do not offset the encounter effect by attacking more members of the group (i.e. dilution effect). In this paper, we propose a novel mechanism by which prey insects could gain by producing decoys that act as multiple targets for predators or parasitoids if these decoys are recognised as preys or hosts and negatively affect the patch foraging strategy of these predators and parasitoids. Such a decoy mechanism could be present in aphid colonies in which aphid exuviae are recognised and attacked by Aphidiine wasps. Results: We conducted a behavioural study to evaluate the effect of exuviae on parasitoid patch residence time and egg allocation in experimental aphid patches with or without exuviae. We showed that exuviae are recognised and attacked at the same level as aphids when both are present in the patch. While parasitism rate was not significantly lower in patches with exuviae when the parasitoid left the patch, the time wasted by parasitoids to handle exuviae did influence the patch residence time. As a consequence, the attack rate on the live aphids was lower in patches that contain exuviae. Conclusions: Aphids had more time available to flee and thus each individual might gain protection against parasitoids by leaving their exuviae near and within the colony. These results demonstrate that the encounter-dilution effect provided by living in a group can be enhanced by extra-materials that acts as decoy for natural enemies.

Background: Thanks to recent high coverage mass-spectrometry studies and reconstructed protein complexes, we are now in an unprecedented position to study the evolution of biological systems. Gene duplications, known to be a major source of innovation in evolution, can now be readily examined in the context of protein complexes. Results: We observe that paralogs operating in the same complex fulfill different roles: mRNA dosage increase for more than a hundred cytosolic ribosomal proteins, mutually exclusive participation of at least 54 paralogs resulting in alternative forms of complexes, and 24 proteins contributing to bona fide structural growth. Inspection of paralogous proteins participating in two independent complexes shows that an ancient, pre-duplication protein functioned in both multi-protein assemblies and a gene duplication event allowed the respective copies to specialize and split their roles. Conclusions: Variants with conditionally assembled, paralogous subunits likely have played a role in yeast's adaptation to anaerobic conditions. In a number of cases the gene duplication has given rise to one duplicate that is no longer part of a protein complex and shows an accelerated rate of evolution. Such genes could provide the raw material for the evolution of new functions.

Background: Based on the observation of an increased number of paralogous genes in teleost fishes compared with other vertebrates and on the conserved synteny between duplicated copies, it has been shown that a whole genome duplication (WGD) occurred during the evolution of Actinopterygian fish. Comparative phylogenetic dating of this duplication event suggests that it occurred early on, specifically in teleosts. It has been proposed that this event might have facilitated the evolutionary radiation and the phenotypic diversification of the teleost fish, notably by allowing the sub or neo functionalization of many duplicated genes. Results: In this paper, we studied in a wide range of Actinopterygians the duplication and fate of the androgen receptor (AR, NR3C4), a nuclear receptor known to play a key role in sex determination in vertebrates. The pattern of AR gene duplication is consistent with an early WGD event: it has been duplicated into two genes AR A and AR B after the split of the Acipenseriformes from the lineage leading to teleost fish but before the divergence of Osteoglossiformes. Genomic and syntenic analyses in addition to lack of PCR amplification show that one of the duplicated copies, AR B, was lost in several basal Clupeocephala such as Cypriniformes (including the model species zebrafish), Siluriformes, Characiformes and Salmoniformes. Interestingly, we also found that, in basal teleost fish (Osteoglossiformes and Anguilliformes), the two copies remain very similar, whereas, specifically in Percomorphs, one of the copies, AR B, has accumulated substitutions in both the ligand binding domain (LBD) and the DNA binding domain (DBD). Conclusions: The comparison of the mutations present in these divergent AR B with those known in human to be implicated in complete, partial or mild androgen insensitivity syndrome suggests that the existence of two distinct AR duplicates may be correlated to specific functional differences that may be connected to the well-known plasticity of sex determination in fish. This suggests that three specific events have shaped the present diversity of ARs in Actinopterygians: (i) early WGD, (ii) parallel loss of one duplicate in several lineages and (iii) putative neofunctionalization of the same duplicate in percomorphs, which occurred a long time after the WGD.

Background: Viruses that have spent most of their evolutionary time associated with a single host lineage should have sequences that reflect codivergence of virus and host. Several examples for RNA viruses of host-virus tree congruence are being challenged. DNA viruses, such as mastreviruses, are more likely than RNA viruses to have maintained a record of host lineage association. Results: The full genomes of 28 isolates of Wheat dwarf virus (WDV), a member of the Mastrevirus genus, from different regions of China were sequenced. The analysis of these 28 entire genomes and 18 entire genome sequences of cereal mastreviruses from other countries support the designation of wheat, barley and oat mastrevirus isolates as separate species. They revealed that relative divergence times for the viruses WDV, Barley dwarf virus (BDV), Oat dwarf virus (ODV) and Maize streak virus (MSV) are proportional to divergence times of their hosts, suggesting codivergence. Considerable diversity among Chinese isolates was found and was concentrated in hot spots in the Rep intron, SIR, LIR, and intron regions in WDV genomes. Two probable recombination events were detected in Chinese WDV isolates. Analysis including further Mastrevirus genomes concentrated on coding regions to avoid difficulties due to recombination and hyperdiversity. The analysis demonstrated congruence of trees in two branches of the genus, but not in the third. Assuming codivergence, an evolutionary rate of 10-8 substitutions per site per year was calculated. The low rate implies stronger constraints against change than are obtained by other methods of estimating the rate. Conclusions: We report tests of the hypothesis that mastreviruses have codiverged with their monocotyledonous hosts over 50 million years of evolution. The tests support the hypothesis for WDV, BDV and ODV, but not for MSV and other African streak viruses.

Background: Adaptive protein evolution is common in several Drosophila species investigated. Some studies point to very weak selection operating on amino-acid mutations, with average selection intensities on the order of Nes ~ 5 in D. melanogaster and D. simulans. Species with lower effective population sizes should undergo less adaptation since they generate fewer mutations and selection is ineffective on a greater proportion of beneficial mutations. Results: Here I study patterns of polymorphism and divergence at 91 X-linked loci in D. miranda, a species with a roughly 5-fold smaller effective population size than D. melanogaster. Surprisingly, I find a similar fraction of amino-acid mutations being driven to fixation by positive selection in D. miranda and D. melanogaster. Genes with higher rates of amino-acid evolution show lower levels of neutral diversity, a pattern predicted by recurrent adaptive protein evolution. I fit a hitchhiking model to patterns of polymorphism in D. miranda and D. melanogaster and estimate an order of magnitude higher selection coefficients for beneficial mutations in D. miranda. Conclusions: This analysis suggests that effective population size may not be a major determinant in rates of protein adaptation. Instead, adaptation may not be mutation-limited, or the distribution of fitness effects for beneficial mutations might differ vastly between different species or populations. Alternative explanation such as biases in estimating the fraction of beneficial mutations or slightly deleterious mutation models are also discussed.

Background: For survival, scorpions depend on a wide array of short neurotoxic polypeptides. The venoms of scorpions from the most studied group, the Buthida, are a rich source of small, 23-78 amino acid-long peptides, well packed by either three or four disulfide bridges that affect ion channel function in excitable and non-excitable cells. Results: In this work, by constructing a toxin transcripts data set from the venom gland of the scorpion Buthus occitanus israelis, we were able to follow the evolutionary path leading to mature toxin diversification and suggest a mechanism for leader peptide hyper-conservation. Toxins from each family were more closely related to one another than to toxins from other species, implying that fixation of duplicated genes followed speciation, suggesting early gene conversion events. Upon fixation, the mature toxin-coding domain was subjected to diversifying selection resulting in a significantly higher substitution rate that can be explained solely by diversifying selection. In contrast to the mature peptide, the leader peptide sequence was hyper-conserved and characterized by an atypical sub-neutral synonymous substitution rate. We interpret this as resulting from purifying selection acting on both the peptide and, as reported here for the first time, the DNA sequence, to create a toxin family-specific codon bias. Conclusions: We thus propose that scorpion toxin genes were shaped by selective forces acting at three levels, namely (1) diversifying the mature toxin, (2) conserving the leader peptide amino acid sequence and intriguingly, (3) conserving the leader DNA sequences.

Background: A fundamental challenge in evolutionary biology is to resolve the mechanisms that maintain paternity a hypervariable fitness component. Because females are often sexually promiscuous, this challenge hinges on establishing the mechanisms through which the ejaculates of different males compete for fertilisation (sperm competition). The competitive quality of an ejaculate is mediated by the relative number of live sperm and their motile performance. The differential rate at which rival ejaculates loose their fertilising efficiency over time is therefore expected to influence the outcome of sperm competition. Results: Here, we artificially inseminated into sets of replicate domestic hens, Gallus gallus domesticus, experimentally engineered heterospermic ejaculates containing a large number of low-quality sperm from one male, and a lower number of high-quality sperm from another male. Large, low-quality ejaculates fertilised the first eggs produced after insemination, but small, high-quality ejaculates prevailed in the long run despite their numerical disadvantage. Conclusions: Together, these results provide the first experimental demonstration that the relative competitive value of an ejaculate changes drastically over the time during which competing ejaculates are stored within the reproductive tract of a female, resulting in a marked temporal pattern of variation in paternity. A high level of replication makes these results robust. Our study was restricted to few males of a well characterised study population, and future work should explore the generality of these results.

Background: Widely used substitution models for proteins, such as the Jones-Taylor-Thornton (JTT) or Whelan and Goldman (WAG) models, are based on empirical amino acid interchange matrices estimated from databases of protein alignments that incorporate the average amino acid frequencies of the data set under examination (e.g JTT + F). Variation in the evolutionary process between sites is typically modeled using a rates-across-sites distribution such as the gamma (Gamma) distribution. However, sites in proteins also vary in the kinds of amino acid interchanges that are favoured, a feature that is ignored by standard empirical substitution matrices. Here we examine the degree to which the pattern of evolution at sites differs from that expected based on empirical amino acid substitution models and evaluate the impact of this phenomenon on phylogenetic estimation. Results: We analyzed 21 large protein alignments with two statistical tests designed to detect deviation of site-specific amino acid distributions from data simulated under the standard empirical substitution model: JTT+ F + Gamma. We found that the number of states at a given site is, on average, smaller and the frequencies of these states are less uniform than expected based on a JTT + F + Gamma substitution model. With a four-taxon example, we show that phylogenetic estimation under the JTT + F + Gamma model is seriously biased by a long-branch attraction artefact if the data are simulated under a model utilizing the observed site-specific amino acid frequencies from an alignment. Principal components analyses indicate the existence of at least four major site-specific frequency classes in these 21 protein alignments. Using a mixture model with these four separate classes of site-specific state frequencies plus a fifth class of global frequencies (the JTT + cF + Gamma model), significant improvements in model fit over the standard JTT + F + Gamma model for real data sets can be achieved. This simple mixture model also reduces the long-branch attraction problem, as shown by simulations and analyses of a real phylogenomic data set. Conclusions: Protein families display site-specific evolutionary dynamics, departing from standard protein phylogenetic models. Accurate estimation of protein phylogenies requires models that accommodate the heterogeneity in the evolutionary process across sites. To this end, we have implemented the class frequency mixture model proposed here in a freely-available program called QmmRAxML for phylogenetic estimation.

Background: Brain-expressed genes that were created in primate lineage represent obvious candidates to investigate molecular mechanisms that contributed to neural reorganization and emergence of new behavioural functions in Homo sapiens. PMCHL1 arose from retroposition of a pro-melanin-concentrating hormone (PMCH) antisense mRNA on the ancestral human chromosome 5p14 when platyrrhines and catarrhines diverged. Mutations before divergence of hylobatidae led to creation of new exons and finally PMCHL1 duplicated in an ancestor of hominids to generate PMCHL2 at the human chromosome 5q13. A complex pattern of spliced and unspliced PMCHL RNAs were found in human brain and testis. Results: Several novel spliced PMCHL transcripts have been characterized in human testis and fetal brain, identifying an additional exon and novel splice sites. Sequencing of PMCHL genes in several non-human primates allowed to carry out phylogenetic analyses revealing that the initial retroposition event took place within an intron of the brain cadherin (CDH12) gene, soon after platyrrhine/catarrhine divergence, i.e. 30-35 Mya, and was concomitant with the insertion of an AluSg element. Sequence analysis of the spliced PMCHL transcripts identified only short ORFs of less than 300 bp, with low (VMCH-p8 and protein variants) or no evolutionary conservation. Western blot analyses of human and macaque tissues expressing PMCHL RNA failed to reveal any protein corresponding to VMCH-p8 and protein variants encoded by spliced transcripts. Conclusions: Our present results improve our knowledge of the gene structure and the evolutionary history of the primate-specific chimeric PMCHL genes. These genes produce multiple spliced transcripts, bearing short, non-conserved and apparently non-translated ORFs that may function as mRNA-like non-coding RNAs.

Background: The braconid subfamily Rogadinae is a large, cosmopolitan group of endoparasitoid wasps characterised by 'mummifying' their lepidopteran host larvae, from which the adult subsequently emerges. Rogadines attack a variety of both macro- and microlepidopteran taxa, although the speciose genus Aleiodes almost exclusively attacks macrolepidopterans. Here, we investigate the phylogenetic history of the Rogadinae, revise their higher-level classification and assess the evolution of their host ranges and mummy types. We also assess the divergence times within the subfamily and discuss the reasons for the extraordinary evolutionary diversification of Aleiodes. Results: Our Bayesian analyses weakly support the monophyly of the subfamily. A clade comprising all Aleiodes species and some other taxa is not nested within the tribe Rogadini as previously supposed, but instead is recovered as sister to the Yeliconini, with the remaining Rogadini genera being recovered as sister to the Stiropiini. The Rogadinae is estimated to have originated during the mid to late Eocene, 36.1-51.62 MYA. Molecular dating gives a more recent origin for the Aleiodes clade (17.98-41.76 MYA) compared to the origins proposed for two of its principal lepidopteran host groups (Noctuidae: 60.7-113.4 MYA; Geometridae 48-62 MYA). The Bayesian ancestral reconstruction of the emergence habits from the mummified hosts weakly recovered an anterior emergence as the ancestral condition for the subfamily. Producing a hard mummy has evolved at various times independently, though most of the species with this biology belong to the Aleiodes clade. Conclusion: Based on our results, we erect the tribe Aleiodini nov. to include Aleiodes and Heterogamus stat. rev. Cordylorhogas, Pholichora and Hemigyroneuron are synonymised with Aleiodes. The molecular dating of clades and the ancestral reconstruction of host ranges support the hypothesis that radiation within Aleiodes s. s. was due to host recruitment leading to host range expansion followed by speciation, and not to parasitoid-host coevolution. Within the Rogadinae, variation in the site of emergence from the mummified host probably evolved as a consequence of the mummy's site and mode of formation, and the extent of mummy tanning/hardness to the degree of protection needed in relation to the cost of providing it.