Bayesian Relaxed Molecular Clock Analyses Research Articles

2D approach to reconstruct the evolutionary history of clevelandellids (Ciliophora, Armophorea) inhabiting the hindgut of the Panesthiinae cockroaches

AbstractThe ciliate order Clevelandellida unites endosymbionts of the digestive tract of a variety of invertebrates and vertebrates. In the present study, the primary and secondary structures of nuclear and hydrogenosomal rRNA molecules were employed to reconstruct the phylogenetic relationships and to estimate the divergence times of clevelandellids inhabiting the hindgut of the Panesthiinae cockroaches. The secondary structure information was incorporated in phylogenetic analyses using two different strategies, viz., indirectly through 2D‐guided alignments and directly through so‐called pseudo‐protein data. Nuclear and hydrogenosomal markers carried a consistent phylogenetic signal and robustly supported the monophyletic origin of the family Clevelandellidae as well as of its four genera. According to Bayesian relaxed molecular clock analyses, the last common ancestor (LCA) of the family Clevelandellidae very likely emerged during the Late Cretaceous in the Oriental region. Its descendants most likely expanded to Australia in concert with the Neogene colonization and radiation of their host Panesthiinae cockroaches. Taking into account the time‐calibrated phylogenies and the fact that early branching members of the order Clevelandellida inhabit the digestive tracts of amphibians, it is tempting to speculate that the LCA of the Clevelandellida evolved in ectothermic vertebrates. Amphibians could have brought clevelandellids to the land, where they may have been transmitted to the cockroach digestive tract upon feeding on amphibian faeces.

The timescale of early land plant evolution

Establishing the timescale of early land plant evolution is essential for testing hypotheses on the coevolution of land plants and Earth's System. The sparseness of early land plant megafossils and stratigraphic controls on their distribution make the fossil record an unreliable guide, leaving only the molecular clock. However, the application of molecular clock methodology is challenged by the current impasse in attempts to resolve the evolutionary relationships among the living bryophytes and tracheophytes. Here, we establish a timescale for early land plant evolution that integrates over topological uncertainty by exploring the impact of competing hypotheses on bryophyte-tracheophyte relationships, among other variables, on divergence time estimation. We codify 37 fossil calibrations for Viridiplantae following best practice. We apply these calibrations in a Bayesian relaxed molecular clock analysis of a phylogenomic dataset encompassing the diversity of Embryophyta and their relatives within Viridiplantae. Topology and dataset sizes have little impact on age estimates, with greater differences among alternative clock models and calibration strategies. For all analyses, a Cambrian origin of Embryophyta is recovered with highest probability. The estimated ages for crown tracheophytes range from Late Ordovician to late Silurian. This timescale implies an early establishment of terrestrial ecosystems by land plants that is in close accord with recent estimates for the origin of terrestrial animal lineages. Biogeochemical models that are constrained by the fossil record of early land plants, or attempt to explain their impact, must consider the implications of a much earlier, middle Cambrian-Early Ordovician, origin.

Molecular phylogenetic analysis of nuclear genes suggests a Cenozoic over-water dispersal origin for the Cuban solenodon.

The Cuban solenodon (Solenodon cubanus) is one of the most enigmatic mammals and is an extremely rare species with a distribution limited to a small part of the island of Cuba. Despite its rarity, in 2012 seven individuals of S. cubanus were captured and sampled successfully for DNA analysis, providing new insights into the evolutionary origin of this species and into the origins of the Caribbean fauna, which remain controversial. We conducted molecular phylogenetic analyses of five nuclear genes (Apob, Atp7a, Bdnf, Brca1 and Rag1; total, 4,602 bp) from 35 species of the mammalian order Eulipotyphla. Based on Bayesian relaxed molecular clock analyses, the family Solenodontidae diverged from other eulipotyphlan in the Paleocene, after the bolide impact on the Yucatan Peninsula, and S. cubanus diverged from the Hispaniolan solenodon (S. paradoxus) in the Early Pliocene. The strikingly recent divergence time estimates suggest that S. cubanus and its ancestral lineage originated via over-water dispersal rather than vicariance events, as had previously been hypothesised.

Crown group Oxyphotobacteria postdate the rise of oxygen.

The rise of oxygen ca. 2.3billion years ago (Ga) is the most distinct environmental transition in Earth history. This event was enabled by the evolution of oxygenic photosynthesis in the ancestors of Cyanobacteria. However, long-standing questions concern the evolutionary timing of this metabolism, with conflicting answers spanning more than one billion years. Recently, knowledge of the Cyanobacteria phylum has expanded with the discovery of non-photosynthetic members, including a closely related sister group termed Melainabacteria, with the known oxygenic phototrophs restricted to a clade recently designated Oxyphotobacteria. By integrating genomic data from the Melainabacteria, cross-calibrated Bayesian relaxed molecular clock analyses show that crown group Oxyphotobacteria evolved ca. 2.0 billion years ago (Ga), well after the rise of atmospheric dioxygen. We further estimate the divergence between Oxyphotobacteria and Melainabacteria ca. 2.5-2.6Ga, which-if oxygenic photosynthesis is an evolutionary synapomorphy of the Oxyphotobacteria-marks an upper limit for the origin of oxygenic photosynthesis. Together, these results are consistent with the hypothesis that oxygenic photosynthesis evolved relatively close in time to the rise of oxygen.

Origin of marine planktonic cyanobacteria.

Marine planktonic cyanobacteria contributed to the widespread oxygenation of the oceans towards the end of the Pre-Cambrian and their evolutionary origin represents a key transition in the geochemical evolution of the Earth surface. Little is known, however, about the evolutionary events that led to the appearance of marine planktonic cyanobacteria. I present here phylogenomic (135 proteins and two ribosomal RNAs), Bayesian relaxed molecular clock (18 proteins, SSU and LSU) and Bayesian stochastic character mapping analyses from 131 cyanobacteria genomes with the aim to unravel key evolutionary steps involved in the origin of marine planktonic cyanobacteria. While filamentous cell types evolved early on at around 2,600–2,300 Mya and likely dominated microbial mats in benthic environments for most of the Proterozoic (2,500–542 Mya), marine planktonic cyanobacteria evolved towards the end of the Proterozoic and early Phanerozoic. Crown groups of modern terrestrial and/or benthic coastal cyanobacteria appeared during the late Paleoproterozoic to early Mesoproterozoic. Decrease in cell diameter and loss of filamentous forms contributed to the evolution of unicellular planktonic lineages during the middle of the Mesoproterozoic (1,600–1,000 Mya) in freshwater environments. This study shows that marine planktonic cyanobacteria evolved from benthic marine and some diverged from freshwater ancestors during the Neoproterozoic (1,000–542 Mya).

Phylogeography of cliff-dwelling relicts with a highly narrow and disjunct distribution in the western Mediterranean.

The high biodiversity in the Baetic-Rifan hotspot of Mediterranean region is shaped by complex geological and climatic histories and has been a subject of recent intensive studies. However, very little is known about phylogenetic and biogeographic history of three rare and critically endangered cliff-dwelling species of Sonchus in section Pustulati in this region. We investigated the genetic variation and phylogenetic relationships of populations based on nuclear (ITS/ETS) and plastid (3'trnL-ndhJ/psaI-accD) DNA sequences, and amplified fragment length polymorphisms (AFLPs). We performed a Bayesian relaxed molecular clock analysis with ITS data to estimate divergence times for major lineages. ITS/ETS and AFLP phylogenies showed high concordance and contrasted with cpDNA data. The divergence between S. masguindalii and S. fragilis was dated at 5.48 Ma, between S. fragilis and S. pustulatus at 3.89 Ma, and between the Baetic and Rifan S. pustulatus at 1.18 Ma. Within each distribution area, AFLP data showed a relatively high genetic structuring and moderate genetic diversity, the latter being impoverished in the Baetic populations. Our results further confirm the hybrid origin of S. pustulatus, a critically endangered species. The origin and diversification of lineages appear to have occurred on the temporary land bridge that joined Iberian and North Africa during the Messinian Salinity Crisis (5.96-5.33 Ma) and the subsequent Zanclean flood that progressively refilled the Mediterranean Basin (5.33-3.60 Ma). The only Baetic populations of S. pustulatus most likely originated from the Rifan ones.

Molecular phylogeny of Pompilinae (Hymenoptera: Pompilidae): Evidence for rapid diversification and host shifts in spider wasps

Pompilinae is one of the largest subfamilies of spider wasps (Pompilidae). Most pompilines are generalist spider predators at the family level, but some taxa exhibit ecological specificity (i.e., to spider-host guild). Here we present the first molecular phylogenetic analysis of Pompilinae, toward the aim of evaluating the monophyly of tribes and genera. We further test whether changes in the rate of diversification are associated with host-guild shifts. Molecular data were collected from five nuclear loci (28S, EF1-F2, LWRh, Wg, Pol2) for 76 taxa in 39 genera. Data were analyzed using maximum likelihood (ML) and Bayesian inference (BI). The phylogenetic results were compared with previous hypotheses of subfamilial and tribal classification, as well as generic relationships in the subfamily. The classification of Pompilus and Agenioideus is also discussed. A Bayesian relaxed molecular clock analysis was used to examine divergence times. Diversification rate-shift tests accounted for taxon-sampling bias using ML and BI approaches. Ancestral host family and host guild were reconstructed using MP and ML methods. Ancestral host guild for all Pompilinae, for the ancestor at the node where a diversification rate-shift was detected, and two more nodes back in time was inferred using BI. In the resulting phylogenies, Aporini was the only previously proposed monophyletic tribe. Several genera (e.g., Pompilus, Microphadnus and Schistonyx) are also not monophyletic. Dating analyses produced a well-supported chronogram consistent with topologies from BI and ML results. The BI ancestral host-use reconstruction inferred the use of spiders belonging to the guild “other hunters” (frequenting the ground and vegetation) as the ancestral state for Pompilinae. This guild had the highest probability for the ML reconstruction and was equivocal for the MP reconstruction; various switching events to other guilds occurred throughout the evolution of the group. The diversification of Pompilinae shows one main rate-shift coinciding with a shift to ground-hunter spiders, as reconstructed by the BI ancestral character-state analysis.

Phylogenetic trait conservation in the partner choice of a group of ectomycorrhizal trees.

Ecological interactions are frequently conserved across evolutionary time. In the case of mutualisms, these conserved interactions may play a large role in structuring mutualist communities. We hypothesized that phylogenetic trait conservation could play a key role in determining patterns of association in the ectomycorrhizal symbiosis, a globally important trophic mutualism. We used the association between members of the pantropical plant tribe Pisonieae and its fungal mutualist partners as a model system to test the prediction that Pisonieae-associating ectomycorrhizal fungi will be more closely related than expected by chance, reflecting a conserved trait. We tested this prediction using previously published and newly generated sequences in a Bayesian framework incorporating phylogenetic uncertainty. We report that phylogenetic trait conservation does exist in this association. We generated a five-marker phylogeny of members of the Pisonieae and used this phylogeny in a Bayesian relaxed molecular clock analysis. We established that the most recent common ancestors of Pisonieae species and Pisonieae-associating fungi sharing phylogenetic conservation of their patterns of ectomycorrhizal association occurred no more recently than 14.2 Ma. We therefore suggest that phylogenetic trait conservation in the Pisonieae ectomycorrhizal mutualism association represents an inherited syndrome which has existed for at least 14 Myr.

Phylogenetic relationships and species delimitation in pinus section trifoliae inferrred from plastid DNA.

Recent diversification followed by secondary contact and hybridization may explain complex patterns of intra- and interspecific morphological and genetic variation in the North American hard pines (Pinus section Trifoliae), a group of approximately 49 tree species distributed in North and Central America and the Caribbean islands. We concatenated five plastid DNA markers for an average of 3.9 individuals per putative species and assessed the suitability of the five regions as DNA bar codes for species identification, species delimitation, and phylogenetic reconstruction. The ycf1 gene accounted for the greatest proportion of the alignment (46.9%), the greatest proportion of variable sites (74.9%), and the most unique sequences (75 haplotypes). Phylogenetic analysis recovered clades corresponding to subsections Australes, Contortae, and Ponderosae. Sequences for 23 of the 49 species were monophyletic and sequences for another 9 species were paraphyletic. Morphologically similar species within subsections usually grouped together, but there were exceptions consistent with incomplete lineage sorting or introgression. Bayesian relaxed molecular clock analyses indicated that all three subsections diversified relatively recently during the Miocene. The general mixed Yule-coalescent method gave a mixed model estimate of only 22 or 23 evolutionary entities for the plastid sequences, which corresponds to less than half the 49 species recognized based on morphological species assignments. Including more unique haplotypes per species may result in higher estimates, but low mutation rates, recent diversification, and large effective population sizes may limit the effectiveness of this method to detect evolutionary entities.

A divergence dating analysis of turtles using fossil calibrations: an example of best practices

Turtles have served as a model system for molecular divergence dating studies using fossil calibrations. However, because some parts of the fossil record of turtles are very well known, divergence age estimates from molecular phylogenies often do not differ greatly from those observed directly from the fossil record alone. Also, the phylogenetic position and age of turtle fossil calibrations used in previous studies have not been adequately justified. We provide the first explicitly justified minimum and soft maximum age constraints on 22 clades of turtles following best practice protocols. Using these data we undertook a Bayesian relaxed molecular clock analysis establishing a timescale for the evolution of crown Testudines that we exploit in attempting to address evolutionary questions that cannot be resolved with fossils alone. Some of these questions, such as whether the turtle crown originated in the Triassic or Jurassic, cannot be resolved by our analysis. However, our results generate novel age-of-origination estimates for clades within crown Testudines. Finally, we compare our fossil calibrations and posterior age estimates to those from other studies, revealing substantial differences in results and interpretation.

Phylogenetic relationships of the mycoheterotrophic genus Voyria and the implications for the biogeographic history of Gentianaceae

The angiosperm family Gentianaceae comprises over 1700 species in 91 genera. Gentianaceae are distributed worldwide, but most species occur in temperate zones. Phylogenetic studies demonstrate that the family consists of six monophyletic tribes. However, the phylogenetic position of the mycoheterotrophic genus Voyria, with a remarkable trans-Atlantic distribution, remained unknown. We used nuclear ITS and 18S rDNA and mitochondrial apt1 and matR data to infer the phylogenetic position of Voyria in Gentianaceae. In addition, with Bayesian relaxed molecular clock analyses we obtained age estimates for the diversification of Voyria and Gentianaceae in general and used these results to reconstruct the ancestral areas associated with the early diversification events in Gentianaceae. Our results demonstrate that Voyria is an early diverging lineage within Gentianaceae with no close relationships to other mycoheterotrophic Gentianaceae lineages. Voyria originated in the neotropics during the Early Eocene but only reached its current transoceanic distribution around the end of the Oligocene. The neotropics were an important area for the early diversification events in Gentianaceae, most of which occurred during the Eocene. Voyria is an old, phylogenetically isolated lineage within Gentianaceae, and the current distribution of the genus is indicative of the ancestral area in which the early diversification events of Gentianaceae occurred. In parallel with many other pantropical families, our results suggest that migration of tropical taxa through Laurasia during the Early Eocene has played an important role in shaping the current global distribution of Gentianaceae.

Trapped in desert springs: phylogeography of Australian desert spring snails

AbstractAim We investigate the phylogeographical history and determine the time‐scale of population divergence of hydrobiid freshwater snails (genus Trochidrobia) inhabiting groundwater springs in the Australian desert. We test the hypothesis that divergence between geographically distinct snail populations occurred simultaneously due to their isolation in hydrologically discrete spring systems, i.e. ‘trapped in desert springs’.Location Groundwater springs of the Great Artesian Basin (GAB) in central Australia.Methods DNA sequence data from the mitochondrial cytochrome c oxidase subunit I gene and the nuclear 28S and internal transcribed spacer rRNA genes were used to reconstruct phylogenetic relationships within and among three species of Trochidrobia (Hydrobiidae): T. punicea (13 spring groups, n = 90), T. smithi (12 spring groups, n = 62) and T. minuta (2 spring groups, n = 4). Bayesian relaxed molecular clock analyses and approximate Bayesian computation were used to date lineage divergence and distinguish between alternative biogeographical scenarios.Results The diversification of the three Trochidrobia species probably occurred between 2.54 and 9.3 Ma, prior to the formation of the springs c. 1 Ma. Intraspecific divergences within the two widespread species occurred after the formation and colonization of the springs. Coalescent modelling and molecular clock analyses supported a simultaneous radiation of five allopatric intraspecific snail lineages within T. punicea (two lineages) and T. smithi (three lineages) across the GAB springs examined.Main conclusions The analyses support the ‘trapped in desert springs’ hypothesis for the diversification of intraspecific lineages within the species T. punicea and T. smithi. This hypothesis suggests that the formation of deserts around Lake Eyre in the early Pleistocene led to the hydrological isolation of spring complexes in the GAB, resulting in significant molecular divergence, but no morphological divergence, of Trochidrobia snail populations.

Understanding the Evolutionary Processes of Fungal Fruiting Bodies: Correlated Evolution and Divergence Times in the Psathyrellaceae

Fruiting body evolution is one of the central topics in fungal evolutionary biology. A number of hypotheses have been developed to explain the contemporary diversity of fruiting body forms, but their evaluation has been hampered by the lack of well-sampled data sets and suitable statistical methods. Phylogenetic evidence of the physiological changes that accompany switches in fruiting body type is lacking, and very little is known about the age of major events of fruiting body evolution. Based on a new multigene phylogeny, by using Bayesian methods, we demonstrate the existence of correlation between a number of morphological features and switches from nondeliquescent to deliquescent (autodigesting) fruiting bodies in the mushroom family Psathyrellaceae. Our results show that switches in the anatomy of two types of spacer cells (cystidia and pseudoparaphyses) and basidia (bimorphic or monomorphic) as well as the structure of the mushroom cap follow the evolution of deliquescent fruiting bodies, which suggests strong functional linkage between these traits. We performed Bayes factor-based tests, referred hereafter to as evolutionary pathway test (EPT), to decide which of the correlated characters were gained first during evolution. The EPTs strongly suggest that deliquescence was gained first, followed after short waiting times by the other morphological features. Bayesian relaxed molecular clock analyses suggest that the various events of switching between fruiting body types occurred independently at various ages during the history of the family. The utility of two mushroom fossils (Archaemarasmius and Protomycena), the only ones with unambiguous taxonomic positions, for the calibration of agaric trees were also examined. Based on our results, we suggest that the evolutionary benefit of deliquescence may be prevention against desiccation via accelerated ontogeny of the fruiting body. Hypotheses regarding the functional significance of the correlated evolution are presented and discussed. Further, we argue that the changes in fruiting body types in mushrooms in general can be attributed to independent events (e.g., dispersal and adaptation) and not to particular geologic ages.

A multi-locus time-calibrated phylogeny of the brown algae (Heterokonta, Ochrophyta, Phaeophyceae): Investigating the evolutionary nature of the “brown algal crown radiation”

The most conspicuous feature in previous phaeophycean phylogenies is a large polytomy known as the brown algal crown radiation (BACR). The BACR encompasses 10 out of the 17 currently recognized brown algal orders. A recent study has been able to resolve a few nodes of the BACR, suggesting that it may be a soft polytomy caused by a lack of signal in molecular markers. The present work aims to refine relationships within the BACR and investigate the nature and timeframe of the diversification in question using a dual approach. A multi-marker phylogeny of the brown algae was built from 10 mitochondrial, plastid and nuclear loci (>10,000 nt) of 72 phaeophycean taxa, resulting in trees with well-resolved inter-ordinal relationships within the BACR. Using Bayesian relaxed molecular clock analysis, it is shown that the BACR is likely to represent a gradual diversification spanning most of the Lower Cretaceous rather than a sudden radiation. Non-molecular characters classically used in ordinal delimitation were mapped on the molecular topology to study their evolutionary history.

Molecular evolution of glutamine synthetase II: Phylogenetic evidence of a non-endosymbiotic gene transfer event early in plant evolution.

BackgroundGlutamine synthetase (GS) is essential for ammonium assimilation and the biosynthesis of glutamine. The three GS gene families (GSI, GSII, and GSIII) are represented in both prokaryotic and eukaryotic organisms. In this study, we examined the evolutionary relationship of GSII from eubacterial and eukaryotic lineages and present robust phylogenetic evidence that GSII was transferred from γ-Proteobacteria (Eubacteria) to the Chloroplastida.ResultsGSII sequences were isolated from four species of green algae (Trebouxiophyceae), and additional green algal (Chlorophyceae and Prasinophytae) and streptophyte (Charales, Desmidiales, Bryophyta, Marchantiophyta, Lycopodiophyta and Tracheophyta) sequences were obtained from public databases. In Bayesian and maximum likelihood analyses, eubacterial (GSIIB) and eukaryotic (GSIIE) GSII sequences formed distinct clades. Both GSIIB and GSIIE were found in chlorophytes and early-diverging streptophytes. The GSIIB enzymes from these groups formed a well-supported sister clade with the γ-Proteobacteria, providing evidence that GSIIB in the Chloroplastida arose by horizontal gene transfer (HGT). Bayesian relaxed molecular clock analyses suggest that GSIIB and GSIIE coexisted for an extended period of time but it is unclear whether the proposed HGT happened prior to or after the divergence of the primary endosymbiotic lineages (the Archaeplastida). However, GSIIB genes have not been identified in glaucophytes or red algae, favoring the hypothesis that GSIIB was gained after the divergence of the primary endosymbiotic lineages. Duplicate copies of the GSIIB gene were present in Chlamydomonas reinhardtii, Volvox carteri f. nagariensis, and Physcomitrella patens. Both GSIIB proteins in C. reinhardtii and V. carteri f. nagariensis had N-terminal transit sequences, indicating they are targeted to the chloroplast or mitochondrion. In contrast, GSIIB proteins of P. patens lacked transit sequences, suggesting a cytosolic function. GSIIB sequences were absent in vascular plants where the duplication of GSIIE replaced the function of GSIIB.ConclusionsPhylogenetic evidence suggests GSIIB in Chloroplastida evolved by HGT, possibly after the divergence of the primary endosymbiotic lineages. Thus while multiple GS isoenzymes are common among members of the Chloroplastida, the isoenzymes may have evolved via different evolutionary processes. The acquisition of essential enzymes by HGT may provide rapid changes in biochemical capacity and therefore be favored by natural selection.

Molecular phylogeny of Solms‐laubachia (Brassicaceae) s.l., based on multiple nuclear and plastid DNA sequences, and its biogeographic implications

Abstract The Hengduan Mountains region of south‐west China is a noted biodiversity hotspot, but the geographic origins and historical assembly of its rich endemic flora, including the sky‐island species of Solms‐laubachia Muschl. (Brassicaceae), have been little studied. Previous molecular studies on the phylogeny of Solms‐laubachia showed it to be paraphyletic, leading to considerable expansion not only of its taxonomic limits, but also its geographic range, with the inclusion of taxa from outside the Hengduan region. However, these studies provided little resolution of interspecific relationships, preventing inferences about historical biogeography within the clade. In the present study, new sequence data from two nuclear genes (LEAFY and G3pdh) and two chloroplast intergenic spacers (petN–psbM and psbM–trnD) were combined with existing markers to increase phylogenetic signals. Phaeonychium villosum (Maxim.) Al‐Shehbaz was found to be nested within Solms‐laubachia s.l. In general, phylogenetic relationships appear to be a good predictor of geography, with the Hengduan Mountain endemics embedded in a paraphyletic grade of species from the western Himalayas and central Asia, but they also imply morphological homoplasy. Incongruence was detected between the nuclear and chloroplast gene trees, perhaps resulting from incomplete lineage sorting of ancestral polymorphisms. The crown age of Solms‐laubachia s.l. was estimated to be approximately 1.42–3.68 mya, using Bayesian relaxed molecular clock analysis. Historical biogeographic analysis using a parametric dispersal–extinction–cladogenesis model inferred central Asia and the western Himalayas as most probable ancestral range of Solms‐laubachia s.l., and estimated higher rates of eastward expansion than westward during the diversification of descendant lineages. In summary, our results suggest that Solms‐laubachia s.l. originated during the Pliocene in central Asia, and subsequently migrated eastward into the Hengduan Mountains, colonizing sky‐island, alpine scree‐slope habitats that may have provided novel ecological opportunity and accelerated speciation, ultimately establishing this region as the present center of diversity of the genus.

Deciphering and dating the red panda’s ancestry and early adaptive radiation of Musteloidea

Few species have been of more disputed affinities than the red or lesser panda ( Ailurus fulgens), an endangered endemic Southeast Asian vegetarian member of the placental mammalian order Carnivora. This peculiar carnivoran has mostly been classified with raccoons (Procyonidae) or bears (Ursidae), grouped with the giant panda ( Ailuropoda melanoleuca) in their own family, or considered a separate lineage of equivocal ancestry. Recent molecular studies have indicated a close affinity of the red panda to a clade of procyonids and mustelids (weasels, otters, martens, badgers, and allies), but have failed to unambiguously resolve the position of this species relative to mephitids (skunks and stink badgers). We examined the relationship of the red panda to other extant species of the carnivoran suborder Caniformia using a set of concatenated ∼5.5-kb sequences from protein-coding exons of five nuclear genes. Bayesian, maximum likelihood, and parsimony phylogenetic analyses strongly supported the red panda as the closest living relative of a clade containing Procyonidae and Mustelidae to the exclusion of Mephitidae. These three families together with the red panda (which is classified here as a single extant species of a distinct family, Ailuridae) compose the superfamily Musteloidea, a clade strongly supported by all our phylogenetic analyses as sister to the monophyletic Pinnipedia (seals, sea lions, walruses). The approximately unbiased, Kishino–Hasegawa, and Templeton topology tests rejected ( P < 0.05) each of all possible alternative hypotheses about the relationships among the red panda and mephitids, procyonids, and mustelids. We also estimated divergence times for the red panda’s lineage and ones of other caniform taxa, as well as the ages of the first appearance datums for the crown and total clades of musteloids and the total clades of the red panda, mephitids, procyonids, and mustelids. Bayesian relaxed molecular-clock analysis using combined information from all sampled genes yielded a ∼42-Myr timescale to caniform evolution and provided evidence of five periods of increased diversification. The red panda’s lineage and those of other extant musteloid families are estimated to have diverged during a 3-Myr interval from the mid-Early Oligocene to near the Early/Late Oligocene boundary. We present fossil evidence that extends the early adaptive radiation of the total clade of musteloids to the Eocene–Oligocene transition and also suggests Asia as a center of this radiation.

Divergence time of the two regional medaka populations in Japan as a new time scale for comparative genomics of vertebrates

The southern and northern Japanese populations of the medaka fish provide useful tools to gain insights into the comparative genomics and speciation of vertebrates, because they can breed to produce healthy and fertile offspring despite their highly divergent genetic backgrounds compared with those of human-chimpanzee. Comparative genomics analysis has suggested that such large genetic differences between the two populations are caused by higher molecular evolutionary rates among the medakas than those of the hominids. The argument, however, was based on the assumption that the two Japanese populations diverged approximately at the same time (4.0-4.7 Myr ago) as the human-chimpanzee lineage (5.0-6.0 Myr ago). This can be misleading, because the divergence time of the two populations was calculated based on estimated, extremely higher molecular evolutionary rates of other fishes with an implicit assumption of a global molecular clock. Here we show that our estimate, based on a Bayesian relaxed molecular-clock analysis of whole mitogenome sequences from 72 ray-finned fishes (including 14 medakas), is about four times older than that of the previous study (18 Myr). This remarkably older estimate can be reconciled with the vicariant events of the Japanese archipelago, and the resulting rates of molecular evolution are almost identical between the medaka and hominid lineages. Our results further highlight the fact that reproductive isolation may not evolve despite a long period of geographical isolation.

The relative ages of ectomycorrhizal mushrooms and their plant hosts estimated using Bayesian relaxed molecular clock analyses.

BackgroundEctomycorrhizae (ECM) are symbioses formed by polyphyletic assemblages of fungi (mostly Agaricomycetes) and plants (mostly Pinaceae and angiosperms in the rosid clade). Efforts to reconstruct the evolution of the ECM habit in Agaricomycetes have yielded vastly different results, ranging from scenarios with many relatively recent origins of the symbiosis and no reversals to the free-living condition; a single ancient origin of ECM and many subsequent transitions to the free-living condition; or multiple gains and losses of the association. To test the plausibility of these scenarios, we performed Bayesian relaxed molecular clock analyses including fungi, plants, and other eukaryotes, based on the principle that a symbiosis cannot evolve prior to the origin of both partners. As we were primarily interested in the relative ages of the plants and fungi, we did not attempt to calibrate the molecular clock using the very limited fossil record of Agaricomycetes.ResultsTopologically constrained and unconstrained analyses suggest that the root node of the Agaricomycetes is much older than either the rosids or Pinaceae. The Agaricomycetidae, a large clade containing the Agaricales and Boletales (collectively representing 70% of Agaricomycetes), is also significantly older than the rosids. The relative age of Agaricomycetidae and Pinaceae, however, is sensitive to tree topology, and the inclusion or exclusion of the gnetophyte Welwitschia mirabilis.ConclusionThe ancestor of the Agaricomycetes could not have been an ECM species because it existed long before any of its potential hosts. Within more derived clades of Agaricomycetes, there have been at least eight independent origins of ECM associations involving angiosperms, and at least six to eight origins of associations with gymnosperms. The first ECM symbioses may have involved Pinaceae, which are older than rosids, but several major clades of Agaricomycetes, such as the Boletales and Russulales, are young enough to have been plesiomorphically associated with either rosids or Pinaceae, suggesting that some contemporary ECM partnerships could be of very ancient origin.

Comparative phylogeography of livebearing fishes in the genera Poeciliopsis and Poecilia (Poeciliidae: Cyprinodontiformes) in central Mexico

AbstractAim To test the hypothesis that the vicariant event responsible for north–south divergences in two clades of the fish genus Poeciliopsis Regan was also responsible for north–south divergences in the fish Poecilia butleri Jordan.Location Central Mexico.Methods Parsimony, distance, maximum likelihood and Bayesian phylogenetic analyses of a mitochondrial gene. Molecular clock test and Bayesian analyses of divergence time.Results We report concordant phylogeographical patterns between two clades in the genus Poeciliopsis (i.e. the other formed by P. latidens Garman and P. fasciata Meek, and the other formed by P. presidionis Jordan and P. turneri Miller) and the clade of Poecilia butleri, with northern and southern individuals within each clade grouping into separate lineages. There is also evidence for slower substitution rates in Poecilia compared with Poeciliopsis. After taking into account these substitution rate discrepancies with Bayesian relaxed molecular clock analyses, north–south divergences in Poecilia butleri were equivalent to those reported for Poeciliopsis latidens‐fasciata and P. presidionis‐turneri.Main conclusions The same Plio‐Pleistocene vicariant event associated with geological activity of the Trans‐Mexican Volcanic Belt (TMVB) appears to have caused divergence in three different freshwater fish lineages. This study is an example of how comparative phylogeography can strengthen inferences about vicariant events in regions of high biological diversity and complex geological history such as the TMVB.