The Radiation of Darwin’s Giant Daisies in the Galápagos Islands

Genomes of Meniocus linifolius and Tetracme quadricornis reveal the ancestral karyotype and genomic features of core Brassicaceae

Brassicaceae represents an important plant family from both a scientific and economic perspective. However, genomic features related to the early diversification of this family have not been fully characterized, especially upon the uplift of the Tibetan Plateau, which was followed by increasing aridity in the Asian interior, intensifying monsoons in Eastern Asia, and significantly fluctuating daily temperatures. Here, we reveal the genomic architecture that accompanied early Brassicaceae diversification by analyzing two high-quality chromosome-level genomes for Meniocus linifolius (Arabodae; clade D) and Tetracme quadricornis (Hesperodae; clade E), together with genomes representing all major Brassicaceae clades and the basal Aethionemeae. We reconstructed an ancestral core Brassicaceae karyotype (CBK) containing 9 pseudochromosomes with 65 conserved syntenic genomic blocks and identified 9702 conserved genes in Brassicaceae. We detected pervasive conflicting phylogenomic signals accompanied by widespread ancient hybridization events, which correlate well with the early divergence of core Brassicaceae. We identified a successive Brassicaceae-specific expansion of the class I TREHALOSE-6-PHOSPHATE SYNTHASE 1 (TPS1) gene family, which encodes enzymes with essential regulatory roles in flowering time and embryo development. The TPS1s were mainly randomly amplified, followed by expression divergence. Our results provide fresh insights into historical genomic features coupled with Brassicaceae evolution and offer a potential model for broad-scale studies of adaptive radiation under an ever-changing environment.

Phylogenetic and phylogeographic evidence for the colonization success of the diplochorous Astydamia latifolia across the Canary Islands

Astydamia latifolia is the only species of the genus Astydamia, which forms an early-diverging lineage of Apiaceae, subfamily Apioideae. This species is subendemic to the Canary Islands and one of the most representative species of the coastal environments of this archipelago. Astydamia displays diplochory, that is, diaspores with two long-distance dispersal (LDD) syndromes. In particular, A. latifolia has both anemochorous and thalassochorous fruit traits (corky and winged mericarps). Although we expect this species to have a high dispersal capacity, there is no genetic study addressing it. The divergence time of this species from its sister taxon is also unknown. In this study, we aimed (i) to unveil the phylogenetic relationships and divergence times of A. latifolia; (ii) to reconstruct its phylogeographic structure across the Canary Islands; and (iii) to estimate the number of inter-island colonization events. To these ends, we first sequenced the internal transcribed spacer (ITS) region for A. latifolia, reconstructed the phylogenetic relationships of Astydamia and Apiaceae relatives and estimated divergence times. Then, two plastid DNA regions (psaI-aacD and psbK-trnS) were sequenced for 158 individuals (from 36 Canarian population and one NW African population) to reconstruct a haplotype network. The recently developed method Phylogeographic Analysis of Island Colonization Events (PAICE) was used to estimate the number of inter-island colonization events. Results show that A. latifolia is a phylogenetically isolated lineage that diverged from the most closely related genus (Molopospermum) in the Eocene-Miocene. It displays a low plastid DNA diversity (only four haplotypes detected), which is accompanied by a high degree of haplotype sharing between islands and highly linear rarefaction curves of colonization obtained in PAICE. These findings suggest a high colonization ability for this species, most likely related to the availability of two LDD syndromes.

Parallel evolution, atavism, and extensive introgression explain the radiation of Epimedium sect. Diphyllon (Berberidaceae) in southern East Asia

East Asia is the richest region of plant biodiversity in the northern temperate zone, and its radiation provides key insights for understanding rapid speciation, including evolutionary patterns and processes. However, it is challenging to investigate the recent evolutionary radiation among plants because of the lack of genetic divergence, phenotypic convergence, and interspecific gene flow. Epimedium sect. Diphyllon is a rarely studied plant lineage endemic to East Asia, especially highly diversified in its southern part. In this study, we report a robust phylogenomic analysis based on genotyping-by-sequencing data of this lineage. The results revealed a clear biogeographic pattern for Epimedium sect. Diphyllon with recognition into two major clades corresponding to the Sino-Himalayan and Sino-Japanese subkingdoms of East Asian Flora and rapid diversification of the extant species dated to the Pleistocene. Evolutionary radiation of Epimedium sect. Diphyllon is characterized by recent and predominant parallel evolution and atavism between the two subkingdom regions, with extensive reticulating hybridization within each region during the course of diversification in southern East Asia. A parallel-atavism-introgression hypothesis is referred to in explaining the radiation of plant diversity in southern East Asia, which represents a potential model for the rapid diversification of plants under global climate cooling in the late Tertiary. Our study advances our understanding of the evolutionary processes of plant radiation in East Asia as well as in other biodiversity hotspot regions.

Climatic stability and geological history shape global centers of neo- and paleoendemism in seed plants

Assessing the distribution of geographically restricted and evolutionarily unique species and their underlying drivers is key to understanding biogeographical processes and critical for global conservation prioritization. Here, we quantified the geographic distribution and drivers of phylogenetic endemism for ~320,000 seed plants worldwide and identified centers and drivers of evolutionarily young (neoendemism) and evolutionarily old endemism (paleoendemism). Tropical and subtropical islands as well as tropical mountain regions displayed the world's highest phylogenetic endemism. Most tropical rainforest regions emerged as centers of paleoendemism, while most Mediterranean-climate regions showed high neoendemism. Centers where high neo- and paleoendemism coincide emerged on some oceanic and continental fragment islands, in Mediterranean-climate regions and parts of the Irano-Turanian floristic region. Global variation in phylogenetic endemism was well explained by a combination of past and present environmental factors (79.8 to 87.7% of variance explained) and most strongly related to environmental heterogeneity. Also, warm and wet climates, geographic isolation, and long-term climatic stability emerged as key drivers of phylogenetic endemism. Neo- and paleoendemism were jointly explained by climatic and geological history. Long-term climatic stability promoted the persistence of paleoendemics, while the isolation of oceanic islands and their unique geological histories promoted neoendemism. Mountainous regions promoted both neo- and paleoendemism, reflecting both diversification and persistence over time. Our study provides insights into the evolutionary underpinnings of biogeographical patterns in seed plants and identifies the areas on Earth with the highest evolutionary and biogeographical uniqueness-key information for setting global conservation priorities.

Comparative analysis of Diospyros (Ebenaceae) plastomes: Insights into genomic features, mutational hotspots, and adaptive evolution

Diospyros (Ebenaceae) is a widely distributed genus of trees and shrubs from pantropical to temperate regions, with numerous species valued for their fruits (persimmons), timber, and medicinal values. However, information regarding their plastomes and chloroplast evolution is scarce. The present study performed comparative genomic and evolutionary analyses on plastomes of 45 accepted Diospyros species, including three newly sequenced ones. Our study showed a highly conserved genomic structure across the Diospyros species, with 135-136 encoding genes, including 89 protein-coding genes, 1-2 pseudogenes (Ψycf1 for all, Ψrps19 for a few), 37 tRNA genes and 8 rRNA genes. Comparative analysis of Diospyros identified three intergenic regions (ccsA-ndhD, rps16-psbK and petA-psbJ) and five genes (rpl33, rpl22, petL, psaC and rps15) as the mutational hotspots in these species. Phylogenomic analysis identified the phylogenetic position of three newly sequenced ones and well supported a monophylogenetic (sub)temperate taxa and four clades in the pantropical taxa. The analysis codon usage identified 30 codons with relative synonymous codon usage (RSCU) values >1 and 29 codons ending with A and U bases. A total of three codons (UUA, GCU, and AGA) with highest RSCU values were identified as the optimal codons. Effective number of codons (ENC)-plot indicated the significant role of mutational pressure in shaping codon usage, while most protein-coding genes in Diospyros experienced relaxed purifying selection (d _N/d _S < 1). Additionally, the psbH gene showed positive selection (d _N/d _S > 1) in the (sub)temperate species. Thus, the results provide a meaningful foundation for further elaborating Diospyros's genetic architecture and taxonomy, enriching genetic diversity and conserving genetic resources.

Evolutionary genomics of oceanic island radiations

A recurring feature of oceanic archipelagos is the presence of adaptive radiations that generate endemic, species-rich clades that can offer outstanding insight into the links between ecology and evolution. Recent developments in evolutionary genomics have contributed towards solving long-standing questions at this interface. Using a comprehensive literature search, we identify studies spanning 19 oceanic archipelagos and 110 putative adaptive radiations, but find that most of these radiations have not yet been investigated from an evolutionary genomics perspective. Our review reveals different gaps in knowledge related to the lack of implementation of genomic approaches, as well as undersampled taxonomic and geographic areas. Filling those gaps with the required data will help to deepen our understanding of adaptation, speciation, and other evolutionary processes.

Reference-Based RADseq Unravels the Evolutionary History of Polar Species in 'the Crux Lichenologorum' Genus Usnea (Parmeliaceae, Ascomycota)

Nearly 90% of fungal diversity, one of the most speciose branches in the tree of life, remains undescribed. Lichenized fungi as symbiotic associations are still a challenge for species delimitation, and current species diversity is vastly underestimated. The ongoing democratization of Next-Generation Sequencing is turning the tables. Particularly, reference-based RADseq allows for metagenomic filtering of the symbiont sequence and yields robust phylogenomic trees of closely related species. We implemented reference-based RADseq to disentangle the evolution of neuropogonoid lichens, which inhabit harsh environments and belong to Usnea (Parmeliaceae, Ascomycota), one of the most taxonomically intriguing genera within lichenized fungi. Full taxon coverage of neuropogonoid lichens was sampled for the first time, coupled with phenotype characterizations. More than 20,000 loci of 126 specimens were analyzed through concatenated and coalescent-based methods, including time calibrations. Our analysis addressed the major taxonomic discussions over recent decades. Subsequently, two species are newly described, namely U. aymondiana and U. fibriloides, and three species names are resurrected. The late Miocene and Pliocene-Pleistocene boundary is inferred as the timeframe for neuropogonoid lichen diversification. Ultimately, this study helped fill the gap of fungal diversity by setting a solid backbone phylogeny which raises new questions about which factors may trigger complex evolutionary scenarios.

Assessing the relative performance of fast molecular dating methods for phylogenomic data

Advances in genome sequencing techniques produced a significant growth of phylogenomic datasets. This massive amount of data represents a computational challenge for molecular dating with Bayesian approaches. Rapid molecular dating methods have been proposed over the last few decades to overcome these issues. However, a comparative evaluation of their relative performance on empirical data sets is lacking. We analyzed 23 empirical phylogenomic datasets to investigate the performance of two commonly employed fast dating methodologies: penalized likelihood (PL), implemented in treePL, and the relative rate framework (RRF), implemented in RelTime. They were compared to Bayesian analyses using the closest possible substitution models and calibration settings. We found that RRF was computationally faster and generally provided node age estimates statistically equivalent to Bayesian divergence times. PL time estimates consistently exhibited low levels of uncertainty. Overall, to approximate Bayesian approaches, RelTime is an efficient method with significantly lower computational demand, being more than 100 times faster than treePL. Thus, to alleviate the computational burden of Bayesian divergence time inference in the era of massive genomic data, molecular dating can be facilitated using the RRF, allowing evolutionary hypotheses to be tested more quickly and efficiently.

Climatic niche pre-adaptation facilitated island colonization followed by budding speciation in the Madeiran ivy (Hedera maderensis, Araliaceae)

The path followed by species in the colonization of remote oceanic islands ultimately depends on their phylogenetic constraints and ecological responses. In this study, we aim to evaluate the relative role of geographical and ecological forces in the origin and evolution of the Madeiran ivy (Hedera maderensis), a single-species endemic belonging to the western polyploid clade of Hedera. To determine the phylogenetic placement of H. maderensis within the western polyploid clade, we analyzed 40 populations (92 individuals) using genotyping-by-sequencing and including Hedera helix as outgroup. Climatic niche differences among the study species were evaluated using a database with 867 records representing the entire species ranges. To test species responses to climate, 13 vegetative and reproductive functional traits were examined for 70 populations (335 individuals). Phylogenomic results revealed a nested pattern with H. maderensis embedded within the south-western Iberian H. iberica. Gradual niche differentiation from the coldest and most continental populations of H. iberica to the warm and stable coastal population sister to H. maderensis parallels the geographical pattern observed in the phylogeny. Similarity in functional traits is observed for H. maderensis and H. iberica. The two species show leaves with higher specific leaf area (SLA), lower leaf dry matter content (LDMC) and thickness and fruits with lower pulp fraction than the other western polyploid species H. hibernica. Acquisition of a Macaronesian climatic niche and the associated functional syndrome in mainland European ivies (leaves with high SLA, and low LDMC and thickness, and fruits with less pulp content) was a key step in the colonization of Madeira by the H. iberica/H. maderensis lineage, which points to climatic pre-adaptation as key in the success of island colonization (dispersal and establishment). Once in Madeira, budding speciation was driven by geographical isolation, while ecological processes are regarded as secondary forces with a putative impact in the lack of further in situ diversification.

The genomic basis of the plant island syndrome in Darwin's giant daisies

The repeated, rapid and often pronounced patterns of evolutionary divergence observed in insular plants, or the ‘plant island syndrome’, include changes in leaf phenotypes, growth, as well as the acquisition of a perennial lifestyle. Here, we sequence and describe the genome of the critically endangered, Galápagos-endemic species Scalesia atractyloides Arnot., obtaining a chromosome-resolved, 3.2-Gbp assembly containing 43,093 candidate gene models. Using a combination of fossil transposable elements, k-mer spectra analyses and orthologue assignment, we identify the two ancestral genomes, and date their divergence and the polyploidization event, concluding that the ancestor of all extant Scalesia species was an allotetraploid. There are a comparable number of genes and transposable elements across the two subgenomes, and while their synteny has been mostly conserved, we find multiple inversions that may have facilitated adaptation. We identify clear signatures of selection across genes associated with vascular development, growth, adaptation to salinity and flowering time, thus finding compelling evidence for a genomic basis of the island syndrome in one of Darwin’s giant daisies.

Many island plant species share a syndrome of characteristic phenotype and life history. Cerca et al. find the genomic basis of the plant island syndrome in one of Darwin’s giant daisies, while separating ancestral genomes in a chromosome-resolved polyploid assembly.

Multiplexed shotgun genotyping (MSG) data resolve phylogenetic relationships within and among archipelagos in Macaronesian Tolpis

PREMISE

Plants endemic to oceanic archipelagos are suitable for studying evolution, being isolated on substrates of different ages. Evolution has been recent, rendering traditionally employed sequences insufficiently variable for resolving relationships. This study includes sampling in the genus Tolpis (Asteraceae) from the Azores, Madeira, and Cape Verde, and expands upon an earlier study demonstrating the efficacy of multiplexed shotgun genotyping (MSG) for resolving relationships in Canarian Tolpis.

METHODS

Genomic libraries for 90 accessions of Tolpis and two from the outgroup were generated for genotyping individuals using MSG. Loci were de novo assembled with iPyrad, which clusters MSG loci within and between samples. A maximum likelihood phylogeny was generated with RAxML. Ancestral area reconstruction was inferred using R package BioGeoBEARS.

RESULTS

MSG data recovered a highly resolved phylogeny from population to inter-archipelago levels. Ancestral area reconstruction provided biogeographic hypotheses for the radiation of Macaronesian Tolpis.

CONCLUSIONS

Four major clades were resolved. The Madeiran endemic T. macrorhiza is sister to other Tolpis. Species from the Canaries, Cape Verdes, and the continent are sister to T. succulenta from Madeira, which has a sister subclade of Azorean populations composed of T. succulenta and T. azorica. Population-level resolution suggests unrecognized taxa on several archipelagos. Ancestral reconstruction suggests initial dispersal from the continent to Madeira, with dispersal to the Azores, then dispersal from Madeira to the Canary Islands, with both subsequent dispersal to the Cape Verdes and back-dispersal to the continent. Single-island radiations and inter-island dispersal are implicated in divergence in Macaronesian Tolpis.

exTREEmaTIME: a method for incorporating uncertainty into divergence time estimates

We present a method of divergence time estimation (exTREEmaTIME) that aims to effectively account for uncertainty in divergence time estimates. The method requires a minimal set of assumptions, and, based on these assumptions, estimates the oldest possible divergence times and youngest possible divergence times that are consistent with the assumptions. We use a series of simulations and empirical analyses to illustrate that exTREEmaTIME is effective at representing uncertainty. We then describe how exTREEmaTIME can act as a basis to determine the implications of the more stringent assumptions that are incorporated into other methods of divergence time estimation that produce more precise estimates. This is critically important given that many of the assumptions that are incorporated into these methods are highly complex, difficult to justify biologically, and as such can lead to estimates that are highly inaccurate. This article has an associated First Person interview with the first author of the paper.

Spatiotemporal evolution of the global species diversity of Rhododendron

Evolutionary radiation is a widely recognized mode of species diversification, but its underlying mechanisms have not been unambiguously resolved for species-rich cosmopolitan plant genera. In particular, it remains largely unknown how biological and environmental factors have jointly driven its occurrence in specific regions. Here, we use Rhododendron, the largest genus of woody plants in the Northern Hemisphere, to investigate how geographic and climatic factors, as well as functional traits, worked together to trigger plant evolutionary radiations and shape the global patterns of species richness based on a solid species phylogeny. Using 3,437 orthologous nuclear genes, we reconstructed the first highly supported and dated backbone phylogeny of Rhododendron comprising 200 species that represent all subgenera, sections, and nearly all multispecies subsections, and found that most extant species originated by evolutionary radiations when the genus migrated southward from circumboreal areas to tropical/subtropical mountains, showing rapid increases of both net diversification rate and evolutionary rate of environmental factors in the Miocene. We also found that the geographically uneven diversification of Rhododendron led to a much higher diversity in Asia than in other continents, which was mainly driven by two environmental variables, that is, elevation range and annual precipitation, and were further strengthened by the adaptation of leaf functional traits. Our study provides a good example of integrating phylogenomic and ecological analyses in deciphering the mechanisms of plant evolutionary radiations, and sheds new light on how the intensification of the Asian monsoon has driven evolutionary radiations in large plant genera of the Himalaya-Hengduan Mountains.

Genome-Wide Single Nucleotide Polymorphism Analysis Elucidates the Evolution of Prunus takesimensis in Ulleung Island: The Genetic Consequences of Anagenetic Speciation

Of the two major speciation modes of endemic plants on oceanic islands, cladogenesis and anagenesis, the latter has been recently emphasized as an effective mechanism for increasing plant diversity in isolated, ecologically homogeneous insular settings. As the only flowering cherry occurring on Ulleung Island in the East Sea (concurrently known as Sea of Japan), Prunus takesimensis Nakai has been presumed to be derived through anagenetic speciation on the island. Based on morphological similarities, P. sargentii Rehder distributed in adjacent continental areas and islands has been suggested as a purported continental progenitor. However, the overall genetic complexity and resultant non-monophyly of closely related flowering cherries have hindered the determination of their phylogenetic relationships as well as the establishment of concrete continental progenitors and insular derivative relationships. Based on extensive sampling of wild flowering cherries, including P. takesimensis and P. sargentii from Ulleung Island and its adjacent areas, the current study revealed the origin and evolution of P. takesimensis using multiple molecular markers. The results of phylogenetic reconstruction and population genetic structure analyses based on single nucleotide polymorphisms detected by multiplexed inter-simple sequence repeat genotyping by sequencing (MIG-seq) and complementary cpDNA haplotypes provided evidence for (1) the monophyly of P. takesimensis; (2) clear genetic differentiation between P. takesimensis (insular derivative) and P. sargentii (continental progenitor); (3) uncertain geographic origin of P. takesimensis, but highly likely via single colonization from the source population of P. sargentii in the Korean Peninsula; (4) no significant reduction in genetic diversity in anagenetically derived insular species, i.e., P. takesimensis, compared to its continental progenitor P. sargentii; (5) no strong population genetic structuring or geographical patterns in the insular derivative species; and (6) MIG-seq method as an effective tool to elucidate the complex evolutionary history of plant groups.

The Galápagos Islands: biogeographic patterns and geology

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

How to become a crab: Phenotypic constraints on a recurring body plan

A fundamental question in biology is whether phenotypes can be predicted by ecological or genomic rules. At least five cases of convergent evolution of the crab-like body plan (with a wide and flattened shape, and a bent abdomen) are known in decapod crustaceans, and have, for over 140 years, been known as "carcinization." The repeated loss of this body plan has been identified as "decarcinization." In reviewing the field, we offer phylogenetic strategies to include poorly known groups, and direct evidence from fossils, that will resolve the history of crab evolution and the degree of phenotypic variation within crabs. Proposed ecological advantages of the crab body are summarized into a hypothesis of phenotypic integration suggesting correlated evolution of the carapace shape and abdomen. Our premise provides fertile ground for future studies of the genomic and developmental basis, and the predictability, of the crab-like body form.

Polyploidy on Islands: Its Emergence and Importance for Diversification

Whole genome duplication or polyploidy is widespread among floras globally, but traditionally has been thought to have played a minor role in the evolution of island biodiversity, based on the low proportion of polyploid taxa present. We investigate five island systems (Juan Fernández, Galápagos, Canary Islands, Hawaiian Islands, and New Zealand) to test whether polyploidy (i) enhances or hinders diversification on islands and (ii) is an intrinsic feature of a lineage or an attribute that emerges in island environments. These island systems are diverse in their origins, geographic and latitudinal distributions, levels of plant species endemism (37% in the Galapagos to 88% in the Hawaiian Islands), and ploidy levels, and taken together are representative of islands more generally. We compiled data for vascular plants and summarized information for each genus on each island system, including the total number of species (native and endemic), generic endemicity, chromosome numbers, genome size, and ploidy levels. Dated phylogenies were used to infer lineage age, number of colonization events, and change in ploidy level relative to the non-island sister lineage. Using phylogenetic path analysis, we then tested how the diversification of endemic lineages varied with the direct and indirect effects of polyploidy (presence of polyploidy, time on island, polyploidization near colonization, colonizer pool size) and other lineage traits not associated with polyploidy (time on island, colonizer pool size, repeat colonization). Diploid and tetraploid were the most common ploidy levels across all islands, with the highest ploidy levels (>8x) recorded for the Canary Islands (12x) and New Zealand (20x). Overall, we found that endemic diversification of our focal island floras was shaped by polyploidy in many cases and certainly others still to be detected considering the lack of data in many lineages. Polyploid speciation on the islands was enhanced by a larger source of potential congeneric colonists and a change in ploidy level compared to overseas sister taxa.