Exploring the Relationships between Mutation Rates, Life History, Genome Size, Environment, and Species Richness in Flowering Plants

Effects of urban-induced mutations on ecology, evolution and health

Increasing evidence suggests that urbanization is associated with higher mutation rates, which can affect the health and evolution of organisms that inhabit cities. Elevated pollution levels in urban areas can induce DNA damage, leading to de novo mutations. Studies on mutations induced by urban pollution are most prevalent in humans and microorganisms, whereas studies of non-human eukaryotes are rare, even though increased mutation rates have the potential to affect organisms and their populations in contemporary time. Our Perspective explores how higher mutation rates in urban environments could impact the fitness, ecology and evolution of populations. Most mutations will be neutral or deleterious, and higher mutation rates associated with elevated pollution in urban populations can increase the risk of cancer in humans and potentially other species. We highlight the potential for urban-driven increased deleterious mutational loads in some organisms, which could lead to a decline in population growth of a wide diversity of organisms. Although beneficial mutations are expected to be rare, we argue that higher mutation rates in urban areas could influence adaptive evolution, especially in organisms with short generation times. Finally, we explore avenues for future research to better understand the effects of urban-induced mutations on the fitness, ecology and evolution of city-dwelling organisms.

Uncovering epigenetic and transcriptional regulation of growth in Douglas-fir: identification of differential methylation regions in mega-sized introns

Tree growth performance can be partly explained by genetics, while a large proportion of growth variation is thought to be controlled by environmental factors. However, to what extent DNA methylation, a stable epigenetic modification, contributes to phenotypic plasticity in the growth performance of long-lived trees remains unclear. In this study, a comparative analysis of targeted DNA genotyping, DNA methylation and mRNAseq profiling for needles of 44-year-old Douglas-fir trees (Pseudotsuga menziesii (Mirb.) Franco) having contrasting growth characteristics was performed. In total, we identified 195 differentially expressed genes (DEGs) and 115 differentially methylated loci (DML) that are associated with genes involved in fitness-related processes such as growth, stress management, plant development and energy resources. Interestingly, all four intronic DML were identified in mega-sized (between 100 and 180 kbp in length) and highly expressed genes, suggesting specialized regulation mechanisms of these long intron genes in gymnosperms. DNA repetitive sequences mainly comprising long-terminal repeats of retroelements are involved in growth-associated DNA methylation regulation (both hyper- and hypomethylation) of 99 DML (86.1% of total DML). Furthermore, nearly 14% of the DML was not tagged by single nucleotide polymorphisms, suggesting a unique contribution of the epigenetic variation in tree growth.

Half a Century of Controversy: The Neutralist/Selectionist Debate in Molecular Evolution

The neutral and nearly neutral theories, introduced more than 50 yr ago, have raised and still raise passionate discussion regarding the forces governing molecular evolution and their relative importance. The debate, initially focused on the amount of within-species polymorphism and constancy of the substitution rate, has spread, matured, and now underlies a wide range of topics and questions. The neutralist/selectionist controversy has structured the field and influences the way molecular evolutionary scientists conceive their research.

Evolution of the Araliaceae family involved rapid diversification of the Asian Palmate group and Hydrocotyle specific mutational pressure

The Araliaceae contain many valuable species in medicinal and industrial aspects. We performed intensive phylogenomics using the plastid genome (plastome) and 45S nuclear ribosomal DNA sequences. A total of 66 plastome sequences were used, 13 of which were newly assembled in this study, 12 from new sequences, and one from existing data. While Araliaceae plastomes showed conserved genome structure, phylogenetic reconstructions based on four different plastome datasets revealed phylogenetic discordance within the Asian Palmate group. The divergence time estimation revealed that splits in two Araliaceae subfamilies and the clades exhibiting phylogenetic discordances in the Asian Palmate group occurred at two climatic optima, suggesting that global warming events triggered species divergence, particularly the rapid diversification of the Asian Palmate group during the Middle Miocene. Nucleotide substitution analyses indicated that the Hydrocotyloideae plastomes have undergone accelerated AT-biased mutations (C-to-T transitions) compared with the Aralioideae plastomes, and the acceleration may occur in their mitochondrial and nuclear genomes as well. This implies that members of the genus Hydrocotyle, the only aquatic plants in the Araliaceae, have experienced a distinct evolutionary history from the other species. We also discussed the intercontinental disjunction in the genus Panax and proposed a hypothesis to complement the previously proposed hypothesis. Our results provide the evolutionary trajectory of Araliaceae and advance our current understanding of the evolution of Araliaceae species.

Tensor Decomposition-based Feature Extraction and Classification to Detect Natural Selection from Genomic Data

Inferences of adaptive events are important for learning about traits, such as human digestion of lactose after infancy and the rapid spread of viral variants. Early efforts toward identifying footprints of natural selection from genomic data involved development of summary statistic and likelihood methods. However, such techniques are grounded in simple patterns or theoretical models that limit the complexity of settings they can explore. Due to the renaissance in artificial intelligence, machine learning methods have taken center stage in recent efforts to detect natural selection, with strategies such as convolutional neural networks applied to images of haplotypes. Yet, limitations of such techniques include estimation of large numbers of model parameters under nonconvex settings and feature identification without regard to location within an image. An alternative approach is to use tensor decomposition to extract features from multidimensional data although preserving the latent structure of the data, and to feed these features to machine learning models. Here, we adopt this framework and present a novel approach termed T-REx, which extracts features from images of haplotypes across sampled individuals using tensor decomposition, and then makes predictions from these features using classical machine learning methods. As a proof of concept, we explore the performance of T-REx on simulated neutral and selective sweep scenarios and find that it has high power and accuracy to discriminate sweeps from neutrality, robustness to common technical hurdles, and easy visualization of feature importance. Therefore, T-REx is a powerful addition to the toolkit for detecting adaptive processes from genomic data.

Evolution of the Mutation Spectrum Across a Mammalian Phylogeny

Although evolutionary biologists have long theorized that variation in DNA repair efficacy might explain some of the diversity of lifespan and cancer incidence across species, we have little data on the variability of normal germline mutagenesis outside of humans. Here, we shed light on the spectrum and etiology of mutagenesis across mammals by quantifying mutational sequence context biases using polymorphism data from thirteen species of mice, apes, bears, wolves, and cetaceans. After normalizing the mutation spectrum for reference genome accessibility and k-mer content, we use the Mantel test to deduce that mutation spectrum divergence is highly correlated with genetic divergence between species, whereas life history traits like reproductive age are weaker predictors of mutation spectrum divergence. Potential bioinformatic confounders are only weakly related to a small set of mutation spectrum features. We find that clock-like mutational signatures previously inferred from human cancers cannot explain the phylogenetic signal exhibited by the mammalian mutation spectrum, despite the ability of these signatures to fit each species' 3-mer spectrum with high cosine similarity. In contrast, parental aging signatures inferred from human de novo mutation data appear to explain much of the 1-mer spectrum's phylogenetic signal in combination with a novel mutational signature. We posit that future models purporting to explain the etiology of mammalian mutagenesis need to capture the fact that more closely related species have more similar mutation spectra; a model that fits each marginal spectrum with high cosine similarity is not guaranteed to capture this hierarchy of mutation spectrum variation among species.

A high-quality Bougainvillea genome provides new insights into evolutionary history and pigment biosynthetic pathways in the Caryophyllales

Bougainvillea is a perennial ornamental shrub that is highly regarded in ornamental horticulture around the world. However, the absence of genome data limits our understanding of the pathways involved in bract coloration and breeding. Here, we report a chromosome-level assembly of the giga-genome of Bougainvillea × buttiana 'Mrs Butt', a cultivar thought to be the origin of many other Bougainvillea cultivars. The assembled genome is ~5 Gb with a scaffold N50 of 151 756 278 bp and contains 86 572 genes which have undergone recent whole-genome duplication. We confirmed that multiple rounds of whole-genome multiplication have occurred in the evolutionary history of the Caryophyllales, reconstructed the relationship in the Caryophyllales at whole genome level, and found discordance between species and gene trees as the result of complex introgression events. We investigated betalain and anthocyanin biosynthetic pathways and found instances of independent evolutionary innovations in the nine different Caryophyllales species. To explore the potential formation mechanism of diverse bract colors in Bougainvillea, we analyzed the genes involved in betalain and anthocyanin biosynthesis and found extremely low expression of ANS and DFR genes in all cultivars, which may limit anthocyanin biosynthesis. Our findings indicate that the expression pattern of the betalain biosynthetic pathway did not directly correlate with bract color, and a higher expression level in the betalain biosynthetic pathway is required for colored bracts. This improved understanding of the correlation between gene expression and bract color allows plant breeding outcomes to be predicted with greater certainty.

Fewer chromosomes, more co-occurring species within plant lineages: A likely effect of local survival and colonization

PREMISE

Plant lineages differ markedly in species richness globally, regionally, and locally. Differences in whole-genome characteristics (WGCs) such as monoploid chromosome number, genome size, and ploidy level may explain differences in global species richness through speciation or global extinction. However, it is unknown whether WGCs drive species richness within lineages also in a recent, postglacial regional flora or in local plant communities through local extinction or colonization and regional species turnover.

METHODS

We tested for relationships between WGCs and richness of angiosperm families across the Netherlands/Germany/Czechia as a region, and within 193,449 local vegetation plots.

RESULTS

Families that are species-rich across the region have lower ploidy levels and small monoploid chromosomes numbers or both (interaction terms), but the relationships disappear after accounting for continental and local richness of families. Families that are species-rich within occupied localities have small numbers of polyploidy and monoploid chromosome numbers or both, independent of their own regional richness and the local richness of all other locally co-occurring species in the plots. Relationships between WGCs and family species-richness persisted after accounting for niche characteristics and life histories.

CONCLUSIONS

Families that have few chromosomes, either monoploid or holoploid, succeed in maintaining many species in local communities and across a continent and, as indirect consequence of both, across a region. We suggest evolutionary mechanisms to explain how small chromosome numbers and ploidy levels might decrease rates of local extinction and increase rates of colonization. The genome of a macroevolutionary lineage may ultimately control whether its species can ecologically coexist.

Evolution of genome space occupation in ferns: linking genome diversity and species richness

BACKGROUND AND AIMS

The dynamics of genome evolution caused by whole genome duplications and other processes are hypothesized to shape the diversification of plants and thus contribute to the astonishing variation in species richness among the main lineages of land plants. Ferns, the second most species-rich lineage of land plants, are highly suitable to test this hypothesis because of several unique features that distinguish fern genomes from those of seed plants. In this study, we tested the hypothesis that genome diversity and disparity shape fern species diversity by recording several parameters related to genome size and chromosome number.

METHODS

We conducted de novo measurement of DNA C-values across the fern phylogeny to reconstruct the phylogenetic history of the genome space occupation in ferns by integrating genomic parameters such as genome size, chromosome number and average DNA amount per chromosome into a time-scaled phylogenetic framework. Using phylogenetic generalized least square methods, we determined correlations between chromosome number and genome size, species diversity and evolutionary rates of their transformation.

KEY RESULTS

The measurements of DNA C-values for 233 species more than doubled the taxon coverage from ~2.2 % in previous studies to 5.3 % of extant diversity. The dataset not only documented substantial differences in the accumulation of genomic diversity and disparity among the major lineages of ferns but also supported the predicted correlation between species diversity and the dynamics of genome evolution.

CONCLUSIONS

Our results demonstrated substantial genome disparity among different groups of ferns and supported the prediction that alterations of reproductive modes alter trends of genome evolution. Finally, we recovered evidence for a close link between the dynamics of genome evolution and species diversity in ferns for the first time.

Ten Plastomes of Crassula (Crassulaceae) and Phylogenetic Implications

The genus Crassula is the second-largest genus in the family Crassulaceae, with about 200 species. As an acknowledged super-barcode, plastomes have been extensively utilized for plant evolutionary studies. Here, we first report 10 new plastomes of Crassula. We further focused on the structural characterizations, codon usage, aversion patterns, and evolutionary rates of plastomes. The IR junction patterns-IRb had 110 bp expansion to rps19-were conservative among Crassula species. Interestingly, we found the codon usage patterns of matK gene in Crassula species are unique among Crassulaceae species with elevated ENC values. Furthermore, subgenus Crassula species have specific GC-biases in the matK gene. In addition, the codon aversion motifs from matK, pafI, and rpl22 contained phylogenetic implications within Crassula. The evolutionary rates analyses indicated all plastid genes of Crassulaceae were under the purifying selection. Among plastid genes, ycf1 and ycf2 were the most rapidly evolving genes, whereas psaC was the most conserved gene. Additionally, our phylogenetic analyses strongly supported that Crassula is sister to all other Crassulaceae species. Our findings will be useful for further evolutionary studies within the Crassula and Crassulaceae.

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.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02047-0.

Evolutionary importance of the relationship between cytogeography and climate: New insights on creosote bushes from North and South America

Relationships between genome size and environmental variables suggest that DNA content might be adaptive and of evolutionary importance in plants. The genus Larrea provides an interesting system to test this hypothesis, since it shows both intra- and interspecific variation in genome size. Larrea has an amphitropical distribution in North and South American deserts, where it is most speciose. Larrea tridentata in North America shows a gradient of increasing autopolyploidy; while three of the four studied South American species are diploids, Larrea divaricata, Larrea nitida, Larrea ameghinoi, and the fourth is an allopolyploid, Larrea cuneifolia. We downloaded available focal species' georeferenced records from seven data reservoirs. We used these records to extract biologically relevant environmental variables from WorldClim at 30 arc seconds scale, to have a broad characterization of the variable climatic conditions of both regions, and a climatic envelope for each species. We estimated relative DNA content index and relative monoploid genome values, by flow cytometry, of four most abundant Larrea species throughout their respective ranges. Then we winnow the bioclimatic dataset down to uncorrelated variables and sampled locales, to analyse the degree of association between both intra- and interspecific relative DNA content and climatic variables that are functionally relevant in arid environments using Pearson correlations, general linear and mixed effects models. Within the genus Larrea, relative DNA content increases with rising temperature and decreases with rising precipitation. At the intraspecific level, all four species show relative DNA content variation across climatic conditions. Larrea is a genus that shows genome size variation correlated with climate. Our results are also consistent with the hypothesis that extreme environmental pressures may have facilitated repeated whole genome duplication events in North America, while in South America, reticulate evolution, as allopolyploidization, and speciation might have been climate-dependent since the Oligocene.

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.

Chromosome-level assembly, annotation and phylome of Pelobates cultripes, the western spadefoot toad

Genomic resources for amphibians are still hugely under-represented in vertebrate genomic research, despite being a group of major interest for ecology, evolution and conservation. Amphibians constitute a highly threatened group of vertebrates, present a vast diversity in reproductive modes, are extremely diverse in morphology, occupy most ecoregions of the world, and present the widest range in genome sizes of any major group of vertebrates. We combined Illumina, Nanopore and Hi-C sequencing technologies to assemble a chromosome-level genome sequence for an anuran with a moderate genome size (assembly span 3.09 Gb); Pelobates cultripes, the western spadefoot toad. The genome has an N50 length of 330 Mb with 98.6% of the total sequence length assembled into 14 super scaffolds, and 87.7% complete BUSCO genes. We use published transcriptomic data to provide annotations, identifying 32,684 protein-coding genes. We also reconstruct the P. cultripes phylome and identify 2,527 gene expansions. We contribute the first draft of the genome of the western spadefoot toad, P. cultripes. This species represents a relatively basal lineage in the anuran tree with an interesting ecology and a high degree of developmental plasticity, and thus is an important resource for amphibian genomic research.

Investigating the reliability of molecular estimates of evolutionary time when substitution rates and speciation rates vary

Background

An accurate timescale of evolutionary history is essential to testing hypotheses about the influence of historical events and processes, and the timescale for evolution is increasingly derived from analysis of DNA sequences. But variation in the rate of molecular evolution complicates the inference of time from DNA. Evidence is growing for numerous factors, such as life history and habitat, that are linked both to the molecular processes of mutation and fixation and to rates of macroevolutionary diversification. However, the most widely used methods rely on idealised models of rate variation, such as the uncorrelated and autocorrelated clocks, and molecular dating methods are rarely tested against complex models of rate change. One relationship that is not accounted for in molecular dating is the potential for interaction between molecular substitution rates and speciation, a relationship that has been supported by empirical studies in a growing number of taxa. If these relationships are as widespread as current evidence suggests, they may have a significant influence on molecular dates.

Results

We simulate phylogenies and molecular sequences under three different realistic rate variation models—one in which speciation rates and substitution rates both vary but are unlinked, one in which they covary continuously and one punctuated model in which molecular change is concentrated in speciation events, using empirical case studies to parameterise realistic simulations. We test three commonly used “relaxed clock” molecular dating methods against these realistic simulations to explore the degree of error in molecular dates under each model. We find average divergence time inference errors ranging from 12% of node age for the unlinked model when reconstructed under an uncorrelated rate prior using BEAST 2, to up to 91% when sequences evolved under the punctuated model are reconstructed under an autocorrelated prior using PAML.

Conclusions

We demonstrate the potential for substantial errors in molecular dates when both speciation rates and substitution rates vary between lineages. This study highlights the need for tests of molecular dating methods against realistic models of rate variation generated from empirical parameters and known relationships.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-022-02015-8.

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.

Diversification Rate is Associated with Rate of Molecular Evolution in Ray-Finned Fish (Actinopterygii)

Understanding the factors that drive diversification of taxa across the tree of life is a key focus of macroevolutionary research. While the effects of life history, ecology, climate and geography on diversity have been studied for many taxa, the relationship between molecular evolution and diversification has received less attention. However, correlations between rates of molecular evolution and diversification rate have been detected in a range of taxa, including reptiles, plants and birds. A correlation between rates of molecular evolution and diversification rate is a prediction of several evolutionary theories, including the evolutionary speed hypothesis which links variation in mutation rates to differences in speciation rates. If it is widespread, such correlations could also have significant practical impacts, if they are not adequately accounted for in phylogenetic inference of evolutionary rates and timescales. Ray-finned fish (Actinopterygii) offer a prime target to test for this relationship due to their extreme variation in clade size suggesting a wide range of diversification rates. We employ both a sister-pairs approach and a whole-tree approach to test for correlations between substitution rate and net diversification. We also collect life history and ecological trait data and account for potential confounding factors including body size, latitude, max depth and reef association. We find evidence to support a relationship between diversification and synonymous rates of nuclear evolution across two published backbone phylogenies, as well as weak evidence for a relationship between mitochondrial nonsynonymous rates and diversification at the genus level.

From Genome Sequencing to CRISPR-Based Genome Editing for Climate-Resilient Forest Trees

Due to the economic and ecological importance of forest trees, modern breeding and genetic manipulation of forest trees have become increasingly prevalent. The CRISPR-based technology provides a versatile, powerful, and widely accepted tool for analyzing gene function and precise genetic modification in virtually any species but remains largely unexplored in forest species. Rapidly accumulating genetic and genomic resources for forest trees enabled the identification of numerous genes and biological processes that are associated with important traits such as wood quality, drought, or pest resistance, facilitating the selection of suitable gene editing targets. Here, we introduce and discuss the latest progress, opportunities, and challenges of genome sequencing and editing for improving forest sustainability.

Nucleotide substitutions during speciation may explain substitution rate variation

Although molecular mechanisms associated with the generation of mutations are highly conserved across taxa, there is widespread variation in mutation rates between evolutionary lineages. When phylogenies are reconstructed based on nucleotide sequences, such variation is typically accounted for by the assumption of a relaxed molecular clock, which is a statistical distribution of mutation rates without much underlying biological mechanism. Here, we propose that variation in accumulated mutations may be partly explained by an elevated mutation rate during speciation. Using simulations, we show how shifting mutations from branches to speciation events impacts inference of branching times in phylogenetic reconstruction. Furthermore, the resulting nucleotide alignments are better described by a relaxed than by a strict molecular clock. Thus, elevated mutation rates during speciation potentially explain part of the variation in substitution rates that is observed across the tree of life. [Molecular clock; phylogenetic reconstruction; speciation; substitution rate variation.]

Evolutionary patterns of nucleotide substitution rates in plastid genomes of Quercus

Molecular evolution, including nucleotide substitutions, plays an important role in understanding the dynamics and mechanisms of species evolution. Here, we sequenced whole plastid genomes (plastomes) of Quercus fabri, Quercus semecarpifolia, Quercus engleriana, and Quercus phellos and compared them with 14 other Quercus plastomes to explore their evolutionary relationships using 67 shared protein-coding sequences. While many previously identified evolutionary relationships were found, our findings do not support previous research which retrieve Quercus subg. Cerris sect. Ilex as a monophyletic group, with sect. Ilex found to be polyphyletic and composed of three strongly supported lineages inserted between sections Cerris and Cyclobalanposis. Compared with gymnosperms, Quercus plastomes showed higher evolutionary rates (Dn/Ds = 0.3793). Most protein-coding genes experienced relaxed purifying selection, and the high Dn value (0.1927) indicated that gene functions adjusted to environmental changes effectively. Our findings suggest that gene interval regions play an important role in Quercus evolution. We detected greater variation in the intergenic regions (trnH-psbA, trnK_UUU-rps16, trnfM_CAU-rps14, trnS_GCU-trnG_GCC, and atpF-atpH), intron losses (petB and petD), and pseudogene loss and degradation (ycf15). Additionally, the loss of some genes suggested the existence of gene exchanges between plastid and nuclear genomes, which affects the evolutionary rate of the former. However, the connective mechanism between these two genomes is still unclear.

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.

Evolutionary Responses to Warming

Climate change is predicted to dramatically alter biological diversity and distributions, driving extirpations, extinctions, and extensive range shifts across the globe. Warming can also, however, lead to phenotypic or behavioural plasticity, as species adapt to new conditions. Recent genomic research indicates that some species are capable of rapid evolution as selection favours adaptive responses to environmental change and altered or novel niche spaces. New advances are providing mechanistic insights into how temperature might accelerate evolution in the Anthropocene. These discoveries highlight intriguing new research directions - such as using geothermal and polar systems combined with powerful genomic tools - that will help us to understand the processes underpinning adaptive evolution and better project how ecosystems will change in a warming world.

The Effects of Ecological Traits on the Rate of Molecular Evolution in Ray-Finned Fishes: A Multivariable Approach

Myriad environmental and biological traits have been investigated for their roles in influencing the rate of molecular evolution across various taxonomic groups. However, most studies have focused on a single trait, while controlling for additional factors in an informal way, generally by excluding taxa. This study utilized a dataset of cytochrome c oxidase subunit I (COI) barcode sequences from over 7000 ray-finned fish species to test the effects of 27 traits on molecular evolutionary rates. Environmental traits such as temperature were considered, as were traits associated with effective population size including body size and age at maturity. It was hypothesized that these traits would demonstrate significant correlations with substitution rate in a multivariable analysis due to their associations with mutation and fixation rates, respectively. A bioinformatics pipeline was developed to assemble and analyze sequence data retrieved from the Barcode of Life Data System (BOLD) and trait data obtained from FishBase. For use in phylogenetic regression analyses, a maximum likelihood tree was constructed from the COI sequence data using a multi-gene backbone constraint tree covering 71% of the species. A variable selection method that included both single- and multivariable analyses was used to identify traits that contribute to rate heterogeneity estimated from different codon positions. Our analyses revealed that molecular rates associated most significantly with latitude, body size, and habitat type. Overall, this study presents a novel and systematic approach for integrative data assembly and variable selection methodology in a phylogenetic framework.

Nucleotide substitution rates of diatom plastid encoded protein genes are positively correlated with genome architecture

Diatoms are the largest group of heterokont algae with more than 100,000 species. As one of the single-celled photosynthetic organisms that inhabit marine, aquatic and terrestrial ecosystems, diatoms contribute ~ 45% of global primary production. Despite their ubiquity and environmental significance, very few diatom plastid genomes (plastomes) have been sequenced and studied. This study explored patterns of nucleotide substitution rates of diatom plastids across the entire suite of plastome protein-coding genes for 40 taxa representing the major clades. The highest substitution rate was lineage-specific within the araphid 2 taxon Astrosyne radiata and radial 2 taxon Proboscia sp. Rate heterogeneity was also evident in different functional classes and individual genes. Similar to land plants, proteins genes involved in photosynthetic metabolism have lower synonymous and nonsynonymous substitutions rates than those involved in transcription and translation. Significant positive correlations were identified between substitution rates and measures of genomic rearrangements, including indels and inversions, which is a similar result to what was found in legume plants. This work advances the understanding of the molecular evolution of diatom plastomes and provides a foundation for future studies.

The genetics of evolutionary radiations

With the realization that much of the biological diversity on Earth has been generated by discrete evolutionary radiations, there has been a rapid increase in research into the biotic (key innovations) and abiotic (key environments) circumstances in which such radiations took place. Here we focus on the potential importance of population genetic structure and trait genetic architecture in explaining radiations. We propose a verbal model describing the stages of an evolutionary radiation: first invading a suitable adaptive zone and expanding both spatially and ecologically through this zone; secondly, diverging genetically into numerous distinct populations; and, finally, speciating. There are numerous examples of the first stage; the difficulty, however, is explaining how genetic diversification can take place from the establishment of a, presumably, genetically depauperate population in a new adaptive zone. We explore the potential roles of epigenetics and transposable elements (TEs), of neutral process such as genetic drift in combination with trait genetic architecture, of gene flow limitation through isolation by distance (IBD), isolation by ecology and isolation by colonization, the possible role of intra-specific competition, and that of admixture and hybridization in increasing the genetic diversity of the founding populations. We show that many of the predictions of this model are corroborated. Most radiations occur in complex adaptive zones, which facilitate the establishment of many small populations exposed to genetic drift and divergent selection. We also show that many radiations (especially those resulting from long-distance dispersal) were established by polyploid lineages, and that many radiating lineages have small genome sizes. However, there are several other predictions which are not (yet) possible to test: that epigenetics has played a role in radiations, that radiations occur more frequently in clades with small gene flow distances, or that the ancestors of radiations had large fundamental niches. At least some of these may be testable in the future as more genome and epigenome data become available. The implication of this model is that many radiations may be hard polytomies because the genetic divergence leading to speciation happens within a very short time, and that the divergence history may be further obscured by hybridization. Furthermore, it suggests that only lineages with the appropriate genetic architecture will be able to radiate, and that such a radiation will happen in a meta-population environment. Understanding the genetic architecture of a lineage may be an essential part of accounting for why some lineages radiate, and some do not.

Sexual selection, body mass and molecular evolution interact to predict diversification in birds

Sexual selection is a powerful agent of evolution, driving microevolutionary changes in the genome and macroevolutionary rates of lineage diversification. The mechanisms by which sexual selection might influence macroevolution remain poorly understood. For example, sexual selection might drive positive selection for key adaptations that facilitate diversification. Furthermore, sexual selection might be a general driver of molecular evolutionary rate. We lay out some of the potential mechanisms that create a link between sexual selection and diversification, based on causal effects on other life-history traits such as body mass and the rate of molecular evolution. Birds are ideally suited for testing the importance of these relationships because of their diverse reproductive systems and the multiple evolutionary radiations that have produced their astounding modern diversity. We show that sexual selection (measured as the degree of polygyny) interacts with the rate of molecular evolution and with body mass to predict species richness at the genus level. A high degree of polygyny and rapid molecular evolution are positively associated with the net rate of diversification, with the two factors being especially important for explaining diversification in large-bodied taxa. Our findings further suggest that mutation rates underpin some of the macroevolutionary effects of sexual selection. We synthesize the existing theory on sexual selection as a force for diversity and propose avenues for exploring this association using genome data.

Disentangling the intertwined roles of mutation, selection and drift in the mitochondrial genome

Understanding and quantifying the rates of change in the mitochondrial genome is a major component of many areas of biological inquiry, from phylogenetics to human health. A critical parameter in understanding rates of change is estimating the mitochondrial mutation rate (mtDNA MR). Although the first direct estimates of mtDNA MRs were reported almost 20 years ago, the number of estimates has not grown markedly since that time. This is largely owing to the challenges associated with time- and labour-intensive mutation accumulation (MA) experiments. But even MA experiments do not solve a major problem with estimating mtDNA MRs-the challenge of disentangling the role of mutation from other evolutionary forces acting within the cell. Now that it is widely understood that any newly generated mutant allele in the mitochondria will initially be at very low frequency (1/N, where N is the number of mtDNA molecules in the cell), the importance of understanding the effective population size (Ne) of the mtDNA and the size of genetic bottlenecks during gametogenesis and development has come into the spotlight. In addition to these factors regulating the role of genetic drift, advances in our understanding of mitochondrial replication and turnover allow us to more easily envision how natural selection within the cell might favour or purge mutations in multi-copy organellar genomes. Here, we review the unique features of the mitochondrial genome that pose a challenge for accurate MR estimation and discuss ways to overcome those challenges. Estimates of mtDNA MRs remain one of the most widely used parameters in biology, thus accurate quantification and a deeper understanding of how and why they may vary within and between individuals, populations and species is an important goal. This article is part of the theme issue 'Linking the mitochondrial genotype to phenotype: a complex endeavour'.

The C-Fern (Ceratopteris richardii) genome: insights into plant genome evolution with the first partial homosporous fern genome assembly

Ferns are notorious for possessing large genomes and numerous chromosomes. Despite decades of speculation, the processes underlying the expansive genomes of ferns are unclear, largely due to the absence of a sequenced homosporous fern genome. The lack of this crucial resource has not only hindered investigations of evolutionary processes responsible for the unusual genome characteristics of homosporous ferns, but also impeded synthesis of genome evolution across land plants. Here, we used the model fern species Ceratopteris richardii to address the processes (e.g., polyploidy, spread of repeat elements) by which the large genomes and high chromosome numbers typical of homosporous ferns may have evolved and have been maintained. We directly compared repeat compositions in species spanning the green plant tree of life and a diversity of genome sizes, as well as both short- and long-read-based assemblies of Ceratopteris. We found evidence consistent with a single ancient polyploidy event in the evolutionary history of Ceratopteris based on both genomic and cytogenetic data, and on repeat proportions similar to those found in large flowering plant genomes. This study provides a major stepping-stone in the understanding of land plant evolutionary genomics by providing the first homosporous fern reference genome, as well as insights into the processes underlying the formation of these massive genomes.

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.

Plastid Genomes of Five Species of Riverweeds (Podostemaceae): Structural Organization and Comparative Analysis in Malpighiales

With the advent of next-generation sequencing technologies, whole-plastome data can be obtained as a byproduct of low-coverage sequencing of the plant genomic DNA. This provides an opportunity to study plastid evolution across groups, as well as testing phylogenetic relationships among taxa. Within the order Malpighiales (∼16,000 spp.), the Podostemaceae (∼300 spp.) stand out for their unique habit, living attached to rocks in fast-flowing aquatic habitats, and displaying highly modified morphologies that confound our understanding of their classification, biology, and evolution. In this study, we used genome skimming data to assemble the full plastid genome of 5 species within Podostemaceae. We analyzed our data in a comparative framework within Malpighiales to determine the structure, gene content, and rearrangements in the plastomes of the family. The Podostemaceae have one of the smallest plastid genomes reported so far for the Malpighiales, possibly due to variation in length of inverted repeat (IR) regions, gene loss, and intergenic region variation. We also detected a major inversion in the large single-copy region unique to the family. The uncommon loss or pseudogenization of ycf1 and ycf2 in angiosperms and in land plants in general is also found to be characteristic of Podostemaceae, but the compensatory mechanisms and implications of this and of the pseudogenization of accD, rpl22, and clpP and loss of rps16 remain to be explained in this group. In addition, we estimated a phylogenetic tree among selected species in Malpighiales. Our findings indicate that the Podostemaceae are a distinct lineage with long branches that suggest faster rates of evolution in the plastome of the group, compared with other taxa in the order. This study lays the foundations for future phylogenomic studies in the family.

Evolutionary Rates of Bumblebee Genomes Are Faster at Lower Elevations

The importance of climate in determining biodiversity patterns has been well documented. However, the relationship between climate and rates of genetic evolution remains controversial. Latitude and elevation have been associated with rates of change in genetic markers such as cytochrome b. What is not known, however, is the strength of such associations and whether patterns found among these genes apply across entire genomes. Here, using bumblebee genetic data from seven subgenera of Bombus, we demonstrate that all species occupying warmer elevations have undergone faster genome-wide evolution than those in the same subgenera occupying cooler elevations. Our findings point to a critical biogeographic role in the relative rates of whole species evolution, potentially influencing global biodiversity patterns.

Fire as a key driver of Earth's biodiversity

Many terrestrial ecosystems are fire prone, such that their composition and structure are largely due to their fire regime. Regions subject to regular fire have exceptionally high levels of species richness and endemism, and fire has been proposed as a major driver of their diversity, within the context of climate, resource availability and environmental heterogeneity. However, current fire-management practices rarely take into account the ecological and evolutionary roles of fire in maintaining biodiversity. Here, we focus on the mechanisms that enable fire to act as a major ecological and evolutionary force that promotes and maintains biodiversity over numerous spatiotemporal scales. From an ecological perspective, the vegetation, topography and local weather conditions during a fire generate a landscape with spatial and temporal variation in fire-related patches (pyrodiversity), and these produce the biotic and environmental heterogeneity that drives biodiversity across local and regional scales. There have been few empirical tests of the proposition that 'pyrodiversity begets biodiversity' but we show that biodiversity should peak at moderately high levels of pyrodiversity. Overall species richness is greatest immediately after fire and declines monotonically over time, with postfire successional pathways dictated by animal habitat preferences and varying lifespans among resident plants. Theory and data support the 'intermediate disturbance hypothesis' when mean patch species diversity is correlated with mean fire intervals. Postfire persistence, recruitment and immigration allow species with different life histories to coexist. From an evolutionary perspective, fire drives population turnover and diversification by promoting a wide range of adaptive responses to particular fire regimes. Among 39 comparisons, the number of species in 26 fire-prone lineages is much higher than that in their non-fire-prone sister lineages. Fire and its byproducts may have direct mutagenic effects, producing novel genotypes that can lead to trait innovation and even speciation. A paradigm shift aimed at restoring biodiversity-maintaining fire regimes across broad landscapes is required among the fire research and management communities. This will require ecologists and other professionals to spread the burgeoning fire-science knowledge beyond scientific publications to the broader public, politicians and media.

Somatic mutations substantially increase the per-generation mutation rate in the conifer Picea sitchensis

The rates and biological significance of somatic mutations have long been a subject of debate. Somatic mutations in plants are expected to accumulate with vegetative growth and time, yet rates of somatic mutations are unknown for conifers, which can reach exceptional sizes and ages. We investigated somatic mutation rates in the conifer Sitka spruce (Picea sitchensis (Bong.) Carr.) by analyzing a total of 276 Gb of nuclear DNA from the tops and bottoms of 20 old‐growth trees averaging 76 m in height. We estimate a somatic base substitution rate of 2.7 × 10⁻⁸ per base pair within a generation. To date, this is one of the highest estimated per‐generation rates of mutation among eukaryotes, indicating that somatic mutations contribute substantially to the total per‐generation mutation rate in conifers. Nevertheless, as the sampled trees are centuries old, the per‐year rate is low in comparison with nontree taxa. We argue that although somatic mutations raise genetic load in conifers, they generate important genetic variation and enable selection both among cell lineages within individual trees and among offspring.

Six Impossible Things before Breakfast: Assumptions, Models, and Belief in Molecular Dating

Confidence in molecular dating analyses has grown with the increasing sophistication of the methods. Some problematic cases where molecular dates disagreed with paleontological estimates appear to have been resolved with a growing agreement between molecules and fossils. But we cannot relax just yet. The growing analytical sophistication of many molecular dating methods relies on an increasingly large number of assumptions about evolutionary history and processes. Many of these assumptions are based on statistical tractability rather than being informed by improved understanding of molecular evolution, yet changing the assumptions can influence molecular dates. How can we tell if the answers we get are driven more by the assumptions we make than by the molecular data being analyzed?

Phylogenetic trends and environmental correlates of nuclear genome size variation in Helianthus sunflowers

Flowering plants serve as a powerful model for studying the evolution of nuclear genome size (GS) given the tremendous GS variation that exists both within and across angiosperm lineages. Helianthus sunflowers consist of c. 50 species native to North America that occupy diverse habitats and vary in ploidy level. In the current study, we generated a comprehensive GS database for 49 Helianthus species using flow cytometric approaches. We examined variability across the genus and present a comparative phylogenetic analysis of GS evolution in diploid Helianthus species. Results demonstrated that different clades of diploid Helianthus species showed evolutionary patterns of GS contraction, expansion and relative stasis, with annual diploid species evolving smaller GS with the highest rate of evolution. Phylogenetic comparative analyses of diploids revealed significant negative associations of GS with temperature seasonality and cell production rate, indicating that the evolution of larger GS in Helianthus diploids may be more permissible in habitats with longer growing seasons where selection for more rapid growth may be relaxed. The Helianthus GS database presented here and corresponding analyses of environmental and phenotypic correlates will facilitate ongoing and future research on the ultimate drivers of GS evolution in this well-studied North American plant genus.

Genome Evolution in Outcrossing vs. Selfing vs. Asexual Species

A major current molecular evolution challenge is to link comparative genomic patterns to species' biology and ecology. Breeding systems are pivotal because they affect many population genetic processes and thus genome evolution. We review theoretical predictions and empirical evidence about molecular evolutionary processes under three distinct breeding systems-outcrossing, selfing, and asexuality. Breeding systems may have a profound impact on genome evolution, including molecular evolutionary rates, base composition, genomic conflict, and possibly genome size. We present and discuss the similarities and differences between the effects of selfing and clonality. In reverse, comparative and population genomic data and approaches help revisiting old questions on the long-term evolution of breeding systems.

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

The Deep Past Controls the Phylogenetic Structure of Present, Local Communities

Coexisting species may be evolutionarily proximate or distant, resulting in phylogenetically poor or rich communities. This variation is often considered to result from present assembly processes. We argue that, under certain conditions, deep-past processes might control the phylogenetic diversity of communities. First, deep-past effects involve macroevolutionary processes, such as diversification rate, niche conservatism, or dispersal, in the lineages that constitute communities. Second, deep-past processes in the respective region or in the habitat type play a role, for instance, through age, area, stability, or connectivity. Third, the deep past may affect communities via trophic interactions (i.e., communities of enemies or mutualists or communities of hosts). We suggest that deep-past effects can be identified in local communities by measuring phylogenetic diversity in different species pools. We also show how community phylogenetic diversity results in positive or negative eco-evolutionary feedback, and we identify present-day conservation challenges that may profit from a deep-time perspective.

Is molecular evolution faster in the tropics?

The evolutionary speed hypothesis (ESH) suggests that molecular evolutionary rates are higher among species inhabiting warmer environments. Previously, the ESH has been investigated using small numbers of latitudinally-separated sister lineages; in animals, these studies typically focused on subsets of Chordata and yielded mixed support for the ESH. This study analyzed public DNA barcode sequences from the cytochrome c oxidase subunit I (COI) gene for six of the largest animal phyla (Arthropoda, Chordata, Mollusca, Annelida, Echinodermata, and Cnidaria) and paired latitudinally-separated taxa together informatically. Of 8037 lineage pairs, just over half (51.6%) displayed a higher molecular rate in the lineage inhabiting latitudes closer to the equator, while the remainder (48.4%) displayed a higher rate in the higher-latitude lineage. To date, this study represents the most comprehensive analysis of latitude-related molecular rate differences across animals. While a statistically-significant pattern was detected from our large sample size, our findings suggest that the EHS may not serve as a strong universal mechanism underlying the latitudinal diversity gradient and that COI molecular clocks may generally be applied across latitudes. This study also highlights the merits of using automation to analyze large DNA barcode datasets.

Exploring the Limits and Causes of Plastid Genome Expansion in Volvocine Green Algae

Plastid genomes are not normally celebrated for being large. But researchers are steadily uncovering algal lineages with big and, in rare cases, enormous plastid DNAs (ptDNAs), such as volvocine green algae. Plastome sequencing of five different volvocine species has revealed some of the largest, most repeat-dense plastomes on record, including that of Volvox carteri (∼525 kb). Volvocine algae have also been used as models for testing leading hypotheses on organelle genome evolution (e.g., the mutational hazard hypothesis), and it has been suggested that ptDNA inflation within this group might be a consequence of low mutation rates and/or the transition from a unicellular to multicellular existence. Here, we further our understanding of plastome size variation in the volvocine line by examining the ptDNA sequences of the colonial species Yamagishiella unicocca and Eudorina sp. NIES-3984 and the multicellular Volvox africanus, which are phylogenetically situated between species with known ptDNA sizes. Although V. africanus is closely related and similar in multicellular organization to V. carteri, its ptDNA was much less inflated than that of V. carteri. Synonymous- and noncoding-site nucleotide substitution rate analyses of these two Volvox ptDNAs suggest that there are drastically different plastid mutation rates operating in the coding versus intergenic regions, supporting the idea that error-prone DNA repair in repeat-rich intergenic spacers is contributing to genome expansion. Our results reinforce the idea that the volvocine line harbors extremes in plastome size but ultimately shed doubt on some of the previously proposed hypotheses for ptDNA inflation within the lineage.

Correlates of rate heterogeneity in avian ecomorphological traits

Heterogeneity in rates of trait evolution is widespread, but it remains unclear which processes drive fast and slow character divergence across global radiations. Here, we test multiple hypotheses for explaining rate variation in an ecomorphological trait (beak shape) across a globally distributed group (birds). We find low support that variation in evolutionary rates of species is correlated with life history, environmental mutagenic factors, range size, number of competitors, or living on islands. Indeed, after controlling for the negative effect of species' age, 80% of variation in species‐specific evolutionary rates remains unexplained. At the clade level, high evolutionary rates are associated with unusual phenotypes or high species richness. Taken together, these results imply that macroevolutionary rates of ecomorphological traits are governed by both ecological opportunity in distinct adaptive zones and niche differentiation among closely related species.

RNA-seq highlights parallel and contrasting patterns in the evolution of the nuclear genome of fully mycoheterotrophic plants

BACKGROUND

While photosynthesis is the most notable trait of plants, several lineages of plants (so-called full heterotrophs) have adapted to obtain organic compounds from other sources. The switch to heterotrophy leads to profound changes at the morphological, physiological and genomic levels.

RESULTS

Here, we characterize the transcriptomes of three species representing two lineages of mycoheterotrophic plants: orchids (Epipogium aphyllum and Epipogium roseum) and Ericaceae (Hypopitys monotropa). Comparative analysis is used to highlight the parallelism between distantly related fully heterotrophic plants. In both lineages, we observed genome-wide elimination of nuclear genes that encode proteins related to photosynthesis, while systems associated with protein import to plastids as well as plastid transcription and translation remain active. Genes encoding components of plastid ribosomes that have been lost from the plastid genomes have not been transferred to the nuclear genomes; instead, some of the encoded proteins have been substituted by homologs. The nuclear genes of both Epipogium species accumulated nucleotide substitutions twice as rapidly as their photosynthetic relatives; in contrast, no increase in the substitution rate was observed in H. monotropa.

CONCLUSIONS

Full heterotrophy leads to profound changes in nuclear gene content. The observed increase in the rate of nucleotide substitutions is lineage specific, rather than a universal phenomenon among non-photosynthetic plants.

Pinaceae show elevated rates of gene turnover that are robust to incomplete gene annotation

Gene duplications and gene losses are major determinants of genome evolution and phenotypic diversity. The frequency of gene turnover (gene gains and gene losses combined) is known to vary between organisms. Comparative genomic analyses of gene families can highlight such variation; however, estimates of gene turnover may be biased when using highly fragmented genome assemblies resulting in poor gene annotations. Here, we address potential biases introduced by gene annotation errors in estimates of gene turnover frequencies in a dataset including both well-annotated angiosperm genomes and the incomplete gene sets of four Pinaceae, including two pine species, Norway spruce and Douglas-fir. We show that Pinaceae experienced higher gene turnover rates than angiosperm lineages lacking recent whole-genome duplications. This finding is robust to both known major issues in Pinaceae gene sets: missing gene models and erroneous annotation of pseudogenes. A separate analysis limited to the four Pinaceae gene sets pointed to an accelerated gene turnover rate in pines compared with Norway spruce and Douglas-fir. Our results indicate that gene turnover significantly contributes to genome variation and possibly to speciation in Pinaceae, particularly in pines. Moreover, these findings indicate that reliable estimates of gene turnover frequencies can be discerned in incomplete and potentially inaccurate gene sets. Because gymnosperms are known to exhibit low overall substitution rates compared with angiosperms, our results suggest that the rate of single-base pair mutations is uncoupled from the rate of large DNA duplications and deletions associated with gene turnover in Pinaceae.

Population Size and the Rate of Language Evolution: A Test Across Indo-European, Austronesian, and Bantu Languages

What role does speaker population size play in shaping rates of language evolution? There has been little consensus on the expected relationship between rates and patterns of language change and speaker population size, with some predicting faster rates of change in smaller populations, and others expecting greater change in larger populations. The growth of comparative databases has allowed population size effects to be investigated across a wide range of language groups, with mixed results. One recent study of a group of Polynesian languages revealed greater rates of word gain in larger populations and greater rates of word loss in smaller populations. However, that test was restricted to 20 closely related languages from small Oceanic islands. Here, we test if this pattern is a general feature of language evolution across a larger and more diverse sample of languages from both continental and island populations. We analyzed comparative language data for 153 pairs of closely-related sister languages from three of the world's largest language families: Austronesian, Indo-European, and Niger-Congo. We find some evidence that rates of word loss are significantly greater in smaller languages for the Indo-European comparisons, but we find no significant patterns in the other two language families. These results suggest either that the influence of population size on rates and patterns of language evolution is not universal, or that it is sufficiently weak that it may be overwhelmed by other influences in some cases. Further investigation, for a greater number of language comparisons and a wider range of language features, may determine which of these explanations holds true.

Genome Size Diversity and Its Impact on the Evolution of Land Plants

Genome size is a biodiversity trait that shows staggering diversity across eukaryotes, varying over 64,000-fold. Of all major taxonomic groups, land plants stand out due to their staggering genome size diversity, ranging ca. 2400-fold. As our understanding of the implications and significance of this remarkable genome size diversity in land plants grows, it is becoming increasingly evident that this trait plays not only an important role in shaping the evolution of plant genomes, but also in influencing plant community assemblages at the ecosystem level. Recent advances and improvements in novel sequencing technologies, as well as analytical tools, make it possible to gain critical insights into the genomic and epigenetic mechanisms underpinning genome size changes. In this review we provide an overview of our current understanding of genome size diversity across the different land plant groups, its implications on the biology of the genome and what future directions need to be addressed to fill key knowledge gaps.