An integrative taxonomic analysis reveals a new species of lotic Hynobius salamander from Japan

We examine the phylogenetic structure and morphological differentiation within the Hynobius kimurae–H. boulengeri species complex—a widely-distributed group of stream-breeding hynobiid salamanders, inhabiting montane areas of western, central and eastern parts of Honshu Island, Japan. Phylogenetic relationships were assessed based on analyses of mitochondrial (mtDNA) and nuclear (nuDNA) gene fragments for a total of 51 samples representing 23 localities covering the entire range of the species complex. Morphological study included one-way analysis of variance (ANOVA) and principal components analysis (PCA) analyses of 26 morphometric and six meristic characters for 38 adult specimens of H. kimurae and three adult specimens of H. boulengeri. MtDNA genealogy supported monophyly of the H. kimurae–H. boulengeri complex, which is sister to all other Hynobius except H. retardatus. The complex is subdivided into three major clades, corresponding to the Eastern (Clade I) and Western (Clade II) populations of H. kimurae sensu lato, and to H. boulengeri (Clade III). Monophyly of H. kimurae sensu lato is only moderately supported by mtDNA, while nuDNA suggested that the Western form of H. kimurae is closer to H. boulengeri than to the eastern form. The time of the split of the H. kimurae–H. boulengeri complex is estimated as late Miocene and coincides with intensive crust movement in western Japan. Divergence between Clades I and II took place in early Pliocene and was likely influenced by the uplift of Central Japanese Highlands. All three clades were found to be different in a number of morphological characters, allowing us to describe the eastern form of H. kimurae as a new species, Hynobius fossigenussp. nov.

Osteology of Batrachuperus londongensis (Urodela, Hynobiidae): study of bony anatomy of a facultatively neotenic salamander from Mount Emei, Sichuan Province, China

The Longdong Stream Salamander Batrachuperus londongensis, living in a mountain stream environment at Mt. Emei in Sichuan Province, China, represents a rare species that is facultatively neotenic in the family Hynobiidae. Although the species has been known to science for some 40 years since its initial discovery in the late 1970s, anatomical details of its osteology remain poorly understood and developmental information is still lacking for the species. This study (1) provides a detailed osteological account of B. londongensis based on micro-CT scanning and clearing and staining of multiple specimens from the type locality; (2) provides a discussion of intraspecific variation related to life-history differences; and (3) presents a discussion on limb features related to morphological evolution of limb patterns correlative with ecological adaptation to mountain stream environments. Osteological comparisons with congeneric species has led to recognition of several diagnostic features that are unique to B. londongensis, including: vomers widely separated from one another, lacking a midline contact; presence of uncommon perichondral ossification of the ascending process of the palatoquadrate as part of the suspensorium; and presence of a prominent posterodorsal process of the scapular blade, which serves as a ligamentous insertion of the levator muscle of the scapula. In addition, some but not all neotenic individuals retain the palatine as a discrete element, indicative of its delayed absorption after sexual maturity. Postmetamorphic and neotenic individuals are strikingly different in the complexity of hyobranchial structures. Neotenes display a high degree of ossification of hyobranchial elements, tend to increase ossification of both hypobranchial I and ceratobranchial I during aging, and retain fully ossified ceratobranchial III and IV; in contrast, these elements remain entirely cartilaginous or are totally lost by resorption in postmetamorphic individuals. In addition, all postmetamorphic forms display an inverted "T"-shaped basibranchial II, whereas neotenes show transformation from a "fork"-shaped to the "T"-shaped configuration after sexual maturity. B. londongensis displays a mosaic of apomorphic and plesiomorphic states in its limb ossifications: presence of a single centrale element in both the manus and pes is a derived condition in Hynobiidae and other families as well, whereas retention of a postminimus in the pes is obviously plesiomorphic within Urodela. Reduction in number of digits from five to four in the pes and possession of a cornified sheath covering the terminal phalanges are also derived features shared with some but not all mountain stream salamanders that are adapted to a similar type of environment.

Explosive diversification of marine fishes at the Cretaceous-Palaeogene boundary

The Cretaceous-Palaeogene (K-Pg) mass extinction is linked to the rapid emergence of ecologically divergent higher taxa (for example, families and orders) across terrestrial vertebrates, but its impact on the diversification of marine vertebrates is less clear. Spiny-rayed fishes (Acanthomorpha) provide an ideal system for exploring the effects of the K-Pg on fish diversification, yet despite decades of morphological and molecular phylogenetic efforts, resolution of both early diverging lineages and enormously diverse subclades remains problematic. Recent multilocus studies have provided the first resolved phylogenetic backbone for acanthomorphs and suggested novel relationships among major lineages. However, these new relationships and associated timescales have not been interrogated using phylogenomic approaches. Here, we use targeted enrichment of >1,000 ultraconserved elements in conjunction with a divergence time analysis to resolve relationships among 120 major acanthomorph lineages and provide a new timescale for acanthomorph radiation. Our results include a well-supported topology that strongly resolves relationships along the acanthomorph backbone and the recovery of several new relationships within six major percomorph subclades. Divergence time analyses also reveal that crown ages for five of these subclades, and for the bulk of the species diversity in the sixth, coincide with the K-Pg boundary, with divergences between anatomically and ecologically distinctive suprafamilial clades concentrated in the first 10 million years of the Cenozoic.

Mitogenome Sequencing in the Genus Camelus Reveals Evidence for Purifying Selection and Long-term Divergence between Wild and Domestic Bactrian Camels

The genus Camelus is an interesting model to study adaptive evolution in the mitochondrial genome, as the three extant Old World camel species inhabit hot and low-altitude as well as cold and high-altitude deserts. We sequenced 24 camel mitogenomes and combined them with three previously published sequences to study the role of natural selection under different environmental pressure, and to advance our understanding of the evolutionary history of the genus Camelus. We confirmed the heterogeneity of divergence across different components of the electron transport system. Lineage-specific analysis of mitochondrial protein evolution revealed a significant effect of purifying selection in the concatenated protein-coding genes in domestic Bactrian camels. The estimated dN/dS < 1 in the concatenated protein-coding genes suggested purifying selection as driving force for shaping mitogenome diversity in camels. Additional analyses of the functional divergence in amino acid changes between species-specific lineages indicated fixed substitutions in various genes, with radical effects on the physicochemical properties of the protein products. The evolutionary time estimates revealed a divergence between domestic and wild Bactrian camels around 1.1 [0.58-1.8] million years ago (mya). This has major implications for the conservation and management of the critically endangered wild species, Camelus ferus.

Evolutionary history of the extant amphioxus lineage with shallow-branching diversification

Amphioxus or lancelets have been regarded as a key animal in understanding the origin of vertebrates. However, the evolutionary history within this lineage remains unexplored. As the amphioxus lineage has likely been separated from other chordates for a very long time and displays a marked left-right asymmetry, its evolutionary history is potentially helpful in better understanding chordate and vertebrate origins. We studied the phylogenetic relationships within the extant amphioxus lineage based on mitochondrial genomes incorporating new Asymmetron and Epigonichthys populations, and based on previously reported nuclear transcriptomes. The resulting tree patterns are consistent, showing the Asymmetron clade diverging first, followed by the Epigonichthys and Branchiostoma clades splitting. Divergence time estimates based on nuclear transcriptomes with vertebrate calibrations support a shallow diversification of the extant amphioxus lineage in the Tertiary. These estimates fit well with the closure of seaways between oceans by continental drift, ocean currents, and present geographical distributions, and suggest a long cryptic history from the origin of amphioxus to its most recent diversification. Deduced character polarities based on phylogenetic analyses suggest that the common ancestor of the extant amphioxus existed in a tiny epibenthic state with larva-like appearance of extant amphioxus, likely with ciliate epidermis.

Colonization and diversification of aquatic insects on three Macaronesian archipelagos using 59 nuclear loci derived from a draft genome

The study of processes driving diversification requires a fully sampled and well resolved phylogeny, although a lack of phylogenetic markers remains a limitation for many non-model groups. Multilocus approaches to the study of recent diversification provide a powerful means to study the evolutionary process, but their application remains restricted because multiple unlinked loci with suitable variation for phylogenetic or coalescent analysis are not available for most non-model taxa. Here we identify novel, putative single-copy nuclear DNA (nDNA) phylogenetic markers to study the colonization and diversification of an aquatic insect species complex, Cloeon dipterum L. 1761 (Ephemeroptera: Baetidae), in Macaronesia. Whole-genome sequencing data from one member of the species complex were used to identify 59 nDNA loci (32,213 base pairs), followed by Sanger sequencing of 29 individuals sampled from 13 islands of three Macaronesian archipelagos. Multispecies coalescent analyses established six putative species. Three island species formed a monophyletic clade, with one species occurring on the Azores, Europe and North America. Ancestral state reconstruction indicated at least two colonization events from the mainland (to the Canaries, respectively Azores) and one within the archipelago (between Madeira and the Canaries). Random subsets of the 59 loci showed a positive linear relationship between number of loci and node support. In contrast, node support in the multispecies coalescent tree was negatively correlated with mean number of phylogenetically informative sites per locus, suggesting a complex relationship between tree resolution and marker variability. Our approach highlights the value of combining genomics, coalescent-based phylogeography, species delimitation, and phylogenetic reconstruction to resolve recent diversification events in an archipelago species complex.

What explains patterns of species richness? The relative importance of climatic-niche evolution, morphological evolution, and ecological limits in salamanders

A major goal of evolutionary biology and ecology is to understand why species richness varies among clades. Previous studies have suggested that variation in richness among clades might be related to variation in rates of morphological evolution among clades (e.g., body size and shape). Other studies have suggested that richness patterns might be related to variation in rates of climatic-niche evolution. However, few studies, if any, have tested the relative importance of these variables in explaining patterns of richness among clades. Here, we test their relative importance among major clades of Plethodontidae, the most species-rich family of salamanders. Earlier studies have suggested that climatic-niche evolution explains patterns of diversification among plethodontid clades, whereas rates of morphological evolution do not. A subsequent study stated that rates of morphological evolution instead explained patterns of species richness among plethodontid clades (along with "ecological limits" on richness of clades, leading to saturation of clades with species, given limited resources). However, they did not consider climatic-niche evolution. Using phylogenetic multiple regression, we show that rates of climatic-niche evolution explain most variation in richness among plethodontid clades, whereas rates of morphological evolution do not. We find little evidence that ecological limits explain patterns of richness among plethodontid clades. We also test whether rates of morphological and climatic-niche evolution are correlated, and find that they are not. Overall, our results help explain richness patterns in a major amphibian group and provide possibly the first test of the relative importance of climatic niches and morphological evolution in explaining diversity patterns.

Population Structure and Evolution after Speciation of the Hokkaido Salamander (Hynobius retardatus)

The Hokkaido salamander (Hynobius retardatus) is endemic to Hokkaido Island, Japan, and shows intriguing flexible phenotypic plasticity and regional morphological diversity. However, to date, allozymes and partial mitochondria DNA sequences have provided only an outline of its demographic histories and the pattern of its genetic diversification. To understand the finer details of the population structure of this species and its evolution since speciation, we genotyped five regional populations by using 12 recently developed microsatellite polymorphic markers. We found a clear population structure with low gene flow among the five populations, but a close genetic relationship between the Teshio and Kitami populations. Our demographic analysis suggested that Teshio and Erimo had the largest effective population sizes among the five populations. These findings regarding the population structure and demography of H. retardatus improve our understanding of the faunal phylogeography on Hokkaido Island and also provide fundamental genetic information that will be useful for future studies.

Molecular and morphological phylogenetics of chelonine parasitoid wasps (Hymenoptera: Braconidae), with a critical assessment of divergence time estimations

Parasitoid wasps of the subfamily Cheloninae are both species rich and poorly known. Although the taxonomy of Cheloninae appears to be relatively stable, there is no clear understanding of relationships among higher-level taxa. We here applied molecular phylogenetic analyses using three markers (COI, EF1α, 28S) and 37 morphological characters to elucidate the evolution and systematics of these wasps. Analyses were based on 83 specimens representing 13 genera. All genera except Ascogaster, Phanerotoma, and Pseudophanerotoma formed monophyletic groups; Furcidentia (stat. rev.) is raised to generic rank. Neither Chelonus (Chelonus) nor Chelonus (Microchelonus) were recovered as monophyletic, but together formed a monophyletic lineage. The tribes Chelonini and Odontosphaeropygini formed monophyletic groups, but the Phanerotomini sensu Zettel and Pseudophanerotomini were retrieved as either para- or polyphyletic. The genera comprising the former subfamily Adeliinae were confirmed as being nested within the Cheloninae. To estimate the age of the subfamily, we used 16 fossil taxa. Three approaches were compared: fixed-rate dating, node dating, and total-evidence dating, with age estimates differing greatly between the three methods. Shortcomings of each approach in relation to our dataset are discussed, and none of the age estimates is deemed sufficiently reliable. Given that most dating studies use a single method only, in most cases without presenting analyses on the sensitivity to priors, it is likely that numerous age estimates in the literature suffer from a similar lack of robustness. We argue for a more rigorous approach to dating analyses and for a faithful presentation of uncertainties in divergence time estimates. Given the results of the phylogenetic analysis the following taxonomic changes are proposed: Furcidentia Zettel (stat. rev.), previously treated as a subgenus of Pseudophanerotoma Zettel is raised to generic rank; Microchelonus Szépligeti (syn. nov.), variously treated by previous authors, is proposed as a junior synonym of Chelonus Jurine; the following subgenera of Microchelonus - Baculonus Braet & van Achterberg (syn. nov.), Carinichelonus Tobias (syn. nov.) and Scabrichelonus He, Chen & van Achterberg (syn. nov.), are proposed as junior synonyms of Chelonus; a number of new species names are proposed due to homonyms resulting from the above changes and these are listed in the paper.

Consequences of Secondary Calibrations on Divergence Time Estimates

Secondary calibrations (calibrations based on the results of previous molecular dating studies) are commonly applied in divergence time analyses in groups that lack fossil data; however, the consequences of applying secondary calibrations in a relaxed-clock approach are not fully understood. I tested whether applying the posterior estimate from a primary study as a prior distribution in a secondary study results in consistent age and uncertainty estimates. I compared age estimates from simulations with 100 randomly replicated secondary trees. On average, the 95% credible intervals of node ages for secondary estimates were significantly younger and narrower than primary estimates. The primary and secondary age estimates were significantly different in 97% of the replicates after Bonferroni corrections. Greater error in magnitude was associated with deeper than shallower nodes, but the opposite was found when standardized by median node age, and a significant positive relationship was determined between the number of tips/age of secondary trees and the total amount of error. When two secondary calibrated nodes were analyzed, estimates remained significantly different, and although the minimum and median estimates were associated with less error, maximum age estimates and credible interval widths had greater error. The shape of the prior also influenced error, in which applying a normal, rather than uniform, prior distribution resulted in greater error. Secondary calibrations, in summary, lead to a false impression of precision and the distribution of age estimates shift away from those that would be inferred by the primary analysis. These results suggest that secondary calibrations should not be applied as the only source of calibration in divergence time analyses that test time-dependent hypotheses until the additional error associated with secondary calibrations is more properly modeled to take into account increased uncertainty in age estimates.

Stop and Go - Waves of Tarsier Dispersal Mirror the Genesis of Sulawesi Island

The Indonesian island of Sulawesi harbors a highly endemic and diverse fauna sparking fascination since long before Wallace’s contemplation of biogeographical patterns in the region. Allopatric diversification driven by geological or climatic processes has been identified as the main mechanism shaping present faunal distribution on the island. There is both consensus and conflict among range patterns of terrestrial species pointing to the different effects of vicariant events on once co-distributed taxa. Tarsiers, small nocturnal primates with possible evidence of an Eocene fossil record on the Asian mainland, are at present exclusively found in insular Southeast Asia. Sulawesi is hotspot of tarsier diversity, whereby island colonization and subsequent radiation of this old endemic primate lineage remained largely enigmatic. To resolve the phylogeographic history of Sulawesi tarsiers we analyzed an island-wide sample for a set of five approved autosomal phylogenetic markers (ABCA1, ADORA3, AXIN1, RAG1, and TTR) and the paternally inherited SRY gene. We constructed ML and Bayesian phylogenetic trees and estimated divergence times between tarsier populations. We found that their arrival at the Proto-Sulawesi archipelago coincided with initial Miocene tectonic uplift and hypothesize that tarsiers dispersed over the region in distinct waves. Intra-island diversification was spurred by land emergence and a rapid succession of glacial cycles during the Plio-Pleistocene. Some tarsier range boundaries concur with spatial limits in other taxa backing the notion of centers of faunal endemism on Sulawesi. This congruence, however, has partially been superimposed by taxon-specific dispersal patterns.

Phylogeny and Evolution of Multiple Common Carp (Cyprinus carpio L.) Populations Clarified by Phylogenetic Analysis Based on Complete Mitochondrial Genomes

Common carp (Cyprinus carpio L.) is one of the oldest, most widely farmed commercially important freshwater fish in the world. However, many undetermined phylogenetic relationships and origins of common carp lineages remain, which are obstacles to conservation and genetic breeding of this species. Phylogenetic analyses based on molecular tools are helpful to distinguish the origin of species, understand and clarify their evolutionary history, and provide a genetic basis for selective breeding. In this study, we demonstrated a method to extract complete mitochondrial genome sequences from whole-genome resequencing data using the Illumina platform. The complete mitochondrial genome sequences of 26 individuals representing nine strains were obtained and subjected to a phylogenetic analysis. We reconstructed the phylogenetic topologies of the nine strains and analyzed the haplotypes. Results from both analyses suggested that the genome sequences belonged to two distinct subspecies from Europe and East Asia. We also estimated the time of divergence of the nine strains, which was up to 100 KYA. The phylogenetic results clarified the breeding history of Songpu mirror carp and suggest that this species may be hybrid of paternal European mirror carp and maternal Xingguo red carp. The results also support a previous hypothesis that koi may have originated from or have close ancestry with Oujiang color carp in China.

Enlarged Multilocus Data set Provides Surprisingly Younger Time of Origin for the Plethodontidae, the Largest Family of Salamanders

Deep phylogenetic relationships of the largest salamander family Plethodontidae have been difficult to resolve, probably reflecting a rapid diversification early in their evolutionary history. Here, data from 50 independent nuclear markers (total 48,582 bp) are used to reconstruct the phylogeny and divergence times for plethodontid salamanders, using both concatenation and coalescence-based species tree analyses. Our results robustly resolve the position of the enigmatic eastern North American four-toed salamander (Hemidactylium) as the sister taxon of Batrachoseps + Tribe Bolitoglossini, thus settling a long-standing question. Furthermore, we statistically reject sister taxon status of Karsenia and Hydromantes, the only plethodontids to occur outside the Americas, leading us to new biogeographic hypotheses. Contrary to previous long-standing arguments that plethodontid salamanders are an old lineage originating in the Cretaceous (more than 90 Ma), our analyses lead to the hypothesis that these salamanders are much younger, arising close to the K-T boundary (~66 Ma). These time estimates are highly stable using alternative calibration schemes and dating methods. Our data simulation highlights the potential risk of making strong arguments about phylogenetic timing based on inferences from a handful of nuclear genes, a common practice. Based on the newly obtained timetree and ancestral area reconstruction results, we argue that (i) the classic "Out of Appalachia" hypothesis of plethodontid origins is problematic; (ii) the common ancestor of extant plethodontids may have originated in northwestern North America in the early Paleocene; (iii) origins of Eurasian plethodontids likely result from two separate dispersal events from western North America via Beringia in the late Eocene (~42 Ma) and the early Miocene (~23 Ma), respectively.

Is geographic variation within species related to macroevolutionary patterns between species?

The relationship between microevolution and macroevolution is a central topic in evolutionary biology. An aspect of this relationship that remains very poorly studied in modern evolutionary biology is the relationship between within-species geographic variation and among-species patterns of trait variation. Here, we tested the relationship between climate and morphology among and within species in the salamander genus Plethodon. We focus on a discrete colour polymorphism (presence and absence of a red dorsal stripe) that appears to be related to climatic distributions in a common, wide-ranging species (Plethodon cinereus). We find that this trait has been variable among (and possibly within) species for >40 million years. Furthermore, we find a strong relationship among species between climatic variation and within-species morph frequencies. These between-species patterns are similar (but not identical) to those in the broadly distributed Plethodon cinereus. Surprisingly, there are no significant climate-morphology relationships within most other polymorphic species, despite the strong between-species patterns. Overall, our study provides an initial exploration of how within-species geographic variation and large-scale macroevolutionary patterns of trait variation may be related.

Potential for bias and low precision in molecular divergence time estimation of the Canopy of Life: an example from aquatic bird families

Uncertainty in divergence time estimation is frequently studied from many angles but rarely from the perspective of phylogenetic node age. If appropriate molecular models and fossil priors are used, a multi-locus, partitioned analysis is expected to equally minimize error in accuracy and precision across all nodes of a given phylogeny. In contrast, if available models fail to completely account for rate heterogeneity, substitution saturation and incompleteness of the fossil record, uncertainty in divergence time estimation may increase with node age. While many studies have stressed this concern with regard to deep nodes in the Tree of Life, the inference that molecular divergence time estimation of shallow nodes is less sensitive to erroneous model choice has not been tested explicitly in a Bayesian framework. Because of available divergence time estimation methods that permit fossil priors across any phylogenetic node and the present increase in efficient, cheap collection of species-level genomic data, insight is needed into the performance of divergence time estimation of shallow (<10 MY) nodes. Here, we performed multiple sensitivity analyses in a multi-locus data set of aquatic birds with six fossil constraints. Comparison across divergence time analyses that varied taxon and locus sampling, number and position of fossil constraint and shape of prior distribution showed various insights. Deviation from node ages obtained from a reference analysis was generally highest for the shallowest nodes but determined more by temporal placement than number of fossil constraints. Calibration with only the shallowest nodes significantly underestimated the aquatic bird fossil record, indicating the presence of saturation. Although joint calibration with all six priors yielded ages most consistent with the fossil record, ages of shallow nodes were overestimated. This bias was found in both mtDNA and nDNA regions. Thus, divergence time estimation of shallow nodes may suffer from bias and low precision, even when appropriate fossil priors and best available substitution models are chosen. Much care must be taken to address the possible ramifications of substitution saturation across the entire Tree of Life.

A Diverse Tetrapod Fauna at the Base of 'Romer's Gap'

The lack of fossil tetrapod bearing deposits in the earliest Carboniferous (‘Romer’s Gap’) has provoked some recent discussions regarding the proximal cause, with three explanations being offered: environmental, taphonomic, and collection failure. One of the few, and earliest, windows into this time is the locality of Blue Beach exposed in the Tournaisian deposits at Horton Bluff lying along the Avon River near Hantsport, Nova Scotia, Canada. This locality has long been known but, because the fossils were deposited in high energy settings they are almost always disarticulated, so the fauna has not been described in detail. Recent intensive collection has revealed a diverse assemblage of material, including for the first time associated elements, which permits an evaluation of the faunal constituents at the locality. Although not diagnosable to a fine taxonomic level, sufficient apomorphies are present to identify representatives from numerous clades known from more complete specimens elsewhere. The evidence suggests a diverse fauna was present, including whatcheeriids and embolomeres. A single humerus previously had been attributed to a colosteid, but there is some uncertainty with this identification. Additional elements suggest the presence of taxa otherwise only known from the late Devonian. Depositional biases at the locality favor tetrapod fossils from larger individuals, but indirect evidence from trackways and tantalizing isolated bones evidences the presence of small taxa that remain to be discovered. The fossils from Blue Beach demonstrate that when windows into the fauna of ‘Romer’s Gap’ are found a rich diversity of tetrapods will be shown to be present, contra arguments that suggested this hiatus in the fossil record was due to extrinsic factors such as atmospheric oxygen levels. They also show that the early tetrapod fauna is not easily divisible into Devonian and Carboniferous faunas, suggesting that some tetrapods passed through the end Devonian extinction event unaffected.

Hynobiidae origin in middle Cretaceous corroborated by the new mitochondrial genome of Hynobius chinensis

Hynobius chinensis was first described by Günther in the nineteenth century. At present, the origins of the extinct Hynobius chinensis on the Zhoushan Island (Hynobius chinensis-ZI) remain a mystery. It is the only species of family Hynobiidae on the Zhoushan Island. However, there is very little empirical evidence regarding Hynobius chinensis-ZI phylogenetic relationship, and when or how did its ancestors colonized the island. Here, we used mitochondrial genome data to recover the phylogeny of family Hynobiidae. Results suggested that the origin of Hynobiidae was most likely in Middle Cretaceous (~112.9 Mya), and some Hynobius species of Taiwan and Japan diverged earlier than that of the mainland of China. Hynobius chinensis-ZI diverged from its closest living relative (Hynobius yiwuensis) around 6.5 Mya, and Hynobius chinensis-ZI was isolated on Zhoushan Island since the postglacial transgression in Holocene period.

Low levels of LTR retrotransposon deletion by ectopic recombination in the gigantic genomes of salamanders

Across the tree of life, species vary dramatically in nuclear genome size. Mutations that add or remove sequences from genomes-insertions or deletions, or indels-are the ultimate source of this variation. Differences in the tempo and mode of insertion and deletion across taxa have been proposed to contribute to evolutionary diversity in genome size. Among vertebrates, most of the largest genomes are found within the salamanders, an amphibian clade with genome sizes ranging from ~14 to ~120 Gb. Salamander genomes have been shown to experience slower rates of DNA loss through small (i.e., <30 bp) deletions than do other vertebrate genomes. However, no studies have addressed DNA loss from salamander genomes resulting from larger deletions. Here, we focus on one type of large deletion-ectopic-recombination-mediated removal of LTR retrotransposon sequences. In ectopic recombination, double-strand breaks are repaired using a "wrong" (i.e., ectopic, or non-allelic) template sequence-typically another locus of similar sequence. When breaks occur within the LTR portions of LTR retrotransposons, ectopic-recombination-mediated repair can produce deletions that remove the internal transposon sequence and the equivalent of one of the two LTR sequences. These deletions leave a signature in the genome-a solo LTR sequence. We compared levels of solo LTRs in the genomes of four salamander species with levels present in five vertebrates with smaller genomes. Our results demonstrate that salamanders have low levels of solo LTRs, suggesting that ectopic-recombination-mediated deletion of LTR retrotransposons occurs more slowly than in other vertebrates with smaller genomes.

Phylogenetic relationships of Amazonian anole lizards (Dactyloa): taxonomic implications, new insights about phenotypic evolution and the timing of diversification

The ecology and evolution of Caribbean anoles are well described, yet little is known about mainland anole species. Lack of phylogenetic information limits our knowledge about species boundaries, morphological evolution, and the biogeography of anoles in South America. To help fill this gap, we provide an updated molecular phylogeny of the Dactyloa (Dactyloidae), with emphasis on the punctata species group. By sampling understudied Amazonian taxa, we (i) assess the phylogenetic placement of the 'odd anole', D. dissimilis; (ii) infer the relationships of the proboscis-bearing D. phyllorhina, testing the hypothesis of independent nasal appendage evolution within the anole radiation; and (iii) examine genetic and dewlap color variation in D. punctata and D. philopunctata. Combining multiple nuclear loci with a review of the fossil record, we also (iv) estimate divergence times within the pleurodont iguanian clade of lizards, including Amazonian representatives of Dactyloa and Norops (Dactyloidae) and of Polychrus (Polychrotidae). We recover the five Dactyloa clades previously referred to as the aequatorialis, heteroderma, latifrons, punctata and roquet species groups, as well as a sixth clade composed of D. dissimilis and the non-Amazonian D. neblinina and D. calimae. We find D. phyllorhina to be nested within the punctata group, suggesting independent evolution of the anole proboscis. We consistently recover D. philopunctata nested within D. punctata, and report limited genetic divergence between distinct dewlap phenotypes. The most recent common ancestor of Dactyloa, Anolis and Norops dates back to the Eocene. Most Amazonian taxa within both Dactyloa and Norops diverged in the Miocene, but some diversification events were as old as the late Eocene and late Oligocene. Amazonian Polychrus diverged in the Pliocene. Our findings have broad implications for anole biogeography, disputing recent suggestions that modern dactyloid genera were present in the Caribbean region during the Cretaceous.

Comparative and phylogenetic perspectives of the cleavage process in tailed amphibians

The order Caudata includes about 660 species and displays a variety of important developmental traits such as cleavage pattern and egg size. However, the cleavage process of tailed amphibians has never been analyzed within a phylogenetic framework. We use published data on the embryos of 36 species concerning the character of the third cleavage furrow (latitudinal, longitudinal or variable) and the magnitude of synchronous cleavage period (up to 3–4 synchronous cell divisions in the animal hemisphere or a considerably longer series of synchronous divisions followed by midblastula transition). Several species from basal caudate families Cryptobranchidae (Andrias davidianus and Cryptobranchus alleganiensis) and Hynobiidae (Onychodactylus japonicus) as well as several representatives from derived families Plethodontidae (Desmognathus fuscus and Ensatina eschscholtzii) and Proteidae (Necturus maculosus) are characterized by longitudinal furrows of the third cleavage and the loss of synchrony as early as the 8-cell stage. By contrast, many representatives of derived families Ambystomatidae and Salamandridae have latitudinal furrows of the third cleavage and extensive period of synchronous divisions. Our analysis of these ontogenetic characters mapped onto a phylogenetic tree shows that the cleavage pattern of large, yolky eggs with short series of synchronous divisions is an ancestral trait for the tailed amphibians, while the data on the orientation of third cleavage furrows seem to be ambiguous with respect to phylogeny. Nevertheless, the midblastula transition, which is characteristic of the model species Ambystoma mexicanum (Caudata) and Xenopus laevis (Anura), might have evolved convergently in these two amphibian orders.

Phylogenetic informativeness reconciles ray-finned fish molecular divergence times

Background

Discordance among individual molecular age estimates, or between molecular age estimates and the fossil record, is observed in many clades across the Tree of Life. This discordance is attributed to a variety of variables including calibration age uncertainty, calibration placement, nucleotide substitution rate heterogeneity, or the specified molecular clock model. However, the impact of changes in phylogenetic informativeness of individual genes over time on phylogenetic inferences is rarely analyzed. Using nuclear and mitochondrial sequence data for ray-finned fishes (Actinopterygii) as an example, we extend the utility of phylogenetic informativeness profiles to predict the time intervals when nucleotide substitution saturation results in discordance among molecular ages estimated.

Results

We demonstrate that even with identical calibration regimes and molecular clock methods, mitochondrial based molecular age estimates are systematically older than those estimated from nuclear sequences. This discordance is most severe for highly nested nodes corresponding to more recent (i.e., Jurassic-Recent) divergences. By removing data deemed saturated, we reconcile the competing age estimates and highlight that the older mtDNA based ages were driven by nucleotide saturation.

Conclusions

Homoplasious site patterns in a DNA sequence alignment can systematically bias molecular divergence time estimates. Our study demonstrates that PI profiles can provide a non-arbitrary criterion for data exclusion to mitigate the influence of homoplasy on time calibrated branch length estimates. Analyses of actinopterygian molecular clocks demonstrate that scrutiny of the time scale on which sequence data is informative is a fundamental, but generally overlooked, step in molecular divergence time estimation.

Mitogenomic sequences effectively recover relationships within brush-footed butterflies (Lepidoptera: Nymphalidae)

BACKGROUND

Mitogenomic phylogenies have revealed well-supported relationships for many eukaryote groups. In the order Lepidoptera, 113 species mitogenomes had been sequenced (May 14, 2014). However, these data are restricted to ten of the forty-three recognised superfamilies, while it has been challenging to recover large numbers of mitogenomes due to the time and cost required for primer design and sequencing. Nuclear rather than mitochondrial genes have been preferred to reconstruct deep-level lepidopteran phylogenies, without seriously evaluating the potential of entire mitogenomes. Next-generation sequencing methods remove these limitations by providing efficiently massive amounts of sequence data. In the present study, we simultaneously obtained a large number of nymphalid butterfly mitogenomes to evaluate the utility of mitogenomic phylogenies by comparing reconstructions to the now quite well established phylogeny of Nymphalidae.

RESULTS

We newly obtained 30 nymphalid mitogenomes via pyrosequencing on the Roche 454 GS Junior system, and combined these sequences with publicly accessible data to provide a 70-taxa dataset covering 37 genes for a 15,495 bp alignment. Polymorphic sites were not homogeneously distributed across the gene. Two gene regions, nad6 and 3' end of nad5, were most variable, whereas the cox1 and 5' ends of rrnL were most conserved. Phylogenetic relationships inferred by two likelihood methods were congruent and strongly supported (>0.95 posterior probability; ML bootstrap >85%), across the majority of nodes for multiple partitioning strategies and substitution models. Bayes factor results showed that the most highly partitioned dataset is the preferred strategy among different partitioning schemes. The most striking phylogenetic findings were that the subfamily Danainae not Libytheinae was sister of the remaining brush-footed butterflies and that, within Limenitidini, the genus Athyma was clearly polyphyletic. None of the single-gene phylogenies recovered the highly supported topologies generated on the basis of the whole mitogenomic data.

CONCLUSIONS

Thirty mitogenomes were assembled with 89% completeness from the contigs of pyrosequencing-derived reads. Entire mitogenomes or higher-quality sequences could be obtained by increasing pyrosequencing read coverage or by additional Sanger sequencing. Our mitogenomic phylogenies provide robust nodal support at a range of levels, demonstrating that mitogenomes are both accurate and efficient molecular markers for inferring butterfly phylogeny.

Molecular phylogeny and biogeography of air sac catfishes of the Heteropneustes fossilis species complex (Siluriformes: Heteropneustidae)

The air sac catfish, Heteropneustes fossilis (Siluriformes: Heteropneustidae), is widely distributed in freshwaters of the Indian subcontinent and mainland southeast Asia. No comprehensive molecular studies that cover the broad distributional areas have been carried out to date. Here, we conducted molecular phylogenetic analyses using both mitochondrial and nuclear gene sequences to suggest that the Heteropneustes fossilis species complex consists of three clades that may potentially be separate species with distinct geographical distribution (southeast Asia, northeastern India, and southwestern India). The first and second clades are more closely related to each other than they are to the third clade. Within the first clade there is a basal divergence of a subclade consisting of individuals from the Upper Irrawaddy River basin of Myanmar, which share some morphological traits with members of the Indian clades. Our molecular and morphological data are congruent with hypotheses that the Early-Middle Miocene disconnection between the paleo-Tsangpo River and the Irrawaddy River caused the vicariant divergence between southeast Asian and northeastern Indian clades, and that the southeast Asian Heteropneustes originated from the Upper Irrawaddy.

The complete mitochondrial genome of the Hynobius maoershanensis (Caudata: hynobiidae)

The Hynobius maoershanensis is a member of hynobiidae, endemic to Mountain Maoer in Guangxi province, China. It was first found and reported in 2006 and so far there is a little molecular research about it. The complete mitochondrial genome of H. maoershanensis has been obtained for the first time in this study. The circular genome (16,412 bp in length) consisted of 37 typical animal mitochondrial genes (13 protein-coding genes, 22 transfer RNA genes, 2 ribosomal RNA genes) and 1 control region. Overall base composition of the complete mitochondrial DNA was 33% A, 32% T, 21% C, and 14% G with AT (65%).

Genetic variation in the green anole lizard (Anolis carolinensis) reveals island refugia and a fragmented Florida during the quaternary

The green anole lizard (Anolis carolinensis) is a model organism for behavior and genomics that is native to the southeastern United States. It is currently thought that the ancestors of modern green anoles dispersed to peninsular Florida from Cuba. However, the climatic changes and geological features responsible for the early diversification of A. carolinensis in North America have remained largely unexplored. This is because previous studies (1) differ in their estimates of the divergence times of populations, (2) are based on a single genetic locus or (3) did not test specific hypotheses regarding the geologic and topographic history of Florida. Here we provide a multi-locus study of green anole genetic diversity and find that the Florida peninsula contains a larger number of genetically distinct populations that are more diverse than those on the continental mainland. As a test of the island refugia hypothesis in Pleistocene Florida, we use a coalescent approach to estimate the divergence times of modern green anole lineages. We find that all demographic events occurred during or after the Upper Pliocene and suggest that green anole diversification was driven by population divergence on interglacial island refugia in Florida during the Lower Pleistocene, while the region was often separated from continental North America. When Florida reconnected to the mainland, two separate dispersal events led to the expansion of green anole populations across the Atlantic Seaboard and Gulf Coastal Plain.

Molecular and fossil evidence place the origin of cichlid fishes long after Gondwanan rifting

Cichlid fishes are a key model system in the study of adaptive radiation, speciation and evolutionary developmental biology. More than 1600 cichlid species inhabit freshwater and marginal marine environments across several southern landmasses. This distributional pattern, combined with parallels between cichlid phylogeny and sequences of Mesozoic continental rifting, has led to the widely accepted hypothesis that cichlids are an ancient group whose major biogeographic patterns arose from Gondwanan vicariance. Although the Early Cretaceous (ca 135 Ma) divergence of living cichlids demanded by the vicariance model now represents a key calibration for teleost molecular clocks, this putative split pre-dates the oldest cichlid fossils by nearly 90 Myr. Here, we provide independent palaeontological and relaxed-molecular-clock estimates for the time of cichlid origin that collectively reject the antiquity of the group required by the Gondwanan vicariance scenario. The distribution of cichlid fossil horizons, the age of stratigraphically consistent outgroup lineages to cichlids and relaxed-clock analysis of a DNA sequence dataset consisting of 10 nuclear genes all deliver overlapping estimates for crown cichlid origin centred on the Palaeocene (ca 65-57 Ma), substantially post-dating the tectonic fragmentation of Gondwana. Our results provide a revised macroevolutionary time scale for cichlids, imply a role for dispersal in generating the observed geographical distribution of this important model clade and add to a growing debate that questions the dominance of the vicariance paradigm of historical biogeography.

Multilocus phylogeny and cryptic diversity in Asian shrew-like moles (Uropsilus, Talpidae): implications for taxonomy and conservation

BACKGROUND

The genus Uropsilus comprises a group of terrestrial, montane mammals endemic to the Hengduan and adjacent mountains. These animals are the most primitive living talpids. The taxonomy has been primarily based on cursory morphological comparisons and the evolutionary affinities are little known. To provide insight into the systematics of this group, we estimated the first multi-locus phylogeny and conducted species delimitation, including taxon sampling throughout their distribution range.

RESULTS

We obtained two mitochondrial genes (~1, 985 bp) and eight nuclear genes (~4, 345 bp) from 56 specimens. Ten distinct evolutionary lineages were recovered from the three recognized species, eight of which were recognized as species/putative species. Five of these putative species were found to be masquerading as the gracile shrew mole. The divergence time estimation results indicated that climate change since the last Miocene and the uplift of the Himalayas may have resulted in the diversification and speciation of Uropsilus.

CONCLUSIONS

The cryptic diversity found in this study indicated that the number of species is strongly underestimated under the current taxonomy. Two synonyms of gracilis (atronates and nivatus) should be given full species status, and the taxonomic status of another three potential species should be evaluated using extensive taxon sampling, comprehensive morphological, and morphometric approaches. Consequently, the conservation status of Uropsilus spp. should also be re-evaluated, as most of the species/potential species have very limited distribution.

Efficient sequencing of Anuran mtDNAs and a mitogenomic exploration of the phylogeny and evolution of frogs

Anura (frogs and toads) constitute over 88% of living amphibian diversity but many important questions about their phylogeny and evolution remain unresolved. For this study, we developed an efficient method for sequencing anuran mitochondrial DNAs (mtDNAs) by amplifying the mitochondrial genome in 12 overlapping fragments using frog-specific universal primer sets. Based on this method, we generated 47 nearly complete, new anuran mitochondrial genomes and discovered nine novel gene arrangements. By combining the new data and published anuran mitochondrial genomes, we assembled a large mitogenomic data set (11,007 nt) including 90 frog species, representing 39 of 53 recognized anuran families, to investigate their phylogenetic relationships and evolutionary history. The resulting tree strongly supported a paraphyletic arrangement of archaeobatrachian (=nonneobatrachian) frogs, with Leiopelmatoidea branching first, followed by Discoglossoidea, Pipoidea, and Pelobatoidea. Within Neobatrachia, the South African Heleophrynidae is the sister-taxon to all other neobatrachian frogs and the Seychelles-endemic Sooglossidae is recovered as the sister-taxon to Ranoidea. These phylogenetic relationships agree with many nuclear gene studies. The chronogram derived from two Bayesian relaxed clock methods (MultiDivTime and BEAST) suggests that modern frogs (Anura) originated in the early Triassic about 244 Ma and the appearance of Neobatrachia took place in the late Jurassic about 163 Ma. The initial diversifications of two species-rich superfamilies Hyloidea and Ranoidea commenced 110 and 133 Ma, respectively. These times are older than some other estimates by approximately 30-40 My. Compared with nuclear data, mtDNA produces compatible time estimates for deep nodes (>150 Ma), but apparently older estimates for more shallow nodes. Our study shows that, although it evolves relatively rapidly and behaves much as a single locus, mtDNA performs well for both phylogenetic and divergence time inferences and will provide important reference hypotheses for the phylogeny and evolution of frogs.

Estimating divergence dates and evaluating dating methods using phylogenomic and mitochondrial data in squamate reptiles

Recently, phylogenetics has expanded to routinely include estimation of clade ages in addition to their relationships. Various dating methods have been used, but their relative performance remains understudied. Here, we generate and assemble an extensive phylogenomic data set for squamate reptiles (lizards and snakes) and evaluate two widely used dating methods, penalized likelihood in r8s (r8s-PL) and Bayesian estimation with uncorrelated relaxed rates among lineages (BEAST). We obtained sequence data from 25 nuclear loci (∼500-1,000 bp per gene; 19,020bp total) for 64 squamate species and nine outgroup taxa, estimated the phylogeny, and estimated divergence dates using 14 fossil calibrations. We then evaluated how well each method approximated these dates using random subsets of the nuclear loci (2, 5, 10, 15, and 20; replicated 10 times each), and using ∼1 kb of the mitochondrial ND2 gene. We find that estimates from r8s-PL based on 2, 5, or 10 loci can differ considerably from those based on 25 loci (mean absolute value of differences between 2-locus and 25-locus estimates were 9.0 Myr). Estimates from BEAST are somewhat more consistent given limited sampling of loci (mean absolute value of differences between 2 and 25-locus estimates were 5.0 Myr). Most strikingly, age estimates using r8s-PL for ND2 were ∼68-82 Myr older (mean=73.1) than those using 25 nuclear loci with r8s-PL. These results show that dates from r8s-PL with a limited number of loci (and especially mitochondrial data) can differ considerably from estimates derived from a large number of nuclear loci, whereas estimates from BEAST derived from fewer nuclear loci or mitochondrial data alone can be surprisingly similar to those from many nuclear loci. However, estimates from BEAST using relatively few loci and mitochondrial data could still show substantial deviations from the full data set (>50 Myr), suggesting the benefits of sampling many nuclear loci. Finally, we found that confidence intervals on ages from BEAST were not significantly different when sampling 2 vs. 25 loci, suggesting that adding loci decreased errors but did not increase confidence in those estimates.

How is the rate of climatic-niche evolution related to climatic-niche breadth?

The rate of climatic-niche evolution is important to many research areas in ecology, evolution, and conservation biology, including responses of species to global climate change, spread of invasive species, speciation, biogeography, and patterns of species richness. Previous studies have implied that clades with higher rates of climatic-niche evolution among species should have species with narrower niche breadths, but there is also evidence suggesting the opposite pattern. However, the relationships between rate and breadth have not been explicitly analyzed. Here, we examine the relationships between the rate of climatic-niche evolution and climatic-niche breadth using phylogenetic and climatic data for 250 species in the salamander family Plethodontidae, a group showing considerable variation in both rates of climatic-niche evolution and climatic-niche breadths. Contrary to some expectations, we find no general relationship between climatic-niche breadth and the rate of climatic-niche evolution. Climatic-niche breadths for some ecologically important climatic variables considered separately (temperature seasonality and annual precipitation) do show significant relationships with the rate of climatic-niche evolution, but rates are faster in clades in which species have broader (not narrower) niche breadths. In summary, our results show that narrower niche breadths are not necessarily associated with faster rates of niche evolution.

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.

Winding up the molecular clock in the genus Carabus (Coleoptera: Carabidae): assessment of methodological decisions on rate and node age estimation

BACKGROUND

Rates of molecular evolution are known to vary across taxa and among genes, and this requires rate calibration for each specific dataset based on external information. Calibration is sensitive to evolutionary model parameters, partitioning schemes and clock model. However, the way in which these and other analytical aspects affect both the rates and the resulting clade ages from calibrated phylogenies are not yet well understood. To investigate these aspects we have conducted calibration analyses for the genus Carabus (Coleoptera, Carabidae) on five mitochondrial and four nuclear DNA fragments with 7888 nt total length, testing different clock models and partitioning schemes to select the most suitable using Bayes Factors comparisons.

RESULTS

We used these data to investigate the effect of ambiguous character and outgroup inclusion on both the rates of molecular evolution and the TMRCA of Carabus. We found considerable variation in rates of molecular evolution depending on the fragment studied (ranging from 5.02% in cob to 0.26% divergence/My in LSU-A), but also on analytical conditions. Alternative choices of clock model, partitioning scheme, treatment of ambiguous characters, and outgroup inclusion resulted in rate increments ranging from 28% (HUWE1) to 1000% (LSU-B and ITS2) and increments in the TMRCA of Carabus ranging from 8.4% (cox1-A) to 540% (ITS2). Results support an origin of the genus Carabus during the Oligocene in the Eurasian continent followed by a Miocene differentiation that originated all main extant lineages.

CONCLUSIONS

The combination of several genes is proposed as the best strategy to minimise both the idiosyncratic behaviors of individual markers and the effect of analytical aspects in rate and age estimations. Our results highlight the importance of estimating rates of molecular evolution for each specific dataset, selecting for optimal clock and partitioning models as well as other methodological issues potentially affecting rate estimation.

LTR retrotransposons contribute to genomic gigantism in plethodontid salamanders

Among vertebrates, most of the largest genomes are found within the salamanders, a clade of amphibians that includes 613 species. Salamander genome sizes range from ~14 to ~120 Gb. Because genome size is correlated with nucleus and cell sizes, as well as other traits, morphological evolution in salamanders has been profoundly affected by genomic gigantism. However, the molecular mechanisms driving genomic expansion in this clade remain largely unknown. Here, we present the first comparative analysis of transposable element (TE) content in salamanders. Using high-throughput sequencing, we generated genomic shotgun data for six species from the Plethodontidae, the largest family of salamanders. We then developed a pipeline to mine TE sequences from shotgun data in taxa with limited genomic resources, such as salamanders. Our summaries of overall TE abundance and diversity for each species demonstrate that TEs make up a substantial portion of salamander genomes, and that all of the major known types of TEs are represented in salamanders. The most abundant TE superfamilies found in the genomes of our six focal species are similar, despite substantial variation in genome size. However, our results demonstrate a major difference between salamanders and other vertebrates: salamander genomes contain much larger amounts of long terminal repeat (LTR) retrotransposons, primarily Ty3/gypsy elements. Thus, the extreme increase in genome size that occurred in salamanders was likely accompanied by a shift in TE landscape. These results suggest that increased proliferation of LTR retrotransposons was a major molecular mechanism contributing to genomic expansion in salamanders.