Taller plants have lower rates of molecular evolution

Drivers of avian genomic change revealed by evolutionary rate decomposition

Modern birds have diversified into a striking array of forms, behaviours and ecological roles. Analyses of molecular evolutionary rates can reveal the links between genomic and phenotypic change1-4, but disentangling the drivers of rate variation at the whole-genome scale has been difficult. Using comprehensive estimates of traits and evolutionary rates across a family-level phylogeny of birds5,6, we find that genome-wide mutation rates across lineages are predominantly explained by clutch size and generation length, whereas rate variation across genes is driven by the content of guanine and cytosine. Here, to find the subsets of genes and lineages that dominate evolutionary rate variation in birds, we estimated the influence of individual lineages on decomposed axes of gene-specific evolutionary rates. We find that most of the rate variation occurs along recent branches of the tree, associated with present-day families of birds. Additional tests on axes of rate variation show rapid changes in microchromosomes immediately after the Cretaceous-Palaeogene transition. These apparent pulses of evolution are consistent with major changes in the genetic machineries for meiosis, heart performance, and RNA splicing, surveillance and translation, and correlate with the ecological diversity reflected in increased tarsus length. Collectively, our analyses paint a nuanced picture of avian evolution, revealing that the ancestors of the most diverse lineages of birds underwent major genomic changes related to mutation, gene usage and niche expansion in the early Palaeogene period.

Allometric scaling of somatic mutation and epimutation rates in trees

How long-lived trees escape "mutational meltdown" despite centuries of continuous growth remains puzzling. Here we integrate recent studies to show that the yearly rate of somatic mutations and epimutations (μY) scales inversely with generation time (G), and follows the same allometric power law found in mammals (μY ∝ G-1). Deeper insights into the scaling function may permit predictions of somatic (epi)mutation rates from life-history traits without the need for genomic data.

Combining Molecular, Macroevolutionary, and Macroecological Perspectives on the Generation of Diversity

Charles Darwin presented a unified process of diversification driven by the gradual accumulation of heritable variation. The growth in DNA databases and the increase in genomic sequencing, combined with advances in molecular phylogenetic analyses, gives us an opportunity to realize Darwin's vision, connecting the generation of variation to the diversification of lineages. The rate of molecular evolution is correlated with the rate of diversification across animals and plants, but the relationship between genome change and speciation is complex: Mutation rates evolve in response to life history and niche; substitution rates are influenced by mutation, selection, and population size; rates of acquisition of reproductive isolation vary between populations; and traits, niches, and distribution can influence diversification rates. The connection between mutation rate and diversification rate is one part of the complex and varied story of speciation, which has theoretical importance for understanding the generation of biodiversity and also practical impacts on the use of DNA to understand the dynamics of speciation over macroevolutionary timescales.

Fast Bayesian Inference of Phylogenies from Multiple Continuous Characters

Time-scaled phylogenetic trees are an ultimate goal of evolutionary biology and a necessary ingredient in comparative studies. The accumulation of genomic data has resolved the tree of life to a great extent, yet timing evolutionary events remain challenging if not impossible without external information such as fossil ages and morphological characters. Methods for incorporating morphology in tree estimation have lagged behind their molecular counterparts, especially in the case of continuous characters. Despite recent advances, such tools are still direly needed as we approach the limits of what molecules can teach us. Here, we implement a suite of state-of-the-art methods for leveraging continuous morphology in phylogenetics, and by conducting extensive simulation studies we thoroughly validate and explore our methods' properties. While retaining model generality and scalability, we make it possible to estimate absolute and relative divergence times from multiple continuous characters while accounting for uncertainty. We compile and analyze one of the most data-type diverse data sets to date, comprised of contemporaneous and ancient molecular sequences, and discrete and continuous morphological characters from living and extinct Carnivora taxa. We conclude by synthesizing lessons about our method's behavior, and suggest future research venues.

The chromosome-level genome of the submerged plant Cryptocoryne crispatula provides insights into the terrestrial-freshwater transition in Araceae

Plant terrestrialization (i.e. the transition to a terrestrial environment) is a significant evolutionary event that has been intensively studied. While certain plant lineages, particularly in angiosperms, have re-adapted to freshwater habitats after colonizing terrene, however, the molecular mechanism of the terrestrial-freshwater (T-F) transition remains limited. Here, the basal monocot Araceae was selected as the study object to explore the T-F transition adaptation mechanism by comparative genomic analysis. Our findings revealed that the substitution rates significantly increased in the lineage of freshwater Araceae, which may promote their adaptation to the freshwater habitat. Additionally, 20 gene sets across all four freshwater species displayed signs of positive selection contributing to tissue development and defense responses in freshwater plants. Comparative synteny analysis showed that genes specific to submerged plants were enriched in cellular respiration and photosynthesis. In contrast, floating plants were involved in regulating gene expression, suggesting that gene and genome duplications may provide the original material for plants to adapt to the freshwater environment. Our study provides valuable insights into the genomic aspects of the transition from terrestrial to aquatic environments in Araceae, laying the groundwork for future research in the angiosperm.

A reference genome of Commelinales provides insights into the commelinids evolution and global spread of water hyacinth (Pontederia crassipes)

Commelinales belongs to the commelinids clade, which also comprises Poales that includes the most important monocot species, such as rice, wheat, and maize. No reference genome of Commelinales is currently available. Water hyacinth (Pontederia crassipes or Eichhornia crassipes), a member of Commelinales, is one of the devastating aquatic weeds, although it is also grown as an ornamental and medical plant. Here, we present a chromosome-scale reference genome of the tetraploid water hyacinth with a total length of 1.22 Gb (over 95% of the estimated size) across 8 pseudochromosome pairs. With the representative genomes, we reconstructed a phylogeny of the commelinids, which supported Zingiberales and Commelinales being sister lineages of Arecales and shed lights on the controversial relationship of the orders. We also reconstructed ancestral karyotypes of the commelinids clade and confirmed the ancient commelinids genome having 8 chromosomes but not 5 as previously reported. Gene family analysis revealed contraction of disease-resistance genes during polyploidization of water hyacinth, likely a result of fitness requirement for its role as a weed. Genetic diversity analysis using 9 water hyacinth lines from 3 continents (South America, Asia, and Europe) revealed very closely related nuclear genomes and almost identical chloroplast genomes of the materials, as well as provided clues about the global dispersal of water hyacinth. The genomic resources of P. crassipes reported here contribute a crucial missing link of the commelinids species and offer novel insights into their phylogeny.

Compositional shifts associated with major evolutionary transitions in plants

Heterogeneity in gene trees, morphological characters, and composition has been associated with several major plant clades. Here, we examine heterogeneity in composition across a large transcriptomic dataset of plants to better understand whether locations of shifts in composition are shared across gene regions and whether directions of shifts within clades are shared across gene regions. We estimate mixed models of composition for both nucleotide and amino acids across a recent large-scale transcriptomic dataset for plants. We find shifts in composition across both nucleotide and amino acid datasets, with more shifts detected in nucleotides. We find that Chlorophytes and lineages within experience the most shifts. However, many shifts occur at the origins of land, vascular, and seed plants. While genes in these clades do not typically share the same composition, they tend to shift in the same direction. We discuss potential causes of these patterns. Compositional heterogeneity has been highlighted as a potential problem for phylogenetic analysis, but the variation presented here highlights the need to further investigate these patterns for the signal of biological processes.

CALANGO: A phylogeny-aware comparative genomics tool for discovering quantitative genotype-phenotype associations across species

Living species vary significantly in phenotype and genomic content. Sophisticated statistical methods linking genes with phenotypes within a species have led to breakthroughs in complex genetic diseases and genetic breeding. Despite the abundance of genomic and phenotypic data available for thousands of species, finding genotype-phenotype associations across species is challenging due to the non-independence of species data resulting from common ancestry. To address this, we present CALANGO (comparative analysis with annotation-based genomic components), a phylogeny-aware comparative genomics tool to find homologous regions and biological roles associated with quantitative phenotypes across species. In two case studies, CALANGO identified both known and previously unidentified genotype-phenotype associations. The first study revealed unknown aspects of the ecological interaction between Escherichia coli, its integrated bacteriophages, and the pathogenicity phenotype. The second identified an association between maximum height in angiosperms and the expansion of a reproductive mechanism that prevents inbreeding and increases genetic diversity, with implications for conservation biology and agriculture.

Weaker selection on genes with treatment-specific expression consistent with a limit on plasticity evolution in Arabidopsis thaliana

Differential gene expression between environments often underlies phenotypic plasticity. However, environment-specific expression patterns are hypothesized to relax selection on genes, and thus limit plasticity evolution. We collated over 27 terabases of RNA-sequencing data on Arabidopsis thaliana from over 300 peer-reviewed studies and 200 treatment conditions to investigate this hypothesis. Consistent with relaxed selection, genes with more treatment-specific expression have higher levels of nucleotide diversity and divergence at nonsynonymous sites but lack stronger signals of positive selection. This result persisted even after controlling for expression level, gene length, GC content, the tissue specificity of expression, and technical variation between studies. Overall, our investigation supports the existence of a hypothesized trade-off between the environment specificity of a gene's expression and the strength of selection on said gene in A. thaliana. Future studies should leverage multiple genome-scale datasets to tease apart the contributions of many variables in limiting plasticity evolution.

Rate accelerations in plastid and mitochondrial genomes of Cyperaceae occur in the same clades

Cyperaceae, the second largest family in the monocot order Poales, comprises >5500 species and includes the genus Eleocharis with ∼ 250 species. A previous study of complete plastomes of two Eleocharis species documented extensive structural heteroplasmy, gene order changes, high frequency of dispersed repeats along with gene losses and duplications. To better understand the phylogenetic distribution of gene and intron content as well as rates and patterns of sequence evolution within and between mitochondrial and plastid genomes of Eleocharis and Cyperaceae, an additional 29 Eleocharis organelle genomes were sequenced and analyzed. Eleocharis experienced extensive gene loss in both genomes while loss of introns was mitochondria-specific. Eleocharis has higher rates of synonymous (dS) and nonsynonymous (dN) substitutions in the plastid and mitochondrion than most sampled angiosperms, and the pattern was distinct from other eudicot lineages with accelerated rates. Several clades showed higher dS and dN in mitochondrial genes than in plastid genes. Furthermore, nucleotide substitution rates of mitochondrial genes were significantly accelerated on the branch leading to Cyperaceae compared to most angiosperms. Mitochondrial genes of Cyperaceae exhibited dramatic loss of RNA editing sites and a negative correlation between RNA editing and dS values was detected among angiosperms. Mutagenic retroprocessing and dysfunction of DNA replication, repair and recombination genes are the most likely cause of striking rate accelerations and loss of edit sites and introns in Eleocharis and Cyperaceae organelle genomes.

Multifeature analysis of age-related microbiome structures reveals defense mechanisms of Populus tomentosa trees

Root microbiota composition shifts during the development of most annual plants. Although some perennial plants can live for centuries, the host-microbiome partnerships and interaction mechanisms underlying their longevity remain unclear. To address this gap, we investigated age-related changes in the root metabolites, transcriptomes, and microbiome compositions of 1- to 35-yr-old Populus tomentosa trees. Ten co-response clusters were obtained according to their accumulation patterns, and members of each cluster displayed a uniform and clear pattern of abundance. Multi-omics network analysis demonstrated that the increased abundance of Actinobacteria with tree age was strongly associated with the flavonoid biosynthesis. Using genetic approaches, we demonstrate that the flavonoid biosynthesis regulator gene Transparent Testa 8 is associated with the recruitment of flavonoid-associated Actinobacteria. Further inoculation experiments of Actinobacteria isolates indicated that their colonization could significantly improve the host's phenotype. Site-directed mutagenesis revealed that the hyBl gene cluster, involved in biosynthesis of an aminocyclitol hygromycin B analog in Streptomyces isolate bj1, is associated with disease suppression. We hypothesize that interactions between perennial plants and soil microorganisms lead to gradual enrichment of a subset of microorganisms that may harbor a wealth of currently unknown functional traits.

Opposing effects of plant traits on diversification

Species diversity can vary dramatically across lineages due to differences in speciation and extinction rates. Here, we explore the effects of several plant traits on diversification, finding that most traits have opposing effects on diversification. For example, outcrossing may increase the efficacy of selection and adaptation but also decrease mate availability, two processes with contrasting effects on lineage persistence. Such opposing trait effects can manifest as differences in diversification rates that depend on ecological context, spatiotemporal scale, and associations with other traits. The complexity of pathways linking traits to diversification suggests that the mechanistic underpinnings behind their correlations may be difficult to interpret with any certainty, and context dependence means that the effects of specific traits on diversification are likely to differ across multiple lineages and timescales. This calls for taxonomically and context-controlled approaches to studies that correlate traits and diversification.

Genome structure-based Juglandaceae phylogenies contradict alignment-based phylogenies and substitution rates vary with DNA repair genes

In lineages of allopolyploid origin, sets of homoeologous chromosomes may coexist that differ in gene content and syntenic structure. Presence or absence of genes and microsynteny along chromosomal blocks can serve to differentiate subgenomes and to infer phylogenies. We here apply genome-structural data to infer relationships in an ancient allopolyploid lineage, the walnut family (Juglandaceae), by using seven chromosome-level genomes, two of them newly assembled. Microsynteny and gene-content analyses yield identical topologies that place Platycarya with Engelhardia as did a 1980s morphological-cladistic study. DNA-alignment-based topologies here and in numerous earlier studies instead group Platycarya with Carya and Juglans, perhaps misled by past hybridization. All available data support a hybrid origin of Juglandaceae from extinct or unsampled progenitors nested within, or sister to, Myricaceae. Rhoiptelea chiliantha, sister to all other Juglandaceae, contains proportionally more DNA repair genes and appears to evolve at a rate 2.6- to 3.5-times slower than the remaining species.

Origin and evolutionary history of Populus (Salicaceae): Further insights based on time divergence and biogeographic analysis

INTRODUCTION

Populus (Salicaceae) species harbour rich biodiversity and are widely distributed throughout the Northern Hemisphere. However, the origin and biogeography of Populus remain poorly understood.

METHODS

We infer the divergence times and the historical biogeography of the genus Populus through phylogenetic analysis of 34 chloroplast fragments based on a large sample.

RESULTS AND DISCUSSION

Eurasia is the likely location of the early divergences of Salicaceae after the Cretaceous-Paleogene (K-Pg) mass extinction, followed by recurrent spread to the remainder of the Old World and the New World beginning in the Eocene; the extant Populus species began to diversity during the early Oligocene (approximately 27.24 Ma), climate changes during the Oligocene may have facilitated the diversification of modern poplar species; three separate lineages of Populus from Eurasia colonized North America in the Cenozoic via the Bering Land Bridges (BLB); We hypothesize that the present day disjunction in Populus can be explained by two scenarios: (i) Populus likely originated in Eurasia and subsequently colonized other regions, including North America; and (ii) the fact that the ancestor of the genus Populus that was once widely distributed in the Northern Hemisphere and eventually wiped out due to the higher extinction rates in North America, similar to the African Rand flora. We hypothesize that disparities in extinction across the evolutionary history of Populus in different regions shape the modern biogeography of Populus. Further studies with dense sampling and more evidence are required to test these hypotheses. Our research underscores the significance of combining phylogenetic analyses with biogeographic interpretations to enhance our knowledge of the origin, divergence, and distribution of biodiversity in temperate plant floras.

A common whole-genome paleotetraploidization in Cucurbitales

Cucurbitales are an important order of flowering plants known for encompassing edible plants of economic and medicinal value and numerous ornamental plants of horticultural value. By reanalyzing the genomes of two representative families (Cucurbitaceae and Begoniaceae) in Cucurbitales, we found that the previously identified Cucurbitaceae common paleotetraploidization that occurred shortly after the core-eudicot-common hexaploidization event is shared by Cucurbitales, including Begoniaceae. We built a multigenome alignment framework for Cucurbitales by identifying orthologs and paralogs and systematically redating key evolutionary events in Cucurbitales. Notably, characterizing the gene retention levels and genomic fractionation patterns between subgenomes generated from different polyploidizations in Cucurbitales suggested the autopolyploid nature of the Begoniaceae common tetraploidization and the allopolyploid nature of the Cucurbitales common tetraploidization and the Cucurbita-specific tetraploidization. Moreover, we constructed the ancestral Cucurbitales karyotype comprising 17 proto-chromosomes, confirming that the most recent common ancestor of Cucurbitaceae contained 15 proto-chromosomes and rejecting the previous hypothesis for an ancestral Cucurbitaceae karyotype with 12 proto-chromosomes. In addition, we found that the polyploidization and tandem duplication events promoted the expansion of gene families involved in the cucurbitacin biosynthesis pathway; however, gene loss and chromosomal rearrangements likely limited the expansion of these gene families.

A Multigene Phylogeny of Native American Hawkweeds (Hieracium Subgen. Chionoracium, Cichorieae, Asteraceae): Origin, Speciation Patterns, and Migration Routes

Native American hawkweeds are mainly mountainous species that are distributed all over the New World. They are severely understudied with respect to their origin, colonization of the vast distribution area, and species relationships. Here, we attempt to reconstruct the evolutionary history of the group by applying seven molecular markers (plastid, nuclear ribosomal and low-copy genes). Phylogenetic analyses revealed that Chionoracium is a subgenus of the mainly Eurasian genus Hieracium, which originated from eastern European hawkweeds about 1.58-2.24 million years ago. Plastid DNA suggested a single origin of all Chionoracium species. They colonized the New World via Beringia and formed several distinct lineages in North America. Via one Central American lineage, the group colonized South America and radiated into more than a hundred species within about 0.8 million years, long after the closure of the Isthmus of Panama and the most recent uplift of the Andes. Despite some incongruences shown by different markers, most of them revealed the same crown groups of closely related taxa, which were, however, largely in conflict with traditional sectional classifications. We provide a basic framework for further elucidation of speciation patterns. A thorough taxonomic revision of Hieracium subgen. Chionoracium is recommended.

Environmental niche and flight intensity are associated with molecular evolutionary rates in a large avian radiation

BACKGROUND

Metabolic activity and environmental energy are two of the most studied putative drivers of molecular evolutionary rates. Their extensive study, however, has resulted in mixed results and has rarely included the exploration of interactions among various factors impacting molecular evolutionary rates across large clades. Taking the diverse avian family Furnariidae as a case study, we examined the association between several estimates of molecular evolutionary rates with proxies of metabolic demands imposed by flight (wing loading and wing shape) and proxies of environmental energy across the geographic ranges of species (temperature and UV radiation).

RESULTS

We found weak evidence of a positive effect of environmental and morphological variables on mitochondrial substitution rates. Additionally, we found that temperature and UV radiation interact to explain molecular rates at nucleotide sites affected by selection and population size (non-synonymous substitutions), contrary to the expectation of their impact on sites associated with mutation rates (synonymous substitutions). We also found a negative interaction between wing shape (as described by the hand-wing index) and body mass explaining mitochondrial molecular rates, suggesting molecular signatures of positive selection or reduced population sizes in small-bodied species with greater flight activity.

CONCLUSIONS

Our results suggest that the demands of flight and environmental energy pose multiple evolutionary pressures on the genome either by driving mutation rates or via their association with natural selection or population size. Data from whole genomes and detailed physiology across taxa will bring a more complete picture of the impact of metabolism, population size, and the environment on avian genome evolution.

Temperature-dependent evolutionary speed shapes the evolution of biodiversity patterns across tetrapod radiations

Biodiversity varies predictably with environmental energy around the globe, but the underlaying mechanisms remain incompletely understood. The evolutionary speed hypothesis predicts that environmental kinetic energy shapes variation in speciation rates through temperature- or life history-dependent rates of evolution. To test whether variation in evolutionary speed can explain the relationship between energy and biodiversity in birds, mammals, amphibians, and reptiles, we simulated diversification over 65 million years of geological and climatic change with a spatially explicit eco-evolutionary simulation model. We modelled four distinct evolutionary scenarios in which speciation-completion rates were dependent on temperature (M1), life history (M2), temperature and life history (M3), or were independent of temperature and life-history (M0). To assess the agreement between simulated and empirical data, we performed model selection by fitting supervised machine learning models to multidimensional biodiversity patterns. We show that a model with temperature-dependent rates of speciation (M1) consistently had the strongest support. In contrast to statistical inferences, which showed no general relationships between temperature and speciation rates in tetrapods, we demonstrate how process-based modelling can disentangle the causes behind empirical biodiversity patterns. Our study highlights how environmental energy has played a fundamental role in the evolution of biodiversity over deep time.

Phylogenomic approaches untangle early divergences and complex diversifications of the olive plant family

BACKGROUND

Deep-branching phylogenetic relationships are often difficult to resolve because phylogenetic signals are obscured by the long history and complexity of evolutionary processes, such as ancient introgression/hybridization, polyploidization, and incomplete lineage sorting (ILS). Phylogenomics has been effective in providing information for resolving both deep- and shallow-scale relationships across all branches of the tree of life. The olive family (Oleaceae) is composed of 25 genera classified into five tribes with tribe Oleeae consisting of four subtribes. Previous phylogenetic analyses showed that ILS and/or hybridization led to phylogenetic incongruence in the family. It was essential to distinguish phylogenetic signal conflicts, and explore mechanisms for the uncertainties concerning relationships of the olive family, especially at the deep-branching nodes.

RESULTS

We used the whole plastid genome and nuclear single nucleotide polymorphism (SNP) data to infer the phylogenetic relationships and to assess the variation and rates among the main clades of the olive family. We also used 2608 and 1865 orthologous nuclear genes to infer the deep-branching relationships among tribes of Oleaceae and subtribes of tribe Oleeae, respectively. Concatenated and coalescence trees based on the plastid genome, nuclear SNPs and multiple nuclear genes suggest events of ILS and/or ancient introgression during the diversification of Oleaceae. Additionally, there was extreme heterogeneity in the substitution rates across the tribes. Furthermore, our results supported that introgression/hybridization, rather than ILS, is the main factor for phylogenetic discordance among the five tribes of Oleaceae. The tribe Oleeae is supported to have originated via ancient hybridization and polyploidy, and its most likely parentages are the ancestral lineage of Jasmineae or its sister group, which is a "ghost lineage," and Forsythieae. However, ILS and ancient introgression are mainly responsible for the phylogenetic discordance among the four subtribes of tribe Oleeae.

CONCLUSIONS

This study showcases that using multiple sequence datasets (plastid genomes, nuclear SNPs and thousands of nuclear genes) and diverse phylogenomic methods such as data partition, heterogeneous models, quantifying introgression via branch lengths (QuIBL) analysis, and species network analysis can facilitate untangling long and complex evolutionary processes of ancient introgression, paleopolyploidization, and ILS.

ksrates: positioning whole-genome duplications relative to speciation events in KS distributions

SUMMARY

We present ksrates, a user-friendly command-line tool to position ancient whole-genome duplication events with respect to speciation events in a phylogeny by comparing paralog and ortholog KS distributions derived from genomic or transcriptomic sequences, while adjusting for substitution rate differences among the lineages involved.

AVAILABILITY AND IMPLEMENTATION

ksrates is implemented in Python 3 and as a Nextflow pipeline. The source code, Singularity and Docker containers, documentation and tutorial are available via https://github.com/VIB-PSB/ksrates.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Should we hail the Red King? Evolutionary consequences of a mutualistic lifestyle in genomes of lichenized ascomycetes

The Red Queen dynamic is often brought into play for antagonistic relationships. However, the coevolutionary effects of mutualistic interactions, which predict slower evolution for interacting organisms (Red King), have been investigated to a lesser extent. Lichens are a stable, mutualistic relationship of fungi and cyanobacteria and/or algae, which originated several times independently during the evolution of fungi. Therefore, they represent a suitable system to investigate the coevolutionary effect of mutualism on the fungal genome. We measured substitution rates and selective pressure of about 2000 protein-coding genes (plus the rDNA region) in two different classes of Ascomycota, each consisting of closely related lineages of lichenized and non-lichenized fungi. Our results show that independent lichenized clades are characterized by significantly slower rates for both synonymous and non-synonymous substitutions. We hypothesize that this evolutionary pattern is connected to the lichen life cycle (longer generation time of lichenized fungi) rather than a result of different selection strengths, which is described as the main driver for the Red Kind dynamic. This first empirical evidence of slower evolution in lichens provides an important insight on how biotic cooperative interactions are able to shape the evolution of symbiotic organisms.

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein-coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein-coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non-conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non-conifers. It is more obvious that in non-conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6- to 3.1-fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2- to 3.6-fold of non-conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein-coding genes due to gene function, plant habitat, or newly acquired IRs.

Distinct Responses of Leaf Traits to Environment and Phylogeny Between Herbaceous and Woody Angiosperm Species in China

Leaf traits play key roles in plant resource acquisition and ecosystem processes; however, whether the effects of environment and phylogeny on leaf traits differ between herbaceous and woody species remains unclear. To address this, in this study, we collected data for five key leaf traits from 1,819 angiosperm species across 530 sites in China. The leaf traits included specific leaf area, leaf dry matter content, leaf area, leaf N concentration, and leaf P concentration, all of which are closely related to trade-offs between resource uptake and leaf construction. We quantified the relative contributions of environment variables and phylogeny to leaf trait variation for all species, as well as for herbaceous and woody species separately. We found that environmental factors explained most of the variation (44.4-65.5%) in leaf traits (compared with 3.9-23.3% for phylogeny). Climate variability and seasonality variables, in particular, mean temperature of the warmest and coldest seasons of a year (MTWM/MTWQ and MTCM/MTCQ) and mean precipitation in the wettest and driest seasons of a year (MPWM/MPWQ and MPDM/MPDQ), were more important drivers of leaf trait variation than mean annual temperature (MAT) and mean annual precipitation (MAP). Furthermore, the responses of leaf traits to environment variables and phylogeny differed between herbaceous and woody species. Our study demonstrated the different effects of environment variables and phylogeny on leaf traits among different plant growth forms, which is expected to advance the understanding of plant adaptive strategies and trait evolution under different environmental conditions.

Population-level consequences of inheritable somatic mutations and the evolution of mutation rates in plants

Inbreeding depression, that is the decrease in fitness of inbred relative to outbred individuals, was shown to increase strongly as life expectancy increases in plants. Because plants are thought to not have a separated germline, it was proposed that this pattern could be generated by somatic mutations accumulating during growth, since larger and more long-lived species have more opportunities for mutations to accumulate. A key determinant of the role of somatic mutations is the rate at which they occur, which probably differs between species because mutation rates may evolve differently in species with constrasting life histories. In this paper, I study the evolution of the mutation rates in plants, and consider the population-level consequences of inheritable somatic mutations given this evolution. I show that despite substantially lower somatic and meiotic mutation rates, more long-lived species still tend to accumulate larger amounts of deleterious mutations because of the increased number of opportunities they have to acquire mutations during growth, leading to higher levels of inbreeding depression in these species. However, the magnitude of this increase depends strongly on how mutagenic meiosis is relative to growth, to the point of being close to non-existent in some situations.

Rates and patterns of molecular evolution in bryophyte genomes, with focus on complex thalloid liverworts, Marchantiopsida

Plants commonly referred to as "bryophytes" belong to three major lineages of non-vascular plants: the liverworts, the hornworts and the mosses. They are unique among land plants in having a dominant haploid generation and a short-lived diploid sporophytic generation. The dynamics of selection acting on a haploid genome differs from those acting on a diploid genome: new mutations are directly exposed to selection. The general aim of this paper is to investigate the diversification rateof bryophytes - measured as silent site substitution rate representing neutral evolution (mutation rate) and the nonsynonymous to synonymous substitution rate ratio (dN/dS) representing selective evolution - and compare it with earlier studies on vascular plants. Results show that the silent site substitution rate is lower for liverworts as compared to angiosperms, but not as low as for gymnosperms. The selection pressure, measured as dN/dS, isnot remarkably lower for bryophytes as compared to other diploid dominant plants as would be expected by the masking hypothesis, indicating that other factors are more important than ploidy.

The challenge and promise of estimating the de novo mutation rate from whole-genome comparisons among closely related individuals

Germline mutations are the raw material for natural selection, driving species evolution and the generation of earth's biodiversity. Without this driver of genetic diversity, life on earth would stagnate. Yet, it is a double-edged sword. An excess of mutations can have devastating effects on fitness and population viability. It is therefore one of the great challenges of molecular ecology to determine the rate and mechanisms by which these mutations accrue across the tree of life. Advances in high-throughput sequencing technologies are providing new opportunities for characterizing the rates and mutational spectra within species and populations thus informing essential evolutionary parameters such as the timing of speciation events, the intricacies of historical demography, and the degree to which lineages are subject to the burdens of mutational load. Here, we will focus on both the challenge and promise of whole-genome comparisons among parents and their offspring from known pedigrees for the detection of germline mutations as they arise in a single generation. The potential of these studies is high, but the field is still in its infancy and much uncertainty remains. Namely, the technical challenges are daunting given that pedigree-based genome comparisons are essentially searching for needles in a haystack given the very low signal to noise ratio. Despite the challenges, we predict that rapidly developing methods for whole-genome comparisons hold great promise for integrating empirically derived estimates of de novo mutation rates and mutation spectra across many molecular ecological applications.

Rapid Accumulation of Mutations in Growing Mycelia of a Hypervariable Fungus Schizophyllum commune

The basidiomycete Schizophyllum commune has the highest level of genetic polymorphism known among living organisms. In a previous study, it was also found to have a moderately high per-generation mutation rate of 2×10⁻⁸, likely contributing to its high polymorphism. However, this rate has been measured only in an experiment on Petri dishes, and it is unclear how it translates to natural populations. Here, we used an experimental design that measures the rate of accumulation of de novo mutations in a linearly growing mycelium. We show that S. commune accumulates mutations at a rate of 1.24×10⁻⁷ substitutions per nucleotide per meter of growth, or ∼2.04×10⁻¹¹ per nucleotide per cell division. In contrast to what has been observed in a number of species with extensive vegetative growth, this rate does not decline in the course of propagation of a mycelium. As a result, even a moderate per-cell-division mutation rate in S. commune can translate into a very high per-generation mutation rate when the number of cell divisions between consecutive meiosis is large.

Effects of Climate, Plant Height, and Evolutionary Age on Geographical Patterns of Fruit Type

Fruit type is a key reproductive trait associated with plant evolution and adaptation. However, large-scale geographical patterns in fruit type composition and the mechanisms driving these patterns remain to be established. Contemporary environment, plant functional traits and evolutionary age may all influence fruit type composition, while their relative importance remains unclear. Here, using data on fruit types, plant height and distributions of 28,222 (∼ 90.1%) angiosperm species in China, we analyzed the geographical patterns in the proportion of fleshy-fruited species for all angiosperms, trees, shrubs, and herbaceous species separately, and compared the relative effects of contemporary climate, ecosystem primary productivity, plant height, and evolutionary age on these patterns. We found that the proportion of fleshy-fruited species per grid cell for all species and different growth forms all showed significant latitudinal patterns, being the highest in southeastern China. Mean plant height per grid cell and actual evapotranspiration (AET) representing ecosystem primary productivity were the strongest drivers of geographical variations in the proportion of fleshy-fruited species, but their relative importance varied between growth forms. From herbaceous species to shrubs and trees, the relative effects of mean plant height decreased. Mean genus age had significant yet consistently weaker effects on proportion of fleshy-fruited species than mean plant height and AET, and environmental temperature and precipitation contributed to those of only trees and shrubs. These results suggest that biotic and environmental factors and evolutionary age of floras jointly shape the pattern in proportion of fleshy-fruited species, and improve our understanding of the mechanisms underlying geographical variations in fruit type composition. Our study also demonstrates the need of integrating multiple biotic and abiotic factors to fully understand the drivers of large-scale patterns of plant reproductive traits.

Correlation Analysis Reveals an Important Role of GC Content in Accumulation of Deletion Mutations in the Coding Region of Angiosperm Plastomes

Variation in GC content is assumed to correlate with various processes, including mutation biases, recombination, and environmental parameters. To date, most genomic studies exploring the evolution of GC content have focused on nuclear genomes, but relatively few have concentrated on organelle genomes. We explored the mechanisms maintaining the GC content in angiosperm plastomes, with a particular focus on the hypothesis of phylogenetic dependence and the correlation with deletion mutations. We measured three genetic traits, namely, GC content, A/T tracts, and G/C tracts, in the coding region of plastid genomes for 1382 angiosperm species representing 350 families and 64 orders, and tested the phylogenetic signal. Then, we performed correlation analyses and revealed the variation in evolutionary rate of selected traits using RRphylo. The plastid GC content in the coding region varied from 28.10% to 43.20% across angiosperms, with a few non-photosynthetic species showing highly reduced values, highlighting the significance of functional constraints. We found strong phylogenetic signal in A/T tracts, but weak ones in GC content and G/C tracts, indicating adaptive potential. GC content was positively and negatively correlated with G/C and A/T tracts, respectively, suggesting a trade-off between these two deletion events. GC content evolved at various rates across the phylogeny, with significant increases in monocots and Lamiids, and a decrease in Fabids, implying the effects of some other factors. We hypothesize that variation in plastid GC content might be a mixed strategy of species to optimize fitness in fluctuating climates, partly through influencing the trade-off between AT → GC and GC → AT mutations.

The Implications of Lineage-Specific Rates for Divergence Time Estimation

Rate variation adds considerable complexity to divergence time estimation in molecular phylogenies. Here, we evaluate the impact of lineage-specific rates—which we define as among-branch-rate-variation that acts consistently across the entire genome. We compare its impact to residual rates—defined as among-branch-rate-variation that shows a different pattern of rate variation at each sampled locus, and gene-specific rates—defined as variation in the average rate across all branches at each sampled locus. We show that lineage-specific rates lead to erroneous divergence time estimates, regardless of how many loci are sampled. Further, we show that stronger lineage-specific rates lead to increasing error. This contrasts to residual rates and gene-specific rates, where sampling more loci significantly reduces error. If divergence times are inferred in a Bayesian framework, we highlight that error caused by lineage-specific rates significantly reduces the probability that the 95% highest posterior density includes the correct value, and leads to sensitivity to the prior. Use of a more complex rate prior—which has recently been proposed to model rate variation more accurately—does not affect these conclusions. Finally, we show that the scale of lineage-specific rates used in our simulation experiments is comparable to that of an empirical data set for the angiosperm genus Ipomoea. Taken together, our findings demonstrate that lineage-specific rates cause error in divergence time estimates, and that this error is not overcome by analyzing genomic scale multilocus data sets. [Divergence time estimation; error; rate variation.]

Uncertainty in Divergence Time Estimation

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

The rise of angiosperms pushed conifers to decline during global cooling

Competition among species and entire clades can impact species diversification and extinction, which can shape macroevolutionary patterns. The fossil record shows successive biotic turnovers such that a dominant group is replaced by another. One striking example involves the decline of gymnosperms and the rapid diversification and ecological dominance of angiosperms in the Cretaceous. It is generally believed that angiosperms outcompeted gymnosperms, but the macroevolutionary processes and alternative drivers explaining this pattern remain elusive. Using extant time trees and vetted fossil occurrences for conifers, we tested the hypotheses that clade competition or climate change led to the decline of conifers at the expense of angiosperms. Here, we find that both fossil and molecular data show high congruence in revealing 1) low diversification rates, punctuated by speciation pulses, during warming events throughout the Phanerozoic and 2) that conifer extinction increased significantly in the Mid-Cretaceous (100 to 110 Ma) and remained high ever since. Their extinction rates are best explained by the rise of angiosperms, rejecting alternative models based on either climate change or time alone. Our results support the hypothesis of an active clade replacement, implying that direct competition with angiosperms increased the extinction of conifers by pushing their remaining species diversity and dominance out of the warm tropics. This study illustrates how entire branches on the Tree of Life may actively compete for ecological dominance under changing climates.

The Longevity of Colonies of Fungus-Growing Termites and the Stability of the Symbiosis

The agricultural mutualistic symbiosis between macrotermitine termites and Termitomyces fungi is obligate for both partners. The termites provide a protective growth environment for the fungus by cultivating it inside their colony and providing it with foraged plant material. The termites use the fungus for plant substrate degradation, and the production of asexual fruiting bodies for nourishment and re-inoculation of the fungus garden. The termite colony can reach an age of up to several decades, during which time it is believed that a single fungal monoculture is asexually propagated by the offspring of a single founding royal pair. The termite-fungus mutualism has a long evolutionary history dating back more than 30 million years. Both on the time-scale of a termite colony lifespan and that of the mutualistic symbiosis, questions arise about stability. We address the physical stability of the mound, the termite colony and the monoculture fungal garden during a colony's lifetime. On the long-term evolutionary scale, we address the stability of the symbiosis, where horizontal transmission of the symbiotic fungus raises the question of how the mutualistic interaction between host and symbiont persists over generations.

Spatial Patterns and Drivers of Angiosperm Sexual Systems in China Differ Between Woody and Herbaceous Species

Plant sexual systems play an important role in the evolution of angiosperm diversity. However, large-scale patterns in the frequencies of sexual systems (i.e. dioecy, monoecy, and hermaphroditism) and their drivers for species with different growth forms remain poorly known. Here, using a newly compiled database on the sexual systems and distributions of 19780 angiosperm species in China, we map the large-scale geographical patterns in frequencies of the sexual systems of woody and herbaceous species separately. We use these data to test the following two hypotheses: (1) the prevalence of sexual systems differs between woody and herbaceous assemblies because woody plants have taller canopies and are found in warm and humid climates; (2) the relative contributions of different drivers (specifically climate, evolutionary age, and mature plant height) to these patterns differ between woody and herbaceous species. We show that geographical patterns in proportions of different sexual systems (especially dioecy) differ between woody and herbaceous species. Geographical variations in sexual systems of woody species were influenced by climate, evolutionary age and plant height. In contrast, these have only weakly significant effects on the patterns of sexual systems of herbaceous species. We suggest that differences between species with woody and herbaceous growth forms in terms of biogeographic patterns of sexual systems, and their drivers, may reflect their differences in physiological and ecological adaptions, as well as the coevolution of sexual system with vegetative traits in response to environmental changes.

Molecular Clocks without Rocks: New Solutions for Old Problems

Molecular data have been used to date species divergences ever since they were described as documents of evolutionary history in the 1960s. Yet, an inadequate fossil record and discordance between gene trees and species trees are persistently problematic. We examine how, by accommodating gene tree discordance and by scaling branch lengths to absolute time using mutation rate and generation time, multispecies coalescent (MSC) methods can potentially overcome these challenges. We find that time estimates can differ - in some cases, substantially - depending on whether MSC methods or traditional phylogenetic methods that apply concatenation are used, and whether the tree is calibrated with pedigree-based mutation rates or with fossils. We discuss the advantages and shortcomings of both approaches and provide practical guidance for data analysis when using these methods.

Heterogeneity in the rate of molecular sequence evolution substantially impacts the accuracy of detecting shifts in diversification rates

As species richness varies along the tree of life, there is a great interest in identifying factors that affect the rates by which lineages speciate or go extinct. To this end, theoretical biologists have developed a suite of phylogenetic comparative methods that aim to identify where shifts in diversification rates had occurred along a phylogeny and whether they are associated with some traits. Using these methods, numerous studies have predicted that speciation and extinction rates vary across the tree of life. In this study, we show that asymmetric rates of sequence evolution lead to systematic biases in the inferred phylogeny, which in turn lead to erroneous inferences regarding lineage diversification patterns. The results demonstrate that as the asymmetry in sequence evolution rates increases, so does the tendency to select more complicated models that include the possibility of diversification rate shifts. These results thus suggest that any inference regarding shifts in diversification pattern should be treated with great caution, at least until any biases regarding the molecular substitution rate have been ruled out.

Mutation-rate plasticity and the germline of unicellular organisms

The mutation rate is a fundamental factor in evolutionary genetics. Recently, mutation rates were found to be strongly reduced at high density in a wide range of unicellular organisms, prokaryotic and eukaryotic. Independently, cell division was found to become more asymmetrical at increasing density in diverse organisms; some 'mother' cells continue dividing, while their 'offspring' cells do not divide further. Here, we investigate how this increased asymmetry in cell division at high density can be reconciled with reduced mutation-rate estimates. We calculated the expected number of mutant cells due to replication errors under various modes of segregation of template-DNA strands and copy-DNA strands, both under symmetrical (exponential) and asymmetrical (linear) growth. We show that the observed reduction in the mutation rate at high density can be explained if mother cells preferentially retain the template-DNA strands, since new mutations are then confined to non-dividing daughter cells, thus reducing the spread of mutant cells. Any other inheritance mode results in an increase in the number of mutant cells at higher density. The proposed hypothesis that patterns of DNA-strand segregation are density-dependent fundamentally challenges our current understanding of mutation-rate estimates and extends the distinction between germline and soma to unicellular organisms.

Oak symbolism in the light of genomics

Throughout the Northern Hemisphere, human societies, political systems, and religions have appropriated oaks in symbolic representations. In this review, we explore the possible associations between recent genetic and genomic findings and the symbolic representations of oaks. We first consider the ways in which evolutionary history during the Holocene has tightened links between humans and oaks in Europe, and how this may have led to symbolic representations. We then show how recent findings concerning the structure and evolution of the oak genome have provided additional knowledge about symbolic representations, such as longevity, cohesiveness, and robustness.

A phylogenomic approach reveals a low somatic mutation rate in a long-lived plant

Somatic mutations can have important effects on the life history, ecology, and evolution of plants, but the rate at which they accumulate is poorly understood and difficult to measure directly. Here, we develop a method to measure somatic mutations in individual plants and use it to estimate the somatic mutation rate in a large, long-lived, phenotypically mosaic Eucalyptus melliodora tree. Despite being 100 times larger than Arabidopsis, this tree has a per-generation mutation rate only ten times greater, which suggests that this species may have evolved mechanisms to reduce the mutation rate per unit of growth. This adds to a growing body of evidence that illuminates the correlated evolutionary shifts in mutation rate and life history in plants.

-Complete plastid genome sequences of two species of the Neotropical genus Brunellia (Brunelliaceae)

Here we present the first two complete plastid genomes for Brunelliaceae, a Neotropical family with a single genus, Brunellia. We surveyed the entire plastid genome in order to find variable cpDNA regions for further phylogenetic analyses across the family. We sampled morphologically different species, B. antioquensis and B. trianae, and found that the plastid genomes are 157,685 and 157,775 bp in length and display the typical quadripartite structure found in angiosperms. Despite the clear morphological distinction between both species, the molecular data show a very low level of divergence. The amount of nucleotide substitutions per site is one of the lowest reported to date among published congeneric studies (π = 0.00025). The plastid genomes have gene order and content coincident with other COM (Celastrales, Oxalidales, Malpighiales) relatives. Phylogenetic analyses of selected superrosid representatives show high bootstrap support for the ((C,M)O) topology. The N-fixing clade appears as the sister group of the COM clade and Zygophyllales as the sister to the rest of the fabids group.

Meta-analysis Shows That Rapid Phenotypic Change in Angiosperms in Response to Environmental Change Is Followed by Stasis

The amount and rate of phenotypic change at ecological timescales varies widely, but there has not been a comprehensive quantitative synthesis of the patterns and causes of such variation for plants. Present knowledge is based predominantly on animals, whose differences with plants in the origin of germ cells and the level of modularity (among others) could make it invalid for plants. We synthesized data on contemporary phenotypic responses of angiosperms to environmental change and show that if extinction does not occur, quantitative traits change quickly in the first few years following the environmental novelty and then remain stable. This general pattern is independent from life span, growth form, spatial scale, or the type of trait. Our work shows that high amounts and rates of phenotypic change at contemporary timescales observed in plants are consistent with the pattern of stasis and bounded evolution previously observed over longer time frames. We also found evidence that may contradict some common ideas about phenotypic evolution: (1) the total amount of phenotypic change observed does not differ significantly according to growth form or life span; (2) greater and faster divergence tends to occur between populations connected at the local scale, where gene flow could be intense, rather than between distant populations; and (3) traits closely related to fitness change as much and as fast as other traits.

Somatic Mutation and Evolution in Plants

Somatic mutations are common in plants, and they may accumulate and be passed on to gametes. The determinants of somatic mutation accumulation include the intraorganismal selective effect of mutations, the number of cell divisions that separate the zygote from the formation of gametes, and shoot apical meristem structure and branching. Somatic mutations can promote the evolution of diploidy, polyploidy, sexual recombination, outcrossing, clonality, and separate sexes, and they may contribute genetic variability in many other traits. The amplification of beneficial mutations via intraorganismal selection may relax selection to reduce the genomic mutation rate or to protect the germline in plants. The total rate of somatic mutation, the distribution of selective effects and fates in the plant body, and the degree to which the germline is sheltered from somatic mutations are still poorly understood. Our knowledge can be improved through empirical estimates of mutation rates and effects on cell lineages and whole organisms, such as estimates of the reduction in fitness of progeny produced by within- versus between-flower crosses on the same plant, mutation coalescent studies within the canopy, and incorporation of somatic mutation into theoretical models of plant evolutionary genetics.

Maintenance of High Genome Integrity over Vegetative Growth in the Fairy-Ring Mushroom Marasmius oreades

Most mutations in coding regions of the genome are deleterious, causing selection to favor mechanisms that minimize the mutational load over time [1-5]. DNA replication during cell division is a major source of new mutations. It is therefore important to limit the number of cell divisions between generations, particularly for large and long-lived organisms [6-9]. The germline cells of animals and the slowly dividing cells in plant meristems are adaptations to control the number of mutations that accumulate over generations [9-11]. Fungi lack a separated germline while harboring species with very large and long-lived individuals that appear to maintain highly stable genomes within their mycelia [8, 12, 13]. Here, we studied genomic mutation accumulation in the fairy-ring mushroom Marasmius oreades. We generated a chromosome-level genome assembly using a combination of cutting-edge DNA sequencing technologies and re-sequenced 40 samples originating from six individuals of this fungus. The low number of mutations recovered in the sequencing data suggests the presence of an unknown mechanism that works to maintain extraordinary genome integrity over vegetative growth in M. oreades. The highly structured growth pattern of M. oreades allowed us to estimate the number of cell divisions leading up to each sample [14, 15], and from this data, we infer an incredibly low per mitosis mutation rate (3.8 × 10^-12 mutations per site and cell division) as one of several possible explanations for the low number of identified mutations.

Dates and rates in grape's plastomes: evolution in slow motion

The family Vitaceae includes the domesticated grapevine (Vitis vinifera), one of the most economically important crops in the world. Despite the importance of Vitaceae, there is still considerable controversy surrounding their phylogenetic relationships and evolutionary timescales. Moreover, variation in rates of molecular evolution among Vitaceae remains mostly unexplored. The present research aims to fill these knowledge gaps through the analysis of plastome sequences. Thirteen newly sequenced grape plastomes are presented and their phylogenetic relationships examined. Divergence times and absolute substitution rates are inferred under different molecular clocks by the analysis of 95 non-coding plastid regions and 43 representative accessions of the major lineages of Vitaceae. Furthermore, the phylogenetic informativeness of non-coding plastid regions is investigated. We find strong evidence in favor of the random local clock model and rate heterogeneity within Vitaceae. Substitution rates decelerate in Ampelocissus, Ampelopsis, Nekemias, Parthenocissus, Rhoicissus, and Vitis, with genus Vitis showing the lowest values up to a minimum of ~ 4.65 × 10^-11 s/s/y. We suggest that liana-like species of Vitaceae evolve slower than erect growth habit plants and we invoke the "rate of mitosis hypothesis" to explain the observed pattern of the substitution rates. We identify a reduced set of 20 non-coding regions able to accurately reconstruct the phylogeny of Vitaceae and we provide a detailed description of all 152 non-coding regions identified in the plastomes of subg. Vitis. These polymorphic regions will find their applications in phylogenetics, phylogeography, and population genetics as well in grapes identification through DNA barcoding techniques.

Optimal Rates for Phylogenetic Inference and Experimental Design in the Era of Genome-Scale Data Sets

With the rise of genome-scale data sets, there has been a call for increased data scrutiny and careful selection of loci that are appropriate to use in an attempt to resolve a phylogenetic problem. Such loci should maximize phylogenetic information content while minimizing the risk of homoplasy. Theory posits the existence of characters that evolve at an optimum rate, and efforts to determine optimal rates of inference have been a cornerstone of phylogenetic experimental design for over two decades. However, both theoretical and empirical investigations of optimal rates have varied dramatically in their conclusions: spanning no relationship to a tight relationship between the rate of change and phylogenetic utility. Herein, we synthesize these apparently contradictory views, demonstrating both empirical and theoretical conditions under which each is correct. We find that optimal rates of characters-not genes-are generally robust to most experimental design decisions. Moreover, consideration of site rate heterogeneity within a given locus is critical to accurate predictions of utility. Factors such as taxon sampling or the targeted number of characters providing support for a topology are additionally critical to the predictions of phylogenetic utility based on the rate of character change. Further, optimality of rates and predictions of phylogenetic utility are not equivalent, demonstrating the need for further development of comprehensive theory of phylogenetic experimental design. [Divergence time; GC bias; homoplasy; incongruence; information content; internode length; optimal rates; phylogenetic informativeness; phylogenetic theory; phylogenetic utility; phylogenomics; signal and noise; subtending branch length; state space; taxon and character sampling.].