Using target sequence capture to improve the phylogenetic resolution of a rapid radiation in New Zealand Veronica.

PREMISE: Recent, rapid radiations present a challenge for phylogenetic reconstruction. Fast successive speciation events typically lead to low sequence divergence and poorly resolved relationships with standard phylogenetic markers. Target sequence capture of many independent nuclear loci has the potential to improve phylogenetic resolution for rapid radiations. METHODS: Here we applied target sequence capture with 353 protein-coding genes (Angiosperms353 bait kit) to Veronica sect. Hebe (common name hebe) to determine its utility for improving the phylogenetic resolution of rapid radiations. Veronica section Hebe originated 5-10 million years ago in New Zealand, forming a monophyletic radiation of ca 130 extant species. RESULTS: We obtained approximately 150 kbp of 353 protein-coding exons and an additional 200 kbp of flanking noncoding sequences for each of 77 hebe and two outgroup species. When comparing coding, noncoding, and combined data sets, we found that the latter provided the best overall phylogenetic resolution. While some deep nodes in the radiation remained unresolved, our phylogeny provided broad and often improved support for subclades identified by both morphology and standard markers in previous studies. Gene-tree discordance was nonetheless widespread, indicating that additional methods are needed to disentangle fully the history of the radiation. CONCLUSIONS: Phylogenomic target capture data sets both increase phylogenetic signal and deliver new insights into the complex evolutionary history of rapid radiations as compared with traditional markers. Improving methods to resolve remaining discordance among loci from target sequence capture is now important to facilitate the further study of rapid radiations.

Angiosperm speciation cools down in the tropics.

Recent evidence has questioned whether the Latitudinal Diversity Gradient (LDG), whereby species richness increases towards the Equator, results in higher rates of speciation in the tropics. Allowing for time heterogeneity in speciation rate estimates for over 60,000 angiosperm species, we found that the LDG does not arise from variation in speciation rates because lineages do not speciate faster in the tropics. These results were consistently retrieved using two other methods to test the association between occupancy of tropical habitats and speciation rates. Our speciation rate estimates were robust to the effects of both undescribed species and missing taxa. Overall, our results show that speciation rates follow an opposite pattern to global variation in species richness. Greater ecological opportunity in the temperate zones, stemming from less saturated communities, higher species turnover or greater environmental change, may ultimately explain these results.

Does Evolutionary History Correlate with Contemporary Extinction Risk by Influencing Range Size Dynamics?

Extinction threatens many species yet is predicted by few factors across the plant tree of life (ToL). Taxon age is one factor that may associate with extinction if occupancy of geographic and adaptive zones varies with time, but evidence for such an association has been equivocal. Age-dependent occupancy can also influence diversification rates and thus extinction risk where new taxa have small range and population sizes. To test how age, diversification, and range size were correlated with extinction, we analyzed 639 well-sampled genera representing 8,937 species from across the plant ToL. We found a greater proportion of species were threatened by contemporary extinction in younger and faster-diversifying genera. When we directly tested how range size mediated this pattern in two large, well-sampled groups, our results varied. In conifers, potential range size was smaller in older species and was correlated with higher extinction risk. Age on its own had no direct effect on extinction when accounting for its influence on range size. In palm species, age was neither directly nor indirectly correlated with extinction risk. Our results suggest that range size dynamics may explain differing patterns of extinction risk across the ToL, with consequences for biodiversity conservation.

Seed size and its rate of evolution correlate with species diversification across angiosperms.

Species diversity varies greatly across the different taxonomic groups that comprise the Tree of Life (ToL). This imbalance is particularly conspicuous within angiosperms, but is largely unexplained. Seed mass is one trait that may help clarify why some lineages diversify more than others because it confers adaptation to different environments, which can subsequently influence speciation and extinction. The rate at which seed mass changes across the angiosperm phylogeny may also be linked to diversification by increasing reproductive isolation and allowing access to novel ecological niches. However, the magnitude and direction of the association between seed mass and diversification has not been assessed across the angiosperm phylogeny. Here, we show that absolute seed size and the rate of change in seed size are both associated with variation in diversification rates. Based on the largest available angiosperm phylogenetic tree, we found that smaller-seeded plants had higher rates of diversification, possibly due to improved colonisation potential. The rate of phenotypic change in seed size was also strongly positively correlated with speciation rates, providing rare, large-scale evidence that rapid morphological change is associated with species divergence. Our study now reveals that variation in morphological traits and, importantly, the rate at which they evolve can contribute to explaining the extremely uneven distribution of diversity across the ToL.

Ten-year assessment of the 100 priority questions for global biodiversity conservation.

In 2008, a group of conservation scientists compiled a list of 100 priority questions for the conservation of the world's biodiversity. However, now almost a decade later, no one has yet published a study gauging how much progress has been made in addressing these 100 high-priority questions in the peer-reviewed literature. We took a first step toward reexamining the 100 questions to identify key knowledge gaps that remain. Through a combination of a questionnaire and a literature review, we evaluated each question on the basis of 2 criteria: relevance and effort. We defined highly relevant questions as those that - if answered - would have the greatest impact on global biodiversity conservation and quantified effort based on the number of review publications addressing a particular question, which we used as a proxy for research effort. Using this approach, we identified a set of questions that, despite being perceived as highly relevant, have been the focus of relatively few review publications over the past 10 years. These questions covered a broad range of topics but predominantly tackled 3 major themes: conservation and management of freshwater ecosystems, role of societal structures in shaping interactions between people and the environment, and impacts of conservation interventions. We believe these questions represent important knowledge gaps that have received insufficient attention and may need to be prioritized in future research.

Evidence of positive selection associated with placental loss in tiger sharks.

BACKGROUND: All vertebrates initially feed their offspring using yolk reserves. In some live-bearing species these yolk reserves may be supplemented with extra nutrition via a placenta. Sharks belonging to the Carcharhinidae family are all live-bearing, and with the exception of the tiger shark (Galeocerdo cuvier), develop placental connections after exhausting yolk reserves. Phylogenetic relationships suggest the lack of placenta in tiger sharks is due to secondary loss. This represents a dramatic shift in reproductive strategy, and is likely to have left a molecular footprint of positive selection within the genome. RESULTS: We sequenced the transcriptome of the tiger shark and eight other live-bearing shark species. From this data we constructed a time-calibrated phylogenetic tree estimating the tiger shark lineage diverged from the placental carcharhinids approximately 94 million years ago. Along the tiger shark lineage, we identified five genes exhibiting a signature of positive selection. Four of these genes have functions likely associated with brain development (YWHAE and ARL6IP5) and sexual reproduction (VAMP4 and TCTEX1D2). CONCLUSIONS: Our results indicate the loss of placenta in tiger sharks may be associated with subsequent adaptive changes in brain development and sperm production.

Multilocus species trees and species delimitation in a temporal context: application to the water shrews of the genus Neomys.

BACKGROUND: Multilocus data are becoming increasingly important in determining the phylogeny of closely related species and delimiting species. In species complexes where unequivocal fossil calibrations are not available, rigorous dating of the coalescence-based species trees requires accurate mutation rates of the loci under study but, generally, these rates are unknown. Here, we obtained lineage-specific mutation rates of these loci from a higher-level phylogeny with a reliable fossil record and investigated how different choices of mutation rates and species tree models affected the split time estimates. We implemented this strategy with a genus of water shrews, Neomys, whose taxonomy has been contentious over the last century. RESULTS: We sequenced 13 introns and cytochrome b from specimens of the three species currently recognized in this genus including two subspecies of N. anomalus that were originally described as species. A Bayesian multilocus species delimitation method and estimation of gene flow supported that these subspecies are distinct evolutionary lineages that should be treated as distinct species: N. anomalus (sensu stricto), limited to part of the Iberian Peninsula, and N. milleri, with a larger Eurasian range. We then estimated mutation rates from a Bayesian relaxed clock analysis of the mammalian orthologues with several fossil calibrations. Next, using the estimated Neomys-specific rates for each locus in an isolation-with-migration model, the split time for these sister taxa was dated at 0.40 Myr ago (with a 95 % confidence interval of 0.26 - 0.86 Myr), likely coinciding with one of the major glaciations of the Middle Pleistocene. We also showed that the extrapolation of non-specific rates or the use of simpler models would lead to very different split time estimates. CONCLUSIONS: We showed that the estimation of rigorous lineage-specific mutation rates for each locus allows the inference of robust split times in a species tree framework. These times, in turn, afford a better understanding of the timeframe required to achieve isolation and, eventually, speciation in sister lineages. The application of species delimitation methods and an accurate dating strategy to the genus Neomys helped to clarify its controversial taxonomy.

Phylogeography and postglacial expansion of the endangered semi-aquatic mammal Galemys pyrenaicus.

BACKGROUND: Species with strict ecological requirements may provide new insights into the forces that shaped the geographic variation of genetic diversity. The Pyrenean desman, Galemys pyrenaicus, is a small semi-aquatic mammal that inhabits clean streams of the northern half of the Iberian Peninsula and is endangered in most of its geographic range, but its genetic structure is currently unknown. While the stringent ecological demands derived from its aquatic habitat might have caused a partition of the genetic diversity among river basins, Pleistocene glaciations would have generated a genetic pattern related to glacial refugia. RESULTS: To study the relative importance of historical and ecological factors in the genetic structure of G. pyrenaicus, we used mitochondrial and intronic sequences of specimens covering most of the species range. We show, first, that the Pyrenean desman has very low levels of genetic diversity compared to other mammals. In addition, phylogenetic and dating analyses of the mitochondrial sequences reveal a strong phylogeographic structure of a Middle Pleistocene origin, suggesting that the main lineages arose during periods of glacial isolation. Furthermore, both the spatial distribution of nuclear and mitochondrial diversity and the results of species distribution modeling suggest the existence of a major glacial refugium in the northwestern part of the Iberian Peninsula. Finally, the main mitochondrial lineages show a striking parapatric distribution without any apparent exchange of mitochondrial haplotypes between the lineages that came into secondary contact (although with certain permeability to nuclear genes), indicating incomplete mixing after the post-glacial recolonization. On the other hand, when we analyzed the partition of the genetic diversity among river basins, the Pyrenean desman showed a lower than expected genetic differentiation among main rivers. CONCLUSIONS: The analysis of mitochondrial and intronic markers in G. pyrenaicus showed the predominant effects of Pleistocene glaciations on the genetic structure of this species, while the distribution of the genetic diversity was not greatly influenced by the main river systems. These results and, particularly, the discovery of a marked phylogeographic structure, may have important implications for the conservation of the Pyrenean desman.

Global topographic uplift has elevated speciation in mammals and birds over the last 3 million years.

Topographic change shapes the evolution of biodiversity by influencing both habitat connectivity and habitat diversity as well as abiotic factors like climate. However, its role in creating global biodiversity gradients remains poorly characterized because geology, climate and evolutionary data have rarely been integrated across concordant timescales. Here we show that topographic uplift over the last 3 million years explains more spatial variation in the speciation of all mammals and birds than do the direct effects of palaeoclimate change and both present-day elevation and present-day temperature. By contrast, the effects of topographic changes are much smaller than those of present-day temperatures in eroded areas. Together, our results stress that historical geological processes rather than traditionally studied macroecological gradients may ultimately generate much of the world's biodiversity. More broadly, as the Earth's surface continues to rise and fall, topography will remain an important driver of evolutionary change and novelty.

Polyploidy on Islands: Its Emergence and Importance for Diversification.

Whole genome duplication or polyploidy is widespread among floras globally, but traditionally has been thought to have played a minor role in the evolution of island biodiversity, based on the low proportion of polyploid taxa present. We investigate five island systems (Juan Fernández, Galápagos, Canary Islands, Hawaiian Islands, and New Zealand) to test whether polyploidy (i) enhances or hinders diversification on islands and (ii) is an intrinsic feature of a lineage or an attribute that emerges in island environments. These island systems are diverse in their origins, geographic and latitudinal distributions, levels of plant species endemism (37% in the Galapagos to 88% in the Hawaiian Islands), and ploidy levels, and taken together are representative of islands more generally. We compiled data for vascular plants and summarized information for each genus on each island system, including the total number of species (native and endemic), generic endemicity, chromosome numbers, genome size, and ploidy levels. Dated phylogenies were used to infer lineage age, number of colonization events, and change in ploidy level relative to the non-island sister lineage. Using phylogenetic path analysis, we then tested how the diversification of endemic lineages varied with the direct and indirect effects of polyploidy (presence of polyploidy, time on island, polyploidization near colonization, colonizer pool size) and other lineage traits not associated with polyploidy (time on island, colonizer pool size, repeat colonization). Diploid and tetraploid were the most common ploidy levels across all islands, with the highest ploidy levels (>8x) recorded for the Canary Islands (12x) and New Zealand (20x). Overall, we found that endemic diversification of our focal island floras was shaped by polyploidy in many cases and certainly others still to be detected considering the lack of data in many lineages. Polyploid speciation on the islands was enhanced by a larger source of potential congeneric colonists and a change in ploidy level compared to overseas sister taxa.

Novel intron markers to study the phylogeny of closely related mammalian species.

BACKGROUND: Multilocus phylogenies can be used to infer the species tree of a group of closely related species. In species trees, the nodes represent the actual separation between species, thus providing essential information about their evolutionary history. In addition, multilocus phylogenies can help in analyses of species delimitation, gene flow and genetic differentiation within species. However, few adequate markers are available for such studies. RESULTS: In order to develop nuclear markers that can be useful in multilocus studies of mammals, we analyzed the mammalian genomes of human, chimpanzee, macaque, dog and cow. Rodents were excluded due to their unusual genomic features. Introns were extracted from the mammalian genomes because of their greater genetic variability and ease of amplification from the flanking exons. To an initial set of more than 10,000 one-to-one orthologous introns we applied several filters to select introns that belong to single-copy genes, show neutral evolutionary rates and have an adequate length for their amplification. This analysis led to a final list of 224 intron markers randomly distributed along the genome. To experimentally test their validity, we amplified twelve of these introns in a panel of six mammalian species. The result was that seven of these introns gave rise to a PCR band of the expected size in all species. In addition, we sequenced these bands and analyzed the accumulation of substitutions in these introns in five pairs of closely related species. The results showed that the estimated genetic distances in the five species pairs was quite variable among introns and that this divergence cannot be directly predicted from the overall intron divergence in mammals. CONCLUSIONS: We have designed a new set of 224 nuclear introns with optimal features for the phylogeny of closely related mammalian species. A large proportion of the introns tested experimentally showed a perfect amplification and enough variability in most species, indicating that this marker set can be very helpful in multilocus phylogenetics of mammals. Due to the lower variability and stronger stochasticity of nuclear markers with respect to mitochondrial genes, studies should be designed to make use of several markers like the ones designed here.

Ecological speciation in sympatric palms: 4. Demographic analyses support speciation of Howea in the face of high gene flow.

The idea that populations must be geographically isolated (allopatric) to evolve into separate species has persisted for a long time. It is now clear that new species can also diverge despite ongoing genetic exchange, but few accepted cases of speciation in sympatry have held up when scrutinized using modern approaches. Here, we examined evidence for speciation of the Howea palms of Lord Howe Island, Australia, in light of new genomic data. We used coalescence-based demographic models combined with double digest restriction site associated DNA sequencing of multiple individuals and provide support for previous claims by Savolainen et al. that speciation in Howea did occur in the face of gene flow.

Ecological speciation in sympatric palms: 3. Genetic map reveals genomic islands underlying species divergence in Howea.

Although it is now widely accepted that speciation can occur in the face of continuous gene flow, with little or no spatial separation, the mechanisms and genomic architectures that permit such divergence are still debated. Here, we examined speciation in the face of gene flow in the Howea palms of Lord Howe Island, Australia. We built a genetic map using a novel method applicable to long-lived tree species, combining it with double digest restriction site-associated DNA sequencing of multiple individuals. Based upon various metrics, we detected 46 highly differentiated regions throughout the genome, four of which contained genes with functions that are particularly relevant to the speciation scenario for Howea, specifically salt and drought tolerance.

The K tree score: quantification of differences in the relative branch length and topology of phylogenetic trees.

SUMMARY: We introduce a new phylogenetic comparison method that measures overall differences in the relative branch length and topology of two phylogenetic trees. To do this, the algorithm first scales one of the trees to have a global divergence as similar as possible to the other tree. Then, the branch length distance, which takes differences in topology and branch lengths into account, is applied to the two trees. We thus obtain the minimum branch length distance or K tree score. Two trees with very different relative branch lengths get a high K score whereas two trees that follow a similar among-lineage rate variation get a low score, regardless of the overall rates in both trees. There are several applications of the K tree score, two of which are explained here in more detail. First, this score allows the evaluation of the performance of phylogenetic algorithms, not only with respect to their topological accuracy, but also with respect to the reproduction of a given branch length variation. In a second example, we show how the K score allows the selection of orthologous genes by choosing those that better follow the overall shape of a given reference tree. AVAILABILITY: http://molevol.ibmb.csic.es/Ktreedist.html

Viviparity stimulates diversification in an order of fish.

Species richness is distributed unevenly across the tree of life and this may be influenced by the evolution of novel phenotypes that promote diversification. Viviparity has originated ∼150 times in vertebrates and is considered to be an adaptation to highly variable environments. Likewise, possessing an annual life cycle is common in plants and insects, where it enables the colonization of seasonal environments, but rare in vertebrates. The extent to which these reproductive life-history traits have enhanced diversification and their relative importance in the process remains unknown. We show that convergent evolution of viviparity causes bursts of diversification in fish. We built a phylogenetic tree for Cyprinodontiformes, an order in which both annualism and viviparity have arisen, and reveal that while both traits have evolved multiple times, only viviparity played a major role in shaping the patterns of diversity. These results demonstrate that changes in reproductive life-history strategy can stimulate diversification.

A comparative analysis of island floras challenges taxonomy-based biogeographical models of speciation.

Speciation on islands, and particularly the divergence of species in situ, has long been debated. Here, we present one of the first, complete assessments of the geographic modes of speciation for the flora of a small oceanic island. Cocos Island (Costa Rica) is pristine; it is located 550 km off the Pacific coast of Central America. It harbors 189 native plant species, 33 of which are endemic. Using phylogenetic data from insular and mainland congeneric species, we show that all of the endemic species are derived from independent colonization events rather than in situ speciation. This is in sharp contrast to the results of a study carried out in a comparable system, Lord Howe Island (Australia), where as much as 8.2% of the plant species were the product of sympatric speciation. Differences in physiography and age between the islands may be responsible for the contrasting patterns of speciation observed. Importantly, comparing phylogenetic assessments of the modes of speciation with taxonomy-based measures shows that widely used island biogeography approaches overestimate rates of in situ speciation.

Development of rapidly evolving intron markers to estimate multilocus species trees of rodents.

One of the major challenges in the analysis of closely related species, speciation and phylogeography is the identification of variable sequence markers that allow the determination of genealogical relationships in multiple genomic regions using coalescent and species tree approaches. Rodent species represent nearly half of the mammalian diversity, but so far no systematic study has been carried out to detect suitable informative markers for this group. Here, we used a bioinformatic pipeline to extract intron sequences from rodent genomes available in databases and applied a series of filters that allowed the identification of 208 introns that adequately fulfilled several criteria for these studies. The main required characteristics of the introns were that they had the maximum possible mutation rates, that they were part of single-copy genes, that they had an appropriate sequence length for amplification, and that they were flanked by exons with suitable regions for primer design. In addition, in order to determine the validity of this approach, we chose ten of these introns for primer design and tested them in a panel of eleven rodent species belonging to different representative families. We show that all these introns can be amplified in the majority of species and that, overall, 79% of the amplifications worked with minimum optimization of the annealing temperature. In addition, we confirmed for a pair of sister species the relatively high level of sequence divergence of these introns. Therefore, we provide here a set of adequate intron markers that can be applied to different species of Rodentia for their use in studies that require significant sequence variability.