Evaluating fossil calibrations for dating phylogenies in light of rates of molecular evolution: a comparison of three approaches

Diversity and biogeography of South American mud turtles elucidated by multilocus DNA sequencing (Testudines: Kinosternidae)

Kinosternon is the most speciose genus of extant turtles, with 22 currently recognized species, distributed across large parts of the Americas. Most species have small distributions, but K. leucostomum and K. scorpioides range from Mexico to South America. Previous studies have found discordance between mitochondrial and nuclear phylogenies in some kinosternid groups, with the current taxonomy following the nuclear-based results. Herein, based on extended molecular, geographic, and taxonomic sampling, we explore the phylogeographic structure and taxonomic limits for K. leucostomum and the K. scorpioides group and present a fossil-calibrated nuclear time tree for Kinosternon. Our results reveal contrasting differentiation patterns for the K. scorpioides group and K. leucostomum, despite overlapping distributions. Kinosternon leucostomum shows only shallow geographic divergence, whereas the K. scorpioides group is polyphyletic with up to 10 distinct taxa, some of them undescribed. We support the elevation of K. s. albogulare and K. s. cruentatum to species level. Given the deep divergence within the genus Kinosternon, we propose the recognition of three subgenera, Kinosternon, Cryptochelys and Thyrosternum, and the abandonment of the group-based classification, at least for the K. leucostomum and K. scorpioides groups. Our results show an initial split in Kinosternon that gave rise to two main radiations, one Nearctic and one mainly Neotropical. Most speciation events in Kinosternon occurred during the Quaternary and we hypothesize that they were mediated by both climatic and geological events. Additionally, our data imply that at least three South American colonizations occurred, two in the K. leucostomum group, and one in the K. scorpioides group. Additionally, we hypothesize that discordance between mitochondrial and nuclear phylogenetic signal is due to mitochondrial capture from an extinct kinosternine lineage.

A comprehensive molecular phylogeny of Cephalotrichum and Microascus provides novel insights into their systematics and evolutionary history

The genera Cephalotrichum and Microascus contain ecologically, morphologically and lifestyle diverse fungi in Microascaceae (Microascales, Sordariomycetes) with a world-wide distribution. Despite previous studies having elucidated that Cephalotrichum and Microascus are highly polyphyletic, the DNA phylogeny of many traditionally morphology-defined species is still poorly resolved, and a comprehensive taxonomic overview of the two genera is lacking. To resolve this issue, we integrate broad taxon sampling strategies and the most comprehensive multi-gene (ITS, LSU, tef1 and tub2) datasets to date, with fossil calibrations to address the phylogenetic relationships and divergence times among major lineages of Microascaceae. Two previously recognised main clades, Cephalotrichum (24 species) and Microascus (49 species), were re-affirmed based on our phylogenetic analyses, as well as the phylogenetic position of 15 genera within Microascaceae. In this study, we provide an up-to-date overview on the taxonomy and phylogeny of species belonging to Cephalotrichum and Microascus, as well as detailed descriptions and illustrations of 21 species of which eight are newly described. Furthermore, the divergence time estimates indicate that the crown age of Microascaceae was around 210.37 Mya (95 % HPD: 177.18-246.96 Mya) in the Late Triassic, and that Cephalotrichum and Microascus began to diversify approximately 27.07 Mya (95 % HPD: 20.47-34.37 Mya) and 70.46 Mya (95 % HPD: 56.96-86.24 Mya), respectively. Our results also demonstrate that multigene sequence data coupled with broad taxon sampling can help elucidate previously unresolved clade relationships. Citation: Wei TP, Wu YM, Zhang X, et al. 2024. A comprehensive molecular phylogeny of Cephalotrichum and Microascus provides novel insights into their systematics and evolutionary history. Persoonia 52: 119-160. https://doi.org/10.3767/persoonia.2024.52.05 .

Phylogenetic conservation in plant phenological traits varies between temperate and subtropical climates in China

Phenological traits, such as leaf and flowering dates, are proven to be phylogenetically conserved. The relationship between phylogenetic conservation, plant phenology, and climatic factors remains unknown. Here, we assessed phenological features among flowering plants as evidence for phylogenetic conservatism, the tendency for closely related species to share similar ecological and biological attributes. We use spring phenological traits data from 1968-2018 of 65 trees and 49 shrubs in Xi'an (temperate climate) and Guiyang (subtropical climate) to understand plant phenological traits' relationship with phylogeny. Molecular datasets are employed in evolutionary models to test the phylogenetic conservatism in spring phenological characteristics in response to climate-sensitive phenological features. Significant phylogenetic conservation was found in the Xi'an plant's phenological traits, while there was a non-significant conservation in the Guiyang plant species. Phylogenetic generalized least squares (PGLS) models correlate with phenological features significantly in Xi'an while non-significantly in Guiyang. Based on the findings of molecular dating, it was suggested that the Guiyang species split off from their relatives around 46.0 mya during the middle Eocene of the Tertiary Cenozoic Era, while Xi'an species showed a long evolutionary history and diverged from their relatives around 95 mya during the late Cretaceous Mesozoic Era. First leaf dates (FLD) indicative of spring phenology, show that Xi'an adjourned the case later than Guiyang. Unlike FLD, first flower dates (FFD) yield different results as Guiyang flowers appear later than Xi'an's. Our research revealed that various factors, including phylogeny, growth form, and functional features, influenced the diversity of flowering phenology within species in conjunction with local climate circumstances. These results are conducive to understanding evolutionary conservation mechanisms in plant phenology concerning evolutionary processes in different geographical and climate zones.

Is it inappropriate to ask for your age? Evaluating parameter impact on tree dating in a challenging clade (Macroscelidea)

Sengis (order Macroscelidea) are small mammals endemic to Africa. The taxonomy and phylogeny of sengis has been difficult to resolve due to a lack of clear morphological apomorphies. Molecular phylogenies have already significantly revised sengi systematics, but until now no molecular phylogeny has included all 20 extant species. In addition, the age of origin of the sengi crown clade and the divergence age of its two extant families remain unclear. Two recently published studies based on different datasets and age-calibration parameters (DNA type, outgroup selection, fossil calibration points) proposed highly different divergent age estimates and evolutionary scenarios. We obtained nuclear and mitochondrial DNA from mainly museum specimens using target enrichment of single-stranded DNA libraries to generate the first phylogeny of all extant macroscelidean species. We then explored the effects of different parameters (type of DNA, ratio of ingroup to outgroup sampling, number and type of fossil calibration points) and their resulting impacts on age estimates for the origin and initial diversification of Macroscelidea. We show that, even after correcting for substitution saturation, both using mitochondrial DNA in conjunction with nuclear DNA or alone results in much older ages and different branch lengths than when using nuclear DNA alone. We further show that the former effect can be attributed to insufficient amounts of nuclear data. If multiple calibration points are included, the age of the sengi crown group fossil prior has minimal impact on the estimated time frame of sengi evolution. In contrast, the inclusion or exclusion of outgroup fossil priors has a major effect on the resulting node ages. We also find that a reduced sampling of ingroup species does not significantly affect overall age estimates and that terminal specific substitution rates can serve as a means to evaluate the biological likeliness of the produced temporal estimates. Our study demonstrates how commonly varied parameters in temporal calibration of phylogenies affect age estimates. Dated phylogenies should therefore always be seen in the context of the dataset which was used to produce them.

Molecular phylogeny of the tropical wandering spiders (Araneae, Ctenidae) and the evolution of eye conformation in the RTA clade

Tropical wandering spiders (Ctenidae) are a diverse group of cursorial predators with its greatest species richness in the tropics. Traditionally, Ctenidae are diagnosed based on the presence of eight eyes arranged in three rows (a 2-4-2 pattern). We present a molecular phylogeny of Ctenidae, including for the first time representatives of all of its subfamilies. The molecular phylogeny was inferred using five nuclear (histone H3, 28S, 18S, Actin and ITS-2) and four mitochondrial (NADH, COI, 12S and 16S) markers. The final matrix includes 259 terminals, 103 of which belong to Ctenidae and represent 28 of the current 49 described genera. We estimated divergence times by including fossils as calibration points and biogeographic events, and used the phylogenetic hypothesis obtained to reconstruct the evolution of the eye conformation in the retrolateral tibial apophysis (RTA) clade. Ctenidae and its main lineages originated during the Paleocene-Eocene and have diversified in the tropics since then. However, in some analyses Ctenidae was recovered as polyphyletic as the genus Ancylometes Bertkau, 1880 was placed as sister to Oxyopidae. Except for Acantheinae, in which the type genus Acantheis Thorell, 1891 is placed inside Cteninae, the four recognized subfamilies of Ctenidae are monophyletic in most analyses. The ancestral reconstruction of the ocular conformation in the retrolateral tibial apophysis clade suggests that the ocular pattern of Ctenidae has evolved convergently seven times and that it has originated from ocular conformations of two rows of four eyes (4-4) and the ocular pattern of lycosids (4-2-2). We also synonymize the monotypic genus Parabatinga Polotov & Brescovit, 2009 with Centroctenus Mello-Leitão, 1929. We discuss some of the putative morphological synapomorphies of the main ctenid lineages within the phylogenetic framework offered by the molecular phylogenetic results of the study.

Genomic diversity and Demographic History of the Dromiciops genus (Marsupialia: Microbiotheriidae)

Three orders represent the South American fauna of marsupials. Of these, Microbiotheria was until recently known as a monotypic genus with the only surviving species Dromiciops gliroides (monito del monte). The recent proposal of a new Dromiciops species (Dromiciops bicinovici), together with new information on the origin and diversification of living microbioterians has changed the prevailing paradigm around the evolutionary history of these emblematic marsupials. Here, we used a RADseq approach to test for evidence of admixture and past or current gene flow among both species of Dromiciops and evaluate the genetic structure within D. gliroides. We analyzed 127 samples of Dromiciops distributed across the known distribution range of both species. We also inferred the joint demographic history of these lineages, thus corroborating the status of D. bozinovici as a distinct species. Demographic history reconstruction indicated that D. bozinovici diverged from D. gliroides around 4my ago and has remained isolated and demographically stable ever since. In contrast, D. gliroides is subdivided into three subclades that experienced recent expansions and moderate gene flow among them (mostly from north to south). Furthermore, genetic distances among populations within D. gliroides were significantly correlated with geographic distances. These results suggest that some of the D. gliroides populations would have survived in glacial refuges, with posterior expansions after ice retreat. Our results have important implications for the systematics of the genus and have profound conservation consequences for the new species, especially considering the fragmentation level of the temperate rainforest.

The phylogeny, phylogeography, and diversification history of the westernmost Asian cobra (Serpentes: Elapidae: Naja oxiana) in the Trans-Caspian region

We conducted a comprehensive analysis of the phylogenetic, phylogeographic, and demographic relationships of Caspian cobra (Naja oxiana; Eichwald, 1831) populations based on a concatenated dataset of two mtDNA genes (cyt b and ND4) across the species' range in Iran, Afghanistan, and Turkmenistan, along with other members of Asian cobras (i.e., subgenus Naja Laurenti, 1768). Our results robustly supported that the Asiatic Naja are monophyletic, as previously suggested by other studies. Furthermore, N. kaouthia and N. sagittifera were recovered as sister taxa to each other, and in turn sister clades to N. oxiana. Our results also highlighted the existence of a single major evolutionary lineage for populations of N. oxiana in the Trans-Caspian region, suggesting a rapid expansion of this cobra from eastern to western Asia, coupled with a rapid range expansion from east of Iran toward the northeast. However, across the Iranian range of N. oxiana, subdivision of populations was not supported, and thus, a single evolutionary significant unit is proposed for inclusion in future conservation plans in this region.

Climatic and topographic changes since the Miocene influenced the diversification and biogeography of the tent tortoise (Psammobates tentorius) species complex in Southern Africa

BACKGROUND

Climatic and topographic changes function as key drivers in shaping genetic structure and cladogenic radiation in many organisms. Southern Africa has an exceptionally diverse tortoise fauna, harbouring one-third of the world's tortoise genera. The distribution of Psammobates tentorius (Kuhl, 1820) covers two of the 25 biodiversity hotspots in the world, the Succulent Karoo and Cape Floristic Region. The highly diverged P. tentorius represents an excellent model species for exploring biogeographic and radiation patterns of reptiles in Southern Africa.

RESULTS

We investigated genetic structure and radiation patterns against temporal and spatial dimensions since the Miocene in the Psammobates tentorius species complex, using multiple types of DNA markers and niche modelling analyses. Cladogenesis in P. tentorius started in the late Miocene (11.63-5.33 Ma) when populations dispersed from north to south to form two geographically isolated groups. The northern group diverged into a clade north of the Orange River (OR), followed by the splitting of the group south of the OR into a western and an interior clade. The latter divergence corresponded to the intensification of the cold Benguela current, which caused western aridification and rainfall seasonality. In the south, tectonic uplift and subsequent exhumation, together with climatic fluctuations seemed responsible for radiations among the four southern clades since the late Miocene. We found that each clade occurred in a habitat shaped by different climatic parameters, and that the niches differed substantially among the clades of the northern group but were similar among clades of the southern group.

CONCLUSION

Climatic shifts, and biome and geographic changes were possibly the three major driving forces shaping cladogenesis and genetic structure in Southern African tortoise species. Our results revealed that the cladogenesis of the P. tentorius species complex was probably shaped by environmental cooling, biome shifts and topographic uplift in Southern Africa since the late Miocene. The Last Glacial Maximum (LGM) may have impacted the distribution of P. tentorius substantially. We found the taxonomic diversify of the P. tentorius species complex to be highest in the Greater Cape Floristic Region. All seven clades discovered warrant conservation attention, particularly Ptt-B-Ptr, Ptt-A and Pv-A.

Phylogenomic Reconstruction Sheds Light on New Relationships and Timescale of Rails (Aves: Rallidae) Evolution

The integration of state-of-the-art molecular techniques and analyses, together with a broad taxonomic sampling, can provide new insights into bird interrelationships and divergence. Despite their evolutionary significance, the relationships among several rail lineages remain unresolved as does the general timescale of rail evolution. Here, we disentangle the deep phylogenetic structure of rails using anchored phylogenomics. We analysed a set of 393 loci from 63 species, representing approximately 40% of the extant familial diversity. Our phylogenomic analyses reconstruct the phylogeny of rails and robustly infer several previously contentious relationships. Concatenated maximum likelihood and coalescent species-tree approaches recover identical topologies with strong node support. The results are concordant with previous phylogenetic studies using small DNA datasets, but they also supply an additional resolution. Our dating analysis provides contrasting divergence times using fossils and Bayesian and non-Bayesian approaches. Our study refines the evolutionary history of rails, offering a foundation for future evolutionary studies of birds.

Paleomicrobiology: Diagnosis and Evolution of Ancient Pathogens

The last century has witnessed progress in the study of ancient infectious disease from purely medical descriptions of past ailments to dynamic interpretations of past population health that draw upon multiple perspectives. The recent adoption of high-throughput DNA sequencing has led to an expanded understanding of pathogen presence, evolution, and ecology across the globe. This genomic revolution has led to the identification of disease-causing microbes in both expected and unexpected contexts, while also providing for the genomic characterization of ancient pathogens previously believed to be unattainable by available methods. In this review we explore the development of DNA-based ancient pathogen research, the specialized methods and tools that have emerged to authenticate and explore infectious disease of the past, and the unique challenges that persist in molecular paleopathology. We offer guidelines to mitigate the impact of these challenges, which will allow for more reliable interpretations of data in this rapidly evolving field of investigation.

East African cichlid lineages (Teleostei: Cichlidae) might be older than their ancient host lakes: new divergence estimates for the east African cichlid radiation

BACKGROUND

Cichlids are a prime model system in evolutionary research and several of the most prominent examples of adaptive radiations are found in the East African Lakes Tanganyika, Malawi and Victoria, all part of the East African cichlid radiation (EAR). In the past, great effort has been invested in reconstructing the evolutionary and biogeographic history of cichlids (Teleostei: Cichlidae). In this study, we present new divergence age estimates for the major cichlid lineages with the main focus on the EAR based on a dataset encompassing representative taxa of almost all recognized cichlid tribes and ten mitochondrial protein genes. We have thoroughly re-evaluated both fossil and geological calibration points, and we included the recently described fossil †Tugenchromis pickfordi in the cichlid divergence age estimates.

RESULTS

Our results estimate the origin of the EAR to Late Eocene/Early Oligocene (28.71 Ma; 95% HPD: 24.43-33.15 Ma). More importantly divergence ages of the most recent common ancestor (MRCA) of several Tanganyika cichlid tribes were estimated to be substantially older than the oldest estimated maximum age of the Lake Tanganyika: Trematocarini (16.13 Ma, 95% HPD: 11.89-20.46 Ma), Bathybatini (20.62 Ma, 95% HPD: 16.88-25.34 Ma), Lamprologini (15.27 Ma; 95% HPD: 12.23-18.49 Ma). The divergence age of the crown haplochromine H-lineage is estimated to 22.8 Ma (95% HPD: 14.40-26.32 Ma) and of the Lake Malawi radiation to 4.07 Ma (95% HDP: 2.93-5.26 Ma). In addition, we recovered a novel lineage within the Lamprologini tribe encompassing only Lamprologus of the lower and central Congo drainage with its divergence estimated to the Late Miocene or early Pliocene. Furthermore we recovered two novel mitochondrial haplotype lineages within the Haplochromini tribe: 'Orthochromis' indermauri and 'Haplochormis' vanheusdeni.

CONCLUSIONS

Divergence time estimates of the MRCA of several Tanganyika cichlid tribes predate the age of the extant Lake Tanganyika basin, and hence are in line with the recently formulated "Melting-Pot Tanganyika" hypothesis. The radiation of the 'Lower Congo Lamprologus clade' might be linked with the Pliocene origin of the modern lower Congo rapids as has been shown for other Lower Congo cichlid assemblages. Finally, the age of origin of the Lake Malawi cichlid flock agrees well with the oldest age estimate for lacustrine conditions in Lake Malawi.

Lizards of the lost arcs: mid-Cenozoic diversification, persistence and ecological marginalization in the West Pacific

Regions with complex geological histories often have diverse and highly endemic biotas, yet inferring the ecological and historical processes shaping this relationship remains challenging. Here, in the context of the taxon cycle model of insular community assembly, we investigate patterns of lineage diversity and habitat usage in a newly characterized vertebrate radiation centred upon the world's most geologically complex insular region: island arcs spanning from the Philippines to Fiji. On island arcs taxa are ecologically widespread, and provide evidence to support one key prediction of the taxon cycle, specifically that interior habitats (lowland rainforests, montane habitats) are home to a greater number of older or relictual lineages than are peripheral habitats (coastal and open forests). On continental fringes, however, the clade shows a disjunct distribution away from lowland rainforest, occurring in coastal, open or montane habitats. These results are consistent with a role for biotic interactions in shaping disjunct distributions (a central tenant of the taxon cycle), but we find this pattern most strongly on continental fringes not islands. Our results also suggest that peripheral habitats on islands, and especially island arcs, may be important for persistence and diversification, not just dispersal and colonization. Finally, new phylogenetic evidence for subaerial island archipelagos (with an associated biota) east of present-day Wallace's Line since the Oligocene has important implications for understanding long-term biotic interchange and assembly across Asia and Australia.

The complex study of complexes: The first well-supported phylogeny of two species complexes within genus Caridina (Decapoda: Caridea: Atyidae) sheds light on evolution, biogeography, and habitat

Atyid shrimps, a key component of tropical freshwater ecosystems, face multiple anthropogenic threats and thus need special attention. With more than 300 described species, the genus Caridina is the most speciose of all the Caridea infra-order. Caridina spp. occupy diverse habitats in tropical freshwaters of the Indo-West Pacific region. Several species complexes have been recognized, based on common morphological features, but little is known about how well these morphological characteristics align with phylogenetic characteristics. Furthermore, no phylogeny of the genus Caridina published so far has provided well-resolved and supported relationships among different species, thus impeding the possibility of proposing evolutionary hypotheses. In this study we used next generation sequencing (NGS) to provide new insights into the phylogenetic relationships among the genus Caridina, focusing on two complexes: 'Caridina nilotica' and 'Caridina weberi'. We collected 92 specimens belonging to these two groups from most of their known geographical range, representing 50 species, for which we sequenced seven mitochondrial genes and two nuclear markers using ion torrent NGS. We performed a phylogenetic analysis, which yielded the first well-supported tree for the genus Caridina. On this tree were mapped the geographic ranges and the habitats used by the different species, and a time calibration was tested. We found the driving factors that most likely account for separation of clades are differences in habitat and to a lesser extent geography. This work provides new insights into the taxonomy of this group and identifies opportunities for further studies in order to fill knowledge gaps that currently impede the management and conservation of atyid species.

Species delimitation of the North American orchard-spider Leucauge venusta (Walckenaer, 1841) (Araneae, Tetragnathidae)

The orchard spider, Leucauge venusta (Walckenaer, 1841) is one of the most common and abundant orb-weavers in North America. This species has a broad geographic distribution extending across tropical and temperate regions of the Americas from Canada to Brazil. Guided by a preliminary observation of the barcode gap between sequences from specimens of L. venusta collected in Florida and other North American localities, we collected across a transect through the southeastern USA to investigate the observed genetic divide. The dataset, complemented with additional samples from Mexico, and Brazil was analyzed for species delimitation using STACEY and bGMYC based on sequences from one nuclear (ITS2) and one mitochondrial marker (COI). The analyses clearly separate USA samples into two deeply divergent and geographically structured groups (north-south) which we interpret as two different species. We generated ecological niche models for these two groups rejecting a niche equivalence hypothesis for these lineages. Taxonomic changes are proposed based on these findings, Leucauge venusta is restricted to denote the northern clade, and its known distribution restricted to the USA. Leucauge argyrobapta (White, 1841) is removed from synonymy to denote the populations in Florida, Mexico and Brazil. Although the delimitation analyses suggest each of these geographic clusters within the L. argyrobapta samples represent different species, more specimens from Central and South America are needed to properly test the cohesion of L. argyrobapta populations.

Molecular Phylogenies indicate a Paleo-Tibetan Origin of Himalayan Lazy Toads (Scutiger)

The Himalaya presents an outstanding geologically active orogen and biodiversity hotspot. However, our understanding of the historical biogeography of its fauna is far from comprehensive. Many taxa are commonly assumed to have originated from China-Indochina and dispersed westward along the Himalayan chain. Alternatively, the “Tibetan-origin hypothesis” suggests primary diversification of lineages in Paleo-Tibet, and secondary diversification along the slopes of the later uplifted Greater Himalaya. We test these hypotheses in high-mountain megophryid anurans (Scutiger). Extensive sampling from High Asia, and analyses of mitochondrial (2839 bp) and nuclear DNA (2208 bp), using Bayesian and Maximum likelihood phylogenetics, suggest that the Himalayan species form a distinct clade, possibly older than those from the eastern Himalaya-Tibet orogen. While immigration from China-Indochina cannot be excluded, our data may indicate that Himalayan Scutiger originated to the north of the Himalaya by colonization from Paleo-Tibet and then date back to the Oligocene. High intraspecific diversity of Scutiger implies limited migration across mountains and drainages along the Himalaya. While our study strengthens support for a “Tibetan-origin hypothesis”, current sampling (10/22 species; 1 revalidated: S. occidentalis) remains insufficient to draw final conclusions on Scutiger but urges comparative phylogeographers to test alternative, geologically supported hypotheses for a true future understanding of Himalayan biogeography.

Phylogenetic evidence that both ancient vicariance and dispersal have contributed to the biogeographic patterns of anchialine cave shrimps

Cave shrimps from the genera Typhlatya, Stygiocaris and Typhlopatsa (Atyidae) are restricted to specialised coastal subterranean habitats or nearby freshwaters and have a highly disconnected distribution (Eastern Pacific, Caribbean, Atlantic, Mediterranean, Madagascar, Australia). The combination of a wide distribution and a limited dispersal potential suggests a large-scale process has generated this geographic pattern. Tectonic plates that fragment ancestral ranges (vicariance) has often been assumed to cause this process, with the biota as passive passengers on continental blocks. The ancestors of these cave shrimps are believed to have inhabited the ancient Tethys Sea, with three particular geological events hypothesised to have led to their isolation and divergence; (1) the opening of the Atlantic Ocean, (2) the breakup of Gondwana, and (3) the closure of the Tethys Seaway. We test the relative contribution of vicariance and dispersal in the evolutionary history of this group using mitochondrial genomes to reconstruct phylogenetic and biogeographic scenarios with fossil-based calibrations. Given that the Australia/Madagascar shrimp divergence postdates the Gondwanan breakup, our results suggest both vicariance (the Atlantic opening) and dispersal. The Tethys closure appears not to have been influential, however we hypothesise that changing marine currents had an important early influence on their biogeography.

Life on the rocks: Multilocus phylogeography of rock hyrax (Procavia capensis) from southern Africa

Understanding the role of geography and climatic cycles in determining patterns of biodiversity is important in comparative and evolutionary biology and conservation. We studied the phylogeographic pattern and historical demography of a rock-dwelling small mammal species from southern Africa, the rock hyrax Procavia capensis capensis. Using a multilocus coalescent approach, we assessed the influence of strong habitat dependence and fluctuating regional climates on genetic diversity. We sequenced a mitochondrial gene (cytochrome b) and two nuclear introns (AP5, PRKC1) supplemented with microsatellite genotyping, in order to assess evolutionary processes over multiple temporal scales. In addition, distribution modelling was used to investigate the current and predicted distribution of the species under different climatic scenarios. Collectively, the data reveal a complex history of isolation followed by secondary contact shaping the current intraspecific diversity. The cyt b sequences confirmed the presence of two previously proposed geographically and genetically distinct lineages distributed across the southern African Great Escarpment and north-western mountain ranges. Molecular dating suggests Miocene divergence of the lineages, yet there are no discernible extrinsic barriers to gene flow. The nuclear markers reveal incomplete lineage sorting or ongoing mixing of the two lineages. Although the microsatellite data lend some support to the presence of two subpopulations, there is weak structuring within and between lineages. These data indicate the presence of gene flow from the northern into the southern parts of the southern African sub-region likely following the secondary contact. The distribution modelling predictably reveal the species' preference for rocky areas, with stable refugia through time in the northern mountain ranges, the Great Escarpment, as well as restricted areas of the Northern Cape Province and the Cape Fold Mountains of South Africa. Different microclimatic variables appear to determine the distributional range of the species. Despite strong habitat preference, the micro-habitat offered by rocky crevices and unique life history traits likely promoted the adaptability of P. capensis, resulting in the widespread distribution and persistence of the species over a long evolutionary period. Spatio-temporal comparison of the evolutionary histories of other co-distributed species across the rocky landscapes of southern Africa will improve our understanding of the regional patterns of biodiversity and local endemism.

Fossils matter: improved estimates of divergence times in Pinus reveal older diversification

BACKGROUND

The taxonomy of pines (genus Pinus) is widely accepted and a robust gene tree based on entire plastome sequences exists. However, there is a large discrepancy in estimated divergence times of major pine clades among existing studies, mainly due to differences in fossil placement and dating methods used. We currently lack a dated molecular phylogeny that makes use of the rich pine fossil record, and this study is the first to estimate the divergence dates of pines based on a large number of fossils (21) evenly distributed across all major clades, in combination with applying both node and tip dating methods.

RESULTS

We present a range of molecular phylogenetic trees of Pinus generated within a Bayesian framework. We find the origin of crown Pinus is likely up to 30 Myr older (Early Cretaceous) than inferred in most previous studies (Late Cretaceous) and propose generally older divergence times for major clades within Pinus than previously thought. Our age estimates vary significantly between the different dating approaches, but the results generally agree on older divergence times. We present a revised list of 21 fossils that are suitable to use in dating or comparative analyses of pines.

CONCLUSIONS

Reliable estimates of divergence times in pines are essential if we are to link diversification processes and functional adaptation of this genus to geological events or to changing climates. In addition to older divergence times in Pinus, our results also indicate that node age estimates in pines depend on dating approaches and the specific fossil sets used, reflecting inherent differences in various dating approaches. The sets of dated phylogenetic trees of pines presented here provide a way to account for uncertainties in age estimations when applying comparative phylogenetic methods.

Inferring Trees

Molecular evolution can reveal the relationship between sets of homologous sequences and the patterns of change that occur during their evolution. An important aspect of these studies is the inference of a phylogenetic tree, which explicitly describes evolutionary relationships between homologous sequences. This chapter provides an introduction to evolutionary trees and how to infer them from sequence data using some commonly used inferential methodology. It focuses on statistical methods for inferring trees and how to assess the confidence one should have in any resulting tree, with a particular emphasis on the underlying assumptions of the methods and how they might affect the tree estimate. There is also some discussion of the underlying algorithms used to perform tree search and recommendations regarding the performance of different algorithms. Finally, there are a few practical guidelines, including how to combine multiple software packages to improve inference, and a comparison between Bayesian and Maximum likelihood phylogenetics.

A Time-Calibrated Mitogenome Phylogeny of Catfish (Teleostei: Siluriformes)

A very significant part of the world's freshwater ichthyofauna is represented by ancient, exceptionally diverse and cosmopolitan ray-finned teleosts of the order Siluriformes. Over the years, catfish have been established as an exemplary model for probing historical biogeography at various scales. Yet, several tantalizing gaps still exist in their phylogenetic history, timeline and mode of diversification. Here, we re-examine the phylogeny of catfish by assembling and analyzing almost all publicly available mitogenome data. We constructed an ingroup matrix of 62 full-length mitogenome sequences from 20 catfish families together with four cypriniform outgroups, spanning 15,557 positions in total. Partitioned maximum likelihood analyses and Bayesian relaxed clock dating using fossil age constraints provide some useful and novel insights into the evolutionary history of this group. Loricarioidei are recovered as the first siluriform group to diversify, rendering Neotropics the cradle of the order. The next deepest clade is the South American Diplomystoidei placed as a sister group to all the remaining Siluroidei. The two multifamilial clades of "Big Asia" and "Big Africa" are also recovered, albeit nodal support for the latter is poor. Within "Big Asia", Bagridae are clearly polyphyletic. Other interfamilial relationships, including Clariidae + Heteropneustidae, Doradidae + Auchenipteridae and Ictaluridae + Cranoglanididae are robustly resolved. Our chronogram shows that siluriforms have a Pangaean origin, at least as far back as the Early Cretaceous. The inferred timeline of the basal splits corroborates the "Out-of-South America" hypothesis and accords well with the fossil record. The divergence of Siluroidei most likely postdated the final separation of Africa and South America. An appealing case of phylogenetic affinity elaborated by biogeographic dispersal is exemplified by the Early Paleogene split between the Southeast Asian Cranoglanididae and Ictaluridae, with the latter radiating into North America's freshwater realm by Eocene. The end of Cretaceous probably concludes the major bout of diversification at the family level while with the dawn of the Cenozoic a prolific radiation is evident at the generic level.

Phylogeny and diversification of mountain vipers (Montivipera, Nilson et al., 2001) triggered by multiple Plio-Pleistocene refugia and high-mountain topography in the Near and Middle East

The Near and Middle East is a hotspot of biodiversity, but the region remains underexplored at the level of genetic biodiversity. Here, we present an extensive molecular phylogeny of the viperid snake genus Montivipera, including all known taxa. Based on nuclear and mitochondrial data, we present novel insights into the phylogeny of the genus and review the status of its constituent species. Maximum likelihood methods revealed a montane origin of Montivipera at 12.3Mya. We then analyzed factors of mountain viper diversity. Our data support substantial changes in effective population size through Plio-Pleistocene periods. We conclude that climatic oscillations were drivers of allopatric speciation, and that mountain systems of the Near and Middle East have strongly influenced the evolution and survival of taxa, because climatic and topographical heterogeneities induced by mountains have played a crucial role as filters for dispersal and as multiple refugia. The wide diversity of montane microhabitats enabled mountain vipers to retain their ecological niche during climatic pessima. In consequence the varied geological and topographical conditions between refugia favoured genetic isolation and created patterns of species richness resulting in the formation of neoendemic taxa. Our data support high concordance between geographic distributions of Montivipera haplotypes with putative plant refugia.

Diversification rates and phenotypic evolution in venomous snakes (Elapidae)

The relationship between rates of diversification and of body size change (a common proxy for phenotypic evolution) was investigated across Elapidae, the largest radiation of highly venomous snakes. Time-calibrated phylogenetic trees for 175 species of elapids (more than 50% of known taxa) were constructed using seven mitochondrial and nuclear genes. Analyses using these trees revealed no evidence for a link between speciation rates and changes in body size. Two clades (Hydrophis, Micrurus) show anomalously high rates of diversification within Elapidae, yet exhibit rates of body size evolution almost identical to the general elapid 'background' rate. Although correlations between speciation rates and rates of body size change exist in certain groups (e.g. ray-finned fishes, passerine birds), the two processes appear to be uncoupled in elapid snakes. There is also no detectable shift in diversification dynamics associated with the colonization of Australasia, which is surprising given that elapids appear to be the first clade of venomous snakes to reach the continent.

Danger in the reef: Proteome, toxicity, and neutralization of the venom of the olive sea snake, Aipysurus laevis

Four specimens of the olive sea snake, Aipysurus laevis, were collected off the coast of Western Australia, and the venom proteome was characterized and quantitatively estimated by RP-HPLC, SDS-PAGE, and MALDI-TOF-TOF analyses. A. laevis venom is remarkably simple and consists of phospholipases A2 (71.2%), three-finger toxins (3FTx; 25.3%), cysteine-rich secretory proteins (CRISP; 2.5%), and traces of a complement control module protein (CCM; 0.2%). Using a Toxicity Score, the most lethal components were determined to be short neurotoxins. Whole venom had an intravenous LD50 of 0.07 mg/kg in mice and showed a high phospholipase A2 activity, but no proteinase activity in vitro. Preclinical assessment of neutralization and ELISA immunoprofiling showed that BioCSL Sea Snake Antivenom was effective in cross-neutralizing A. laevis venom with an ED50 of 821 μg venom per mL antivenom, with a binding preference towards short neurotoxins, due to the high degree of conservation between short neurotoxins from A. laevis and Enhydrina schistosa venom. Our results point towards the possibility of developing recombinant antibodies or synthetic inhibitors against A. laevis venom due to its simplicity.

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.

Alectorioid Morphologies in Paleogene Lichens: New Evidence and Re-Evaluation of the Fossil Alectoria succini Mägdefrau

One of the most important issues in molecular dating studies concerns the incorporation of reliable fossil taxa into the phylogenies reconstructed from DNA sequence variation in extant taxa. Lichens are symbiotic associations between fungi and algae and/or cyanobacteria. Several lichen fossils have been used as minimum age constraints in recent studies concerning the diversification of the Ascomycota. Recent evolutionary studies of Lecanoromycetes, an almost exclusively lichen-forming class in the Ascomycota, have utilized the Eocene amber inclusion Alectoria succinic as a minimum age constraint. However, a re-investigation of the type material revealed that this inclusion in fact represents poorly preserved plant remains, most probably of a root. Consequently, this fossil cannot be used as evidence of the presence of the genus Alectoria (Parmeliaceae, Lecanorales) or any other lichens in the Paleogene. However, newly discovered inclusions from Paleogene Baltic and Bitterfeld amber verify that alectorioid morphologies in lichens were in existence by the Paleogene. The new fossils represent either a lineage within the alectorioid group or belong to the genus Oropogon.

Examining phylogenetic relationships among gibbon genera using whole genome sequence data using an approximate bayesian computation approach

Gibbons are believed to have diverged from the larger great apes ∼16.8 MYA and today reside in the rainforests of Southeast Asia. Based on their diploid chromosome number, the family Hylobatidae is divided into four genera, Nomascus, Symphalangus, Hoolock, and Hylobates. Genetic studies attempting to elucidate the phylogenetic relationships among gibbons using karyotypes, mitochondrial DNA (mtDNA), the Y chromosome, and short autosomal sequences have been inconclusive . To examine the relationships among gibbon genera in more depth, we performed second-generation whole genome sequencing (WGS) to a mean of ∼15× coverage in two individuals from each genus. We developed a coalescent-based approximate Bayesian computation (ABC) method incorporating a model of sequencing error generated by high coverage exome validation to infer the branching order, divergence times, and effective population sizes of gibbon taxa. Although Hoolock and Symphalangus are likely sister taxa, we could not confidently resolve a single bifurcating tree despite the large amount of data analyzed. Instead, our results support the hypothesis that all four gibbon genera diverged at approximately the same time. Assuming an autosomal mutation rate of 1 × 10(-9)/site/year this speciation process occurred ∼5 MYA during a period in the Early Pliocene characterized by climatic shifts and fragmentation of the Sunda shelf forests. Whole genome sequencing of additional individuals will be vital for inferring the extent of gene flow among species after the separation of the gibbon genera.

A species-level phylogeny of all extant and late Quaternary extinct mammals using a novel heuristic-hierarchical Bayesian approach

Across large clades, two problems are generally encountered in the estimation of species-level phylogenies: (a) the number of taxa involved is generally so high that computation-intensive approaches cannot readily be utilized and (b) even for clades that have received intense study (e.g., mammals), attention has been centered on relatively few selected species, and most taxa must therefore be positioned on the basis of very limited genetic data. Here, we describe a new heuristic-hierarchical Bayesian approach and use it to construct a species-level phylogeny for all extant and late Quaternary extinct mammals. In this approach, species with large quantities of genetic data are placed nearly freely in the mammalian phylogeny according to these data, whereas the placement of species with lower quantities of data is performed with steadily stricter restrictions for decreasing data quantities. The advantages of the proposed method include (a) an improved ability to incorporate phylogenetic uncertainty in downstream analyses based on the resulting phylogeny, (b) a reduced potential for long-branch attraction or other types of errors that place low-data taxa far from their true position, while maintaining minimal restrictions for better-studied taxa, and (c) likely improved placement of low-data taxa due to the use of closer outgroups.

An evaluation of fossil tip-dating versus node-age calibrations in tetraodontiform fishes (Teleostei: Percomorphaceae)

Time-calibrated phylogenies based on molecular data provide a framework for comparative studies. Calibration methods to combine fossil information with molecular phylogenies are, however, under active development, often generating disagreement about the best way to incorporate paleontological data into these analyses. This study provides an empirical comparison of the most widely used approach based on node-dating priors for relaxed clocks implemented in the programs BEAST and MrBayes, with two recently proposed improvements: one using a new fossilized birth-death process model for node dating (implemented in the program DPPDiv), and the other using a total-evidence or tip-dating method (implemented in MrBayes and BEAST). These methods are applied herein to tetraodontiform fishes, a diverse group of living and extinct taxa that features one of the most extensive fossil records among teleosts. Previous estimates of time-calibrated phylogenies of tetraodontiforms using node-dating methods reported disparate estimates for their age of origin, ranging from the late Jurassic to the early Paleocene (ca. 150-59Ma). We analyzed a comprehensive dataset with 16 loci and 210 morphological characters, including 131 taxa (95 extant and 36 fossil species) representing all families of fossil and extant tetraodontiforms, under different molecular clock calibration approaches. Results from node-dating methods produced consistently younger ages than the tip-dating approaches. The older ages inferred by tip dating imply an unlikely early-late Jurassic (ca. 185-119Ma) origin for this order and the existence of extended ghost lineages in their fossil record. Node-based methods, by contrast, produce time estimates that are more consistent with the stratigraphic record, suggesting a late Cretaceous (ca. 86-96Ma) origin. We show that the precision of clade age estimates using tip dating increases with the number of fossils analyzed and with the proximity of fossil taxa to the node under assessment. This study suggests that current implementations of tip dating may overestimate ages of divergence in calibrated phylogenies. It also provides a comprehensive phylogenetic framework for tetraodontiform systematics and future comparative studies.

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.

Methods for the quantitative comparison of molecular estimates of clade age and the fossil record

Approaches quantifying the relative congruence, or incongruence, of molecular divergence estimates and the fossil record have been limited. Previously proposed methods are largely node specific, assessing incongruence at particular nodes for which both fossil data and molecular divergence estimates are available. These existing metrics, and other methods that quantify incongruence across topologies including entirely extinct clades, have so far not taken into account uncertainty surrounding both the divergence estimates and the ages of fossils. They have also treated molecular divergence estimates younger than previously assessed fossil minimum estimates of clade age as if they were the same as cases in which they were older. However, these cases are not the same. Recovered divergence dates younger than compared oldest known occurrences require prior hypotheses regarding the phylogenetic position of the compared fossil record and standard assumptions about the relative timing of morphological and molecular change to be incorrect. Older molecular dates, by contrast, are consistent with an incomplete fossil record and do not require prior assessments of the fossil record to be unreliable in some way. Here, we compare previous approaches and introduce two new descriptive metrics. Both metrics explicitly incorporate information on uncertainty by utilizing the 95% confidence intervals on estimated divergence dates and data on stratigraphic uncertainty concerning the age of the compared fossils. Metric scores are maximized when these ranges are overlapping. MDI (minimum divergence incongruence) discriminates between situations where molecular estimates are younger or older than known fossils reporting both absolute fit values and a number score for incompatible nodes. DIG range (divergence implied gap range) allows quantification of the minimum increase in implied missing fossil record induced by enforcing a given set of molecular-based estimates. These metrics are used together to describe the relationship between time trees and a set of fossil data, which we recommend be phylogenetically vetted and referred on the basis of apomorphy. Differences from previously proposed metrics and the utility of MDI and DIG range are illustrated in three empirical case studies from angiosperms, ostracods, and birds. These case studies also illustrate the ways in which MDI and DIG range may be used to assess time trees resultant from analyses varying in calibration regime, divergence dating approach or molecular sequence data analyzed.

The impact of calibration and clock-model choice on molecular estimates of divergence times

Phylogenetic estimates of evolutionary timescales can be obtained from nucleotide sequence data using the molecular clock. These estimates are important for our understanding of evolutionary processes across all taxonomic levels. The molecular clock needs to be calibrated with an independent source of information, such as fossil evidence, to allow absolute ages to be inferred. Calibration typically involves fixing or constraining the age of at least one node in the phylogeny, enabling the ages of the remaining nodes to be estimated. We conducted an extensive simulation study to investigate the effects of the position and number of calibrations on the resulting estimate of the timescale. Our analyses focused on Bayesian estimates obtained using relaxed molecular clocks. Our findings suggest that an effective strategy is to include multiple calibrations and to prefer those that are close to the root of the phylogeny. Under these conditions, we found that evolutionary timescales could be estimated accurately even when the relaxed-clock model was misspecified and when the sequence data were relatively uninformative. We tested these findings in a case study of simian foamy virus, where we found that shallow calibrations caused the overall timescale to be underestimated by up to three orders of magnitude. Finally, we provide some recommendations for improving the practice of molecular-clock calibration.

Assignment of Calibration Information to Deeper Phylogenetic Nodes is More Effective in Obtaining Precise and Accurate Divergence Time Estimates

Divergence time estimation has become an essential tool for understanding macroevolutionary events. Molecular dating aims to obtain reliable inferences, which, within a statistical framework, means jointly increasing the accuracy and precision of estimates. Bayesian dating methods exhibit the propriety of a linear relationship between uncertainty and estimated divergence dates. This relationship occurs even if the number of sites approaches infinity and places a limit on the maximum precision of node ages. However, how the placement of calibration information may affect the precision of divergence time estimates remains an open question. In this study, relying on simulated and empirical data, we investigated how the location of calibration within a phylogeny affects the accuracy and precision of time estimates. We found that calibration priors set at median and deep phylogenetic nodes were associated with higher precision values compared to analyses involving calibration at the shallowest node. The results were independent of the tree symmetry. An empirical mammalian dataset produced results that were consistent with those generated by the simulated sequences. Assigning time information to the deeper nodes of a tree is crucial to guarantee the accuracy and precision of divergence times. This finding highlights the importance of the appropriate choice of outgroups in molecular dating.

Estimating the Phanerozoic history of the Ascomycota lineages: combining fossil and molecular data

The phylum Ascomycota is by far the largest group in the fungal kingdom. Ecologically important mutualistic associations such as mycorrhizae and lichens have evolved in this group, which are regarded as key innovations that supported the evolution of land plants. Only a few attempts have been made to date the origin of Ascomycota lineages by using molecular clock methods, which is primarily due to the lack of satisfactory fossil calibration data. For this reason we have evaluated all of the oldest available ascomycete fossils from amber (Albian to Miocene) and chert (Devonian and Maastrichtian). The fossils represent five major ascomycete classes (Coniocybomycetes, Dothideomycetes, Eurotiomycetes, Laboulbeniomycetes, and Lecanoromycetes). We have assembled a multi-gene data set (18SrDNA, 28SrDNA, RPB1 and RPB2) from a total of 145 taxa representing most groups of the Ascomycota and utilized fossil calibration points solely from within the ascomycetes to estimate divergence times of Ascomycota lineages with a Bayesian approach. Our results suggest an initial diversification of the Pezizomycotina in the Ordovician, followed by repeated splits of lineages throughout the Phanerozoic, and indicate that this continuous diversification was unaffected by mass extinctions. We suggest that the ecological diversity within each lineage ensured that at least some taxa of each group were able to survive global crises and rapidly recovered.

The emergence of lobsters: phylogenetic relationships, morphological evolution and divergence time comparisons of an ancient group (decapoda: achelata, astacidea, glypheidea, polychelida)

Lobsters are a ubiquitous and economically important group of decapod crustaceans that include the infraorders Polychelida, Glypheidea, Astacidea and Achelata. They include familiar forms such as the spiny, slipper, clawed lobsters and crayfish and unfamiliar forms such as the deep-sea and "living fossil" species. The high degree of morphological diversity among these infraorders has led to a dynamic classification and conflicting hypotheses of evolutionary relationships. In this study, we estimated phylogenetic relationships among the major groups of all lobster families and 94% of the genera using six genes (mitochondrial and nuclear) and 195 morphological characters across 173 species of lobsters for the most comprehensive sampling to date. Lobsters were recovered as a non-monophyletic assemblage in the combined (molecular + morphology) analysis. All families were monophyletic, with the exception of Cambaridae, and 7 of 79 genera were recovered as poly- or paraphyletic. A rich fossil history coupled with dense taxon coverage allowed us to estimate and compare divergence times and origins of major lineages using two drastically different approaches. Age priors were constructed and/or included based on fossil age information or fossil discovery, age, and extant species count data. Results from the two approaches were largely congruent across deep to shallow taxonomic divergences across major lineages. The origin of the first lobster-like decapod (Polychelida) was estimated in the Devonian (∼409-372 Ma) with all infraorders present in the Carboniferous (∼353-318 Ma). Fossil calibration subsampling studies examined the influence of sampling density (number of fossils) and placement (deep, middle, and shallow) on divergence time estimates. Results from our study suggest including at least 1 fossil per 10 operational taxonomic units (OTUs) in divergence dating analyses. [Dating; decapods; divergence; lobsters; molecular; morphology; phylogenetics.].

Darwinian shortfalls in biodiversity conservation

If we were to describe all the species on Earth and determine their distributions, we would solve the popularly termed 'Linnean' and 'Wallacean' shortfalls in biodiversity conservation. Even so, we would still be hindered by a 'Darwinian shortfall', that is, the lack of relevant phylogenetic information for most organisms. Overall, there are too few comprehensive phylogenies, large uncertainties in the estimation of divergence times, and, most critically, unknown evolutionary models linking phylogenies to relevant ecological traits and life history variation. Here, we discuss these issues and offer suggestions for further research to support evolutionary-based conservation planning.

From β- to α-proteobacteria: the origin and evolution of rhizobial nodulation genes nodIJ

Although many α- and some β-proteobacterial species are symbiotic with legumes, the evolutionary origin of nitrogen-fixing nodulation remains unclear. We examined α- and β-proteobacteria whose genomes were sequenced using large-scale phylogenetic profiling and revealed the evolutionary origin of two nodulation genes. These genes, nodI and nodJ (nodIJ), play key roles in the secretion of Nod factors, which are recognized by legumes during nodulation. We found that only the nodulating β-proteobacteria, including the novel strains isolated in this study, possess both nodIJ and their paralogous genes (DRA-ATPase/permease genes). Contrary to the widely accepted scenario of the α-proteobacterial origin of rhizobia, our exhaustive phylogenetic analysis showed that the entire nodIJ clade is included in the clade of Burkholderiaceae DRA-ATPase/permease genes, that is, the nodIJ genes originated from gene duplication in a lineage of the β-proteobacterial family. After duplication, the evolutionary rates of nodIJ were significantly accelerated relative to those of homologous genes, which is consistent with their novel function in nodulation. The likelihood analyses suggest that this accelerated evolution is not associated with changes in either nonsynonymous/synonymous substitution rates or transition/transversion rates, but rather, in the GC content. Although the low GC content of the nodulation genes has been assumed to reflect past horizontal transfer events from donor rhizobial genomes with low GC content, no rhizobial genome with such low GC content has yet been found. Our results encourage a reconsideration of the origin of nodulation and suggest new perspectives on the role of the GC content of bacterial genes in functional adaptation.

A multi-calibrated mitochondrial phylogeny of extant Bovidae (Artiodactyla, Ruminantia) and the importance of the fossil record to systematics

BACKGROUND

Molecular phylogenetics has provided unprecedented resolution in the ruminant evolutionary tree. However, molecular age estimates using only one or a few (often misapplied) fossil calibration points have produced a diversity of conflicting ages for important evolutionary events within this clade. I here identify 16 fossil calibration points of relevance to the phylogeny of Bovidae and Ruminantia and use these, individually and together, to construct a dated molecular phylogeny through a reanalysis of the full mitochondrial genome of over 100 ruminant species.

RESULTS

The new multi-calibrated tree provides ages that are younger overall than found in previous studies. Among these are young ages for the origin of crown Ruminantia (39.3-28.8 Ma), and crown Bovidae (17.3-15.1 Ma). These are argued to be reasonable hypotheses given that many basal fossils assigned to these taxa may in fact lie on the stem groups leading to the crown clades, thus inflating previous age estimates. Areas of conflict between molecular and fossil dates do persist, however, especially with regard to the base of the rapid Pecoran radiation and the sister relationship of Moschidae to Bovidae. Results of the single-calibrated analyses also show that a very wide range of molecular age estimates are obtainable using different calibration points, and that the choice of calibration point can influence the topology of the resulting tree. Compared to the single-calibrated trees, the multi-calibrated tree exhibits smaller variance in estimated ages and better reflects the fossil record.

CONCLUSIONS

The use of a large number of vetted fossil calibration points with soft bounds is promoted as a better approach than using just one or a few calibrations, or relying on internal-congruency metrics to discard good fossil data. This study also highlights the importance of considering morphological and ecological characteristics of clades when delimiting higher taxa. I also illustrate how phylogeographic and paleoenvironmental hypotheses inferred from a tree containing only extant taxa can be problematic without consideration of the fossil record. Incorporating the fossil record of Ruminantia is a necessary step for future analyses aiming to reconstruct the evolutionary history of this clade.

A multi-locus molecular timescale for the origin and diversification of eels (Order: Anguilliformes)

Anguilliformes are an ecologically diverse group of predominantly marine fishes whose members are easily recognized by their extremely elongate bodies, and universal lack of pelvic fins. Recent studies based on mitochondrial loci, including full mitogenomes, have called into question the monophyly of both the Anguilliformes, which appear to be paraphyletic without the inclusion of the Saccopharyngiformes (gulper eels and allies), as well as other more commonly known eel families (e.g., Congridae, Serrivomeridae). However, no study to date has investigated anguilliform interrelationships using nuclear loci. Here we present a new phylogenetic hypothesis for the Anguilliformes based on five markers (the nuclear loci Early Growth Hormone 3, Myosin Heavy Polypeptide 6 and Recombinase Activating Gene 1, as well as the mitochondrial genes Cytochrome b and Cytochrome Oxidase I). Our sampling spans 148 species and includes 19 of the 20 extant families of anguilliforms and saccopharyngiforms. Maximum likelihood analysis reveals that saccopharyngiform eels are deeply nested within the anguilliforms, and supports the non-monophyly of Congridae and Nettastomatidae, as well as that of Derichthyidae and Chlopsidae. Our analyses suggest that Protanguilla may be the sister group of the Synaphobranchidae, though the recent hypothesis that this species is the sister group to all other anguilliforms cannot be rejected. The molecular phylogeny, time-calibrated using a Bayesian relaxed clock approach and seven fossil calibration points, reveals a Late Cretaceous origin of this expanded anguilliform clade (stem age ~116 Ma, crown age ~99 Ma). Most major (family level) lineages originated between the end of the Cretaceous and Early Eocene, suggesting that anguilliform radiation may have been facilitated by the recovery of marine ecosystems following the KP extinction.