Pleistocene diversification of unifoliolate-leaved Lupinus (Leguminosae: Papilionoideae) in Florida.

The importance and prevalence of recent ice-age and post-glacial speciation and species diversification during the Pleistocene across many organismal groups and physiographic settings are well established. However, the extent to which Pleistocene diversification can be attributed to climatic oscillations and their effects on distribution ranges and population structure remains debatable. In this study, we use morphologic, geographic and genetic (RADseq) data to document Pleistocene speciation and intra-specific diversification of the unifoliolate-leaved clade of Florida Lupinus, a small group of species largely restricted to inland and coastal sand ridges across the Florida peninsula and panhandle. Phylogenetic and demographic analyses alongside morphological and geographic evidence suggest that recent speciation and intra-specific divergence within this clade were driven by a combination of non-adaptive allopatric divergence caused by edaphic niche conservatism and opportunities presented by the emergence of new post-glacial sand ridge habitats. These results highlight the central importance of even modest geographic isolation and short periods of allopatric divergence following range expansion in the emergence of new taxa and add to the growing evidence that Pleistocene climatic oscillations may contribute to rapid diversification in a myriad of physiographic settings. Furthermore, our results shed new light on long-standing taxonomic debate surrounding the number of species in the Florida unifoliate Lupinus clade providing support for recognition of five species and a set of intra-specific variants. The important conservation implications for the narrowly restricted, highly endangered species Lupinus aridorum, which we show to be genetically distinct from its sister species Lupinus westianus, are discussed.

Rapid evolution of hybrid breakdown following recent divergence with gene flow in Senecio species on Mount Etna, Sicily.

How do nascent species evolve reproductive isolation during speciation with on-going gene flow? How do hybrid lineages become stabilised hybrid species? While commonly used genomic approaches provide an indirect way to identify species incompatibility factors, synthetic hybrids generated from interspecific crosses allow direct pinpointing of phenotypic traits involved in incompatibilities and the traits that are potentially adaptive in hybrid species. Here we report the analysis of phenotypic variation and hybrid breakdown in crosses between closely-related Senecio aethnensis and S. chrysanthemifolius, and their homoploid hybrid species, S. squalidus. The two former species represent a likely case of recent (<200 ky) speciation with gene flow driven by adaptation to contrasting conditions of high- and low-elevations on Mount Etna, Sicily. As these species form viable and fertile hybrids, it remains unclear whether they have started to evolve reproductive incompatibility. Our analysis represents the first study of phenotypic variation and hybrid breakdown involving multiple Senecio hybrid families. It revealed wide range of variation in multiple traits, including the traits previously unrecorded in synthetic hybrids. Leaf shape, highly distinct between S. aethnensis and S. chrysanthemifolius, was extremely variable in F2 hybrids, but more consistent in S. squalidus. Our study demonstrates that interspecific incompatibilities can evolve rapidly despite on-going gene flow between the species. Further work is necessary to understand the genetic bases of these incompatibilities and their role in speciation with gene flow.

Low Spontaneous Mutation Rate and Pleistocene Radiation of Pea Aphids.

Accurate estimates of divergence times are essential to understand the evolutionary history of species. It allows linking evolutionary histories of the diverging lineages with past geological, climatic, and other changes in environment and shed light on the processes involved in speciation. The pea aphid radiation includes multiple host races adapted to different legume host plants. It is thought that diversification in this system occurred very recently, over the past 8,000-16,000 years. This young age estimate was used to link diversification in pea aphids to the onset of human agriculture, and led to the establishment of the pea aphid radiation as a model system in the study of speciation with gene flow. Here, we re-examine the age of the pea aphid radiation, by combining a mutation accumulation experiment with a genome-wide estimate of divergence between distantly related pea aphid host races. We estimate the spontaneous mutation rate for pea aphids as 2.7×10-10 per haploid genome per parthenogenic generation. Using this estimate of mutation rate and the genome-wide genetic differentiation observed between pea aphid host races, we show that the pea aphid radiation is much more ancient than assumed previously, predating Neolithic agriculture by several hundreds of thousands of years. Our results rule out human agriculture as the driver of diversification of the pea aphid radiation, and call for re-assessment of the role of allopatric isolation during Pleistocene climatic oscillations in divergence of the pea aphid complex.

Senecio as a model system for integrating studies of genotype, phenotype and fitness.

Two major developments have made it possible to use examples of ecological radiations as model systems to understand evolution and ecology. First, the integration of quantitative genetics with ecological experiments allows detailed connections to be made between genotype, phenotype, and fitness in the field. Second, dramatic advances in molecular genetics have created new possibilities for integrating field and laboratory experiments with detailed genetic sequencing. Combining these approaches allows evolutionary biologists to better study the interplay between genotype, phenotype, and fitness to explore a wide range of evolutionary processes. Here, we present the genus Senecio (Asteraceae) as an excellent system to integrate these developments, and to address fundamental questions in ecology and evolution. Senecio is one of the largest and most phenotypically diverse genera of flowering plants, containing species ranging from woody perennials to herbaceous annuals. These Senecio species exhibit many growth habits, life histories, and morphologies, and they occupy a multitude of environments. Common within the genus are species that have hybridized naturally, undergone polyploidization, and colonized diverse environments, often through rapid phenotypic divergence and adaptive radiation. These diverse experimental attributes make Senecio an attractive model system in which to address a broad range of questions in evolution and ecology.

Rapid homoploid hybrid speciation in British gardens: The origin of Oxford ragwort (Senecio squalidus).

Hybridisation can lead to homoploid hybrid speciation, i.e., the origin of new species without change in chromosome number between parents and offspring. Central to homoploid hybrid speciation is the role of hybridisation in the establishment of reproductive isolation between the hybrid and the parental species in the early stages of speciation, when typically all species occur at least partly in sympatry. In this work we analyse genome-wide polymorphism data obtained by transcriptome sequencing of the British hybrid species Oxford ragwort (Senecio squalidus, Asteraceae), its two Italian parental species (S. aethnensis and S. chrysanthemifolius) and their naturally occurring hybrids on Mt Etna (Italy). We show that Oxford ragwort most likely originated from de novo hybridisation between its two Italian parental species whilst they were in cultivation in British gardens at the turn of the 18th century. Reproductive isolation between the new hybrid species and its parental species probably resulted from inheritance of genetic incompatibilities between the two parental species and subsequent ecological segregation - both of which have been shown in previous studies. Our results imply that S. squalidus meets the most stringent criteria set forth to identify homoploid hybrid speciation, and call attention to the creative role of hybridisation in responding to novel environmental conditions.

Strong divergent selection at multiple loci in two closely related species of ragworts adapted to high and low elevations on Mount Etna.

Recently diverged species present particularly informative systems for studying speciation and maintenance of genetic divergence in the face of gene flow. We investigated speciation in two closely related Senecio species, S. aethnensis and S. chrysanthemifolius, which grow at high and low elevations, respectively, on Mount Etna, Sicily and form a hybrid zone at intermediate elevations. We used a newly generated genome-wide single nucleotide polymorphism (SNP) dataset from 192 individuals collected over 18 localities along an elevational gradient to reconstruct the likely history of speciation, identify highly differentiated SNPs, and estimate the strength of divergent selection. We found that speciation in this system involved heterogeneous and bidirectional gene flow along the genome, and species experienced marked population size changes in the past. Furthermore, we identified highly-differentiated SNPs between the species, some of which are located in genes potentially involved in ecological differences between species (such as photosynthesis and UV response). We analysed the shape of these SNPs' allele frequency clines along the elevational gradient. These clines show significantly variable coincidence and concordance, indicative of the presence of multifarious selective forces. Selection against hybrids is estimated to be very strong (0.16-0.78) and one of the highest reported in literature. The combination of strong cumulative selection across the genome and previously identified intrinsic incompatibilities probably work together to maintain the genetic and phenotypic differentiation between these species - pointing to the importance of considering both intrinsic and extrinsic factors when studying divergence and speciation.

Resolving the backbone of the Brassicaceae phylogeny for investigating trait diversity.

The Brassicaceae family comprises c. 4000 species including economically important crops and the model plant Arabidopsis thaliana. Despite their importance, the relationships among major lineages in the family remain unresolved, hampering comparative research. Here, we inferred a Brassicaceae phylogeny using newly generated targeted enrichment sequence data of 1827 exons (> 940 000 bases) representing 63 species, as well as sequenced genome data of 16 species, together representing 50 of the 52 currently recognized Brassicaceae tribes. A third of the samples were derived from herbarium material, facilitating broad taxonomic coverage of the family. Six major clades formed successive sister groups to the rest of Brassicaceae. We also recovered strong support for novel relationships among tribes, and resolved the position of 16 taxa previously not assigned to a tribe. The broad utility of these phylogenetic results is illustrated through a comparative investigation of genome-wide expression signatures that distinguish simple from complex leaves in Brassicaceae. Our study provides an easily extendable dataset for further advances in Brassicaceae systematics and a timely higher-level phylogenetic framework for a wide range of comparative studies of multiple traits in an intensively investigated group of plants.

Pleistocene glacial cycles drive isolation, gene flow and speciation in the high-elevation Andes.

Mountain ranges are amongst the most species-rich habitats, with many large and rapid evolutionary radiations. The tempo and mode of diversification in these systems are key unanswered questions in evolutionary biology. Here we study the Andean Lupinus radiation to understand the processes driving very rapid diversification in montane systems. We use genomic and transcriptomic data of multiple species and populations, and apply phylogenomic and demographic analyses to test whether diversification proceeded without interspecific gene flow - as expected if Andean orogeny and geographic isolation were the main drivers of diversification - or if diversification was accompanied by gene flow, in which case other processes were probably involved. We uncover several episodes of gene flow between species, including very recent events likely to have been prompted by changes in habitat connectivity during Pleistocene glacial cycles. Furthermore, we find that gene flow between species was heterogeneously distributed across the genome. We argue that exceptionally fast diversification of Andean Lupinus was partly a result of Late Pleistocene glacial cycles, with associated cycles of expansion and contraction driving geographic isolation or secondary contact of species. Furthermore, heterogeneous gene flow across the genome suggests a role for selection and ecological speciation in rapid diversification in this system.

The western Mediterranean region provided the founder population of domesticated narrow-leafed lupin.

KEY MESSAGE: This study revealed that the western Mediterranean provided the founder population for domesticated narrow-leafed lupin and that genetic diversity decreased significantly during narrow-leafed lupin domestication. The evolutionary history of plants during domestication profoundly shaped the genome structure and genetic diversity of today's crops. Advances in next-generation sequencing technologies allow unprecedented opportunities to understand genome evolution in minor crops, which constitute the majority of plant domestications. A diverse set of 231 wild and domesticated narrow-leafed lupin (Lupinus angustifolius L.) accessions were subjected to genotyping-by-sequencing using diversity arrays technology. Phylogenetic, genome-wide divergence and linkage disequilibrium analyses were applied to identify the founder population of domesticated narrow-leafed lupin and the genome-wide effect of domestication on its genome. We found wild western Mediterranean population as the founder of domesticated narrow-leafed lupin. Domestication was associated with an almost threefold reduction in genome diversity in domesticated accessions compared to their wild relatives. Selective sweep analysis identified no significant footprints of selection around domestication loci. A genome-wide association study identified single nucleotide polymorphism markers associated with pod dehiscence. This new understanding of the genomic consequences of narrow-leafed lupin domestication along with molecular marker tools developed here will assist plant breeders more effectively access wild genetic diversity for crop improvement.

Exploring the genetic and adaptive diversity of a pan-Mediterranean crop wild relative: narrow-leafed lupin.

KEY MESSAGE: This first pan-Mediterranean analysis of genetic diversity in wild narrow-leafed lupin revealed strong East-West genetic differentiation of populations, an historic eastward migration, and signatures of genetic adaptation to climatic variables. Most grain crops suffer from a narrow genetic base, which limits their potential for adapting to new challenges such as increased stresses associated with climate change. Plant breeders are returning to the wild ancestors of crops and their close relatives to broaden the genetic base of their crops. Understanding the genetic adaptation of these wild relatives will help plant breeders most effectively use available wild diversity. Here, we took narrow-leafed lupin (Lupinus angustifolius L.) as a model to understand adaptation in a wild crop ancestor. A set of 142 wild accessions of narrow-leafed lupin from across the Mediterranean basin were subjected to genotyping-by-sequencing using Diversity Arrays Technology. Phylogenetic, linkage disequilibrium and demographic analyses were employed to explore the history of narrow-leafed lupin within the Mediterranean region. We found strong genetic differentiation between accessions from the western and eastern Mediterranean, evidence of an historic West to East migration, and that eastern Mediterranean narrow-leafed lupin experienced a severe and recent genetic bottleneck. We showed that these two populations differ for flowering time as a result of local adaptation, with the West flowering late while the East flowers early. A genome-wide association study identified single nucleotide polymorphism markers associated with climatic adaptation. Resolving the origin of wild narrow-leafed lupin and how its migration has induced adaptation to specific regions of the Mediterranean serves as a useful resource not only for developing narrow-leafed lupin cultivars with greater resilience to a changing climate, but also as a model which can be applied to other legumes.

Intrinsic pre-zygotic reproductive isolation of distantly related pea aphid host races.

Human activities may weaken or destroy reproductive isolation between young taxa, leading to their fusion with consequences for population and community ecology. Pea aphid host races are adapted to different legume taxa, providing a degree of pre-mating isolation mediated by habitat choice. Yet, all races can feed and reproduce on the broad bean (Vicia faba), a major crop which represents a 'universal host plant', which can promote hybridization between races. Here, we ask if pea aphid host races have reproductive barriers which prevent or reduce gene flow when they co-occur on the universal host plant. We observed mating behaviour, female survival, number of eggs and egg fertilization rates for three types of crosses: among individuals of the same host race, between closely related host races and between distantly related host races. We did not find significant differences in mating behaviour and female survival among the three types of crosses. However, we observed a drastic reduction in the number of eggs laid, and in the number of fertilized eggs, in distant crosses. We conclude that widespread broad bean cultivation in agriculture may predispose closely related-but not distantly related-host races to hybridize, disrupting reproductive isolation between incipient species.

Lost crops of the Incas: Origins of domestication of the Andean pulse crop tarwi, Lupinus mutabilis.

PREMISE OF THE STUDY: The Andean highlands are a hotspot of domestication, yet our understanding of the origins of early Andean agriculture remains fragmentary. Key questions of where, when, how many times, and from what progenitors many Andean crops were domesticated remain unanswered. The Andean lupine crop tarwi (Lupinus mutabilis) is a regionally important pulse crop with exceptionally high seed protein and oil content and is the focus of modern breeding efforts, but its origins remain obscure. METHODS: A large genome-wide DNA polymorphism data set was generated using nextRADseq to infer relationships among more than 200 accessions of Andean Lupinus species, including 24 accessions of L. mutabilis and close relatives. Phylogenetic and demographic analyses were used to identify the likely progenitor of tarwi and elucidate the area and timing of domestication in combination with archaeological evidence. KEY RESULTS: We infer that tarwi was domesticated once in northern Peru, most likely in the Cajamarca region within, or adjacent to the extant distribution of L. piurensis, which is the most likely wild progenitor. Demographic analyses suggest that tarwi split from L. piurensis around 2600 BP and suffered a classical domestication bottleneck. The earliest unequivocal archaeological evidence of domesticated tarwi seeds is from the Mantaro Valley, central Peru ca. 1800 BP. CONCLUSIONS: A single origin of tarwi from L. piurensis in northern Peru provides a robust working hypothesis for the domestication of this regionally important crop and is one of the first clear-cut examples of a crop originating in the highlands of northern Peru.

Genomic changes and stabilization following homoploid hybrid speciation of the Oxford ragwort Senecio squalidus.

Oxford ragwort (Senecio squalidus) is one of only two homoploid hybrid species known to have originated very recently, so it is a unique model for determining genomic changes and stabilization following homoploid hybrid speciation. Here, we provide a chromosome-level genome assembly of S. squalidus with 95% of the assembly contained in the 10 longest scaffolds, corresponding to its haploid chromosome number. We annotated 30,249 protein-coding genes and estimated that ∼62% of the genome consists of repetitive elements. We then characterized genome-wide patterns of linkage disequilibrium, polymorphism, and divergence in S. squalidus and its two parental species, finding that (1) linkage disequilibrium is highly heterogeneous, with a region on chromosome 4 showing increased values across all three species but especially in S. squalidus; (2) regions harboring genetic incompatibilities between the two parental species tend to be large, show reduced recombination, and have lower polymorphism in S. squalidus; (3) the two parental species have an unequal contribution (70:30) to the genome of S. squalidus, with long blocks of parent-specific ancestry supporting a very rapid stabilization of the hybrid lineage after hybrid formation; and (4) genomic regions with major parent ancestry exhibit an overrepresentation of loci with evidence for divergent selection occurring between the two parental species on Mount Etna. Our results show that both genetic incompatibilities and natural selection play a role in determining genome-wide reorganization following hybrid speciation and that patterns associated with homoploid hybrid speciation-typically seen in much older systems-can evolve very quickly following hybridization.

SNP calling by sequencing pooled samples.

BACKGROUND: Performing high throughput sequencing on samples pooled from different individuals is a strategy to characterize genetic variability at a small fraction of the cost required for individual sequencing. In certain circumstances some variability estimators have even lower variance than those obtained with individual sequencing. SNP calling and estimating the frequency of the minor allele from pooled samples, though, is a subtle exercise for at least three reasons. First, sequencing errors may have a much larger relevance than in individual SNP calling: while their impact in individual sequencing can be reduced by setting a restriction on a minimum number of reads per allele, this would have a strong and undesired effect in pools because it is unlikely that alleles at low frequency in the pool will be read many times. Second, the prior allele frequency for heterozygous sites in individuals is usually 0.5 (assuming one is not analyzing sequences coming from, e.g. cancer tissues), but this is not true in pools: in fact, under the standard neutral model, singletons (i.e. alleles of minimum frequency) are the most common class of variants because P(f) ∝ 1/f and they occur more often as the sample size increases. Third, an allele appearing only once in the reads from a pool does not necessarily correspond to a singleton in the set of individuals making up the pool, and vice versa, there can be more than one read - or, more likely, none - from a true singleton. RESULTS: To improve upon existing theory and software packages, we have developed a Bayesian approach for minor allele frequency (MAF) computation and SNP calling in pools (and implemented it in a program called snape): the approach takes into account sequencing errors and allows users to choose different priors. We also set up a pipeline which can simulate the coalescence process giving rise to the SNPs, the pooling procedure and the sequencing. We used it to compare the performance of snape to that of other packages. CONCLUSIONS: We present a software which helps in calling SNPs in pooled samples: it has good power while retaining a low false discovery rate (FDR). The method also provides the posterior probability that a SNP is segregating and the full posterior distribution of f for every SNP. In order to test the behaviour of our software, we generated (through simulated coalescence) artificial genomes and computed the effect of a pooled sequencing protocol, followed by SNP calling. In this setting, snape has better power and False Discovery Rate (FDR) than the comparable packages samtools, PoPoolation, Varscan : for N = 50 chromosomes, snape has power ≈ 35%and FDR ≈ 2.5%. snape is available at http://code.google.com/p/snape-pooled/ (source code and precompiled binaries).

Repeated unidirectional introgression of nuclear and mitochondrial DNA between four congeneric Tanganyikan cichlids.

With an increasing number of reported cases of hybridization and introgression, interspecific gene flow between animals has recently become a widely accepted and broadly studied phenomenon. In this study, we examine patterns of hybridization and introgression in Ophthalmotilapia spp., a genus of cichlid fish from Lake Tanganyika, using mitochondrial and nuclear DNA from all four species in the genus and including specimens from over 800 km of shoreline. These four species have very different, partially overlapping distribution ranges, thus allowing us to study in detail patterns of gene flow between sympatric and allopatric populations of the different species. We show that a significant proportion of individuals of the lake-wide distributed O. nasuta carry mitochondrial and/or nuclear DNA typical of other Ophthalmotilapia species. Strikingly, all such individuals were found in populations living in sympatry with each of the other Ophthalmotilapia species, strongly suggesting that this pattern originated by repeated and independent episodes of genetic exchange in different parts of the lake, with unidirectional introgression occurring into O. nasuta. Our analysis rejects the hypotheses that unidirectional introgression is caused by natural selection favoring heterospecific DNA, by skewed abundances of Ophthalmotilapia species or by hybridization events occurring during a putative spatial expansion in O. nasuta. Instead, cytonuclear incompatibilities or asymmetric behavioral reproductive isolation seem to have driven repeated, unidirectional introgression of nuclear and mitochondrial DNA into O. nasuta in different parts of the lake.

Adaptive divergence generates distinct plastic responses in two closely related Senecio species.

The evolution of plastic responses to external cues allows species to maintain fitness in response to the environmental variations they regularly experience. However, it remains unclear how plasticity evolves during adaptation. To test whether distinct patterns of plasticity are associated with adaptive divergence, we quantified plasticity for two closely related but ecologically divergent Sicilian daisy species (Senecio, Asteraceae). We sampled 40 representative genotypes of each species from their native range on Mt. Etna and then reciprocally transplanted multiple clones of each genotype into four field sites along an elevational gradient that included the native elevational range of each species, and two intermediate elevations. At each elevation, we quantified survival and measured leaf traits that included investment (specific leaf area), morphology, chlorophyll fluorescence, pigment content, and gene expression. Traits and differentially expressed genes that changed with elevation in one species often showed little changes in the other species, or changed in the opposite direction. As evidence of adaptive divergence, both species performed better at their native site and better than the species from the other habitat. Adaptive divergence is, therefore, associated with the evolution of distinct plastic responses to environmental variation, despite these two species sharing a recent common ancestor.

Comparative performances of DNA barcoding across insect orders.

BACKGROUND: Previous studies on insect DNA barcoding provide contradictory results and suggest not consistent performances across orders. This work aims at providing a general evaluation of insect DNA barcoding and "mini-barcoding" by performing simulations on a large database of 15,948 DNA barcodes. We compared the proportions of correctly identified queries across a) six insect orders (Coleoptera, Diptera, Hemiptera, Hymenoptera, Lepidoptera and Orthoptera), b) four identification criteria (Best Match: BM; Best Close Match: BCM; All Species Barcodes: ASB; tree-based identification: NJT), and c) reference databases with different taxon coverage (100, 500, 1,000, 1,500 and 1,995 insect species). RESULTS: Analysis of variance revealed highly significant differences among ID criteria and insect orders. A posteriori comparisons of means showed that NJT had always a significantly lower identification success (NJT = 0.656, S.D. = 0.118) compared to both BM and BCM (BM = 0.948, S.D. = 0.026; BCM = 0.946, S.D. = 0.031). NJT showed significant variations among orders, with the highest proportion of correctly identified queries in Hymenoptera and Orthoptera and the lowest in Diptera. Conversely, the proportions of correct matches of BM and BCM were consistent across orders but a progressive increase in false identification was observed when larger reference databases were used. CONCLUSIONS: Regardless the relatively low proportion of Type I errors (misidentification of queries which are represented in the reference database) of BM and BCM, the lack of reference DNA barcodes for 98% of the known insect species implies that insect DNA barcoding is heavily biased by Type II errors (misidentification of queries without conspecifics in the database). The detrimental effects of Type II errors could be circumvented if insect DNA barcoding is used to verify the lack of correspondence between a query and a list of properly referenced target species (e.g. insect pests). This "negative identification" would only be subjected to Type I errors and could be profitably adopted in insect quarantine procedures.

When environmental changes do not cause geographic separation of fauna: differential responses of Baikalian invertebrates.

BACKGROUND: While the impact of climate fluctuations on the demographic histories of species caused by changes in habitat availability is well studied, populations of species from systems without geographic isolation have received comparatively little attention. Using CO1 mitochondrial sequences, we analysed phylogeographic patterns and demographic histories of populations of five species (four gastropod and one amphipod species) co-occurring in the southwestern shore of Lake Baikal, an area where environmental oscillations have not resulted in geographical isolation of habitats. RESULTS: Species with stronger habitat preferences (gastropods B. turriformis, B. carinata and B. carinatocostata) exhibit rather stable population sizes through their evolutionary history, and their phylogeographic pattern indicates moderate habitat fragmentation. Conversely, species without strong habitat preference (gastropod M. herderiana and amphipod G. fasciatus) exhibit haplotype networks with a very abundant and widespread central haplotype and a big number of singleton haplotypes, while their reconstructed demographic histories show a population expansion starting about 25-50 thousand years ago, a period marked by climate warming and increase in diatom abundance as inferred from bottom-lake sedimentary cores. CONCLUSIONS: In agreement with previous studies, we found that species reacted differently to the same environmental changes. Our results highlight the important role of dispersal ability and degree of ecological specialization in defining a species' response to environmental changes.

Biogeographic origin and radiation of Cuban Eleutherodactylus frogs of the auriculatus species group, inferred from mitochondrial and nuclear gene sequences.

We studied phylogenetic relationships of the Eleutherodactylus auriculatus species group to infer colonization and diversification patterns in this endemic radiation of terrestrial frogs of the genus Eleutherodactylus in the largest of the Greater Antilles, Cuba. An initial screening of genetic diversity based on partial sequences of the 16S rRNA gene in almost 100 individuals of all species of the group and covering the complete known geographic range of their occurrence found most species endemic to small ranges in the eastern Cuban mountains while a single species was widespread over most of Cuba. Our molecular phylogeny, based on 3731 bp of four mitochondrial and one nuclear gene, suggests that most cladogenetic events within the group occurred among clades restricted to the eastern mountains, which acted as refugia and facilitated the diversification in this group. Our results reveal two separate colonization events of Central and Western Cuba and allow inferring the timing of the subsequent diversification events that occurred between 11 and 2 Mya. Because populations previously assigned to E. auriculatus represent four genetically strongly divergent lineages that also differ in their advertisement calls, we propose that E. auriculatus as currently recognized comprises four species. The difficulties in assigning the name auriculatus to any of these four species, and the fact that E. principalis is nested within one of them, stress the need for a thorough taxonomic revision of this group.

Adaptive Evolution Is Common in Rapid Evolutionary Radiations.

One of the most long-standing and important mysteries in evolutionary biology is why biological diversity is so unevenly distributed across space and taxonomic lineages. Nowhere is this disparity more evident than in the multitude of rapid evolutionary radiations found on oceanic islands and mountain ranges across the globe [1-5]. The evolutionary processes driving these rapid diversification events remain unclear [6-8]. Recent genome-wide studies suggest that natural selection may be frequent during rapid evolutionary radiations, as inferred from work in cichlid fish [9], white-eye birds [10], new world lupins [11], and wild tomatoes [12]. However, whether frequent adaptive evolution is a general feature of rapid evolutionary radiations remains untested. Here we show that adaptive evolution is significantly more frequent in rapid evolutionary radiations compared to background levels in more slowly diversifying lineages. This result is consistent across a wide range of angiosperm lineages analyzed: 12 evolutionary radiations, which together comprise 1,377 described species, originating from some of the most biologically diverse systems on Earth. In addition, we find a significant negative correlation between population size and frequency of adaptive evolution in rapid evolutionary radiations. A possible explanation for this pattern is that more frequent adaptive evolution is at least partly driven by positive selection for advantageous mutations that compensate for the fixation of slightly deleterious mutations in smaller populations.