Land plant evolutionary timeline: gene effects are secondary to fossil constraints in relaxed clock estimation of age and substitution rates

Hunting the Snark: the flawed search for mythical Jurassic angiosperms

Several recent palaeobotanical studies claim to have found and described pre-Cretaceous angiosperm macrofossils. With rare exceptions, these papers fail to define a flower, do not acknowledge that fossils require character-based rather than group-based classification, do not explicitly state which morphological features would unambiguously identify a fossil as angiospermous, ignore the modern conceptual framework of phylogeny reconstruction, and infer features in the fossils in question that are interpreted differently by (or even invisible to) other researchers. This unfortunate situation is compounded by the relevant fossils being highly disarticulated two-dimensional compression-impressions lacking anatomical preservation. Given current evidence, all supposed pre-Cretaceous angiosperms are assignable to other major clades among the gymnosperms sensu lato. By any workable morphological definition, flowers are not confined to, and therefore cannot delimit, the angiosperm clade. More precisely defined character states that are potentially diagnostic of angiosperms must by definition originate on the phylogenetic branch that immediately precedes the angiosperm crown group. Although the most reliable candidates for diagnostic characters (triploid endosperm reflecting double fertilization, closed carpel, bitegmic ovule, and phloem companion cells) are rarely preserved and/or difficult to detect unambiguously, similar characters have occasionally been preserved in high-quality permineralized non-angiosperm fossils. The angiosperm radiation documented by Early Cretaceous fossils involves only lineages closely similar to extant taxonomic families, lacks obvious morphological gaps, and (as agreed by both the fossil record and molecular phylogenies) was relatively rapid-all features that suggest a primary radiation. It is unlikely that ancestors of the crown group common ancestor would have fulfilled a character-based definition of (and thereby required expansion of the concept of) an angiosperm; they would instead form a new element of the non-angiosperm members of the 'anthophyte' grade, competing with Caytonia to be viewed as morphologically determined sister group for angiosperms. Conclusions drawn from molecular phylogenetics should not be allowed to routinely constrain palaeobotanical inferences; reciprocal illumination between different categories of data offers greater explanatory power than immediately resorting to Grand Syntheses. The Jurassic angiosperm-essentially a product of molecular phylogenetics-may have become the holy grail of palaeobotany but it appears equally mythical.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution

PREMISE

Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution.

METHODS

We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements.

RESULTS

Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses.

CONCLUSIONS

A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

Large-scale phylogenomic analysis suggests three ancient superclades of the WUSCHEL-RELATED HOMEOBOX transcription factor family in plants

The adaptation of plants to land required multiple morphological innovations. Among these include a variety of lateral organs that are initiated from apical meristems, in which the mantainance of undifferentiated stem cells is regulated by the homeodomain WUSCHEL-RELATED (WOX) transcription factors. Expansion of the WOX gene family has been associated with whole genome duplication (WGD) events and postulated to have been pivotal to the evolution of morphological complexity in land plants. Previous studies have classified the WOX gene family into three superclades (e.g., the ancient clade, the intermediate clade, and the modern clade). In order to improve our understanding of the evolution of the WOX gene family, we surveyed the WOX gene sequences from 38 genomes and 440 transcriptomes spanning the Viridiplantae and Rhodophyta. The WOX phylogeny inferred from 1039 WOX proteins drawn from 267 species with improved support along the backbone of the phylogeny suggests that the plant-specific WOX family contains three ancient superclades, which we term Type 1 (T1WOX, the WOX10/13/14 clade), Type 2 (T2WOX, the WOX8/9 and WOX11/12 clades), and Type 3 (T3WOX, the WUS, WOX1/6, WOX2, WOX3, WOX4 and WOX5/7 clades). Divergence of the T1WOX and T2WOX superclades may predate the diversification of vascular plants. Synteny analysis suggests contribution of WGD to expansion of the WOX family. Promoter analysis finds that the capacity of the WOX genes to be regulated by the auxin and cytokinin signaling pathways may be deeply conserved in the Viridiplantae. This study improves our phylogenetic context for elucidating functional evolution of the WOX gene family, which has likely contributed to the morphological complexity of land plants.

Recalibration of the insect evolutionary time scale using Monte San Giorgio fossils suggests survival of key lineages through the End-Permian Extinction

Insects are a highly diverse group of organisms and constitute more than half of all known animal species. They have evolved an extraordinary range of traits, from flight and complete metamorphosis to complex polyphenisms and advanced eusociality. Although the rich insect fossil record has helped to chart the appearance of many phenotypic innovations, data are scarce for a number of key periods. One such period is that following the End-Permian Extinction, recognized as the most catastrophic of all extinction events. We recently discovered several 240-million-year-old insect fossils in the Mount San Giorgio Lagerstätte (Switzerland-Italy) that are remarkable for their state of preservation (including internal organs and soft tissues), and because they extend the records of their respective taxa by up to 200 million years. By using these fossils as calibrations in a phylogenomic dating analysis, we present a revised time scale for insect evolution. Our date estimates for several major lineages, including the hyperdiverse crown groups of Lepidoptera, Hemiptera: Heteroptera and Diptera, are substantially older than their currently accepted post-Permian origins. We found that major evolutionary innovations, including flight and metamorphosis, appeared considerably earlier than previously thought. These results have numerous implications for understanding the evolution of insects and their resilience in the face of extreme events such as the End-Permian Extinction.

Fossil evidence reveals how plants responded to cooling during the Cretaceous-Paleogene transition

BACKGROUND

Around the Cretaceous-Paleogene (K-Pg) boundary, an obvious global cooling occurred, which resulted in dramatic changes in terrestrial ecosystems and the evolutionary trends of numerous organisms. However, how plant lineages responded to the cooling has remained unknown until now. Between ca. 70-60 Ma Mesocyparis McIver & Basinger (Cupressaceae), an extinct conifer genus, was distributed from eastern Asia to western North America and provides an excellent opportunity to solve this riddle.

RESULTS

Here we report a new species, Mesocyparis sinica from the early Paleocene of Jiayin, Heilongjiang, northeastern China. By integrating lines of evidence from phylogeny and comparative morphology of Mesocyparis, we found that during ca.70-60 Ma, the size of seed cone of Mesocyparis more than doubled, probably driven by the cooling during the K-Pg transition, which might be an effective adaptation for seed dispersal by animals. More importantly, we discovered that the northern limit of this genus, as well as those of two other arboreal taxa Metasequoia Miki ex Hu et Cheng (gymnosperm) and Nordenskioldia Heer (angiosperm), migrated ca.4-5° southward in paleolatitude during this time interval.

CONCLUSIONS

Our results suggest that the cooling during the K-Pg transition may have been responsible for the increase in size of the seed cone of Mesocyparis and have driven the migration of plants southwards.

Battistuzzi F, Kumar S. A Machine Learning Method for Detecting Autocorrelation of Evolutionary Rates in Large Phylogenies

New species arise from pre-existing species and inherit similar genomes and environments. This predicts greater similarity of the tempo of molecular evolution between direct ancestors and descendants, resulting in autocorrelation of evolutionary rates in the tree of life. Surprisingly, molecular sequence data have not confirmed this expectation, possibly because available methods lack the power to detect autocorrelated rates. Here, we present a machine learning method, CorrTest, to detect the presence of rate autocorrelation in large phylogenies. CorrTest is computationally efficient and performs better than the available state-of-the-art method. Application of CorrTest reveals extensive rate autocorrelation in DNA and amino acid sequence evolution of mammals, birds, insects, metazoans, plants, fungi, parasitic protozoans, and prokaryotes. Therefore, rate autocorrelation is a common phenomenon throughout the tree of life. These findings suggest concordance between molecular and nonmolecular evolutionary patterns, and they will foster unbiased and precise dating of the tree of life.

Assembly of modern mammal community structure driven by Late Cretaceous dental evolution, rise of flowering plants, and dinosaur demise

The long-standing view that Mesozoic mammaliaforms living in dinosaur-dominated ecosystems were ecologically constrained to small size and insectivory has been challenged by astonishing fossil discoveries over the last three decades. By studying these well-preserved early mammaliaform specimens, paleontologists now agree that mammaliaforms underwent ecomorphological diversification during the Mesozoic Era. This implies that Mesozoic mammaliaform communities had ecological structure and breadth that were comparable to today's small-bodied mammalian communities. However, this hypothesis remains untested in part because the primary focus of most studies is on individual taxa. Here, we present a study quantifying the ecological structure of Mesozoic mammaliaform communities with the aim of identifying evolutionary and ecological drivers that influenced the deep-time assembly of small-bodied mammaliaform communities. We used body size, dietary preference, and locomotor mode to establish the ecospace occupation of 98 extant, small-bodied mammalian communities from diverse biomes around the world. We calculated ecological disparity and ecological richness to measure the magnitude of ecological differences among species in a community and the number of different eco-cells occupied by species of a community, respectively. This modern dataset served as a reference for analyzing five exceptionally preserved, extinct mammaliaform communities (two Jurassic, two Cretaceous, one Eocene) from Konservat-Lagerstätten. Our results indicate that the interplay of at least three factors, namely the evolution of the tribosphenic molar, the ecological rise of angiosperms, and potential competition with other vertebrates, may have been critical in shaping the ecological structure of small-bodied mammaliaform communities through time.

Origin of angiosperms and the puzzle of the Jurassic gap

Angiosperms are by far the most species-rich clade of land plants, but their origin and early evolutionary history remain poorly understood. We reconstructed angiosperm phylogeny based on 80 genes from 2,881 plastid genomes representing 85% of extant families and all orders. With a well-resolved plastid tree and 62 fossil calibrations, we dated the origin of the crown angiosperms to the Upper Triassic, with major angiosperm radiations occurring in the Jurassic and Lower Cretaceous. This estimated crown age is substantially earlier than that of unequivocal angiosperm fossils, and the difference is here termed the 'Jurassic angiosperm gap'. Our time-calibrated plastid phylogenomic tree provides a highly relevant framework for future comparative studies of flowering plant evolution.

Phylogenomics resolves the deep phylogeny of seed plants and indicates partial convergent or homoplastic evolution between Gnetales and angiosperms

After decades of molecular phylogenetic studies, the deep phylogeny of gymnosperms has not been resolved, and the phylogenetic placement of Gnetales remains one of the most controversial issues in seed plant evolution. To resolve the deep phylogeny of seed plants and to address the sources of phylogenetic conflict, we conducted a phylotranscriptomic study with a sampling of all 13 families of gymnosperms and main lineages of angiosperms. Multiple datasets containing up to 1 296 042 sites across 1308 loci were analysed, using concatenation and coalescence approaches. Our study generated a consistent and well-resolved phylogeny of seed plants, which places Gnetales as sister to Pinaceae and thus supports the Gnepine hypothesis. Cycads plus Ginkgo is sister to the remaining gymnosperms. We also found that Gnetales and angiosperms have similar molecular evolutionary rates, which are much higher than those of other gymnosperms. This implies that Gnetales and angiosperms might have experienced similar selective pressures in evolutionary histories. Convergent molecular evolution or homoplasy is partially responsible for the phylogenetic conflicts in seed plants. Our study provides a robustly reconstructed backbone phylogeny that is important for future molecular and morphological studies of seed plants, in particular gymnosperms, in the light of evolution.

Nuclear protein phylogenies support the monophyly of the three bryophyte groups (Bryophyta Schimp.)

Unraveling the phylogenetic relationships between the four major lineages of terrestrial plants (mosses, liverworts, hornworts, and vascular plants) is essential for an understanding of the evolution of traits specific to land plants, such as their complex life cycles, and the evolutionary development of stomata and vascular tissue. Well supported phylogenetic hypotheses resulting from different data and methods are often incongruent due to processes of nucleotide evolution that are difficult to model, for example substitutional saturation and composition heterogeneity. We reanalysed a large published dataset of nuclear data and modelled these processes using degenerate-codon recoding and tree-heterogeneous composition substitution models. Our analyses resolved bryophytes as a monophyletic group and showed that the nonnonmonophyly of the clade that is supported by the analysis of nuclear nucleotide data is due solely to fast-evolving synonymous substitutions. The current congruence among phylogenies of both nuclear and chloroplast analyses lent considerable support to the conclusion that the bryophytes are a monophyletic group. An initial split between bryophytes and vascular plants implies that the bryophyte life cycle (with a dominant gametophyte nurturing an unbranched sporophyte) may not be ancestral to all land plants and that stomata are likely to be a symplesiomorphy among embryophytes.

Lycophyte plastid genomics: extreme variation in GC, gene and intron content and multiple inversions between a direct and inverted orientation of the rRNA repeat

Lycophytes are a key group for understanding vascular plant evolution. Lycophyte plastomes are highly distinct, indicating a dynamic evolutionary history, but detailed evaluation is hindered by the limited availability of sequences. Eight diverse plastomes were sequenced to assess variation in structure and functional content across lycophytes. Lycopodiaceae plastomes have remained largely unchanged compared with the common ancestor of land plants, whereas plastome evolution in Isoetes and especially Selaginella is highly dynamic. Selaginella plastomes have the highest GC content and fewest genes and introns of any photosynthetic land plant. Uniquely, the canonical inverted repeat was converted into a direct repeat (DR) via large-scale inversion in some Selaginella species. Ancestral reconstruction identified additional putative transitions between an inverted and DR orientation in Selaginella and Isoetes plastomes. A DR orientation does not disrupt the activity of copy-dependent repair to suppress substitution rates within repeats. Lycophyte plastomes include the most archaic examples among vascular plants and the most reconfigured among land plants. These evolutionary trends correlate with the mitochondrial genome, suggesting shared underlying mechanisms. Copy-dependent repair for DR-localized genes indicates that recombination and gene conversion are not inhibited by the DR orientation. Gene relocation in lycophyte plastomes occurs via overlapping inversions rather than transposase/recombinase-mediated processes.

Thirty clues to the exceptional diversification of flowering plants

BACKGROUND AND AIMS

As angiosperms became one of the megadiverse groups of macroscopic eukaryotes, they forged modern ecosystems and promoted the evolution of extant terrestrial biota. Unequal distribution of species among lineages suggests that diversification, the process that ultimately determines species richness, acted differentially through angiosperm evolution.

METHODS

We investigate how angiosperms became megadiverse by identifying the phylogenetic and temporal placement of exceptional radiations, by combining the most densely fossil-calibrated molecular clock phylogeny with a Bayesian model that identifies diversification shifts among evolutionary lineages and through time. We evaluate the effect of the prior number of expected shifts in the phylogenetic tree.

KEY RESULTS

Major diversification increases took place over 100 Ma, from the Early Cretaceous to the end of the Paleogene, and are distributed across the angiosperm phylogeny. The long-term diversification trajectory of angiosperms shows moderate rate variation, but is underlain by increasing speciation and extinction, and results from temporally overlapping, independent radiations and depletions in component lineages.

CONCLUSIONS

The identified deep time diversification shifts are clues to the identification of ultimate drivers of angiosperm megadiversity, which probably involve multivariate interactions among intrinsic traits and extrinsic forces. An enhanced understanding of angiosperm diversification will involve a more precise phylogenetic location of diversification shifts, and integration of fossil information.

Contemporaneous radiations of fungi and plants linked to symbiosis

Interactions between fungi and plants, including parasitism, mutualism, and saprotrophy, have been invoked as key to their respective macroevolutionary success. Here we evaluate the origins of plant-fungal symbioses and saprotrophy using a time-calibrated phylogenetic framework that reveals linked and drastic shifts in diversification rates of each kingdom. Fungal colonization of land was associated with at least two origins of terrestrial green algae and preceded embryophytes (as evidenced by losses of fungal flagellum, ca. 720 Ma), likely facilitating terrestriality through endomycorrhizal and possibly endophytic symbioses. The largest radiation of fungi (Leotiomyceta), the origin of arbuscular mycorrhizae, and the diversification of extant embryophytes occurred ca. 480 Ma. This was followed by the origin of extant lichens. Saprotrophic mushrooms diversified in the Late Paleozoic as forests of seed plants started to dominate the landscape. The subsequent diversification and explosive radiation of Agaricomycetes, and eventually of ectomycorrhizal mushrooms, were associated with the evolution of Pinaceae in the Mesozoic, and establishment of angiosperm-dominated biomes in the Cretaceous.

Plants and fungi interact widely and in diverse ways, from mutualism to parasitism and decomposition. Here, Lutzoni et al. analyse the timing of plant and fungal evolutionary radiations and identify four major periods in which plant-fungal interactions likely drove lineage diversification.

Characterization of chloroplast iridescence in Selaginella erythropus

Iridescence in shade-dwelling plants has previously been described in only a few plant groups, and even fewer where the structural colour is produced by intracellular structures. In contrast with other Selaginella species, this work reports the first example in the genus of structural colour originating from modified chloroplasts. Characterization of these structures determines that they form one-dimensional photonic multilayers. The Selaginella bizonoplasts present an analogous structure to recently reported Begonia iridoplasts; however, unlike Begonia species that produce iridoplasts, this Selaginella species was not previously described as iridescent. This therefore raises the possibility of widespread but unobserved and uncharacterized photonic structures in plants.

Ant-plant interactions evolved through increasing interdependence

Ant-plant interactions are diverse and abundant and include classic models in the study of mutualism and other biotic interactions. By estimating a time-scaled phylogeny of more than 1,700 ant species and a time-scaled phylogeny of more than 10,000 plant genera, we infer when and how interactions between ants and plants evolved and assess their macroevolutionary consequences. We estimate that ant-plant interactions originated in the Mesozoic, when predatory, ground-inhabiting ants first began foraging arboreally. This served as an evolutionary precursor to the use of plant-derived food sources, a dietary transition that likely preceded the evolution of extrafloral nectaries and elaiosomes. Transitions to a strict, plant-derived diet occurred in the Cenozoic, and optimal models of shifts between strict predation and herbivory include omnivory as an intermediate step. Arboreal nesting largely evolved from arboreally foraging lineages relying on a partially or entirely plant-based diet, and was initiated in the Mesozoic, preceding the evolution of domatia. Previous work has suggested enhanced diversification in plants with specialized ant-associated traits, but it appears that for ants, living and feeding on plants does not affect ant diversification. Together, the evidence suggests that ants and plants increasingly relied on one another and incrementally evolved more intricate associations with different macroevolutionary consequences as angiosperms increased their ecological dominance.

Facets of phylodiversity: evolutionary diversification, divergence and survival as conservation targets

Biodiversity is often described as having multiple facets, including species richness, functional diversity and phylogenetic diversity. In this paper, we argue that phylogenetic diversity itself has three distinct facets-lineage diversification, character divergence and survival time-that can be quantified using distinct branch length metrics on an evolutionary tree. Each dimension is related to different processes of macroevolution, has different spatial patterns and is tied to distinct goals for conserving biodiversity and protecting its future resilience and evolutionary potential. We compared the landscapes identified as top conservation priorities by each of these three metrics in a conservation gap analysis for California, a world biodiversity hotspot, using herbarium data on the biogeography and evolutionary relationships of more than 5000 native plant species. Our analysis incorporated a novel continuous metric of current land protection status, fine-scale data on landscape intactness and an optimization algorithm used to identify complementary priority sites containing concentrations of taxa that are evolutionarily unique, vulnerable due to small range size and/or poorly protected across their ranges. Top conservation priorities included pockets of coastal and northern California that ranked highly for all three phylodiversity dimensions and for species richness, as well as sites uniquely identified by each metric whose value may depend on whether properties such as genetic divergence, high net diversification or independent survival experience are most desirable in an Anthropocene flora.This article is part of the theme issue 'Biological collections for understanding biodiversity in the Anthropocene'.

Spatial phylogenetics reveals evolutionary constraints on the assembly of a large regional flora

PREMISE OF THE STUDY

We used spatial phylogenetics to analyze the assembly of the Wisconsin flora, linking processes of dispersal and niche evolution to spatial patterns of floristic and phylogenetic diversity and testing whether phylogenetic niche conservatism can account for these patterns.

METHODS

We used digitized records and a new molecular phylogeny for 93% of vascular plants in Wisconsin to estimate spatial variation in species richness and phylogenetic α and β diversity in a native flora shaped mainly by postglacial dispersal and response to environmental gradients. We developed distribution models for all species and used these to infer fine-scale variation in potential diversity, phylogenetic distance, and interspecific range overlaps. We identified 11 bioregions based on floristic composition, mapped areas of neo- and paleo-endemism to establish new conservation priorities and predict how community-assembly patterns should shift with climatic change.

KEY RESULTS

Spatial phylogenetic turnover most strongly reflects differences in temperature and spatial distance. For all vascular plants, assemblages shift from phylogenetically clustered to overdispersed northward, contrary to most other studies. This pattern is lost for angiosperms alone, illustrating the importance of phylogenetic scale.

CONCLUSIONS

Species ranges and assemblage composition appear driven primarily by phylogenetic niche conservatism. Closely related species are ecologically similar and occupy similar territories. The average level and geographic structure of plant phylogenetic diversity within Wisconsin are expected to greatly decline over the next half century, while potential species richness will increase throughout the state. Our methods can be applied to allochthonous communities throughout the world.

Community phylogeny of the globally critically imperiled pine rockland ecosystem

PREMISE OF THE STUDY

Community phylogenetic methods incorporate information on evolutionary relationships into studies of organismal assemblages. We used a community phylogenetic framework to investigate relationships and biogeographic affinities and to calculate phylogenetic signal of endemism and invasiveness for the flora of the pine rocklands-a globally critically imperiled ecosystem with a significant portion of its distribution in South Florida, United States.

METHODS

We reconstructed phylogenetic relationships of 538 vascular plant taxa, which represent 92.28% of the vascular flora of the pine rocklands. We estimated phylogenetic signal for endemism and invasiveness using phylogenetic generalized linear mixed models. We determined the native range for each species in the data set and calculated the total number of species sourced from each region and all possible combinations of these regions.

KEY RESULTS

The pine rockland flora includes representatives of all major vascular plant lineages, and most species have native ranges in the New World. There was strong phylogenetic signal for endemism, but not for invasiveness.

CONCLUSIONS

Community phylogenetics has high potential value for conservation planning, particularly for fragmented and endangered ecosystems like the pine rockland. Strong phylogenetic signal for endemic species in our data set, which also tend to be threatened or endangered, can help to identify species at risk, as well as fragments where those species occur, highlighting conservation priorities. Our results indicate, at least in the pine rockland ecosystem, no phylogenetic signal for invasive species, and thus other information must be used to predict the potential for invasiveness.

The Vegetation Red Edge Biosignature Through Time on Earth and Exoplanets

The high reflection of land vegetation in the near-infrared, the vegetation red edge (VRE), is often cited as a spectral biosignature for surface vegetation on exoplanets. The VRE involves only a few percentage change in reflectivity for a disk-integrated observation of present-day Earth. Here we show that the strength of Earth's VRE has increased over the past ∼500 million years of land plant evolution and may continue to increase as solar luminosity increases and the planet warms, until either vegetation coverage is reduced, or the planet's atmosphere becomes opaque to light reflected off the surface. Early plants such as mosses and liverworts, which dominated the land 500-400 million years ago, produce a weaker VRE, approximately half as strong as that of modern vegetation. We explore how the changes in land plants, as well as geological changes such as ice coverage during ice ages and interglacial periods, influence the detectability of the VRE through Earth's geological past. Our results show that the VRE has varied through the evolutionary history of land plants on Earth and could continue to change into the future if hotter climate conditions became dominant, encouraging the spread of vegetation. Our findings suggest that older and hotter Earth-like planets are good targets for the search for a VRE signature. In addition, hot exoplanets and dry exoplanets with some water could be the best targets for a successful vegetation biosignature detection. As well as a strong red edge, lower cloud fractions and low levels of atmospheric water vapor on such planets could make it easier to detect surface features in general.

Biosynthesis of riccionidins and marchantins is regulated by R2R3-MYB transcription factors in Marchantia polymorpha

R2R3-MYB transcription factors constitute the largest gene family among plant transcription factor families. They became largely divergent during the evolution of land plants and regulate various biological processes. The functions of R2R3-MYBs are mostly characterized in seed plants but are poorly understood in non-seed plants. Here, we examined the function of two R2R3-MYB genes of Marchantia polymorpha (Mapoly0073s0038 and Mapoly0006s0226) that are closely related to subgroup 4 of the R2R3-MYB family. We performed LC/MS/MS metabolomics, RNA-seq analysis and expression analysis in overexpressors and knockout mutants of MpMYB14 and MpMYB02. Overexpression of MpMYB14 remarkably increased the amount of riccionidins, which are specific anthocyanins in liverworts and a few flowering plants. In contrast, overexpression of MpMYB02 increased the amount of several marchantins, which are characteristic cyclic bis (bibenzyl ether) compounds in M. polymorpha and related liverworts. Knockouts of MpMYB14 and MpMYB02 abolished the accumulation of riccionidins and marchantins, respectively. The expression of MpMYB14 was up-regulated by UV-B irradiation, N deficiency, and NaCl treatment, whereas the expression of MpMYB02 was down-regulated by NaCl treatment. Our results suggest that the regulatory framework of phenolic metabolism by R2R3-MYB was already established in early land plants.

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.

Tangled history of a multigene family: The evolution of ISOPENTENYLTRANSFERASE genes

ISOPENTENYLTRANSFERASE (IPT) genes play important roles in the initial steps of cytokinin synthesis, exist in plant and pathogenic bacteria, and form a multigene family in plants. Protein domain searches revealed that bacteria and plant IPT proteins were to assigned to different protein domains families in the Pfam database, namely Pfam IPT (IPT^Pfam) and Pfam IPPT (IPPT^Pfam) families, both are closely related in the P-loop NTPase clan. To understand the origin and evolution of the genes, a species matrix was assembled across the tree of life and intensively in plant lineages. The IPT^Pfam domain was only found in few bacteria lineages, whereas IPPT^Pfam is common except in Archaea and Mycoplasma bacteria. The bacterial IPPT^Pfam domain miaA genes were shown as ancestral of eukaryotic IPPT^Pfam domain genes. Plant IPTs diversified into class I, class II tRNA-IPTs, and Adenosine-phosphate IPTs; the class I tRNA-IPTs appeared to represent direct successors of miaA genes were found in all plant genomes, whereas class II tRNA-IPTs originated from eukaryotic genes, and were found in prasinophyte algae and in euphyllophytes. Adenosine-phosphate IPTs were only found in angiosperms. Gene duplications resulted in gene redundancies with ubiquitous expression or diversification in expression. In conclusion, it is shown that IPT genes have a complex history prior to the protein family split, and might have experienced losses or HGTs, and gene duplications that are to be likely correlated with the rise in morphological complexity involved in fine tuning cytokinin production.

A well-resolved fern nuclear phylogeny reveals the evolution history of numerous transcription factor families

Ferns account for 80% of nonflowering vascular plant species and are the sister lineage of seed plants. Recent molecular phylogenetics have greatly advanced understanding of fern tree of life, but relationships among some major lineages remain unclear. To better resolve the phylogenetic relationships of ferns, we generated transcriptomes from 125 ferns and two lycophytes, with three additional public datasets, to represent all 11 orders and 85% of families of ferns. Our nuclear phylogeny provides strong supports for the monophyly of all four subclasses and nearly all orders and families, and for relationships among these lineages. The only exception is Gleicheniales, which was highly supported as being paraphyletic with Dipteridaceae sister to a clade with Gleicheniaceae + Hymenophyllales. In addition, new and strongly supported phylogenetic relationships are found for suborders and families in Polypodiales. We provide the first dated fern phylogenomic tree using many nuclear genes from a large majority of families, with an estimate for separation of the ancestors of ferns and seed plants in early Devonian at ∼400 Mya and subsequent gradual divergences of fern orders from ∼380 to 200 Mya. Moreover, the newly obtained fern phylogeny provides a framework for gene family analyses, which indicate that the vast majority of transcription factor families found in seed plants were already present in the common ancestor of extant vascular plants. In addition, fern transcription factor genes show similar duplication patterns to those in seed plants, with some showing stable copy number and others displaying independent expansions in both ferns and seed plants. This study provides a robust phylogenetic and gene family evolution framework, as well as rich molecular resources for understanding the morphological and functional evolution in ferns.

Animal origins and the Tonian Earth system

The Neoproterozoic Era (1000-541 million years ago, Ma) was characterized by dramatic environmental and evolutionary change, including at least two episodes of extensive, low-latitude glaciation, potential changes in the redox structure of the global ocean, and the origin and diversification of animal life. How these different events related to one another remains an active area of research, particularly how these environmental changes influenced, and were influenced by, the earliest evolution of animals. Animal multicellularity is estimated to have evolved in the Tonian Period (1000-720 Ma) and represents one of at least six independent acquisitions of complex multicellularity, characterized by cellular differentiation, three-dimensional body plans, and active nutrient transport. Compared with the other instances of complex multicellularity, animals represent the only clade to have evolved from wall-less, phagotrophic flagellates, which likely placed unique cytological and trophic constraints on the evolution of animal multicellularity. Here, we compare recent molecular clock estimates with compilations of the chromium isotope, micropaleontological, and organic biomarker records, suggesting that, as of now, the origin of animals was not obviously correlated to any environmental-ecological change in the Tonian Period. This lack of correlation is consistent with the idea that the evolution of animal multicellularity was primarily dictated by internal, developmental constraints and occurred independently of the known environmental-ecological changes that characterized the Neoproterozoic Era.

Evolutionary and genomic analysis of the caleosin/peroxygenase (CLO/PXG) gene/protein families in the Viridiplantae

Bioinformatics analyses of caleosin/peroxygenases (CLO/PXG) demonstrated that these genes are present in the vast majority of Viridiplantae taxa for which sequence data are available. Functionally active CLO/PXG proteins with roles in abiotic stress tolerance and lipid droplet storage are present in some Trebouxiophycean and Chlorophycean green algae but are absent from the small number of sequenced Prasinophyceaen genomes. CLO/PXG-like genes are expressed during dehydration stress in Charophyte algae, a sister clade of the land plants (Embryophyta). CLO/PXG-like sequences are also present in all of the >300 sequenced Embryophyte genomes, where some species contain as many as 10–12 genes that have arisen via selective gene duplication. Angiosperm genomes harbour at least one copy each of two distinct CLO/PX isoforms, termed H (high) and L (low), where H-forms contain an additional C-terminal motif of about 30–50 residues that is absent from L-forms. In contrast, species in other Viridiplantae taxa, including green algae, non-vascular plants, ferns and gymnosperms, contain only one (or occasionally both) of these isoforms per genome. Transcriptome and biochemical data show that CLO/PXG-like genes have complex patterns of developmental and tissue-specific expression. CLO/PXG proteins can associate with cytosolic lipid droplets and/or bilayer membranes. Many of the analysed isoforms also have peroxygenase activity and are involved in oxylipin metabolism. The distribution of CLO/PXG-like genes is consistent with an origin >1 billion years ago in at least two of the earliest diverging groups of the Viridiplantae, namely the Chlorophyta and the Streptophyta, after the Viridiplantae had already diverged from other Archaeplastidal groups such as the Rhodophyta and Glaucophyta. While algal CLO/PXGs have roles in lipid packaging and stress responses, the Embryophyte proteins have a much wider spectrum of roles and may have been instrumental in the colonisation of terrestrial habitats and the subsequent diversification as the major land flora.

Genomic Changes Associated with the Evolutionary Transitions of Nostoc to a Plant Symbiont

Cyanobacteria belonging to the genus Nostoc comprise free-living strains and also facultative plant symbionts. Symbiotic strains can enter into symbiosis with taxonomically diverse range of host plants. Little is known about genomic changes associated with evolutionary transition of Nostoc from free-living to plant symbiont. Here, we compared the genomes derived from 11 symbiotic Nostoc strains isolated from different host plants and infer phylogenetic relationships between strains. Phylogenetic reconstructions of 89 Nostocales showed that symbiotic Nostoc strains with a broad host range, entering epiphytic and intracellular or extracellular endophytic interactions, form a monophyletic clade indicating a common evolutionary history. A polyphyletic origin was found for Nostoc strains which enter only extracellular symbioses, and inference of transfer events implied that this trait was likely acquired several times in the evolution of the Nostocales. Symbiotic Nostoc strains showed enriched functions in transport and metabolism of organic sulfur, chemotaxis and motility, as well as the uptake of phosphate, branched-chain amino acids, and ammonium. The genomes of the intracellular clade differ from that of other Nostoc strains, with a gain/enrichment of genes encoding proteins to generate l-methionine from sulfite and pathways for the degradation of the plant metabolites vanillin and vanillate, and of the macromolecule xylan present in plant cell walls. These compounds could function as C-sources for members of the intracellular clade. Molecular clock analysis indicated that the intracellular clade emerged ca. 600 Ma, suggesting that intracellular Nostoc symbioses predate the origin of land plants and the emergence of their extant hosts.

Ginkgo biloba's footprint of dynamic Pleistocene history dates back only 390,000 years ago

BACKGROUND

At the end of the Pliocene and the beginning of Pleistocene glaciation and deglaciation cycles Ginkgo biloba went extinct all over the world, and only few populations remained in China in relict areas serving as sanctuary for Tertiary relict trees. Yet the status of these regions as refuge areas with naturally existing populations has been proven not earlier than one decade ago. Herein we elaborated the hypothesis that during the Pleistocene cooling periods G. biloba expanded its distribution range in China repeatedly. Whole plastid genomes were sequenced, assembled and annotated, and sequence data was analyzed in a phylogenetic framework of the entire gymnosperms to establish a robust spatio-temporal framework for gymnosperms and in particular for G. biloba Pleistocene evolutionary history.

RESULTS

Using a phylogenetic approach, we identified that Ginkgoatae stem group age is about 325 million years, whereas crown group radiation of extant Ginkgo started not earlier than 390,000 years ago. During repeated warming phases, Gingko populations were separated and isolated by contraction of distribution range and retreated into mountainous regions serving as refuge for warm-temperate deciduous forests. Diversification and phylogenetic splits correlate with the onset of cooling phases when Ginkgo expanded its distribution range and gene pools merged.

CONCLUSIONS

Analysis of whole plastid genome sequence data representing the entire spatio-temporal genetic variation of wild extant Ginkgo populations revealed the deepest temporal footprint dating back to approximately 390,000 years ago. Present-day directional West-East admixture of genetic diversity is shown to be the result of pronounced effects of the last cooling period. Our evolutionary framework will serve as a conceptual roadmap for forthcoming genomic sequence data, which can then provide deep insights into the demographic history of Ginkgo.

The Gastrodia elata genome provides insights into plant adaptation to heterotrophy

We present the 1.06 Gb sequenced genome of Gastrodia elata, an obligate mycoheterotrophic plant, which contains 18,969 protein-coding genes. Many genes conserved in other plant species have been deleted from the G. elata genome, including most of those for photosynthesis. Additional evidence of the influence of genome plasticity in the adaptation of this mycoheterotrophic lifestyle is evident in the large number of gene families that are expanded in G. elata, including glycoside hydrolases and urease that likely facilitate the digestion of hyphae are expanded, as are genes associated with strigolactone signaling, and ATPases that may contribute to the atypical energy metabolism. We also find that the plastid genome of G. elata is markedly smaller than that of green plant species while its mitochondrial genome is one of the largest observed to date. Our report establishes a foundation for studying adaptation to a mycoheterotrophic lifestyle.

The Gastrodia elata genome provides insights into plant adaptation to heterotrophy

We present the 1.06 Gb sequenced genome of Gastrodia elata, an obligate mycoheterotrophic plant, which contains 18,969 protein-coding genes. Many genes conserved in other plant species have been deleted from the G. elata genome, including most of those for photosynthesis. Additional evidence of the influence of genome plasticity in the adaptation of this mycoheterotrophic lifestyle is evident in the large number of gene families that are expanded in G. elata, including glycoside hydrolases and urease that likely facilitate the digestion of hyphae are expanded, as are genes associated with strigolactone signaling, and ATPases that may contribute to the atypical energy metabolism. We also find that the plastid genome of G. elata is markedly smaller than that of green plant species while its mitochondrial genome is one of the largest observed to date. Our report establishes a foundation for studying adaptation to a mycoheterotrophic lifestyle.

Enhancing (crop) plant photosynthesis by introducing novel genetic diversity

Although some elements of the photosynthetic light reactions might appear to be ideal, the overall efficiency of light conversion to biomass has not been optimized during evolution. Because crop plants are depleted of genetic diversity for photosynthesis, efforts to enhance its efficiency with respect to light conversion to yield must generate new variation. In principle, three sources of natural variation are available: (i) rare diversity within extant higher plant species, (ii) photosynthetic variants from algae, and (iii) reconstruction of no longer extant types of plant photosynthesis. Here, we argue for a novel approach that outsources crop photosynthesis to a cyanobacterium that is amenable to adaptive evolution. This system offers numerous advantages, including a short generation time, virtually unlimited population sizes and high mutation rates, together with a versatile toolbox for genetic manipulation. On such a synthetic bacterial platform, 10 000 years of (crop) plant evolution can be recapitulated within weeks. Limitations of this system arise from its unicellular nature, which cannot reproduce all aspects of crop photosynthesis. But successful establishment of such a bacterial host for crop photosynthesis promises not only to enhance the performance of eukaryotic photosynthesis but will also reveal novel facets of the molecular basis of photosynthetic flexibility.This article is part of the themed issue 'Enhancing photosynthesis in crop plants: targets for improvement'.

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.

The timescale of early land plant evolution

Establishing the timescale of early land plant evolution is essential for testing hypotheses on the coevolution of land plants and Earth's System. The sparseness of early land plant megafossils and stratigraphic controls on their distribution make the fossil record an unreliable guide, leaving only the molecular clock. However, the application of molecular clock methodology is challenged by the current impasse in attempts to resolve the evolutionary relationships among the living bryophytes and tracheophytes. Here, we establish a timescale for early land plant evolution that integrates over topological uncertainty by exploring the impact of competing hypotheses on bryophyte-tracheophyte relationships, among other variables, on divergence time estimation. We codify 37 fossil calibrations for Viridiplantae following best practice. We apply these calibrations in a Bayesian relaxed molecular clock analysis of a phylogenomic dataset encompassing the diversity of Embryophyta and their relatives within Viridiplantae. Topology and dataset sizes have little impact on age estimates, with greater differences among alternative clock models and calibration strategies. For all analyses, a Cambrian origin of Embryophyta is recovered with highest probability. The estimated ages for crown tracheophytes range from Late Ordovician to late Silurian. This timescale implies an early establishment of terrestrial ecosystems by land plants that is in close accord with recent estimates for the origin of terrestrial animal lineages. Biogeochemical models that are constrained by the fossil record of early land plants, or attempt to explain their impact, must consider the implications of a much earlier, middle Cambrian-Early Ordovician, origin.

Deconstructing the evolutionary complexity between rust fungi (Pucciniales) and their plant hosts

The rust fungi (Pucciniales) are the most speciose natural group of plant pathogens, members of which possess the most complex lifecycles in Fungi. How natural selection works on the Pucciniales has been the subject of several hypotheses in mycology. This study uses molecular age estimation using sequence data from multiple loci, and cophylogeny reconciliation analyses to test hypotheses regarding how the aecial and telial stages in the lifecycle of rust fungi may have differentially impacted their diversification. Molecular age estimates show that the timing of diversification in the Pucciniales correlates with the diversification of their gymnosperm and angiosperm hosts. Host reconciliation analyses suggest that systematic relationships of hosts from the aecial stage of the Pucciniales lifecycle better reflect the systematic relationships among the Pucciniales. The results demonstrate the relative importance of this stage on the overall evolution of the Pucciniales and supports hypotheses made by Leppik over half a century ago. This study represents the first evaluation of how different life stages in the Pucciniales shape the evolution of these fungi.

An Evaluation of Different Partitioning Strategies for Bayesian Estimation of Species Divergence Times

The explosive growth of molecular sequence data has made it possible to estimate species divergence times under relaxed-clock models using genome-scale data sets with many gene loci. In order to improve both model realism and to best extract information about relative divergence times in the sequence data, it is important to account for the heterogeneity in the evolutionary process across genes or genomic regions. Partitioning is a commonly used approach to achieve those goals. We group sites that have similar evolutionary characteristics into the same partition and those with different characteristics into different partitions, and then use different models or different values of model parameters for different partitions to account for the among-partition heterogeneity. However, how to partition data in practical phylogenetic analysis, and in particular in relaxed-clock dating analysis, is more art than science. Here, we use computer simulation and real data analysis to study the impact of the partition scheme on divergence time estimation. The partition schemes had relatively minor effects on the accuracy of posterior time estimates when the prior assumptions were correct and the clock was not seriously violated, but showed large differences when the clock was seriously violated, when the fossil calibrations were in conflict or incorrect, or when the rate prior was mis-specified. Concatenation produced the widest posterior intervals with the least precision. Use of many partitions increased the precision, as predicted by the infinite-sites theory, but the posterior intervals might fail to include the true ages because of the conflicting fossil calibrations or mis-specified rate priors. We analyzed a data set of 78 plastid genes from 15 plant species with serious clock violation and showed that time estimates differed significantly among partition schemes, irrespective of the rate drift model used. Multiple and precise fossil calibrations reduced the differences among partition schemes and were important to improving the precision of divergence time estimates. While the use of many partitions is an important approach to reducing the uncertainty in posterior time estimates, we do not recommend its general use for the present, given the limitations of current models of rate drift for partitioned data and the challenges of interpreting the fossil evidence to construct accurate and informative calibrations.

The Interrelationships of Land Plants and the Nature of the Ancestral Embryophyte

The evolutionary emergence of land plant body plans transformed the planet. However, our understanding of this formative episode is mired in the uncertainty associated with the phylogenetic relationships among bryophytes (hornworts, liverworts, and mosses) and tracheophytes (vascular plants). Here we attempt to clarify this problem by analyzing a large transcriptomic dataset with models that allow for compositional heterogeneity between sites. Zygnematophyceae is resolved as sister to land plants, but we obtain several distinct relationships between bryophytes and tracheophytes. Concatenated sequence analyses that can explicitly accommodate site-specific compositional heterogeneity give more support for a mosses-liverworts clade, "Setaphyta," as the sister to all other land plants, and weak support for hornworts as the sister to all other land plants. Bryophyte monophyly is supported by gene concatenation analyses using models explicitly accommodating lineage-specific compositional heterogeneity and analyses of gene trees. Both maximum-likelihood analyses that compare the fit of each gene tree to proposed species trees and Bayesian supertree estimation based on gene trees support bryophyte monophyly. Of the 15 distinct rooted relationships for embryophytes, we reject all but three hypotheses, which differ only in the position of hornworts. Our results imply that the ancestral embryophyte was more complex than has been envisaged based on topologies recognizing liverworts as the sister lineage to all other embryophytes. This requires many phenotypic character losses and transformations in the liverwort lineage, diminishes inconsistency between phylogeny and the fossil record, and prompts re-evaluation of the phylogenetic affinity of early land plant fossils, the majority of which are considered stem tracheophytes.

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.].

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.

Historical Biogeography of endemic seed plant genera in the Caribbean: Did GAARlandia play a role?

The Caribbean archipelago is a region with an extremely complex geological history and an outstanding plant diversity with high levels of endemism. The aim of this study was to better understand the historical assembly and evolution of endemic seed plant genera in the Caribbean, by first determining divergence times of endemic genera to test whether the hypothesized Greater Antilles and Aves Ridge (GAARlandia) land bridge played a role in the archipelago colonization and second by testing South America as the main colonization source as expected by the position of landmasses and recent evidence of an asymmetrical biotic interchange. We reconstructed a dated molecular phylogenetic tree for 625 seed plants including 32 Caribbean endemic genera using Bayesian inference and ten calibrations. To estimate the geographic range of the ancestors of endemic genera, we performed a model selection between a null and two complex biogeographic models that included timeframes based on geological information, dispersal probabilities, and directionality among regions. Crown ages for endemic genera ranged from Early Eocene (53.1 Ma) to Late Pliocene (3.4 Ma). Confidence intervals for divergence times (crown and/or stem ages) of 22 endemic genera occurred within the GAARlandia time frame. Contrary to expectations, the Antilles appears as the main ancestral area for endemic seed plant genera and only five genera had a South American origin. In contrast to patterns shown for vertebrates and other organisms and based on our sampling, we conclude that GAARlandia did not act as a colonization route for plants between South America and the Antilles. Further studies on Caribbean plant dispersal at the species and population levels will be required to reveal finer-scale biogeographic patterns and mechanisms.

Spatial phylogenetics of the native California flora

BACKGROUND

California is a world floristic biodiversity hotspot where the terms neo- and paleo-endemism were first applied. Using spatial phylogenetics, it is now possible to evaluate biodiversity from an evolutionary standpoint, including discovering significant areas of neo- and paleo-endemism, by combining spatial information from museum collections and DNA-based phylogenies. Here we used a distributional dataset of 1.39 million herbarium specimens, a phylogeny of 1083 operational taxonomic units (OTUs) and 9 genes, and a spatial randomization test to identify regions of significant phylogenetic diversity, relative phylogenetic diversity, and phylogenetic endemism (PE), as well as to conduct a categorical analysis of neo- and paleo-endemism (CANAPE).

RESULTS

We found (1) extensive phylogenetic clustering in the South Coast Ranges, southern Great Valley, and deserts of California; (2) significant concentrations of short branches in the Mojave and Great Basin Deserts and the South Coast Ranges and long branches in the northern Great Valley, Sierra Nevada foothills, and the northwestern and southwestern parts of the state; (3) significant concentrations of paleo-endemism in Northwestern California, the northern Great Valley, and western Sonoran Desert, and neo-endemism in the White-Inyo Range, northern Mojave Desert, and southern Channel Islands. Multiple analyses were run to observe the effects on significance patterns of using different phylogenetic tree topologies (uncalibrated trees versus time-calibrated ultrametric trees) and using different representations of OTU ranges (herbarium specimen locations versus species distribution models).

CONCLUSIONS

These analyses showed that examining the geographic distributions of branch lengths in a statistical framework adds a new dimension to California floristics that, in comparison with climatic data, helps to illuminate causes of endemism. In particular, the concentration of significant PE in more arid regions of California extends previous ideas about aridity as an evolutionary stimulus. The patterns seen are largely robust to phylogenetic uncertainty and time calibration but are sensitive to the use of occurrence data versus modeled ranges, indicating that special attention toward improving geographic distributional data should be top priority in the future for advancing understanding of spatial patterns of biodiversity.

Contrasting Rates of Molecular Evolution and Patterns of Selection among Gymnosperms and Flowering Plants

The majority of variation in rates of molecular evolution among seed plants remains both unexplored and unexplained. Although some attention has been given to flowering plants, reports of molecular evolutionary rates for their sister plant clade (gymnosperms) are scarce, and to our knowledge differences in molecular evolution among seed plant clades have never been tested in a phylogenetic framework. Angiosperms and gymnosperms differ in a number of features, of which contrasting reproductive biology, life spans, and population sizes are the most prominent. The highly conserved morphology of gymnosperms evidenced by similarity of extant species to fossil records and the high levels of macrosynteny at the genomic level have led scientists to believe that gymnosperms are slow-evolving plants, although some studies have offered contradictory results. Here, we used 31,968 nucleotide sites obtained from orthologous genes across a wide taxonomic sampling that includes representatives of most conifers, cycads, ginkgo, and many angiosperms with a sequenced genome. Our results suggest that angiosperms and gymnosperms differ considerably in their rates of molecular evolution per unit time, with gymnosperm rates being, on average, seven times lower than angiosperm species. Longer generation times and larger genome sizes are some of the factors explaining the slow rates of molecular evolution found in gymnosperms. In contrast to their slow rates of molecular evolution, gymnosperms possess higher substitution rate ratios than angiosperm taxa. Finally, our study suggests stronger and more efficient purifying and diversifying selection in gymnosperm than in angiosperm species, probably in relation to larger effective population sizes.

The Apostasia genome and the evolution of orchids

WebComparing the whole genome sequence of Apostasia shenzhenica with transcriptome and genome data from five orchid subfamilies permits the reconstruction of an ancestral gene toolkit, providing insight into orchid origins, evolution and diversification.

Around 10 per cent of flowering plant species are orchids, with a broad diversity in both morphology and lifestyle. Apostasia is one of the earliest-diverging genera of Orchidaceae. To study the evolution and diversity of Orchidaceae, Zhong-Jian Liu, Yves Van de Peer and colleagues sequenced the genome of Apostasia shenzhenica, a self-pollinating species found in southeast China. The authors also report improved genomes for two species of Epidendroideae, Phalaenopsis equestris and Dendrobium catenatum, as well as transcriptome analysis of representatives of subfamilies of Orchidaceae. Their analyses provide insights into orchid origins, genome evolution, adaptation and diversification.

Constituting approximately 10% of flowering plant species, orchids (Orchidaceae) display unique flower morphologies, possess an extraordinary diversity in lifestyle, and have successfully colonized almost every habitat on Earth^1,2,3. Here we report the draft genome sequence of Apostasia shenzhenica⁴, a representative of one of two genera that form a sister lineage to the rest of the Orchidaceae, providing a reference for inferring the genome content and structure of the most recent common ancestor of all extant orchids and improving our understanding of their origins and evolution. In addition, we present transcriptome data for representatives of Vanilloideae, Cypripedioideae and Orchidoideae, and novel third-generation genome data for two species of Epidendroideae, covering all five orchid subfamilies. A. shenzhenica shows clear evidence of a whole-genome duplication, which is shared by all orchids and occurred shortly before their divergence. Comparisons between A. shenzhenica and other orchids and angiosperms also permitted the reconstruction of an ancestral orchid gene toolkit. We identify new gene families, gene family expansions and contractions, and changes within MADS-box gene classes, which control a diverse suite of developmental processes, during orchid evolution. This study sheds new light on the genetic mechanisms underpinning key orchid innovations, including the development of the labellum and gynostemium, pollinia, and seeds without endosperm, as well as the evolution of epiphytism; reveals relationships between the Orchidaceae subfamilies; and helps clarify the evolutionary history of orchids within the angiosperms.

Evidence for an ancient whole genome duplication in the cycad lineage

Contrary to the many whole genome duplication events recorded for angiosperms (flowering plants), whole genome duplications in gymnosperms (non-flowering seed plants) seem to be much rarer. Although ancient whole genome duplications have been reported for most gymnosperm lineages as well, some are still contested and need to be confirmed. For instance, data for ginkgo, but particularly cycads have remained inconclusive so far, likely due to the quality of the data available and flaws in the analysis. We extracted and sequenced RNA from both the cycad Encephalartos natalensis and Ginkgo biloba. This was followed by transcriptome assembly, after which these data were used to build paralog age distributions. Based on these distributions, we identified remnants of an ancient whole genome duplication in both cycads and ginkgo. The most parsimonious explanation would be that this whole genome duplication event was shared between both species and had occurred prior to their divergence, about 300 million years ago.

Evolutionary Cell Biology of Proteins from Protists to Humans and Plants

During evolution, the cell as a fine-tuned machine had to undergo permanent adjustments to match changes in its environment, while "closed for repair work" was not possible. Evolution from protists (protozoa and unicellular algae) to multicellular organisms may have occurred in basically two lineages, Unikonta and Bikonta, culminating in mammals and angiosperms (flowering plants), respectively. Unicellular models for unikont evolution are myxamoebae (Dictyostelium) and increasingly also choanoflagellates, whereas for bikonts, ciliates are preferred models. Information accumulating from combined molecular database search and experimental verification allows new insights into evolutionary diversification and maintenance of genes/proteins from protozoa on, eventually with orthologs in bacteria. However, proteins have rarely been followed up systematically for maintenance or change of function or intracellular localization, acquirement of new domains, partial deletion (e.g. of subunits), and refunctionalization, etc. These aspects are discussed in this review, envisaging "evolutionary cell biology." Protozoan heritage is found for most important cellular structures and functions up to humans and flowering plants. Examples discussed include refunctionalization of voltage-dependent Ca²⁺ channels in cilia and replacement by other types during evolution. Altogether components serving Ca²⁺ signaling are very flexible throughout evolution, calmodulin being a most conservative example, in contrast to calcineurin whose catalytic subunit is lost in plants, whereas both subunits are maintained up to mammals for complex functions (immune defense and learning). Domain structure of R-type SNAREs differs in mono- and bikonta, as do Ca²⁺ -dependent protein kinases. Unprecedented selective expansion of the subunit a which connects multimeric base piece and head parts (V0, V1) of H⁺ -ATPase/pump may well reflect the intriguing vesicle trafficking system in ciliates, specifically in Paramecium. One of the most flexible proteins is centrin when its intracellular localization and function throughout evolution is traced. There are many more examples documenting evolutionary flexibility of translation products depending on requirements and potential for implantation within the actual cellular context at different levels of evolution. From estimates of gene and protein numbers per organism, it appears that much of the basic inventory of protozoan precursors could be transmitted to highest eukaryotic levels, with some losses and also with important additional "inventions."

The ancestral flower of angiosperms and its early diversification

Recent advances in molecular phylogenetics and a series of important palaeobotanical discoveries have revolutionized our understanding of angiosperm diversification. Yet, the origin and early evolution of their most characteristic feature, the flower, remains poorly understood. In particular, the structure of the ancestral flower of all living angiosperms is still uncertain. Here we report model-based reconstructions for ancestral flowers at the deepest nodes in the phylogeny of angiosperms, using the largest data set of floral traits ever assembled. We reconstruct the ancestral angiosperm flower as bisexual and radially symmetric, with more than two whorls of three separate perianth organs each (undifferentiated tepals), more than two whorls of three separate stamens each, and more than five spirally arranged separate carpels. Although uncertainty remains for some of the characters, our reconstruction allows us to propose a new plausible scenario for the early diversification of flowers, leading to new testable hypotheses for future research on angiosperms.