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.

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.

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.

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.

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.

Multiple macroevolutionary routes to becoming a biodiversity hotspot

Why is species diversity so unevenly distributed across different regions on Earth? Regional differences in biodiversity may stem from differences in rates of speciation and dispersal and colonization times, but these hypotheses have rarely been tested simultaneously at a global scale. Our study reveals the macroevolutionary routes that have generated hotspots of mammal and bird biodiversity by analyzing the tempo and mode of diversification and dispersal within major biogeographic realms. Hotspots in tropical realms had higher rates of speciation, whereas those in temperate realms received more immigrant species from their surrounding regions. We also found that hotspots had higher spatial complexity and energy availability, providing a link between the environment and macroevolutionary history. Our study highlights how assessing differences in macroevolutionary history can help to explain why biodiversity varies so much worldwide.

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.

Ecological speciation in sympatric palms: 1. Gene expression, selection and pleiotropy

Ecological speciation requires divergent selection, reproductive isolation and a genetic mechanism to link the two. We examined the role of gene expression and coding sequence evolution in this process using two species of Howea palms that have diverged sympatrically on Lord Howe Island, Australia. These palms are associated with distinct soil types and have displaced flowering times, representing an ideal candidate for ecological speciation. We generated large amounts of RNA‐Seq data from multiple individuals and tissue types collected on the island from each of the two species. We found that differentially expressed loci as well as those with divergent coding sequences between Howea species were associated with known ecological and phenotypic differences, including response to salinity, drought, pH and flowering time. From these loci, we identified potential ‘ecological speciation genes’ and further validate their effect on flowering time by knocking out orthologous loci in a model plant species. Finally, we put forward six plausible ecological speciation loci, providing support for the hypothesis that pleiotropy could help to overcome the antagonism between selection and recombination during speciation with gene flow.

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.

Electronic supplementary material

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

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.

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.

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.

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

Posthemipelvectomy hernia

We report the case of a white male who underwent a classic hemipelvectomy due to a femur fibrosarcoma with inguinal metastases, which 33 years later, developed into a posthemipelvectomy hernia in the amputation stump that impaired the use of his Canadian prosthesis. The hernia was repaired with a polypropylene mesh in a subaponeurotic position. A seroma was drained in the postoperative and it was only 2 months after the operation that he could use his prosthesis with any difficulty. A year after the operation, the hernia had not recurred. Only seven similar cases have been published, and there are only four cases with details of their correction, two with a mesh as was our case, and the rest with a primary suture of the aponeurotic borders. A brief review of the bibliography is given on this subject.

[Other indications for growth hormone treatment: chronic kidney failure and Turner's syndrome]

Recently, it has been demonstrated that treatment with growth hormone (GH) can accelerate height velocity in children with chronic renal insufficiency (CRI), after kidney transplantation and in Turner syndrome. The pathogenesis of growth retardation in CRI is complex. Possibly the most important factor involved is the presence of peripheral resistance to insulin-like growth factors. The administration of GH in high doses may restore catch-up growth, both in children with CRI and in kidney post-transplant patients. After 12 months of treatment height velocity increased from 4.3 to 6.6 cm/year in 4 children with CRI, and from 2.5 to 7.4 cm/year in 3 patients with kidney transplant. GH therapy alone, but even more in combination with oxandrolone increases the growth rate in Turner syndrome and may increase the final height, but the long-term results of this treatment are awaited. No undesirable collateral effects have been reported.

[Diabetic ketoacidosis: results with two different patterns of insulin administration (author's transl)]

The response to insulin treatment in 15 children with diabetic ketoacidosis is studied. Insulin continous infusion was administered to nine patients. The other six patients received insulin subcutaneously except one half of the first dose injected by vein as a bolus. Both patterns of insulin administration proved to be equally effective. Anyway the insulin continuous infusion appears as a simple and easier method for a good general control and to prevent hypoglycemic episodes. To prevent a hyperglycemic rebound it is suggested that a dose of subcutaneous insulin must be injected immediately after the insulin infusion is discontinued. A discordance between hyperglycemia correction and the degree of acidosis has been noticed in some patients; this can be corrected decreasing insulin infusion rate.

[Stimulation of growth hormone (HGH) secretion by physical exercise (author's transl)]

Authors study forty eight children with severe lineal growth retardation. The more frequent causes of a stunted growth were ruled out previously (systemic and chronic diseases, metabolic disorders, genetic alterations, etc.). As a screening test, physical exercise was used to detect growth hormone (HGH) deficiency. A normal response was considered when HGH levels rose more than 7,0 ng./ml. after 20 minutes of continuous exercise. In those cases in which this level was not attained. insulin-hypoglycemia alone and followed by insulin plus an arginine infusion were used as stimuli, to confirm the lack of HGH response. A correct response to the exercise was found in 81% of the cases. It is concluded that the physical exercise as HGH stimulus is a good screening test. It presents a high degree of confidence with some other qualities, mainly: little disturbances to the patients, a rapid performance, lack of unpleasant colateral effects, and the possibility that it can be carried out by sanitary non-medical personnel.

[Infantile cushing's disease with paradoxical response to dexamethasona due to a probable periodic hormogenesis (author's transl)]

A case of Cushing's disease in a ten year old girl with bilateral diffuse hyperplasia of the adrenal cortex and postoperative enlargment of the sella turcica is presented. The administration of dexamethasone elicited a paradoxical response with a clear elevation of the already high excretion of 17-hydroxycorticoids. The possible mechanisms for this previously described, but infrequent, response to dexamethasone are discussed. It is concluded that: 1) Interpretation of the classical dexamethasone suppression test can occasionally be misleading; and 2) Periodic hormonogenesis may account for this type of paradoxical response.