Convergence, adaptation, and constraint

The joint evolution of traits and habitat: ontogenetic shifts in leaf morphology and wetland specialization in Lasthenia

The interplay between functional traits and habitat associations drives species' evolutionary responses to environmental heterogeneity, including processes such as adaptation, ecological speciation, and niche evolution. Seasonal variation is an aspect of the environment that varies across habitats, and could result in adaptive shifts in trait values across the life cycle of a plant. Here, we use phylogenetic comparative methods to evaluate the joint evolution of plant traits and habitat associations in Lasthenia (Asteraceae), a small clade of predominantly annual plants that have differentiated into an ecologically diverse range of habitats, including seasonal ephemeral wetlands known as vernal pools. Our results support the hypothesis that there is a link between the evolution of leaf morphology and the ecohydrological niche in Lasthenia, and, in the formation of aerenchyma (air space), differentiation between vernal pool and terrestrial taxa is fine-tuned to specific stages of plant ontogeny that reflects the evolution of heterophylly. Our findings demonstrate how the relationships between traits and habitat type can vary across the development of an organism, while highlighting a carefully considered comparative approach for examining correlated trait and niche evolution in a recently diversified and ecologically diverse plant clade.

On the evolutionary relationship between chondrocytes and osteoblasts

Vertebrates are the only animals that produce bone, but the molecular genetic basis for this evolutionary novelty remains obscure. Here, we synthesize information from traditional evolutionary and modern molecular genetic studies in order to generate a working hypothesis on the evolution of the gene regulatory network (GRN) underlying bone formation. Since transcription factors are often core components of GRNs (i.e., kernels), we focus our analyses on Sox9 and Runx2. Our argument centers on three skeletal tissues that comprise the majority of the vertebrate skeleton: immature cartilage, mature cartilage, and bone. Immature cartilage is produced during early stages of cartilage differentiation and can persist into adulthood, whereas mature cartilage undergoes additional stages of differentiation, including hypertrophy and mineralization. Functionally, histologically, and embryologically, these three skeletal tissues are very similar, yet unique, suggesting that one might have evolved from another. Traditional studies of the fossil record, comparative anatomy and embryology demonstrate clearly that immature cartilage evolved before mature cartilage or bone. Modern molecular approaches show that the GRNs regulating differentiation of these three skeletal cell fates are similar, yet unique, just like the functional and histological features of the tissues themselves. Intriguingly, the Sox9 GRN driving cartilage formation appears to be dominant to the Runx2 GRN of bone. Emphasizing an embryological and evolutionary transcriptomic view, we hypothesize that the Runx2 GRN underlying bone formation was co-opted from mature cartilage. We discuss how modern molecular genetic experiments, such as comparative transcriptomics, can test this hypothesis directly, meanwhile permitting levels of constraint and adaptation to be evaluated quantitatively. Therefore, comparative transcriptomics may revolutionize understanding of not only the clade-specific evolution of skeletal cells, but also the generation of evolutionary novelties, providing a modern paradigm for the evolutionary process.

Extent of QTL Reuse During Repeated Phenotypic Divergence of Sympatric Threespine Stickleback

How predictable is the genetic basis of phenotypic adaptation? Answering this question begins by estimating the repeatability of adaptation at the genetic level. Here, we provide a comprehensive estimate of the repeatability of the genetic basis of adaptive phenotypic evolution in a natural system. We used quantitative trait locus (QTL) mapping to discover genomic regions controlling a large number of morphological traits that have diverged in parallel between pairs of threespine stickleback (Gasterosteus aculeatus species complex) in Paxton and Priest lakes, British Columbia. We found that nearly half of QTL affected the same traits in the same direction in both species pairs. Another 40% influenced a parallel phenotypic trait in one lake but not the other. The remaining 10% of QTL had phenotypic effects in opposite directions in the two species pairs. Similarity in the proportional contributions of all QTL to parallel trait differences was about 0.4. Surprisingly, QTL reuse was unrelated to phenotypic effect size. Our results indicate that repeated use of the same genomic regions is a pervasive feature of parallel phenotypic adaptation, at least in sticklebacks. Identifying the causes of this pattern would aid prediction of the genetic basis of phenotypic evolution.

Past climate change on Sky Islands drives novelty in a core developmental gene network and its phenotype

BACKGROUND

A fundamental and enduring problem in evolutionary biology is to understand how populations differentiate in the wild, yet little is known about what role organismal development plays in this process. Organismal development integrates environmental inputs with the action of gene regulatory networks to generate the phenotype. Core developmental gene networks have been highly conserved for millions of years across all animals, and therefore, organismal development may bias variation available for selection to work on. Biased variation may facilitate repeatable phenotypic responses when exposed to similar environmental inputs and ecological changes. To gain a more complete understanding of population differentiation in the wild, we integrated evolutionary developmental biology with population genetics, morphology, paleoecology and ecology. This integration was made possible by studying how populations of the ant species Monomorium emersoni respond to climatic and ecological changes across five 'Sky Islands' in Arizona, which are mountain ranges separated by vast 'seas' of desert. Sky Islands represent a replicated natural experiment allowing us to determine how repeatable is the response of M. emersoni populations to climate and ecological changes at the phenotypic, developmental, and gene network levels.

RESULTS

We show that a core developmental gene network and its phenotype has kept pace with ecological and climate change on each Sky Island over the last ~90,000 years before present (BP). This response has produced two types of evolutionary change within an ant species: one type is unpredictable and contingent on the pattern of isolation of Sky lsland populations by climate warming, resulting in slight changes in gene expression, organ growth, and morphology. The other type is predictable and deterministic, resulting in the repeated evolution of a novel wingless queen phenotype and its underlying gene network in response to habitat changes induced by climate warming.

CONCLUSION

Our findings reveal dynamics of developmental gene network evolution in wild populations. This holds important implications: (1) for understanding how phenotypic novelty is generated in the wild; (2) for providing a possible bridge between micro- and macroevolution; and (3) for understanding how development mediates the response of organisms to past, and potentially, future climate change.

Seahorse Brood Pouch Transcriptome Reveals Common Genes Associated with Vertebrate Pregnancy

Viviparity (live birth) has evolved more than 150 times in vertebrates, and represents an excellent model system for studying the evolution of complex traits. There are at least 23 independent origins of viviparity in fishes, with syngnathid fishes (seahorses and pipefish) unique in exhibiting male pregnancy. Male seahorses and pipefish have evolved specialized brooding pouches that provide protection, gas exchange, osmoregulation, and limited nutrient provisioning to developing embryos. Pouch structures differ widely across the Syngnathidae, offering an ideal opportunity to study the evolution of reproductive complexity. However, the physiological and genetic changes facilitating male pregnancy are largely unknown. We used transcriptome profiling to examine pouch gene expression at successive gestational stages in a syngnathid with the most complex brood pouch morphology, the seahorse Hippocampus abdominalis. Using a unique time-calibrated RNA-seq data set including brood pouch at key stages of embryonic development, we identified transcriptional changes associated with brood pouch remodeling, nutrient and waste transport, gas exchange, osmoregulation, and immunological protection of developing embryos at conception, development and parturition. Key seahorse transcripts share homology with genes of reproductive function in pregnant mammals, reptiles, and other live-bearing fish, suggesting a common toolkit of genes regulating pregnancy in divergent evolutionary lineages.

The significance and scope of evolutionary developmental biology: a vision for the 21st century

Evolutionary developmental biology (evo-devo) has undergone dramatic transformations since its emergence as a distinct discipline. This paper aims to highlight the scope, power, and future promise of evo-devo to transform and unify diverse aspects of biology. We articulate key questions at the core of eleven biological disciplines-from Evolution, Development, Paleontology, and Neurobiology to Cellular and Molecular Biology, Quantitative Genetics, Human Diseases, Ecology, Agriculture and Science Education, and lastly, Evolutionary Developmental Biology itself-and discuss why evo-devo is uniquely situated to substantially improve our ability to find meaningful answers to these fundamental questions. We posit that the tools, concepts, and ways of thinking developed by evo-devo have profound potential to advance, integrate, and unify biological sciences as well as inform policy decisions and illuminate science education. We look to the next generation of evolutionary developmental biologists to help shape this process as we confront the scientific challenges of the 21st century.

Evolutionary determinants of cancer

Our understanding of cancer is being transformed by exploring clonal diversity, drug resistance, and causation within an evolutionary framework. The therapeutic resilience of advanced cancer is a consequence of its character as a complex, dynamic, and adaptive ecosystem engendering robustness, underpinned by genetic diversity and epigenetic plasticity. The risk of mutation-driven escape by self-renewing cells is intrinsic to multicellularity but is countered by multiple restraints, facilitating increasing complexity and longevity of species. But our own species has disrupted this historical narrative by rapidly escalating intrinsic risk. Evolutionary principles illuminate these challenges and provide new avenues to explore for more effective control.

SIGNIFICANCE

Lifetime risk of cancer now approximates to 50% in Western societies. And, despite many advances, the outcome for patients with disseminated disease remains poor, with drug resistance the norm. An evolutionary perspective may provide a clearer understanding of how cancer clones develop robustness and why, for us as a species, risk is now off the scale. And, perhaps, of what we might best do to achieve more effective control.

Evolutionary determinants of cancer

Our understanding of cancer is being transformed by exploring clonal diversity, drug resistance, and causation within an evolutionary framework. The therapeutic resilience of advanced cancer is a consequence of its character as a complex, dynamic, and adaptive ecosystem engendering robustness, underpinned by genetic diversity and epigenetic plasticity. The risk of mutation-driven escape by self-renewing cells is intrinsic to multicellularity but is countered by multiple restraints, facilitating increasing complexity and longevity of species. But our own species has disrupted this historical narrative by rapidly escalating intrinsic risk. Evolutionary principles illuminate these challenges and provide new avenues to explore for more effective control.

SIGNIFICANCE

Lifetime risk of cancer now approximates to 50% in Western societies. And, despite many advances, the outcome for patients with disseminated disease remains poor, with drug resistance the norm. An evolutionary perspective may provide a clearer understanding of how cancer clones develop robustness and why, for us as a species, risk is now off the scale. And, perhaps, of what we might best do to achieve more effective control.

Repeated evolution and the impact of evolutionary history on adaptation

BACKGROUND

Whether natural selection can erase the imprint of past evolutionary history from phenotypes has been a topic of much debate. A key source of evidence that present-day selection can override historically contingent effects comes from the repeated evolution of similar adaptations in different taxa. Yet classic examples of repeated evolution are often among closely related taxa, suggesting the likelihood that similar adaptations evolve is contingent on the length of time separating taxa. To resolve this, we performed a meta-analysis of published reports of repeated evolution.

RESULTS

Overall, repeated evolution was far more likely to be documented among closely related than distantly related taxa. However, not all forms of adaptation seemed to exhibit the same pattern. The evolution of similar behavior and physiology seemed frequent in distantly related and closely related taxa, while the repeated evolution of morphology was heavily skewed towards closely related taxa. Functionally redundant characteristics-alternative phenotypes that achieve the same functional outcome-also appeared less contingent.

CONCLUSIONS

If the literature provides a reasonable reflection of the incidence of repeated evolution in nature, our findings suggest that natural selection can overcome contingent effects to an extent, but it depends heavily on the aspect of the phenotype targeted by selection.

Adaptive radiation versus 'radiation' and 'explosive diversification': why conceptual distinctions are fundamental to understanding evolution

Adaptive radiation is the rise of a diversity of ecological roles and role-specific adaptations within a lineage. Recently, some researchers have begun to use 'adaptive radiation' or 'radiation' as synonymous with 'explosive species diversification'. This essay aims to clarify distinctions between these concepts, and the related ideas of geographic speciation, sexual selection, key innovations, key landscapes and ecological keys. Several examples are given to demonstrate that adaptive radiation and explosive diversification are not the same phenomenon, and that focusing on explosive diversification and the analysis of phylogenetic topology ignores much of the rich biology associated with adaptive radiation, and risks generating confusion about the nature of the evolutionary forces driving species diversification. Some 'radiations' involve bursts of geographic speciation or sexual selection, rather than adaptive diversification; some adaptive radiations have little or no effect on speciation, or even a negative effect. Many classic examples of 'adaptive radiation' appear to involve effects driven partly by geographic speciation, species' dispersal abilities, and the nature of extrinsic dispersal barriers; partly by sexual selection; and partly by adaptive radiation in the classical sense, including the origin of traits and invasion of adaptive zones that result in decreased diversification rates but add to overall diversity.

Implications of Nubian-Like Core Reduction Systems in Southern Africa for the Identification of Early Modern Human Dispersals

Lithic technologies have been used to trace dispersals of early human populations within and beyond Africa. Convergence in lithic systems has the potential to confound such interpretations, implying connections between unrelated groups. Due to their reductive nature, stone artefacts are unusually prone to this chance appearance of similar forms in unrelated populations. Here we present data from the South African Middle Stone Age sites Uitpanskraal 7 and Mertenhof suggesting that Nubian core reduction systems associated with Late Pleistocene populations in North Africa and potentially with early human migrations out of Africa in MIS 5 also occur in southern Africa during early MIS 3 and with no clear connection to the North African occurrence. The timing and spatial distribution of their appearance in southern and northern Africa implies technological convergence, rather than diffusion or dispersal. While lithic technologies can be a critical guide to human population flux, their utility in tracing early human dispersals at large spatial and temporal scales with stone artefact types remains questionable.

Iterative adaptive radiations of fossil canids show no evidence for diversity-dependent trait evolution

A long-standing hypothesis in adaptive radiation theory is that ecological opportunity constrains rates of phenotypic evolution, generating a burst of morphological disparity early in clade history. Empirical support for the early burst model is rare in comparative data, however. One possible reason for this lack of support is that most phylogenetic tests have focused on extant clades, neglecting information from fossil taxa. Here, I test for the expected signature of adaptive radiation using the outstanding 40-My fossil record of North American canids. Models implying time- and diversity-dependent rates of morphological evolution are strongly rejected for two ecologically important traits, body size and grinding area of the molar teeth. Instead, Ornstein-Uhlenbeck processes implying repeated, and sometimes rapid, attraction to distinct dietary adaptive peaks receive substantial support. Diversity-dependent rates of morphological evolution seem uncommon in clades, such as canids, that exhibit a pattern of replicated adaptive radiation. Instead, these clades might best be thought of as deterministic radiations in constrained Simpsonian subzones of a major adaptive zone. Support for adaptive peak models may be diagnostic of subzonal radiations. It remains to be seen whether early burst or ecological opportunity models can explain broader adaptive radiations, such as the evolution of higher taxa.

More than one way to evolve a weed: parallel evolution of US weedy rice through independent genetic mechanisms

Many different crop species were selected for a common suite of 'domestication traits', which facilitates their use for studies of parallel evolution. Within domesticated rice (Oryza sativa), there has also been independent evolution of weedy strains from different cultivated varieties. This makes it possible to examine the genetic basis of parallel weed evolution and the extent to which this process occurs through shared genetic mechanisms. We performed comparative QTL mapping of weediness traits using two recombinant inbred line populations derived from crosses between an indica crop variety and representatives of each of the two independently evolved weed strains found in US rice fields, strawhull (S) and blackhull awned (B). Genotyping-by-sequencing provided dense marker coverage for linkage map construction (average marker interval <0.25 cM), with 6016 and 13 730 SNPs mapped in F5 lines of the S and B populations, respectively. For some weediness traits (awn length, hull pigmentation and pericarp pigmentation), QTL mapping and sequencing of underlying candidate genes confirmed that trait variation was largely attributable to individual loci. However, for more complex quantitative traits (including heading date, panicle length and seed shattering), we found multiple QTL, with little evidence of shared genetic bases between the S and B populations or across previous studies of weedy rice. Candidate gene sequencing revealed causal genetic bases for 8 of 27 total mapped QTL. Together these findings suggest that despite the genetic bottleneck that occurred during rice domestication, there is ample genetic variation in this crop to allow agricultural weed evolution through multiple genetic mechanisms.

RAD-QTL Mapping Reveals Both Genome-Level Parallelism and Different Genetic Architecture Underlying the Evolution of Body Shape in Lake Whitefish (Coregonus clupeaformis) Species Pairs

Parallel changes in body shape may evolve in response to similar environmental conditions, but whether such parallel phenotypic changes share a common genetic basis is still debated. The goal of this study was to assess whether parallel phenotypic changes could be explained by genetic parallelism, multiple genetic routes, or both. We first provide evidence for parallelism in fish shape by using geometric morphometrics among 300 fish representing five species pairs of Lake Whitefish. Using a genetic map comprising 3438 restriction site−associated DNA sequencing single-nucleotide polymorphisms, we then identified quantitative trait loci underlying body shape traits in a backcross family reared in the laboratory. A total of 138 body shape quantitative trait loci were identified in this cross, thus revealing a highly polygenic architecture of body shape in Lake Whitefish. Third, we tested for evidence of genetic parallelism among independent wild populations using both a single-locus method (outlier analysis) and a polygenic approach (analysis of covariation among markers). The single-locus approach provided limited evidence for genetic parallelism. However, the polygenic analysis revealed genetic parallelism for three of the five lakes, which differed from the two other lakes. These results provide evidence for both genetic parallelism and multiple genetic routes underlying parallel phenotypic evolution in fish shape among populations occupying similar ecological niches.

Genetic and developmental basis for parallel evolution and its significance for hominoid evolution

Greater understanding of ape comparative anatomy and evolutionary history has brought a general appreciation that the hominoid radiation is characterized by substantial homoplasy.(1-4) However, little consensus has been reached regarding which features result from repeated evolution. This has important implications for reconstructing ancestral states throughout hominoid evolution, including the nature of the Pan-Homo last common ancestor (LCA). Advances from evolutionary developmental biology (evo-devo) have expanded the diversity of model organisms available for uncovering the morphogenetic mechanisms underlying instances of repeated phenotypic change. Of particular relevance to hominoids are data from adaptive radiations of birds, fish, and even flies demonstrating that parallel phenotypic changes often use similar genetic and developmental mechanisms. The frequent reuse of a limited set of genes and pathways underlying phenotypic homoplasy suggests that the conserved nature of the genetic and developmental architecture of animals can influence evolutionary outcomes. Such biases are particularly likely to be shared by closely related taxa that reside in similar ecological niches and face common selective pressures. Consideration of these developmental and ecological factors provides a strong theoretical justification for the substantial homoplasy observed in the evolution of complex characters and the remarkable parallel similarities that can occur in closely related taxa. Thus, as in other branches of the hominoid radiation, repeated phenotypic evolution within African apes is also a distinct possibility. If so, the availability of complete genomes for each of the hominoid genera makes them another model to explore the genetic basis of repeated evolution.

Space use variation in co-occurring sister species: response to environmental variation or competition?

Coexistence often involves niche differentiation either as the result of environmental divergence, or in response to competition. Disentangling the causes of such divergence requires that environmental variation across space is taken into account, which is rarely done in empirical studies. We address the role of environmental variation versus competition in coexistence between two rodent species: Rhabdomys bechuanae (bechuanae) and Rhabdomys dilectus dilectus (dilectus) comparing their habitat preference and home range (HR) size in areas with similar climates, where their distributions abut (allopatry) or overlap (sympatry). Using Outlying Mean Index analyses, we test whether habitat characteristics of the species deviate significantly from a random sample of available habitats. In allopatry, results suggest habitat selection: dilectus preferring grasslands with little bare soil while bechuanae occurring in open shrublands. In sympatry, shrubland type habitats dominate and differences are less marked, yet dilectus selects habitats with more cover than bechuanae. Interestingly, bechuanae shows larger HRs than dilectus, and both species display larger HRs in sympatry. Further, HR overlaps between species are lower than expected. We discuss our results in light of data on the phylogeography of the genus and propose that evolution in allopatry resulted in adaptation leading to different habitat preferences, even at their distribution margins, a divergence expected to facilitate coexistence. However, since sympatry occurs in sites where environmental characteristics do not allow complete species separation, competition may explain reduced inter-species overlap and character displacement in HR size. This study reveals that both environmental variation and competition may shape species coexistence.

Two developmentally temporal quantitative trait loci underlie convergent evolution of increased branchial bone length in sticklebacks

In convergent evolution, similar phenotypes evolve repeatedly in independent populations, often reflecting adaptation to similar environments. Understanding whether convergent evolution proceeds via similar or different genetic and developmental mechanisms offers insight towards the repeatability and predictability of evolution. Oceanic populations of threespine stickleback fish, Gasterosteus aculeatus, have repeatedly colonized countless freshwater lakes and streams, where new diets lead to morphological adaptations related to feeding. Here, we show that heritable increases in branchial bone length have convergently evolved in two independently derived freshwater stickleback populations. In both populations, an increased bone growth rate in juveniles underlies the convergent adult phenotype, and one population also has a longer cartilage template. Using F2 crosses from these two freshwater populations, we show that two quantitative trait loci (QTL) control branchial bone length at distinct points in development. In both populations, a QTL on chromosome 21 controls bone length throughout juvenile development, and a QTL on chromosome 4 controls bone length only in adults. In addition to these similar developmental profiles, these QTL show similar chromosomal locations in both populations. Our results suggest that sticklebacks have convergently evolved longer branchial bones using similar genetic and developmental programmes in two independently derived populations.

Parallel evolution in Ugandan crater lakes: repeated evolution of limnetic body shapes in haplochromine cichlid fish

Background

The enormous diversity found in East African cichlid fishes in terms of morphology, coloration, and behavior have made them a model for the study of speciation and adaptive evolution. In particular, haplochromine cichlids, by far the most species-rich lineage of cichlids, are a well-known textbook example for parallel evolution. Southwestern Uganda is an area of high tectonic activity, and is home to numerous crater lakes. Many Ugandan crater lakes were colonized, apparently independently, by a single lineage of haplochromine cichlids. Thereby, this system could be considered a natural experiment in which one can study the interaction between geographical isolation and natural selection promoting phenotypic diversification.

Results

We sampled 13 crater lakes and six potentially-ancestral older lakes and, using both mitochondrial and microsatellite markers, discovered strong genetic and morphological differentiation whereby (a) geographically close lakes tend to be genetically more similar and (b) three different geographic areas seem to have been colonized by three independent waves of colonization from the same source population. Using a geometric morphometric approach, we found that body shape elongation (i.e. a limnetic morphology) evolved repeatedly from the ancestral deeper-bodied benthic morphology in the clear and deep crater lake habitats.

Conclusions

A pattern of strong genetic and morphological differentiation was observed in the Ugandan crater lakes. Our data suggest that body shape changes have repeatedly evolved into a more limnetic-like form in several Ugandan crater lakes after independent waves of colonization from the same source population. The observed morphological changes in crater lake cichlids are likely to result from a common selective regime.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0287-3) contains supplementary material, which is available to authorized users.

Correlated evolution of microhabitat, morphology, and behavior in West Indian Anolis lizards: a test of the habitat matrix model

The habitat matrix model (HMM) explains convergence among arboreal animals as a result of the correlated evolution of morphology, locomotor mode, and habitat use. Although the HMM has generated important insights into the ecology of arboreal species, these tests have left a gap in the habitat-behavior-morphology story by focusing primarily on locomotor performance in lab and field experiments and thus failing to include data on locomotor behavior of undisturbed animals in the wild. We combined data on undisturbed locomotion, habitat use, and morphology for 31 species of arboreal lizard in the genusAnolisand used these data to test nine specific predictions arising from the HMM. We find strong support for nearly all aspects of this model. The addition of data on locomotion by undisturbed wild animals offers a more direct and compelling case for the HMM than most previous tests.

Co-option of alternate sperm activation programs in the evolution of self-fertile nematodes

Self-fertility evolved independently in three species of Caenorhabditis, yet the underlying genetic changes remain unclear. This transition required that XX animals acquire the ability to produce sperm and then signal those sperm to activate and fertilise oocytes. Here, we show that all genes that regulate sperm activation in C. elegans are conserved throughout the genus, even in male/female species. By using gene editing, we show that C. elegans and C. briggsae hermaphrodites use the SPE-8 tyrosine kinase pathway to activate sperm, whereas C. tropicalis hermaphrodites use a TRY-5 serine protease pathway. Finally, our analysis of double mutants shows that these pathways were redundant in ancestral males. Thus, newly evolving hermaphrodites became self-fertile by co-opting either of the two redundant male programs. The existence of these alternatives helps explain the frequent origin of self-fertility in nematode lineages. This work also demonstrates that the new genome-editing techniques allow unprecedented power and precision in evolutionary studies.

Evolutionary Pathways to Break Down the Resistance of Allelic Versions of the PVY Resistance Gene va

Emergence of viral genotypes can make control strategies based on resistance genes ineffective. A few years after the deployment of tobacco genotypes carrying alleles of the Potato virus Y (PVY) recessive resistance gene va, virulent PVY isolates have been reported, suggesting the low durability of va. To have a broader view of the evolutionary processes involved in PVY adaptation to va, we studied mutational pathways leading to the emergence of PVY resistance-breaking populations. The viral genome-linked protein (VPg) has been described to be potentially involved in va adaptation. Analyses of the VPg sequence of PVY isolates sampled from susceptible and resistant tobacco allowed us to identify mutations in the central part of the VPg. Analysis of the virulence of wild-type isolates with known VPg sequences and of mutated versions of PVY infectious clones allowed us to (i) validate VPg as the PVY virulence factor corresponding to va, (ii) highlight the fact that virulence gain in PVY occurs rapidly and preferentially by substitution at position AA₁₀₅ in the VPg, and (iii) show that the 101G substitution in the VPg of a PVY^C isolate is responsible for cross-virulence toward two resistance sources. Moreover, it appears that the evolutionary pathway of PVY adaptation to va depends on both virus and host genetic backgrounds.

Stochastic simulation in systems biology

Natural systems are, almost by definition, heterogeneous: this can be either a boon or an obstacle to be overcome, depending on the situation. Traditionally, when constructing mathematical models of these systems, heterogeneity has typically been ignored, despite its critical role. However, in recent years, stochastic computational methods have become commonplace in science. They are able to appropriately account for heterogeneity; indeed, they are based around the premise that systems inherently contain at least one source of heterogeneity (namely, intrinsic heterogeneity). In this mini-review, we give a brief introduction to theoretical modelling and simulation in systems biology and discuss the three different sources of heterogeneity in natural systems. Our main topic is an overview of stochastic simulation methods in systems biology. There are many different types of stochastic methods. We focus on one group that has become especially popular in systems biology, biochemistry, chemistry and physics. These discrete-state stochastic methods do not follow individuals over time; rather they track only total populations. They also assume that the volume of interest is spatially homogeneous. We give an overview of these methods, with a discussion of the advantages and disadvantages of each, and suggest when each is more appropriate to use. We also include references to software implementations of them, so that beginners can quickly start using stochastic methods for practical problems of interest.

Widespread signals of convergent adaptation to high altitude in Asia and america

Living at high altitude is one of the most difficult challenges that humans had to cope with during their evolution. Whereas several genomic studies have revealed some of the genetic bases of adaptations in Tibetan, Andean, and Ethiopian populations, relatively little evidence of convergent evolution to altitude in different continents has accumulated. This lack of evidence can be due to truly different evolutionary responses, but it can also be due to the low power of former studies that have mainly focused on populations from a single geographical region or performed separate analyses on multiple pairs of populations to avoid problems linked to shared histories between some populations. We introduce here a hierarchical Bayesian method to detect local adaptation that can deal with complex demographic histories. Our method can identify selection occurring at different scales, as well as convergent adaptation in different regions. We apply our approach to the analysis of a large SNP data set from low- and high-altitude human populations from America and Asia. The simultaneous analysis of these two geographic areas allows us to identify several candidate genome regions for altitudinal selection, and we show that convergent evolution among continents has been quite common. In addition to identifying several genes and biological processes involved in high-altitude adaptation, we identify two specific biological pathways that could have evolved in both continents to counter toxic effects induced by hypoxia.

Origins and functional diversification of salinity-responsive Na(+) , K(+) ATPase α1 paralogs in salmonids

The Salmoniform whole-genome duplication is hypothesized to have facilitated the evolution of anadromy, but little is known about the contribution of paralogs from this event to the physiological performance traits required for anadromy, such as salinity tolerance. Here, we determined when two candidate, salinity-responsive paralogs of the Na(+) , K(+) ATPase α subunit (α1a and α1b) evolved and studied their evolutionary trajectories and tissue-specific expression patterns. We found that these paralogs arose during a small-scale duplication event prior to the Salmoniform, but after the teleost, whole-genome duplication. The 'freshwater paralog' (α1a) is primarily expressed in the gills of Salmoniformes and an unduplicated freshwater sister species (Esox lucius) and experienced positive selection in the freshwater ancestor of Salmoniformes and Esociformes. Contrary to our predictions, the 'saltwater paralog' (α1b), which is more widely expressed than α1a, did not experience positive selection during the evolution of anadromy in the Coregoninae and Salmonine. To determine whether parallel mutations in Na(+) , K(+) ATPase α1 may contribute to salinity tolerance in other fishes, we studied independently evolved salinity-responsive Na(+) , K(+) ATPase α1 paralogs in Anabas testudineus and Oreochromis mossambicus. We found that a quarter of the mutations occurring between salmonid α1a and α1b in functionally important sites also evolved in parallel in at least one of these species. Together, these data argue that paralogs contributing to salinity tolerance evolved prior to the Salmoniform whole-genome duplication and that strong selection and/or functional constraints have led to parallel evolution in salinity-responsive Na(+) , K(+) ATPase α1 paralogs in fishes.

Ancient homology underlies adaptive mimetic diversity across butterflies

Convergent evolution provides a rare, natural experiment with which to test the predictability of adaptation at the molecular level. Little is known about the molecular basis of convergence over macro-evolutionary timescales. Here we use a combination of positional cloning, population genomic resequencing, association mapping and developmental data to demonstrate that positionally orthologous nucleotide variants in the upstream region of the same gene, WntA, are responsible for parallel mimetic variation in two butterfly lineages that diverged >65 million years ago. Furthermore, characterization of spatial patterns of WntA expression during development suggests that alternative regulatory mechanisms underlie wing pattern variation in each system. Taken together, our results reveal a strikingly predictable molecular basis for phenotypic convergence over deep evolutionary time.

Little is known about the genetic basis of convergent evolution in deeply diverged species. Here, the authors show that variation in the WntA gene is associated with parallel wing pattern variation in two butterflies that diverged more than 65 million years ago.

Interspecific tests of allelism reveal the evolutionary timing and pattern of accumulation of reproductive isolation mutations

Despite extensive theory, little is known about the empirical accumulation and evolutionary timing of mutations that contribute to speciation. Here we combined QTL (Quantitative Trait Loci) analyses of reproductive isolation, with information on species evolutionary relationships, to reconstruct the order and timing of mutations contributing to reproductive isolation between three plant (Solanum) species. To evaluate whether reproductive isolation QTL that appear to coincide in more than one species pair are homologous, we used cross-specific tests of allelism and found evidence for both homologous and lineage-specific (non-homologous) alleles at these co-localized loci. These data, along with isolation QTL unique to single species pairs, indicate that >85% of isolation-causing mutations arose later in the history of divergence between species. Phylogenetically explicit analyses of these data support non-linear models of accumulation of hybrid incompatibility, although the specific best-fit model differs between seed (pairwise interactions) and pollen (multi-locus interactions) sterility traits. Our findings corroborate theory that predicts an acceleration ('snowballing') in the accumulation of isolation loci as lineages progressively diverge, and suggest different underlying genetic bases for pollen versus seed sterility. Pollen sterility in particular appears to be due to complex genetic interactions, and we show this is consistent with a snowball model where later arising mutations are more likely to be involved in pairwise or multi-locus interactions that specifically involve ancestral alleles, compared to earlier arising mutations.

Repeated evolutionary changes of leaf morphology caused by mutations to a homeobox gene

Elucidating the genetic basis of morphological changes in evolution remains a major challenge in biology. Repeated independent trait changes are of particular interest because they can indicate adaptation in different lineages or genetic and developmental constraints on generating morphological variation. In animals, changes to "hot spot" genes with minimal pleiotropy and large phenotypic effects underlie many cases of repeated morphological transitions. By contrast, only few such genes have been identified from plants, limiting cross-kingdom comparisons of the principles of morphological evolution. Here, we demonstrate that the REDUCED COMPLEXITY (RCO) locus underlies more than one naturally evolved change in leaf shape in the Brassicaceae. We show that the difference in leaf margin dissection between the sister species Capsella rubella and Capsella grandiflora is caused by cis-regulatory variation in the homeobox gene RCO-A, which alters its activity in the developing lobes of the leaf. Population genetic analyses in the ancestral C. grandiflora indicate that the more-active C. rubella haplotype is derived from a now rare or lost C. grandiflora haplotype via additional mutations. In Arabidopsis thaliana, the deletion of the RCO-A and RCO-B genes has contributed to its evolutionarily derived smooth leaf margin, suggesting the RCO locus as a candidate for an evolutionary hot spot. We also find that temperature-responsive expression of RCO-A can explain the phenotypic plasticity of leaf shape to ambient temperature in Capsella, suggesting a molecular basis for the well-known negative correlation between temperature and leaf margin dissection.

The evolution of bipedal running in lizards suggests a consequential origin may be exploited in later lineages

The origin of bipedal locomotion in lizards is unclear. Modeling studies have suggested that bipedalism may be an exaptation, a byproduct of features originally designed to increase maneuverability, which were only later exploited. Measurement of the body center of mass (BCOM) in 124 species of lizards confirms a significant rearward shift among bipedal lineages. Further racetrack trials showed a significant acceleration threshold between bipedal and quadrupedal runs. These suggest good general support for a passive bipedal model, in which the combination of these features lead to passive lifting of the front of the body. However, variation in morphology could only account for 56% of the variation in acceleration thresholds, suggesting that dynamics have a significant influence on bipedalism. Deviation from the passive bipedal model was compared with node age, supporting an increase in the influence of dynamics over time. Together, these results show that bipedalism may have first arisen as a consequence of acceleration and a rearward shift in the BCOM, but subsequent linages have exploited this consequence to become bipedal more often, suggesting that bipedalism in lizards may convey some advantage. Exploitation of bipedalism was also associated with increased rates of phenotypic diversity, suggesting exploiting bipedalism may promote adaptive radiation.