The phylogeny and biogeography of Hakea (Proteaceae) reveals the role of biome shifts in a continental plant radiation

A Phylogenetic Model of Established and Enabled Biome Shifts

Where each species actually lives is distinct from where it could potentially survive and persist. This suggests that it may be important to distinguish established from enabled biome affinities when considering how ancestral species moved and evolved among major habitat types. We introduce a new phylogenetic method, called RFBS, to model how anagenetic and cladogenetic events cause established and enabled biome affinities (or, more generally, other discrete realized versus fundamental niche states) to shift over evolutionary timescale. We provide practical guidelines for how to assign established and enabled biome affinity states to extant taxa, using the flowering plant clade Viburnum as a case study. Through a battery of simulation experiments, we show that RFBS performs well, even when we have realistically imperfect knowledge of enabled biome affinities for most analyzed species. We also show that RFBS reliably discerns established from enabled affinities, with similar accuracy to standard competing models that ignore the existence of enabled biome affinities. Lastly, we apply RFBS to Viburnum to infer ancestral biomes throughout the tree and to highlight instances where repeated shifts between established affinities for warm and cold temperate forest biomes were enabled by a stable and slowly-evolving enabled affinity for both temperate biomes.

Transition of survival strategies under global climate shifts in the grape family

Faced with environmental changes, plants may either move to track their ancestral niches or evolve to adapt to new niches. Vitaceae, the grape family, has evolved diverse adaptive traits facilitating a global expansion in wide-ranging habitats, making it ideal for investigating transition between move and evolve strategies and exploring the underlying mechanisms. Here we inferred the patterns of biogeographic diversification and trait evolution in Vitaceae based on a robust phylogeny with dense sampling including 495 species (~52% of Vitaceae species). Vitaceae probably originated from Asia-the diversity centre of extant genera and the major source of dispersals. Boundaries of the Eocene, Oligocene and Miocene were identified as turning points in shifting strategies. A significant decrease in move strategy was identified during the Oligocene, followed by increases in move and evolve. After the Miocene, evolve began to dominate, during which increased niche opportunities and key trait innovations played important roles.

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

Invasiveness, Monitoring and Control of Hakea sericea: A Systematic Review

Solutions for ecological and economic problems posed by Hakea sericea invasions rely on scientific knowledge. We conducted a systematic review to analyze and synthesize the past and current scientific knowledge concerning H. sericea invasion processes and mechanisms, as well as monitoring and control techniques. We used ISI Web of Science, Scopus, and CAPES Periodicals to look for publications on the ecological and environmental factors involved in H. sericea establishment (question 1); responses of H. sericea to fire in native and invaded ecosystems (question 2); and H. sericea monitoring and control methods (question 3). We identified 207 publications, 47.4% of which related to question 1, mainly from Australia and South Africa, with an increasing trend in the number of publications on monitoring and modeling. The traits identified in our systematic review, such as adaptations to dystrophic environments, drought resistance, sclerophylly, low transpiration rates, high nutrient use efficiency, stomatal conductance and photosynthetic rates, strong serotiny, proteoid roots and high post-fire seed survival and seedling recruitment, highlighted that H. sericea is a successful invader species due to its long adaptive history mediated by an arsenal of ecophysiological mechanisms that place it at a superior competitive level, especially in fire-prone ecosystems. Integrated cost-effective control methods in selected areas and the incorporation of information on the temporal invasion dynamics can significantly improve invasion control and mitigate H. sericea impacts while maintaining the supply of ecosystem services in invaded areas.

Capturing Genetic Diversity in Seed Collections: An Empirical Study of Two Congeners with Contrasting Mating Systems

Plant mating systems shape patterns of genetic diversity and impact the long-term success of populations. As such, they are relevant to the design of seed collections aiming to maximise genetic diversity (e.g., germplasm conservation, ecological restoration). However, for most species, little is known empirically about how variation in mating systems and genetic diversity is distributed. We investigated the relationship between genetic diversity and mating systems in two functionally similar, co-occurring species of Hakea (Proteaceae), and evaluated the extent to which genetic diversity was captured in seeds. We genotyped hundreds of seedlings and mother plants via DArTseq, and developed novel implementations of two approaches to inferring the mating system from SNP data. A striking contrast in patterns of genetic diversity between H. sericea and H. teretifolia was revealed, consistent with a contrast in their mating systems. While both species had mixed mating systems, H. sericea was found to be habitually selfing, while H. teretifolia more evenly employed both selfing and outcrossing. In both species, seed collection schemes maximised genetic diversity by increasing the number of maternal lines and sites sampled, but twice as many sites were needed for the selfing species to capture equivalent levels of genetic variation at a regional scale.

Phylogenetic biome conservatism as a key concept for an integrative understanding of evolutionary history: Galliformes and Falconiformes as study cases

Biomes are climatically and biotically distinctive macroecological units that formed over geological time scales. Their features consolidate them as ‘evolutionary scenarios’, with their own diversification dynamics. Under the concept of phylogenetic niche conservatism, we assessed, for the first time, the evolution of biome occupation in birds. We aimed to analyse patterns of adaptation to different climatic regimes and the determinant factors for colonization of emerging biomes by clades from different ancestral biomes. In this work, we reconstructed the biome occupation history of two clades of birds (Galliformes and Falconiformes) under an integrative perspective through a comprehensive review of ecological, phylogenetic, palaeontological and biogeographical evidence. Our findings for both groups are consistent with a scenario of phylogenetic biome conservatism and highlight the importance of changes in climate during the Miocene in the adaptation and evolution of climatic niches. In particular, our results indicate high biome conservatism associated with biomes situated in some of the extremes of the global climate gradient (evergreen tropical rainforest, steppe and tundra) for both bird taxa. Finally, the historical dynamics of tropical seasonal biomes, such as tropical deciduous woodlands and savannas, appear to have played a preponderant role during the diversification processes of these bird lineages.

Pioneering polyploids: the impact of whole-genome duplication on biome shifting in New Zealand Coprosma (Rubiaceae) and Veronica (Plantaginaceae)

The role of whole-genome duplication (WGD) in facilitating shifts into novel biomes remains unknown. Focusing on two diverse woody plant groups in New Zealand, Coprosma (Rubiaceae) and Veronica (Plantaginaceae), we investigate how biome occupancy varies with ploidy level, and test the hypothesis that WGD increases the rate of biome shifting. Ploidy levels and biome occupancy (forest, open and alpine) were determined for indigenous species in both clades. The distribution of low-ploidy (Coprosma: 2x, Veronica: 6x) versus high-ploidy (Coprosma: 4-10x, Veronica: 12-18x) species across biomes was tested statistically. Estimation of the phylogenetic history of biome occupancy and WGD was performed using time-calibrated phylogenies and the R package BioGeoBEARS. Trait-dependent dispersal models were implemented to determine support for an increased rate of biome shifting among high-ploidy lineages. We find support for a greater than random portion of high-ploidy species occupying multiple biomes. We also find strong support for high-ploidy lineages showing a three- to eightfold increase in the rate of biome shifts. These results suggest that WGD promotes ecological expansion into new biomes.

Evolution of crassulacean acid metabolism (CAM) as an escape from ecological niche conservatism in Malagasy Bulbophyllum (Orchidaceae)

Despite growing evidence that niche shifts are more common in flowering plants than previously thought, little is known of whether such shifts are promoted by changes in photosynthetic pathways. Here we combine the most complete phylogeny for epiphytic Malagasy Bulbophyllum orchids (c. 210 spp.) with climatic niche and carbon isotope ratios to infer the group's spatial-temporal history, and the role of strongly expressed crassulacean acid metabolism (CAM) in facilitating niche shifts and diversification. We find that most extant species still retain niche (Central Highland) and photosynthesis (C₃ ) states as present in the single mid-Miocene (c. 12.70 million yr ago (Ma)) ancestor colonizing Madagascar. However, we also infer a major transition to CAM, linked to a late Miocene (c. 7.36 Ma) invasion of species from the sub-humid highland first into the island's humid eastern coastal, and then into the seasonally dry 'Northwest Sambirano' rainforests, yet without significant effect on diversification rates. These findings indicate that CAM in tropical epiphytes may be selectively advantageous even in high rainfall habitats, rather than presenting a mere adaptation to dry environments or epiphytism per se. Overall, our study qualifies CAM as an evolutionary 'gateway' trait that considerably widened the spatial-ecological amplitude of Madagascar's most species-rich orchid genus.

Vulnerability to xylem cavitation of Hakea species (Proteaceae) from a range of biomes and life histories predicted by climatic niche

BACKGROUND AND AIMS

Extreme drought conditions across the globe are impacting biodiversity, with serious implications for the persistence of native species. However, quantitative data on physiological tolerance are not available for diverse flora to inform conservation management. We quantified physiological resistance to cavitation in the diverse Hakea genus (Proteaceae) to test predictions based on climatic origin, life history and functional traits.

METHODS

We sampled terminal branches of replicate plants of 16 species in a common garden. Xylem cavitation was induced in branches under varying water potentials (tension) in a centrifuge, and the tension generating 50 % loss of conductivity (stem P50) was characterized as a metric for cavitation resistance. The same branches were used to estimate plant functional traits, including wood density, specific leaf area and Huber value (sap flow area to leaf area ratio).

KEY RESULTS

There was significant variation in stem P50 among species, which was negatively associated with the species climate origin (rainfall and aridity). Cavitation resistance did not differ among life histories; however, a drought avoidance strategy with terete leaf form and greater Huber value may be important for species to colonize and persist in the arid biome.

CONCLUSIONS

This study highlights climate (rainfall and aridity), rather than life history and functional traits, as the key predictor of variation in cavitation resistance (stem P50). Rainfall for species origin was the best predictor of cavitation resistance, explaining variation in stem P50, which appears to be a major determinant of species distribution. This study also indicates that stem P50 is an adaptive trait, genetically determined, and hence reliable and robust for predicting species vulnerability to climate change. Our findings will contribute to future prediction of species vulnerability to drought and adaptive management under climate change.

Ecological interactions shape the evolution of flower color in communities across a temperate biodiversity hotspot

Processes driving the divergence of floral traits may be integral to the extraordinary richness of flowering plants and the assembly of diverse plant communities. Several models of pollinator-mediated floral evolution have been proposed; floral divergence may (i) be directly involved in driving speciation or may occur after speciation driven by (ii) drift or local adaptation in allopatry or (iii) negative interactions between species in sympatry. Here, we generate predictions for patterns of trait divergence and community assembly expected under these three models, and test these predictions in Hakea (Proteaceae), a diverse genus in the Southwest Australian biodiversity hotspot. We quantified functional richness for two key floral traits (pistil length and flower color), as well as phylogenetic distances between species, across ecological communities, and compared these to patterns generated from null models of community assembly. We also estimated the statistical relationship between rates of trait evolution and lineage diversification across the phylogeny. Patterns of community assembly suggest that flower color, but not floral phenology or morphology, or phylogenetic relatedness, is more divergent in communities than expected. Rates of lineage diversification and flower color evolution were negatively correlated across the phylogeny and rates of flower colour evolution were positively related to branching times. These results support a role for diversity-dependent species interactions driving floral divergence during the Hakea radiation, contributing to the development of the extraordinary species richness of southwest Australia.

The living heart: Climate gradients predict desert mountain endemism

Mountain regions are centers of biodiversity endemism at a global scale but the role of arid-zone mountain ranges in shaping biodiversity patterns is poorly understood. Focusing on three guilds of taxa from a desert upland refugium in Australia, we sought to determine: (a) the relative extent to which climate, terrain or geological substrate predict endemism, and (b) whether patterns of endemism are complimentary across broad taxonomic guilds. We mapped regional endemism for plants, land snails, and vertebrates using combined Species Distribution Models (SDMs) for all endemic taxa (n = 82). We then modelled predictors of endemism using Generalised Additive Models (GAMs) and geology, terrain, and climate variables. We tested for the presence of inter- and intraguild hotspots of endemism. Many individual plant and land snail taxa were tightly linked with geology, corresponding to small distributions. Conversely, most vertebrate taxa were not constrained to specific geological substrates and occurred over larger areas. However, across all three guilds climate was the strongest predictor of regional endemism, particularly for plants wherein discrete hotspots of endemism were buffered from extreme summer temperatures. Land snail and vertebrate endemism peaked in areas with highest precipitation in the driest times of the year. Hotspots of endemism within each guild poorly predicted endemism in other guilds. We found an overarching signal that climatic gradients play a dominant role in the persistence of endemic taxa in an arid-zone mountain range system. An association with higher rainfall and cooler temperatures indicates that continuing trends toward hotter and drier climates may lead to range contractions in this, and potentially other, arid-zone mountain biotas. Contrasting patterns of endemism across guilds highlight the need to couple comprehensive regional planning for the protection of climate refugia, with targeted management of more localized and habitat specialist taxa.

Coordination of stem and leaf traits define different strategies to regulate water loss and tolerance ranges to aridity

Adaptation to drought involves complex interactions of traits that vary within and among species. To date, few data are available to quantify within-species variation in functional traits and they are rarely integrated into mechanistic models to improve predictions of species response to climate change. We quantified intraspecific variation in functional traits of two Hakea species growing along an aridity gradient in southeastern Australia. Measured traits were later used to parameterise the model SurEau to simulate a transplantation experiment to identify the limits of drought tolerance. Embolism resistance varied between species but not across populations. Instead, populations adjusted to drier conditions via contrasting sets of trait trade-offs that facilitated homeostasis of plant water status. The species from relatively mesic climate, Hakea dactyloides, relied on tight stomatal control whereas the species from xeric climate, Hakea leucoptera dramatically increased Huber value and leaf mass per area, while leaf area index (LAI) and epidermal conductance (g_min ) decreased. With trait variability, SurEau predicts the plasticity of LAI and g_min buffers the impact of increasing aridity on population persistence. Knowledge of within-species variability in multiple drought tolerance traits will be crucial to accurately predict species distributional limits.

The effect of single biome occupancy on the estimation of biome shifts and the detection of biome conservatism

Biome conservatism is often regarded as common in diversifying lineages, based on the detection of low biome shift rates or high phylogenetic signal. However, many studies testing biome conservatism utilise a single-biome-per-species approach, which may influence the detection of biome conservatism. Meta-analyses show that biome shift rates are significantly lower (less than a tenth), when single biome occupancy approaches are adopted. Using New Zealand plant lineages, estimated biome shifts were also significantly lower (14-67% fewer biome shifts) when analysed under the assumption of a single biome per species. Although a single biome approach consistently resulted in lower biome shifts, it detected fewer instances of biome conservatism. A third of clades (3 out of 9) changed status in biome conservatism tests between single and multiple biome occupancy approaches, with more instances of significant biome conservatism when using a multiple biome occupancy approach. A single biome approach may change the likelihood of finding biome conservatism because it assumes biome specialisation within species, falsely recognises some biome shift types and fails to include other biome shift types. Our results indicate that the degree of biome fidelity assumed has a strong influence on analyses assessing biome shift rates, and biome conservatism testing. We advocate analyses that allow species to occupy multiple biomes.

The implications of interrelated assumptions on estimates of divergence times and rates of diversification

Phylogenies are increasingly being used as a basis to provide insight into macroevolutionary history. Here, we use simulation experiments and empirical analyses to evaluate methods that use phylogenies as a basis to make estimates of divergence times and rates of diversification. This is the first study to present a comprehensive assessment of the key variables that underpin analyses in this field - including substitution rates, speciation rates, and extinction, plus character sampling and taxon sampling. We show that in unrealistically simplistic cases (where substitution rates and speciation rates are constant, and where there is no extinction), increased character and taxon sampling lead to more accurate and precise parameter estimates. By contrast, in more complex but realistic cases (where substitution rates, speciation rates, and extinction rates vary), gains in accuracy and precision from increased character and taxon sampling are far more limited. The lack of accuracy and precision even occurs when using methods that are designed to account for more complex cases, such as relaxed clocks, fossil calibrations, and models that allow speciation rates and extinction rates to vary. The problem also persists when analysing genomic scale datasets. These results suggest two interrelated problems that occur when the processes that generated the data are more complex. First, methodological assumptions are more likely to be violated. Second, limitations in the information content of the data become more important.

Historical biogeography of Pomaderris (Rhamnaceae): Continental vicariance in Australia and repeated independent dispersals to New Zealand

AIM

Gondwanan biogeographic patterns include a combination of old vicariance events following the breakup of the supercontinent, and more recent long-distance dispersals across the southern landmasses. Floristic relationships between Australia and New Zealand have mostly been attributed to recent dispersal events rather than vicariance. We assessed the biogeographic history of Pomaderris (Rhamnaceae), which occurs in both Australia and New Zealand, by constructing a time-calibrated molecular phylogeny to infer (1) phylogenetic relationships and (2) the relative contributions of vicariance and dispersal events in the biogeographic history of the genus.

LOCATION

Australia and New Zealand.

METHODS

Using hybrid capture and high throughput sequencing, we generated nuclear and plastid data sets to estimate phylogenetic relationships and fossil calibrated divergence time estimates for Pomaderris. BioGeoBEARS and biogeographical stochastic mapping (BSM) were used to assess the ancestral area of the genus and the relative contributions of vicariance vs dispersal, and the directionality of dispersal events.

RESULTS

Our analyses indicate that Pomaderris originated in the Oligocene and had a widespread Australian distribution. Vicariance of western and eastern Australian clades coincides with the uplift of the Nullarbor Plain c. 14 Ma, followed by subsequent in-situ and within-biome diversification with little exchange across regions. A rapid radiation of southeastern Australian taxa beginning c. 10 Ma was the source for at least six independent long-distance dispersal events to New Zealand during the Pliocene-Pleistocene.

MAIN CONCLUSIONS

Our study demonstrates the importance of dispersal in explaining not only the current cross-Tasman distributions of Pomaderris, but for the New Zealand flora more broadly. The pattern of multiple independent long-distance dispersal events for Pomaderris, without significant radiation within New Zealand, is congruent with other lowland plant groups, suggesting that this biome has a different evolutionary history compared with the younger alpine flora of New Zealand, which exhibits extensive radiations often following single long distance dispersal events.

Uncertainty in Divergence Time Estimation

Understanding and representing uncertainty is crucial in academic research because it enables studies to build on the conclusions of previous studies, leading to robust advances in a particular field. Here, we evaluate the nature of uncertainty and the manner by which it is represented in divergence time estimation, a field that is fundamental to many aspects of macroevolutionary research, and where there is evidence that uncertainty has been seriously underestimated. We address this issue in the context of methods used in divergence time estimation, and with respect to the manner by which time-calibrated phylogenies are interpreted. With respect to methods, we discuss how the assumptions underlying different methods may not adequately reflect uncertainty about molecular evolution, the fossil record, or diversification rates. Therefore, divergence time estimates may not adequately reflect uncertainty and may be directly contradicted by subsequent findings. For the interpretation of time-calibrated phylogenies, we discuss how the use of time-calibrated phylogenies for reconstructing general evolutionary timescales leads to inferences about macroevolution that are highly sensitive to methodological limitations in how uncertainty is accounted for. By contrast, we discuss how the use of time-calibrated phylogenies to test specific hypotheses leads to inferences about macroevolution that are less sensitive to methodological limitations. Given that many biologists wish to use time-calibrated phylogenies to reconstruct general evolutionary timescales, we conclude that the development of methods of divergence time estimation that adequately account for uncertainty is necessary. [Divergence time estimation; macroevolution; uncertainty.].

Modeling Phylogenetic Biome Shifts on a Planet with a Past

The spatial distribution of biomes has changed considerably over deep time, so the geographical opportunity for an evolutionary lineage to shift into a new biome may depend on how the availability and connectivity of biomes has varied temporally. To better understand how lineages shift between biomes in space and time, we developed a phylogenetic biome shift model in which each lineage shifts between biomes and disperses between regions at rates that depend on the lineage's biome affinity and location relative to the spatial distribution of biomes at any given time. To study the behavior of the biome shift model in an empirical setting, we developed a literature-based representation of paleobiome structure for three mesic forest biomes, six regions, and eight time strata, ranging from the Late Cretaceous (100 Ma) through the present. We then fitted the model to a time-calibrated phylogeny of 119 Viburnum species to compare how the results responded to various realistic or unrealistic assumptions about paleobiome structure. Ancestral biome estimates that account for paleobiome dynamics reconstructed a warm temperate (or tropical) origin of Viburnum, which is consistent with previous fossil-based estimates of ancestral biomes. Imposing unrealistic paleobiome distributions led to ancestral biome estimates that eliminated support for tropical origins, and instead inflated support for cold temperate ancestry throughout the warmer Paleocene and Eocene. The biome shift model we describe is applicable to the study of evolutionary systems beyond Viburnum, and the core mechanisms of our model are extensible to the design of richer phylogenetic models of historical biogeography and/or lineage diversification. We conclude that biome shift models that account for dynamic geographical opportunities are important for inferring ancestral biomes that are compatible with our understanding of Earth history.[Ancestral states; biome shifts; historical biogeography; niche conservatism; phylogenetics].

Out of Sight, Out of Mind: Widespread Nuclear and Plastid-Nuclear Discordance in the Flowering Plant Genus Polemonium (Polemoniaceae) Suggests Widespread Historical Gene Flow Despite Limited Nuclear Signal

Phylogenomic data from a rapidly increasing number of studies provide new evidence for resolving relationships in recently radiated clades, but they also pose new challenges for inferring evolutionary histories. Most existing methods for reconstructing phylogenetic hypotheses rely solely on algorithms that only consider incomplete lineage sorting (ILS) as a cause of intra- or intergenomic discordance. Here, we utilize a variety of methods, including those to infer phylogenetic networks, to account for both ILS and introgression as a cause for nuclear and cytoplasmic-nuclear discordance using phylogenomic data from the recently radiated flowering plant genus Polemonium (Polemoniaceae), an ecologically diverse genus in Western North America with known and suspected gene flow between species. We find evidence for widespread discordance among nuclear loci that can be explained by both ILS and reticulate evolution in the evolutionary history of Polemonium. Furthermore, the histories of organellar genomes show strong discordance with the inferred species tree from the nuclear genome. Discordance between the nuclear and plastid genome is not completely explained by ILS, and only one case of discordance is explained by detected introgression events. Our results suggest that multiple processes have been involved in the evolutionary history of Polemonium and that the plastid genome does not accurately reflect species relationships. We discuss several potential causes for this cytoplasmic-nuclear discordance, which emerging evidence suggests is more widespread across the Tree of Life than previously thought. [Cyto-nuclear discordance, genomic discordance, phylogenetic networks, plastid capture, Polemoniaceae, Polemonium, reticulations.].

Spatiophylogenetic modelling of extinction risk reveals evolutionary distinctiveness and brief flowering period as threats in a hotspot plant genus

Comparative models used to predict species threat status can help identify the diagnostic features of species at risk. Such models often combine variables measured at the species level with spatial variables, causing multiple statistical challenges, including phylogenetic and spatial non-independence. We present a novel Bayesian approach for modelling threat status that simultaneously deals with both forms of non-independence and estimates their relative contribution, and we apply the approach to modelling threat status in the Australian plant genus Hakea. We find that after phylogenetic and spatial effects are accounted for, species with greater evolutionary distinctiveness and a shorter annual flowering period are more likely to be threatened. The model allows us to combine information on evolutionary history, species biology and spatial data, calculate latent extinction risk (potential for non-threatened species to become threatened), estimate the most important drivers of risk for individual species and map spatial patterns in the effects of different predictors on extinction risk. This could be of value for proactive conservation decision-making based on the early identification of species and regions of potential conservation concern.

Uncovering the genomic signature of ancient introgression between white oak lineages (Quercus)

Botanists have long recognised interspecific gene flow as a common occurrence within white oaks (Quercus section Quercus). Historical allele exchange, however, has not been fully characterised and the complex genomic signals resulting from the combination of vertical and horizontal gene transmission may confound phylogenetic inference and obscure our ability to accurately infer the deep evolutionary history of oaks. Using anchored enrichment, we obtained a phylogenomic dataset consisting of hundreds of single-copy nuclear loci. Concatenation, species-tree and network analyses were carried out in an attempt to uncover the genomic signal of ancient introgression and infer the divergent phylogenetic topology for the white oak clade. Locus and site-level likelihood comparisons were then conducted to further explore the introgressed signal within our dataset. Historical, intersectional gene flow is suggested to have occurred between an ancestor of the Eurasian Roburoid lineage and Quercus pontica and North American Dumosae and Prinoideae lineages. Despite extensive time past, our approach proved successful in detecting the genomic signature of ancient introgression. Our results, however, highlight the importance of sampling and the use of a plurality of analytical tools and methods to sufficiently explore genomic datasets, uncover this signal, and accurately infer evolutionary history.

Joint Phylogenetic Estimation of Geographic Movements and Biome Shifts during the Global Diversification of Viburnum

Phylogeny, molecular sequences, fossils, biogeography, and biome occupancy are all lines of evidence that reflect the singular evolutionary history of a clade, but they are most often studied separately, by first inferring a fossil-dated molecular phylogeny, then mapping on ancestral ranges and biomes inferred from extant species. Here we jointly model the evolution of biogeographic ranges, biome affinities, and molecular sequences, while incorporating fossils to estimate a dated phylogeny for all of the 163 extant species of the woody plant clade Viburnum (Adoxaceae) that we currently recognize in our ongoing worldwide monographic treatment of the group. Our analyses indicate that while the major Viburnum lineages evolved in the Eocene, the majority of extant species originated since the Miocene. Viburnum radiated first in Asia, in warm, broad-leaved evergreen (lucidophyllous) forests. Within Asia, we infer several early shifts into more tropical forests, and multiple shifts into forests that experience prolonged freezing. From Asia, we infer two early movements into the New World. These two lineages probably first occupied warm temperate forests and adapted later to spreading cold climates. One of these lineages (Porphyrotinus) occupied cloud forests and moved south through the mountains of the Neotropics. Several other movements into North America took place more recently, facilitated by prior adaptations to freezing in the Old World. We also infer four disjunctions between Asia and Europe: the Tinus lineage is the oldest and probably occupied warm forests when it spread, whereas the other three were more recent and in cold-adapted lineages. These results variously contradict published accounts, especially the view that Viburnum radiated initially in cold forests and, accordingly, maintained vessel elements with scalariform perforations. We explored how the location and biome assignments of fossils affected our inference of ancestral areas and biome states. Our results are sensitive to, but not entirely dependent upon, the inclusion of fossil biome data. It will be critical to take advantage of all available lines of evidence to decipher events in the distant past. The joint estimation approach developed here provides cautious hope even when fossil evidence is limited. [Biogeography; biome; combined evidence; fossil pollen; phylogeny; Viburnum.].

Increased diversification rates are coupled with higher rates of climate space exploration in Australian Acacia (Caesalpinioideae)

Australia is an excellent setting to explore relationships between climate change and diversification dynamics. Aridification since the Eocene has resulted in spectacular radiations within one or more Australian biomes. Acacia is the largest plant genus on the Australian continent, with around 1000 species, and is present in all biomes. We investigated the macroevolutionary dynamics of Acacia within climate space. We analysed phylogenetic and climatic data for 503 Acacia species to estimate a time-calibrated phylogeny and central climatic tendencies for BioClim layers from 132 000 herbarium specimens. Diversification rate heterogeneity and rates of climate space exploration were tested. We inferred two diversification rate increases, both associated with significantly higher rates of climate space exploration. Observed spikes in climate disparity within the Pleistocene correspond with onset of Pleistocene glacial-interglacial cycling. Positive time dependency in environmental disparity applies in the basal grade of Acacia, though climate space exploration rates were lower. Incongruence between rates of climate space exploration and disparity suggests different Acacia lineages have experienced different macroevolutionary processes. The second diversification rate increase is associated with a south-east Australian mesic lineage, suggesting adaptations to progressively aridifying environments and ability to transition into mesic environments contributed to Acacia's dominance across Australia.

The temporal dynamics of evolutionary diversification in Ipomoea

Molecular phylogenies are used as a basis for making inferences about macroevolutionary history. However, a robust phylogeny does not contain the information that is necessary to make many of these inferences. Complex methodologies that incorporate important assumptions about the nature of evolutionary history are therefore required. Here, we explore the implications of these assumptions for making inferences about the macroevolutionary history of Ipomoea - a large pantropical genus of flowering plants that contains the sweet potato (Ipomoea batatas), a crop of global economic importance. We focus on assumptions that underlie inferences of divergence times, and diversification parameters (speciation rates, extinction rates, and net diversification rates). These are among the most fundamental variables in macroevolutionary research. We use a series of novel approaches to explore the implications of these assumptions for inferring the age of Ipomoea, the ages of major clades within Ipomoea, whether there are significant differences in diversification parameters among clades within Ipomoea, and whether the storage root of I. batatas evolved in pre-human times. We show that inferring an age estimate for Ipomoea and major clades within Ipomoea is highly problematic. Inferred divergence times are sensitive to uncertain fossil calibrations and differing assumptions about among-branch-substitution-rate-variation. Despite this uncertainty, we are able to make robust inferences about patterns of variation in diversification parameters within Ipomoea, and that the storage root of I. batatas evolved in pre-human times. Taken together, this study presents novel and generalizable insights into the implications of methodological assumptions for making inferences about macroevolutionary history. Further, by presenting novel findings relating to the temporal dynamics of evolution in Ipomoea, as well as more specifically to I. batatas, this study makes a valuable contribution to our understanding of tropical plant evolution, and the evolutionary context in which economically important crops evolve.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

Insights from Empirical Analyses and Simulations on Using Multiple Fossil Calibrations with Relaxed Clocks to Estimate Divergence Times

Relaxed clock methods account for among-branch-rate-variation when estimating divergence times by inferring different rates for individual branches. In order to infer different rates for individual branches, important assumptions are required. This is because molecular sequence data do not provide direct information about rates but instead provide direct information about the total number of substitutions along any branch, which is a product of the rate and time for that branch. Often, the assumptions required for estimating rates for individual branches depend heavily on the implementation of multiple fossil calibrations in a single phylogeny. Here, we show that the basis of these assumptions is often critically undermined. First, we highlight that the temporal distribution of the fossil record often violates key assumptions of methods that use multiple fossil calibrations with relaxed clocks. With respect to “node calibration” methods, this conclusion is based on our inference that different fossil calibrations are unlikely to reflect the relative ages of different clades. With respect to the fossilized birth–death process, this conclusion is based on our inference that the fossil recovery rate is often highly heterogeneous. We then demonstrate that methods of divergence time estimation that use multiple fossil calibrations are highly sensitive to assumptions about the fossil record and among-branch-rate-variation. Given the problems associated with these assumptions, our results highlight that using multiple fossil calibrations with relaxed clocks often does little to improve the accuracy of divergence time estimates.

Digest: The evolutionary relationship between environment and size in Australian rodents

Biome shifts are thought to be relatively rare, but some clades shift among starkly different environments with relative ease. What causes these shifts, and how do they shape phenotypic evolution? Roycroft et al. found that biome shifts in the Pseudomys Division of murid rodents were repeatedly accompanied by body size evolution in accordance with Bergmann's rule, suggesting adaptive evolution in response to changing climate conditions.

Environment predicts repeated body size shifts in a recent radiation of Australian mammals

Closely related species that occur across steep environmental gradients often display clear body size differences, and examining this pattern is crucial to understanding how environmental variation shapes diversity. Australian endemic rodents in the Pseudomys Division (Muridae: Murinae) have repeatedly colonized the arid, monsoon, and mesic biomes over the last 5 million years. Using occurrence records, body mass data, and Bayesian phylogenetic models, we test whether body mass of 31 species in the Pseudomys Division can be predicted by their biome association. We also model the effect of eight environmental variables on body mass. Despite high phylogenetic signal in body mass evolution across the phylogeny, we find that mass predictably increases in the mesic biome and decreases in arid and monsoon biomes. As per Bergmann's rule, temperature is strongly correlated with body mass, as well as several other variables. Our results highlight two important findings. First, body size in Australian rodents has tracked with climate through the Pleistocene, likely due to several environmental variables rather than a single factor. Second, support for both Brownian motion and predictable change at different taxonomic levels in the Pseudomys Division phylogeny demonstrates how the level at which we test hypotheses can alter interpretation of evolutionary processes.

Island woodiness underpins accelerated disparification in plant radiations

The evolution of secondary (insular) woodiness and the rapid disparification of plant growth forms associated with island radiations show intriguing parallels between oceanic islands and tropical alpine sky islands. However, the evolutionary significance of these phenomena remains poorly understood and the focus of debate. We explore the evolutionary dynamics of species diversification and trait disparification across evolutionary radiations in contrasting island systems compared with their nonisland relatives. We estimate rates of species diversification, growth form evolution and phenotypic space saturation for the classical oceanic island plant radiations - the Hawaiian silverswords and Macaronesian Echium - and the well-studied sky island radiations of Lupinus and Hypericum in the Andes. We show that secondary woodiness is associated with dispersal to islands and with accelerated rates of species diversification, accelerated disparification of plant growth forms and occupancy of greater phenotypic trait space for island clades than their nonisland relatives, on both oceanic and sky islands. We conclude that secondary woodiness is a prerequisite that could act as a key innovation, manifest as the potential to occupy greater trait space, for plant radiations on island systems in general, further emphasizing the importance of combinations of clade-specific traits and ecological opportunities in driving adaptive radiations.

Phylogeny of Neotropical Sicarius sand spiders suggests frequent transitions from deserts to dry forests despite antique, broad-scale niche conservatism

Phylogenetic niche conservatism (PNC) shapes the distribution of organisms by constraining lineages to particular climatic conditions. Conversely, if areas with similar climates are geographically isolated, diversification may also be limited by dispersal. Neotropical xeric habitats provide an ideal system to test the relative roles of climate and geography on diversification, as they occur in disjunct areas with similar biotas. Sicariinae sand spiders are intimately associated with these xeric environments, particularly seasonally dry tropical forests (SDTFs) and subtropical deserts/scrublands in Africa (Hexophthalma) and the Neotropics (Sicarius). We explore the role of PNC, geography and biome shifts in their evolution and timing of diversification. We estimated a time-calibrated, total-evidence phylogeny of Sicariinae, and used published distribution records to estimate climatic niche and biome occupancy. Topologies were used for estimating ancestral niches and biome shifts. We used variation partitioning methods to test the relative importance of climate and spatially autocorrelated factors in explaining the spatial variation in phylogenetic structure of Sicarius across the Neotropics. Neotropical Sicarius are ancient and split from their African sister-group around 90 (57-131) million years ago. Most speciation events took place in the Miocene. Sicariinae records can be separated in two groups corresponding to temperate/dry and tropical/seasonally dry climates. The ancestral climatic niche of Sicariinae are temperate/dry areas, with 2-3 shifts to tropical/seasonally dry areas in Sicarius. Similarly, ancestral biomes occupied by the group are temperate and dry (deserts, Mediterranean scrub, temperate grasslands), with 2-3 shifts to tropical, seasonally dry forests and grasslands. Most of the variation in phylogenetic structure is explained by long-distance dispersal limitation that is independent of the measured climatic conditions. Sicariinae have an ancient association to arid lands, suggesting that PNC prevented them from colonizing mesic habitats. However, niches are labile at a smaller scale, with several shifts from deserts to SDTFs. This suggests that PNC and long-distance dispersal limitation played major roles in confining lineages to isolated areas of SDTF/desert over evolutionary history, although shifts between xeric biomes occurred whenever geographical opportunities were presented.

Global Succulent Biome phylogenetic conservatism across the pantropical Caesalpinia Group (Leguminosae)

The extent to which phylogenetic biome conservatism vs biome shifting determines global patterns of biodiversity remains poorly understood. To address this question, we investigated the biogeography and trajectories of biome and growth form evolution across the Caesalpinia Group (Leguminosae), a clade of 225 species of trees, shrubs and lianas distributed across the Rainforest, Succulent, Temperate and Savanna Biomes. We focused especially on the little-known Succulent Biome, an assemblage of succulent-rich, grass-poor, seasonally dry tropical vegetation distributed disjunctly across the Neotropics, Africa, Arabia and Madagascar. We reconstructed a time-calibrated phylogeny, assembled species occurrence data and assigned species to areas, biomes and growth forms. These data are used to estimate the frequency of transcontinental disjunctions, biome shifts and evolutionary transitions between growth forms and test for phylogenetic biome conservatism and correlated evolution of growth forms and biome shifts. We uncovered a pattern of strong phylogenetic Succulent Biome conservatism. We showed that transcontinental disjunctions confined within the Succulent Biome are frequent and that biome shifts to the Savanna, Rainforest and Temperate Biomes are infrequent and closely associated with shifts in plant growth forms. Our results suggest that the Succulent Biome comprises an ecologically constrained evolutionary arena spanning large geographical disjunctions across the tropics.

Tracking temporal shifts in area, biomes, and pollinators in the radiation of Salvia (sages) across continents: leveraging anchored hybrid enrichment and targeted sequence data

PREMISE OF THE STUDY

A key question in evolutionary biology is why some clades are more successful by being widespread geographically, biome diverse, or species-rich. To extend understanding of how shifts in area, biomes, and pollinators impact diversification in plants, we examined the relationships of these shifts to diversification across the mega-genus Salvia.

METHODS

A chronogram was developed from a supermatrix of anchored hybrid enrichment genomic data and targeted sequence data for over 500 of the nearly 1000 Salvia species. Ancestral areas and biomes were reconstructed using BioGeoBEARS. Pollinator guilds were scored, ancestral pollinators determined, shifts in pollinator guilds identified, and rates of pollinator switches compared.

KEY RESULTS

A well-resolved phylogenetic backbone of Salvia and updated subgeneric designations are presented. Salvia originated in Southwest Asia in the Oligocene and subsequently dispersed worldwide. Biome shifts are frequent from a likely ancestral lineage utilizing broadleaf and/or coniferous forests and/or arid shrublands. None of the four species diversification shifts are correlated to shifts in biomes. Shifts in pollination system are not correlated to species diversification shifts, except for one hummingbird shift that precedes a major shift in diversification near the crown of New World subgen. Calosphace. Multiple reversals back to bee pollination occurred within this hummingbird clade.

CONCLUSIONS

Salvia diversified extensively in different continents, biomes, and with both bee and bird pollinators. The lack of tight correlation of area, biome, and most pollinator shifts to the four documented species diversification shifts points to other important drivers of speciation in Salvia.

Monocot plastid phylogenomics, timeline, net rates of species diversification, the power of multi-gene analyses, and a functional model for the origin of monocots

PREMISE OF THE STUDY

We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification.

METHODS

We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth.

KEY RESULTS

Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches.

CONCLUSIONS

We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

Miocene biome turnover drove conservative body size evolution across Australian vertebrates

On deep time scales, changing climatic trends can have a predictable influence on macroevolution. From evidence of mass extinctions, we know that rapid climatic oscillations can indirectly open niche space and precipitate adaptive radiation, changing the course of ecological diversification. These dramatic shifts in the global climate, however, are rare events relative to extended periods of protracted climate change and biome turnover. It remains unclear whether during gradually changing periods, shifting habitats may instead promote non-adaptive speciation by facilitating allopatry and phenotypic conservatism. Using fossil-calibrated, species-level phylogenies for five Australian radiations comprising more than 800 species, we investigated temporal trends in biogeography and body size evolution. Here, we demonstrate that gradual Miocene cooling and aridification correlates with the restricted phenotypic diversification of multiple ecologically diverse vertebrate groups. This probably occurred as species ranges became fractured and isolated during continental biome restructuring, encouraging a shift towards conservatism in body size evolution. Our results provide further evidence that abiotic changes, not only biotic interactions, may act as selective forces influencing phenotypic macroevolution.