Quantifying the Error of Secondary vs. Distant Primary Calibrations in a Simulated Environment

Using calibrations to obtain absolute divergence times is standard practice in molecular clock studies. While the use of primary (e.g., fossil) calibrations is preferred, this approach can be limiting because of their rarity in fast-growing datasets. Thus, alternatives need to be explored, such as the use of secondary (molecularly-derived) calibrations that can anchor a timetree in a larger number of nodes. However, the use of secondary calibrations has been discouraged in the past because of concerns in the error rates of the node estimates they produce with an apparent high precision. Here, we quantify the amount of errors in estimates produced by the use of secondary calibrations relative to true times and primary calibrations placed on distant nodes. We find that, overall, the inaccuracies in estimates based on secondary calibrations are predictable and mirror errors associated with primary calibrations and their confidence intervals. Additionally, we find comparable error rates in estimated times from secondary calibrations and distant primary calibrations, although the precision of estimates derived from distant primary calibrations is roughly twice as good as that of estimates derived from secondary calibrations. This suggests that increasing dataset size to include primary calibrations may produce divergence times that are about as accurate as those from secondary calibrations, albeit with a higher precision. Overall, our results suggest that secondary calibrations may be useful to explore the parameter space of plausible evolutionary scenarios when compared to time estimates obtained with distant primary calibrations.

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.

Accounting for Uncertainty in the Evolutionary Timescale of Green Plants Through Clock-Partitioning and Fossil Calibration Strategies

Establishing an accurate evolutionary timescale for green plants (Viridiplantae) is essential to understanding their interaction and coevolution with the Earth's climate and the many organisms that rely on green plants. Despite being the focus of numerous studies, the timing of the origin of green plants and the divergence of major clades within this group remain highly controversial. Here, we infer the evolutionary timescale of green plants by analyzing 81 protein-coding genes from 99 chloroplast genomes, using a core set of 21 fossil calibrations. We test the sensitivity of our divergence-time estimates to various components of Bayesian molecular dating, including the tree topology, clock models, clock-partitioning schemes, rate priors, and fossil calibrations. We find that the choice of clock model affects date estimation and that the independent-rates model provides a better fit to the data than the autocorrelated-rates model. Varying the rate prior and tree topology had little impact on age estimates, with far greater differences observed among calibration choices and clock-partitioning schemes. Our analyses yield date estimates ranging from the Paleoproterozoic to Mesoproterozoic for crown-group green plants, and from the Ediacaran to Middle Ordovician for crown-group land plants. We present divergence-time estimates of the major groups of green plants that take into account various sources of uncertainty. Our proposed timeline lays the foundation for further investigations into how green plants shaped the global climate and ecosystems, and how embryophytes became dominant in terrestrial environments.

Priors and Posteriors in Bayesian Timing of Divergence Analyses: The Age of Butterflies Revisited

The need for robust estimates of times of divergence is essential for downstream analyses, yet assessing this robustness is still rare. We generated a time-calibrated genus-level phylogeny of butterflies (Papilionoidea), including 994 taxa, up to 10 gene fragments and an unprecedented set of 12 fossils and 10 host-plant node calibration points. We compared marginal priors and posterior distributions to assess the relative importance of the former on the latter. This approach revealed a strong influence of the set of priors on the root age but for most calibrated nodes posterior distributions shifted from the marginal prior, indicating significant information in the molecular data set. Using a very conservative approach we estimated an origin of butterflies at 107.6 Ma, approximately equivalent to the latest Early Cretaceous, with a credibility interval ranging from 89.5 Ma (mid Late Cretaceous) to 129.5 Ma (mid Early Cretaceous). In addition, we tested the effects of changing fossil calibration priors, tree prior, different sets of calibrations and different sampling fractions but our estimate remained robust to these alternative assumptions. With 994 genera, this tree provides a comprehensive source of secondary calibrations for studies on butterflies.

The Prevalence and Impact of Model Violations in Phylogenetic Analysis

In phylogenetic inference, we commonly use models of substitution which assume that sequence evolution is stationary, reversible, and homogeneous (SRH). Although the use of such models is often criticized, the extent of SRH violations and their effects on phylogenetic inference of tree topologies and edge lengths are not well understood. Here, we introduce and apply the maximal matched-pairs tests of homogeneity to assess the scale and impact of SRH model violations on 3,572 partitions from 35 published phylogenetic data sets. We show that roughly one-quarter of all the partitions we analyzed (23.5%) reject the SRH assumptions, and that for 25% of data sets, tree topologies inferred from all partitions differ significantly from topologies inferred using the subset of partitions that do not reject the SRH assumptions. This proportion increases when comparing trees inferred using the subset of partitions that rejects the SRH assumptions, to those inferred from partitions that do not reject the SRH assumptions. These results suggest that the extent and effects of model violation in phylogenetics may be substantial. They highlight the importance of testing for model violations and possibly excluding partitions that violate models prior to tree reconstruction. Our results also suggest that further effort in developing models that do not require SRH assumptions could lead to large improvements in the accuracy of phylogenomic inference. The scripts necessary to perform the analysis are available in https://github.com/roblanf/SRHtests, and the new tests we describe are available as a new option in IQ-TREE (http://www.iqtree.org).

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.

The effect of alignment uncertainty, substitution models and priors in building and dating the mammal tree of life

BACKGROUND

The flood of genomic data to help build and date the tree of life requires automation at several critical junctures, most importantly during sequence assembly and alignment. It is widely appreciated that automated alignment protocols can yield inaccuracies, but the relative impact of various sources error on phylogenomic analysis is not yet known. This study employs an updated mammal data set of 5162 coding loci sampled from 90 species to evaluate the effects of alignment uncertainty, substitution models, and fossil priors on gene tree, species tree, and divergence time estimation. Additionally, a novel coalescent likelihood ratio test is introduced for comparing competing species trees against a given set of gene trees.

RESULTS

The aligned DNA sequences of 5162 loci from 90 species were trimmed and filtered using trimAL and two filtering protocols. The final dataset contains 4 sets of alignments - before trimming, after trimming, filtered by a recently proposed pipeline, and further filtered by comparing ML gene trees for each locus with the concatenation tree. Our analyses suggest that the average discordance among the coalescent trees is significantly smaller than that among the concatenation trees estimated from the 4 sets of alignments or with different substitution models. There is no significant difference among the divergence times estimated with different substitution models. However, the divergence dates estimated from the alignments after trimming are more recent than those estimated from the alignments before trimming.

CONCLUSIONS

Our results highlight that alignment uncertainty of the updated mammal data set and the choice of substitution models have little impact on tree topologies yielded by coalescent methods for species tree estimation, whereas they are more influential on the trees made by concatenation. Given the choice of calibration scheme and clock models, divergence time estimates are robust to the choice of substitution models, but removing alignments deemed problematic by trimming algorithms can lead to more recent dates. Although the fossil prior is important in divergence time estimation, Bayesian estimates of divergence times in this data set are driven primarily by the sequence data.

Evolutionary Insight into the Trypanosomatidae Using Alignment-Free Phylogenomics of the Kinetoplast

Advancements in next-generation sequencing techniques have led to a substantial increase in the genomic information available for analyses in evolutionary biology. As such, this data requires the exponential growth in bioinformatic methods and expertise required to understand such vast quantities of genomic data. Alignment-free phylogenomics offer an alternative approach for large-scale analyses that may have the potential to address these challenges. The evolutionary relationships between various species within the trypanosomatid family, specifically members belonging to the genera Leishmania and Trypanosoma have been extensively studies over the last 30 years. However, there is a need for a more exhaustive analysis of the Trypanosomatidae, summarising the evolutionary patterns amongst the entire family of these important protists. The mitochondrial DNA of the trypanosomatids, better known as the kinetoplast, represents a valuable taxonomic marker given its unique presence across all kinetoplastid protozoans. The aim of this study was to validate the reliability and robustness of alignment-free approaches for phylogenomic analyses and its applicability to reconstruct the evolutionary relationships between the trypanosomatid family. In the present study, alignment-free analyses demonstrated the strength of these methods, particularly when dealing with large datasets compared to the traditional phylogenetic approaches. We present a maxicircle genome phylogeny of 46 species spanning the trypanosomatid family, demonstrating the superiority of the maxicircle for the analysis and taxonomic resolution of the Trypanosomatidae.

Molecular phylogeny and historical biogeography of Parnara butterflies (Lepidoptera: Hesperiidae)

The butterfly genus Parnara (Hesperiinae: Baorini), of which some are major pests of economic crops (e.g., rice, wild rice stems and sugarcane), currently consists of 10 species and several subspecies and has a highly disjunct distribution in Australia, Africa, and Asia. We determined the systematic relationships and biogeographical history of the genus by reconstructing the phylogeny based on eight genes and 101 specimens representing all 10 recognized species. Four species delimitation methods (ABGD, bPTP, GMYC and BPP) were also employed to assess the taxonomic status of each species. Based on these results and analyses, we recognize 11 extant species in the genus. The status of the taxon P. naso poutieri (Boisduval, 1833) from Madagascar is revised as a distinct species, Parnara poutieri (Boisduval, 1833) stat. rev. The subspecies P. guttata mangala (Moore, 1866) syn. nov. is synonymized with P. guttata guttata (Bremer & Grey, 1853), while P. bada (Moore, 1878) is provisionally treated as a complex of two species, namely P. bada and P. apostata (Snellen, 1886). The monophyly of Parnara is strongly supported, with the following relationships: P. amalia + ((P. monasi + (P. poutieri + P. naso)) + ((P. kawazoei + P. bada complex) + (P. ganga + (P. ogasawarensis + (P. guttata + P. batta))))). Divergence time and ancestral range estimates indicate that the common ancestor of Parnara originated in an implausible area of Australia, Africa, and Oriental region in the mid-Oligocene and then differentiated in the late Miocene-late Pliocene. Dispersal and range expansion have played an important role in diversification of the genus in Asia and Afica. Relatively stable geotectonic plates at the time when most extant lineages appeared during the late Miocene-early Pliocene might have been the factor responsible for the relatively constant low dynamic rate of diversification within the group.

Changing Ecological Opportunities Facilitated the Explosive Diversification of New Caledonian Oxera (Lamiaceae)

Phylogenies recurrently demonstrate that oceanic island systems have been home to rapid clade diversification and adaptive radiations. The existence of adaptive radiations posits a central role of natural selection causing ecological divergence and speciation, and some plant radiations have been highlighted as paradigmatic examples of such radiations. However, neutral processes may also drive speciation during clade radiations, with ecological divergence occurring following speciation. Here, we document an exceptionally rapid and unique radiation of Lamiaceae within the New Caledonian biodiversity hotspot. Specifically, we investigated various biological, ecological, and geographical drivers of species diversification within the genus Oxera. We found that Oxera underwent an initial process of rapid cladogenesis likely triggered by a dramatic period of aridity during the early Pliocene. This early diversification of Oxera was associated with an important phase of ecological diversification triggered by significant shifts of pollination syndromes, dispersal modes, and life forms. Finally, recent diversification of Oxera appears to have been further driven by the interplay of allopatry and habitat shifts likely related to climatic oscillations. This suggests that Oxera could be regarded as an adaptive radiation at an early evolutionary stage that has been obscured by more recent joint habitat diversification and neutral geographical processes. Diversification within Oxera has perhaps been triggered by varied ecological and biological drivers acting in a leapfrog pattern, but geographic processes may have been an equally important driver. We suspect that strictly adaptive radiations may be rare in plants and that most events of rapid clade diversification may have involved a mixture of geographical and ecological divergence.

Six Impossible Things before Breakfast: Assumptions, Models, and Belief in Molecular Dating

Confidence in molecular dating analyses has grown with the increasing sophistication of the methods. Some problematic cases where molecular dates disagreed with paleontological estimates appear to have been resolved with a growing agreement between molecules and fossils. But we cannot relax just yet. The growing analytical sophistication of many molecular dating methods relies on an increasingly large number of assumptions about evolutionary history and processes. Many of these assumptions are based on statistical tractability rather than being informed by improved understanding of molecular evolution, yet changing the assumptions can influence molecular dates. How can we tell if the answers we get are driven more by the assumptions we make than by the molecular data being analyzed?

Data partitioning and correction for ascertainment bias reduce the uncertainty of placental mammal divergence times inferred from the morphological clock

Bayesian estimates of divergence times based on the molecular clock yield uncertainty of parameter estimates measured by the width of posterior distributions of node ages. For the relaxed molecular clock, previous works have reported that some of the uncertainty inherent to the variation of rates among lineages may be reduced by partitioning data. Here we test this effect for the purely morphological clock, using placental mammals as a case study. We applied the uncorrelated lognormal relaxed clock to morphological data of 40 extant mammalian taxa and 4,533 characters, taken from the largest published matrix of discrete phenotypic characters. The morphologically derived timescale was compared to divergence times inferred from molecular and combined data. We show that partitioning data into anatomical units significantly reduced the uncertainty of divergence time estimates for morphological data. For the first time, we demonstrate that ascertainment bias has an impact on the precision of morphological clock estimates. While analyses including molecular data suggested most divergences between placental orders occurred near the K-Pg boundary, the partitioned morphological clock recovered older interordinal splits and some younger intraordinal ones, including significantly later dates for the radiation of bats and rodents, which accord to the short-fuse hypothesis.

Estimating divergence times and ancestral breeding systems in Ficus and Moraceae

Background and Aims

Although dioecy, which characterizes only 6 % of angiosperm species, has been considered an evolutionary dead end, recent studies have demonstrated that this is not necessarily the case. Moraceae (40 genera, 1100 spp., including Ficus, 750 spp.) are particularly diverse in breeding systems (including monoecy, gynodioecy, androdioecy and dioecy) and thus represent a model clade to study macroevolution of dioecy.

Methods

Ancestral breeding systems of Ficus and Moraceae were inferred. To do so, a new dated phylogenetic tree of Ficus and Moraceae was first reconstructed by combining a revised 12-fossil calibration set and a densely sampled molecular data set of eight markers and 320 species. Breeding system evolution was then reconstructed using both parsimony and model-based (maximum likelihood and Bayesian) approaches with this new time scale.

Key Results

The crown group ages of Ficus and Moraceae were estimated in the Eocene (40.6-55.9 Ma) and Late Cretaceous (73.2-84.7 Ma), respectively. Strong support was found for ancestral dioecy in Moraceae. Although the ancestral state of Ficus remained particularly sensitive to model selection, the results show that monoecy and gynodioecy evolved from dioecy in Moraceae, and suggest that gynodioecy probably evolved from monoecy in Ficus.

Conclusions

Dioecy was found not to be an evolutionary dead end in Moraceae. This study provides a new time scale for the phylogeny and a new framework of breeding system evolution in Ficus and Moraceae.

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

BACKGROUND AND AIMS

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

METHODS

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

RESULTS

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

CONCLUSIONS

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

Accounting for Calibration Uncertainty: Bayesian Molecular Dating as a "Doubly Intractable" Problem

This study introduces a new Bayesian technique for molecular dating that explicitly accommodates for uncertainty in the phylogenetic position of calibrated nodes derived from the analysis of fossil data. The proposed approach thus defines an adequate framework for incorporating expert knowledge and/or prior information about the way fossils were collected in the inference of node ages. Although it belongs to the class of "node-dating" approaches, this method shares interesting properties with "tip-dating" techniques. Yet, it alleviates some of the computational and modeling difficulties that hamper tip-dating approaches. The influence of fossil data on the probabilistic distribution of trees is the crux of the matter considered here. More specifically, among all the phylogenies that a tree model (e.g., the birth-death process) generates, only a fraction of them "agree" with the fossil data. Bayesian inference under the new model requires taking this fraction into account. However, evaluating this quantity is difficult in practice. A generic solution to this issue is presented here. The proposed approach relies on a recent statistical technique, the so-called exchange algorithm, dedicated to drawing samples from "doubly intractable" distributions. A small example illustrates the problem of interest and the impact of uncertainty in the placement of calibration constraints in the phylogeny given fossil data. An analysis of land plant sequences and multiple fossils further highlights the pertinence of the proposed approach.

Climatic Changes and Orogeneses in the Late Miocene of Eurasia: The Main Triggers of an Expansion at a Continental Scale?

Migrations from the Qinghai-Tibetan Plateau (QTP) to other temperate regions represent one of the main biogeographical patterns for the Northern Hemisphere. However, the ages and routes of these migrations are largely not known. We aimed to reconstruct a well-resolved and dated phylogeny of Hippophae L. (Elaeagnaceae) and test hypothesis of a westward migration of this plant out of the QTP across Eurasian mountains in the Miocene. We produced two data matrices of five chloroplast DNA (cpDNA) and five nuclear DNA markers for all distinct taxa of Hippophae. These matrices were used to reconstruct phylogenetic relationships in the genus. In dating analyses, we first estimated the stem node age of Elaeagnaceae using five fossil records evenly distributed across a tree of Rosales. We used this estimate and two fossil records to calibrate the cpDNA and nDNA phylogenies of Hippophae. The same phylogenies were used to reconstruct ancestral areas within the genus. The monophyly of Hippophae, all five species, and most of subspecies was strongly supported by both plastid and nuclear data sets. Diversification of Hippophae likely started in central Himalayas/southern Tibet in the early Miocene and all extant distinct species had probably originated by the middle Miocene. Diversification of Hippophae rhamnoides likely started in the late Miocene east of the QTP from where this species rapidly expanded to central and western Eurasia. Our findings highlight the impact of different stages in uplift of the QTP and Eurasian mountains and climatic changes in the Neogene on diversification and range shifts in the highland flora on the continent. The results provide support to the idea of an immigration route for some European highland plants from their ancestral areas on the QTP across central and western mountain ranges of Eurasia in the late Miocene.

Rapid Intraspecific Diversification of the Alpine Species Saxifraga sinomontana (Saxifragaceae) in the Qinghai-Tibetan Plateau and Himalayas

An increasing number of phylogeographic studies have been conducted for plant species in the Qinghai-Tibetan Plateau (QTP) and its flanking mountains. However, these studies have mainly focused on the determination of glacial refugia and routes of inter-/post-glacial expansions. Rapid intraspecific diversification of plants in this region have not been thoroughly discussed. Herein, we investigate the effects of the Quaternary climate changes on population genetic structure and diversifications of a herbaceous alpine species, Saxifraga sinomontana, which may have an evolutionary time scale <5 million years in the QTP and Himalayan regions. Using a total of 350 individuals from 29 populations, we studied the evolutionary history of S. sinomontana by analyzing cpDNA trnL-trnF, rpl16 and nrDNA ITS sequences. A total of 89 haplotypes and 158 genotypes were detected for cpDNA and ITS sequences, respectively. Only a few haplotypes/genotypes were widespread, while an extremely large number of haplotypes/genotypes were restricted to single populations, which were scattered throughout the current geographical range of S. sinomontana. This suggests the existence of microrefugia of this species during the Quaternary glaciations. In addition, the relationships of the haplotypes/genotypes were almost completely not resolved by phylogenetic reconstruction. Combining characteristics in terms of high haplotype richness, large proportion of private haplotypes, and shallow haplotype divergence, we speculate that recent intraspecific diversification has occurred in S. sinomontana. Molecular clock analysis estimated that the onset diversification within S. sinomontana to be 1.09 Ma (95% HPD = 0.80–1.45), coinciding with the extensive Quaternary glaciations on the QTP which started ca. 1.17 Ma. The Quaternary climatic oscillations may have triggered rapid intraspecific diversification in this QTP-Himalayan species. However, large niche breadth, as well as introgression/hybridization between the studied species and its closely related sympatric saxifrages, may also played a role to some extent on the current genetic structure of S. sinomontana, which need to be further studied.

An overview of extant conifer evolution from the perspective of the fossil record

PREMISE OF THE STUDY

Conifers are an important living seed plant lineage with an extensive fossil record spanning more than 300 million years. The group therefore provides an excellent opportunity to explore congruence and conflict between dated molecular phylogenies and the fossil record.

METHODS

We surveyed the current state of knowledge in conifer phylogenetics to present a new time-calibrated molecular tree that samples ~90% of extant species diversity. We compared phylogenetic relationships and estimated divergence ages in this new phylogeny with the paleobotanical record, focusing on clades that are species-rich and well known from fossils.

KEY RESULTS

Molecular topologies and estimated divergence ages largely agree with the fossil record in Cupressaceae, conflict with it in Araucariaceae, and are ambiguous in Pinaceae and Podocarpaceae. Molecular phylogenies provide insights into some fundamental questions in conifer evolution, such as the origin of their seed cones, but using them to reconstruct the evolutionary history of specific traits can be challenging.

CONCLUSIONS

Molecular phylogenies are useful for answering deep questions in conifer evolution if they depend on understanding relationships among extant lineages. Because of extinction, however, molecular datasets poorly sample diversity from periods much earlier than the Late Cretaceous. This fundamentally limits their utility for understanding deep patterns of character evolution and resolving the overall pattern of conifer phylogeny.

Key questions and challenges in angiosperm macroevolution

Contents Summary 1170 I. Introduction 1170 II. Six key questions 1172 III. Three key challenges 1177 IV. Conclusions 1181 Acknowledgements 1182 References 1183 SUMMARY: The origin and rapid diversification of angiosperms (flowering plants) represent one of the most intriguing topics in evolutionary biology. Despite considerable progress made in complementary fields over the last two decades (paleobotany, phylogenetics, ecology, evo-devo, genomics), many important questions remain. For instance, what has been the impact of mass extinctions on angiosperm diversification? Are the angiosperms an adaptive radiation? Has morphological evolution in angiosperms been gradual or pulsed? We propose that the recent and ongoing revolution in macroevolutionary methods provides an unprecedented opportunity to explore long-standing questions that probably hold important clues to understand present-day biodiversity. We present six key questions that explore the origin and diversification of angiosperms. We also identify three key challenges to address these questions: (1) the development of new integrative models that include diversification, multiple intrinsic and environmental traits, biogeography and the fossil record all at once, whilst accounting for sampling bias and heterogeneity of macroevolutionary processes through time and among lineages; (2) the need for large and standardized synthetic databases of morphological variation; and (3) continuous effort on sampling the fossil record, but with a revolution in current paleobotanical practice.

Evolution of Climatic Related Leaf Traits in the Family Nothofagaceae

The current relationship between leaf traits and environmental variables has been widely used as a proxy for climate estimates. However, it has been observed that the phylogenetic relationships between taxa also influence the evolution of climatic related leaf traits, implying that the direct use of the physiognomy-climate relation should be corrected by their ancestor-descendant relations. Here, we analyze the variation of 20 leaf traits during the evolution of 27 species in the Gondwana family Nothofagaceae. We evaluate whether the evolution of these traits is exclusively associated with past climate variations or whether they are restricted by phylogenetic relationships. Our results indicate that four leaf traits, associated with size and shape, had consistently a phylogenetic independent evolution, suggesting adaptive variation with the environment. While three of the traits, presented consistently phylogenetic signal and fit a Brownian motion or Ornstein-Uhlenbeck model of evolution, suggesting that the evolution of these traits is restrained by phylogenetic relationships and implying that phylogenetic corrections should be made for the family Nothofagaceae to use them as climatic proxy. Finally, this study highlights the importance of evaluating the evolutionary history of climatic related leaf traits before conducting paleoclimate estimates.

Fossils know it best: Using a new set of fossil calibrations to improve the temporal phylogenetic framework of murid rodents (Rodentia: Muridae)

Murid rodents (Rodentia: Muridae) represent the most diverse and abundant mammalian family. In this study, we provide a refined set of fossil calibrations which is used to reconstruct a dated phylogeny of the family using a multilocus dataset (six nuclear and nine mitochondrial gene fragments) encompassing 161 species representing 82 murid genera from four extant subfamilies (Deomyinae, Gerbillinae, Lophiomyinae and Murinae). In comparison with previous studies on murid or muroid rodents, our work stands out for the implementation of nine robust fossil constraints within the Muridae thanks to a thorough review of the fossil record. Before being assigned to specific nodes of the phylogeny, all potential fossil constraints were carefully assessed; they were also subjected to several cross-validation analyses. The resulting phylogeny is consistent with previous phylogenetic studies on murids, and recovers the monophyly of all sampled murid subfamilies and tribes. Based on nine controlled fossil calibrations, our inferred temporal timeframe indicates that the murid family likely originated in the course of the Early Miocene, 22.0-17.0 million years ago (Ma), and that most major lineages (i.e. tribes) started diversifying ca. 10 Ma. Historical biogeography analyses support the tropical origin for the family, with an initial internal split (vicariance event) between Afrotropical and Oriental (Indomalaya and Philippines) lineages. During the course of their diversification, the biogeographic pattern of murids is marked by several dispersal events toward the Australasian and the Palearctic regions. The Afrotropical region was also secondarily colonized at least three times from the Indomalaya, indicating that the latter region has acted as a major centre of diversification for the family.

Strategies for Partitioning Clock Models in Phylogenomic Dating: Application to the Angiosperm Evolutionary Timescale

Evolutionary timescales can be inferred from molecular sequence data using a Bayesian phylogenetic approach. In these methods, the molecular clock is often calibrated using fossil data. The uncertainty in these fossil calibrations is important because it determines the limiting posterior distribution for divergence-time estimates as the sequence length tends to infinity. Here, we investigate how the accuracy and precision of Bayesian divergence-time estimates improve with the increased clock-partitioning of genome-scale data into clock-subsets. We focus on a data set comprising plastome-scale sequences of 52 angiosperm taxa. There was little difference among the Bayesian date estimates whether we chose clock-subsets based on patterns of among-lineage rate heterogeneity or relative rates across genes, or by random assignment. Increasing the degree of clock-partitioning usually led to an improvement in the precision of divergence-time estimates, but this increase was asymptotic to a limit presumably imposed by fossil calibrations. Our clock-partitioning approaches yielded highly precise age estimates for several key nodes in the angiosperm phylogeny. For example, when partitioning the data into 20 clock-subsets based on patterns of among-lineage rate heterogeneity, we inferred crown angiosperms to have arisen 198–178 Ma. This demonstrates that judicious clock-partitioning can improve the precision of molecular dating based on phylogenomic data, but the meaning of this increased precision should be considered critically.

Plastome sequences and exploration of tree-space help to resolve the phylogeny of riceflowers (Thymelaeaceae: Pimelea)

Data sets comprising small numbers of genetic markers are not always able to resolve phylogenetic relationships. This has frequently been the case in molecular systematic studies of plants, with many analyses being based on sequence data from only two or three chloroplast genes. An example of this comes from the riceflowers Pimelea Banks & Sol. ex Gaertn. (Thymelaeaceae), a large genus of flowering plants predominantly distributed in Australia. Despite the considerable morphological variation in the genus, low sequence divergence in chloroplast markers has led to the phylogeny of Pimelea remaining largely uncertain. In this study, we resolve the backbone of the phylogeny of Pimelea in comprehensive Bayesian and maximum-likelihood analyses of plastome sequences from 41 taxa. However, some relationships received only moderate to poor support, and the Pimelea clade contained extremely short internal branches. By using topology-clustering analyses, we demonstrate that conflicting phylogenetic signals can be found across the trees estimated from individual chloroplast protein-coding genes. A relaxed-clock dating analysis reveals that Pimelea arose in the mid-Miocene, with most divergences within the genus occurring during a subsequent rapid diversification. Our new phylogenetic estimate offers better resolution and is more strongly supported than previous estimates, providing a platform for future taxonomic revisions of both Pimelea and the broader subfamily. Our study has demonstrated the substantial improvements in phylogenetic resolution that can be achieved using plastome-scale data sets in plant molecular systematics.

Constraining uncertainty in the timescale of angiosperm evolution and the veracity of a Cretaceous Terrestrial Revolution

Through the lens of the fossil record, angiosperm diversification precipitated a Cretaceous Terrestrial Revolution (KTR) in which pollinators, herbivores and predators underwent explosive co-diversification. Molecular dating studies imply that early angiosperm evolution is not documented in the fossil record. This mismatch remains controversial. We used a Bayesian molecular dating method to analyse a dataset of 83 genes from 644 taxa and 52 fossil calibrations to explore the effect of different interpretations of the fossil record, molecular clock models, data partitioning, among other factors, on angiosperm divergence time estimation. Controlling for different sources of uncertainty indicates that the timescale of angiosperm diversification is much less certain than previous molecular dating studies have suggested. Discord between molecular clock and purely fossil-based interpretations of angiosperm diversification may be a consequence of false precision on both sides. We reject a post-Jurassic origin of angiosperms, supporting the notion of a cryptic early history of angiosperms, but this history may be as much as 121 Myr, or as little as 23 Myr. These conclusions remain compatible with palaeobotanical evidence and a more general KTR in which major groups of angiosperms diverged later within the Cretaceous, alongside the diversification of pollinators, herbivores and their predators.

The Past Sure is Tense: On Interpreting Phylogenetic Divergence Time Estimates

Divergence time estimation-the calibration of a phylogeny to geological time-is an integral first step in modeling the tempo of biological evolution (traits and lineages). However, despite increasingly sophisticated methods to infer divergence times from molecular genetic sequences, the estimated age of many nodes across the tree of life contrast significantly and consistently with timeframes conveyed by the fossil record. This is perhaps best exemplified by crown angiosperms, where molecular clock (Triassic) estimates predate the oldest (Early Cretaceous) undisputed angiosperm fossils by tens of millions of years or more. While the incompleteness of the fossil record is a common concern, issues of data limitation and model inadequacy are viable (if underexplored) alternative explanations. In this vein, Beaulieu et al. (2015) convincingly demonstrated how methods of divergence time inference can be misled by both (i) extreme state-dependent molecular substitution rate heterogeneity and (ii) biased sampling of representative major lineages. These results demonstrate the impact of (potentially common) model violations. Here, we suggest another potential challenge: that the configuration of the statistical inference problem (i.e., the parameters, their relationships, and associated priors) alone may preclude the reconstruction of the paleontological timeframe for the crown age of angiosperms. We demonstrate, through sampling from the joint prior (formed by combining the tree (diversification) prior with the calibration densities specified for fossil-calibrated nodes) that with no data present at all, that an Early Cretaceous crown angiosperms is rejected (i.e., has essentially zero probability). More worrisome, however, is that for the 24 nodes calibrated by fossils, almost all have indistinguishable marginal prior and posterior age distributions when employing routine lognormal fossil calibration priors. These results indicate that there is inadequate information in the data to over-rule the joint prior. Given that these calibrated nodes are strategically placed in disparate regions of the tree, they act to anchor the tree scaffold, and so the posterior inference for the tree as a whole is largely determined by the pseudodata present in the (often arbitrary) calibration densities. We recommend, as for any Bayesian analysis, that marginal prior and posterior distributions be carefully compared to determine whether signal is coming from the data or prior belief, especially for parameters of direct interest. This recommendation is not novel. However, given how rarely such checks are carried out in evolutionary biology, it bears repeating. Our results demonstrate the fundamental importance of prior/posterior comparisons in any Bayesian analysis, and we hope that they further encourage both researchers and journals to consistently adopt this crucial step as standard practice. Finally, we note that the results presented here do not refute the biological modeling concerns identified by Beaulieu et al. (2015). Both sets of issues remain apposite to the goals of accurate divergence time estimation, and only by considering them in tandem can we move forward more confidently.

Phylogeographic analysis reveals high genetic structure with uniform phenotypes in the paper wasp Protonectarina sylveirae (Hymenoptera: Vespidae)

Swarm-founding wasps are endemic and common representatives of neotropical fauna and compose an interesting social tribe of vespids, presenting both complex social characteristics and uncommon traits for a eusocial group, such as the absence of castes with distinct morphology. The paper wasp Protonectarina sylveirae (Saussure) presents a broad distribution from Brazil, Argentina and Paraguay, occurring widespread in the Atlantic rainforest and arboreal Caatinga, being absent in the Amazon region. Given the peculiar distribution among swarm-founding wasps, an integrative approach to reconstruct the evolutionary history of P. sylveirae in a spatial-temporal framework was performed to investigate: the presence of genetic structure and its relationship with the geography, the evolution of distinct morphologic lineages and the possible historical event(s) in Neotropical region, which could explain the observed phylogeographic pattern. Individuals of P. sylveirae were obtained from populations of 16 areas throughout its distribution for DNA extraction and amplification of mitochondrial genes 12S, 16S and COI. Analysis of genetic diversity, construction of haplotype net, analysis of population structure and dating analysis of divergence time were performed. A morphometric analysis was also performed using 8 measures of the body of the adult (workers) to test if there are morphological distinction among populations. Thirty-five haplotypes were identified, most of them exclusively of a group and a high population structure was found. The possibility of genetic divergence because of isolation by distance was rejected. Morphological analysis pointed to a great uniformity in phenotypes, with only a small degree of differentiation between populations of south and the remaining. Divergence time analysis showed a Middle/Late Miocene origin, a period where an extensive marine ingression occurred in South America. Divergence of haplogroups began from the Plio/Pleistocene boundary and the last glacial maximum most likely modeled the current distribution of species, even though it was not the cause of genetic breaks.

Evaluating the Impact of Genomic Data and Priors on Bayesian Estimates of the Angiosperm Evolutionary Timescale

The evolutionary timescale of angiosperms has long been a key question in biology. Molecular estimates of this timescale have shown considerable variation, being influenced by differences in taxon sampling, gene sampling, fossil calibrations, evolutionary models, and choices of priors. Here, we analyze a data set comprising 76 protein-coding genes from the chloroplast genomes of 195 taxa spanning 86 families, including novel genome sequences for 11 taxa, to evaluate the impact of models, priors, and gene sampling on Bayesian estimates of the angiosperm evolutionary timescale. Using a Bayesian relaxed molecular-clock method, with a core set of 35 minimum and two maximum fossil constraints, we estimated that crown angiosperms arose 221 (251-192) Ma during the Triassic. Based on a range of additional sensitivity and subsampling analyses, we found that our date estimates were generally robust to large changes in the parameters of the birth-death tree prior and of the model of rate variation across branches. We found an exception to this when we implemented fossil calibrations in the form of highly informative gamma priors rather than as uniform priors on node ages. Under all other calibration schemes, including trials of seven maximum age constraints, we consistently found that the earliest divergences of angiosperm clades substantially predate the oldest fossils that can be assigned unequivocally to their crown group. Overall, our results and experiments with genome-scale data suggest that reliable estimates of the angiosperm crown age will require increased taxon sampling, significant methodological changes, and new information from the fossil record. [Angiospermae, chloroplast, genome, molecular dating, Triassic.].

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.

Toward a Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa

Rapidly growing biological data—including molecular sequences and fossils—hold an unprecedented potential to reveal how evolutionary processes generate and maintain biodiversity. However, researchers often have to develop their own idiosyncratic workflows to integrate and analyze these data for reconstructing time-calibrated phylogenies. In addition, divergence times estimated under different methods and assumptions, and based on data of various quality and reliability, should not be combined without proper correction. Here we introduce a modular framework termed SUPERSMART (Self-Updating Platform for Estimating Rates of Speciation and Migration, Ages, and Relationships of Taxa), and provide a proof of concept for dealing with the moving targets of evolutionary and biogeographical research. This framework assembles comprehensive data sets of molecular and fossil data for any taxa and infers dated phylogenies using robust species tree methods, also allowing for the inclusion of genomic data produced through next-generation sequencing techniques. We exemplify the application of our method by presenting phylogenetic and dating analyses for the mammal order Primates and for the plant family Arecaceae (palms). We believe that this framework will provide a valuable tool for a wide range of hypothesis-driven research questions in systematics, biogeography, and evolution. SUPERSMART will also accelerate the inference of a “Dated Tree of Life” where all node ages are directly comparable.

A continuous morphological approach to study the evolution of pollen in a phylogenetic context: An example with the order Myrtales

The study of pollen morphology has historically allowed evolutionary biologists to assess phylogenetic relationships among Angiosperms, as well as to better understand the fossil record. During this process, pollen has mainly been studied by discretizing some of its main characteristics such as size, shape, and exine ornamentation. One large plant clade in which pollen has been used this way for phylogenetic inference and character mapping is the order Myrtales, composed by the small families Alzateaceae, Crypteroniaceae, and Penaeaceae (collectively the "CAP clade"), as well as the large families Combretaceae, Lythraceae, Melastomataceae, Myrtaceae, Onagraceae and Vochysiaceae. In this study, we present a novel way to study pollen evolution by using quantitative size and shape variables. We use morphometric and morphospace methods to evaluate pollen change in the order Myrtales using a time-calibrated, supermatrix phylogeny. We then test for conservatism, divergence, and morphological convergence of pollen and for correlation between the latitudinal gradient and pollen size and shape. To obtain an estimate of shape, Myrtales pollen images were extracted from the literature, and their outlines analyzed using elliptic Fourier methods. Shape and size variables were then analyzed in a phylogenetic framework under an Ornstein-Uhlenbeck process to test for shifts in size and shape during the evolutionary history of Myrtales. Few shifts in Myrtales pollen morphology were found which indicates morphological conservatism. Heterocolpate, small pollen is ancestral with largest pollen in Onagraceae. Convergent shifts in shape but not size occurred in Myrtaceae and Onagraceae and are correlated to shifts in latitude and biogeography. A quantitative approach was applied for the first time to examine pollen evolution across a large time scale. Using phylogenetic based morphometrics and an OU process, hypotheses of pollen size and shape were tested across Myrtales. Convergent pollen shifts and position in the latitudinal gradient support the selective role of harmomegathy, the mechanism by which pollen grains accommodate their volume in response to water loss.

Contrasting evolutionary processes during Quaternary climatic changes and historical orogenies: a case study of the Japanese endemic primroses Primula sect. Reinii

BACKGROUND AND AIMS

Recent biogeographic studies have illustrated that Quaternary climatic changes and historical orogenies have driven the development of high levels of biodiversity. In this context, phylogenetic niche conservatism may play a role as a major precursor of allopatric speciation. However, the effects of niche evolution on the diversification patterns of plant species under rapid habitat changes are still unknown. Here, Primula section Reinii, one of the few primroses endemic to the Japanese Archipelago, was investigated. This study aimed to clarify the phylogenetic position and relationships of section Reinii, interpret the biogeographic and diversification patterns of this group and gain a better understanding of the role of climatic niche evolution in the Japanese endemic primroses.

METHODS

Dated phylogeny for Primula section Reinii is presented based on the sequences of six chloroplast genes and one nuclear gene. Biogeographic history was reconstructed using statistical dispersal-vicariance analysis. Macro-evolutionary modelling of the climatic niche was combined with biogeographic inferences.

KEY RESULTS

Section Reinii was shown to be monophyletic based on chloroplast and nuclear sequences. Fossil-calibrated dating analysis estimated that this section diverged from its sister taxon, sect. Cortusoides, around 1.82 million years ago, and intraspecific diversification occurred within the last million years. This time frame was characterized by Quaternary climatic oscillations and uplift of high mountains in Japan. Biogeographic inference suggested that this section originated at the northern end of the Japanese Archipelago and then dispersed southward to other islands. Models of climatic niche evolution indicated that the closely related species P. reinii and P. tosanensis have contrasting niche optima and rates of niche evolution.

CONCLUSION

Our results highlight that spatiotemporal heterogeneity in the Japanese islands may play a significant role in the biogeographic history of Japanese endemic primroses. Contrasting evolutionary processes found in closely related species illustrate the effects of climatic niche evolution on species' diversification patterns.

Both temperature fluctuations and East Asian monsoons have driven plant diversification in the karst ecosystems from southern China

Karst ecosystems in southern China are species-rich and have high levels of endemism, yet little is known regarding the evolutionary processes responsible for the origin and diversification of karst biodiversity. The genus Primulina (Gesneriaceae) comprises ca. 170 species endemic to southern China with high levels of ecological (edaphic) specialization, providing an exceptional model to study the plant diversification in karsts. We used molecular data from nine chloroplast and 11 nuclear regions and macroevolutionary analyses to assess the origin and cause of species diversification due to palaeoenvironmental changes and edaphic specialization in Primulina. We found that speciation was positively associated with changes in past temperatures and East Asian monsoons through the evolutionary history of Primulina. Climatic change around the mid-Miocene triggered an early burst followed by a slowdown of diversification rate towards the present with the climate cooling. We detected different speciation rates among edaphic types, and transitions among soil types were infrequently and did not impact the overall speciation rate. Our findings suggest that both global temperature changes and East Asian monsoons have played crucial roles in floristic diversification within the karst ecosystems in southern China, such that speciation was higher when climate was warmer and wetter. This is the first study to directly demonstrate that past monsoon activity is positively correlated with speciation rate in East Asia. This case study could motivate further investigations to assess the impacts of past environmental changes on the origin and diversification of biodiversity in global karst ecosystems, most of which are under threat.

Genomic evidence reveals a radiation of placental mammals uninterrupted by the KPg boundary

The timing of the diversification of placental mammals relative to the Cretaceous-Paleogene (KPg) boundary mass extinction remains highly controversial. In particular, there have been seemingly irreconcilable differences in the dating of the early placental radiation not only between fossil-based and molecular datasets but also among molecular datasets. To help resolve this discrepancy, we performed genome-scale analyses using 4,388 loci from 90 taxa, including representatives of all extant placental orders and transcriptome data from flying lemurs (Dermoptera) and pangolins (Pholidota). Depending on the gene partitioning scheme, molecular clock model, and genic deviation from molecular clock assumptions, extensive sensitivity analyses recovered widely varying diversification scenarios for placental mammals from a given gene set, ranging from a deep Cretaceous origin and diversification to a scenario spanning the KPg boundary, suggesting that the use of suboptimal molecular clock markers and methodologies is a major cause of controversies regarding placental diversification timing. We demonstrate that reconciliation between molecular and paleontological estimates of placental divergence times can be achieved using the appropriate clock model and gene partitioning scheme while accounting for the degree to which individual genes violate molecular clock assumptions. A birth-death-shift analysis suggests that placental mammals underwent a continuous radiation across the KPg boundary without apparent interruption by the mass extinction, paralleling a genus-level radiation of multituberculates and ecomorphological diversification of both multituberculates and therians. These findings suggest that the KPg catastrophe evidently played a limited role in placental diversification, which, instead, was likely a delayed response to the slightly earlier radiation of angiosperms.

Evolutionary priority effects persist in anthropogenically created habitats, but not through nonnative plant invasion

Evolutionary priority effects, where early-arriving lineages occupy niche space via diversification and preclude dominance of later arrivals, have been observed in alpine and forest communities. However, the potential for evolutionary priority effects to persist in an era of rapid global change remains unclear. Here, we use a natural experiment of historical disturbance in New Zealand to test whether anthropogenic changes in available habitat and nonnative invasion eliminate the role of evolutionary priority effects in community assembly. We also test whether evolutionary priority effects diminish with decreasing resource availability. Older plant clades, as estimated by clade crown age, were relatively more abundant in both primary and secondary grassland. Relative abundance in primary grassland decreased with clade stem age, but only weakly. However, for both clade age estimates, relative abundance decreased with age when nonnative biomass was high and soil moisture was low. Our data show that patterns in community structure consistent with evolutionary priority effects can occur in both primary and secondary grasslands, the latter created by anthropogenic disturbance. However, nonnative invasion may overwhelm the effect of immigration timing on community dominance, possibly as a result of high immigration rates and preadaptation to anthropogenically modified environments.

Updating the Phylogenetic Dating of New Caledonian Biodiversity with a Meta-analysis of the Available Evidence

For a long time, New Caledonia was considered a continental island, a fragment of Gondwana harbouring old clades that originated by vicariance and so were thought to be locally ancient. Recent molecular phylogenetic studies dating diversification and geological data indicating important events of submergence during the Paleocene and Eocene (until 37 Ma) brought evidence to dismiss this old hypothesis. In spite of this, some authors still insist on the idea of a local permanence of a Gondwanan biota, justifying this assumption through a complex scenario of survival by hopping to and from nearby and now-vanished islands. Based on a comprehensive review of the literature, we found 40 studies dating regional clades of diverse organisms and we used them to test the hypothesis that New Caledonian and inclusive Pacific island clades are older than 37 Ma. The results of this meta-analysis provide strong evidence for refuting the hypothesis of a Gondwanan refuge with a biota that originated by vicariance. Only a few inclusive Pacific clades (6 out of 40) were older than the oldest existing island. We suggest that these clades could have extinct members either on vanished islands or nearby continents, emphasizing the role of dispersal and extinction in shaping the present-day biota.

Plant phylogeny as a window on the evolution of hyperdiversity in the tropical rainforest biome

I. II. III. IV. V. VI. VII. VIII. IX. References SUMMARY: Tropical rainforest (TRF) is the most species-rich terrestrial biome on Earth, harbouring just under half of the world's plant species in c. 7% of the land surface. Phylogenetic trees provide important insights into mechanisms underpinning TRF hyperdiversity that are complementary to those obtained from the fossil record. Phylogenetic studies of TRF plant diversity have mainly focused on whether this biome is an evolutionary 'cradle' or 'museum', emphasizing speciation and extinction rates. However, other explanations, such as biome age, immigration and ecological limits, must also be considered. We present a conceptual framework for addressing the drivers of TRF diversity, and review plant studies that have tested them with phylogenetic data. Although surprisingly few in number, these studies point to old age of TRF, low extinction and high speciation rates as credible drivers of TRF hyperdiversity. There is less evidence for immigration and ecological limits, but these cannot be dismissed owing to the limited number of studies. Rapid methodological developments in DNA sequencing, macroevolutionary analysis and the integration of phylogenetics with other disciplines may improve our grasp of TRF hyperdiversity in the future. However, such advances are critically dependent on fundamental systematic research, yielding numerous, additional, well-sampled phylogenies of TRF lineages.

Another look at the phylogenetic relationships and intercontinental biogeography of eastern Asian - North American Rhus gall aphids (Hemiptera: Aphididae: Eriosomatinae): Evidence from mitogenome sequences via genome skimming

The Rhus gall aphids are sometimes referred to as subtribe Melaphidina (Aphididae: Eriosomatinae: Fordini) and comprise a unique group that forms galls on the primary host plants, Rhus. We examined the evolutionary relationships within the Melaphidina aphids using sequences of the complete mitochondrial genome and with samples of 11 of the 12 recognized species representing all six genera. Bayesian, maximum likelihood and parsimony analyses of the mitochondrial genome data support five well-supported clades within Melaphidina: (1) Nurudea (except N. ibofushi), (2) Schlechtendalia-Nurudea ibofushi, (3) Meitanaphis-Kaburagia, (4) Floraphis, and (5) Melaphis. Nurudea shiraii and N. yanoniella are sister to each other, but N. ibofushi is nested within Schlechtendalia. The Nurudea shiraii-N. yanoniella clade is sister to the large clade of the remaining taxa of Melaphidina aphids. The Bayesian and maximum likelihood analyses support the North American Melaphis rhois as sister to the clade of Floraphis-Kaburagia-Meitanaphis-Schlechtendalia from eastern Asia, whereas the parsimony analysis suggests Melaphis sister to Floraphis with low support (bootstrap support 38%), and the amino acid data weakly place it sister to Schlechtendalia-Nurudea ibofushi. The Melaphis position needs to be further tested with nuclear data. Meitanaphis flavogallis is sister to Kaburagia species instead of grouping with Meitanaphis elongallis. Using the Bayesian method, the North American Melaphis was estimated to have diverged from its closest Asian relatives around 64.6 (95% HPD 59.4-69.8) Ma, which is in the early Paleocene near the Cretaceous and Paleogene boundary (K/Pg boundary). At the K/Pg boundary, mass extinctions caused many types of insect-plant associations to disappear, and these extinctions may explain some of the difficulties in the phylogenetic placement of Melaphis within the analyses.