Correlated evolution of mating strategy and inbreeding depression within and among populations of the hermaphroditic snail Physa acuta

Inbreeding depression is one of the main forces opposing the evolution of self-fertilization. Of central importance is the hypothesis that inbreeding depression and selfing coevolve antagonistically, generating either low selfing rate and high inbreeding depression or vice versa. However, there is limited evidence for this coevolution within species. We investigated this topic in the hermaphroditic snail Physa acuta. In this species, isolated individuals delay the onset of egg laying compared to individuals having access to mates. Longer delays (''waiting times'') indicate more intense selfing avoidance. We measured inbreeding depression and waiting time in a large quantitative-genetic experiment (281 outbred families derived from 26 natural populations). We observed large genetic variance for both traits and a strong positive genetic covariance between them, most of which resided within rather than among populations. It means that, within populations, individuals with higher mutation load avoided selfing more strongly on average. This genetic covariance may result from pleiotropy and/or linkage disequilibrium. Whatever its genetic architecture, the fact it emerges specifically when individuals are deprived of mates suggests it is not fortuitous and rather reflects the action of natural selection. We conclude that a diversity of mating strategies can arise within populations subjected to variation in inbreeding depression.

Bayesian analysis of correlated evolution of discrete characters by reversible-jump Markov chain Monte Carlo

We describe a Bayesian method for investigating correlated evolution of discrete binary traits on phylogenetic trees. The method fits a continuous-time Markov model to a pair of traits, seeking the best fitting models that describe their joint evolution on a phylogeny. We employ the methodology of reversible-jump ( RJ) Markov chain Monte Carlo to search among the large number of possible models, some of which conform to independent evolution of the two traits, others to correlated evolution. The RJ Markov chain visits these models in proportion to their posterior probabilities, thereby directly estimating the support for the hypothesis of correlated evolution. In addition, the RJ Markov chain simultaneously estimates the posterior distributions of the rate parameters of the model of trait evolution. These posterior distributions can be used to test among alternative evolutionary scenarios to explain the observed data. All results are integrated over a sample of phylogenetic trees to account for phylogenetic uncertainty. We implement the method in a program called RJ Discrete and illustrate it by analyzing the question of whether mating system and advertisement of estrus by females have coevolved in the Old World monkeys and great apes.

Constrained models of evolution lead to improved prediction of functional linkage from correlated gain and loss of genes

Motivation: We compare phylogenetic approaches for inferring functional gene links. The approaches detect independent instances of the correlated gain and loss of pairs of genes from species' genomes. We investigate the effect on results of basing evidence of correlations on two phylogenetic approaches, Dollo parsminony and maximum likelihood (ML). We further examine the effect of constraining the ML model by fixing the rate of gene gain at a low value, rather than estimating it from the data. Results: We detect correlated evolution among a test set of pairs of yeast (Saccharomyces cerevisiae) genes, with a case study of 21 eukaryotic genomes and test data derived from known yeast protein complexes. If the rate at which genes are gained is constrained to be low, ML achieves by far the best results at detecting known functional links. The model then has fewer parameters but it is more realistic by preventing genes from being gained more than once. Availability: BayesTraits by M. Pagel and A. Meade, and a script to configure and repeatedly launch it by D. Barker and M. Pagel, are available at http://www.evolution.reading.ac.uk .

Phylogeny and evolution of anomalous roots in Daviesia (Fabaceae : Mirbelieae)

The phylogeny of the Australian legume genus Daviesia was estimated using sequences of the internal transcribed spacers of nuclear ribosomal DNA. Partial congruence was found with previous analyses using morphology, including strong support for monophyly of the genus and for a sister group relationship between the clade D. pachyloma and the rest of the genus. A previously unplaced bird-pollinated species, anceps + D. D. epiphyllum, was well supported as sister to the only other bird-pollinated species in the genus, D. speciosa, indicating a single origin of bird pollination in their common ancestor. Other morphological groups within Daviesia were not supported and require reassessment. A strong and previously unreported sister clade of Daviesia consists of the two monotypic genera Erichsenia and Viminaria. These share phyllode-like leaves and indehiscent fruits. The evolutionary history of cord roots, which have anomalous secondary thickening, was explored using parsimony. Cord roots are limited to three separate clades but have a complex history involving a small number of gains (most likely 0-3) and losses (0-5). The anomalous structure of cord roots ( adventitious vascular strands embedded in a parenchymatous matrix) may facilitate nutrient storage, and the roots may be contractile. Both functions may be related to a postfire resprouting adaptation. Alternatively, cord roots may be an adaptation to the low-nutrient lateritic soils of Western Australia. However, tests for association between root type, soil type, and growth habit were equivocal, depending on whether the variables were treated as phylogenetically dependent (insignificant) or independent ( significant).

Latitudinal patterns in plant defense: evolution of cardenolides, their toxicity and induction following herbivory.

Attempts over the past 50 years to explain variation in the abundance, distribution and diversity of plant secondary compounds gave rise to theories of plant defense. Remarkably, few phylogenetically robust tests of these long-standing theories have been conducted. Using >50 species of milkweed (Asclepias spp.), we show that variation among plant species in the induction of toxic cardenolides is explained by latitude, with higher inducibility evolving more frequently at lower latitudes. We also found that: (1) the production of cardenolides showed positive-correlated evolution with the diversity of cardenolides, (2) greater cardenolide investment by a species is accompanied by an increase in an estimate of toxicity (measured as chemical polarity) and (3) instead of trading off, constitutive and induced cardenolides were positively correlated. Analyses of root and shoot cardenolides showed concordant patterns. Thus, milkweed species from lower latitudes are better defended with higher inducibility, greater diversity and added toxicity of cardenolides.

The evolution of host associations in the parasitic wasp genus Ichneumon (Hymenoptera: Ichneumonidae): convergent adaptations to host pupation sites

Background: The diversification of organisms with a parasitic lifestyle is often tightly linked to the evolution of their host associations. If a tight host association exists, closely related species tend to attack closely related hosts; host associations are less stable if associations are determined by more plastic traits like parasitoid searching and oviposition behaviour. The pupal-parasitoids of the genus Ichneumon attack a variety of macrolepidopteran hosts.They are either monophagous or polyphagous, and therefore offer a promissing system to investigate the evolution of host associations. Ichneumon was previously divided into two groups based on general body shape; however, a stout shape has been suggested as an adaptation to buried host pupation sites, and might thus not represent a reliable phylogenetic character. Results: We here reconstruct the first molecular phylogeny of the genus Ichneumon using two mitochondrial (CO1 and NADH1) and one nuclear marker (28S). The resulting phylogeny only supports monophyly of Ichneumon when Ichneumon lugens Gravenhorst, 1829 (formerly in Chasmias, stat. rev.) and Ichneumon deliratorius Linnaeus, 1758 (formerly Coelichneumon) are included. Neither parasitoid species that attack hosts belonging to one family nor those attacking butterflies (Rhopalocera) form monophyletic clades. Ancestral state reconstructions suggest multiple transitions between searching for hosts above versus below ground and between a stout versus elongated body shape. A model assuming correlated evolution between the two characters was preferred over independent evolution of host-searching niche and body shape. Conclusions: Host relations, both in terms of phylogeny and ecology, evolved at a high pace in the genus Ichneumon. Numerous switches between hosts of different lepidopteran families have occurred, a pattern that seems to be the rule among idiobiont parasitoids. A stout body and antennal shape in the parasitoid female is confirmed as an ecological adaptation to host pupation sites below ground and has evolved convergently several times. Morphological characters that might be involved in adaptation to hosts should be avoided as diagnostic characters for phylogeny and classification, as they can be expected to show high levels of homoplasy.

The ecological basis of life history variation in marsupials

Our understanding of the diversity of mammalian life histories is based almost exclusively on eutherian mammals, in which the slow-fast continuum persists even after controlling for effects of body size and phylogeny. In this paper, we use modern comparative methods to test the extent to which this eutherian-based framework can be extrapolated to metatherian mammals. First, we examine the pattern of covariation among life history traits, and second, we test for correlations between variation in life history and variation in six candidate ecological variables: type of diet, extent of intraspecific competition, risk of juvenile mortality, diurnal pattern of activity, arboreality, and rainfall pattern. Even when controlling for body size and phylogeny, we observe a slow-fast continuum in metatherian mammals. Some parameters involved are different from those identified by studies of eutherians, but the underlying relationships among longevity, fecundity, and age at maturity persist. We also show that overall variation in a key life history variable, reproductive output (measured by annual reproductive rate and litter size), is significantly related to variation in type of diet, with a foliage-rich diet being associated with low fecundity. This is interesting because, although ecological correlations have been found within some eutherian subgroups, modern comparative approaches have failed to reveal robust ecological correlates of overall life history diversity in eutherians. Copyright ESA. All rights reserved.

Convergent maternal care strategies in ungulates and macropods

Mammals show extensive interspecific variation in the form of maternal care. Among ungulates, there is a dichotomy between species in which offspring follow the mother (following strategy) versus species in which offspring remain concealed (hiding strategy). Here we reveal that the same dichotomy exists among macropods (kangaroos, wallabies and allies). We test three traditional adaptive explanations and one new life history hypothesis. and find very similar patterns among both ungulates and macropods. The three traditional explanations that we tested were that a ''following'' strategy is associated with (1) open habitat, (2) large mothers, and (3) gregariousness. Our new life-history hypothesis is that a following strategy'' is associated with delayed weaning, and thus with the slow end of the slow-fast mammalian life-history continuum, because offspring devote resources to locomotion rather than rapid growth. Our comparative test strongly supports the habitat structure hypothesis and provides some support for this new delayed weaning hypothesis for both ungulates and macropods. We propose that sedentary young in closed habitats benefit energetically by having milk brought to them. In open habitats, predation pressure will select against hiding. Followers will suffer slower growth to independence. Taken together, therefore, our results provide the first quantitative evidence that macropods and ungulates are convergent with respect to interspecific variation in maternal care strategy. In both clades, differences between species in the form of parental care are due to a similar interaction between habitat, social behavior, and life history.

Species richness in agamid lizards: chance, body size, sexual selection or ecology?

Why does species richness vary so greatly across lineages? Traditionally, variation in species richness has been attributed to deterministic processes, although it is equally plausible that it may result from purely stochastic processes. We show that, based on the best available phylogenetic hypothesis, the pattern of cladogenesis among agamid lizards is not consistent with a random model, with some lineages having more species, and others fewer species, than expected by chance. We then use phylogenetic comparative methods to test six types of deterministic explanation for variation in species richness: body size, life history, sexual selection, ecological generalism, range size and latitude. Of eight variables we tested, only sexual size dimorphism and sexual dichromatism predicted species richness. Increases in species richness are associated with increases in sexual dichromatism but reductions in sexual size dimorphism. Consistent with recent comparative studies, we find no evidence that species richness is associated with small body size or high fecundity. Equally, we find no evidence that species richness covaries with ecological generalism, latitude or range size.

Macroecologia de carnívoros do Novo Mundo (Mammalia): envelopes de restrição e análise de padrões filogenéticos

The relationship between body size and geographic range was analyzed for 70 species of terrestrial Carnivora ("fissipeds") of the New World, after the control of phylogenetic patterns in the data using phylogenetic eigenvector regression. The analysis from EcoSim software showed that the variables are related as a triangular envelope. Phylogenetic patterns in data were detected by means of phylogenetic correlograms, and 200 simulations of the phenotypic evolution were also performed over the phylogeny. For body size, the simulations suggested a non-linear relationship for the evolution of this character along the phylogeny. For geographic range size, the correlogram showed no phylogenetic patterns. A phylogenetic eigenvector regression was performed on original data and on data simulated under Ornstein-Uhlenbeck process. Since both characters did not evolve under a simple Brownian motion process, the Type I errors should be around 10%, compatible with other methods to analyze correlated evolution. The significant correlation of the original data (r = 0.38; P < 0.05), as well as the triangular envelope, then indicate ecological and adaptive processes connecting the two variables, such as those proposed in minimum viable population models.

Analyse de la corrélation conditionnelle dérivée de la coévolution d’un système de trois gènes par un modèle du maximum de vraisemblance

Les gènes codant pour des protéines peuvent souvent être regroupés et intégrés en modules fonctionnels par rapport à un organelle. Ces modules peuvent avoir des composantes qui suivent une évolution corrélée pouvant être conditionnelle à un phénotype donné. Les gènes liés à la motilité possèdent cette caractéristique, car ils se suivent en cascade en réponse à des stimuli extérieurs. L’hyperthermophilie, d’autre part, est interreliée à la reverse gyrase, cependant aucun autre élément qui pourrait y être associé avec certitude n’est connu. Ceci peut être dû à un déplacement de gènes non orthologues encore non résolu. En utilisant une approche bio-informatique, une modélisation mathématique d’évolution conditionnelle corrélée pour trois gènes a été développée et appliquée sur des profils phylétiques d’archaea. Ceci a permis d’établir des théories quant à la fonction potentielle du gène du flagelle FlaD/E ainsi que l’histoire évolutive des gènes lui étant liés et ayant contribué à sa formation. De plus, une histoire évolutive théorique a été établie pour une ligase liée à l’hyperthermophilie.

Estimation des corrélations phylogénétiques entre paramètres d'évolution moléculaire et Traits d'histoire de vie

Depuis quelques années, l'évolution moléculaire cherche à caractériser les variations et l'intensité de la sélection grâce au rapport entre taux de substitution synonyme et taux de substitution non-synonyme (dN/dS). Cette mesure, dN/dS, a permis d'étudier l'histoire de la variation de l'intensité de la sélection au cours du temps ou de détecter des épisodes de la sélection positive. Les liens entre sélection et variation de taille efficace interfèrent cependant dans ces mesures. Les méthodes comparatives, quant a elle, permettent de mesurer les corrélations entre caractères quantitatifs le long d'une phylogénie. Elles sont également utilisées pour tester des hypothèses sur l'évolution corrélée des traits d'histoire de vie, mais pour être employées pour étudier les corrélations entre traits d'histoire de vie, masse, taux de substitution ou dN/dS. Nous proposons ici une approche combinant une méthode comparative basée sur le principe des contrastes indépendants et un modèle d'évolution moléculaire, dans un cadre probabiliste Bayésien. Intégrant, le long d'une phylogénie, sur les reconstructions ancestrales des traits et et de dN/dS nous estimons les covariances entre traits ainsi qu'entre traits et paramètres du modèle d'évolution moléculaire. Un modèle hiérarchique, a été implémenté dans le cadre du logiciel coevol, publié au cours de cette maitrise. Ce modèle permet l'analyse simultané de plusieurs gènes sans perdre la puissance donnée par l'ensemble de séquences. Un travail deparallélisation des calculs donne la liberté d'augmenter la taille du modèle jusqu'à l'échelle du génome. Nous étudions ici les placentaires, pour lesquels beaucoup de génomes complets et de mesures phénotypiques sont disponibles. À la lumière des théories sur les traits d'histoire de vie, notre méthode devrait permettre de caractériser l'implication de groupes de gènes dans les processus biologique liés aux phénotypes étudiés.

Assembly rules for protein networks derived from phylogenetic-statistical analysis of whole genomes

Background: We report an analysis of a protein network of functionally linked proteins, identified from a phylogenetic statistical analysis of complete eukaryotic genomes. Phylogenetic methods identify pairs of proteins that co-evolve on a phylogenetic tree, and have been shown to have a high probability of correctly identifying known functional links. Results: The eukaryotic correlated evolution network we derive displays the familiar power law scaling of connectivity. We introduce the use of explicit phylogenetic methods to reconstruct the ancestral presence or absence of proteins at the interior nodes of a phylogeny of eukaryote species. We find that the connectivity distribution of proteins at the point they arise on the tree and join the network follows a power law, as does the connectivity distribution of proteins at the time they are lost from the network. Proteins resident in the network acquire connections over time, but we find no evidence that 'preferential attachment' - the phenomenon of newly acquired connections in the network being more likely to be made to proteins with large numbers of connections - influences the network structure. We derive a 'variable rate of attachment' model in which proteins vary in their propensity to form network interactions independently of how many connections they have or of the total number of connections in the network, and show how this model can produce apparent power-law scaling without preferential attachment. Conclusion: A few simple rules can explain the topological structure and evolutionary changes to protein-interaction networks: most change is concentrated in satellite proteins of low connectivity and small phenotypic effect, and proteins differ in their propensity to form attachments. Given these rules of assembly, power law scaled networks naturally emerge from simple principles of selection, yielding protein interaction networks that retain a high-degree of robustness on short time scales and evolvability on longer evolutionary time scales.

Loss of male secondary sexual structures in allopatry in the Neotropical butterfly genus Arcas (Lycaenidae: Iheclinae: Eumaeini)

Male secondary sexual characters in Lepidoptera may be present or absent in species that otherwise appear to be closely related, an observation that has led to differences of opinion over the taxonomic usefulness of these structures above the species level. An evolutionary issue raised by this debate is whether male secondary sexual characters (1) can be regained after being lost evolutionarily, (2) are not lost after being evolved, or (3) are 'switched on and off' by genes that regulate development. A second evolutionary issue is the conditions under which male secondary sexual characters might be lost or gained evolutionarily. Because these structures are thought to promote species recognition, theory predicts evolutionary losses to be most likely in allopatry; evolutionary gains to be most likely during the process of secondarily establishing sympatry or during sympatric speciation. We updated the species-level taxonomy of the brilliant emerald winged Neotropical lycaenid butterfly genus Arcas and performed an analysis of phylogenetic relations among species to assess these evolutionary issues. We morphologically detail a scent pouch on the ventral hindwing of Areas and report that six species possess the pouch with androconia, one possesses the pouch without androconia, and the remaining two species have neither pouch nor androconia. In addition, eight Areas species have a morphologically species-specific male forewing scent pad, and one lacks a scent pad. This variation appears to be the result of three evolutionary losses and no gains of male secondary sexual organs. The four Areas species lacking a scent pouch or a scent pad are allopatric with their closest phylogenetic relatives while four of five with both of these structures are sympatric. Although Arcas is a small genus, these results are significantly more extreme than predicted by chance. For taxonomy, this study provides a rationale for the evolutionary loss of male secondary sexual structures and suggests that their absence, but itself, does not indicate a lack of relationship above the species level.