Dataset Information

Data concatenation, Bayesian concordance and coalescent-based analyses of the species tree for the rapid radiation of Triturus newts.

ABSTRACT: The phylogenetic relationships for rapid species radiations are difficult to disentangle. Here we study one such case, namely the genus Triturus, which is composed of the marbled and crested newts. We analyze data for 38 genetic markers, positioned in 3-prime untranslated regions of protein-coding genes, obtained with 454 sequencing. Our dataset includes twenty Triturus newts and represents all nine species. Bayesian analysis of population structure allocates all individuals to their respective species. The branching patterns obtained by data concatenation, Bayesian concordance analysis and coalescent-based estimations of the species tree differ from one another. The data concatenation based species tree shows high branch support but branching order is considerably affected by allele choice in the case of heterozygotes in the concatenation process. Bayesian concordance analysis expresses the conflict between individual gene trees for part of the Triturus species tree as low concordance factors. The coalescent-based species tree is relatively similar to a previously published species tree based upon morphology and full mtDNA and any conflicting internal branches are not highly supported. Our findings reflect high gene tree discordance due to incomplete lineage sorting (possibly aggravated by hybridization) in combination with low information content of the markers employed (as can be expected for relatively recent species radiations). This case study highlights the complexity of resolving rapid radiations and we acknowledge that to convincingly resolve the Triturus species tree even more genes will have to be consulted.

SUBMITTER: Wielstra B

PROVIDER: S-EPMC4206468 | biostudies-literature | 2014

REPOSITORIES: biostudies-literature

ACCESS DATA

Publications

Data concatenation, Bayesian concordance and coalescent-based analyses of the species tree for the rapid radiation of Triturus newts.

Wielstra Ben B Arntzen Jan W JW van der Gaag Kristiaan J KJ Pabijan Maciej M Babik Wieslaw W

PloS one 20141022 10

The phylogenetic relationships for rapid species radiations are difficult to disentangle. Here we study one such case, namely the genus Triturus, which is composed of the marbled and crested newts. We analyze data for 38 genetic markers, positioned in 3-prime untranslated regions of protein-coding genes, obtained with 454 sequencing. Our dataset includes twenty Triturus newts and represents all nine species. Bayesian analysis of population structure allocates all individuals to their respective ...[more]

PMID: 25337997

Similar Datasets

Project description:BackgroundThe rapid radiation of crested newts (Triturus cristatus superspecies) comprises four morphotypes: 1) the T. karelinii group, 2) T. carnifex - T. macedonicus, 3) T. cristatus and 4) T. dobrogicus. These vary in body build and the number of rib-bearing pre-sacral vertebrae (NRBV). The phylogenetic relationships of the morphotypes have not yet been settled, despite several previous attempts, employing a variety of molecular markers. We here resolve the crested newt phylogeny by using complete mitochondrial genome sequences.ResultsBayesian inference based on the mitogenomic data yields a fully bifurcating, significantly supported tree, though Maximum Likelihood inference yields low support values. The internal branches connecting the morphotypes are short relative to the terminal branches. Seen from the root of Triturus (NRBV = 13), a basal dichotomy separates the T. karelinii group (NRBV = 13) from the remaining crested newts. The next split divides the latter assortment into T. carnifex - T. macedonicus (NRBV = 14) versus T. cristatus (NRBV = 15) and T. dobrogicus (NRBV = 16 or 17).ConclusionsWe argue that the Bayesian full mitochondrial DNA phylogeny is superior to previous attempts aiming to recover the crested newt species tree. Furthermore, our new phylogeny involves a maximally parsimonious interpretation of NRBV evolution. Calibrating the phylogeny allows us to evaluate potential drivers for crested newt cladogenesis. The split between the T. karelinii group and the three other morphotypes, at ca. 10.4 Ma, is associated with the separation of the Balkan and Anatolian landmasses (12-9 Ma). No currently known vicariant events can be ascribed to the other two splits, first at ca. 9.3 Ma, separating T. carnifex - T. macedonicus, and second at ca. 8.8 Ma, splitting T. cristatus and T. dobrogicus. The crested newt morphotypes differ in the duration of their annual aquatic period. We speculate on the role that this ecological differentiation could have played during speciation.

Project description:We explored intraspecific variation in vertebral formulae, more specifically the variation in the number of thoracic vertebrae and frequencies of transitional sacral vertebrae in Triturus newts (Caudata: Salamandridae). Within salamandrid salamanders this monophyletic group shows the highest disparity in the number of thoracic vertebrae and considerable intraspecific variation in the number of thoracic vertebrae. Triturus species also differ in their ecological preferences, from predominantly terrestrial to largely aquatic. Following Geoffroy St. Hilaire's and Darwin's rule which states that structures with a large number of serially homologous repetitive elements are more variable than structures with smaller numbers, we hypothesized that the variation in vertebral formulae increases in more elongated species with a larger number of thoracic vertebrae. We furthermore hypothesized that the frequency of transitional vertebrae will be correlated with the variation in the number of thoracic vertebrae within the species. We also investigated potential effects of species hybridization on the vertebral formula. The proportion of individuals with a number of thoracic vertebrae different from the modal number and the range of variation in number of vertebrae significantly increased in species with a larger number of thoracic vertebrae. Contrary to our expectation, the frequencies of transitional vertebrae were not correlated with frequencies of change in the complete vertebrae number. The frequency of transitional sacral vertebra in hybrids did not significantly differ from that of the parental species. Such a pattern could be a result of selection pressure against transitional vertebrae and/or a bias towards the development of full vertebrae numbers. Although our data indicate relaxed selection for vertebral count changes in more elongated, aquatic species, more data on different selective pressures in species with different numbers of vertebrae in the two contrasting, terrestrial and aquatic environments are needed to test for causality.

Project description:BackgroundIn amphibians, thyroid hormone (TH) has a profound role in cranial development, especially in ossification of the late-appearing bones and remodeling of the skull. In the present study, we explored the influence of TH deficiency on bone ossification and resulting skull shape during the ontogeny of Triturus newt hybrid larvae obtained from interspecific crosses between T. ivanbureschi and T. macedonicus.MethodsLarvae were treated with two concentrations of thiourea (an endocrine disruptor that chemically inhibits synthesis of TH) during the midlarval and late larval periods. Morphological differences of the cranium were assessed at the end of the midlarval period (ontogenetic stage 62) and the metamorphic stage after treatment during the late larval period.ResultsThere was no difference in the ossification level and shape of the skull between the experimental groups (control and two treatment concentrations) at stage 62. During the late larval period and metamorphosis, TH deficit had a significant impact on the level of bone ossification and skull shape with no differences between the two treatment concentrations of thiourea. The most pronounced differences in bone development were: the palatopterygoid failed to disintegrate into the palatal and pterygoid portions, retardation was observed in development of the maxilla, nasal and prefrontal bones and larval organization of the vomer was retained in thiourea-treated larvae.ConclusionsThis implies that deficiency of TH caused retardation in development and arrested metamorphic cranium skeletal reorganization, which resulted in divergent cranial shape compared to the control group. Our results confirmed that skull remodeling and ossification of late-appearing bones is TH-dependent, as in other studied Urodela species. Also, our results indicate that TH plays an important role in the establishment of skull shape during the ontogeny of Triturus newts, especially during the late larval period and metamorphosis, when TH concentrations reach their maximum.

Project description:Taxonomic identifications in some groups of lichen-forming fungi have been challenge largely due to the scarcity of taxonomically relevant features and limitations of morphological and chemical characters traditionally used to distinguish closely related taxa. Delineating species boundaries in closely related species or species complexes often requires a range of multisource data sets and comprehensive analytical methods. Here we aim to examine species boundaries in a group of saxicolous lichen forming fungi, the Aspiciliella intermutans complex (Megasporaceae), widespread mainly in the Mediterranean. We gathered DNA sequences of the nuclear ribosomal internal transcribed spacer (nuITS), the nuclear large subunit (nuLSU), the mitochondrial small subunit (mtSSU) ribosomal DNA, and the DNA replication licensing factor MCM7 from 80 samples mostly from Iran, Caucasia, Greece and eastern Europe. We used a combination of phylogenetic strategies and a variety of empirical, sequence-based species delimitation approaches to infer species boundaries in this group. The latter included: the automatic barcode gap discovery (ABGD), the multispecies coalescent approach *BEAST and Bayesian Phylogenetics and Phylogeography (BPP) program. Different species delimitation scenarios were compared using Bayes factors species delimitation analysis. Furthermore, morphological, chemical, ecological and geographical features of the sampled specimens were examined. Our study uncovered cryptic species diversity in A. intermutans and showed that morphology-based taxonomy may be unreliable, underestimating species diversity in this group of lichens. We identified a total of six species-level lineages in the A. intermutans complex using inferences from multiple empirical operational criteria. We found little corroboration between morphological and ecological features with our proposed candidate species, while secondary metabolite data do not corroborate tree topology. The present study on the A. intermutans species-complex indicates that the genus Aspiciliella, as currently circumscribed, is more diverse in Eurasia than previously expected.

Project description:A statistical framework of model comparison and model validation is essential to resolving the debates over concatenation and coalescent models in phylogenomic data analysis. A set of statistical tests are here applied and developed to evaluate and compare the adequacy of substitution, concatenation, and multispecies coalescent (MSC) models across 47 phylogenomic data sets collected across tree of life. Tests for substitution models and the concatenation assumption of topologically congruent gene trees suggest that a poor fit of substitution models, rejected by 44% of loci, and concatenation models, rejected by 38% of loci, is widespread. Logistic regression shows that the proportions of GC content and informative sites are both negatively correlated with the fit of substitution models across loci. Moreover, a substantial violation of the concatenation assumption of congruent gene trees is consistently observed across six major groups (birds, mammals, fish, insects, reptiles, and others, including other invertebrates). In contrast, among those loci adequately described by a given substitution model, the proportion of loci rejecting the MSC model is 11%, significantly lower than those rejecting the substitution and concatenation models. Although conducted on reduced data sets due to computational constraints, Bayesian model validation and comparison both strongly favor the MSC over concatenation across all data sets; the concatenation assumption of congruent gene trees rarely holds for phylogenomic data sets with more than 10 loci. Thus, for large phylogenomic data sets, model comparisons are expected to consistently and more strongly favor the coalescent model over the concatenation model. We also found that loci rejecting the MSC have little effect on species tree estimation. Our study reveals the value of model validation and comparison in phylogenomic data analysis, as well as the need for further improvements of multilocus models and computational tools for phylogenetic inference. [Bayes factor; Bayesian model validation; coalescent prior; congruent gene trees; independent prior; Metazoa; posterior predictive simulation.].

Project description:BackgroundIf the geographical displacement of one species by another is accompanied by hybridization, mitochondrial DNA can introgress asymmetrically, from the outcompeted species into the invading species, over a large area. We explore this phenomenon using the two parapatric crested newt species, Triturus macedonicus and T. karelinii, distributed on the Balkan Peninsula in south-eastern Europe, as a model.ResultsWe first delimit a ca. 54,000 km(2) area in which T. macedonicus contains T. karelinii mitochondrial DNA. This introgression zone bisects the range of T. karelinii, cutting off a T. karelinii enclave. The high similarity of introgressed mitochondrial DNA haplotypes with those found in T. karelinii suggests a recent transfer across the species boundary. We then use ecological niche modeling to explore habitat suitability of the location of the present day introgression zone under current, mid-Holocene and Last Glacial Maximum conditions. This area was inhospitable during the Last Glacial Maximum for both species, but would have been habitable at the mid-Holocene. Since the mid-Holocene, habitat suitability generally increased for T. macedonicus, whereas it decreased for T. karelinii.ConclusionThe presence of a T. karelinii enclave suggests that T. karelinii was the first to colonize the area where the present day introgression zone is positioned after the Last Glacial Maximum. Subsequently, we propose T. karelinii was outcompeted by T. macedonicus, which captured T. karelinii mitochondrial DNA via introgressive hybridization in the process. Ecological niche modeling suggests that this replacement was likely facilitated by a shift in climate since the mid-Holocene. We suggest that the northwestern part of the current introgression zone was probably never inhabited by T. karelinii itself, and that T. karelinii mitochondrial DNA spread there through T. macedonicus exclusively. Considering the spatial distribution of the introgressed mitochondrial DNA and the signal derived from ecological niche modeling, we do not favor the hypothesis that foreign mitochondrial DNA was pulled into the T. macedonicus range by natural selection.

Dataset Information

Data concatenation, Bayesian concordance and coalescent-based analyses of the species tree for the rapid radiation of Triturus newts.

Publications

Data concatenation, Bayesian concordance and coalescent-based analyses of the species tree for the rapid radiation of Triturus newts.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets