Resolving deep phylogenetic relationships in salamanders: analyses of mitochondrial and nuclear genomic data

Seventy-eight entire mitochondrial genomes and nuclear rRNA genes provide insight into the phylogeny of the hard ticks, particularly the Haemaphysalis species, Africaniella transversale and Robertsicus elaphensis

Hoogstraal and Kim (1985) proposed from morphology, three groups of Haemaphysalis subgenera: (i) the "structurally advanced"; (ii) the "structurally intermediate"; and (iii) the "structurally primitive" subgenera. Nuclear gene phylogenies, however, did not indicate monophyly of these morphological groups but alas, only two mitochondrial (mt) genomes from the "structurally intermediate" subgenera had been sequenced. The phylogeny of Haemaphysalis has not yet been resolved. We aimed to resolve the phylogeny of the genus Haemaphysalis, with respect to the subgenus Alloceraea. We presented 15 newly sequenced and annotated mt genomes from 15 species of ticks, five species of which have not been sequenced before, and four new 18S rRNA and 28S rRNA nuclear gene sequences. Our datasets were constructed from 10 mt protein-coding genes, cox1, and the 18S and 28S nuclear rRNA genes. We found a 132-bp insertion between tRNA-Glu (E) gene and the nad1 gene in the mt genome of Haemaphysalis (Alloceraea) inermis that resembles insertions in H. (Alloceraea) kitaokai and Rhipicephalus (Boophilus) geigyi. Our mt phylogenies had the three species of Amblyomma (Aponomma) we sequenced embedded in the main clade of Amblyomma: Am. (Aponomma) fimbriatum, Am. (Aponomma) gervaisi and Am. (Aponomma) latum. This is further support for the hypothesis that the evolution of eyes appears to have occurred in the most-recent-common-ancestor of Amblyocephalus (i.e. Amblyomminae plus Rhipicephalinae) and that eyes were subsequently lost in the most-recent-common-ancestor of the subgenus Am. (Aponomma). Either Africaniella transversale or Robertsicus elaphensis, or perhaps Af. transversale plus Ro. elaphensis, appear to be the sister-group to the rest of the metastriate Ixodida. Our cox1 phylogenies did not indicate monophyly of the "structurally primitive", "structurally intermediate" nor the "structurally advanced" groups of Haemaphysalis subgenera. Indeed, the subgenus Alloceraea may be the only monophyletic subgenus of the genus Haemaphysalis sequenced thus far. All of our mt genome and cox1 phylogenies had the subgenus Alloceraea in a clade that was separate from the rest of the Haemaphysalis ticks. If Alloceraea is indeed the sister to the rest of the Haemaphysalis subgenera this would resonate with the argument of Hoogstraal and Kim (1985), that Alloceraea was a subgenus of "primitive" Haemaphysalis. Alectorobius capensis from Japan had a higher genetic-identity to A. sawaii, which was also from Japan, than to the A. capensis from South Africa. This indicates that A. capensis from Japan may be a cryptic species with respect to the A. capensis from South Africa.

Comparative osteology of the hynobiid complex Liua-Protohynobius-Pseudohynobius (Amphibia, Urodela): Ⅰ. Cranial anatomy of Pseudohynobius

Hynobiidae are a clade of salamanders that diverged early within the crown radiation and that retain a considerable number of features plesiomorphic for the group. Their evolutionary history is informed by a fossil record that extends to the Middle Jurassic Bathonian time. Our understanding of the evolution within the total group of Hynobiidae has benefited considerably from recent discoveries of stem hynobiids but is constrained by inadequate anatomical knowledge of some extant forms. Pseudohynobius is a derived hynobiid clade consisting of five to seven extant species living endemic to southwestern China. Although this clade has been recognized for over 37 years, osteological details of these extant hynobiids remain elusive, which undoubtedly has contributed to taxonomic controversies over the hynobiid complex Liua-Protohynobius-Pseudohynobius. Here we provide a bone-by-bone study of the cranium in the five extant species of Pseudohynobius (Ps. flavomaculatus, Ps. guizhouensis, Ps. jinfo, Ps. kuankuoshuiensis and Ps. shuichengensis) based on x-ray computer tomography data for 18 specimens. Our results indicate that the cranium in each of these species has a combination of differences in morphology, proportions and articulation patterns in both dermal and endochondral bones. Our study establishes a range of intraspecific differences that will serve as organizing hypotheses for future studies as more extensive collections of these species become available. Morphological features in the cranium for terrestrial ecological adaptation in Hynobiidae are summarized. Based on the results, we also discuss the evolution and development of several potential synapomorphies of Hynobiidae, including features of the orbitosphenoid and articular.

Genome-wide analyses reveal drivers of penguin diversification

Penguins are the only extant family of flightless diving birds. They currently comprise at least 18 species, distributed from polar to tropical environments in the Southern Hemisphere. The history of their diversification and adaptation to these diverse environments remains controversial. We used 22 new genomes from 18 penguin species to reconstruct the order, timing, and location of their diversification, to track changes in their thermal niches through time, and to test for associated adaptation across the genome. Our results indicate that the penguin crown-group originated during the Miocene in New Zealand and Australia, not in Antarctica as previously thought, and that Aptenodytes is the sister group to all other extant penguin species. We show that lineage diversification in penguins was largely driven by changing climatic conditions and by the opening of the Drake Passage and associated intensification of the Antarctic Circumpolar Current (ACC). Penguin species have introgressed throughout much of their evolutionary history, following the direction of the ACC, which might have promoted dispersal and admixture. Changes in thermal niches were accompanied by adaptations in genes that govern thermoregulation and oxygen metabolism. Estimates of ancestral effective population sizes (N _e ) confirm that penguins are sensitive to climate shifts, as represented by three different demographic trajectories in deeper time, the most common (in 11 of 18 penguin species) being an increased N _e between 40 and 70 kya, followed by a precipitous decline during the Last Glacial Maximum. The latter effect is most likely a consequence of the overall decline in marine productivity following the last glaciation.

The choice of tree prior and molecular clock does not substantially affect phylogenetic inferences of diversification rates

Comparative methods allow researchers to make inferences about evolutionary processes and patterns from phylogenetic trees. In Bayesian phylogenetics, estimating a phylogeny requires specifying priors on parameters characterizing the branching process and rates of substitution among lineages, in addition to others. Accordingly, characterizing the effect of prior selection on phylogenies is an active area of research. The choice of priors may systematically bias phylogenetic reconstruction and, subsequently, affect conclusions drawn from the resulting phylogeny. Here, we focus on the impact of priors in Bayesian phylogenetic inference and evaluate how they affect the estimation of parameters in macroevolutionary models of lineage diversification. Specifically, we simulate trees under combinations of tree priors and molecular clocks, simulate sequence data, estimate trees, and estimate diversification parameters (e.g., speciation and extinction rates) from these trees. When substitution rate heterogeneity is large, diversification rate estimates deviate substantially from those estimated under the simulation conditions when not captured by an appropriate choice of relaxed molecular clock. However, in general, we find that the choice of tree prior and molecular clock has relatively little impact on the estimation of diversification rates insofar as the sequence data are sufficiently informative and substitution rate heterogeneity among lineages is low-to-moderate.

A transcriptome-based resolution for a key taxonomic controversy in Cupressaceae

BACKGROUND AND AIMS

Rapid evolutionary divergence and reticulate evolution may result in phylogenetic relationships that are difficult to resolve using small nucleotide sequence data sets. Next-generation sequencing methods can generate larger data sets that are better suited to solving these puzzles. One major and long-standing controversy in conifers concerns generic relationships within the subfamily Cupressoideae (105 species, approx. 1/6 of all conifers) of Cupressaceae, in particular the relationship between Juniperus, Cupressus and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) clade. Here we attempt to resolve this question using transcriptome-derived data.

METHODS

Transcriptome sequences of 20 species from Cupressoideae were collected. Using MarkerMiner, single-copy nuclear (SCN) genes were extracted. These were applied to estimate phylogenies based on concatenated data, species trees and a phylogenetic network. We further examined the effect of alternative backbone topologies on downstream analyses, including biogeographic inference and dating analysis.

RESULTS

Based on the 73 SCN genes (>200 000 bp total alignment length) we considered, all tree-building methods lent strong support for the relationship (HCX, (Juniperus, Cupressus)); however, strongly supported conflicts among individual gene trees were also detected. Molecular dating suggests that these three lineages shared a most recent common ancestor approx. 60 million years ago (Mya), and that Juniperus and Cupressus diverged about 56 Mya. Ancestral area reconstructions (AARs) suggest an Asian origin for the entire clade, with subsequent dispersal to North America, Europe and Africa.

CONCLUSIONS

Our analysis of SCN genes resolves a controversial phylogenetic relationship in the Cupressoideae, a major clade of conifers, and suggests that rapid evolutionary divergence and incomplete lineage sorting probably acted together as the source for conflicting phylogenetic inferences between gene trees and between our robust results and recently published studies. Our updated backbone topology has not substantially altered molecular dating estimates relative to previous studies; however, application of the latest AAR approaches has yielded a clearer picture of the biogeographic history of Cupressoideae.

Incomplete lineage sorting rather than hybridization explains the inconsistent phylogeny of the wisent

The wisent or European bison is the largest European herbivore and is completely cross-fertile with its American relative. However, mtDNA genome of wisent is similar to that of cattle, which suggests that wisent emerged as a hybrid of bison and an extinct cattle-like species. Here, we analyzed nuclear whole-genome sequences of the bovine species, and found only a minor and recent gene flow between wisent and cattle. Furthermore, we identified an appreciable heterogeneity of the nuclear gene tree topologies of the bovine species. The relative frequencies of various topologies, including the mtDNA topology, were consistent with frequencies of incomplete lineage sorting (ILS) as estimated by tree coalescence analysis. This indicates that ILS has occurred and may well account for the anomalous wisent mtDNA phylogeny as the outcome of a rare event. We propose that ILS is a possible explanation of phylogenomic anomalies among closely related species.

Kun Wang et al. present a genomic analysis identifying incomplete lineage sorting and hybridization in the mitochondrial DNA of the European bison (wisent). They find that incomplete lineage sorting is the most feasible explanation for the phylogenetic heterogeneity observed in Bovidae.

A new hynobiid-like salamander (Amphibia, Urodela) from Inner Mongolia, China, provides a rare case study of developmental features in an Early Cretaceous fossil urodele

A new fossil salamander, Nuominerpeton aquilonaris (gen. et sp. nov.), is named and described based on specimens from the Lower Cretaceous Guanghua Formation of Inner Mongolia, China. The new discovery documents a far northern occurrence of Early Cretaceous salamanders in China, extending the geographic distribution for the Mesozoic fossil record of the group from the Jehol area (40th-45th parallel north) to near the 49th parallel north. The new salamander is characterized by having the orbitosphenoid semicircular in shape; coracoid plate of the scapulocoracoid greatly expanded with a convex ventral and posterior border; ossification of two centralia in carpus and tarsus; and first digit being about half the length of the second digit in both manus and pes. The new salamander appears to be closely related to hynobiids, although this inferred relationship awaits confirmation by research in progress by us on a morphological and molecular combined analysis of cryptobranchoid relationships. Comparison of adult with larval and postmetamorphic juvenile specimens provides insights into developmental patterns of cranial and postcranial skeletons in this fossil species, especially resorption of the palatine and anterior portions of the palatopterygoid in the palate and the coronoid in the mandible during metamorphosis, and postmetamorphic ossification of the mesopodium in both manus and pes. Thus, this study provides a rare case study of developmental features in a Mesozoic salamander.

The first endemic West African vertebrate family - a new anuran family highlighting the uniqueness of the Upper Guinean biodiversity hotspot

BACKGROUND

Higher-level systematics in amphibians is relatively stable. However, recent phylogenetic studies of African torrent-frogs have uncovered high divergence in these phenotypically and ecologically similar frogs, in particular between West African torrent-frogs versus Central (Petropedetes) and East African (Arthroleptides and Ericabatrachus) lineages. Because of the considerable molecular divergence, and external morphology of the single West African torrent-frog species a new genus was erected (Odontobatrachus). In this study we aim to clarify the systematic position of West African torrent-frogs (Odontobatrachus). We determine the relationships of torrent-frogs using a multi-locus, nuclear and mitochondrial, dataset and include genera of all African and Asian ranoid families. Using micro-tomographic scanning we examine osteology and external morphological features of West African torrent-frogs to compare them with other ranoids.

RESULTS

Our analyses reveal Petropedetidae (Arthroleptides, Ericabatrachus, Petropedetes) as the sister taxon of the Pyxicephalidae. The phylogenetic position of Odontobatrachus is clearly outside Petropedetidae, and not closely related to any other ranoid family. According to our time-tree estimation Odontobatrachus has been separated from other frog lineages since the Cretaceous (90.1 Ma; confidence interval: 84.2-97.1 Ma). Along with this molecular evidence, osteological and external diagnostic characters recognize West African torrent-frogs as distinct from other ranoids and provide strong support for the necessity of the recognition of a new family of frogs. This is the only endemic vertebrate family occurring in the Upper Guinea biodiversity hotspot.

CONCLUSION

Based on molecular and morphological distinctiveness, the West African torrent-frog Odontobatrachus natator is allocated to a newly described anuran family. The discovery of an endemic vertebrate family in West Africa highlights the Upper Guinean forests as an outstanding, but highly endangered biodiversity hotspot.

A versatile and highly efficient toolkit including 102 nuclear markers for vertebrate phylogenomics, tested by resolving the higher level relationships of the caudata

Resolving difficult nodes for any part of the vertebrate tree of life often requires analyzing a large number of loci. Developing molecular markers that are workable for the groups of interest is often a bottleneck in phylogenetic research. Here, on the basis of a nested polymerase chain reaction (PCR) strategy, we present a universal toolkit including 102 nuclear protein-coding locus (NPCL) markers for vertebrate phylogenomics. The 102 NPCL markers have a broad range of evolutionary rates, which makes them useful for a wide range of time depths. The new NPCL toolkit has three important advantages compared with all previously developed NPCL sets: 1) the kit is universally applicable across vertebrates, with a PCR success rate of 94.6% in 16 widely divergent tested vertebrate species; 2) more than 90% of PCR reactions produce strong and single bands of the expected sizes that can be directly sequenced; and 3) all cleanup PCR reactions can be sequenced with only two specific universal primers. To test its actual phylogenetic utility, 30 NPCLs from this toolkit were used to address the higher level relationships of living salamanders. Of the 639 target PCR reactions performed on 19 salamanders and several outgroup species, 632 (98.9%) were successful, and 602 (94.1%) were directly sequenced. Concatenation and species-tree analyses on this 30-locus data set produced a fully resolved phylogeny and showed that Cryptobranchoidea (Cryptobranchidae + Hynobiidae) branches first within the salamander tree, followed by Sirenidae. Our experimental tests and our demonstration for a particular case show that our NPCL toolkit is a highly reliable, fast, and cost-effective approach for vertebrate phylogenomic studies and thus has the potential to accelerate the completion of many parts of the vertebrate tree of life.

Mitochondrial genomes and divergence times of crocodile newts: inter-islands distribution of Echinotriton andersoni and the origin of a unique repetitive sequence found in Tylototriton mt genomes

Crocodile newts, which constitute the genera Echinotriton and Tylototriton, are known as living fossils, and these genera comprise many endangered species. To identify mitochondrial (mt) genes suitable for future population genetic analyses for endangered taxa, we determined the complete nucleotide sequences of the mt genomes of the Japanese crocodile newt Echinotriton andersoni and Himalayan crocodile newt Tylototriton verrucosus. Although the control region (CR) is known as the most variable mtDNA region in many animal taxa, the CRs of crocodile newts are highly conservative. Rather, the genes of NADH dehydrogenase subunits and ATPase subunit 6 were found to have high sequence divergences and to be usable for population genetics studies. To estimate the inter-population divergence ages of E. andersoni endemic to the Ryukyu Islands, we performed molecular dating analysis using whole and partial mt genomic data. The estimated divergence ages of the inter-island individuals are older than the paleogeographic segmentation ages of the islands, suggesting that the lineage splits of E. andersoni populations were not caused by vicariant events. Our phylogenetic analysis with partial mt sequence data also suggests the existence of at least two more undescribed species in the genus Tylototriton. We also found unusual repeat sequences containing the 3' region of cytochrome apoenzyme b gene, whole tRNA-Thr gene, and a noncoding region (the T-P noncoding region characteristic in caudate mtDNAs) from T. verrucosus mtDNA. Similar repeat sequences were found in two other Tylototriton species. The Tylototriton taxa with the repeats become a monophyletic group, indicating a single origin of the repeat sequences. The intra-and inter-specific comparisons of the repeat sequences suggest the occurrences of homologous recombination-based concerted evolution among the repeat sequences.

Late Jurassic salamandroid from western Liaoning, China

A Jurassic salamander, Beiyanerpeton jianpingensis (gen. et sp. nov.), from a recently found site in western Liaoning Province, China is the earliest known record of Salamandroidea. As a Late Jurassic record of the group, it extends the range of the clade by ~40 Ma. The Late Jurassic taxon is neotenic and represented by exceptionally preserved specimens, including fully articulated cranial and postcranial skeletons and bony gill structures close to the cheek region. The fossil beds, consisting of dark-brown volcanic ash shales of the Upper Jurassic Tiaojishan (Lanqi) Formation (Oxfordian), underlie trachyandesite rocks that have yielded a SHRIMP zircon U-Pb date of 157 ± 3 Ma. The fossiliferous beds are substantially older than the Jehol Group, including the Yixian Formation ((40)Ar/(39)Ar dates of 122-129 Ma), but slightly younger than the Middle Jurassic Daohugou horizon ((40)Ar/(39)Ar date of 164 ± 4 Ma). The early fossil taxon shares with extant salamandroids derived character states, including: separated nasals lacking a midline contact, angular fused to the prearticular in the lower jaw, and double-headed ribs on the presacral vertebrae. In contrast to extant salamandroids, however, the salamander shows a discrete and tooth-bearing palatine, and unequivocally nonpedicellate and monocuspid marginal teeth in large and presumably mature individuals. The finding provides insights into the evolution of key characters of salamanders, and also provides direct evidence supporting the hypothesis that the split between Cryptobranchoidea and Salamandroidea had taken placed before the Late Jurassic Oxfordian time. In this aspect, both paleontological and molecular data now come to agree.

Analyzing the relationship between sequence divergence and nodal support using Bayesian phylogenetic analyses

Determining the appropriate gene for phylogeny reconstruction can be a difficult process. Rapidly evolving genes tend to resolve recent relationships, but suffer from alignment issues and increased homoplasy among distantly related species. Conversely, slowly evolving genes generally perform best for deeper relationships, but lack sufficient variation to resolve recent relationships. We determine the relationship between sequence divergence and Bayesian phylogenetic reconstruction ability using both natural and simulated datasets. The natural data are based on 28 well-supported relationships within the subphylum Vertebrata. Sequences of 12 genes were acquired and Bayesian analyses were used to determine phylogenetic support for correct relationships. Simulated datasets were designed to determine whether an optimal range of sequence divergence exists across extreme phylogenetic conditions. Across all genes we found that an optimal range of divergence for resolving the correct relationships does exist, although this level of divergence expectedly depends on the distance metric. Simulated datasets show that an optimal range of sequence divergence exists across diverse topologies and models of evolution. We determine that a simple to measure property of genetic sequences (genetic distance) is related to phylogenic reconstruction ability in Bayesian analyses. This information should be useful for selecting the most informative gene to resolve any relationships, especially those that are difficult to resolve, as well as minimizing both cost and confounding information during project design.

Data mining approach identifies research priorities and data requirements for resolving the red algal tree of life

BACKGROUND

The assembly of the tree of life has seen significant progress in recent years but algae and protists have been largely overlooked in this effort. Many groups of algae and protists have ancient roots and it is unclear how much data will be required to resolve their phylogenetic relationships for incorporation in the tree of life. The red algae, a group of primary photosynthetic eukaryotes of more than a billion years old, provide the earliest fossil evidence for eukaryotic multicellularity and sexual reproduction. Despite this evolutionary significance, their phylogenetic relationships are understudied. This study aims to infer a comprehensive red algal tree of life at the family level from a supermatrix containing data mined from GenBank. We aim to locate remaining regions of low support in the topology, evaluate their causes and estimate the amount of data required to resolve them.

RESULTS

Phylogenetic analysis of a supermatrix of 14 loci and 98 red algal families yielded the most complete red algal tree of life to date. Visualization of statistical support showed the presence of five poorly supported regions. Causes for low support were identified with statistics about the age of the region, data availability and node density, showing that poor support has different origins in different parts of the tree. Parametric simulation experiments yielded optimistic estimates of how much data will be needed to resolve the poorly supported regions (ca. 103 to ca. 104 nucleotides for the different regions). Nonparametric simulations gave a markedly more pessimistic image, some regions requiring more than 2.8 105 nucleotides or not achieving the desired level of support at all. The discrepancies between parametric and nonparametric simulations are discussed in light of our dataset and known attributes of both approaches.

CONCLUSIONS

Our study takes the red algae one step closer to meaningful inclusion in the tree of life. In addition to the recovery of stable relationships, the recognition of five regions in need of further study is a significant outcome of this work. Based on our analyses of current availability and future requirements of data, we make clear recommendations for forthcoming research.

Evolution of gigantism in amphiumid salamanders

The Amphiumidae contains three species of elongate, permanently aquatic salamanders with four diminutive limbs that append one, two, or three toes. Two of the species, Amphiuma means and A. tridactylum, are among the largest salamanders in the world, reaching lengths of more than one meter, whereas the third species (A. pholeter), extinct amphiumids, and closely related salamander families are relatively small. Amphiuma means and A. tridactylum are widespread species and live in a wide range of lowland aquatic habitats on the Coastal Plain of the southeastern United States, whereas A. pholeter is restricted to very specialized organic muck habitats and is syntopic with A. means. Here we present analyses of sequences of mitochondrial and nuclear loci from across the distribution of the three taxa to assess lineage diversity, relationships, and relative timing of divergence in amphiumid salamanders. In addition we analyze the evolution of gigantism in the clade. Our analyses indicate three lineages that have diverged since the late Miocene, that correspond to the three currently recognized species, but the two gigantic species are not each other's closest relatives. Given that the most closely related salamander families and fossil amphiumids from the Upper Cretaceous and Paleocene are relatively small, our results suggest at least two extreme changes in body size within the Amphuimidae. Gigantic body size either evolved once as the ancestral condition of modern amphiumas, with a subsequent strong size reduction in A. pholeter, or gigantism independently evolved twice in the modern species, A. means and A. tridactylum. These patterns are concordant with differences in habitat breadth and range size among lineages, and have implications for reproductive isolation and diversification of amphiumid salamanders.

Island biogeography of Galápagos lava lizards (Tropiduridae: Microlophus): species diversity and colonization of the archipelago

The "lava lizards" (Microlophus) are distributed throughout the Galápagos Archipelago, and consist of radiations derived from two independent colonizations. The "Eastern Radiation" includes M. bivittatus and M. habeli endemic to San Cristobal and Marchena Islands. The "Western Radiation" includes five to seven historically recognized species distributed across almost the entire Archipelago. We combine dense geographic sampling and multilocus sequence data to estimate a phylogenetic hypothesis for the Western Radiation, to delimit species boundaries in this radiation, and to estimate a time frame for colonization events. Our phylogenetic hypothesis rejects two earlier topologies for the Western Radiation and paraphyly of M. albemarlensis, while providing strong support for single colonizations on each island. The colonization history implied by our phylogeny is consistent with general expectations of an east-to-west route predicted by the putative age of island groups, and prevailing ocean currents in the Archipelago. Additionally, combined evidence suggests that M. indefatigabilis from Santa Fe should be recognized as a full species. Finally, molecular divergence estimates suggest that the two colonization events likely occurred on the oldest existing islands, and the Western Radiation represents a recent radiation that, in most cases, has produced species that are considerably younger than the islands they inhabit.

The genus Coleodactylus (Sphaerodactylinae, Gekkota) revisited: a molecular phylogenetic perspective

Nucleotide sequence data from a mitochondrial gene (16S) and two nuclear genes (c-mos, RAG-1) were used to evaluate the monophyly of the genus Coleodactylus, to provide the first phylogenetic hypothesis of relationships among its species in a cladistic framework, and to estimate the relative timing of species divergences. Maximum Parsimony, Maximum Likelihood and Bayesian analyses of the combined data sets retrieved Coleodactylus as a monophyletic genus, although weakly supported. Species were recovered as two genetically and morphological distinct clades, with C. amazonicus populations forming the sister taxon to the meridionalis group (C. brachystoma, C. meridionalis, C. natalensis, and C. septentrionalis). Within this group, C. septentrionalis was placed as the sister taxon to a clade comprising the rest of the species, C. meridionalis was recovered as the sister species to C. brachystoma, and C. natalensis was found nested within C. meridionalis. Divergence time estimates based on penalized likelihood and Bayesian dating methods do not support the previous hypothesis based on the Quaternary rain forest fragmentation model proposed to explain the diversification of the genus. The basal cladogenic event between major lineages of Coleodactylus was estimated to have occurred in the late Cretaceous (72.6+/-1.77 Mya), approximately at the same point in time than the other genera of Sphaerodactylinae diverged from each other. Within the meridionalis group, the split between C. septentrionalis and C. brachystoma+C. meridionalis was placed in the Eocene (46.4+/-4.22 Mya), and the divergence between C. brachystoma and C. meridionalis was estimated to have occurred in the Oligocene (29.3+/-4.33 Mya). Most intraspecific cladogenesis occurred through Miocene to Pliocene, and only for two conspecific samples and for C. natalensis could a Quaternary differentiation be assumed (1.9+/-1.3 Mya).

Development and application of a phylogenomic toolkit: resolving the evolutionary history of Madagascar's lemurs

Lemurs and the other strepsirrhine primates are of great interest to the primate genomics community due to their phylogenetic placement as the sister lineage to all other primates. Previous attempts to resolve the phylogeny of lemurs employed limited mitochondrial or small nuclear data sets, with many relationships poorly supported or entirely unresolved. We used genomic resources to develop 11 novel markers from nine chromosomes, representing approximately 9 kb of nuclear sequence data. In combination with previously published nuclear and mitochondrial loci, this yields a data set of more than 16 kb and adds approximately 275 kb of DNA sequence to current databases. Our phylogenetic analyses confirm hypotheses of lemuriform monophyly and provide robust resolution of the phylogenetic relationships among the five lemuriform families. We verify that the genus Daubentonia is the sister lineage to all other lemurs. The Cheirogaleidae and Lepilemuridae are sister taxa and together form the sister lineage to the Indriidae; this clade is the sister lineage to the Lemuridae. Divergence time estimates indicate that lemurs are an ancient group, with their initial diversification occurring around the Cretaceous-Tertiary boundary. Given the power of this data set to resolve branches in a notoriously problematic area of primate phylogeny, we anticipate that our phylogenomic toolkit will be of value to other studies of primate phylogeny and diversification. Moreover, the methods applied will be broadly applicable to other taxonomic groups where phylogenetic relationships have been notoriously difficult to resolve.

Mitochondrial matR sequences help to resolve deep phylogenetic relationships in rosids

Background

Rosids are a major clade in the angiosperms containing 13 orders and about one-third of angiosperm species. Recent molecular analyses recognized two major groups (i.e., fabids with seven orders and malvids with three orders). However, phylogenetic relationships within the two groups and among fabids, malvids, and potentially basal rosids including Geraniales, Myrtales, and Crossosomatales remain to be resolved with more data and a broader taxon sampling. In this study, we obtained DNA sequences of the mitochondrial matR gene from 174 species representing 72 families of putative rosids and examined phylogenetic relationships and phylogenetic utility of matR in rosids. We also inferred phylogenetic relationships within the "rosid clade" based on a combined data set of 91 taxa and four genes including matR, two plastid genes (rbcL, atpB), and one nuclear gene (18S rDNA).

Results

Comparison of mitochondrial matR and two plastid genes (rbcL and atpB) showed that the synonymous substitution rate in matR was approximately four times slower than those of rbcL and atpB; however, the nonsynonymous substitution rate in matR was relatively high, close to its synonymous substitution rate, indicating that the matR has experienced a relaxed evolutionary history. Analyses of our matR sequences supported the monophyly of malvids and most orders of the rosids. However, fabids did not form a clade; instead, the COM clade of fabids (Celastrales, Oxalidales, Malpighiales, and Huaceae) was sister to malvids. Analyses of the four-gene data set suggested that Geraniales and Myrtales were successively sister to other rosids, and that Crossosomatales were sister to malvids.

Conclusion

Compared to plastid genes such as rbcL and atpB, slowly evolving matR produced less homoplasious but not less informative substitutions. Thus, matR appears useful in higher-level angiosperm phylogenetics. Analysis of matR alone identified a novel deep relationship within rosids, the grouping of the COM clade of fabids and malvids, which was not resolved by any previous molecular analyses but recently suggested by floral structural features. Our four-gene analysis supported the placements of Geraniales, Myrtales at basal nodes of the rosid clade and placed Crossosomatales as sister to malvids. We also suggest that the core part of rosids should include fabids, malvids and Crossosomatales.

Data congruence, paedomorphosis and salamanders

BACKGROUND

The retention of ancestral juvenile characters by adult stages of descendants is called paedomorphosis. However, this process can mislead phylogenetic analyses based on morphological data, even in combination with molecular data, because the assessment if a character is primary absent or secondary lost is difficult. Thus, the detection of incongruence between morphological and molecular data is necessary to investigate the reliability of simultaneous analyses. Different methods have been proposed to detect data congruence or incongruence. Five of them (PABA, PBS, NDI, LILD, DRI) are used herein to assess incongruence between morphological and molecular data in a case study addressing salamander phylogeny, which comprises several supposedly paedomorphic taxa. Therefore, previously published data sets were compiled herein. Furthermore, two strategies ameliorating effects of paedomorphosis on phylogenetic studies were tested herein using a statistical rigor. Additionally, efficiency of the different methods to assess incongruence was analyzed using this empirical data set. Finally, a test statistic is presented for all these methods except DRI.

RESULTS

The addition of morphological data to molecular data results in both different positions of three of the four paedomorphic taxa and strong incongruence, but treating the morphological data using different strategies ameliorating the negative impact of paedomorphosis revokes these changes and minimizes the conflict. Of these strategies the strategy to just exclude paedomorphic character traits seem to be most beneficial. Of the three molecular partitions analyzed herein the RAG1 partition seems to be the most suitable to resolve deep salamander phylogeny. The rRNA and mtDNA partition are either too conserved or too variable, respectively. Of the different methods to detect incongruence, the NDI and PABA approaches are more conservative in the indication of incongruence than LILD and PBS.

CONCLUSION

Paedomorphosis induces strong conflicts and can mislead the phylogenetic analyses even in combined analyses. However, different strategies are efficiently minimizing these problems. Though the exploration of different methods to detect incongruence is preferable NDI and PABA are more conservative than the others and NDI is computational less extensive than PABA.

Global Patterns of Diversification and Species Richness in Amphibians

Geographic patterns of species richness ultimately arise through the processes of speciation, extinction, and dispersal, but relatively few studies consider evolutionary and biogeographic processes in explaining these diversity patterns. One explanation for high tropical species richness is that many species-rich clades originated in tropical regions and spread to temperate regions infrequently and more recently, leaving little time for species richness to accumulate there (assuming similar rates of diversification in temperate and tropical regions). However, the major clades of anurans (frogs) and salamanders may offer a compelling counterexample. Most salamander families are predominately temperate in distribution, but the one primarily tropical clade (Bolitoglossinae) contains nearly half of all salamander species. Similarly, most basal clades of anurans are predominately temperate, but one largely tropical clade (Neobatrachia) contains approximately 96% of anurans. In this article, I examine patterns of diversification in frogs and salamanders and their relationship to large-scale patterns of species richness in amphibians. I find that diversification rates in both frogs and salamanders increase significantly with decreasing latitude. These results may shed light on both the evolutionary causes of the latitudinal diversity gradient and the dramatic but poorly explained disparities in the diversity of living amphibian clades.

Profiling phylogenetic informativeness

The resolution of four controversial topics in phylogenetic experimental design hinges upon the informativeness of characters about the historical relationships among taxa. These controversies regard the power of different classes of phylogenetic character, the relative utility of increased taxonomic versus character sampling, the differentiation between lack of phylogenetic signal and a historical rapid radiation, and the design of taxonomically broad phylogenetic studies optimized by taxonomically sparse genome-scale data. Quantification of the informativeness of characters for resolution of phylogenetic hypotheses during specified historical epochs is key to the resolution of these controversies. Here, such a measure of phylogenetic informativeness is formulated. The optimal rate of evolution of a character to resolve a dated four-taxon polytomy is derived. By scaling the asymptotic informativeness of a character evolving at a nonoptimal rate by the derived asymptotic optimum, and by normalizing so that net phylogenetic informativeness is equivalent for all rates when integrated across all of history, an informativeness profile across history is derived. Calculation of the informativeness per base pair allows estimation of the cost-effectiveness of character sampling. Calculation of the informativeness per million years allows comparison across historical radiations of the utility of a gene for the inference of rapid adaptive radiation. The theory is applied to profile the phylogenetic informativeness of the genes BRCA1, RAG1, GHR, and c-myc from a muroid rodent sequence data set. Bounded integrations of the phylogenetic profile of these genes over four epochs comprising the diversifications of the muroid rodents, the mammals, the lobe-limbed vertebrates, and the early metazoans demonstrate the differential power of these genes to resolve the branching order among ancestral lineages. This measure of phylogenetic informativeness yields a new kind of information for evaluation of phylogenetic experiments. It conveys the utility of the addition of characters a phylogenetic study and it provides a basis for deciding whether appropriate phylogenetic power has been applied to a polytomy that is proposed to be a rapid radiation. Moreover, it provides a quantitative measure of the capacity of a gene to resolve soft polytomies.

A Bayesian approach on molecules and behavior: reconsidering phylogenetic and evolutionary patterns of the Salamandridae with emphasis on Triturus newts

The monophyly of European newts of the genus Triturus within the family Salamandridae has for decades rested on presumably homologous behavioral and morphological characters. Molecular data challenge this hypothesis, but the phylogenetic position of Triturus within the Salamandridae has not yet been convincingly resolved. We addressed this issue and the temporal divergence of Triturus within the Salamandridae with novel Bayesian approaches applied to DNA sequence data from three mitochondrial genes (12S, 16S and cytb). We included 38 salamandrid species comprising all 13 recognized species of Triturus and 16 out of 17 salamandrid genera. A clade comprising all the "Newts" can be separated from the "True Salamanders" and Salamandrina clades. Within the "Newts" well-supported clades are: Tylototriton-Pleurodeles, the "New World Newts" (Notophthalmus-Taricha), and the "Modern Eurasian Newts" (Cynops, Pachytriton, Paramesotriton=together the "Modern Asian Newts", Calotriton, Euproctus, Neurergus and Triturus species). We found that Triturus is a non-monophyletic species assemblage, which includes four groups that are themselves monophyletic: (i) the "Large-Bodied Triturus" (six species), (ii) the "Small-Bodied Triturus" (five species), (iii) T. alpestris and (iv) T. vittatus. We estimated that the last common ancestor of Triturus existed around 64 million years ago (mya) while the root of the Salamandridae dates back to 95 mya. This was estimated using a fossil-based molecular dating approach and an explicit framework to select calibration points that least underestimated their corresponding nodes. Using the molecular phylogeny we mapped the evolution of life history and courtship traits in Triturus and found that several Triturus-specific courtship traits evolved independently.

Phylogenetic position of tetraodontiform fishes within the higher teleosts: Bayesian inferences based on 44 whole mitochondrial genome sequences

Tetraodontiformes includes approximately 350 species assigned to nine families, sharing several reduced morphological features of higher teleosts. The order has been accepted as a monophyletic group by many authors, although several alternative hypotheses exist regarding its phylogenetic position within the higher teleosts. To date, acanthuroids, zeiforms, and lophiiforms have been proposed as sister-groups of the tetraodontiforms. The monophyly and sister-group status was investigated using whole mitochondrial genome (mitogenome) sequences from 44 purposefully-chosen species (26 sequences newly-determined during the study) that fully represent the major tetraodontiform lineages plus all the groups that have been hypothesized as being close relatives. Partitioned Bayesian analyses were conducted with the three datasets that comprised concatenated nucleotide sequences from 13 protein-coding genes (with and without, or with RY-coding, 3rd codon positions), plus 22 transfer RNA and two ribosomal RNA genes. The resultant trees were well resolved and largely congruent, with most internal branches being supported by high posterior probabilities. Mitogenomic data strongly supported the monophyly of tetraodontiform fishes, placing them as a sister-group of either Lophiiformes plus Caproidei or Caproidei only. The sister-group relationship between Acanthuroidei and Tetraodontiformes was statistically rejected using Bayes factors. These results were confirmed by a reanalysis of the previously published nuclear RAG1 gene sequences using the Bayesian method. Within the Tetraodontiformes, however, monophylies of the three superfamilies were not recovered and further taxonomic sampling and subsequent efforts should clarify these relationships.

A molecular assessment of phylogenetic relationships and lineage accumulation rates within the family Salamandridae (Amphibia, Caudata)

We examine phylogenetic relationships among salamanders of the family Salamandridae using approximately 2700 bases of new mtDNA sequence data (the tRNALeu, ND1, tRNAIle, tRNAGln, tRNAMet, ND2, tRNATrp, tRNAAla, tRNAAsn, tRNACys, tRNATyr, and COI genes and the origin for light-strand replication) collected from 96 individuals representing 61 of the 66 recognized salamandrid species and outgroups. Phylogenetic analyses using maximum parsimony and Bayesian analysis are performed on the new data alone and combined with previously reported sequences from other parts of the mitochondrial genome. The basal phylogenetic split is a polytomy of lineages ancestral to (1) the Italian newt Salamandrina terdigitata, (2) a strongly supported clade comprising the "true" salamanders (genera Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), and (3) a strongly supported clade comprising all newts except S. terdigitata. Strongly supported clades within the true salamanders include monophyly of each genus and grouping Chioglossa and Mertensiella as the sister taxon to a clade comprising Lyciasalamandra and Salamandra. Among newts, genera Echinotriton, Pleurodeles, and Tylototriton form a strongly supported clade whose sister taxon comprises the genera Calotriton, Cynops, Euproctus, Neurergus, Notophthalmus, Pachytriton, Paramesotriton, Taricha, and Triturus. Our results strongly support monophyly of all polytypic newt genera except Paramesotriton and Triturus, which appear paraphyletic, and Calotriton, for which only one of the two species is sampled. Other well-supported clades within newts include (1) Asian genera Cynops, Pachytriton, and Paramesotriton, (2) North American genera Notophthalmus and Taricha, (3) the Triturus vulgaris species group, and (4) the Triturus cristatus species group; some additional groupings appear strong in Bayesian but not parsimony analyses. Rates of lineage accumulation through time are evaluated using this nearly comprehensive sampling of salamandrid species-level lineages. Rate of lineage accumulation appears constant throughout salamandrid evolutionary history with no obvious fluctuations associated with origins of morphological or ecological novelties.

Molecular phylogenetics of the Anolis onca series: a case history in retrograde evolution revisited

Anoles of the Anolis onca series represent a dramatic case of retrograde evolution, exhibiting great reduction (A. annectens) and loss (A. onca) of the subdigital pads considered a key innovation for the evolutionary radiation of anoles in arboreal environments. We present a molecular phylogenetic analysis of these anoles and their closest known relatives (A. auratus, A. lineatus, A. meridionalis, and A. nitens) using new mitochondrial DNA sequence data from the ND2 gene, five tRNA genes (tRNA(Trp), tRNA(Ala), tRNA(Asn), tRNA(Cys), tRNA(Tyr)), the origin of light-strand replication, and a portion of the CO1 gene (1,446 aligned base positions, 612 parsimony informative). Our results confirm monophyly of the A. onca series and suggest an evolutionary separation of approximately 10 million years between A. annectens and A. onca. Evolution of subdigital structure in this series illustrates ectopic expression of developmental programs that replace flexible subdigital lamellae of the toepad with rigid, keeled scales resembling dorsal digital scales. Our phylogenetic results indicate that narrowing of the toepad in A. auratus evolved separately from toepad reduction in the A. onca series. Expansion of the subdigital lamellae along the phalanges in A. auratus appears to compensate constriction of lamellae by digital narrowing, maintaining greater climbing capability in this species. Toepad evolution in the lineage ancestral to A. auratus features changes of the same developmental modules as the A. onca series but in the opposite direction. Large molecular distances between geographic populations of A. auratus indicate that its derived toepad structure is at least 9 million years old.

Molecular phylogeny and phylogeography of the Australian Diplodactylus stenodactylus (Gekkota; Reptilia) species-group based on mitochondrial and nuclear genes reveals an ancient split between Pilbara and non-Pilbara D. stenodactylus

There is a paucity of research on intra-specific morphological and genetic diversity in Australian arid-zone reptiles, and a number of Australian reptile species have for many years been regarded as "species complexes" that classical morphological analyses could not resolve. We conducted a phylogenetic and phylogeographic study of a widespread species group of Australian geckonid lizards to address two main aims. First, based on a large mitochondrial and nuclear gene data set, we have generated the first molecular phylogeny for the Diplodactylus stenodactylus species group (D. alboguttatus, D. damaeus, D. maini, and D. squarrosus, D. stenodactylus) and multiple outgroups to examine the evolutionary relationships among these arid-zone species and phylogenetic patterns within some species. The edited alignment of 41 individuals comprises 2485 characters (1163 ND2+tRNAs; 490 16s; 832 RAG-1), and of these 717 (29%) were variable and parsimony informative (568 ND2+tRNAs; 89 16s; 60 RAG-1). This broad-scale, multi-gene phylogeny has supported previous conjectures on the higher-level phylogenetic relationships among members of the D. stenodactylus species-group based on morphology, but also has uncovered hidden diversity within the group with two new species identified. Analysis at this broad level has identified patterns associated with the distribution of the D. stenodactylus species group that appear to be influenced by environmental processes operating at large geographic scales. Two major clades within the species group were associated with broad differences in habitat types, with one group largely restricted to the temperate zone of the Southwest Province and another largely restricted to central and northern Western Australia north of Kalgoorlie, in line with the Eremaean Province of the Eremaean Zone and the Northern Province of the Tropical Zone. Second, we have assembled phylogeographic data based on a mitochondrial gene (ND2+tRNAs) for five species (Rhynchoedura ornata, Diplodactylus maini, D. pulcher, D. squarrosus, D. stenodactylus) where larger sampling is available, with particular focus on D. stenodactylus, which is distributed both in the iconic but little-known Pilbara area of endemism in north-western Australia as well as in other parts of the Australian arid zone. The edited alignment of 95 individuals comprises 1142 characters and of these 601 (53%) are variable and parsimony informative. We found significant intra-specific genetic variation in all five species, highlighting the need for large-scale screening of cryptic species, with sampling specifically targeted at determining the geographic limits of such taxa. In particular, within D. stenodactylus, a deep and ancient phylogenetic split distinguishes populations in the Pilbara region from non-Pilbara populations. This split may be the result of broad differences in underlying geological substrate, with the Pilbara clade generally preferring harder soils and the non-Pilbara clade adapted to sand.

Phylogeny, evolution, and biogeography of Asiatic Salamanders (Hynobiidae)

We sequenced 15 complete mitochondrial genomes and performed comprehensive molecular phylogenetic analyses to study the origin and phylogeny of the Hynobiidae, an ancient lineage of living salamanders. Our phylogenetic analyses show that the Hynobiidae is a clade with well resolved relationships, and our results contrast with a morphology-based phylogenetic hypothesis. These salamanders have low vagility and are limited in their distribution primarily by deserts, mountains, and oceans. Our analysis suggests that the relationships among living hynobiids have been shaped primarily by geography. We show that four-toed species assigned to Batrachuperus do not form a monophyletic group, and those that occur in Afghanistan and Iran are transferred to the resurrected Paradactylodon. Convergent morphological characters in different hynobiid lineages are likely produced by similar environmental selective pressures. Clock-independent molecular dating suggests that hynobiids originated in the Middle Cretaceous [ approximately 110 million years ago (Mya)]. We propose an "out of North China" hypothesis for hynobiid origins and hypothesize an ancestral stream-adapted form. Given the particular distributional patterns and our molecular dating estimates, we hypothesize that: (i) the interior desertification from Mongolia to Western Asia began approximately 50 Mya; (ii) the Tibetan plateau (at least on the eastern fringe) experienced rapid uplift approximately 40 Mya and reached an altitude of at least 2,500 m; and (iii) the Ailao-Red River shear zone underwent the most intense orogenic movement approximately 24 Mya.

Structure and evolution of the mitochondrial control region and flanking sequences in the European cave salamander Proteus anguinus

The European cave salamander Proteus anguinus Laurenti 1768 is one of the best-known subterranean animals, yet its evolutionary history and systematic relationships remain enigmatic. This is the first comprehensive study on molecular evolution within the taxon, using an mtDNA segment containing the control region (CR) and adjacent sequences. Two to seven tandem repeats of 24-32 bp were found in the intergenic spacer region (VNTR1), and three, four or six repeats, 59-77 bp each, in the 3' end of the CR (VNTR2). Different molecular mechanisms account for VNTR2 formation in different lineages of Proteus. The overall CR variation was lower than that of the spacer region, the 3' end of the cytb gene, or the tRNA genes. Individual genes and the concatenated non-repetitive sequences produced similar, well resolved maximum likelihood, Bayesian inference and parsimony trees. The numbers of repeat elements as well as the genealogy of the VNTR2 repeat units were mostly inconsistent with the groupings of the non-repetitive sequences. Different degrees of repeat array homogenization were detected in all major groups. Orthology was established for the first and the second VNTR2 elements of some populations. These two copies may therefore be used for analyses at the population level. The pattern of CR sequence variation points to strong genetic isolation of hydrographically separated populations. Genetic separation of the major groups of populations is incongruent with the current division into subspecies.

Evolutionary Rates, Divergence Dates, and the Performance of Mitochondrial Genes in Bayesian Phylogenetic Analysis

The mitochondrial genome is one of the most frequently used loci in phylogenetic and phylogeographic analyses, and it is becoming increasingly possible to sequence and analyze this genome in its entirety from diverse taxa. However, sequencing the entire genome is not always desirable or feasible. Which genes should be selected to best infer the evolutionary history of the mitochondria within a group of organisms, and what properties of a gene determine its phylogenetic performance? The current study addresses these questions in a Bayesian phylogenetic framework with reference to a phylogeny of plethodontid and related salamanders derived from 27 complete mitochondrial genomes; this topology is corroborated by nuclear DNA and morphological data. Evolutionary rates for each mitochondrial gene and divergence dates for all nodes in the plethodontid mitochondrial genome phylogeny were estimated in both Bayesian and maximum likelihood frameworks using multiple fossil calibrations, multiple data partitions, and a clock-independent approach. Bayesian analyses of individual genes were performed, and the resulting trees compared against the reference topology. Ordinal logistic regression analysis of molecular evolution rate, gene length, and the G-shape parameter a demonstrated that slower rate of evolution and longer gene length both increased the probability that a gene would perform well phylogenetically. Estimated rates of molecular evolution vary 84-fold among different mitochondrial genes and different salamander lineages, and mean rates among genes vary 15-fold. Despite having conserved amino acid sequences, cox1, cox2, cox3, and cob have the fastest mean rates of nucleotide substitution, and the greatest variation in rates, whereas rrnS and rrnL have the slowest rates. Reasons underlying this rate variation are discussed, as is the extensive rate variation in cox1 in light of its proposed role in DNA barcoding.