Search

Genitalia are among the fastest evolving morphological traits in arthropods. Among the many hypotheses aimed at explaining this observation, some explicitly or implicitly predict concomitant male and female changes of genital traits that interact during copulation (i.e., lock and key, sexual conflict, cryptic female choice and pleiotropy). Testing these hypotheses requires insights into whether male and female copulatory structures that physically interact during mating also affect each other's evolution and patterns of diversification. Here we compare and contrast size and shape evolution of male and female structures that are known to interact tightly during copulation using two model systems: (a) the sister species O. taurus (1 native, 3 recently established populations) and O. illyricus, and (b) the species-complex O. fracticornis-similis-opacicollis. Partial Least Squares analyses indicated very little to no correlation between size and shape of copulatory structures, both in males and females. Accordingly, comparing shape and size diversification patterns of genitalia within each sex showed that the two components diversify readily - though largely independently of each other - within and between species. Similarly, comparing patterns of divergence across sexes showed that relative sizes of male and female copulatory organs diversify largely independent of each other. However, performing this analysis for genital shape revealed a signature of parallel divergence. Our results therefore suggest that male and female copulatory structures that are linked mechanically during copulation may diverge in concert with respect to their shapes. Furthermore, our results suggest that genital divergence in general, and co-divergence of male and female genital shape in particular, can evolve over an extraordinarily short time frame. Results are discussed in the framework of the hypotheses that assume or predict concomitant evolutionary changes in male and female copulatory organs.

BackgroundThe King Island Emu (Dromaius ater) of Australia is one of several extinct emu taxa whose taxonomic relationship to the modern Emu (D. novaehollandiae) is unclear. King Island Emu were mainly distinguished by their much smaller size and a reported darker colour compared to modern Emu.Methodology and resultsWe investigated the evolutionary relationships between the King Island and modern Emu by the recovery of both nuclear and mitochondrial DNA sequences from sub-fossil remains. The complete mitochondrial control (1,094 bp) and cytochrome c oxidase subunit I (COI) region (1,544 bp), as well as a region of the melanocortin 1 receptor gene (57 bp) were sequenced using a multiplex PCR approach. The results show that haplotypes for King Island Emu fall within the diversity of modern Emu.ConclusionsThese data show the close relationship of these emu when compared to other congeneric bird species and indicate that the King Island and modern Emu share a recent common ancestor. King Island emu possibly underwent insular dwarfism as a result of phenotypic plasticity. The close relationship between the King Island and the modern Emu suggests it is most appropriate that the former should be considered a subspecies of the latter. Although both taxa show a close genetic relationship they differ drastically in size. This study also suggests that rates of morphological and neutral molecular evolution are decoupled.

Background: Antagonistic species often interact via matching of phenotypes, and interactions between brood parasitic common cuckoos (Cuculus canorus) and their hosts constitute classic examples. The outcome of a parasitic event is often determined by the match between host and cuckoo eggs, giving rise to potentially strong associations between fitness and egg phenotype. Yet, empirical efforts aiming to document and understand the resulting evolutionary outcomes are in short supply. Methods/Principal Findings: We used avian color space models to analyze patterns of egg color variation within and between the cuckoo and two closely related hosts, the nomadic brambling (Fringilla montifringilla) and the site fidelic chaffinch (F. coelebs). We found that there is pronounced opportunity for disruptive selection on brambling egg coloration. The corresponding cuckoo host race has evolved egg colors that maximize fitness in both sympatric and allopatric brambling populations. By contrast, the chaffinch has a more bimodal egg color distribution consistent with the evolutionary direction predicted for the brambling. Whereas the brambling and its cuckoo host race show little geographical variation in their egg color distributions, the chaffinch's distribution becomes increasingly dissimilar to the brambling's distribution towards the core area of the brambling cuckoo host race. Conclusion: High rates of brambling gene flow is likely to cool down coevolutionary hot spots by cancelling out the selection imposed by a patchily distributed cuckoo host race, thereby promoting a matching equilibrium. By contrast, the site fidelic chaffinch is more likely to respond to selection from adapting cuckoos, resulting in a markedly more bimodal egg color distribution. The geographic variation in the chaffinch's egg color distribution could reflect a historical gradient in parasitism pressure. Finally, marked cuckoo egg polymorphisms are unlikely to evolve in these systems unless the hosts evolve even more exquisite egg recognition capabilities than currently possessed. © 2011 Vikan et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Coupled behavioural observations and acoustical recordings of aggressive dyadic contests showed that the mudskipper Periophthalmodon septemradiatus communicates acoustically while out of water. An analysis of intraspecific variability showed that specific acoustic components may act as tags for individual recognition, further supporting the sounds' communicative value. A correlative analysis amongst acoustical properties and video-acoustical recordings in slow-motion supported first hypotheses on the emission mechanism. Acoustic transmission through the wet exposed substrate was also discussed. These observations were used to support an "exaptation hypothesis", i.e. the maintenance of key adaptations during the first stages of water-to-land vertebrate eco-evolutionary transitions (based on eco-evolutionary and palaeontological considerations), through a comparative bioacoustic analysis of aquatic and semiterrestrial gobiid taxa. In fact, a remarkable similarity was found between mudskipper vocalisations and those emitted by gobioids and other soniferous benthonic fishes.

The second parity rule states that, if there is no bias in mutation or selection, then within each strand of DNA complementary bases are present at approximately equal frequencies. In bacteria, however, there is commonly an excess of G (over C) and, to a lesser extent, T (over A) in the replicatory leading strand. The low G+C Firmicutes, such as Staphylococcus aureus, are unusual in displaying an excess of A over T on the leading strand. As mutation has been established as a major force in the generation of such skews across various bacterial taxa, this anomaly has been assumed to reflect unusual mutation biases in Firmicute genomes. Here we show that this is not the case and that mutation bias does not explain the atypical AT skew seen in S. aureus. First, recently arisen intergenic SNPs predict the classical replication-derived equilibrium enrichment of T relative to A, contrary to what is observed. Second, sites predicted to be under weak purifying selection display only weak AT skew. Third, AT skew is primarily associated with largely non-synonymous first and second codon sites and is seen with respect to their sense direction, not which replicating strand they lie on. The atypical AT skew we show to be a consequence of the strong bias for genes to be co-oriented with the replicating fork, coupled with the selective avoidance of both stop codons and costly amino acids, which tend to have T-rich codons. That intergenic sequence has more A than T, while at mutational equilibrium a preponderance of T is expected, points to a possible further unresolved selective source of skew., Author Summary When considering a single strand of DNA, it is not necessarily the case that the frequency of each base should equal its complementary partner, such that A = T and G = C. For the leading strand, it is typically the case that Gs are more common than Cs, and Ts more common than As. This bias is widely thought to arise due to different mutational biases during replication. The Firmicutes exhibit an atypical preference for A over T on the leading strand, and here we show that selection, rather than mutation, can explain this exception. For those bases within coding regions, selection acts to inflate the frequency of A over T in order to avoid stop codons and to use metabolically cheap amino acids. Because genes are not orientated randomly, this manifests as an overall enrichment of A on the leading strand. Furthermore, a direct examination of mutational patterns is inconsistent with the observed enrichment of As. Curiously, our data also point to an unresolved source of selection on synonymous and intergenic sites, which are widely assumed to be neutral.

Tetrapod fossil tracks are known from the Middle Devonian (Eifelian at ca. 397 million years ago - MYA), and their earliest bony remains from the Upper Devonian (Frasnian at 375–385 MYA). Tetrapods are now generally considered to have colonized land during the Carboniferous (i.e., after 359 MYA), which is considered to be one of the major events in the history of life. Our analysis on tetrapod evolution was performed using molecular data consisting of 13 proteins from 17 species and different paleontological data. The analysis on the molecular data was performed with the program TreeSAAP and the results were analyzed to see if they had implications on the paleontological data collected. The results have shown that tetrapods evolved from marine environments during times of higher oxygen levels. The change in environmental conditions played a major role in their evolution. According to our analysis this evolution occurred at about 397–416 MYA during the Early Devonian unlike previously thought. This idea is supported by various environmental factors such as sea levels and oxygen rate, and biotic factors such as biodiversity of arthropods and coral reefs. The molecular data also strongly supports lungfish as tetrapod's closest living relative.

Closely related species that show clear phenotypic divergence, but without obvious geographic barriers, can provide opportunities to study how diversification can occur when opportunities for allopatric speciation are limited. We examined genetic divergence in the coral reef fish genus Hypoplectrus (family: Serranidae), which comprises of 10-14 morphotypes that are distinguished solely by their distinct colour patterns, but which show little genetic differentiation. Our goal was to detect loci that show clear disequilibrium between morphotypes and across geographical locations. We conducted Amplified Fragment Length Polymorphism molecular analysis to quantify genetic differentiation among, and selection between, morphotypes. Three loci were consistently divergent beyond neutral expectations in repeated pair-wise morphotype comparisons using two different methods. These loci provide the first evidence for genes that may be associated with colour morphotype in the genus Hypoplectrus., This work was funded by the National Environment Research Council (studentship S NER/SJ2004/13064 and LSMSF application no. EK76-02/05), with the additional funding by SMGF and the Sydney L. Wright Fellowship (contribution number 1700).

Most viruses express during infection products that prevent or neutralize the effect of the host dsRNA activated protein kinase (PKR). Translation of Sindbis virus (SINV) mRNA escapes to PKR activation and eIF2 phosphorylation in infected cells by a mechanism that requires a stem loop structure in viral 26S mRNA termed DLP to initiate translation in the absence of functional eIF2. Unlike the rest of viruses tested, we found that Alphavirus infection allowed a strong PKR activation and eIF2α phosphorylation in vitro and in infected animals so that the presence of DLP structure in mRNA was critical for translation and replication of SINV. Interestingly, infection of MEFs with some viruses that express PKR inhibitors prevented eIF2α phosphorylation after superinfection with SINV, suggesting that viral anti-PKR mechanisms could be exchangeable. Thus, translation of SINV mutant lacking the DLP structure (ΔDLP) in 26S mRNA was partially rescued in cells expressing vaccinia virus (VV) E3 protein, a known inhibitor of PKR. This case of heterotypic complementation among evolutionary distant viruses confirmed experimentally a remarkable case of convergent evolution in viral anti-PKR mechanisms. Our data reinforce the critical role of PKR in regulating virus-host interaction and reveal the versatility of viruses to find different solutions to solve the same conflict., This work was supported in part from the VIRUS-HOST interaction programme (Comunidad de Madrid) and by grants from the Ministerio de Ciencia e Innovación (SAF2006-09810; SAF2008-02036) and the Fundación Mutua Madrileña (FMM 2008). The institutional support from Fundación Ramón Areces is also acknowledged. E.D was a recipient of VIRUS-HOST programme postdoctoral contract. R.T. was a recipient of the SAF2006-09810 contract and I.V. was a researcher of Ramón y Cajal Programme.

Leaf-cutter ants are one of the most important herbivorous insects in the Neotropics, harvesting vast quantities of fresh leaf material. The ants use leaves to cultivate a fungus that serves as the colony's primary food source. This obligate ant-fungus mutualism is one of the few occurrences of farming by non-humans and likely facilitated the formation of their massive colonies. Mature leaf-cutter ant colonies contain millions of workers ranging in size from small garden tenders to large soldiers, resulting in one of the most complex polymorphic caste systems within ants. To begin uncovering the genomic underpinnings of this system, we sequenced the genome of Atta cephalotes using 454 pyrosequencing. One prediction from this ant's lifestyle is that it has undergone genetic modifications that reflect its obligate dependence on the fungus for nutrients. Analysis of this genome sequence is consistent with this hypothesis, as we find evidence for reductions in genes related to nutrient acquisition. These include extensive reductions in serine proteases (which are likely unnecessary because proteolysis is not a primary mechanism used to process nutrients obtained from the fungus), a loss of genes involved in arginine biosynthesis (suggesting that this amino acid is obtained from the fungus), and the absence of a hexamerin (which sequesters amino acids during larval development in other insects). Following recent reports of genome sequences from other insects that engage in symbioses with beneficial microbes, the A. cephalotes genome provides new insights into the symbiotic lifestyle of this ant and advances our understanding of host–microbe symbioses., Author Summary Leaf-cutter ant workers forage for and cut leaves that they use to support the growth of a specialized fungus, which serves as the colony's primary food source. The ability of these ants to grow their own food likely facilitated their emergence as one of the most dominant herbivores in New World tropical ecosystems, where leaf-cutter ants harvest more plant biomass than any other herbivore species. These ants have also evolved one of the most complex forms of division of labor, with colonies composed of different-sized workers specialized for different tasks. To gain insight into the biology of these ants, we sequenced the first genome of a leaf-cutter ant, Atta cephalotes. Our analysis of this genome reveals characteristics reflecting the obligate nutritional dependency of these ants on their fungus. These findings represent the first genetic evidence of a reduced capacity for nutrient acquisition in leaf-cutter ants, which is likely compensated for by their fungal symbiont. These findings parallel other nutritional host–microbe symbioses, suggesting convergent genomic modifications in these types of associations.

Cytochrome c (cyt c) participates in two crucial cellular processes, energy production and apoptosis, and unsurprisingly is a highly conserved protein. However, previous studies have reported for the primate lineage (i) loss of the paralogous testis isoform, (ii) an acceleration and then a deceleration of the amino acid replacement rate of the cyt c somatic isoform, and (iii) atypical biochemical behavior of human cyt c. To gain insight into the cause of these major evolutionary events, we have retraced the history of cyt c loci among primates. For testis cyt c, all primate sequences examined carry the same nonsense mutation, which suggests that silencing occurred before the primates diversified. For somatic cyt c, maximum parsimony, maximum likelihood, and Bayesian phylogenetic analyses yielded the same tree topology. The evolutionary analyses show that a fast accumulation of non-synonymous mutations (suggesting positive selection) occurred specifically on the anthropoid lineage root and then continued in parallel on the early catarrhini and platyrrhini stems. Analysis of evolutionary changes using the 3D structure suggests they are focused on the respiratory chain rather than on apoptosis or other cyt c functions. In agreement with previous biochemical studies, our results suggest that silencing of the cyt c testis isoform could be linked with the decrease of primate reproduction rate. Finally, the evolution of cyt c in the two sister anthropoid groups leads us to propose that somatic cyt c evolution may be related both to COX evolution and to the convergent brain and body mass enlargement in these two anthropoid clades.

This study investigates the genetic diversity, population structure and demographic history of the afrotropical butterfly Bicyclus anynana using mitochondrial DNA (mtDNA). Samples from six wild populations covering most of the species range from Uganda to South Africa were compared for the cytochrome c oxidase subunit gene (COI). Molecular diversity indices show overall high mtDNA diversity for the populations, but low nucleotide divergence between haplotypes. Our results indicate relatively little geographic population structure among the southern populations, especially given the extensive distributional range and an expectation of limited gene flow between populations. We implemented neutrality tests to assess signatures of recent historical demographic events. Tajima's D test and Fu's F(S) test both suggested recent population growth for the populations. The results were only significant for the southernmost populations when applying Tajima's D, but Fu's F(S) indicated significant deviations from neutrality for all populations except the one closest to the equator. Based on our own findings and those from pollen and vegetation studies, we hypothesize that the species range of B. anynana was reduced to equatorial refugia during the last glacial period, and that the species expanded southwards during the past 10.000 years. These results provide crucial background information for studies of phenotypic and molecular adaptation in wild populations of B. anynana.

It is valuable to extend genotyping studies of Helicobacter pylori to strains from indigenous communities across the world to better define adaption, evolution, and associated diseases. We aimed to genetically characterize both human individuals and their infecting H. pylori from indigenous communities of Mexico, and to compare them with those from other human groups. We studied individuals from three indigenous groups, Tarahumaras from the North, Huichols from the West and Nahuas from the center of Mexico. Volunteers were sampled at their community site, DNA was isolated from white blood cells and mtDNA, Y-chromosome, and STR alleles were studied. H. pylori was cultured from gastric juice, and DNA extracted for genotyping of virulence and housekeeping genes. We found Amerindian mtDNA haplogroups (A, B, C, and D), Y-chromosome DYS19T, and Amerindian STRs alleles frequent in the three groups, confirming Amerindian ancestry in these Mexican groups. Concerning H.pylori cagA phylogenetic analyses, although most isolates were of the Western type, a new Amerindian cluster neither Western nor Asian, was formed by some indigenous Mexican, Colombian, Peruvian and Venezuelan isolates. Similarly, vacA phylogenetic analyses showed the existence of a novel Amerindian type in isolates from Alaska, Mexico and Colombia. With hspA strains from Mexico and other American groups clustered within the three major groups, Asian, African or European. Genotyping of housekeeping genes confirmed that Mexican strains formed a novel Asian-related Amerindian group together with strains from remote Amazon Aborigines. This study shows that Mexican indigenous people with Amerindian markers are colonized with H. pylori showing admixture of Asian, European and African strains in genes known to interact with the gastric mucosa. We present evidence of novel Amerindian cagA and vacA alleles in indigenous groups of North and South America.

Variation within a population is a key feature in evolution, because it can increase or impede response to selection, depending on whether or not the intrapopulational variance is correlated to the change under selection. Hence, main directions of genetic variance have been proposed to constitute "lines of least resistance to evolution'' along which evolution would be facilitated. Yet, the screening of selection occurs at the phenotypic level, and the phenotypic variance is not only the product of the underlying genetic variance, but also of developmental processes. It is thus a key issue for interpreting short and long term evolutionary patterns to identify whether main directions of phenotypic variance indeed constitute direction of facilitated evolution, and whether this is favored by developmental processes preferably generating certain phenotypes. We tackled these questions by a morphometric quantification of the directions of variance, compared to the direction of evolution of the first upper and lower molars of wild continental and insular house mice. The main phenotypic variance indeed appeared as channeling evolution between populations. The upper molar emerged as highly evolvable, because a strong allometric component contributed to its variance. This allometric relationship drove a repeated but independent evolution of a peculiar upper molar shape whenever size increased. This repeated evolution, together with knowledge about the molar development, suggest that the main direction of phenotypic variance correspond here to a "line of least developmental resistance'' along which evolution between population is channeled.

International audience; Members of the genus Xenorhabdus are entomopathogenic bacteria that associate with nematodes. The nematode-bacteria pair infects and kills insects, with both partners contributing to insect pathogenesis and the bacteria providing nutrition to the nematode from available insect-derived nutrients. The nematode provides the bacteria with protection from predators, access to nutrients, and a mechanism of dispersal. Members of the bacterial genus Photorhabdus also associate with nematodes to kill insects, and both genera of bacteria provide similar services to their different nematode hosts through unique physiological and metabolic mechanisms. We posited that these differences would be reflected in their respective genomes. To test this, we sequenced to completion the genomes of Xenorhabdus nematophila ATCC 19061 and Xenorhabdus bovienii SS-2004. As expected, both Xenorhabdus genomes encode many anti-insecticidal compounds, commensurate with their entomopathogenic lifestyle. Despite the similarities in lifestyle between Xenorhabdus and Photorhabdus bacteria, a comparative analysis of the Xenorhabdus, Photorhabdus luminescens, and P. asymbiotica genomes suggests genomic divergence. These findings indicate that evolutionary changes shaped by symbiotic interactions can follow different routes to achieve similar end points.

Ascosphaera apis causes chalkbrood in honeybees, a chronic disease that reduces the number of viable offspring in the nest. Although lethal for larvae, the disease normally has relatively low virulence at the colony level. A recent study showed that there is genetic variation for host susceptibility, but whether Ascosphaera apis strains differ in virulence is unknown. We exploited a recently modified in vitro rearing technique to infect honeybee larvae from three colonies with naturally mated queens under strictly controlled laboratory conditions, using four strains from two distinct A. apis clades. We found that both strain and colony of larval origin affected mortality rates. The strains from one clade caused 12-14% mortality while those from the other clade induced 71-92% mortality. Larvae from one colony showed significantly higher susceptibility to chalkbrood infection than larvae from the other two colonies, confirming the existence of genetic variation in susceptibility across colonies. Our results are consistent with antagonistic coevolution between a specialized fungal pathogen and its host, and suggest that beekeeping industries would benefit from more systematic monitoring of this chronic stress factor of their colonies.

Group B of the Sox transcription factor family is crucial in embryo development in the insects and vertebrates. Sox group B, unlike the other Sox groups, has an unusually enlarged functional repertoire in insects, but the timing and mechanism of the expansion of this group were unclear. We collected and analyzed data for Sox group B from 36 species of 12 phyla representing the major metazoan clades, with an emphasis on arthropods, to reconstruct the evolutionary history of SoxB in bilaterians and to date the expansion of Sox group B in insects. We found that the genome of the bilaterian last common ancestor probably contained one SoxB1 and one SoxB2 gene only and that tandem duplications of SoxB2 occurred before the arthropod diversification but after the arthropod-nematode divergence, resulting in the basal repertoire of Sox group B in diverse arthropod lineages. The arthropod Sox group B repertoire expanded differently from the vertebrate repertoire, which resulted from genome duplications. The parallel increases in the Sox group B repertoires of the arthropods and vertebrates are consistent with the parallel increases in the complexity and diversification of these two important organismal groups.

To elucidate the speciation mechanisms prevalent within hotspots of biodiversity, and the evolutionary processes behind the rise of their species-rich and endemic biota, we investigated the phylogeny of the giant fire-millipede genus Aphistogoniulus Silvestri, 1897, a Malagasy endemic. This study is the first comprehensive (molecular and morphological) phylogenetic study focusing on millipede (class Diplopoda) speciation on Madagascar. The morphological analysis is based on 35 morphological characters and incorporates ten described as well as two newly described species (A. rubrodorsalis n. sp. and A. jeekeli n. sp.) of Aphistogoniulus. The molecular analysis is based on both mitochondrial (COI and 16S), and nuclear genes (complete 18S rDNA), together comprised of 3031 base pairs, which were successfully sequenced for 31 individual specimens and eight species of Aphistogoniulus. In addition to the null-model (speciation by distance), two diversification models, mountain refugia and ecotone shift, were discovered to play a role in the speciation of soil arthropods on Madagascar. Mountain refugia were important in the speciation of the A. cowani clade, with three species occurring in the Andringitra and Ranomafana Mountains in the southeast (A. cowani), the Ambohijanahary and Ambohitantely Mountains in the mid-west (A. sanguineus), and the Marojejy Mountain in the northeast (A. rubrodorsalis n. sp.). An ecotone shift from the eastern rainforest to the unique subarid spiny forest of Mahavelo was discovered in the A. vampyrus - A. aridus species-pair. In the monophyletic A. diabolicus clade, evidence for divergent evolution of sexual morphology was detected: species with greatly enlarged gonopods are sister-taxa to species with normal sized gonopods. Among the large-bodied Spirobolida genera of Madagascar, Colossobolus and Sanguinobolus were found to be close sister-genera to Aphistogoniulus. Forest destruction has caused forest corridors between populations to disappear, which might limit the possible resolution of biogeographic analyses on Madagascar.

Contact: john.heraty@ucr.edu; International audience; Chalcidoidea (Hymenoptera) are extremely diverse with more than 23,000 species described and over 500,000 species estimated to exist. This is the first comprehensive phylogenetic analysis of the superfamily based on a molecular analysis of 18S and 28S ribosomal gene regions for 19 families, 72 subfamilies, 343 genera and 649 species. The 56 outgroups are comprised of Ceraphronoidea and most proctotrupomorph families, including Mymarommatidae. Data alignment and the impact of ambiguous regions are explored using a secondary structure analysis and automated (MAFFT) alignments of the core and pairing regions and regions of ambiguous alignment. Both likelihood and parsimony approaches are used to analyze the data. Overall there is no impact of alignment method, and few but substantial differences between likelihood and parsimony approaches. Monophyly of Chalcidoidea and a sister group relationship between Mymaridae and the remaining Chalcidoidea is strongly supported in all analyses. Either Mymarommatoidea or Diaprioidea are the sister group of Chalcidoidea depending on the analysis. Likelihood analyses place Rotoitidae as the sister group of the remaining Chalcidoidea after Mymaridae, whereas parsimony nests them within Chalcidoidea. Some traditional family groups are supported as monophyletic (Agaonidae, Eucharitidae, Encyrtidae, Eulophidae, Leucospidae, Mymaridae, Ormyridae, Signiphoridae, Tanaostigmatidae and Trichogrammatidae). Several other families are paraphyletic (Perilampidae) or polyphyletic (Aphelinidae, Chalcididae, Eupelmidae, Eurytomidae, Pteromalidae, Tetracampidae and Torymidae). Evolutionary scenarios discussed for Chalcidoidea include the evolution of phytophagy, egg parasitism, sternorrhynchan parasitism, hypermetamorphic development and heteronomy.

Owing to their long life span and ecological dominance in many communities, forest trees are subject to attack from a diverse array of herbivores throughout their range, and have therefore developed a large number of both constitutive and inducible defenses. We used molecular population genetics methods to examine the evolution of eight genes in European aspen, Populus tremula, that are all associated with defensive responses against pests and/or pathogens, and have earlier been shown to become strongly up-regulated in poplars as a response to wounding and insect herbivory. Our results show that the majority of these defense genes show patterns of intraspecific polymorphism and site-frequency spectra that are consistent with a neutral model of evolution. However, two of the genes, both belonging to a small gene family of polyphenol oxidases, show multiple deviations from the neutral model. The gene PPO1 has a 600 bp region with a highly elevated K(A)/K(S) ratio and reduced synonymous diversity. PPO1 also shows a skew toward intermediate frequency variants in the SFS, and a pronounced fixation of non-synonymous mutations, all pointing to the fact that PPO1 has been subjected to recurrent selective sweeps. The gene PPO2 shows a marked excess of high frequency, derived variants and shows many of the same trends as PPO1 does, even though the pattern is less pronounced, suggesting that PPO2 might have been the target of a recent selective sweep. Our results supports data from both Populus and other species which have found that the the majority of defense-associated genes show few signs of selection but that a number of genes involved in mediating defense against herbivores show signs of adaptive evolution.

Herpes simplex virus type 1 (HSV-1) and varicella zoster virus (VZV) are closely related viruses causing lifelong infections. They are typically associated with mucocutaneous or skin lesions, but may also cause severe neurological or ophthalmic diseases, possibly due to viral- and/or host-genetic factors. Although these viruses are well characterized, genome-wide evolutionary studies have hitherto only been presented for VZV. Here, we present a genome-wide study on HSV-1. We also compared the evolutionary characteristics of HSV-1 with those for VZV. We demonstrate that, in contrast to VZV for which only a few ancient recombination events have been suggested, all HSV-1 genomes contain mosaic patterns of segments with different evolutionary origins. Thus, recombination seems to occur extremely frequent for HSV-1. We conclude by proposing a timescale for HSV-1 evolution, and by discussing putative underlying mechanisms for why these otherwise biologically similar viruses have such striking evolutionary differences.

Interacting proteins may often experience similar selection pressures. Thus, we may expect that neighbouring proteins in biological interaction networks evolve at similar rates. This has been previously shown for protein-protein interaction networks. Similarly, we find correlated rates of evolution of neighbours in networks based on co-expression, metabolism, and synthetic lethal genetic interactions. While the correlations are statistically significant, their magnitude is small, with network effects explaining only between 2% and 7% of the variation. The strongest known predictor of the rate of protein evolution remains expression level. We confirmed the previous observation that similar expression levels of neighbours indeed explain their similar evolution rates in protein-protein networks, and showed that the same is true for metabolic networks. In co-expression and synthetic lethal genetic interaction networks, however, neighbouring genes still show somewhat similar evolutionary rates even after simultaneously controlling for expression level, gene essentiality and gene length. Thus, similar expression levels and related functions (as inferred from co-expression and synthetic lethal interactions) seem to explain correlated evolutionary rates of network neighbours across all currently available types of biological networks.

We performed proteomics studies of salivas from the genome mouse (C57BL/6 strain) and the genome rat (BN/SsNHsd/Mcwi strain). Our goal was to identify salivary proteins with one or more of three characteristics that may indicate that they have been involved in adaptation: 1) rapid expansion of their gene families; 2) footprints of positive selection; and/or 3) sex-limited expression. The results of our proteomics studies allow direct comparison of the proteins expressed and their levels between the sexes of the two rodent species. Twelve members of the Mus musculus species-specific kallikrein subfamily Klk1b showed sex-limited expression in the mouse saliva proteomes. By contrast, we did not find any of the Rattus norvegicus species-specific kallikrein subfamily Klk1c proteins in male or female genome rat, nor transcripts in their submandibular glands. On the other hand, we detected expression of this family as transcripts in the submandibular glands of both sexes of Sprague-Dawley rats. Using the CODEML program in the PAML package, we demonstrate that the two rodent kallikrein subfamilies have apparently evolved rapidly under the influence of positive selection that continually remodeled the amino acid sites on the same face in the members of the subfamilies. Thus, although their kallikrein subfamily expansions were independent, this evolutionary pattern has occurred in parallel in the two rodent species, suggesting a form of convergent evolution at the molecular level. On the basis of this new data, we suggest that the previous speculative function of the species-specific rodent kallikreins as important solely in wound healing in males be investigated further. In addition to or instead of that function, we propose that their sex-limited expression, coupled with their rapid evolution may be clues to an as-yet-undetermined interaction between the sexes.

International audience; Class A G-protein-coupled receptors (GPCRs) constitute the largest family of transmembrane receptors in the human genome. Understanding the mechanisms which drove the evolution of such a large family would help understand the specificity of each GPCR sub-family with applications to drug design. To gain evolutionary information on class A GPCRs, we explored their sequence space by metric multidimensional scaling analysis (MDS). Three-dimensional mapping of human sequences shows a non-uniform distribution of GPCRs, organized in clusters that lay along four privileged directions. To interpret these directions, we projected supplementary sequences from different species onto the human space used as a reference. With this technique, we can easily monitor the evolutionary drift of several GPCR sub-families from cnidarians to humans. Results support a model of radiative evolution of class A GPCRs from a central node formed by peptide receptors. The privileged directions obtained from the MDS analysis are interpretable in terms of three main evolutionary pathways related to specific sequence determinants. The first pathway was initiated by a deletion in transmembrane helix 2 (TM2) and led to three sub-families by divergent evolution. The second pathway corresponds to the differentiation of the amine receptors. The third pathway corresponds to parallel evolution of several sub-families in relation with a covarion process involving proline residues in TM2 and TM5. As exemplified with GPCRs, the MDS projection technique is an important tool to compare orthologous sequence sets and to help decipher the mutational events that drove the evolution of protein families.

Traits that converge in appearance under similar environmental conditions among phylogenetically independent lineages are thought to represent adaptations to local environments. We tested for convergence in nest morphology and composition of birds breeding in two ecologically different locations in Canada: Churchill in northern Manitoba and Elgin in southeastern Ontario. We examined nests from four families of passerine birds (Turdidae: Turdus, Parulidae: Dendroica, Emberizidae: Passerculus and Fringillidae: Carduelis) where closely related populations or species breed in both locations. Nests of American Robins, Yellow Warblers, and Carduelis finches had heavier nest masses, and tended to have thicker nest-walls, in northern Manitoba compared with conspecifics or congenerics breeding in southeastern Ontario. Together, all species showed evidence for wider internal and external nest-cup diameters in northern Manitoba, while individual species showed varying patterns for internal nest-cup and external nest depths. American Robins, Yellow Warblers, and Carduelis finches in northern Manitoba achieved heavier nest masses in different ways. American Robins increased all materials in similar proportions, and Yellow Warblers and Common Redpolls used greater amounts of select materials. While changes in nest composition vary uniquely for each species, the pattern of larger nests in northern Manitoba compared to southeastern Ontario in three of our four phylogenetically-independent comparisons suggests that birds are adapting to similar selective pressures between locations.

BackgroundIntrogression plays an important role in evolution of plant species via its influences on genetic diversity and differentiation. Outcrossing determines the level of introgression but little is known about the relationships of outcrossing rates, genetic diversity, and differentiation particularly in a weedy taxon that coexists with its conspecific crop.Methodology/principal findingsEleven weedy rice (Oryza sativa f. spontanea) populations from China were analyzed using microsatellite (SSR) fingerprints to study outcrossing rate and its relationship with genetic variability and differentiation. To estimate outcrossing, six highly polymorphic SSR loci were used to analyze >5500 progeny from 216 weedy rice families, applying a mixed mating model; to estimate genetic diversity and differentiation, 22 SSR loci were analyzed based on 301 weedy individuals. Additionally, four weed-crop shared SSR loci were used to estimate the influence of introgression from rice cultivars on weedy rice differentiation. Outcrossing rates varied significantly (0.4~11.7%) among weedy rice populations showing relatively high overall Nei's genetic diversity (0.635). The observed heterozygosity was significantly correlated with outcrossing rates among populations (r² = 0.783; PConclusions/significanceOutcrossing rates can significantly affect heterozygosity of populations, which may shape the evolutionary potential of weedy rice. Introgression from the conspecific crop rice can influence the genetic differentiation and possibly evolution of its coexisting weedy rice populations.

Recent findings suggest a convergence of time scales between ecological and evolutionary processes which is usually explained in terms of rapid micro evolution resulting in evolution on ecological time scales. A similar convergence, however, can also emerge when slow ecological processes take place on evolutionary time scales. A good example of such a slow ecological process is the colonization of remote aquatic habitats by passively dispersed zooplankton. Using variation at the protein coding mitochondrial COI gene, we investigated the balance between mutation and migration as drivers of genetic diversity in two Branchipodopsis fairy shrimp species (Crustacea, Anostraca) endemic to remote temporary rock pool clusters at the summit of isolated mountaintops in central South Africa. We showed that both species colonized the region almost simultaneously c. 0.8 My ago, but exhibit contrasting patterns of regional genetic diversity and demographic history. The haplotype network of the common B. cf. wolfi showed clear evidence of 11 long distance dispersal events (up to 140 km) with five haplotypes that are shared among distant inselbergs, as well as some more spatially isolated derivates. Similar patterns were not observed for B. drakensbergensis presumably since this rarer species experienced a genetic bottleneck. We conclude that the observed genetic patterns reflect rare historic colonization events rather than frequent ongoing gene flow. Moreover, the high regional haplotype diversity combined with a high degree of haplotype endemicity indicates that evolutionary- (mutation) and ecological (migration) processes in this system operate on similar time scales. ispartof: PLoS One vol:6 issue:11 ispartof: location:United States status: published

The tiger is one of the most iconic extant animals, and its origin and evolution have been intensely debated. Fossils attributable to extant pantherine species-lineages are less than 2 MYA and the earliest tiger fossils are from the Calabrian, Lower Pleistocene. Molecular studies predict a much younger age for the divergence of modern tiger subspecies at

Our current understanding of sympatric speciation is that it occurs primarily through disruptive selection on ecological genes driven by competition, followed by reproductive isolation through reinforcement-like selection against inferior intermediates/heterozygotes. Our evolutionary model of selection on resource recognition and preference traits suggests a new mechanism for sympatric speciation. We find speciation can occur in three phases. First a polymorphism of functionally different phenotypes is established through evolution of specialization. On the gene level, regulatory functions have evolved in which some alleles are conditionally switched off (i.e. are silent). These alleles accumulate harmful mutations that potentially may be expressed in offspring through recombination. Second mating associated with resource preference invades because harmful mutations in parents are not expressed in the offspring when mating assortatively, thereby dividing the population into two pre-zygotically isolated resource-specialist lineages. Third, silent alleles that evolved in phase one now accumulate deleterious mutations over the following generations in a Bateson-Dobzhansky-Muller fashion, establishing a post-zygotic barrier to hybridization.

A prophylactic vaccine is needed to slow the spread of HIV-1 infection. Optimization of the wild-type envelope glycoproteins to create immunogens that can elicit effective neutralizing antibodies is a high priority. Starting with ten genes encoding subtype B HIV-1 gp120 envelope glycoproteins and using in vitro homologous DNA recombination, we created chimeric gp120 variants that were screened for their ability to bind neutralizing monoclonal antibodies. Hundreds of variants were identified with novel antigenic phenotypes that exhibit considerable sequence diversity. Immunization of rabbits with these gp120 variants demonstrated that the majority can induce neutralizing antibodies to HIV-1. One novel variant, called ST-008, induced significantly improved neutralizing antibody responses when assayed against a large panel of primary HIV-1 isolates. Further study of various deletion constructs of ST-008 showed that the enhanced immunogenicity results from a combination of effective DNA priming, an enhanced V3-based response, and an improved response to the constant backbone sequences.

A fundamental challenge of morphology is to identify the underlying evolutionary and developmental mechanisms leading to correlated phenotypic characters. Patterns and magnitudes of morphological integration and their association with environmental variables are essential for understanding the evolution of complex phenotypes, yet the nature of the relevant selective pressures remains poorly understood. In this study, the adaptive significance of morphological integration was evaluated through the association between feeding mechanics, ingestive behavior and craniofacial variation. Five capuchin species were examined, Cebus apella sensu stricto, Cebus libidinosus, Cebus nigritus, Cebus olivaceus and Cebus albifrons. Twenty three-dimensional landmarks were chosen to sample facial regions experiencing high strains during feeding, characteristics affecting muscular mechanical advantage and basicranial regions. Integration structure and magnitude between and within the oral and zygomatic subunits, between and within blocks maximizing modularity and within the face, the basicranium and the cranium were examined using partial-least squares, eigenvalue variance, integration indices compared inter-specifically at a common level of sampled population variance and cluster analyses. Results are consistent with previous findings reporting a relative constancy of facial and cranial correlation patterns across mammals, while covariance magnitudes vary. Results further suggest that food material properties structure integration among functionally-linked facial elements and possibly integration between the face and the basicranium. Hard-object-feeding capuchins, especially C. apella s.s., whose faces experience particularly high biomechanical loads are characterized by higher facial and cranial integration especially compared to C. albifrons, likely because morphotypes compromising feeding performance are selected against in species relying on obdurate fallback foods. This is the first study to report a link between food material properties and facial and cranial integration. Furthermore, results do not identify the consistent presence of cranial modules yielding support to suggestions that despite the distinct embryological imprints of its elements the cranium of placental mammals is not characterized by a modular architecture.

Local adaptation to different environments can promote mating isolation – either as an incidental by-product of trait divergence, or as a result of selection to avoid maladaptive mating. Numerous recent empirical examples point to the common influence of divergent natural selection on speciation based largely on evidence of strong pre-mating isolation between populations from different habitat types. Accumulating evidence for natural selection's influence on speciation is therefore no longer a challenge. The difficulty, rather, is in determining the mechanisms involved in the progress of adaptive divergence to speciation once barriers to gene flow are already present. Here, we present results of both laboratory and field experiments with Trinidadian guppies (Poecilia reticulata) from different environments, who do not show complete reproductive isolation despite adaptive divergence. We investigate patterns of mating isolation between populations that do and do not exchange migrants and show evidence for both by-product and reinforcement mechanisms depending on female ecology. Specifically, low-predation females discriminate against all high-predation males thus implying a by-product mechanism, whereas high-predation females only discriminate against low-predation males from further upstream in the same river, implying selection to avoid maladaptive mating. Our study thus confirms that mechanisms of adaptive speciation are not necessarily mutually exclusive and uncovers the complex ecology-geography interactions that underlie the evolution of mating isolation in nature.

Using a phylogenetic approach, the examination of 33 meiosis/meiosis-related genes in 12 Drosophila species, revealed nine independent gene duplications, involving the genes cav, mre11, meiS332, polo and mtrm. Evidence is provided that at least eight out of the nine gene duplicates are functional. Therefore, the rate at which Drosophila meiosis/meiosis-related genes are duplicated and retained is estimated to be 0.0012 per gene per million years, a value that is similar to the average for all Drosophila genes. It should be noted that by using a phylogenetic approach the confounding effect of concerted evolution, that is known to lead to overestimation of the duplication and retention rate, is avoided. This is an important issue, since even in our moderate size sample, evidence for long-term concerted evolution (lasting for more than 30 million years) was found for the meiS332 gene pair in species of the Drosophila subgenus. Most striking, in contrast to theoretical expectations, is the finding that genes that encode proteins that must follow a close stoichiometric balance, such as polo, mtrm and meiS332 have been found duplicated. The duplicated genes may be examples of gene neofunctionalization. It is speculated that meiosis duration may be a trait that is under selection in Drosophila and that it has different optimal values in different species.

Human respiratory syncytial virus (HRSV) has three surface glycoproteins: small hydrophobic (SH), attachment (G) and fusion (F), encoded by three consecutive genes (SH-G-F). A 270-nt fragment of the G gene is used to genotype HRSV isolates. This study genotyped and investigated the variability of the gene and amino acid sequences of the three surface proteins of HRSV strains collected from 1987 to 2005 from one center. Sixty original clinical isolates and 5 prototype strains were analyzed. Sequences containing SH, F and G genes were generated, and multiple alignments and phylogenetic trees were analyzed. Genetic variability by protein domains comparing virus genotypes was assessed. Complete sequences of the SH-G-F genes were obtained for all 65 samples: HRSV-A = 35; HRSV-B = 30. In group A strains, genotypes GA5 and GA2 were predominant. For HRSV-B strains, the genotype GB4 was predominant from 1992 to 1994 and only genotype BA viruses were detected in 2004–2005. Different genetic variability at nucleotide level was detected between the genes, with G gene being the most variable and the highest variability detected in the 270-nt G fragment that is frequently used to genotype the virus. High variability (>10%) was also detected in the signal peptide and transmembrane domains of the F gene of HRSV A strains. Variability among the HRSV strains resulting in non-synonymous changes was detected in hypervariable domains of G protein, the signal peptide of the F protein, a not previously defined domain in the F protein, and the antigenic site Ø in the pre-fusion F. Divergent trends were observed between HRSV -A and -B groups for some functional domains. A diverse population of HRSV -A and -B genotypes circulated in Houston during an 18 year period. We hypothesize that diverse sequence variation of the surface protein genes provide HRSV strains a survival advantage in a partially immune-protected community.

Current evolutionary theory suggests that many natural signaling systems evolved from preexisting cues. In aposematic systems, prey warning signals benefit both predator and prey. When the signal is highly beneficial, a third species often evolves to mimic the unpalatable species, exploiting the signaling system for its own protection. We investigated the evolutionary development of predator cue utilization and prey signaling in a digital predator-prey system in which mimicking prey could evolve to alter their appearance to resemble poison-free or poisonous prey. In predators, we observed rapid evolution of cue recognition (i.e. active behavioral responses) when presented with sufficiently poisonous prey. In addition, active signaling (i.e. mimicry) evolved in prey under all conditions that led to cue utilization. Thus we show that despite imperfect and dishonest signaling, given a high cost of consuming poisonous prey, complex systems of interspecific communication can evolve via predator cue recognition and prey signal manipulation. This provides evidence supporting hypotheses that cues may serve as stepping-stones in the evolution of more advanced communication systems and signals incorporating information about the environment.

Understanding of relationships between morphology and ecological performance can help to reveal how natural selection drives biological diversification. We investigate relationships between feeding behavior, foraging performance and morphology within a diverse group of teleost fishes, and examine the extent to which associations can be explained by evolutionary relatedness. Morphological adaptation associated with sediment sifting was examined using a phylogenetic linear discriminant analysis on a set of ecomorphological traits from 27 species of Neotropical cichlids. For most sifting taxa, feeding behavior could be effectively predicted by a linear discriminant function of ecomorphology across multiple clades of sediment sifters, and this pattern could not be explained by shared evolutionary history alone. Additionally, we tested foraging efficiency in seven Neotropical cichlid species, five of which are specialized benthic feeders with differing head morphology. Efficiency was evaluated based on the degree to which invertebrate prey could be retrieved at different depths of sediment. Feeding performance was compared both with respect to feeding mode and species using a phylogenetic ANCOVA, with substrate depth as a covariate. Benthic foraging performance was constant across sediment depths in non-sifters but declined with depth in sifters. The non-sifting Hypsophrys used sweeping motions of the body and fins to excavate large pits to uncover prey; this tactic was more efficient for consuming deeply buried invertebrates than observed among sediment sifters. Findings indicate that similar feeding performance among sediment-sifting cichlids extracting invertebrate prey from shallow sediment layers reflects constraints associated with functional morphology and, to a lesser extent, phylogeny.

The prevalence of sexual conflict in nature, as well as the supposedly arbitrary direction of the resulting coevolutionary trajectories, suggests that it may be an important driver of phenotypic divergence even in a constant environment. However, natural selection has long been central to the operation of sexual conflict within populations and may therefore constrain or otherwise direct divergence among populations. Ecological context may therefore matter with respect to the diversification of traits involved in sexual conflict, and if natural selection is sufficiently strong, such traits may evolve in correlation with environment, generating a pattern of ecologically-dependent parallel evolution. In this study we assess among-population divergence both within and between environments for several traits involved in sexual conflict. Using eight replicate populations of Drosophila melanogaster from a long-term evolution experiment, we measured remating rates and subsequent offspring production of females when housed with two separate males in sequence. We found no evidence of any variation in male reproductive traits (offense or defense). However, the propensity of females to remate diverged significantly among the eight populations with no evidence of any environmental effect, consistent with sexual conflict promoting diversification even in the absence of ecological differences. On the other hand, females adapted to one environment (ethanol) tended to produce a higher proportion of offspring sired by their first mate as compared to those adapted to the other (cadmium) environment, suggesting ecologically-based divergence of this conflict phenotype. Because we find evidence for both stochastic population divergence operating outside of an ecological context and environment-dependent divergence of traits under sexual conflict, the interaction of these two processes is an important topic for future work.

Lichen symbioses in the Pannariaceae associate an ascomycete and either cyanobacteria alone (usually Nostoc; bipartite thalli) or green algae and cyanobacteria (cyanobacteria being located in dedicated structures called cephalodia; tripartite thalli) as photosynthetic partners (photobionts). In bipartite thalli, cyanobacteria can either be restricted to a well-delimited layer within the thallus ('pannarioid' thalli) or spread over the thallus that becomes gelatinous when wet ('collematoid' thalli). We studied the collematoid genera Kroswia and Physma and an undescribed tripartite species along with representatives of the pannarioid genera Fuscopannaria, Pannaria and Parmeliella. Molecular inferences from 4 loci for the fungus and 1 locus for the photobiont and statistical analyses within a phylogenetic framework support the following: (a) several switches from pannarioid to collematoid thalli occured and are correlated with photobiont switches; the collematoid genus Kroswia is nested within the pannarioid genus Fuscopannaria and the collematoid genus Physma is sister to the pannarioid Parmeliella mariana group; (b) Nostoc associated with collematoid thalli in the Pannariaceae are related to that of the Collemataceae (which contains only collematoid thalli), and never associated with pannarioid thalli; Nostoc associated with pannarioid thalli also associate in other families with similar morphology; (c) ancestors of several lineages in the Pannariaceae developed tripartite thalli, bipartite thalli probably resulting from cephalodia emancipation from tripartite thalli which eventually evolved and diverged, as suggested by the same Nostoc present in the collematoid genus Physma and in the cephalodia of a closely related tripartite species; Photobiont switches and cephalodia emancipation followed by divergence are thus suspected to act as evolutionary drivers in the family Pannariaceae.

Contemporary patterns of genetic diversity and population connectivity within species can be influenced by both historical and contemporary barriers to gene flow. In the marine environment, present day oceanographic features such as currents, fronts and upwelling systems can influence dispersal of eggs/larvae and/juveniles/adults, shaping population substructuring. The Benguela Current system in the southeastern Atlantic is one of the oldest upwelling systems in the world, and provides a unique opportunity to investigate the relative influence of contemporary and historical mechanisms shaping the evolutionary history of warm-temperate fish species. Using the genetic variation in the mitochondrial DNA Control Region and eight nuclear microsatellite DNA loci, we identified the presence of two highly divergent populations in a vagile and warm-temperate fish species, Atractoscion aequidens, across the Benguela region. The geographical distributions of the two populations, on either side of the perennial upwelling cell, suggest a strong correlation between the oceanographic features of the system and the breakdown of gene flow within this species. Genetic divergence (mtDNA φ ST = 0.902, microsatellite F ST = 0.055: probability of genetic homogeneity for either marker = p

During long-term evolution of Escherichia coli, nutrient specialization is primarily driven by the accumulation of neutral mutations, rather than by tradeoffs, and can also be accompanied by general catabolic improvements., Evolutionary adaptation to a constant environment is often accompanied by specialization and a reduction of fitness in other environments. We assayed the ability of the Lenski Escherichia coli populations to grow on a range of carbon sources after 50,000 generations of adaptation on glucose. Using direct measurements of growth rates, we demonstrated that declines in performance were much less widespread than suggested by previous results from Biolog assays of cellular respiration. Surprisingly, there were many performance increases on a variety of substrates. In addition to the now famous example of citrate, we observed several other novel gains of function for organic acids that the ancestral strain only marginally utilized. Quantitative growth data also showed that strains with a higher mutation rate exhibited significantly more declines, suggesting that most metabolic erosion was driven by mutation accumulation and not by physiological tradeoffs. These reductions in growth by mutator strains were ameliorated by growth at lower temperature, consistent with the hypothesis that this metabolic erosion is largely caused by destabilizing mutations to the associated enzymes. We further hypothesized that reductions in growth rate would be greatest for substrates used most differently from glucose, and we used flux balance analysis to formulate this question quantitatively. To our surprise, we found no significant relationship between decreases in growth and dissimilarity to glucose metabolism. Taken as a whole, these data suggest that in a single resource environment, specialization does not mainly result as an inevitable consequence of adaptive tradeoffs, but rather due to the gradual accumulation of disabling mutations in unused portions of the genome., Author Summary Adaptation to a single constant environment is commonly expected to result in decreased performance in alternative conditions, or specialization. It has been proposed that, rather than occurring through the neutral accumulation of mutations in unused alternative pathways, this happens because loss of these pathways enhances fitness in the constant environment via “tradeoffs.” We examined growth rates across a variety of nutrients for 12 independent lineages of Escherichia coli that had evolved in the laboratory for decades in a glucose-containing medium. Surprisingly, after 20,000 generations there were actually widespread improvements in the use of alternative nutrients, rather than the expected declines. After 50,000 generations, however, we find that this trend reversed for those populations that evolved a much higher mutation rate. This indicates that high mutation rate, and not adaptive tradeoffs per se (as had been previously proposed), is the primary driver of specialization. These results caution against general assumptions about the importance of adaptive tradeoffs during evolution, and emphasize the key role that newly evolved changes in mutation rate can play in promoting niche specialization.

The Krüppel-associated box (KRAB) domain interacts with the nuclear hub protein TRIM28 to initiate or mediate chromatin-dependent processes like transcriptional repression, imprinting or suppression of endogenous retroviruses. The prototype KRAB domain initially identified in ZNF10/KOX1 encompasses two subdomains A and B that are found in hundreds of zinc finger transcription factors studied in human and murine genomes. Here we demonstrate for the first time transcriptional repressor activity of an amphibian KRAB domain. After sequence correction, the updated KRAB-AB domain of zinc finger protein XFIN from the frog Xenopus laevis was found to confer transcriptional repression in reporter assays in Xenopus laevis A6 kidney cells as well as in human HeLa, but not in the minnow Pimephales promelas fish cell line EPC. Binding of the XFIN KRAB-AB domain to human TRIM28 was demonstrated in a classical co-immunoprecipitation approach and visualized in a single-cell compartmentalization assay. XFIN-AB displayed reduced potency in repression as well as lower strength of interaction with TRIM28 compared to ZNF10 KRAB-AB. KRAB-B subdomain swapping between the two KRAB domains indicated that it was mainly the KRAB-B subdomain of XFIN that was responsible for its lower capacity in repression and binding to human TRIM28. In EPC fish cells, ZNF10 and XFIN KRAB repressor activity could be partially restored to low levels by adding exogenous human TRIM28. In contrast to XFIN, we did not find any transcriptional repression activity for the KRAB-like domain of human PRDM9 in HeLa cells. PRDM9 is thought to harbor an evolutionary older domain related to KRAB whose homologs even occur in invertebrates. Our results support the notion that functional bona fide KRAB domains which confer transcriptional repression and interact with TRIM28 most likely co-evolved together with TRIM28 at the beginning of tetrapode evolution.

Evolutionary algorithms are widespread heuristic methods inspired by natural evolution to solve difficult problems for which analytical approaches are not suitable. In many domains experimenters are not only interested in discovering optimal solutions, but also in finding the largest number of different solutions satisfying minimal requirements. However, the formulation of an effective performance measure describing these requirements, also known as fitness function, represents a major challenge. The difficulty of combining and weighting multiple problem objectives and constraints of possibly varying nature and scale into a single fitness function often leads to unsatisfactory solutions. Furthermore, selective reproduction of the fittest solutions, which is inspired by competition-based selection in nature, leads to loss of diversity within the evolving population and premature convergence of the algorithm, hindering the discovery of many different solutions. Here we present an alternative abstraction of artificial evolution, which does not require the formulation of a composite fitness function. Inspired from viability theory in dynamical systems, natural evolution and ethology, the proposed method puts emphasis on the elimination of individuals that do not meet a set of changing criteria, which are defined on the problem objectives and constraints. Experimental results show that the proposed method maintains higher diversity in the evolving population and generates more unique solutions when compared to classical competition-based evolutionary algorithms. Our findings suggest that incorporating viability principles into evolutionary algorithms can significantly improve the applicability and effectiveness of evolutionary methods to numerous complex problems of science and engineering, ranging from protein structure prediction to aircraft wing design.

Divergence of acoustic signals in a geographic scale results from diverse evolutionary forces acting in parallel and affecting directly inter-male vocal interactions among disjunct populations. Pleurodema thaul is a frog having an extensive latitudinal distribution in Chile along which males' advertisement calls exhibit an important variation. Using the playback paradigm we studied the evoked vocal responses of males of three populations of P. thaul in Chile, from northern, central and southern distribution. In each population, males were stimulated with standard synthetic calls having the acoustic structure of local and foreign populations. Males of both northern and central populations displayed strong vocal responses when were confronted with the synthetic call of their own populations, giving weaker responses to the call of the southern population. The southern population gave stronger responses to calls of the northern population than to the local call. Furthermore, males in all populations were stimulated with synthetic calls for which the dominant frequency, pulse rate and modulation depth were varied parametrically. Individuals from the northern and central populations gave lower responses to a synthetic call devoid of amplitude modulation relative to stimuli containing modulation depths between 30-100%, whereas the southern population responded similarly to all stimuli in this series. Geographic variation in the evoked vocal responses of males of P. thaul underlines the importance of inter-male interactions in driving the divergence of the acoustic traits and contributes evidence for a role of intra-sexual selection in the evolution of the sound communication system of this anuran.

Reassortments and point mutations are two major contributors to diversity of Influenza A virus; however, the link between these two processes is unclear. It has been suggested that reassortments provoke a temporary increase in the rate of amino acid changes as the viral proteins adapt to new genetic environment, but this phenomenon has not been studied systematically. Here, we use a phylogenetic approach to infer the reassortment events between the 8 segments of influenza A H3N2 virus since its emergence in humans in 1968. We then study the amino acid replacements that occurred in genes encoded in each segment subsequent to reassortments. In five out of eight genes (NA, M1, HA, PB1 and NS1), the reassortment events led to a transient increase in the rate of amino acid replacements on the descendant phylogenetic branches. In NA and HA, the replacements following reassortments were enriched with parallel and/or reversing replacements; in contrast, the replacements at sites responsible for differences between antigenic clusters (in HA) and at sites under positive selection (in NA) were underrepresented among them. Post-reassortment adaptive walks contribute to adaptive evolution in Influenza A: in NA, an average reassortment event causes at least 2.1 amino acid replacements in a reassorted gene, with, on average, 0.43 amino acid replacements per evolving post-reassortment lineage; and at least ∼9% of all amino acid replacements are provoked by reassortments., Author Summary Influenza A is a rapidly evolving virus with genome composed of eight distinct RNA molecules called segments. This genetic structure allows formation of new combinations of segments when a cell is coinfected by multiple viral strains, in a process called reassortment. While “antigenic drift” – the process of continuous accumulation of point mutations that change the antigenic properties of the viral proteins – is mainly responsible for the seasonal flu, the heaviest pandemics were caused by spread of novel reassortant strains and the associated radical “antigenic shift”. However, the association between these two types of processes has not been studied systematically. Here, we use the extensive available complete-genome sequencing data for Influenza A H3N2 subtype to infer the evolutionary timings of within-subtype reassortment events, and study the patterns of point amino acid-changing replacements that followed reassortments. We find that reassortments were often rapidly followed by replacements, which possibly compensated for the loss of fitness associated with the reassortment or explored newly accessible fitness peaks. These findings may be relevant for prediction of future pandemic strains of Influenza A.

The exon junction complex (EJC) plays important roles in RNA metabolisms and the development of eukaryotic organisms. MAGO (short form of MAGO NASHI) and Y14 (also Tsunagi or RBM8) are the EJC core components. Their biological roles have been well investigated in various species, but the evolutionary patterns of the two gene families and their protein-protein interactions are poorly known. Genome-wide survey suggested that the MAGO and Y14 two gene families originated in eukaryotic organisms with the maintenance of a low copy. We found that the two protein families evolved slowly; however, the MAGO family under stringent purifying selection evolved more slowly than the Y14 family that was under relative relaxed purifying selection. MAGO and Y14 were obliged to form heterodimer in a eukaryotic organism, and this obligate mode was plesiomorphic. Lack of binding of MAGO to Y14 as functional barrier was observed only among distantly species, suggesting that a slow co-evolution of the two protein families. Inter-protein co-evolutionary signal was further quantified in analyses of the Tol-MirroTree and co-evolution analysis using protein sequences. About 20% of the 41 significantly correlated mutation groups (involving 97 residues) predicted between the two families was clade-specific. Moreover, around half of the predicted co-evolved groups and nearly all clade-specific residues fell into the minimal interaction domains of the two protein families. The mutagenesis effects of the clade-specific residues strengthened that the co-evolution is required for obligate MAGO-Y14 heterodimerization mode. In turn, the obliged heterodimerization in an organism serves as a strong functional constraint for the co-evolution of the MAGO and Y14 families. Such a co-evolution allows maintaining the interaction between the proteins through large evolutionary time scales. Our work shed a light on functional evolution of the EJC genes in eukaryotes, and facilitates to understand the co-evolutionary processes among protein families.

Large-scale phylogenies provide a valuable source to study background diversification rates and investigate if the rates have changed over time. Unfortunately most large-scale, dated phylogenies are sparsely sampled (fewer than 5% of the described species) and taxon sampling is not uniform. Instead, taxa are frequently sampled to obtain at least one representative per subgroup (e.g. family) and thus to maximize diversity (diversified sampling). So far, such complications have been ignored, potentially biasing the conclusions that have been reached. In this study I derive the likelihood of a birth-death process with non-constant (time-dependent) diversification rates and diversified taxon sampling. Using simulations I test if the true parameters and the sampling method can be recovered when the trees are small or medium sized (fewer than 200 taxa). The results show that the diversification rates can be inferred and the estimates are unbiased for large trees but are biased for small trees (fewer than 50 taxa). Furthermore, model selection by means of Akaike's Information Criterion favors the true model if the true rates differ sufficiently from alternative models (e.g. the birth-death model is recovered if the extinction rate is large and compared to a pure-birth model). Finally, I applied six different diversification rate models – ranging from a constant-rate pure birth process to a decreasing speciation rate birth-death process but excluding any rate shift models – on three large-scale empirical phylogenies (ants, mammals and snakes with respectively 149, 164 and 41 sampled species). All three phylogenies were constructed by diversified taxon sampling, as stated by the authors. However only the snake phylogeny supported diversified taxon sampling. Moreover, a parametric bootstrap test revealed that none of the tested models provided a good fit to the observed data. The model assumptions, such as homogeneous rates across species or no rate shifts, appear to be violated.

Examining the geological past of our planet allows us to study periods of severe climatic and biological crises and recoveries, biotic and abiotic ecosystem fluctuations, and faunal and floral turnovers through time. Furthermore, the recovery dynamics of large predators provide a key for evaluation of the pattern and tempo of ecosystem recovery because predators are interpreted to react most sensitively to environmental turbulences. The end-Permian mass extinction was the most severe crisis experienced by life on Earth, and the common paradigm persists that the biotic recovery from the extinction event was unusually slow and occurred in a step-wise manner, lasting up to eight to nine million years well into the early Middle Triassic (Anisian) in the oceans, and even longer in the terrestrial realm. Here we survey the global distribution and size spectra of Early Triassic and Anisian marine predatory vertebrates (fishes, amphibians and reptiles) to elucidate the height of trophic pyramids in the aftermath of the end-Permian event. The survey of body size was done by compiling maximum standard lengths for the bony fishes and some cartilaginous fishes, and total size (estimates) for the tetrapods. The distribution and size spectra of the latter are difficult to assess because of preservation artifacts and are thus mostly discussed qualitatively. The data nevertheless demonstrate that no significant size increase of predators is observable from the Early Triassic to the Anisian, as would be expected from the prolonged and stepwise trophic recovery model. The data further indicate that marine ecosystems characterized by multiple trophic levels existed from the earliest Early Triassic onwards. However, a major change in the taxonomic composition of predatory guilds occurred less than two million years after the end-Permian extinction event, in which a transition from fish/amphibian to fish/reptile-dominated higher trophic levels within ecosystems became apparent.

Here we report the detection and localisation of chitin in the cuticle of the spinning ducts of both the spider Nephila edulis and the silkworm Bombyx mori. Our observations demonstrate that the duct walls of both animals contain chitin notwithstanding totally independent evolutionary pathways of the systems. We conclude that chitin may well be an essential component for the construction of spinning ducts; we further conclude that in both species chitin may indicate the evolutionary origin of the spinning ducts.

After the surprisingly low number of genes identified in the human genome, alternative splicing emerged as a major mechanism to generate protein diversity in higher eukaryotes. However, it is still not known if its prevalence along the genome evolution has contributed to the overall functional protein diversity or if it simply reflects splicing noise. The (βα)8 barrel or TIM barrel is one of the most frequent, versatile, and ancient fold encountered among enzymes. Here, we analyze the structural modifications present in TIM barrel proteins from the human genome product of alternative splicing events. We found that 87% of all splicing events involved deletions; most of these events resulted in protein fragments that corresponded to the (βα)2, (βα)4, (βα)5, (βα)6, and (βα)7 subdomains of TIM barrels. Because approximately 7% of all the splicing events involved internal β-strand substitutions, we decided, based on the genomic data, to design β-strand and α-helix substitutions in a well-studied TIM barrel enzyme. The biochemical characterization of one of the chimeric variants suggests that some of the splice variants in the human genome with β-strand substitutions may be evolving novel functions via either the oligomeric state or substrate specificity. We provide results of how the splice variants represent subdomains that correlate with the independently folding and evolving structural units previously reported. This work is the first to observe a link between the structural features of the barrel and a recurrent genetic mechanism. Our results suggest that it is reasonable to expect that a sizeable fraction of splice variants found in the human genome represent structurally viable functional proteins. Our data provide additional support for the hypothesis of the origin of the TIM barrel fold through the assembly of smaller subdomains. We suggest a model of how nature explores new proteins through alternative splicing as a mechanism to diversify the proteins encoded in the human genome.

Patterns of intraspecific geographic variation of signaling systems provide insight into the microevolutionary processes driving phenotypic divergence. The acoustic calls of bats are sensitive to diverse evolutionary forces, but processes that shape call variation are largely unexplored. In China, Rhinolophus ferrumequinum displays a diverse call frequency and inhabits a heterogeneous landscape, presenting an excellent opportunity for this kind of research. We quantified geographic variation in resting frequency (RF) of echolocation calls, estimated genetic structure and phylogeny of R. ferrumequinum populations, and combined this with climatic factors to test three hypotheses to explain acoustic variation: genetic drift, cultural drift, and local adaptation. Our results demonstrated significant regional divergence in frequency and phylogeny among the bat populations in China's northeast (NE), central-east (CE) and southwest (SW) regions. The CE region had higher frequencies than the NE and SW regions. Drivers of RF divergence were estimated in the entire range and just the CE/NE region (since these two regions form a clade). In both cases, RF divergence was not correlated with mtDNA or nDNA genetic distance, but was significantly correlated with geographic distance and mean annual temperature, indicating cultural drift and ecological selection pressures are likely important in shaping RF divergence among different regions in China.

Knowledge of the mechanisms that lead to reproductive isolation is essential for understanding population structure and speciation. While several models have been advanced to explain post-mating reproductive isolation, experimental data supporting most are indirect. Laboratory investigations of this phenomenon are typically carried out under benign conditions, which result in low rates of genetic change unlikely to initiate reproductive isolation. Previously, we described an experimental system using the yeast Saccharomyces cerevisiae where starvation served as a proxy to any stress that decreases reproduction and/or survivorship. We showed that novel lineages with restructured genomes quickly emerged in starved populations, and that these survivors were more fit than their ancestors when re-starved. Here we show that certain yeast lineages that survive starvation have become reproductively isolated from their ancestor. We further demonstrate that reproductive isolation arises from genomic rearrangements, whose frequency in starving yeast is several orders of magnitude greater than an unstarved control. By contrast, the frequency of point mutations is less than 2-fold greater. In a particular case, we observe that a starved lineage becomes reproductively isolated as a direct result of the stress-related accumulation of a single chromosome. We recapitulate this result by demonstrating that introducing an extra copy of one or several chromosomes into naïve, i.e. unstarved, yeast significantly diminishes their fertility. This type of reproductive barrier, whether arising spontaneously or via genetic manipulation, can be removed by making a lineage euploid for the altered chromosomes. Our model provides direct genetic evidence that reproductive isolation can arise frequently in stressed populations via genome restructuring without the precondition of geographic isolation.