Biogeographic calibrations for the molecular clock

Phylogeny and biogeography of harmochirine jumping spiders (Araneae: Salticidae)

We use ultraconserved elements (UCE) and Sanger data to study the phylogeny, age, and biogeographical history of harmochirine jumping spiders, a group that includes the species-rich genus Habronattus, whose remarkable courtship has made it the focus of studies of behaviour, sexual selection, and diversification. We recovered 1947 UCE loci from 43 harmochirine taxa and 4 outgroups, yielding a core dataset of 193 UCEs with at least 50 % occupancy. Concatenated likelihood and ASTRAL analyses confirmed the separation of harmochirines into two major clades, here designated the infratribes Harmochirita and Pellenita. Most are African or Eurasian with the notable exception of a clade of pellenites containing Habronattus and Pellenattus of the Americas and Havaika and Hivanua of the Pacific Islands. Biogeographical analysis using the DEC model favours a dispersal of the clade's ancestor from Eurasia to the Americas, from which Havaika's ancestor dispersed to Hawaii and Hivanua's ancestor to the Marquesas Islands. Divergence time analysis on 32 loci with 85 % occupancy, calibrated by fossils and island age, dates the dispersal to the Americas at approximately 4 to 6 million years ago. The explosive radiation of Habronattus perhaps began only about 4 mya. The phylogeny clarifies both the evolution of sexual traits (e.g., the terminal apophyses was enlarged in Pellenes and not subsequently lost) and the taxonomy. Habronattus is confirmed as monophyletic. Pellenattus is raised to the status of genus, and 13 species moved into it as new combinations. Bianor stepposus Logunov, 1991 is transferred to Sibianor, and Pellenes bulawayoensis Wesołowska, 1999 is transferred to Neaetha. A molecular clock rate estimate for spider UCEs is presented and its utility to inform prior distributions is discussed.

Fossil-calibrated molecular clock data enable reconstruction of steps leading to differentiated multicellularity and anisogamy in the Volvocine algae

BACKGROUND

Throughout its nearly four-billion-year history, life has undergone evolutionary transitions in which simpler subunits have become integrated to form a more complex whole. Many of these transitions opened the door to innovations that resulted in increased biodiversity and/or organismal efficiency. The evolution of multicellularity from unicellular forms represents one such transition, one that paved the way for cellular differentiation, including differentiation of male and female gametes. A useful model for studying the evolution of multicellularity and cellular differentiation is the volvocine algae, a clade of freshwater green algae whose members range from unicellular to colonial, from undifferentiated to completely differentiated, and whose gamete types can be isogamous, anisogamous, or oogamous. To better understand how multicellularity, differentiation, and gametes evolved in this group, we used comparative genomics and fossil data to establish a geologically calibrated roadmap of when these innovations occurred.

RESULTS

Our ancestral-state reconstructions, show that multicellularity arose independently twice in the volvocine algae. Our chronograms indicate multicellularity evolved during the Carboniferous-Triassic periods in Goniaceae + Volvocaceae, and possibly as early as the Cretaceous in Tetrabaenaceae. Using divergence time estimates we inferred when, and in what order, specific developmental changes occurred that led to differentiated multicellularity and oogamy. We find that in the volvocine algae the temporal sequence of developmental changes leading to differentiated multicellularity is much as proposed by David Kirk, and that multicellularity is correlated with the acquisition of anisogamy and oogamy. Lastly, morphological, molecular, and divergence time data suggest the possibility of cryptic species in Tetrabaenaceae.

CONCLUSIONS

Large molecular datasets and robust phylogenetic methods are bringing the evolutionary history of the volvocine algae more sharply into focus. Mounting evidence suggests that extant species in this group are the result of two independent origins of multicellularity and multiple independent origins of cell differentiation. Also, the origin of the Tetrabaenaceae-Goniaceae-Volvocaceae clade may be much older than previously thought. Finally, the possibility of cryptic species in the Tetrabaenaceae provides an exciting opportunity to study the recent divergence of lineages adapted to live in very different thermal environments.

Validating a molecular clock for nudibranchs-No fossils to the rescue

Time calibrated phylogenies are typically reconstructed with fossil information but for soft-bodied marine invertebrates that lack hard parts, a fossil record is lacking. In these cases, biogeographic calibrations or the rates of divergence for related taxa are often used. Although nudibranch phylogenies have advanced with the input of molecular data, no study has derived a divergence rate for this diverse group of invertebrates. Here, we use an updated closure date for the Isthmus of Panama (2.8 Ma) to derive the first divergence rates for chromodorid nudibranchs using multigene data from a geminate pair with broad phylogeographic sampling. Examining the species Chromolaichma sedna (Marcus & Marcus, 1967), we uncover deep divergences among eastern Pacific and western Atlantic clades and we erect a new species designation for the latter (Chromolaichma hemera sp. nov.). Next, we discover extensive phylogeographic structure within C. hemera sp. nov. sensu lato, thereby refuting the hypothesis of a recent introduction. Lastly, we derive divergence rates for mitochondrial and nuclear loci that exceed known rates for other gastropods and we highlight significant rate heterogeneity both among markers and taxa. Together, these findings improve understanding of nudibranch systematics and provide rates useful to apply to divergence scenarios in this diverse group.

Parallel signatures of Mycobacterium tuberculosis and human Y-chromosome phylogeography support the Two Layer model of East Asian population history

The Two Layer hypothesis is fast becoming the favoured narrative describing East Asian population history. Under this model, hunter-gatherer groups who initially peopled East Asia via a route south of the Himalayas were assimilated by agriculturalist migrants who arrived via a northern route across Eurasia. A lack of ancient samples from tropical East Asia limits the resolution of this model. We consider insight afforded by patterns of variation within the human pathogen Mycobacterium tuberculosis (Mtb) by analysing its phylogeographic signatures jointly with the human Y-chromosome. We demonstrate the Y-chromosome lineages enriched in the traditionally hunter-gatherer groups associated with East Asia's first layer of peopling to display deep roots, low long-term effective population size, and diversity patterns consistent with a southern entry route. These characteristics mirror those of the evolutionarily ancient Mtb lineage 1. The remaining East Asian Y-chromosome lineage is almost entirely absent from traditionally hunter-gatherer groups and displays spatial and temporal characteristics which are incompatible with a southern entry route, and which link it to the development of agriculture in modern-day China. These characteristics mirror those of the evolutionarily modern Mtb lineage 2. This model paves the way for novel host-pathogen coevolutionary research hypotheses in East Asia.

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.

Origins, phylogenetic relationships and host-parasite interactions of Troglotrematoidea since the cretaceous

In the current study, we raise the issue concerning origins and historical relationships of the trematodes from the families Troglotrematidae and Paragonimidae using phylogenetic analysis and molecular-clock method for estimating evolutionary rates. For the first time we provided 28S rRNA gene fragment (1764 bp) for the type species Troglotrema acutum - zoonotic trematodes that cause cranial lesions (troglotremiasis) in mustelid and canid mammals of the Central Europe, Iberian Peninsula, and North-West Caucasus. Molecular genetic analysis revealed that T. acutum belongs to the monophyletic family Troglotrematidae sister with the family Paragonimidae. The family Troglotrematidae includes five genera: Nanophyetus, Troglotrema, Skrjabinophyetus, Nephrotrema, and Macroorchis; and the family Paragonimidae is monotypic including the only genus Paragonimus. We recover the superfamily Troglotrematoidea for these two families. Divergence of the common ancestor of the superfamily Troglotrematoidea (common troglotrematoid ancestor) likely occurred during the Cretaceous period of the Mesozoic Era and potentially originated in the Asiatic region. The lineage of the family Troglotrematidae is much closer to the common troglotrematoid ancestor than the species of the family Paragonimidae. The radiation time of the common troglotrematoid ancestor (126 Ma, the Early Cretaceous), and formation of the families Troglotrematidae and Paragonimidae (96 Ma and 73 Ma respectively, the Late Cretaceous) corresponds to the time of settling in East Asia by many species of mammaliaforms (about 130-70 Ma).

Variation in the substitution rates among the human mitochondrial haplogroup U sublineages

Resolving the absolute timescale of phylogenetic trees stipulates reliable estimates for the rate of DNA sequence evolution. For this end, various calibration methods have been developed and studied intensively. Intraspecific rate variation among distinct genetic lineages, however, has gained less attention. Here, we have assessed lineage-specific molecular rates of human mitochondrial DNA (mtDNA) by performing tip-calibrated Bayesian phylogenetic analyses. Tip-calibration, as opposed to traditional nodal time stamps from dated fossil evidence or geological events, is based on sample ages and becoming ever more feasible as ancient DNA data from radiocarbon-dated samples accumulate. We focus on subhaplogroups U2, U4, U5a, and U5b, the data including ancient mtDNA genomes from ¹⁴C-dated samples (n = 234), contemporary genomes (n = 301), and two outgroup sequences from haplogroup R. The obtained molecular rates depended on the data sets (with or without contemporary sequences), suggesting time-dependency. More notable was the rate variation between haplogroups: U4 and U5a stand out having a substantially higher rate than U5b. This is also reflected in the divergence times obtained (U5a: 17,700 years and U5b: 29,700 years), a disparity not reported previously. After ruling out various alternative causes (e.g., selection, sampling, and sequence quality), we propose that the substitution rates have been influenced by demographic histories, widely different among populations where U4/U5a or U5b are frequent. As with the Y-chromosomal subhaplogroup R1b, the mitochondrial U4 and U5a have been associated with remarkable range extensions of the Yamnaya culture in the Bronze Age.

An Update of Orthopoxvirus Molecular Evolution

Although variola virus (VARV) has been eradicated through widespread vaccination, other orthopoxviruses pathogenic for humans circulate in nature. Recently, new orthopoxviruses, including some able to infect humans, have been found and their complete genomes have been sequenced. Questions about the orthopoxvirus mutation rate and the emergence of new threats to humankind as a result of the evolution of circulating orthopoxviruses remain open. Based on contemporary data on ancient VARV DNA and DNA of new orthopoxvirus species, an analysis of the molecular evolution of orthopoxviruses was carried out and the timescale of their emergence was estimated. It was calculated that the orthopoxviruses of the Old and New Worlds separated approximately 40,000 years ago; the recently discovered Akhmeta virus and Alaskapox virus separated from other orthopoxviruses approximately 10,000–20,000 years ago; the rest of modern orthopoxvirus species originated from 1700 to 6000 years ago, with the exception of VARV, which emerged in approximately 300 AD. Later, there was a separation of genetic variants of some orthopoxvirus species, so the monkeypox virus West African subtype originated approximately 600 years ago, and the VARV minor alastrim subtype emerged approximately 300 years ago.

Not that young: combining plastid phylogenomic, plate tectonic and fossil evidence indicates a Palaeogene diversification of Cycadaceae

BACKGROUND AND AIMS

Previous molecular dating studies revealed historical mass extinctions and recent radiations of extant cycads, but debates still exist between palaeobotanists and evolutionary biologists regarding the origin and evolution of Cycadaceae.

METHODS

Using whole plastomic data, we revisited the phylogeny of this family and found the Palawan endemic Cycas clade was strongly related to all lineages from Southeast Eurasia, coinciding with a plate drift event occurring in the Early Oligocene. By integrating fossil and biogeographical calibrations as well as molecular data from protein-coding genes, we established different calibration schemes and tested competing evolutionary timelines of Cycadaceae.

KEY RESULTS

We found recent dispersal cannot explain the distribution of Palawan Cycas, yet the scenario including the tectonic calibration yielded a mean crown age of extant Cycadaceae of ~69-43 million years ago by different tree priors, consistent with multiple Palaeogene fossils assigned to this family. Biogeographical analyses incorporating fossil distributions revealed East Asia as the ancestral area of Cycadaceae.

CONCLUSIONS

Our findings challenge the previously proposed Middle-Late Miocene diversification of cycads and an Indochina origin for Cycadaceae and highlight the importance of combining phylogenetic clades, tectonic events and fossils for rebuilding the evolutionary history of lineages that have undergone massive extinctions.

Coastal dunefields maintain pre-Holocene genetic structure in a rocky shore red alga

Most intertidal algae have limited dispersal potential, and areas that lack hard substratum suitable for attachment are thus expected to isolate regional populations from each other. Here, we used nuclear and mitochondrial genetic data to compare genetic structure in two co-distributed intertidal red algae with different dispersal potential along the South African coastline. Gelidium pristoides is divided into a south-eastern and a south-western evolutionary lineage separated by extensive, continuous sandy shoreline habitat adjacent to coastal dunefields. In contrast, Hypnea spicifera is genetically homogeneous throughout its range. In G. pristoides, the genetic breaks are associated with contemporary coastal dunefields. The age of the divergence event suggests that this may reflect the effect of older dispersal barriers, and that genetic structure was subsequently maintained by the formation of contemporary coastal dunefields.

Base-substitution mutation rate across the nuclear genome of Alpheus snapping shrimp and the timing of isolation by the Isthmus of Panama

Background

The formation of the Isthmus of Panama and final closure of the Central American Seaway (CAS) provides an independent calibration point for examining the rate of DNA substitutions. This vicariant event has been widely used to estimate the substitution rate across mitochondrial genomes and to date evolutionary events in other taxonomic groups. Nuclear sequence data is increasingly being used to complement mitochondrial datasets for phylogenetic and evolutionary investigations; these studies would benefit from information regarding the rate and pattern of DNA substitutions derived from the nuclear genome.

Results

To estimate the genome-wide neutral mutation rate (µ), genotype-by-sequencing (GBS) datasets were generated for three transisthmian species pairs in Alpheus snapping shrimp. A range of bioinformatic filtering parameters were evaluated in order to minimize potential bias in mutation rate estimates that may result from SNP filtering. Using a Bayesian coalescent approach (G-PhoCS) applied to 44,960 GBS loci, we estimated µ to be 2.64E−9 substitutions/site/year, when calibrated with the closure of the CAS at 3 Ma. Post-divergence gene flow was detected in one species pair. Failure to account for this post-split migration inflates our substitution rate estimates, emphasizing the importance of demographic methods that can accommodate gene flow.

Conclusions

Results from our study, both parameter estimates and bioinformatic explorations, have broad-ranging implications for phylogeographic studies in other non-model taxa using reduced representation datasets. Our best estimate of µ that accounts for coalescent and demographic processes is remarkably similar to experimentally derived mutation rates in model arthropod systems. These results contradicted recent suggestions that the closure of the Isthmus was completed much earlier (around 10 Ma), as mutation rates based on an early calibration resulted in uncharacteristically low genomic mutation rates. Also, stricter filtering parameters resulted in biased datasets that generated lower mutation rate estimates and influenced demographic parameters, serving as a cautionary tale for the adherence to conservative bioinformatic strategies when generating reduced-representation datasets at the species level. To our knowledge this is the first use of transisthmian species pairs to calibrate the rate of molecular evolution from GBS data.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01836-3.

Fossil-calibrated time tree of Podarcis wall lizards provides limited support for biogeographic calibration models

Podarcis wall lizards are endemic to the Mediterranean Basin where they represent the predominant reptile group. Despite being extensively used as model organisms in evolutionary and ecological studies their phylogeny and historical biogeography are still incompletely resolved. Moreover, molecular clock calibrations used in wall lizard phylogeography are based on the assumption of vicariant speciation triggered by the abrupt Mediterranean Sea level rise at the end of the Messinian salinity crisis (MSC). However, the validity of this biogeographic calibration remains untested. In this study we inferred a robust time tree based on multilocus data and fossil calibrations using both gene concatenation and species-tree approaches and including models with gene-flow. We found five deeply divergent, geographically coherent, and well-supported clades comprising species from i) Iberian Peninsula and North Africa; ii) Western Mediterranean islands, iii) Sicilian and Maltese islands; and iv-v) Balkan region and Aegean islands. The mitochondrial tree shows some inconsistencies with the species tree that warrant future investigation. Diversification of main clades is estimated in a short time frame during the Middle Miocene and might have been associated with a period of global climate cooling with the establishment of a marked climatic zonation in Europe. Cladogenetic events within the main clades are scattered throughout the time tree, from the Late Miocene to the Early Pleistocene, suggesting that speciation events in wall lizards reflect a complex interplay between regional topography, climate and geological history rather than a shared major climatic or paleogeographic event. Our absolute time estimates, as well as a relative dating approach, demonstrate that the assumption of a causal link between sea-level rise at the end of the MSC and the diversification of many island endemics is not justified. This study reinforces the notion that multiple dispersal and vicariant events, at different time frames, are required to explain current allopatric distributions and to account for the historical assembly of Mediterranean biota, and cautions against the use of biogeographic calibrations based on the assumption of vicariance.

The Galápagos Islands: biogeographic patterns and geology

In the traditional biogeographic model, the Galápagos Islands appeared a few million years ago in a sea where no other islands existed and were colonized from areas outside the region. However, recent work has shown that the Galápagos hotspot is 139 million years old (Early Cretaceous), and so groups are likely to have survived at the hotspot by dispersal of populations onto new islands from older ones. This process of metapopulation dynamics means that species can persist indefinitely in an oceanic region, as long as new islands are being produced. Metapopulations can also undergo vicariance into two metapopulations, for example at active island arcs that are rifted by transform faults. We reviewed the geographic relationships of Galápagos groups and found 10 biogeographic patterns that are shared by at least two groups. Each of the patterns coincides spatially with a major tectonic structure; these structures include: the East Pacific Rise; west Pacific and American subduction zones; large igneous plateaus in the Pacific; Alisitos terrane (Baja California), Guerrero terrane (western Mexico); rifting of North and South America; formation of the Caribbean Plateau by the Galápagos hotspot, and its eastward movement; accretion of Galápagos hotspot tracks; Andean uplift; and displacement on the Romeral fault system. All these geological features were active in the Cretaceous, suggesting that geological change at that time caused vicariance in widespread ancestors. The present distributions are explicable if ancestors survived as metapopulations occupying both the Galápagos hotspot and other regions before differentiating, more or less in situ.

Settling taxonomic and nomenclatural problems in brine shrimps, Artemia (Crustacea: Branchiopoda: Anostraca), by integrating mitogenomics, marker discordances and nomenclature rules

High morphological plasticity in populations of brine shrimp subjected to different environmental conditions, mainly salinity, hindered for centuries the identification of the taxonomic entities encompassed within Artemia. In addition, the mismatch between molecular and morphological evolution rates complicates the characterization of evolutionary lineages, generating taxonomic problems. Here, we propose a phylogenetic hypothesis for Artemia based on two new complete mitogenomes, and determine levels of congruence in the definition of evolutionary units using nuclear and mtDNA data. We used a fossil of Artemia to calibrate the molecular clock and discuss divergence times within the genus. The hypothesis proposed herein suggests a more recent time frame for lineage splitting than previously considered. Phylogeographic analyses were performed using GenBank available mitochondrial and nuclear markers. Evidence of gen e flow, identified through discordances between nuclear and mtDNA markers, was used to reconsider the specific status of some taxa. As a result, we consider Artemia to be represented by five evolutionary units: Southern Cone, Mediterranean-South African, New World, Western Asian, and Eastern Asian Lineages. After an exhaustive bibliographical revision, unavailable names for nomenclatural purposes were discarded. The remaining available names have been assigned to their respective evolutionary lineage. The proper names for the evolutionary units in which brine shrimps are structured remain as follows: Artemia persimilis Piccinelli & Prosdocimi, 1968 for the Southern Cone Lineage, Artemia salina (Linnaeus, 1758) for the Mediterranean-SouthAfrican Lineage, Artemia urmiana Günther, 1899 for the Western Asian Lineage, and Artemia sinica Cai, 1989 for the Eastern Asian Lineage. The name Artemia monica Verrill, 1869 has nomenclatural priority over A. franciscana Kellogg, 1906 for naming the New World Lineage. New synonymies are proposed for A. salina (= C. dybowskii Grochowski, 1896 n. syn., and A. tunisiana Bowen & Sterling, 1978 n. syn.), A. monica (= A. franciscana Kellogg, 1906 n. syn., and A. salina var. pacifica Sars, 1904 n. syn.); A. urmiana (= B. milhausenii Fischer de Waldheim, 1834 n. syn., A. koeppeniana Fischer, 1851 n. syn., A. proxima King, 1855 n. syn., A. s. var. biloba Entz, 1886 n. syn., A. s. var. furcata Entz, 1886 n. syn., A. asiatica Walter, 1887 n. syn., A. parthenogenetica Bowen & Sterling, 1978 n. syn., A. ebinurica Qian & Wang, 1992 n. syn., A. murae Naganawa, 2017 n. syn., and A. frameshifta Naganawa & Mura, 2017 n. syn.). Internal deep nuclear structuring within the A. monica and A. salina clades, might suggest the existence of additional evolutionary units within these taxa.

Genetic Diversity of Ligidium Isopods in Hokkaido and Niigata, Northern Japan, Based on Mitochondrial DNA Analysis

The genetic diversity of the genus Ligidium in Hokkaido and Niigata, northern Japan, was investigated by analyzing the cytochrome c oxidase subunit 1 (CO1) region in the mitochondrial DNA (mtDNA). The genetic diversity in Hokkaido was much lower than that in Niigata. Nine different operational taxonomic units (OTUs) were identified. Only a single OTU, most likely Ligidium japonicum, was found in Hokkaido, whereas all nine OTUs were found in Niigata. Using the mtDNA evolutionary rate determined for the marine invertebrate Haptosquilla pulchella (Miers, 1880), population expansion for OTU1 in Hokkaido was estimated to have occurred at 12,600 years BP, suggesting that Ligidium underwent a bottleneck due to glacial cooling, and the population then expanded after postglacial warming. Assuming that the expansion of the OTU1 population occurred at 9600 years BP, when the sea surface temperature rose offshore of Tokachi in the Northwestern Pacific, the evolutionary rate (µ) of the mtDNA CO1 region in Ligidium is calculated as: 0.087 (95% confidence intervals: min: 0.042-max: 0.12) (substitutions/site/million years). The presence of a haplotype common to Hokkaido and Niigata implies that the haplotype migrated across the Tsugaru Strait. Considering that geological evidence indicates that the Tsugaru Strait was continuously present even during the last glacial maximum when the sea level was at its lowest, accidental transport by human beings or animals might have been critical to the migration of Ligidium.

Geographic isolation and human-assisted dispersal in land snails: a Mediterranean story of Helix borealis and its relatives (Gastropoda: Stylommatophora: Helicidae)

The Mediterranean basin is a major centre for land-snail diversity, with many localized endemics, but there are also species widely spread by humans. Both endemics and introductions can be found in the snail genus Helix, which comprises many large-bodied species used for human consumption in the past and present. The Mediterranean clade of Helix is currently distributed throughout this region, but the phylogenetic and biogeographic relationships among its forms from different parts of the basin remain enigmatic. The reasons include insufficient sampling, taxa with unclear taxonomy and a significant impact of human-assisted transport obscuring the natural distribution of phylogenetic lineages. We provide evidence that European and Anatolian populations of H. cincta and its relatives are not native to those regions, but originate from the northern Levant. These results have implications for taxonomy of the genus, but also for the understanding of its evolutionary history. We posit that the Mediterranean clade consists of four geographically separated groups, which diversified in Northern Africa, the Apennine Peninsula and Corsica, the Aegean and Greece, and the northern Levant. This geographic pattern has been subsequently blurred by multiple instances of human-assisted dispersal. However, revealing the founding populations with certainty requires thorough sampling in currently inaccessible countries.

Global Diversification Dynamics Since the Jurassic: Low Dispersal and Habitat-Dependent Evolution Explain Hotspots of Diversity and Shell Disparity in River Snails (Viviparidae)

The Viviparidae, commonly known as River Snails, is a dominant group of freshwater snails with a nearly worldwide distribution that reaches its highest taxonomic and morphological diversity in Southeast Asia. The rich fossil record is indicative of a probable Middle Jurassic origin on the Laurasian supercontinent where the group started to diversify during the Cretaceous. However, it remains uncertain when and how the biodiversity hotspot in Southeast Asia was formed. Here, we used a comprehensive genetic data set containing both mitochondrial and nuclear markers and comprising species representing 24 out of 28 genera from throughout the range of the family. To reconstruct the spatiotemporal evolution of viviparids on a global scale, we reconstructed a fossil-calibrated phylogeny. We further assessed the roles of cladogenetic and anagenetic events in range evolution. Finally, we reconstructed the evolution of shell features by estimating ancestral character states to assess whether the appearance of sculptured shell morphologies was driven by major habitat shifts. The molecular phylogeny supports the monophyly of the three subfamilies, the Bellamyinae, Lioplacinae, and Viviparinae, but challenges the currently accepted genus-level classification in several cases. The almost global distribution of River Snails has been influenced both by comparatively ancient vicariance and more recent founder events. In Southeast Asia, Miocene dispersal was a main factor in shaping the modern species distributions. A recurrent theme across different viviparid taxa is that many species living in lentic waters exhibit sculptured shells, whereas only one strongly sculptured species is known from lotic environments. We show that such shell sculpture is habitat-dependent and indeed evolved several times independently in lentic River Snails. Considerably high transition rates between shell types in lentic habitats probably caused the co-occurrence of morphologically distinct shell types in several lakes. In contrast, directional evolution toward smooth shells in lotic habitats, as identified in the present analyses, explains why sculptured shells are rarely found in these habitats. However, the specific factors that promoted changes in shell morphology require further work. [biogeographical analyses; fossil-calibrated phylogeny; fossil-constrained analyses; Southeast Asia; stochastic character mapping.].

Influence of niche breadth and position on the historical biogeography of seafaring scincid lizards

Niche breadth and position can influence diversification among closely related species or populations, yet limited empirical data exist concerning the predictability of the outcomes. We explored the effects of these factors on the evolution of the Emoia atrocostata species group, an insular radiation of lizards in the western Pacific Ocean and Indo-Australasia composed of both endemic and widespread species that differ in niche occupancy. We used molecular data and phylogeographical diffusion models to estimate the timing and patterns of range expansion, and ancestral reconstruction methods to infer shifts in ecology. We show evidence of multidirectional spread from a centre of origin in western Micronesia, and that the phyletic diversity of the group is derived from a putative habitat specialist that survives in the littoral zone. This species is composed of paraphyletic lineages that represent stages or possible endpoints in the continuum toward speciation. Several descendant species have transitioned to either strand or interior forest habitat, but only on remote islands with depauperate terrestrial faunas. Our results suggest that the atrocostata group might be in the early phases of a Wilsonian taxon cycle and that the capacity to tolerate salt stress has promoted dispersal and colonization of remote oceanic islands. Divergence itself, however, is largely driven by geographical isolation rather than shifts in ecology.

Inferring Evolutionary Timescales without Independent Timing Information: An Assessment of "Universal" Insect Rates to Calibrate a Collembola (Hexapoda) Molecular Clock

Previous estimates of nucleotide substitution rates are routinely applied as secondary or “universal” molecular clock calibrations for estimating evolutionary timescales in groups that lack independent timing information. A major limitation of this approach is that rates can vary considerably among taxonomic groups, but the assumption of rate constancy is rarely evaluated prior to using secondary rate calibrations. Here I evaluate whether an insect mitochondrial DNA clock is appropriate for estimating timescales in Collembola—a group of insect-like arthropods characterized by high levels of cryptic diversity. Relative rates of substitution in cytochrome oxidase subunit 1 (COI) were inferred via Bayesian analysis across a topologically constrained Hexapod phylogeny using a relaxed molecular clock model. Rates for Collembola did not differ significantly from the average rate or from the rates estimated for most other groups (25 of 30), suggesting that (1) their apparent cryptic diversity cannot be explained by accelerated rates of molecular evolution and (2) clocks calibrated using “universal” insect rates may be appropriate for estimating evolutionary timescales in this group. However, of the 31 groups investigated, 10 had rates that deviated significantly from the average (6 higher, 4 lower), underscoring the need for caution and careful consideration when applying secondary insect rate calibrations. Lastly, this study exemplifies a relatively simple approach for evaluating rate constancy within a taxonomic group to determine whether the use of secondary rates are appropriate for molecular clock calibrations.

Rates and Rocks: Strengths and Weaknesses of Molecular Dating Methods

I present here an in-depth, although non-exhaustive, review of two topics in molecular dating. Clock models, which describe the evolution of the rate of evolution, are considered first. Some of the shortcomings of popular approaches-uncorrelated clock models in particular-are presented and discussed. Autocorrelated models are shown to be more reasonable from a biological perspective. Some of the most recent autocorrelated models also rely on a coherent treatment of instantaneous and average substitution rates while previous models are based on implicit approximations. Second, I provide a brief overview of the processes involved in collecting and preparing fossil data. I then review the main techniques that use this data for calibrating the molecular clock. I argue that, in its current form, the fossilized birth-death process relies on assumptions about the mechanisms underlying fossilization and the data collection process that may negatively impact the date estimates. Node-dating approaches make better use of the data available, even though they rest on paleontologists' intervention to prepare raw fossil data. Altogether, this study provides indications that may help practitioners in selecting appropriate methods for molecular dating. It will also hopefully participate in defining the contour of future methodological developments in the field.

Phylogeny of the Subtribe Monoctonina (Hymenoptera, Braconidae, Aphidiinae)

Members of the Monoctonina subtribe have long been neglected in applied studies of the subfamily Aphidiinae, due to their low economic importance, as they do not parasitize pests of cultivated plants. Consequently, data about this group are scarce, including its taxonomy and phylogeny. In the present study, we explore inter- and intraspecific genetic variation of Monoctonina species, including genera Monoctonus Haliday 1833, Monoctonia Starý 1962, Falciconus Mackauer 1959 and Harkeria Cameron 1900. We employ two molecular markers, the barcode region of the mitochondrial cytochrome c oxidase subunit I (COI) and the D2 region of the 28S nuclear gene (28S rDNA), to analyze genetic structuring and phylogeny of all available Monoctonina species, and combine them with morphological data for an integrative approach. We report one new species, and three potentially new species which can be formally described when further specimens are available. Analysis of phylogenetic relationships within the subtribe shows a basal position for the genera Falciconus and Monoctonia, and the close relatedness of Harkeria and Monoctonus.

A Simulation-Based Evaluation of Tip-Dating Under the Fossilized Birth-Death Process

Bayesian molecular dating is widely used to study evolutionary timescales. This procedure usually involves phylogenetic analysis of nucleotide sequence data, with fossil-based calibrations applied as age constraints on internal nodes of the tree. An alternative approach is tip-dating, which explicitly includes fossil data in the analysis. This can be done, for example, through the joint analysis of molecular data from present-day taxa and morphological data from both extant and fossil taxa. In the context of tip-dating, an important development has been the fossilized birth-death process, which allows non-contemporaneous tips and sampled ancestors while providing a model of lineage diversification for the prior on the tree topology and internal node times. However, tip-dating with fossils faces a number of considerable challenges, especially, those associated with fossil sampling and evolutionary models for morphological characters. We conducted a simulation study to evaluate the performance of tip-dating using the fossilized birth-death model. We simulated fossil occurrences and the evolution of nucleotide sequences and morphological characters under a wide range of conditions. Our analyses of these data show that the number and the maximum age of fossil occurrences have a greater influence than the degree of among-lineage rate variation or the number of morphological characters on estimates of node times and the tree topology. Tip-dating with the fossilized birth-death model generally performs well in recovering the relationships among extant taxa but has difficulties in correctly placing fossil taxa in the tree and identifying the number of sampled ancestors. The method yields accurate estimates of the ages of the root and crown group, although the precision of these estimates varies with the probability of fossil occurrence. The exclusion of morphological characters results in a slight overestimation of node times, whereas the exclusion of nucleotide sequences has a negative impact on inference of the tree topology. Our results provide an overview of the performance of tip-dating using the fossilized birth-death model, which will inform further development of the method and its application to key questions in evolutionary biology.

Using the Fossil Record to Evaluate Timetree Timescales

The fossil and geologic records provide the primary data used to established absolute timescales for timetrees. For the paleontological evaluation of proposed timetree timescales, and for node-based methods for constructing timetrees, the fossil record is used to bracket divergence times. Minimum brackets (minimum ages) can be established robustly using well-dated fossils that can be reliably assigned to lineages based on positive morphological evidence. Maximum brackets are much harder to establish, largely because it is difficult to establish definitive evidence that the absence of a taxon in the fossil record is real and not just due to the incompleteness of the fossil and rock records. Five primary methods have been developed to estimate maximum age brackets, each of which is discussed. The fact that the fossilization potential of a group typically decreases the closer one approaches its time of origin increases the challenge of estimating maximum age brackets. Additional complications arise: 1) because fossil data actually bracket the time of origin of the first relevant fossilizable morphology (apomorphy), not the divergence time itself; 2) due to the phylogenetic uncertainty in the placement of fossils; 3) because of idiosyncratic temporal and geographic gaps in the rock and fossil records; and 4) if the preservation potential of a group changed significantly during its history. In contrast, uncertainties in the absolute ages of fossils are typically relatively unimportant, even though the vast majority of fossil cannot be dated directly. These issues and relevant quantitative methods are reviewed, and their relative magnitudes assessed, which typically correlate with the age of the group, its geographic range, and species richness.

Metapopulation Vicariance, Age of Island Taxa and Dispersal: A Case Study Using the Pacific Plant Genus Planchonella (Sapotaceae)

Oceanic islands originate from volcanism or tectonic activity without connections to continental landmasses, are colonized by organisms, and eventually vanish due to erosion and subsidence. Colonization of oceanic islands occurs through long-distance dispersals (LDDs) or metapopulation vicariance, the latter resulting in lineages being older than the islands they inhabit. If metapopulation vicariance is valid, island ages cannot be reliably used to provide maximum age constraints for molecular dating. We explore the relationships between the ages of members of a widespread plant genus (Planchonella, Sapotaceae) and their host islands across the Pacific to test various assumptions of dispersal and metapopulation vicariance. We sampled three nuclear DNA markers from 156 accessions representing some 100 Sapotaceae taxa, and analyzed these in BEAST with a relaxed clock to estimate divergence times and with a phylogeographic diffusion model to estimate range expansions over time. The phylogeny was calibrated with a secondary point (the root) and fossils from New Zealand. The dated phylogeny reveals that the ages of Planchonella species are, in most cases, consistent with the ages of the islands they inhabit. Planchonella is inferred to have originated in the Sahul Shelf region, to which it back-dispersed multiple times. Fiji has been an important source for range expansion in the Pacific for the past 23 myr. Our analyses reject metapopulation vicariance in all cases tested, including between oceanic islands, evolution of an endemic Fiji-Vanuatu flora, and westward rollback vicariance between Vanuatu and the Loyalty Islands. Repeated dispersal is the only mechanism able to explain the empirical data. The longest (8900 km) identified dispersal is between Palau in the Pacific and the Seychelles in the Indian Ocean, estimated at 2.2 Ma (0.4-4.8 Ma). The first split in a Hawaiian lineage (P. sandwicensis) matches the age of Necker Island (11.0 Ma), when its ancestor diverged into two species that are distinguished by purple and yellow fruits. Subsequent establishment across the Hawaiian archipelago supports, in part, progression rule colonization. In summary, we found no explanatory power in metapopulation vicariance and conclude that Planchonella has expanded its range across the Pacific by LDD. We contend that this will be seen in many other groups when analyzed in detail.

Diet-driven ecological radiation and allopatric speciation result in high species diversity in a temperate-cold water marine genus Dendronotus (Gastropoda: Nudibranchia)

While the majority nudibranch clades are more species rich in the tropics, the genus Dendronotus is mainly represented in Arctic and boreal regions. This distribution pattern remains poorly understood. An integrative approach and novel data provided valuable insights into processes driving Dendronotus radiation and speciation. We propose an evolutionary scenario based on molecular phylogenetics and morphological, ecological, ontogenetic data, combined with data on complex geology and paleoclimatology of this region. Estimated phylogenetic relationships based on four molecular markers (COI, 16S, H3 and 28S) shows strong correlation with radular morphology, diet and biogeographical pattern. Ancestral area reconstruction (AAR) provides evidence for a tropical Pacific origin of the genus. Based on AAR and divergence time estimates we conclude that the evolution of Dendronotus has been shaped by different processes: initial migration out of the tropics, diet-driven adaptive radiation in the North Pacific influenced by Miocene climate change, and subsequent allopatric speciation resulting from successive closings of the Bering strait and cooling of the Arctic Ocean during the Pliocene-Pleistocene. At the same time, contemporary amphiboreal species appear to have dispersed into the Atlantic fairly recently.

Insights into the molecular phylogeny of Rhaphidophoridae, an ancient, worldwide lineage of Orthoptera

We investigated the molecular phylogenetic divergence and historical biogeography of cave crickets belonging to the family Rhaphidophoridae (Orthoptera, Ensifera). We used taxa representative of most of the regions embraced by the family, considering samples of Macropathinae from Gondwana land (i.e., Tasmania, Australia, New Zealand, South Africa, and South America); Aemodogryllinae and Rhaphidophorinae from Southern-eastern Asia (i.e., India, Bhutan, China, Philippines and the Sulawesi islands); Dolichopodainae and Troglophilinae from the Mediterranean region and Ceuthophilinae from North America. Based on previous papers, we carried out an analysis of both mitochondrial and nuclear DNA sequences considering the ribosomal RNA units 12S, 16S, 18S, and 28S. To reconstruct phylogeny, we use cladistics, Maximum Likelihood (ML), and Bayesian analyses. All phylogenetic analyses showed the same highly supported topology generally congruent with the classical systematic arrangement at the level of each sub-family but strongly disagree with previous affinity hypotheses between sub-families based on morphological characters. Our results reveal a close affinity between Asiatic and Gondwanian taxa from one hand and between North American and Mediterranean ones from the other hand. Dating estimates indicated that Rhaphidophoridae originated in the Cretaceous period during the Mesozoic era with the ancestral area located both in the northern and southern hemisphere. A possible biogeographic scenario, reconstructed using S-DEC with RASP software, suggested that the current distribution of Rhaphidophoridae might be explained by a combination of both dispersal and vicariance events occurred especially in the ancestral populations. The radiation of Rhaphidophoridae started within the Pangaea, where the ancestor of Rhaphidophoridae occurred throughout an ancestral area including Australia, North America, and the Mediterranean region. The opening of the Atlantic Ocean promoted the divergence of North American and Mediterranean lineages while the differentiation of the southern lineages, spread from Australia, appears to be related to the fragmentation of Gondwana land.

The evolutionary history of the Cape hare (Lepus capensis sensu lato): insights for systematics and biogeography

Inferring the phylogeography of species with large distributions helps deciphering major diversification patterns that may occur in parallel across taxa. Here, we infer the evolutionary history of the Cape hare, Lepus capensis sensu lato, a species distributed from southern Africa to Asia, by analyzing variation at 18 microsatellites and 9 DNA (1 mitochondrial and 8 nuclear) sequenced loci, from field and museum-collected samples. Using a combination of assignment and coalescent-based methods, we show that the Cape hare is composed of five evolutionary lineages, distributed in distinct biogeographic regions-north-western Africa, eastern Africa, southern Africa, the Near East and the Arabian Peninsula. A deep phylogenetic break possibly dating to the Early Pleistocene was inferred between the African and Asian L. capensis groups, and the latter appear more closely related to other Eurasian hare species than to African Cape hares. The inferred phylogeographic structure is shared by numerous taxa distributed across the studied range, suggesting that environmental changes, such as the progressive aridification of the Saharo-Arabian desert and the fluctuations of savannah habitats in Sub-Saharan Africa, had comparable impacts across species. Fine-scale analyses of the western Sahara-Sahel populations showed rich fragmentation patterns for mitochondrial DNA but not for microsatellites, compatible with the environmental heterogeneity of the region and female philopatry. The complex evolutionary history of L. capensis sensu lato, which possibly includes interspecific gene flow, is not reflected by taxonomy. Integrating evolutionary inference contributes to an improved characterization of biodiversity, which is fundamental to foster the conservation of relevant evolutionary units.

Estimation of Evolutionary Rates of Mitochondrial DNA in Two Japanese Wood Mouse Species Based on Calibrations with Quaternary Environmental Changes

Reliable estimates of evolutionary rates of mitochondrial DNA might allow us to build realistic evolutionary scenarios covering broad time scales based on phylogenetic inferences. In the present study, we sought to obtain estimates of evolutionary rates in murine rodents using calibrations against historical biogeographic events. We first assumed that land-bridge-like structures that appeared intermittently at glacial maxima with 100,000-year intervals shaped the divergence patterns of cytochrome b (Cytb) sequences (1140 bp) of the larger Japanese wood mouse Apodemus speciosus. The comparison of sequences from peripheral remote islands that are separated from one another by deep straits allowed us to estimate mitochondrial DNA evolutionary rates (substitutions/site/million years) to be 0.027 to 0.036, with presumed calibrations from 140,000, 250,000, 350,000, and 440,000 years ago. Second, we addressed rapid expansion events inferred from analyses of the Cytb sequences of the lesser Japanese wood mouse A. argenteus. We detected five expansion signals in the dataset and established three categories based on the expansion parameter tau values: 3.9, 5.6-5.7, and 7.8-8.1. Considering that the climate became warmer 15,000, 53,000, and 115,000 years ago after preceding periods of rapid cooling, we calculated evolutionary rates to be 0.114, 0.047, and 0.031, respectively. This preliminary concept of the evolutionary rates on a time scale from 15,000 to 440,000 years ago for the wood mouse should be refined and tested in other species of murine rodents, including mice and rats.

Ancestral range reconstruction of remote oceanic island species of Plantago (Plantaginaceae) reveals differing scales and modes of dispersal

AIM

The aim of this study was to resolve the phylogenetic placement of island taxa, reconstruct ancestral origins and resolve competing hypotheses of dispersal patterns and biogeographical histories for oceanic island endemic taxa within subgenus Plantago (Plantaginaceae).

LOCATION

Juan Fernández Islands, the Auckland Islands, Lord Howe Island, New Amsterdam Island, New Zealand, Tasmania, Falkland Islands, Rapa Iti and the Hawaiian Islands.

TAXON

Island endemics within Plantago (Plantaginaceae), a globally distributed taxonomic group comprising approximately 250 species.

METHODS

We use Bayesian phylogenetic and divergence time analyses and historical biogeographical analysis of molecular sequence data to infer the ancestral origins of the oceanic island species in Plantago.

RESULTS

Taxa within subgenus Plantago form clades based on geographic proximities and challenge previous phylogenetic relationships and classification based on morphology. We infer that biogeographic histories of oceanic island taxa from multiple islands were shaped by dispersal at different scales and possibly by different types of birds. The highly remote Hawaiian Islands and Rapa Iti were colonized from North American taxa in a pattern corresponding to known migration routes of large marine birds, rather than from New Zealand as previously hypothesized. The island endemics of Juan Fernández, the Falkland Islands, Lord Howe, Auckland Islands and New Zealand are found to have sources in the nearest continental areas. The analyses confirm recent speciation within subgenus Plantago - which is particularly heightened in island lineages in Hawaii and Rapa Iti - but show slightly older divergence times than previous molecular dating studies.

MAIN CONCLUSIONS

Using molecular data to infer ancestral ranges for plants with uncertain taxonomic relationships can greatly improve our understanding of biogeographical histories and help elucidate origins, dispersal modes and routes in widespread lineages with complex distribution patterns such as Plantago. We improve understanding of important floristic exchange areas between continents and islands as a result of long-distance dispersal. We infer that a combination of both stepping stone dispersal and extreme long-distance dispersal can shape insular floras, and that multiple floristic areas can be the sources of closely related island taxa. However, despite the successful dispersal of Plantago, radiation in island archipelagos is generally limited suggesting specific traits may limit diversification.

A Not-So-Long Introduction to Computational Molecular Evolution

In this chapter, we give a not-so-long and self-contained introduction to computational molecular evolution. In particular, we present the emergence of the use of likelihood-based methods, review the standard DNA substitution models, and introduce how model choice operates. We also present recent developments in inferring absolute divergence times and rates on a phylogeny, before showing how state-of-the-art models take inspiration from diffusion theory to link population genetics, which traditionally focuses at a taxonomic level below that of the species, and molecular evolution. Although this is not a cookbook chapter, we try and point to popular programs and implementations along the way.

Recent advances in New Caledonian biogeography

The biota of New Caledonia is one of the most unusual in the world. It displays high diversity and endemism, many peculiar absences, and far-flung biogeographic affinities. For example, New Caledonia is the only place on Earth with both main clades of flowering plants - the endemic Amborella and 'all the rest', and it also has the highest concentration of diversity in conifers. The discovery of Amborella's phylogenetic position led to a surge of interest in New Caledonian biogeography, and new studies are appearing at a rapid rate. This paper reviews work on the topic (mainly molecular studies) published since 2013. One current debate is focused on whether any biota survived the marine transgressions of the Paleocene and Eocene. Total submersion would imply that the entire fauna was derived by long-distance dispersal from continental areas since the Eocene, but only if no other islands (now submerged) were emergent. A review of the literature suggests there is little actual evidence in geology for complete submersion. An alternative explanation for New Caledonia's diversity is that the archipelago acted as a refugium, and that the biota avoided the extinctions that occurred in Australia. However, this is contradicted by the many groups that are anomalously absent or depauperate in New Caledonia, although represented there by a sister group. The anomalous absences, together with the unusual levels of endemism, can both be explained by vicariance at breaks in and around New Caledonia. New Caledonia has always been situated at or near a plate boundary, and its complex geological history includes the addition of new terranes (by accretion), orogeny, and rifting. New Caledonia comprises 'basement' terranes that were part of Gondwana, as well as island arc and forearc terranes that accreted to the basement after it separated from Gondwana. The regional tectonic history helps explain the regional biogeography, as well as distribution patterns within New Caledonia. These include endemics on the basement terranes (for example, the basal angiosperm, Amborella), disjunctions at the West Caledonian fault zone, and great biotic differences between Grande Terre and the Loyalty Islands.

Origin and macroevolution of micro-moths on sunken Hawaiian Islands

The origins and evolution of Hawaiian biodiversity are a matter of controversy, and the mechanisms of lineage diversification for many organisms on this remote archipelago remain unclear. Here we focus on the poorly known endemic leaf-mining moth genus Philodoria (Lepidoptera, Gracillariidae), whose species feed on a diversity of Hawaiian plant lineages, many of which are critically endangered. We use anchored hybrid enrichment to assemble the first phylogenomic dataset (507 loci) for any Hawaiian animal taxon. To uncover the timing and pattern of diversification of these moths, we apply two frequently used dating calibration strategies, biogeographic calibrations and secondary calibrations. Island calibrations on their own resulted in much younger and unrealistic dates compared to strategies that relied on secondary calibrations. Philodoria probably originated on the now partially sunken islands of Laysan or Lisianski, approximately 21 Ma, and were associated with host plants in the families Ebenaceae, Malvaceae or Primulaceae. Major feeding groups associated with specific host-plant families originated soon after the plants colonized the islands. Allopatric isolation and host shifts, in concert and independently, probably play major roles in the diversification of Philodoria Our dating results indicate that Philodoria is among the oldest known Hawaiian arthropod lineages, and that island calibrations alone can lead to unrealistically young dates.

Temperature gradient affects differentiation of gene expression and SNP allele frequencies in the dominant Lake Baikal zooplankton species

Local adaptation and phenotypic plasticity are main mechanisms of organisms' resilience in changing environments. Both are affected by gene flow and are expected to be weak in zooplankton populations inhabiting large continuous water bodies and strongly affected by currents. Lake Baikal, the deepest and one of the coldest lakes on Earth, experienced epilimnion temperature increase during the last 100 years, exposing Baikal's zooplankton to novel selective pressures. We obtained a partial transcriptome of Epischura baikalensis (Copepoda: Calanoida), the dominant component of Baikal's zooplankton, and estimated SNP allele frequencies and transcript abundances in samples from regions of Baikal that differ in multiyear average surface temperatures. The strongest signal in both SNP and transcript abundance differentiation is the SW-NE gradient along the 600+ km long axis of the lake, suggesting isolation by distance. SNP differentiation is stronger for nonsynonymous than synonymous SNPs and is paralleled by differential survival during a laboratory exposure to increased temperature, indicating directional selection operating on the temperature gradient. Transcript abundance, generally collinear with the SNP differentiation, shows samples from the warmest, less deep location clustering together with the southernmost samples. Differential expression is more frequent among transcripts orthologous to candidate thermal response genes previously identified in model arthropods, including genes encoding cytoskeleton proteins, heat-shock proteins, proteases, enzymes of central energy metabolism, lipid and antioxidant pathways. We conclude that the pivotal endemic zooplankton species in Lake Baikal exists under temperature-mediated selection and possesses both genetic variation and plasticity to respond to novel temperature-related environmental pressures.

A revised phylogeny of macropathine cave crickets (Orthoptera: Rhaphidophoridae) uncovers a paraphyletic Australian fauna

Australian cave crickets are members of the subfamily Macropathinae (Orthoptera: Rhaphidophoridae). The subfamily is thought to have originated prior to the tectonic separation of the supercontinent Gondwana based on distributions of extant lineages and molecular phylogenetic evidence, although the Australian fauna have been underrepresented in previous studies. The current study augments existing multigene data (using 12S, 16S, and 28S rRNA genes) to investigate the placement of the Australian representatives within the Macropathinae and to assess divergence dates of select clades. Results suggest that the endemic Tasmanian genus Parvotettix is the sister lineage to the remaining members of the subfamily, an outcome that presents a paraphyletic Australian fauna in contrast to previous studies. All other Australian taxa represented in this study (Micropathus and Novotettix) emerged as a sister group to the New Zealand and South American macropathine lineages. Estimation of phylogenetic divergence ages among the aforementioned clades were calibrated using two methods, in absence of suitable fossil records: (i) tectonic events depicting the fragmentation of Gondwanan landmasses that invoke vicariant scenarios of present day geographic distributions; and (ii) molecular evolutionary rates. Geological calibrations place the median age of the most recent common ancestor of extant macropathines at ∼125 to ∼165 Ma, whereas analyses derived from molecular substitution rates suggest a considerably younger origin of ∼32 Ma. This phylogenetic study represents the most rigorous taxonomic sampling of the Australian cave cricket fauna to date and stresses the influence of lineage representation on biogeographic inference.

Phylogeography of Daphnia magna Straus (Crustacea: Cladocera) in Northern Eurasia: Evidence for a deep longitudinal split between mitochondrial lineages

Species with a large geographic distributions present a challenge for phylogeographic studies due to logistic difficulties of obtaining adequate sampling. For instance, in most species with a Holarctic distribution, the majority of studies has concentrated on the European or North American part of the distribution, with the Eastern Palearctic region being notably understudied. Here, we study the phylogeography of the freshwater cladoceran Daphnia magna Straus, 1820 (Crustacea: Cladocera), based on partial mitochondrial COI sequences and using specimens from populations spread longitudinally from westernmost Europe to easternmost Asia, with many samples from previously strongly understudied regions in Siberia and Eastern Asia. The results confirm the previously suspected deep split between Eastern and Western mitochondrial haplotype super-clades. We find a narrow contact zone between these two super-clades in the eastern part of Western Siberia, with proven co-occurrence in a single lake in the Novosibirsk region. However, at present there is no evidence suggesting that the two mitochondrial super-clades represent cryptic species. Rather, they may be explained by secondary contact after expansion from different refugia. Interestingly, Central Siberia has previously been found to be an important contact zone also in other cladoceran species, and may thus be a crucial area for understanding the Eurasian phylogeography of freshwater invertebrates. Together, our study provides an unprecedented complete, while still not global, picture of the phylogeography of this important model species.

Evolutionary Patterns and Processes: Lessons from Ancient DNA

Ever since its emergence in 1984, the field of ancient DNA has struggled to overcome the challenges related to the decay of DNA molecules in the fossil record. With the recent development of high-throughput DNA sequencing technologies and molecular techniques tailored to ultra-damaged templates, it has now come of age, merging together approaches in phylogenomics, population genomics, epigenomics, and metagenomics. Leveraging on complete temporal sample series, ancient DNA provides direct access to the most important dimension in evolution—time, allowing a wealth of fundamental evolutionary processes to be addressed at unprecedented resolution. This review taps into the most recent findings in ancient DNA research to present analyses of ancient genomic and metagenomic data.

Phylogeography of the sand dollar genus Encope: implications regarding the Central American Isthmus and rates of molecular evolution

Vicariant events have been widely used to calibrate rates of molecular evolution, the completion of the Central American Isthmus more extensively than any other. Recent studies have claimed that rather than the generally accepted date of ~3 million years ago (Ma), the Isthmus was effectively complete by the middle Miocene, 13 Ma. We present a fossil calibrated phylogeny of the new world sand dollar genus Encope, based on one nuclear and four mitochondrial genes, calibrated with fossils at multiple nodes. Present day distributions of Encope are likely the result of multiple range contractions and extinction events. Most species are now endemic to a single region, but one widely distributed species in each ocean is composed of morphotypes previously described as separate species. The most recent separation between eastern Pacific and Caribbean extant clades occurred at 4.90 Ma, indicating that the Isthmus of Panama allowed genetic exchange until the Pliocene. The rate of evolution of mitochondrial genes in Encope has been ten times slower than in the closely related genera Mellita and Lanthonia. This large difference in rates suggests that splits between eastern Pacific and Caribbean biota, dated on the assumption of a "universal" mitochondrial DNA clock are not valid.

Transoceanic origin of microendemic and flightless New Caledonian weevils

The origin of the astonishing New Caledonian biota continues to fuel a heated debate among advocates of a Gondwanan relict scenario and defenders of late oceanic dispersal. Here, we study the origin of New Caledonian Trigonopterus flightless weevils using a multimarker molecular phylogeny. We infer two independent clades of species found in the archipelago. Our dating estimates suggest a Late Miocene origin of both clades long after the re-emergence of New Caledonia about 37 Ma. The estimation of ancestral ranges supports an ancestral origin of the genus in a combined region encompassing Australia and New Guinea with subsequent colonizations of New Caledonia out of New Guinea in the mid-Miocene. The two New Caledonian lineages have had very different evolutionary trajectories. Colonizers belonging to a clade of foliage dwellers greatly diversified, whereas species inhabiting leaf-litter have been less successful.

Patterns of divergence in fish species separated by the Isthmus of Panama

Background

The Pleistocene closure of Isthmus of Panama, separating the basins of the Eastern Pacific and the Caribbean Sea, created a unique natural experiment that reveals how marine faunas respond to environmental change. To explore how fishes have been affected by this tectonic event, I compare transisthmian patterns in phylogeny and morphology for geminate lineages in two families, Eleotridae (sleepers) and Apogonidae (cardinalfishes).

Results

Time-calibrated phylogenies for these families show different diversification patterns. In Eleotridae, several independent shallow instances of transisthmian divergences occur, with one or a few species on either side of the Isthmus. Among Apogonidae, a single clade of Eastern Pacific species is nested within a broad Caribbean radiation that also includes a species known from the Mediterranean. Divergence time estimates for taxa isolated by closure of the Isthmus are broadly congruent. Hypotheses dated with deeper, fossil-based legacy calibrations put the divergences in the Miocene at 7.4–15.1 Ma, while those estimated with a shallow biogeographic calibration of final Isthmus closure range from 5.1 to 9.9 Ma, in the late Miocene/early Pliocene. Eleotridae are more euryhaline than Apogonidae, but do not exhibit shallower transisthmian divergences. In both families, descendent lineages on either side of the Isthmus of Panama exhibit significant shape differences, although that distinction disappears for Apogonidae when I apply a correction for phylogenetic relationships. To evaluate the tempo and mode of continuous character evolution, I fit several single and multiple rate evolutionary models to morphometric data reconstructed on the Apogonidae phylogeny. I find that the most highly favored model, as estimated on both legacy and isthmus calibrated hypotheses, is a multiple rate Ornstein-Uhlbeck model, with a mosaic of rate shifts postulated for shape changes among fishes in the Caribbean and Eastern Pacific.

Conclusions

Although many transisthmian taxa have been compared and their phylogenies calibrated to estimate the dates associated with population sundering, few studies correlate these timing estimates with morphological change. I show that in transisthmian fish lineages, morphometric distinctions are detectable across the Isthmus, and that rates and patterns of shape change have also shifted, with variable manifestations across the body and between the Caribbean and Eastern Pacific.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0957-4) contains supplementary material, which is available to authorized users.