Search

Accurately reconstructing the reticulate histories of polyploids remains a central challenge for understanding plant evolution. Although phylogenetic networks can provide insights into relationships among polyploid lineages, inferring networks may be hindered by the complexities of homology determination in polyploid taxa. We use simulations to show that phasing alleles from allopolyploid individuals can improve phylogenetic network inference under the multispecies coalescent by obtaining the true network with fewer loci compared with haplotype consensus sequences or sequences with heterozygous bases represented as ambiguity codes. Phased allelic data can also improve divergence time estimates for networks, which is helpful for evaluating allopolyploid speciation hypotheses and proposing mechanisms of speciation. To achieve these outcomes in empirical data, we present a novel pipeline that leverages a recently developed phasing algorithm to reliably phase alleles from polyploids. This pipeline is especially appropriate for target enrichment data, where the depth of coverage is typically high enough to phase entire loci. We provide an empirical example in the North American Dryopteris fern complex that demonstrates insights from phased data as well as the challenges of network inference. We establish that our pipeline (PATÉ: Phased Alleles from Target Enrichment data) is capable of recovering a high proportion of phased loci from both diploids and polyploids. These data may improve network estimates compared with using haplotype consensus assemblies by accurately inferring the direction of gene flow, but statistical nonidentifiability of phylogenetic networks poses a barrier to inferring the evolutionary history of reticulate complexes., (© The Author(s) 2024. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For commercial re-use, please contact reprints@oup.com for reprints and translation rights for reprints. All other permissions can be obtained through our RightsLink service via the Permissions link on the article page on our site—for further information please contact journals.permissions@oup.com.)

Premise: Reticulate evolution, often accompanied by polyploidy, is prevalent in plants, and particularly in the ferns. Resolving the resulting non-bifurcating histories remains a major challenge for plant phylogenetics. Here, we present a phylogenomic investigation into the complex evolutionary history of the vining ferns, Lygodium (Lygodiaceae, Schizaeales)., Methods: Using a targeted enrichment approach with the GoFlag 408 flagellate land plant probe set, we generated large nuclear and plastid sequence datasets for nearly all taxa in the genus and constructed the most comprehensive phylogeny of the family to date using concatenated maximum likelihood and coalescence approaches. We integrated this phylogeny with cytological and spore data to explore karyotype evolution and generate hypotheses about the origins of putative polyploids and hybrids., Results: Our data and analyses support the origins of several putative allopolyploids (e.g., L. cubense, L. heterodoxum) and hybrids (e.g., L. ×fayae) and also highlight the potential prevalence of autopolyploidy in this clade (e.g., L. articulatum, L. flexuosum, and L. longifolium)., Conclusions: Our robust phylogenetic framework provides valuable insights into dynamic reticulate evolution in this clade and demonstrates the utility of target-capture data for resolving these complex relationships., (© 2024 The Author(s). American Journal of Botany published by Wiley Periodicals LLC on behalf of Botanical Society of America.)

A long-standing hypothesis in evolutionary biology is that the evolution of resource specialization can lead to an evolutionary dead end, where specialists have low diversification rates and limited ability to evolve into generalists. In recent years, advances in comparative methods investigating trait-based differences associated with diversification have enabled more robust tests of this idea and have found mixed support. We test the evolutionary dead end hypothesis by estimating net diversification rate differences associated with nest-type specialization among 3224 species of passerine birds. In particular, we test whether the adoption of hole-nesting, a nest-type specialization that decreases predation, results in reduced diversification rates relative to nesting outside of holes. Further, we examine whether evolutionary transitions to the specialist hole-nesting state have been more frequent than transitions out of hole-nesting. Using diversification models that accounted for background rate heterogeneity and different extinction rate scenarios, we found that hole-nesting specialization was not associated with diversification rate differences. Furthermore, contrary to the assumption that specialists rarely evolve into generalists, we found that transitions out of hole-nesting occur more frequently than transitions into hole-nesting. These results suggest that interspecific competition may limit adoption of hole-nesting, but that such competition does not result in limited diversification of hole-nesters. In conjunction with other recent studies using robust comparative methods, our results add to growing evidence that evolutionary dead ends are not a typical outcome of resource specialization. [Cavity nesting; diversification; hidden-state models; passerines; resource specialization.]., (© The Author(s) 2022. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For permissions, please email: journals.permissions@oup.com.)

The Indo-Pacific legume genus Serianthes was recently placed in the Archidendron clade (sensu Koenen et al. 2020), a subclade of the mimosoid clade in subfamily Caesalpinioideae, which also includes Acacia , Archidendron , Archidendropsis , Falcataria , Pararchidendron , Paraserianthes and Wallaceodendron . Serianthes comprises ca. 18 species, five subspecies and two varieties that are characterised by bipinnately compound leaves with alternate sessile leaflets, branched axillary corymbiform panicles and woody indehiscent pods. Generic relationships, as well as species relationships within genera in the Archidendron clade, remain uncertain. While the sister relationship between Serianthes and the genus Falcataria is strongly supported by molecular data, the distinction between Serianthes and the monotypic genus Wallaceodendron has been questioned, based on their similar flower and fruit morphologies. We combined three gene-enriched hybrid capture DNA sequence datasets (generated from the 964 mimobaits v1 probe set, the expanded 997 mimobaits v2 probe set and the GoFlag angiosperm 408 probe set) and used their overlapping markers (77 loci of the target exonic and flanking regions) to test the monophyly of Serianthes and to investigate generic relationships within the Archidendron clade using 55 ingoid plus two outgroup taxa. We show that Serianthes is monophyletic, confirm the Serianthes + Falcataria sister relationship to Wallaceodendron and recognise this combined clade as the Serianthes clade within the Archidendron clade. We also evaluated the use of overlapping loci across datasets in combination with concordance analyses to test generic relationships and further investigate previously unresolved relationships across the wider ingoid clade. Concordance analysis revealed limited gene tree conflicts near the tips of the Archidendron clade, but increased discordance at the base of the clade, which could be attributed to rapid lineage divergence (radiation) and/or incomplete lineage sorting., (Else Demeulenaere, Tom Schils, J. Gordon Burleigh, Jens J. Ringelberg, Erik J. M. Koenen, Stefanie M. Ickert-Bond.)

Phylogenomic data increase the possibilities of resolving the evolutionary and systematic relationships among taxa. This is especially valuable in groups with few and homoplasious morphological characters, in which systematic and taxonomical delimitations have been traditionally difficult. Such is the case of several lineages within Bryophyta, like Orthotrichaceae, the second most diverse family of mosses. Members of tribe Orthotricheae are common in temperate and cold regions, as well as in high tropical mountains. In extratropical areas, they represent one of the main components of epiphytic communities, both in dry and oceanic or hyperoceanic conditions. The epiphytic environment is considered a hostile one for plant development, mainly due to its low capacity of moisture retention. Thus, the diversification of the Orthotrichaceae in this environment could be seen as striking. Over the last two decades, great taxonomic and systematic progresses have led to a rearrangement at the generic level in this tribe, providing a new framework to link environment to patterns of diversification. Here, we use nuclear loci targeted with the GoFlag 408 enrichment probe set to generate a well-sampled phylogeny with well-supported suprageneric taxa and increasing the phylogenetic resolution within the two recognized subtribes. Specifically, we show that several genera with Ulota -like morphology jointly constitute an independent lineage. Within this lineage, the recently described Atlantichella from Macaronesia and Western Europe appears as the sister group of Ulota bellii from Zealandia. This latter species is here segregated in the new genus Rehubryum . Assessment of the ecological and biogeographical affinities of the species within the phylogenetic framework suggests that niche adaptation (including climate and substrate) may be a key evolutionary driver that shaped the high diversification of Orthotricheae., Competing Interests: The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest., (Copyright © 2022 Draper, Villaverde, Garilleti, Burleigh, McDaniel, Mazimpaka, Calleja and Lara.)

Premise: To support large leaves, many woody plant species evolved a cost-effective way to thicken twigs. As an extension of E. J. H. Corner's rule that twig diameter increases with leaf size, we hypothesized that pith width also increases with leaf size. The benefit to the plant from the proposed relationship is that pith is a low-cost tissue that reduces the metabolic cost of large diameter twig production., Methods: Leaf sizes and cross-sectional areas of bark, xylem, and pith of 81 species of trees and shrubs growing in Gainesville, Florida were measured and compared with standardized major axis regressions of pairwise species trait values and phylogenetically independent contrasts., Results: Pith area increases with leaf size with or without accounting for phylogenetic relationships. In agreement with Corner's rule, overall twig diameter as well as bark and wood thickness also increase with leaf size. Thicker twigs showed more variation in relative pith, wood, and bark cross-sectional areas compared to thinner twigs., Conclusions: Investments in pith, a tissue of low density found in the centers of twigs, provides a low-cost way to increase twig circumference and thereby space for attachment of large leaves while increasing the overall second moment of area of twigs, which increases their ability to biomechanically support large leaves., (© 2021 Botanical Society of America.)

Nonrecombining sex chromosomes, like the mammalian Y, often lose genes and accumulate transposable elements, a process termed degeneration. The correlation between suppressed recombination and degeneration is clear in animal XY systems, but the absence of recombination is confounded with other asymmetries between the X and Y. In contrast, UV sex chromosomes, like those found in bryophytes, experience symmetrical population genetic conditions. Here, we generate nearly gapless female and male chromosome-scale reference genomes of the moss Ceratodon purpureus to test for degeneration in the bryophyte UV sex chromosomes. We show that the moss sex chromosomes evolved over 300 million years ago and expanded via two chromosomal fusions. Although the sex chromosomes exhibit weaker purifying selection than autosomes, we find that suppressed recombination alone is insufficient to drive degeneration. Instead, the U and V sex chromosomes harbor thousands of broadly expressed genes, including numerous key regulators of sexual development across land plants., (Copyright © 2021 The Authors, some rights reserved; exclusive licensee American Association for the Advancement of Science. No claim to original U.S. Government Works. Distributed under a Creative Commons Attribution NonCommercial License 4.0 (CC BY-NC).)

Premise: New sequencing technologies facilitate the generation of large-scale molecular data sets for constructing the plant tree of life. We describe a new probe set for target enrichment sequencing to generate nuclear sequence data to build phylogenetic trees with any flagellate land plants, including hornworts, liverworts, mosses, lycophytes, ferns, and all gymnosperms., Methods: We leveraged existing transcriptome and genome sequence data to design the GoFlag 451 probes, a set of 56,989 probes for target enrichment sequencing of 451 exons that are found in 248 single-copy or low-copy nuclear genes across flagellate plant lineages., Results: Our results indicate that target enrichment using the GoFlag451 probe set can provide large nuclear data sets that can be used to resolve relationships among both distantly and closely related taxa across the flagellate land plants. We also describe the GoFlag 408 probes, an optimized probe set covering 408 of the 451 exons from the GoFlag 451 probe set that is commercialized by RAPiD Genomics., Conclusions: A target enrichment approach using the new probe set provides a relatively low-cost solution to obtain large-scale nuclear sequence data for inferring phylogenetic relationships across flagellate land plants., (© 2021 Breinholt et al. Applications in Plant Sciences is published by Wiley Periodicals, Inc. on behalf of the Botanical Society of America.)

Hybrid zones occur as range boundaries for many animal taxa. One model for how hybrid zones form and stabilize is the tension zone model, a version of which predicts that hybrid zone widths are determined by a balance between random dispersal into hybrid zones and selection against hybrids. Here, we examine whether random dispersal and proxies for selection against hybrids (genetic distances between hybridizing pairs) can explain variation in hybrid zone widths across 131 hybridizing pairs of animals. We show that these factors alone can explain ∼40% of the variation in zone width among animal hybrid zones, with dispersal explaining far more of the variation than genetic distances. Patterns within clades were idiosyncratic. Genetic distances predicted hybrid zone widths particularly well for reptiles, while this relationship was opposite tension zone predictions in birds. Last, the data suggest that dispersal and molecular divergence set lower bounds on hybrid zone widths in animals, indicating that there are geographic restrictions on hybrid zone formation. Overall, our analyses reinforce the fundamental importance of dispersal in hybrid zone formation and more generally in the ecology of range boundaries.

Relative brain sizes in birds can rival those of primates, but large-scale patterns and drivers of avian brain evolution remain elusive. Here, we explore the evolution of the fundamental brain-body scaling relationship across the origin and evolution of birds. Using a comprehensive dataset sampling> 2,000 modern birds, fossil birds, and theropod dinosaurs, we infer patterns of brain-body co-variation in deep time. Our study confirms that no significant increase in relative brain size accompanied the trend toward miniaturization or evolution of flight during the theropod-bird transition. Critically, however, theropods and basal birds show weaker integration between brain size and body size, allowing for rapid changes in the brain-body relationship that set the stage for dramatic shifts in early crown birds. We infer that major shifts occurred rapidly in the aftermath of the Cretaceous-Paleogene mass extinction within Neoaves, in which multiple clades achieved higher relative brain sizes because of a reduction in body size. Parrots and corvids achieved the largest brains observed in birds via markedly different patterns. Parrots primarily reduced their body size, whereas corvids increased body and brain size simultaneously (with rates of brain size evolution outpacing rates of body size evolution). Collectively, these patterns suggest that an early adaptive radiation in brain size laid the foundation for subsequent selection and stabilization., Competing Interests: Declaration of Interests The authors declare no competing interests., (Copyright © 2020 Elsevier Inc. All rights reserved.)

The radiation of angiosperms led to the emergence of the vast majority of today's plant species and all our major food crops. Their extraordinary diversification occurred in conjunction with the evolution of a more efficient vascular system for the transport of water, composed of vessel elements. The physical dimensions of these water-conducting specialized cells have played a critical role in angiosperm evolution; they determine resistance to water flow, influence photosynthesis rate, and contribute to plant stature. However, the genetic factors that determine their dimensions are unclear. Here we show that a previously uncharacterized gene, ENLARGED VESSEL ELEMENT ( EVE ) , contributes to the dimensions of vessel elements in Populus , impacting hydraulic conductivity. Our data suggest that EVE is localized in the plasma membrane and is involved in potassium uptake of differentiating xylem cells during vessel development. In plants, EVE first emerged in streptophyte algae, but expanded dramatically among vessel-containing angiosperms. The phylogeny, structure and composition of EVE indicates that it may have been involved in an ancient horizontal gene-transfer event., Competing Interests: Competing interest statement: The use of the gene described in the paper to improve plant growth and productivity is the subject of a patent (US# 9,650,646) and is currently licensed., (Copyright © 2020 the Author(s). Published by PNAS.)

If particular traits consistently affect rates of speciation and extinction, broad macroevolutionary patterns can be interpreted as consequences of selection at high levels of the biological hierarchy. Identifying traits associated with diversification rates is difficult because of the wide variety of characters under consideration and the statistical challenges of testing for associations from comparative phylogenetic data. Ploidy (diploid vs polyploid states) and breeding system (self-incompatible vs self-compatible states) are both thought to be drivers of differential diversification in angiosperms. We fit 29 diversification models to extensive trait and phylogenetic data in Solanaceae and investigate how speciation and extinction rate differences are associated with ploidy, breeding system, and the interaction between these traits. We show that diversification patterns in Solanaceae are better explained by breeding system and an additional unobserved factor, rather than by ploidy. We also find that the most common evolutionary pathway to polyploidy in Solanaceae occurs via direct breakdown of self-incompatibility by whole genome duplication, rather than indirectly via breakdown followed by polyploidization. Comparing multiple stochastic diversification models that include complex trait interactions alongside hidden states enhances our understanding of the macroevolutionary patterns in plant phylogenies., (© 2019 The Authors. New Phytologist © 2019 New Phytologist Trust.)

Genomic data have provided evidence of previously unknown ancient whole genome duplications (WGDs) and highlighted the role of WGDs in the evolution of many eukaryotic lineages. Ancient WGDs often are detected by examining distributions of synonymous substitutions per site (Ks) within a genome, or "Ks plots." For example, WGDs can be detected from Ks plots by using univariate mixture models to identify peaks in Ks distributions. We performed gene family simulation experiments to evaluate the effects of different Ks estimation methods and mixture models on our ability to detect ancient WGDs from Ks plots. The simulation experiments, which accounted for variation in substitution rates and gene duplication and loss rates across gene families, tested the effects of WGD age and gene retention rates following WGD on inferring WGDs from Ks plots. Our simulations reveal limitations of Ks plot analyses. Strict interpretations of mixture model analyses often overestimate the number of WGD events, and Ks plot analyses typically fail to detect WGDs when ≤10% of the duplicated genes are retained following the WGD. However, WGDs can accurately be characterized over an intermediate range of Ks. The simulation results are supported by empirical analyses of transcriptomic data, which also suggest that biases in gene retention likely affect our ability to detect ancient WGDs. Although our results indicate mixture model results should be interpreted with great caution, using node-averaged Ks estimates and applying more appropriate mixture models can improve the accuracy of detecting WGDs.

Organismal traits may evolve either gradually or in rapid pulses, but the relative importance of these modes in the generation of species differences is unclear. Additionally, while pulsed evolution is frequently assumed to be associated with speciation events, few studies have explicitly examined how the tempo of trait divergence varies with respect to different geographical phases of speciation, starting with geographic isolation and ending, in many cases, with spatial overlap (sympatry). Here we address these issues by combining divergence time estimates, trait measurements and geographic range data for 952 avian sister species pairs worldwide to examine the tempo and timing of trait divergence in recent speciation events. We show that patterns of divergence in key ecological traits are not gradual, but instead seem to follow a pattern of relative stasis interspersed with evolutionary pulses of varying magnitude. We also find evidence that evolutionary pulses generally precede sympatry, and that greater trait disparity is associated with sympatry. These findings suggest that early pulses of trait divergence promote subsequent transitions to sympatry, rather than occurring after sympatry has been established. Incorporating models with evolutionary pulses of varying magnitude into speciation theory may explain why some species pairs achieve rapid sympatry whereas others undergo prolonged geographical exclusion.

Premise of the Study: Polyploidy has profound evolutionary consequences for land plants. Despite the availability of large phylogenetic and chromosomal data sets, estimating the rates of polyploidy and chromosomal evolution across the tree of life remains a challenging, computationally complex problem. We introduce the R package chromploid, which allows scientists to perform inference of chromosomal evolution rates across large phylogenetic trees., Methods and Results: chromploid is an open-source package in the R environment that calculates the likelihood function of models of chromosome evolution. Models of discrete character evolution can be customized using chromploid. We demonstrate the performance of the BiChroM model, testing for associations between rates of chromosome doubling (as a proxy for polyploidy) and a binary phenotypic character, within chromploid using simulations and empirical data from Solanum . In simulations, estimated chromosome-doubling rates were unbiased and the variance decreased with larger trees, but distinguishing small differences in rates of chromosome doubling, even from large data sets, remains challenging. In the Solanum data set, a custom model of chromosome number evolution demonstrated higher rates of chromosome doubling in herbaceous species compared to woody., Conclusions: chromploid enables researchers to perform robust likelihood-based inferences using complex models of chromosome number evolution across large phylogenies.

Premise of the Study: Phenotypic data sets are necessary to elucidate the genealogy of life, but assembling phenotypic data for taxa across the tree of life can be technically challenging and prohibitively time consuming. We describe a semi-automated protocol to facilitate and expedite the assembly of phenotypic character matrices of plants from formal taxonomic descriptions. This pipeline uses new natural language processing (NLP) techniques and a glossary of over 9000 botanical terms., Methods and Results: Our protocol includes the Explorer of Taxon Concepts (ETC), an online application that assembles taxon-by-character matrices from taxonomic descriptions, and MatrixConverter, a Java application that enables users to evaluate and discretize the characters extracted by ETC. We demonstrate this protocol using descriptions from Araucariaceae., Conclusions: The NLP pipeline unlocks the phenotypic data found in taxonomic descriptions and makes them usable for evolutionary analyses.

Providing science and society with an integrated, up-to-date, high quality, open, reproducible and sustainable plant tree of life would be a huge service that is now coming within reach. However, synthesizing the growing body of DNA sequence data in the public domain and disseminating the trees to a diverse audience are often not straightforward due to numerous informatics barriers. While big synthetic plant phylogenies are being built, they remain static and become quickly outdated as new data are published and tree-building methods improve. Moreover, the body of existing phylogenetic evidence is hard to navigate and access for non-experts. We propose that our community of botanists, tree builders, and informaticians should converge on a modular framework for data integration and phylogenetic analysis, allowing easy collaboration, updating, data sourcing and flexible analyses. With support from major institutions, this pipeline should be re-run at regular intervals, storing trees and their metadata long-term. Providing the trees to a diverse global audience through user-friendly front ends and application development interfaces should also be a priority. Interactive interfaces could be used to solicit user feedback and thus improve data quality and to coordinate the generation of new data. We conclude by outlining a number of steps that we suggest the scientific community should take to achieve global phylogenetic synthesis., (© 2018 Botanical Society of America.)

Premise of the Study: Many ecological and evolutionary processes shape the assembly of organisms into local communities from a regional pool of species. We analyzed phylogenetic and functional diversity to understand community assembly of the ferns of Florida at two spatial scales., Methods: We built a phylogeny for 125 of the 141 species of ferns in Florida using five chloroplast markers. We calculated mean pairwise dissimilarity (MPD) and mean nearest taxon distance (MNTD) from phylogenetic distances and functional trait data for both spatial scales and compared the results to null models to assess significance., Key Results: Our results for over vs. underdispersion in functional and phylogenetic diversity differed depending on spatial scale and metric considered. At the county scale, MPD revealed evidence for phylogenetic overdispersion, while MNTD revealed phylogenetic and functional underdispersion, and at the conservation area scale, MPD revealed phylogenetic and functional underdispersion while MNTD revealed evidence only of functional underdispersion., Conclusions: Our results are consistent with environmental filtering playing a larger role at the smaller, conservation area scale. The smaller spatial units are likely composed of fewer local habitat types that are selecting for closely related species, with the larger-scale units more likely to be composed of multiple habitat types that bring together a larger pool of species from across the phylogeny. Several aspects of fern biology, including their unique physiology and water relations and the importance of the independent gametophyte stage of the life cycle, make ferns highly sensitive to local, microhabitat conditions., (© 2018 The Authors. American Journal of Botany is published by Wiley Periodicals, Inc. on behalf of the Botanical Society of America.)

Although numerous studies have surveyed the frequency with which different plant characters are associated with polyploidy, few statistical tools are available to identify the factors that potentially facilitate polyploidy. We describe a new probabilistic model, BiChroM, designed to associate the frequency of polyploidy and chromosomal change with a binary phenotypic character in a phylogeny. BiChroM provides a robust statistical framework for testing differences in rates of polyploidy associated with phenotypic characters along a phylogeny while simultaneously allowing for evolutionary transitions between character states. We used BiChroM to test whether polyploidy is more frequent in woody or herbaceous plants, based on tree with 4711 eudicot species. Although polyploidy occurs in woody species, rates of chromosome doubling were over six times higher in herbaceous species. Rates of single chromosome increases or decreases were also far higher in herbaceous than woody species. Simulation experiments indicate that BiChroM performs well with little to no bias and relatively little variance at a wide range of tree depths when trees have at least 500 taxa. Thus, BiChroM provides a first step toward a rigorous statistical framework for assessing the traits that facilitate polyploidy., (© 2017 The Author(s). Evolution © 2017 The Society for the Study of Evolution.)

Phylogenetic trees can reveal the origins of endosymbiotic lineages of bacteria and detect patterns of co-evolution with their hosts. Although taxon sampling can greatly affect phylogenetic and co-evolutionary inference, most hypotheses of endosymbiont relationships are based on few available bacterial sequences. Here we examined how different sampling strategies of Gammaproteobacteria sequences affect estimates of the number of endosymbiont lineages in parasitic sucking lice (Insecta: Phthirapatera: Anoplura). We estimated the number of louse endosymbiont lineages using both newly obtained and previously sequenced 16S rDNA bacterial sequences and more than 42,000 16S rDNA sequences from other Gammaproteobacteria. We also performed parametric and nonparametric bootstrapping experiments to examine the effects of phylogenetic error and uncertainty on these estimates. Sampling of 16S rDNA sequences affects the estimates of endosymbiont diversity in sucking lice until we reach a threshold of genetic diversity, the size of which depends on the sampling strategy. Sampling by maximizing the diversity of 16S rDNA sequences is more efficient than randomly sampling available 16S rDNA sequences. Although simulation results validate estimates of multiple endosymbiont lineages in sucking lice, the bootstrap results suggest that the precise number of endosymbiont origins is still uncertain.

Premise of the Study: Whole-genome duplications (WGDs) can rapidly increase genome size in angiosperms. Yet their mean genome size is not correlated with ploidy. We compared three hypotheses to explain the constancy of genome size means across ploidies. The genome downsizing hypothesis suggests that genome size will decrease by a given percentage after a WGD. The genome size threshold hypothesis assumes that taxa with large genomes or large monoploid numbers will fail to undergo or survive WGDs. Finally, the genome downsizing and threshold hypothesis suggests that both genome downsizing and thresholds affect the relationship between genome size means and ploidy., Methods: We performed nonparametric bootstrap simulations to compare observed angiosperm genome size means among species or genera against simulated genome sizes under the three different hypotheses. We evaluated the hypotheses using a decision theory approach and estimated the expected percentage of genome downsizing., Key Results: The threshold hypothesis improves the approximations between mean genome size and simulated genome size. At the species level, the genome downsizing with thresholds hypothesis best explains the genome size means with a 15% genome downsizing percentage. In the genus level simulations, the monoploid number threshold hypothesis best explains the data., Conclusions: Thresholds of genome size and monoploid number added to genome downsizing at species level simulations explain the observed means of angiosperm genome sizes, and monoploid number is important for determining the genome size mean at the genus level., (© 2016 Botanical Society of America.)

Whole-genome duplications (WGDs) have helped shape the genomes of land plants, and recent evidence suggests that the genomes of all angiosperms have experienced at least two ancient WGDs. In plants, WGDs often are followed by rapid fractionation, in which many homeologous gene copies are lost. Thus, it can be extremely difficult to identify, let alone characterize, ancient WGDs. In this study, we use a new maximum likelihood estimator to test for evidence of ancient WGDs in land plants and estimate the fraction of new genes copies that are retained following a WGD using gene count data, the number of gene copies in gene families. We identified evidence of many putative ancient WGDs in land plants and found that the genome fractionation rates vary tremendously among ancient WGDs. Analyses of WGDs within Brassicales also indicate that background gene duplication and loss rates vary across land plants, and different gene families have different probabilities of being retained following a WGD. Although our analyses are largely robust to errors in duplication and loss rates and the choice of priors, simulations indicate that this method can have trouble detecting multiple WGDs that occur on the same branch, especially when the gene retention rates for ancient WGDs are very low. They also suggest that we should carefully evaluate evidence for some ancient plant WGD hypotheses., (© The Author 2016. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

Reconstructing the phylogenetic relationships that unite all lineages (the tree of life) is a grand challenge. The paucity of homologous character data across disparately related lineages currently renders direct phylogenetic inference untenable. To reconstruct a comprehensive tree of life, we therefore synthesized published phylogenies, together with taxonomic classifications for taxa never incorporated into a phylogeny. We present a draft tree containing 2.3 million tips-the Open Tree of Life. Realization of this tree required the assembly of two additional community resources: (i) a comprehensive global reference taxonomy and (ii) a database of published phylogenetic trees mapped to this taxonomy. Our open source framework facilitates community comment and contribution, enabling the tree to be continuously updated when new phylogenetic and taxonomic data become digitally available. Although data coverage and phylogenetic conflict across the Open Tree of Life illuminate gaps in both the underlying data available for phylogenetic reconstruction and the publication of trees as digital objects, the tree provides a compelling starting point for community contribution. This comprehensive tree will fuel fundamental research on the nature of biological diversity, ultimately providing up-to-date phylogenies for downstream applications in comparative biology, ecology, conservation biology, climate change, agriculture, and genomics.

Background: Although homologous recombination affects the efficacy of selection in populations, the pattern of recombination rate evolution and its effects on genome evolution across plants are largely unknown. Recombination can reduce genome size by enabling the removal of LTR retrotransposons, alter codon usage by GC biased gene conversion, contribute to complex histories of gene duplication and loss through tandem duplication, and enhance purifying selection on genes. Therefore, variation in recombination rate across species may explain some of the variation in genomic architecture as well as rates of molecular evolution. We used phylogenetic comparative methods to investigate the evolution of global meiotic recombination rate in angiosperms and its effects on genome architecture and selection at the molecular level using genetic maps and genome sequences from thirty angiosperm species., Results: Recombination rate is negatively correlated with genome size, which is likely caused by the removal of LTR retrotransposons. After correcting recombination rates for euchromatin content, we also found an association between global recombination rate and average gene family size. This suggests a role for recombination in the preservation of duplicate genes or expansion of gene families. An analysis of the correlation between the ratio of nonsynonymous to synonymous substitution rates (dN/dS) and recombination rate in 3748 genes indicates that higher recombination rates are associated with an increased efficacy of purifying selection, suggesting that global recombination rates affect variation in rates of molecular evolution across distantly related angiosperm species, not just between populations. We also identified shifts in dN/dS for recombination proteins that are associated with shifts in global recombination rate across our sample of angiosperms., Conclusions: Although our analyses only reveal correlations, not mechanisms, and do not include potential covariates of recombination rate, like effective population size, they suggest that global recombination rates may play an important role in shaping the macroevolutionary patterns of gene and genome evolution in plants. Interspecific recombination rate variation is tightly correlated with genome size as well as variation in overall LTR retrotransposon abundances. Recombination may shape gene-to-gene variation in dN/dS between species, which might impact the overall gene duplication and loss rates.

Birds are the most diverse tetrapod class, with about 10,000 extant species that represent a remarkable evolutionary radiation in which most taxa arose during a short period of time. There has been a tremendous increase in the amount of molecular data available from birds, and more than two-thirds of these species have some sequence data available. Here we assembled these available sequence data from birds to estimate a large-scale avian phylogeny. We performed an unconstrained maximum likelihood analysis of a sparse supermatrix comprising 22 nuclear loci and seven mitochondrial regions from 6714 species. We inferred a phylogeny with a backbone remarkably similar to that obtained by detailed analyses of multigene datasets, yet with the addition of thousands of more taxa. All orders were monophyletic with generally high support. While most families and genera were well supported, a number of them, especially within the oscine passerines, had little or no support. This likely reflects problems with the circumscription of these genera and families. Our results indicate that the amount of sequence data currently available is sufficient to produce a robust estimate of the avian tree of life using current methods of inference. The availability of a tree that is unconstrained by prior information, with branch lengths that have a direct connection to the underlying data, should be useful for comparative methods, taxonomic revisions, and prioritizing taxa that should be targeted for additional data collection., (Copyright © 2014 Elsevier Inc. All rights reserved.)

With the availability of genomic sequence data, there is increasing interest in using genes with a possible history of duplication and loss for species tree inference. Here we assess the performance of both nonprobabilistic and probabilistic species tree inference approaches using gene duplication and loss and coalescence simulations. We evaluated the performance of gene tree parsimony (GTP) based on duplication (Only-dup), duplication and loss (Dup-loss), and deep coalescence (Deep-c) costs, the NJst distance method, the MulRF supertree method, and PHYLDOG, which jointly estimates gene trees and species tree using a hierarchical probabilistic model. We examined the effects of gene tree and species sampling, gene tree error, and duplication and loss rates on the accuracy of phylogenetic estimates. In the 10-taxon duplication and loss simulation experiments, MulRF is more accurate than the other methods when the duplication and loss rates are low, and Dup-loss is generally the most accurate when the duplication and loss rates are high. PHYLDOG performs well in 10-taxon duplication and loss simulations, but its run time is prohibitively long on larger data sets. In the larger duplication and loss simulation experiments, MulRF outperforms all other methods in experiments with at most 100 taxa; however, in the larger simulation, Dup-loss generally performs best. In all duplication and loss simulation experiments with more than 10 taxa, all methods perform better with more gene trees and fewer missing sequences, and they are all affected by gene tree error. Our results also highlight high levels of error in estimates of duplications and losses from GTP methods and demonstrate the usefulness of methods based on generic tree distances for large analyses., (© The Author(s) 2014. Published by Oxford University Press, on behalf of the Society of Systematic Biologists. All rights reserved. For Permissions, please email: journals.permissions@oup.com.)

Unlabelled: •, Premise of the Study: While numerous software packages enable scientists to evaluate molecular data and transform them for phylogenetic analyses, few such tools exist for phenomic data. We introduce MatrixConverter, a program that helps expedite and facilitate the transformation of raw phenomic character data into discrete character matrices that can be used in most evolutionary inference programs. •, Methods and Results: MatrixConverter is an open source program written in Java; a platform-independent binary executable, as well as sample data sets and a user's manual, are available at https://github.com/gburleigh/MatrixConverter/tree/master/distribution. MatrixConverter has a simple, intuitive user interface that enables the user to immediately begin scoring phenomic characters. We demonstrate the performance of MatrixConverter on a phenomic data set from cycads. •, Conclusions: New technologies and software make it possible to obtain phenomic data from species across the tree of life, and MatrixConverter helps to transform these new data for evolutionary or ecological inference.

Analysis of large data sets can help resolve difficult nodes in the tree of life and also reveal complex evolutionary histories. The placement of the Celastrales-Oxalidales-Malpighiales (COM) clade within Rosidae remains one of the most confounding phylogenetic questions in angiosperms, with previous analyses placing it with either Fabidae or Malvidae. To elucidate the position of COM, we assembled multi-gene matrices of chloroplast, mitochondrial, and nuclear sequences, as well as large single- and multi-copy nuclear gene data sets. Analyses of multi-gene data sets demonstrate conflict between the chloroplast and both nuclear and mitochondrial data sets, and the results are robust to various character-coding and data-exclusion treatments. Analyses of single- and multi-copy nuclear loci indicate that most loci support the placement of COM with Malvidae, fewer loci support COM with Fabidae, and almost no loci support COM outside a clade of Fabidae and Malvidae. Although incomplete lineage sorting and ancient introgressive hybridization remain as plausible explanations for the conflict among loci, more complete sampling is necessary to evaluate these hypotheses fully. Our results emphasize the importance of genomic data sets for revealing deep incongruence and complex patterns of evolution., (Copyright © 2014 Elsevier Inc. All rights reserved.)

Despite a large and multifaceted effort to understand the vast landscape of phenotypic data, their current form inhibits productive data analysis. The lack of a community-wide, consensus-based, human- and machine-interpretable language for describing phenotypes and their genomic and environmental contexts is perhaps the most pressing scientific bottleneck to integration across many key fields in biology, including genomics, systems biology, development, medicine, evolution, ecology, and systematics. Here we survey the current phenomics landscape, including data resources and handling, and the progress that has been made to accurately capture relevant data descriptions for phenotypes. We present an example of the kind of integration across domains that computable phenotypes would enable, and we call upon the broader biology community, publishers, and relevant funding agencies to support efforts to surmount today's data barriers and facilitate analytical reproducibility.

Reconstructing the origin and evolution of land plants and their algal relatives is a fundamental problem in plant phylogenetics, and is essential for understanding how critical adaptations arose, including the embryo, vascular tissue, seeds, and flowers. Despite advances in molecular systematics, some hypotheses of relationships remain weakly resolved. Inferring deep phylogenies with bouts of rapid diversification can be problematic; however, genome-scale data should significantly increase the number of informative characters for analyses. Recent phylogenomic reconstructions focused on the major divergences of plants have resulted in promising but inconsistent results. One limitation is sparse taxon sampling, likely resulting from the difficulty and cost of data generation. To address this limitation, transcriptome data for 92 streptophyte taxa were generated and analyzed along with 11 published plant genome sequences. Phylogenetic reconstructions were conducted using up to 852 nuclear genes and 1,701,170 aligned sites. Sixty-nine analyses were performed to test the robustness of phylogenetic inferences to permutations of the data matrix or to phylogenetic method, including supermatrix, supertree, and coalescent-based approaches, maximum-likelihood and Bayesian methods, partitioned and unpartitioned analyses, and amino acid versus DNA alignments. Among other results, we find robust support for a sister-group relationship between land plants and one group of streptophyte green algae, the Zygnematophyceae. Strong and robust support for a clade comprising liverworts and mosses is inconsistent with a widely accepted view of early land plant evolution, and suggests that phylogenetic hypotheses used to understand the evolution of fundamental plant traits should be reevaluated.

The 1,000 plants (1KP) project is an international multi-disciplinary consortium that has generated transcriptome data from over 1,000 plant species, with exemplars for all of the major lineages across the Viridiplantae (green plants) clade. Here, we describe how to access the data used in a phylogenomics analysis of the first 85 species, and how to visualize our gene and species trees. Users can develop computational pipelines to analyse these data, in conjunction with data of their own that they can upload. Computationally estimated protein-protein interactions and biochemical pathways can be visualized at another site. Finally, we comment on our future plans and how they fit within this scalable system for the dissemination, visualization, and analysis of large multi-species data sets.

Rust fungi are a group of fungal pathogens that cause some of the world's most destructive diseases of trees and crops. A shared characteristic among rust fungi is obligate biotrophy, the inability to complete a lifecycle without a host. This dependence on a host species likely affects patterns of gene expansion, contraction, and innovation within rust pathogen genomes. The establishment of disease by biotrophic pathogens is reliant upon effector proteins that are encoded in the fungal genome and secreted from the pathogen into the host's cell apoplast or within the cells. This study uses a comparative genomic approach to elucidate putative effectors and determine their evolutionary histories. We used OrthoMCL to identify nearly 20,000 gene families in proteomes of 16 diverse fungal species, which include 15 basidiomycetes and one ascomycete. We inferred patterns of duplication and loss for each gene family and identified families with distinctive patterns of expansion/contraction associated with the evolution of rust fungal genomes. To recognize potential contributors for the unique features of rust pathogens, we identified families harboring secreted proteins that: (i) arose or expanded in rust pathogens relative to other fungi, or (ii) contracted or were lost in rust fungal genomes. While the origin of rust fungi appears to be associated with considerable gene loss, there are many gene duplications associated with each sampled rust fungal genome. We also highlight two putative effector gene families that have expanded in Cqf that we hypothesize have roles in pathogenicity.

The study of which life history traits primarily affect molecular evolutionary rates is often confounded by the covariance of these traits. Scombroid fishes (billfishes, tunas, barracudas, and their relatives) are unusual in that their mass-specific metabolic rate is positively associated with body size. This study exploits this atypical pattern of trait variation, which allows for direct tests of whether mass-specific metabolic rate or body size is the more important factor of molecular evolutionary rates. We inferred a phylogeny for scombroids from a supermatrix of molecular and morphological characters and used new phylogenetic comparative approaches to assess the associations of body size and mass-specific metabolic rate with substitution rate. As predicted by the body size hypothesis, there is a negative correlation between body size and substitution rate. However, unexpectedly, we also find a negative association between mass-specific metabolic and substitution rates. These relationships are supported by analyses of the total molecular data, separate mitochondrial and nuclear genes, and individual loci, and they are robust to phylogenetic uncertainty. The molecular evolutionary rates of scombroids are primarily tied to body size. This study demonstrates that groups with novel patterns of trait variation can be particularly informative for identifying which life history traits are the primary factors of molecular evolutionary rates.

Background: Next-generation sequencing has provided a wealth of plastid genome sequence data from an increasingly diverse set of green plants (Viridiplantae). Although these data have helped resolve the phylogeny of numerous clades (e.g., green algae, angiosperms, and gymnosperms), their utility for inferring relationships across all green plants is uncertain. Viridiplantae originated 700-1500 million years ago and may comprise as many as 500,000 species. This clade represents a major source of photosynthetic carbon and contains an immense diversity of life forms, including some of the smallest and largest eukaryotes. Here we explore the limits and challenges of inferring a comprehensive green plant phylogeny from available complete or nearly complete plastid genome sequence data., Results: We assembled protein-coding sequence data for 78 genes from 360 diverse green plant taxa with complete or nearly complete plastid genome sequences available from GenBank. Phylogenetic analyses of the plastid data recovered well-supported backbone relationships and strong support for relationships that were not observed in previous analyses of major subclades within Viridiplantae. However, there also is evidence of systematic error in some analyses. In several instances we obtained strongly supported but conflicting topologies from analyses of nucleotides versus amino acid characters, and the considerable variation in GC content among lineages and within single genomes affected the phylogenetic placement of several taxa., Conclusions: Analyses of the plastid sequence data recovered a strongly supported framework of relationships for green plants. This framework includes: i) the placement of Zygnematophyceace as sister to land plants (Embryophyta), ii) a clade of extant gymnosperms (Acrogymnospermae) with cycads + Ginkgo sister to remaining extant gymnosperms and with gnetophytes (Gnetophyta) sister to non-Pinaceae conifers (Gnecup trees), and iii) within the monilophyte clade (Monilophyta), Equisetales + Psilotales are sister to Marattiales + leptosporangiate ferns. Our analyses also highlight the challenges of using plastid genome sequences in deep-level phylogenomic analyses, and we provide suggestions for future analyses that will likely incorporate plastid genome sequence data for thousands of species. We particularly emphasize the importance of exploring the effects of different partitioning and character coding strategies.

Phylogenetic inference is a computationally difficult problem, and constructing high-quality phylogenies that can build upon existing phylogenetic knowledge and synthesize insights from new data remains a major challenge. We introduce knowledge-enhanced phylogenetic problems for both supertree and supermatrix phylogenetic analyses. These problems seek an optimal phylogenetic tree that can only be assembled from a user-supplied set of, possibly incompatible, phylogenetic relationships. We describe exact polynomial time algorithms for the knowledge-enhanced versions of the NP-hard Robinson Foulds, gene duplication, duplication and loss, and deep coalescence supertree problems. Further, we demonstrate that our algorithms can rapidly improve upon results of local search heuristics for these problems. Finally, we introduce a knowledge-enhanced search heuristic that can be applied to any discrete character data set using the maximum parsimony (MP) phylogenetic problem. Although this approach is not guaranteed to find exact solutions, we show that it also can improve upon solutions from commonly used MP heuristics.

The phenotype represents a critical interface between the genome and the environment in which organisms live and evolve. Phenotypic characters also are a rich source of biodiversity data for tree building, and they enable scientists to reconstruct the evolutionary history of organisms, including most fossil taxa, for which genetic data are unavailable. Therefore, phenotypic data are necessary for building a comprehensive Tree of Life. In contrast to recent advances in molecular sequencing, which has become faster and cheaper through recent technological advances, phenotypic data collection remains often prohibitively slow and expensive. The next-generation phenomics project is a collaborative, multidisciplinary effort to leverage advances in image analysis, crowdsourcing, and natural language processing to develop and implement novel approaches for discovering and scoring the phenome, the collection of phentotypic characters for a species. This research represents a new approach to data collection that has the potential to transform phylogenetics research and to enable rapid advances in constructing the Tree of Life. Our goal is to assemble large phenomic datasets built using new methods and to provide the public and scientific community with tools for phenomic data assembly that will enable rapid and automated study of phenotypes across the Tree of Life.

Premise of the Study: We sought novel evolutionary insights for the highly diverse Saxifragales by constructing a large phylogenetic tree encompassing 36.8% of the species-level biodiversity. •, Methods: We built a phylogenetic tree for 909 species of Saxifragales and used this hypothesis to examine character evolution for annual or perennial habit, woody or herbaceous habit, ovary position, petal number, carpel number, and stamen to petal ratio. We employed likelihood approaches to investigate the effect of habit and life history on speciation and extinction within this clade. •, Key Results: Two major shifts occurred from a woody ancestor to the herbaceous habit, with multiple secondary changes from herbaceous to woody. Transitions among superior, subinferior, and inferior ovaries appear equiprobable. A major increase in petal number is correlated with a large increase in carpel number; these increases have co-occurred multiple times in Crassulaceae. Perennial or woody lineages have higher rates of speciation than annual or herbaceous ones, but higher probabilities of extinction offset these differences. Hence, net diversification rates are highest for annual, herbaceous lineages and lowest for woody perennials. The shift from annuality to perenniality in herbaceous taxa is frequent. Conversely, woody perennial lineages to woody annual transitions are infrequent; if they occur, the woody annual state is left immediately. •, Conclusions: The large tree provides new insights into character evolution that are not obvious with smaller trees. Our results indicate that in some cases the evolution of angiosperms might be conditioned by constraints that have been so far overlooked.

The origin and maintenance of separate sexes (dioecy) is an enduring evolutionary puzzle. Although both hermaphroditism and dioecy occur in many diverse clades, we know little about the long-term evolutionary consequences of changing sexual system. Here we find evidence for at least 133 transitions between sexual systems in mosses, representing an almost unparalleled lability in the evolution of their sexual systems. Furthermore, in contrast to predictions, the transition rate from hermaphroditism to dioecy was approximately twice as high as the reverse transition. Our results also suggest that hermaphrodites may have higher rates of diversification than dioecious mosses. These results illustrate the utility of mosses for understanding the genomic and macroevolutionary consequences of hermaphroditism and dioecy., (© 2012 The Author(s). Evolution© 2012 The Society for the Study of Evolution.)

The success of model-based methods in phylogenetics has motivated much research aimed at generating new, biologically informative models. This new computer-intensive approach to phylogenetics demands validation studies and sound measures of performance. To date there has been little practical guidance available as to when and why the parameters in a particular model can be identified reliably. Here, we illustrate how Data Cloning (DC), a recently developed methodology to compute the maximum likelihood estimates along with their asymptotic variance, can be used to diagnose structural parameter nonidentifiability (NI) and distinguish it from other parameter estimability problems, including when parameters are structurally identifiable, but are not estimable in a given data set (INE), and when parameters are identifiable, and estimable, but only weakly so (WE). The application of the DC theorem uses well-known and widely used Bayesian computational techniques. With the DC approach, practitioners can use Bayesian phylogenetics software to diagnose nonidentifiability. Theoreticians and practitioners alike now have a powerful, yet simple tool to detect nonidentifiability while investigating complex modeling scenarios, where getting closed-form expressions in a probabilistic study is complicated. Furthermore, here we also show how DC can be used as a tool to examine and eliminate the influence of the priors, in particular if the process of prior elicitation is not straightforward. Finally, when applied to phylogenetic inference, DC can be used to study at least two important statistical questions: assessing identifiability of discrete parameters, like the tree topology, and developing efficient sampling methods for computationally expensive posterior densities.

A better understanding of the factors that govern individual cell lifespan and the replicative capacity of cells (i.e. Hayflick's limit) is important for addressing disease progression and ageing. Estimates of cell lifespan in vivo and the replicative capacity of cell lines in culture vary substantially both within and across species, but the underlying reasons for this variability remain unclear. Here, we address this issue by presenting a quantitative model of cell lifespan and cell replicative capacity. The model is based on the relationship between cell mortality and metabolic rate, which is supported with data for different cell types from ectotherms and endotherms. These data indicate that much of the observed variation in cell lifespan and cell replicative capacity is explained by differences in cellular metabolic rate, and thus by the three primary factors that control metabolic rate: organism size, organism temperature and cell size. Individual cell lifespan increases as a power law with both body mass and cell mass, and decreases exponentially with increasing temperature. The replicative capacity of cells also increases with body mass, but is independent of temperature. These results provide a point of departure for future comparative studies of cell lifespan and replicative capacity in the laboratory and in the field.

Background: We consider the problem of finding the maximum frequent agreement subtrees (MFASTs) in a collection of phylogenetic trees. Existing methods for this problem often do not scale beyond datasets with around 100 taxa. Our goal is to address this problem for datasets with over a thousand taxa and hundreds of trees., Results: We develop a heuristic solution that aims to find MFASTs in sets of many, large phylogenetic trees. Our method works in multiple phases. In the first phase, it identifies small candidate subtrees from the set of input trees which serve as the seeds of larger subtrees. In the second phase, it combines these small seeds to build larger candidate MFASTs. In the final phase, it performs a post-processing step that ensures that we find a frequent agreement subtree that is not contained in a larger frequent agreement subtree. We demonstrate that this heuristic can easily handle data sets with 1000 taxa, greatly extending the estimation of MFASTs beyond current methods., Conclusions: Although this heuristic does not guarantee to find all MFASTs or the largest MFAST, it found the MFAST in all of our synthetic datasets where we could verify the correctness of the result. It also performed well on large empirical data sets. Its performance is robust to the number and size of the input trees. Overall, this method provides a simple and fast way to identify strongly supported subtrees within large phylogenetic hypotheses.

The first analyses of gene sequence data indicated that the eukaryotic tree of life consisted of a long stem of microbial groups "topped" by a crown-containing plants, animals, and fungi and their microbial relatives. Although more recent multigene concatenated analyses have refined the relationships among the many branches of eukaryotes, the root of the eukaryotic tree of life has remained elusive. Inferring the root of extant eukaryotes is challenging because of the age of the group (∼1.7-2.1 billion years old), tremendous heterogeneity in rates of evolution among lineages, and lack of obvious outgroups for many genes. Here, we reconstruct a rooted phylogeny of extant eukaryotes based on minimizing the number of duplications and losses among a collection of gene trees. This approach does not require outgroup sequences or assumptions of orthology among sequences. We also explore the impact of taxon and gene sampling and assess support for alternative hypotheses for the root. Using 20 gene trees from 84 diverse eukaryotic lineages, this approach recovers robust eukaryotic clades and reveals evidence for a eukaryotic root that lies between the Opisthokonta (animals, fungi and their microbial relatives) and all remaining eukaryotes.

A Robinson-Foulds (RF) supertree for a collection of input trees is a tree containing all the species in the input trees that is at minimum total RF distance to the input trees. Thus, an RF supertree is consistent with the maximum number of splits in the input trees. Constructing RF supertrees for rooted and unrooted data is NP-hard. Nevertheless, effective local search heuristics have been developed for the restricted case where the input trees and the supertree are rooted. We describe new heuristics, based on the Edge Contract and Refine (ECR) operation, that remove this restriction, thereby expanding the utility of RF supertrees. Our experimental results on simulated and empirical data sets show that our unrooted local search algorithms yield better supertrees than those obtained from MRP and rooted RF heuristics in terms of total RF distance to the input trees and, for simulated data, in terms of RF distance to the true tree.

Background: To infer a species phylogeny from unlinked genes, phylogenetic inference methods must confront the biological processes that create incongruence between gene trees and the species phylogeny. Intra-specific gene variation in ancestral species can result in deep coalescence, also known as incomplete lineage sorting, which creates incongruence between gene trees and the species tree. One approach to account for deep coalescence in phylogenetic analyses is the deep coalescence problem, which takes a collection of gene trees and seeks the species tree that implies the fewest deep coalescence events. Although this approach is promising for phylogenetics, the consensus properties of this problem are mostly unknown and analyses of large data sets may be computationally prohibitive., Results: We prove that the deep coalescence consensus tree problem satisfies the highly desirable Pareto property for clusters (clades). That is, in all instances, each cluster that is present in all of the input gene trees, called a consensus cluster, will also be found in every optimal solution. Moreover, we introduce a new divide and conquer method for the deep coalescence problem based on the Pareto property. This method refines the strict consensus of the input gene trees, thereby, in practice, often greatly reducing the complexity of the tree search and guaranteeing that the estimated species tree will satisfy the Pareto property., Conclusions: Analyses of both simulated and empirical data sets demonstrate that the divide and conquer method can greatly improve upon the speed of heuristics that do not consider the Pareto consensus property, while also guaranteeing that the proposed solution fulfills the Pareto property. The divide and conquer method extends the utility of the deep coalescence problem to data sets with enormous numbers of taxa.

Background: Gene tree - species tree reconciliation problems infer the patterns and processes of gene evolution within a species tree. Gene tree parsimony approaches seek the evolutionary scenario that implies the fewest gene duplications, duplications and losses, or deep coalescence (incomplete lineage sorting) events needed to reconcile a gene tree and a species tree. While a gene tree parsimony approach can be informative about genome evolution and phylogenetics, error in gene trees can profoundly bias the results., Results: We introduce efficient algorithms that rapidly search local Subtree Prune and Regraft (SPR) or Tree Bisection and Reconnection (TBR) neighborhoods of a given gene tree to identify a topology that implies the fewest duplications, duplication and losses, or deep coalescence events. These algorithms improve on the current solutions by a factor of n for searching SPR neighborhoods and n2 for searching TBR neighborhoods, where n is the number of taxa in the given gene tree. They provide a fast error correction protocol for ameliorating the effects of gene tree error by allowing small rearrangements in the topology to improve the reconciliation cost. We also demonstrate a simple protocol to use the gene rearrangement algorithm to improve gene tree parsimony phylogenetic analyses., Conclusions: The new gene tree rearrangement algorithms provide a fast method to address gene tree error. They do not make assumptions about the underlying processes of genome evolution, and they are amenable to analyses of large-scale genomic data sets. These algorithms are also easily incorporated into gene tree parsimony phylogenetic analyses, potentially producing more credible estimates of reconciliation cost.

Slow-growing juveniles of shade-tolerant plant species are predicted to have tough leaves because of the high cost of leaf replacement in shade relative to potential carbon gain. We assessed the degree of correlated evolution among eight traits associated with leaf toughness and the relationships of those traits with the growth and mortality rates of 197 tree and shrub species from the understory of the 50-ha forest dynamics plot on Barro Colorado Island, Panama. Path analysis with phylogenetically independent contrasts revealed that leaves attained material toughness (resistance to fracture per unit fracture area) through increases in tissue density, percent cellulose per unit dry mass, and vein fracture toughness. Lamina density and cellulose content evolved independently and thus represent different paths to material toughness. Structural toughness (resistance to fracture per unit fracture length) depended on material toughness and lamina thickness. Mortality rates of individuals 1-10 cm in stem diameter were negatively correlated with material toughness and lamina density but were independent of structural toughness and cell wall fiber contents. Leaf toughness traits were uncorrelated with relative growth rates. Results imply that material toughness enhances resistance to natural enemies, which increases survival and offsets the biomass allocation cost of producing tough leaves in the shaded understory.

Phylogenetic analyses using genome-scale data sets must confront incongruence among gene trees, which in plants is exacerbated by frequent gene duplications and losses. Gene tree parsimony (GTP) is a phylogenetic optimization criterion in which a species tree that minimizes the number of gene duplications induced among a set of gene trees is selected. The run time performance of previous implementations has limited its use on large-scale data sets. We used new software that incorporates recent algorithmic advances to examine the performance of GTP on a plant data set consisting of 18,896 gene trees containing 510,922 protein sequences from 136 plant taxa (giving a combined alignment length of >2.9 million characters). The relationships inferred from the GTP analysis were largely consistent with previous large-scale studies of backbone plant phylogeny and resolved some controversial nodes. The placement of taxa that were present in few gene trees generally varied the most among GTP bootstrap replicates. Excluding these taxa either before or after the GTP analysis revealed high levels of phylogenetic support across plants. The analyses supported magnoliids sister to a eudicot + monocot clade and did not support the eurosid I and II clades. This study presents a nuclear genomic perspective on the broad-scale phylogenic relationships among plants, and it demonstrates that nuclear genes with a history of duplication and loss can be phylogenetically informative for resolving the plant tree of life.

Although Pentapetalae (comprising all core eudicots except Gunnerales) include approximately 70% of all angiosperms, the origin of and relationships among the major lineages of this clade have remained largely unresolved. Phylogenetic analyses of 83 protein-coding and rRNA genes from the plastid genome for 86 species of seed plants, including new sequences from 25 eudicots, indicate that soon after its origin, Pentapetalae diverged into three clades: (i) a "superrosid" clade consisting of Rosidae, Vitaceae, and Saxifragales; (ii) a "superasterid" clade consisting of Berberidopsidales, Santalales, Caryophyllales, and Asteridae; and (iii) Dilleniaceae. Maximum-likelihood analyses support the position of Dilleniaceae as sister to superrosids, but topology tests did not reject alternative positions of Dilleniaceae as sister to Asteridae or all remaining Pentapetalae. Molecular dating analyses suggest that the major lineages within both superrosids and superasterids arose in as little as 5 million years. This phylogenetic hypothesis provides a crucial historical framework for future studies aimed at elucidating the underlying causes of the morphological and species diversity in Pentapetalae.

Background: Supertree methods synthesize collections of small phylogenetic trees with incomplete taxon overlap into comprehensive trees, or supertrees, that include all taxa found in the input trees. Supertree methods based on the well established Robinson-Foulds (RF) distance have the potential to build supertrees that retain much information from the input trees. Specifically, the RF supertree problem seeks a binary supertree that minimizes the sum of the RF distances from the supertree to the input trees. Thus, an RF supertree is a supertree that is consistent with the largest number of clusters (or clades) from the input trees., Results: We introduce efficient, local search based, hill-climbing heuristics for the intrinsically hard RF supertree problem on rooted trees. These heuristics use novel non-trivial algorithms for the SPR and TBR local search problems which improve on the time complexity of the best known (naïve) solutions by a factor of Theta(n) and Theta(n2) respectively (where n is the number of taxa, or leaves, in the supertree). We use an implementation of our new algorithms to examine the performance of the RF supertree method and compare it to matrix representation with parsimony (MRP) and the triplet supertree method using four supertree data sets. Not only did our RF heuristic provide fast estimates of RF supertrees in all data sets, but the RF supertrees also retained more of the information from the input trees (based on the RF distance) than the other supertree methods., Conclusions: Our heuristics for the RF supertree problem, based on our new local search algorithms, make it possible for the first time to estimate large supertrees by directly optimizing the RF distance from rooted input trees to the supertrees. This provides a new and fast method to build accurate supertrees. RF supertrees may also be useful for estimating majority-rule(-) supertrees, which are a generalization of majority-rule consensus trees.