Loss of sexual recombination and segregation is associated with increased diversification in evening primroses

The evolution of multi-gene families and metabolic pathways in the evening primroses (Oenothera: Onagraceae): A comparative transcriptomics approach

The plant genus Oenothera has played an important role in the study of plant evolution of genomes and plant defense and reproduction. Here, we build on the 1kp transcriptomic dataset by creating 44 new transcriptomes and analyzing a total of 63 transcriptomes to present a large-scale comparative study across 29 Oenothera species. Our dataset included 30.4 million reads per individual and 2.3 million transcripts on average. We used this transcriptome resource to examine genome-wide evolutionary patterns and functional diversification by searching for orthologous genes and performed gene family evolution analysis. We found wide heterogeneity in gene family evolution across the genus, with section Oenothera exhibiting the most pronounced evolutionary changes. Overall, more significant gene family expansions occurred than contractions. We also analyzed the molecular evolution of phenolic metabolism by retrieving proteins annotated for phenolic enzymatic complexes. We identified 1,568 phenolic genes arranged into 83 multigene families that varied widely across the genus. All taxa experienced rapid phenolic evolution (fast rate of genomic turnover) involving 33 gene families, which exhibited large expansions, gaining about 2-fold more genes than they lost. Upstream enzymes phenylalanine ammonia-lyase (PAL) and 4-coumaroyl: CoA ligase (4CL) accounted for most of the significant expansions and contractions. Our results suggest that adaptive and neutral evolutionary processes have contributed to Oenothera diversification and rapid gene family evolution.

No evidence for accumulation of deleterious mutations and fitness degradation in clonal fish hybrids: Abandoning sex without regrets

Despite its inherent costs, sexual reproduction is ubiquitous in nature, and the mechanisms to protect it from a competitive displacement by asexuality remain unclear. Popular mutation-based explanations, like the Muller's ratchet and the Kondrashov's hatchet, assume that purifying selection may not halt the accumulation of deleterious mutations in the nonrecombining genomes, ultimately leading to their degeneration. However, empirical evidence is scarce and it remains particularly unclear whether mutational degradation proceeds fast enough to ensure the decay of clonal organisms and to prevent them from outcompeting their sexual counterparts. To test this hypothesis, we jointly analysed the exome sequences and the fitness-related phenotypic traits of the sexually reproducing fish species and their clonal hybrids, whose evolutionary ages ranged from F1 generations to 300 ky. As expected, mutations tended to accumulate in the clonal genomes in a time-dependent manner. However, contrary to the predictions, we found no trend towards increased nonsynonymity of mutations acquired by clones, nor higher radicality of their amino acid substitutions. Moreover, there was no evidence for fitness degeneration in the old clones compared with that in the younger ones. In summary, although an efficacy of purifying selection may still be reduced in the asexual genomes, our data indicate that its efficiency is not drastically decreased. Even the oldest investigated clone was found to be too young to suffer fitness consequences from a mutation accumulation. This suggests that mechanisms other than mutation accumulation may be needed to explain the competitive advantage of sex in the short term.

Multi-gene phylogeny of Tetrahymena refreshed with three new histophagous species invading freshwater planarians

Planarians represent an insufficiently explored group of aquatic invertebrates that might serve as hosts of histophagous ciliates belonging to the hymenostome genus Tetrahymena. During our extensive research on freshwater planarians, parasitic tetrahymenas were detected in two of the eight planarian species investigated, namely, in Dugesia gonocephala and Girardia tigrina. Using the 16S and 18S rRNA genes as well as the barcoding cytochrome oxidase subunit I, one ciliate species was identified as T. scolopax and three species were recognized as new forms: T. acanthophora, T. dugesiae, and T. nigricans. Thus, 25% of the examined planarian taxa are positive for Tetrahymena species and three of them represent new taxa, indicating a large undescribed ciliate diversity in freshwater planarians. According to phylogenetic analyses, histophagous tetrahymenas show a low phylogenetic host specificity. Although T. acanthophora, T. dugesiae, and T. scolopax clustered together within the "borealis" clade, the former species has been detected exclusively in G. tigrina, while the two latter species only in D. gonocephala. Tetrahymena nigricans, which has been isolated only from G. tigrina, was classified within the "paravorax" clade along with T. glochidiophila which feeds on glochidia. The present phylogenetic reconstruction of ancestral life strategies suggested that the last common ancestor of the family Tetrahymenidae was free-living, unlike the progenitor of the subclass Hymenostomatia which was very likely parasitic. Consequently, there were at least seven independent shifts back to parasitism/histophagy within Tetrahymena: one each in the "paravorax" and "australis" clades and at least five transfers back to parasitism in the "borealis" clade.

Species Selection Regime and Phylogenetic Tree Shape

Species selection, the effect of heritable traits in generating between-lineage diversification rate differences, provides a valuable conceptual framework for understanding the relationship between traits, diversification, and phylogenetic tree shape. An important challenge, however, is that the nature of real diversification landscapes—curves or surfaces which describe the propensity of species-level lineages to diversify as a function of one or more traits—remains poorly understood. Here, we present a novel, time-stratified extension of the QuaSSE model in which speciation/extinction rate is specified as a static or temporally shifting Gaussian or skewed-Gaussian function of the diversification trait. We then use simulations to show that the generally imbalanced nature of real phylogenetic trees, as well as their generally greater than expected frequency of deep branching events, are typical outcomes when diversification is treated as a dynamic, trait-dependent process. Focusing on four basic models (Gaussian-speciation with and without background extinction; skewed-speciation; Gaussian-extinction), we also show that particular features of the species selection regime produce distinct tree shape signatures and that, consequently, a combination of tree shape metrics has the potential to reveal the species selection regime under which a particular lineage diversified. We evaluate this idea empirically by comparing the phylogenetic trees of plant lineages diversifying within climatically and geologically stable environments of the Greater Cape Floristic Region, with those of lineages diversifying in environments that have experienced major change through the Late Miocene-Pliocene. Consistent with our expectations, the trees of lineages diversifying in a dynamic context are less balanced, show a greater concentration of branching events close to the present, and display stronger diversification rate-trait correlations. We suggest that species selection plays an important role in shaping phylogenetic trees but recognize the need for an explicit probabilistic framework within which to assess the likelihoods of alternative diversification scenarios as explanations of a particular tree shape. [Cape flora; diversification landscape; environmental change; gamma statistic; species selection; time-stratified QuaSSE model; trait-dependent diversification; tree imbalance.]

Hybridization and a loss of sex shape genome-wide diversity and the origin of species in the evening primroses (Oenothera, Onagraceae)

Hybridization is thought to promote speciation in at least two ways - by fixation of heterozygosity from diploid progenitors in allopolyploids, and by generation of transgressive phenotypes and shifting fitness optima during homoploid hybrid speciation. While recent studies support a hybrid origin for a growing number of species, the extent to which hybrid origins shape patterns of diversity in asexual species remains underexplored. Here we employed transcriptome sequencing and population genomic analysis to describe patterns of genomic variation in the 13 species belonging to Oenothera subsection Oenothera. Eight of these species are functionally asexual and arose by hybrid speciation from parents spanning a range of phylogenetic divergence. We showed that genomic divergence between parents has been retained as heterozygosity in functionally asexual species, and that genome-wide levels of heterozygosity in these asexuals largely reflects the divergence of parental haplotypes coupled with a breakdown in recombination and segregation across the genome. These results show that divergence between parental species and loss of sex in hybrids shape patterns of whole-genome diversity and the origin of asexual species.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

Estimating Age-Dependent Extinction: Contrasting Evidence from Fossils and Phylogenies

The estimation of diversification rates is one of the most vividly debated topics in modern systematics, with considerable controversy surrounding the power of phylogenetic and fossil-based approaches in estimating extinction. Van Valen's seminal work from 1973 proposed the "Law of constant extinction," which states that the probability of extinction of taxa is not dependent on their age. This assumption of age-independent extinction has prevailed for decades with its assessment based on survivorship curves, which, however, do not directly account for the incompleteness of the fossil record, and have rarely been applied at the species level. Here, we present a Bayesian framework to estimate extinction rates from the fossil record accounting for age-dependent extinction (ADE). Our approach, unlike previous implementations, explicitly models unobserved species and accounts for the effects of fossil preservation on the observed longevity of sampled lineages. We assess the performance and robustness of our method through extensive simulations and apply it to a fossil data set of terrestrial Carnivora spanning the past 40 myr. We find strong evidence of ADE, as we detect the extinction rate to be highest in young species and declining with increasing species age. For comparison, we apply a recently developed analogous ADE model to a dated phylogeny of extant Carnivora. Although the phylogeny-based analysis also infers ADE, it indicates that the extinction rate, instead, increases with increasing taxon age. The estimated mean species longevity also differs substantially, with the fossil-based analyses estimating 2.0 myr, in contrast to 9.8 myr derived from the phylogeny-based inference. Scrutinizing these discrepancies, we find that both fossil and phylogeny-based ADE models are prone to high error rates when speciation and extinction rates increase or decrease through time. However, analyses of simulated and empirical data show that fossil-based inferences are more robust. This study shows that an accurate estimation of ADE from incomplete fossil data is possible when the effects of preservation are jointly modeled, thus allowing for a reassessment of Van Valen's model as a general rule in macroevolution.

The evolution of floral sonication, a pollen foraging behavior used by bees (Anthophila)

Over 22,000 species of biotically pollinated flowering plants, including some major agricultural crops, depend primarily on bees capable of floral sonication for pollination services. The ability to sonicate ("buzz") flowers is widespread in bees but not ubiquitous. Despite the prevalence of this pollinator behavior and its importance to natural and agricultural systems, the evolutionary history of floral sonication in bees has not been previously studied. Here, we reconstruct the evolutionary history of floral sonication in bees by generating a time-calibrated phylogeny and reconstructing ancestral states for this pollen extraction behavior. We also test the hypothesis that the ability to sonicate flowers and thereby efficiently access pollen from a diverse assemblage of plant species, led to increased diversification among sonicating bee taxa. We find that floral sonication evolved on average 45 times within bees, possibly first during the Early Cretaceous (100-145 million years ago) in the common ancestor of bees. We find that sonicating lineages are significantly more species rich than nonsonicating sister lineages when comparing sister clades, but a probabilistic structured rate permutation on phylogenies approach failed to support the hypothesis that floral sonication is a key driver of bee diversification. This study provides the evolutionary framework needed to further study how floral sonication by bees may have facilitated the spread and common evolution of angiosperm species with poricidal floral morphology.

Repeated evolution and reversibility of self-fertilization in the volvocine green algae

Outcrossing and self-fertilization are fundamental strategies of sexual reproduction, each with different evolutionary costs and benefits. Self-fertilization is thought to be an evolutionary "dead-end" strategy, beneficial in the short term but costly in the long term, resulting in self-fertilizing species that occupy only the tips of phylogenetic trees. Here, we use volvocine green algae to investigate the evolution of self-fertilization. We use ancestral-state reconstructions to show that self-fertilization has repeatedly evolved from outcrossing ancestors and that multiple reversals from selfing to outcrossing have occurred. We use three phylogenetic metrics to show that self-fertilization is not restricted to the tips of the phylogenetic tree, a finding inconsistent with the view of self-fertilization as a dead-end strategy. We also find no evidence for higher extinction rates or lower speciation rates in selfing lineages. We find that self-fertilizing species have significantly larger colonies than outcrossing species, suggesting the benefits of selfing may counteract the costs of increased size. We speculate that our macroevolutionary results on self-fertilization (i.e., non-tippy distribution, no decreased diversification rates) may be explained by the haploid-dominant life cycle that occurs in volvocine algae, which may alter the costs and benefits of selfing.

Fitness consequences of occasional outcrossing in a functionally asexual plant (Oenothera biennis)

Many clonal organisms occasionally outcross, but the long-term consequences of such infrequent events are often unknown. During five years, representing three to five plant generations, we followed 16 experimental field populations of the forb, Oenothera biennis, originally planted with the same 18 original genotypes. Oenothera biennis usually self fertilizes, which, due to its genetic system (permanent translocation heterozygosity), results in seeds that are clones of the maternal plant. However, rare outcrossing produces genetically novel offspring (but without recombination or increased heterozygosity). We sought to understand whether novel genotypes produced through natural outcrossing had greater fecundity or different multigenerational dynamics compared to our original genotypes. We further assessed whether any differences in fitness or abundances through time between original and novel genotypes were exaggerated in the presence vs. absence of insect herbivores. Over the course of the experiment, we genotyped >12,500 plants using microsatellite DNA markers to identify and track the frequency of specific genotypes and estimated fecundity on a subset (>3,000) of plants. The effective outcrossing rate was 7.3% in the first year and ultimately 50% of the plants were of outcrossed origin by the final year of the experiment. Lifetime fruit production per plant was on average 32% higher across all novel genotypes produced via outcrossing compared to the original genotypes, and this fecundity advantage was significantly enhanced in populations lacking herbivores. Among 43 novel genotypes that were abundant enough to phenotype with replication, plants produced nearly 30% more fruits than the average of their specific two parental genotypes, and marginally more fruits (8%) than their most fecund parent. Mean per capita fecundity of novel genotypes predicted their relative frequencies at the end of the experiment. Novel genotypes increased more dramatically in herbivore-present compared to suppressed populations (45% vs. 27% of all plants), countering the increased competition from dandelions (Taraxacum officinale) that resulted from herbivore suppression. Increased interspecific competition likely also lead to the lower realized fitness of novel vs. original genotypes in herbivore-suppressed populations. These results demonstrate that rare outcrossing and the generation of novel genotypes can create high-fecundity progeny, with the biotic environment influencing the dynamical outcome of such advantages.

Low recombination rates in sexual species and sex-asex transitions

In most sexual, diploid eukaryotes, at least one crossover occurs between each pair of homologous chromosomes during meiosis, presumably in order to ensure proper segregation. Well-known exceptions to this rule are species in which one sex does not recombine and specific chromosomes lacking crossover. We review other possible exceptions, including species with chromosome maps of less than 50 cM in one or both sexes. We discuss the idea that low recombination rates may favour sex-asex transitions, or, alternatively may be a consequence of it. We then show that a yet undescribed species of brine shrimp Artemia from Kazakhstan (A sp. Kazakhstan), the closest known relative of the asexual Artemia parthenogenetica, has one of the shortest genetic linkage maps known. Based on a family of 42 individuals and 589 RAD markers, we find that many linkage groups are considerably shorter than 50 cM, suggesting either no obligate crossover or crossovers concentrated at terminal positions with little effect on recombination. We contrast these findings with the published map of the more distantly related sexual congener, A. franciscana, and conclude that the study of recombination in non-model systems is important to understand the evolutionary causes and consequences of recombination.This article is part of the themed issue 'Evolutionary causes and consequences of recombination rate variation in sexual organisms'.

Comparative tests of the role of dewlap size in Anolis lizard speciation

Phenotypic traits may be linked to speciation in two distinct ways: character values may influence the rate of speciation or diversification in the trait may be associated with speciation events. Traits involved in signal transmission, such as the dewlap of Anolis lizards, are often involved in the speciation process. The dewlap is an important visual signal with roles in species recognition and sexual selection, and dewlaps vary among species in relative size as well as colour and pattern. We compile a dataset of relative dewlap size digitized from photographs of 184 anole species from across the genus' geographical range. We use phylogenetic comparative methods to test two hypotheses: that larger dewlaps are associated with higher speciation rates, and that relative dewlap area diversifies according to a speciational model of evolution. We find no evidence of trait-dependent speciation, indicating that larger signals do not enhance any role the dewlap has in promoting speciation. Instead, we find a signal of mixed speciational and gradual trait evolution, with a particularly strong signal of speciational change in the dewlaps of mainland lineages. This indicates that dewlap size diversifies in association with the speciation process, suggesting that divergent selection may play a role in the macroevolution of this signalling trait.

Diversification dynamics of rhynchostomatian ciliates: the impact of seven intrinsic traits on speciation and extinction in a microbial group

Ciliates are a suitable microbial model to investigate trait-dependent diversification because of their comparatively complex morphology and high diversity. We examined the impact of seven intrinsic traits on speciation, extinction, and net-diversification of rhynchostomatians, a group of comparatively large, predatory ciliates with proboscis carrying a dorsal brush (sensoric structure) and toxicysts (organelles used to kill the prey). Bayesian estimates under the binary-state speciation and extinction model indicate that two types of extrusomes and two-rowed dorsal brush raise diversification through decreasing extinction. On the other hand, the higher number of contractile vacuoles and their dorsal location likely increase diversification via elevating speciation rate. Particular nuclear characteristics, however, do not significantly differ in their diversification rates and hence lineages with various macronuclear patterns and number of micronuclei have similar probabilities to generate new species. Likelihood-based quantitative state diversification analyses suggest that rhynchostomatians conform to Cope’s rule in that their diversity linearly grows with increasing body length and relative length of the proboscis. Comparison with other litostomatean ciliates indicates that rhynchostomatians are not among the cladogenically most successful lineages and their survival over several hundred million years could be associated with their comparatively large and complex bodies that reduce the risk of extinction.

Aging asexual lineages and the evolutionary maintenance of sex

Finite populations of asexual and highly selfing species suffer from a reduced efficacy of selection. Such populations are thought to decline in fitness over time due to accumulating slightly deleterious mutations or failing to adapt to changing conditions. These within-population processes that lead nonrecombining species to extinction may help maintain sex and outcrossing through species level selection. Although inefficient selection is proposed to elevate extinction rates over time, previous models of species selection for sex assumed constant diversification rates. For sex to persist, classic models require that asexual species diversify at rates lower than sexual species; the validity of this requirement is questionable, both conceptually and empirically. We extend past models by allowing asexual lineages to decline in diversification rates as they age, that is nonrecombining lineages "senesce" in diversification rates. At equilibrium, senescing diversification rates maintain sex even when asexual lineages, at young ages, diversify faster than their sexual progenitors. In such cases, the age distribution of asexual lineages contains a peak at intermediate values rather than showing the exponential decline predicted by the classic model. Coexistence requires only that the average rate of diversification in asexuals be lower than that of sexuals.

Assortative mating and self-fertilization differ in their contributions to reinforcement, cascade speciation, and diversification

Cascade speciation and reinforcement can evolve rapidly when traits are pleiotropic and act as both signal/cue in nonrandom mating. Here, we examine the contribution of two key traits-assortative mating and self-fertilization-to reinforcement and (by extension) cascade speciation. First, using a population genetic model of reinforcement we find that both assortative mating and self-fertilization can make independent contributions to increased reproductive isolation, consistent with reinforcement. Self-fertilization primarily evolves due to its 2-fold transmission advantage when inbreeding depression (d) is lower (d < 0.45) but evolves as a function of the cost of hybridization under higher inbreeding depression (0.45 < d < 0.48). When both traits can evolve simultaneously, increased self-fertilization often prohibits the evolution of assortative mating. We infer that, under specific conditions, mating system transitions are more likely to lead to increased reproductive isolation and initiate cascade speciation, than assortative mating. Based on the results of our simulations, we hypothesized that transitions to self-fertilization could contribute to clade-wide diversification if reinforcement or cascade speciation is common. We tested this hypothesis with comparative data from two different groups. Consistent with our hypothesis, there was a trend towards uniparental reproduction being associated with increased diversification rate in the Nematode phylum. For the plant genus Mimulus, however, self-fertilization was associated with reduced diversification. Reinforcement driving speciation via transitions to self-fertilization might be short lived or unsustainable across macroevolutionary scales in some systems (some plants), but not others (such as nematodes), potentially due to differences in susceptibility to inbreeding depression and/or the ability to transition between reproductive modes.

Species selection and random drift in macroevolution

Species selection resulting from trait-dependent speciation and extinction is increasingly recognized as an important mechanism of phenotypic macroevolution. However, the recent bloom in statistical methods quantifying this process faces a scarcity of dynamical theory for their interpretation, notably regarding the relative contributions of deterministic versus stochastic evolutionary forces. I use simple diffusion approximations of birth-death processes to investigate how the expected and random components of macroevolutionary change depend on phenotype-dependent speciation and extinction rates, as can be estimated empirically. I show that the species selection coefficient for a binary trait, and selection differential for a quantitative trait, depend not only on differences in net diversification rates (speciation minus extinction), but also on differences in species turnover rates (speciation plus extinction), especially in small clades. The randomness in speciation and extinction events also produces a species-level equivalent to random genetic drift, which is stronger for higher turnover rates. I then show how microevolutionary processes including mutation, organismic selection, and random genetic drift cause state transitions at the species level, allowing comparison of evolutionary forces across levels. A key parameter that would be needed to apply this theory is the distribution and rate of origination of new optimum phenotypes along a phylogeny.

Widespread flower color convergence in Solanaceae via alternate biochemical pathways

Phenotypic convergence is rampant throughout the tree of life. While recent studies have made significant progress in ascertaining the proximate mechanisms underlying convergent phenotypes, less is known about the frequency and predictability with which convergent phenotypes arise via the same or multiple pathways at the macroevolutionary scale. We investigated the proximate causes and evolutionary patterns of red flower color in the tomato family, Solanaceae, using large-scale data mining and new sequence data to reconstruct a megaphylogeny of 1341 species. We then combined spectral and anatomical data to assess how many times red flowers have evolved, the relative contribution of different pathways to independent origins of red, and whether the underlying pathway is predicted by phylogenetic relatedness. We estimated at least 30 relatively recent origins of red flowers using anthocyanins, carotenoids, or a dual production of both pigments, with significant phylogenetic signal in the use of anthocyanins and dual production, indicating that closely related red-flowered species tend to employ the same mechanism for coloration. Our study is the first to test whether developmental pathways exhibit phylogenetic signal and implies that historical contingency strongly influences the evolution of new phenotypes.

Quantifying Age-dependent Extinction from Species Phylogenies

Several ecological factors that could play into species extinction are expected to correlate with species age, i.e., time elapsed since the species arose by speciation. To date, however, statistical tools to incorporate species age into likelihood-based phylogenetic inference have been lacking. We present here a computational framework to quantify age-dependent extinction through maximum likelihood parameter estimation based on phylogenetic trees, assuming species lifetimes are gamma distributed. Testing on simulated trees shows that neglecting age dependence can lead to biased estimates of key macroevolutionary parameters. We then apply this method to two real data sets, namely a complete phylogeny of birds (class Aves) and a clade of self-compatible and -incompatible nightshades (Solanaceae), gaining initial insights into the extent to which age-dependent extinction may help explain macroevolutionary patterns. Our methods have been added to the R package TreePar.

Multiple independent origins of auto-pollination in tropical orchids (Bulbophyllum) in light of the hypothesis of selfing as an evolutionary dead end

BACKGROUND

The transition from outcrossing to selfing has long been portrayed as an 'evolutionary dead end' because, first, reversals are unlikely and, second, selfing lineages suffer from higher rates of extinction owing to a reduced potential for adaptation and the accumulation of deleterious mutations. We tested these two predictions in a clade of Madagascan Bulbophyllum orchids (30 spp.), including eight species where auto-pollinating morphs (i.e., selfers, without a 'rostellum') co-exist with their pollinator-dependent conspecifics (i.e., outcrossers, possessing a rostellum). Specifically, we addressed this issue on the basis of a time-calibrated phylogeny by means of ancestral character reconstructions and within the state-dependent evolution framework of BiSSE (Binary State Speciation and Extinction), which allowed jointly estimating rates of transition, speciation, and extinction between outcrossing and selfing.

RESULTS

The eight species capable of selfing occurred in scattered positions across the phylogeny, with two likely originating in the Pliocene (ca. 4.4-3.1 Ma), one in the Early Pleistocene (ca. 2.4 Ma), and five since the mid-Pleistocene (ca. ≤ 1.3 Ma). We infer that this scattered phylogenetic distribution of selfing is best described by models including up to eight independent outcrossing-to-selfing transitions and very low rates of speciation (and either moderate or zero rates of extinction) associated with selfing.

CONCLUSIONS

The frequent and irreversible outcrossing-to-selfing transitions in Madagascan Bulbophyllum are clearly congruent with the first prediction of the dead end hypothesis. The inability of our study to conclusively reject or support the likewise predicted higher extinction rate in selfing lineages might be explained by a combination of methodological limitations (low statistical power of our BiSSE approach to reliably estimate extinction in small-sized trees) and evolutionary processes (insufficient time elapsed for selfers to go extinct). We suggest that, in these tropical orchids, a simple genetic basis of selfing (via loss of the 'rostellum') is needed to explain the strikingly recurrent transitions to selfing, perhaps reflecting rapid response to parallel and novel selective environments over Late Quaternary (≤ 1.3 Ma) time scales.

No evidence that sex and transposable elements drive genome size variation in evening primroses

Genome size varies dramatically across species, but despite an abundance of attention there is little agreement on the relative contributions of selective and neutral processes in governing this variation. The rate of sex can potentially play an important role in genome size evolution because of its effect on the efficacy of selection and transmission of transposable elements (TEs). Here, we used a phylogenetic comparative approach and whole genome sequencing to investigate the contribution of sex and TE content to genome size variation in the evening primrose (Oenothera) genus. We determined genome size using flow cytometry for 30 species that vary in genetic system and find that variation in sexual/asexual reproduction cannot explain the almost twofold variation in genome size. Moreover, using whole genome sequences of three species of varying genome sizes and reproductive system, we found that genome size was not associated with TE abundance; instead the larger genomes had a higher abundance of simple sequence repeats. Although it has long been clear that sexual reproduction may affect various aspects of genome evolution in general and TE evolution in particular, it does not appear to have played a major role in genome size evolution in the evening primroses.

Effects of functionally asexual reproduction on quantitative genetic variation in the evening primroses (Oenothera, Onagraceae)

PREMISE OF THE STUDY

It has long been predicted that a loss of sexual reproduction leads to decreased heritable variation within populations and increased differentiation between populations. Despite an abundance of theory, there are few empirical tests of how sex affects genetic variation in phenotypic traits, especially for plants. Here we test whether repeated losses of two critical components of sex (recombination and segregation) in the evening primroses (Oenothera L., Onagraceae) affect quantitative genetic variation within and between populations.

METHODS

We sampled multiple genetic families from 3-5 populations from each of eight Oenothera species, which represented four independent transitions between sexual reproduction and a functionally asexual genetic system called "permanent translocation heterozygosity." We used quantitative genetics methods to partition genetic variation within and between populations for eight plant traits related to growth, leaf physiology, flowering, and resistance to herbivores.

KEY RESULTS

Heritability was, on average, 74% higher in sexual Oenothera populations than in functionally asexual populations, with plant growth rate, specific leaf area, and the percentage of leaf water content showing the strongest differences. By contrast, genetic differentiation among populations was 2.8× higher in functionally asexual vs. sexual Oenothera species. This difference was particularly strong for specific leaf area. Sexual populations tended to exhibit higher genetic correlations among traits, but this difference was weakly supported.

CONCLUSIONS

These results support the prediction that sexual reproduction maintains higher genetic variation within populations, which may facilitate adaptive evolution. We also found partial support for the prediction that a loss of sex leads to greater population differentiation, which may elevate speciation rates.

Taxonomy and Biogeography of Apomixis in Angiosperms and Associated Biodiversity Characteristics

Apomixis in angiosperms is asexual reproduction from seed. Its importance to angiospermous evolution and biodiversity has been difficult to assess mainly because of insufficient taxonomic documentation. Thus, we assembled literature reporting apomixis occurrences among angiosperms and transferred the information to an internet database (http://www.apomixis.uni-goettingen.de). We then searched for correlations between apomixis occurrences and well-established measures of taxonomic diversity and biogeography. Apomixis was found to be taxonomically widespread with no clear tendency to specific groups and to occur with sexuality at all taxonomic levels. Adventitious embryony was the most frequent form (148 genera) followed by apospory (110) and diplospory (68). All three forms are phylogenetically scattered, but this scattering is strongly associated with measures of biodiversity. Across apomictic-containing orders and families, numbers of apomict-containing genera were positively correlated with total numbers of genera. In general, apomict-containing orders, families, and subfamilies of Asteraceae, Poaceae, and Orchidaceae were larger, i.e., they possessed more families or genera, than non-apomict-containing orders, families or subfamilies. Furthermore, many apomict-containing genera were found to be highly cosmopolitan. In this respect, 62% occupy multiple geographic zones. Numbers of genera containing sporophytic or gametophytic apomicts decreased from the tropics to the arctic, a trend that parallels general biodiversity. While angiosperms appear to be predisposed to shift from sex to apomixis, there is also evidence of reversions to sexuality. Such reversions may result from genetic or epigenetic destabilization events accompanying hybridization, polyploidy, or other cytogenetic alterations. Because of increased within-plant genetic and genomic heterogeneity, range expansions and diversifications at the species and genus levels may occur more rapidly upon reversion to sexuality. The significantly-enriched representations of apomicts among highly diverse and geographically-extensive taxa, from genera to orders, support this conclusion.

How traits shape trees: new approaches for detecting character state-dependent lineage diversification

Biologists have long sought to understand the processes underlying disparities in clade size across the tree of life and the extent to which such clade size differences can be attributed to the evolution of particular traits. The association of certain character states with species-rich clades suggests that trait evolution can lead to increased diversification, but such a pattern could also arise due other processes, such as directional trait evolution. Recent advances in phylogenetic comparative methods have provided new statistical approaches for distinguishing between these intertwined and potentially confounded macroevolutionary processes. Here, we review the historical development of methods for detecting state-dependent diversification and explore what new methods have revealed about classic examples of traits that affect diversification, including evolutionary dead ends, key innovations and geographic traits. Applications of these methods thus far collectively suggest that trait diversity commonly arises through the complex interplay between transition, speciation and extinction rates and that long hypothesized evolutionary dead ends and key innovations are instead often cases of directional trends in trait evolution.

Macroevolution of plant defenses against herbivores in the evening primroses

Plant species vary greatly in defenses against herbivores, but existing theory has struggled to explain this variation. Here, we test how phylogenetic relatedness, tradeoffs, trait syndromes, and sexual reproduction affect the macroevolution of defense. To examine the macroevolution of defenses, we studied 26 Oenothera (Onagraceae) species, combining chemistry, comparative phylogenetics and experimental assays of resistance against generalist and specialist herbivores. We detected dozens of phenolic metabolites within leaves, including ellagitannins (ETs), flavonoids, and caffeic acid derivatives (CAs). The concentration and composition of phenolics exhibited low to moderate phylogenetic signal. There were clear negative correlations between multiple traits, supporting the prediction of allocation tradeoffs. There were also positively covarying suites of traits, but these suites did not strongly predict resistance to herbivores and thus did not act as defensive syndromes. By contrast, specific metabolites did correlate with the performance of generalist and specialist herbivores. Finally, that repeated losses of sex in Oenothera was associated with the evolution of increased flavonoid diversity and altered phenolic composition. These results show that secondary chemistry has evolved rapidly during the diversification of Oenothera. This evolution has been marked by allocation tradeoffs between traits, some of which are related to herbivore performance. The repeated loss of sex appears also to have constrained the evolution of plant secondary chemistry, which may help to explain variation in defense among plants.

Dioecy is associated with higher diversification rates in flowering plants

In angiosperms, dioecious clades tend to have fewer species than their nondioecious sister clades. This departure from the expected equal species richness in the standard sister clade test has been interpreted as implying that dioecious clades diversify less and has initiated a series of studies suggesting that dioecy might be an 'evolutionary dead end'. However, two of us recently showed that the 'equal species richness' null hypothesis is not valid in the case of derived char acters, such as dioecy, and proposed a new test for sister clade comparisons; preliminary results, using a data set available in the litterature, indicated that dioecious clades migth diversify more than expected. However, it is crucial for this new test to distinguish between ancestral and derived cases of dioecy, a criterion that was not taken into account in the available data set. Here, we present a new data set that was obtained by searching the phylogenetic literature on more than 600 completely dioecious angiosperm genera and identifying 115 sister clade pairs for which dioecy is likely to be derived (including > 50% of the dioecious species). Applying the new sister clade test to this new dataset, we confirm the preliminary result that dioecy is associated with an increased diversification rate, a result that does not support the idea that dioecy is an evolutionary dead end in angiosperms. The traits usually associated with dioecy, that is, an arborescent growth form, abiotic pollination, fleshy fruits or a tropical distribution, do not influence the diversification rate. Rather than a low diversification rate, the observed species richness patterns of dioecious clades seem to be better explained by a low transition rate to dioecy and frequent losses.

Heterostyly accelerates diversification via reduced extinction in primroses

The exceptional species diversity of flowering plants, exceeding that of their sister group more than 250-fold, is especially evident in floral innovations, interactions with pollinators and sexual systems. Multiple theories, emphasizing flower-pollinator interactions, genetic effects of mating systems or high evolvability, predict that floral evolution profoundly affects angiosperm diversification. However, consequences for speciation and extinction dynamics remain poorly understood. Here, we investigate trajectories of species diversification focusing on heterostyly, a remarkable floral syndrome where outcrossing is enforced via cross-compatible floral morphs differing in placement of their respective sexual organs. Heterostyly evolved at least 20 times independently in angiosperms. Using Darwin's model for heterostyly, the primrose family, we show that heterostyly accelerates species diversification via decreasing extinction rates rather than increasing speciation rates, probably owing to avoidance of the negative genetic effects of selfing. However, impact of heterostyly appears to differ over short and long evolutionary time-scales: the accelerating effect of heterostyly on lineage diversification is manifest only over long evolutionary time-scales, whereas recent losses of heterostyly may prompt ephemeral bursts of speciation. Our results suggest that temporal or clade-specific conditions may ultimately determine the net effects of specific traits on patterns of species diversification.

Let us not be unfair to asexuals: their ephemerality may be explained by neutral models without invoking any evolutionary constraints of asexuality

Phylogenetic studies typically demonstrate lower evolutionary ages of clones, relative to their sexual ancestors. This has often been attributed to heightened extinction risk of asexual organisms. We previously criticized such interpretations and demonstrated that the life span of clones is ultimately limited by neutral drift depending on the rate at which new clones are spawned into an asexual community of a finite size. Therefore, it is important to investigate whether the natural rates of such influxes are sufficiently high to account for the relative ephemerality of clones without assuming their increased extinction rate. I applied the neutral clonal turnover model to phylogenies of polyploid asexual ferns and simulated the coalescent trees over a wide range of demographic structures and sampling schemes. On parameterizing the model with biologically relevant estimates of population sizes and plant polyploidization rates, simulated clonal assemblages appeared younger than their sexual counterparts even in the absence of selection against clones. Therefore, differences observed between the ages of sexual and clonal lineages may be explained by the neutral clonal turnover. Researchers should consider the possibility that natural clones may get lost by neutral drift before their fate could eventually be affected by any long-term constraints of asexuality.

Signatures of recombination in clonal lineages of the citrus brown spot pathogen, Alternaria alternata sensu lato

Most Alternaria spp. are considered asexual but recent molecular evolution analyses of Alternaria mating-type genes show that the mating locus is under strong purifying selection, indicating a possible role in sexual reproduction. The objective of this study was to determine the mode of reproduction of an Alternaria alternata sensu lato population causing citrus brown spot in central Florida. Mating type of each isolate was determined, and isolates were sequenced at six putatively unlinked loci. Three genetically distinct subpopulations (SH1, SH4A, and SH4B) were identified using network and Bayesian population structure analyses. Results demonstrate that most subpopulations of A. alternata associated with citrus are clonal but some have the ability to extensively recombine through a cryptic sexual cycle or parasexual cycle. Although isolates were sampled in close physical proximity (≈2,500-m² area), we were able to reject a random mating model using multilocus gametic disequilibrium tests for two subpopulations, SH1 and SH4B, suggesting that these subpopulations were predominantly asexual. However, three recombination events were identified in SH1 and SH4B and localized to individuals of opposite mating type, possibly indicating meiotic recombination. In contrast, in the third subpopulation (SH4A), where only one mating type was present, extensive reticulation was evident in network analyses, and multilocus gametic disequilibrium tests were consistent with recombination. Recombination among isolates of the same mating type suggests that a nonmeiotic mechanism of recombination such as the parasexual cycle may be operating in this subpopulation. The level of gene flow detected among subpopulations does not appear to be sufficient to prevent differentiation, and perhaps future speciation, of these A. alternata subpopulations.

Lineage Selection and the Maintenance of Sex

Sex predominates in eukaryotes, despite its short-term disadvantage when compared to asexuality. Myriad models have suggested that short-term advantages of sex may be sufficient to counterbalance its twofold costs. However, despite decades of experimental work seeking such evidence, no evolutionary mechanism has yet achieved broad recognition as explanation for the maintenance of sex. We explore here, through lineage-selection models, the conditions favouring the maintenance of sex. In the first model, we allowed the rate of transition to asexuality to evolve, to determine whether lineage selection favoured species with the strongest constraints preventing the loss of sex. In the second model, we simulated more explicitly the mechanisms underlying the higher extinction rates of asexual lineages than of their sexual counterparts. We linked extinction rates to the ecological and/or genetic features of lineages, thereby providing a formalisation of the only figure included in Darwin's "The origin of species". Our results reinforce the view that the long-term advantages of sex and lineage selection may provide the most satisfactory explanations for the maintenance of sex in eukaryotes, which is still poorly recognized, and provide figures and a simulation website for training and educational purposes. Short-term benefits may play a role, but it is also essential to take into account the selection of lineages for a thorough understanding of the maintenance of sex.

Speciation in the asexual realm: is the parthenogenetic weevil Naupactus cervinus a complex of species in statu nascendi?

Population genetic theory shows that asexual organisms may evolve into species, which behave as independent evolutionary units. As a result, they form genotypic clusters separated by deep gaps due to geographic isolation and/or divergent selection. Identification of several genetically divergent groups of weevils embodied in the nominal species Naupactus cervinus deserves further study, in order to test if these lineages are evolving independently. In the present paper we tested if the parthenogenetic weevil N. cervinus, native to South America and broadly distributed throughout the world, contains more than one evolutionary unit. For this purpose, we applied three different approaches, a multilocus phylogenetic analysis, the GMYC approach and the K/θ method. We accomplished these analyses through a survey of mitochondrial (COI and COII genes) and nuclear (ITS1 sequence) genetic variation and morphometric analysis in a sample which included individuals from different locations within the native geographic range of N. cervinus. In addition, we compared the divergence accumulated in this species with that in another weevil of the same tribe (Naupactini) showing identical reproductive mode to see if similar levels of morphological variation matches similar levels of genetic divergence. We report the presence of two independent evolutionary units living in sympatry in forest areas. The incongruence between mitochondrial and nuclear datasets analyzed herein reflects incomplete lineage sorting of the nuclear marker and different evolutionary rates between genomes. Ecological divergence driven by natural selection (sympatry) or secondary contact after geographic isolation (allopatry) might explain the deep gaps in mitochondrial phylogenies. Instead, Wolbachia infection was ruled out as a causal factor for such differentiation. We conclude that N. cervinus is probably a species complex with at least two well differentiated lineages that would represent a cluster of species in statu nascendi.

Exploring power and parameter estimation of the BiSSE method for analyzing species diversification

BACKGROUND

There has been a considerable increase in studies investigating rates of diversification and character evolution, with one of the promising techniques being the BiSSE method (binary state speciation and extinction). This study uses simulations under a variety of different sample sizes (number of tips) and asymmetries of rate (speciation, extinction, character change) to determine BiSSE's ability to test hypotheses, and investigate whether the method is susceptible to confounding effects.

RESULTS

We found that the power of the BiSSE method is severely affected by both sample size and high tip ratio bias (one character state dominates among observed tips). Sample size and high tip ratio bias also reduced accuracy and precision of parameter estimation, and resulted in the inability to infer which rate asymmetry caused the excess of a character state. In low tip ratio bias scenarios with appropriate tip sample size, BiSSE accurately estimated the rate asymmetry causing character state excess, avoiding the issue of confounding effects.

CONCLUSIONS

Based on our findings, we recommend that future studies utilizing BiSSE that have fewer than 300 terminals and/or have datasets where high tip ratio bias is observed (i.e., fewer than 10% of species are of one character state) should be extremely cautious with the interpretation of hypothesis testing results.

The Evolutionary Dynamics of Apomixis in Ferns: A Case Study from Polystichoid Ferns

The disparate distribution of apomixis between the major plant lineages is arguably one of the most paradoxical phenomena in plant evolution. Ferns are particularly interesting for addressing this issue because apomixis is more frequent than in any other group of plants. Here, we use a phylogenetic framework to explore some aspects of the evolution of apomixis in ferns and in particular in the polystichoid ferns. Our findings indicate that apomixis evolved several times independently in three different clades of polystichoid ferns. A lineage-wide perspective across ferns indicates a correlation between apomixis and the species richness of lineages; however BiSSE tests did not recover evidence for a correlation of apomixis and diversification rates. Instead, evidence was recovered supporting an association between the establishment of apomixis and reticulate evolution, especially in the establishment of triploid hybrids. Diversification time estimates supported the hypothesis of short living apomictic lineages and indicated a link between the establishment of apomixis and the strengthening of the monsoons caused by the lifting of the Qinghai-Tibetan plateau. In general our results supported the hypothesis for the rare establishment of apomictic lineages, high extinction risks, and low speciation rates.

Adaptive molecular evolution of a defence gene in sexual but not functionally asexual evening primroses

Theory predicts that sexual reproduction provides evolutionary advantages over asexual reproduction by reducing mutational load and increasing adaptive potential. Here, we test the latter prediction in the context of plant defences against pathogens because pathogens frequently reduce plant fitness and drive the evolution of plant defences. Specifically, we ask whether sexual evening primrose plant lineages (Onagraceae) have faster rates of adaptive molecular evolution and altered gene expression of a class I chitinase, a gene implicated in defence against pathogens, than functionally asexual evening primrose lineages. We found that the ratio of amino acid to silent substitutions (K(a) /K(s) = 0.19 vs. 0.11 for sexual and asexual lineages, respectively), the number of sites identified to be under positive selection (four vs. zero for sexual and asexual lineages, respectively) and the expression of chitinase were all higher in sexual than in asexual lineages. Our results are congruent with the conclusion that a loss of sexual recombination and segregation in the Onagraceae negatively affects adaptive structural and potentially regulatory evolution of a plant defence protein.

High genetic diversity in a remote island population system: sans sex

It has been proposed that long-distance dispersal of mosses to the Hawaiian Islands rarely occurs and that the Hawaiian population of the allopolyploid peat moss Sphagnum palustre probably resulted from a single dispersal event. Here, we used microsatellites to investigate whether the Hawaiian population of the dioicous S. palustre had a single founder and to compare its genetic diversity to that found in populations of S. palustre in other regions. The genetic diversity of the Hawaiian population is comparable to that of larger population systems. Several lines of evidence, including a lack of sporophytes and an apparently restricted natural distribution, suggest that sexual reproduction is absent in the Hawaiian plants. In addition, all samples of Hawaiian S. palustre share a genetic trait rare in other populations. Time to most recent ancestor (TMRCA) analysis indicates that the Hawaiian population was probably founded 49-51 kyr ago. It appears that all Hawaiian plants of S. palustre descend from a single founder via vegetative propagation. The long-term viability of this clonal population coupled with the development of significant genetic diversity suggests that vegetative propagation in a moss does not necessarily preclude evolutionary success in the long term.