Modeling Substitution Rate Evolution across Lineages and Relaxing the Molecular Clock

Relaxing the molecular clock using models of how substitution rates change across lineages has become essential for addressing evolutionary problems. The diversity of rate evolution models and their implementations are substantial, and studies have demonstrated their impact on divergence time estimates can be as significant as that of calibration information. In this review, we trace the development of rate evolution models from the proposal of the molecular clock concept to the development of sophisticated Bayesian and non-Bayesian methods that handle rate variation in phylogenies. We discuss the various approaches to modeling rate evolution, provide a comprehensive list of available software, and examine the challenges and advancements of the prevalent Bayesian framework, contrasting them to faster non-Bayesian methods. Lastly, we offer insights into potential advancements in the field in the era of big data.

Pyricularia oryzae: Lab star and field scourge

Pyricularia oryzae (syn. Magnaporthe oryzae), is a filamentous ascomycete that causes a major disease called blast on cereal crops, as well as on a wide variety of wild and cultivated grasses. Blast diseases have a tremendous impact worldwide particularly on rice and on wheat, where the disease emerged in South America in the 1980s, before spreading to Asia and Africa. Its economic importance, coupled with its amenability to molecular and genetic manipulation, have inspired extensive research efforts aiming at understanding its biology and evolution. In the past 40 years, this plant-pathogenic fungus has emerged as a major model in molecular plant-microbe interactions. In this review, we focus on the clarification of the taxonomy and genetic structure of the species and its host range determinants. We also discuss recent molecular studies deciphering its lifecycle.

TAXONOMY

Kingdom: Fungi, phylum: Ascomycota, sub-phylum: Pezizomycotina, class: Sordariomycetes, order: Magnaporthales, family: Pyriculariaceae, genus: Pyricularia.

HOST RANGE

P. oryzae has the ability to infect a wide range of Poaceae. It is structured into different host-specialized lineages that are each associated with a few host plant genera. The fungus is best known to cause tremendous damage to rice crops, but it can also attack other economically important crops such as wheat, maize, barley, and finger millet.

DISEASE SYMPTOMS

P. oryzae can cause necrotic lesions or bleaching on all aerial parts of its host plants, including leaf blades, sheaths, and inflorescences (panicles, spikes, and seeds). Characteristic symptoms on leaves are diamond-shaped silver lesions that often have a brown margin and whose appearance is influenced by numerous factors such as the plant genotype and environmental conditions. USEFUL WEBSITES Resources URL Genomic data repositories http://genome.jouy.inra.fr/gemo/ Genomic data repositories http://openriceblast.org/ Genomic data repositories http://openwheatblast.net/ Genome browser for fungi (including P. oryzae) http://fungi.ensembl.org/index.html Comparative genomics database https://mycocosm.jgi.doe.gov/mycocosm/home T-DNA mutant database http://atmt.snu.kr/ T-DNA mutant database http://www.phi-base.org/ SNP and expression data https://fungidb.org/fungidb/app/.

Evolution of an intermediate C4 photosynthesis in the non-foliar tissues of the Poaceae

Carbon concentrating mechanisms (CCMs) in plants are abaptive features that have evolved to sustain plant growth in unfavorable environments, especially at low atmospheric carbon levels and high temperatures. Uptake of CO₂ and its storage in the aerenchyma tissues of Lycopsids and diurnal acidity fluctuation in aquatic plants during the Palaeozoic era (ca. 300 Ma.) would represent the earliest evolution of a CCM. The CCM parts of the dark reactions of photosynthesis have evolved many times, while the light reactions are conserved across plant lineages. A C₄ type CCM, leaf C₄ photosynthesis is evolved in the PACMAD clade of the Poaceae family. The evolution of C₄ photosynthesis from C₃ photosynthesis was an abaptation. Photosynthesis in reproductive tissues of sorghum and maize (PACMAD clade) has been shown to be of a weaker C₄ type (high CO₂ compensation point, low carbon isotope discrimination, and lack of Rubisco compartmentalization, when compared to the normal C₄ types) than that in the leaves (normal C₄ type). However, this does not fit well with the character polarity concept from an evolutionary perspective. In a recent model proposed for CCM evolution, the development of a rudimentary CCM prior to the evolution of a more efficient CCM (features contrasting to a weaker C₄ type, leading to greater biomass production rate) has been suggested. An intermediate crassulacean acid metabolism (CAM) type of CCM (rudimentary) was reported in the genera, Brassia, Coryanthes, Eriopsis, Peristeria, of the orchids (well-known group of plants that display the CAM pathway). Similarly, we propose here the evolution of a rudimentary CCM (C₄-like type pathway) in the non-foliar tissues of the Poaceae, prior to the evolution of the C₄ pathway as identified in the leaves of the C₄ species of the PACMAD clade.

A multiple alignment workflow shows the effect of repeat masking and parameter tuning on alignment in plants

Alignments of multiple genomes are a cornerstone of comparative genomics, but generating these alignments remains technically challenging and often impractical. We developed the msa_pipeline workflow (https://bitbucket.org/bucklerlab/msa_pipeline) to allow practical and sensitive multiple alignment of diverged plant genomes and calculation of conservation scores with minimal user inputs. As high repeat content and genomic divergence are substantial challenges in plant genome alignment, we also explored the effect of different masking approaches and parameters of the LAST aligner using genome assemblies of 33 grass species. Compared with conventional masking with RepeatMasker, a masking approach based on k-mers (nucleotide sequences of k length) increased the alignment rate of coding sequence and noncoding functional regions by 25 and 14%, respectively. We further found that default alignment parameters generally perform well, but parameter tuning can increase the alignment rate for noncoding functional regions by over 52% compared with default LAST settings. Finally, by increasing alignment sensitivity from the default baseline, parameter tuning can increase the number of noncoding sites that can be scored for conservation by over 76%. Overall, tuning of masking and alignment parameters can generate optimized multiple alignments to drive biological discovery in plants.

New Insights into the Genomic Structure of Avena L.: Comparison of the Divergence of A-Genome and One C-Genome Oat Species

We used next-generation sequencing analysis of the 3′-part of 18S rDNA, ITS1, and a 5′-part of the 5.8S rDNA region to understand genetic variation among seven diploid A-genome Avena species. We used 4−49 accessions per species that represented the As genome (A. atlantica, A. hirtula, and wiestii), Ac genome (A. canariensis), Ad genome (A. damascena), Al genome (A. longiglumis), and Ap genome (A. prostrata). We also took into our analysis one C-genome species, A. clauda, which previously was found to be related to A-genome species. The sequences of 169 accessions revealed 156 haplotypes of which seven haplotypes were shared by two to five species. We found 16 ribotypes that consisted of a unique sequence with a characteristic pattern of single nucleotide polymorphisms and deletions. The number of ribotypes per species varied from one in A. longiglumis to four in A. wiestii. Although most ribotypes were species-specific, we found two ribotypes shared by three species (one for A. damascena, A. hirtula, and A. wiestii, and the second for A. longiglumis, A. atlantica, and A. wiestii), and a third ribotype shared between A. atlantica and A. wiestii. A characteristic feature of the A. clauda ribotype, a diploid C-genome species, is that two different families of ribotypes have been found in this species. Some of these ribotypes are characteristic of Cc-genome species, whereas others are closely related to As-genome ribotypes. This means that A. clauda can be a hybrid between As- and C-genome oats.

A well-supported nuclear phylogeny of Poaceae and implications for the evolution of C4 photosynthesis

Poaceae (the grasses) includes rice, maize, wheat, and other crops, and is the most economically important angiosperm family. Poaceae is also one of the largest plant families, consisting of over 11 000 species with a global distribution that contributes to diverse ecosystems. Poaceae species are classified into 12 subfamilies, with generally strong phylogenetic support for their monophyly. However, many relationships within subfamilies, among tribes and/or subtribes, remain uncertain. To better resolve the Poaceae phylogeny, we generated 342 transcriptomic and seven genomic datasets; these were combined with other genomic and transcriptomic datasets to provide sequences for 357 Poaceae species in 231 genera, representing 45 tribes and all 12 subfamilies. Over 1200 low-copy nuclear genes were retrieved from these datasets, with several subsets obtained using additional criteria, and used for coalescent analyses to reconstruct a Poaceae phylogeny. Our results strongly support the monophyly of 11 subfamilies; however, the subfamily Puelioideae was separated into two non-sister clades, one for each of the two previously defined tribes, supporting a hypothesis that places each tribe in a separate subfamily. Molecular clock analyses estimated the crown age of Poaceae to be ∼101 million years old. Ancestral character reconstruction of C₃/C₄ photosynthesis supports the hypothesis of multiple independent origins of C₄ photosynthesis. These origins are further supported by phylogenetic analysis of the ppc gene family that encodes the phosphoenolpyruvate carboxylase, which suggests that members of three paralogous subclades (ppc-aL1a, ppc-aL1b, and ppc-B2) were recruited as functional C₄ppc genes. This study provides valuable resources and a robust phylogenetic framework for evolutionary analyses of the grass family.

Phylotranscriptomics Resolves the Phylogeny of Pooideae and Uncovers Factors for Their Adaptive Evolution

Adaptation to cool climates has occurred several times in different angiosperm groups. Among them, Pooideae, the largest grass subfamily with ∼3,900 species including wheat and barley, have successfully occupied many temperate regions and play a prominent role in temperate ecosystems. To investigate possible factors contributing to Pooideae adaptive evolution to cooling climates, we performed phylogenetic reconstruction using five gene sets (with 1,234 nuclear genes and their subsets) from 157 transcriptomes/genomes representing all 15 tribes and 24 of 26 subtribes. Our phylogeny supports the monophyly of all tribes (except Diarrheneae) and all subtribes with at least two species, with strongly supported resolution of their relationships. Molecular dating suggests that Pooideae originated in the late Cretaceous, with subsequent divergences under cooling conditions first among many tribes from the early middle to late Eocene and again among genera in the middle Miocene and later periods. We identified a cluster of gene duplications (CGD5) shared by the core Pooideae (with 80% Pooideae species) near the Eocene–Oligocene transition, coinciding with the transition from closed to open habitat and an upshift of diversification rate. Molecular evolutionary analyses homologs of CBF for cold resistance uncovered tandem duplications during the core Pooideae history, dramatically increasing their copy number and possibly promoting adaptation to cold habitats. Moreover, duplication of AP1/FUL-like genes before the Pooideae origin might have facilitated the regulation of the vernalization pathway under cold environments. These and other results provide new insights into factors that likely have contributed to the successful adaptation of Pooideae members to temperate regions.

The rise of grasslands is linked to atmospheric CO2 decline in the late Palaeogene

Grasslands are predicted to experience a major biodiversity change by the year 2100. A better understanding of how grasslands have responded to past environmental changes will help predict the outcome of current and future environmental changes. Here, we explore the relationship between past atmospheric CO₂ and temperature fluctuations and the shifts in diversification rate of Poaceae (grasses) and Asteraceae (daisies), two exceptionally species-rich grassland families (~11,000 and ~23,000 species, respectively). To this end, we develop a Bayesian approach that simultaneously estimates diversification rates through time from time-calibrated phylogenies and correlations between environmental variables and diversification rates. Additionally, we present a statistical approach that incorporates the information of the distribution of missing species in the phylogeny. We find strong evidence supporting a simultaneous increase in diversification rates for grasses and daisies after the most significant reduction of atmospheric CO₂ in the Cenozoic (~34 Mya). The fluctuations of paleo-temperatures, however, appear not to have had a significant relationship with the diversification of these grassland families. Overall, our results shed new light on our understanding of the origin of grasslands in the context of past environmental changes.

A better understanding of how grasslands have responded to past environmental changes will help predict the outcomes of future changes. This study explores past climatic fluctuations and shifts in the diversification rate of grasses and daisies, finding strong evidence for a simultaneous increase in their diversification rates following a reduction of atmospheric CO₂ in the Cenozoic.

Sample Sequence Analysis Uncovers Recurrent Horizontal Transfers of Transposable Elements among Grasses

Limited genome resources are a bottleneck in the study of horizontal transfer (HT) of DNA in plants. To solve this issue, we tested the usefulness of low-depth sequencing data generated from 19 previously uncharacterized panicoid grasses for HT investigation. We initially searched for horizontally transferred LTR-retrotransposons by comparing the 19 sample sequences to 115 angiosperm genome sequences. Frequent HTs of LTR-retrotransposons were identified solely between panicoids and rice (Oryza sativa). We consequently focused on additional Oryza species and conducted a nontargeted investigation of HT involving the panicoid genus Echinochloa, which showed the most HTs in the first set of analyses. The comparison of nine Echinochloa samples and ten Oryza species identified recurrent HTs of diverse transposable element (TE) types at different points in Oryza history, but no confirmed cases of HT for sequences other than TEs. One case of HT was observed from one Echinochloa species into one Oryza species with overlapping geographic distributions. Variation among species and data sets highlights difficulties in identifying all HT, but our investigations showed that sample sequence analyses can reveal the importance of HT for the diversification of the TE repertoire of plants.

Widespread lateral gene transfer among grasses

Lateral gene transfer (LGT) occurs in a broad range of prokaryotes and eukaryotes, occasionally promoting adaptation. LGT of functional nuclear genes has been reported among some plants, but systematic studies are needed to assess the frequency and facilitators of LGT. We scanned the genomes of a diverse set of 17 grass species that span more than 50 Ma of divergence and include major crops to identify grass-to-grass protein-coding LGT. We identified LGTs in 13 species, with significant variation in the amount each received. Rhizomatous species acquired statistically more genes, probably because this growth habit boosts opportunities for transfer into the germline. In addition, the amount of LGT increases with phylogenetic relatedness, which might reflect genomic compatibility among close relatives facilitating successful transfers. However, genetic exchanges among highly divergent species indicates that transfers can occur across almost the entire family. Overall, we showed that LGT is a widespread phenomenon in grasses that has moved functional genes across the grass family into domesticated and wild species alike. Successful LGTs appear to increase with both opportunity and compatibility.

The Pharus latifolius genome bridges the gap of early grass evolution

The grass family (Poaceae) includes all commercial cereal crops and is a major contributor to biomass in various terrestrial ecosystems. The ancestry of all grass genomes includes a shared whole-genome duplication (WGD), named rho (ρ) WGD, but the evolutionary significance of ρ-WGD remains elusive. We sequenced the genome of Pharus latifolius, a grass species (producing a true spikelet) in the subfamily Pharoideae, a sister lineage to the core Poaceae including the (Panicoideae, Arundinoideae, Chloridoideae, Micrairoideae, Aristidoideae, and Danthonioideae (PACMAD) and Bambusoideae, Oryzoideae, and Pooideae (BOP) clades. Our results indicate that the P. latifolius genome has evolved slowly relative to cereal grass genomes, as reflected by moderate rates of molecular evolution, limited chromosome rearrangements and a low rate of gene loss for duplicated genes. We show that the ρ-WGD event occurred approximately 98.2 million years ago (Ma) in a common ancestor of the Pharoideae and the PACMAD and BOP grasses. This was followed by contrasting patterns of diploidization in the Pharus and core Poaceae lineages. The presence of two FRIZZY PANICLE-like genes in P. latifolius, and duplicated MADS-box genes, support the hypothesis that the ρ-WGD may have played a role in the origin and functional diversification of the spikelet, an adaptation in grasses related directly to cereal yields. The P. latifolius genome sheds light on the origin and early evolution of grasses underpinning the biology and breeding of cereals.

Root Digger: a root placement program for phylogenetic trees

BACKGROUND

In phylogenetic analysis, it is common to infer unrooted trees. However, knowing the root location is desirable for downstream analyses and interpretation. There exist several methods to recover a root, such as molecular clock analysis (including midpoint rooting) or rooting the tree using an outgroup. Non-reversible Markov models can also be used to compute the likelihood of a potential root position.

RESULTS

We present a software called RootDigger which uses a non-reversible Markov model to compute the most likely root location on a given tree and to infer a confidence value for each possible root placement. We find that RootDigger is successful at finding roots when compared to similar tools such as IQ-TREE and MAD, and will occasionally outperform them. Additionally, we find that the exhaustive mode of RootDigger is useful in quantifying and explaining uncertainty in rooting positions.

CONCLUSIONS

RootDigger can be used on an existing phylogeny to find a root, or to asses the uncertainty of the root placement. RootDigger is available under the MIT licence at https://www.github.com/computations/root_digger .

Developmental and biophysical determinants of grass leaf size worldwide

One of the most notable ecological trends-described more than 2,300  years ago by Theophrastus-is the association of small leaves with dry and cold climates, which has recently been recognized for eudicotyledonous plants at a global scale^1-3. For eudicotyledons, this pattern has been attributed to the fact that small leaves have a thinner boundary layer that helps to avoid extreme leaf temperatures⁴ and their leaf development results in vein traits that improve water transport under cold or dry climates^5,6. However, the global distribution of leaf size and its adaptive basis have not been tested in the grasses, which represent a diverse lineage that is distinct in leaf morphology and that contributes 33% of terrestrial primary productivity (including the bulk of crop production)⁷. Here we demonstrate that grasses have shorter and narrower leaves under colder and drier climates worldwide. We show that small grass leaves have thermal advantages and vein development that contrast with those of eudicotyledons, but that also explain the abundance of small leaves in cold and dry climates. The worldwide distribution of leaf size in grasses exemplifies how biophysical and developmental processes result in convergence across major lineages in adaptation to climate globally, and highlights the importance of leaf size and venation architecture for grass performance in past, present and future ecosystems.

A 313 plastome phylogenomic analysis of Pooideae: Exploring relationships among the largest subfamily of grasses

In this study, we analyzed 313 plastid genomes (plastomes) of Poaceae with a focus on expanding our current knowledge of relationships among the subfamily Pooideae, which represented over half the dataset (164 representatives). In total, 47 plastomes were sequenced and assembled for this study. This is the largest study of its kind to include plastome-level data, to not only increase sampling at both the taxonomic and molecular levels with the aim of resolving complex and reticulate relationships, but also to analyze the effects of alignment gaps in large-scale analyses, as well as explore divergences in the subfamily with an expanded set of 14 accepted grass fossils for more accurate calibrations and dating. Incorporating broad systematic assessments of Pooideae taxa conducted by authors within the last five years, we produced a robust phylogenomic reconstruction for the subfamily, which included all but two supergeneric taxa (Calothecinae and Duthieeae). We further explored how including alignment gaps in plastome analyses oftentimes can produce incorrect or misinterpretations of complex or reticulate relationships among taxa of Pooideae. This presented itself as consistently changing relationships at specific nodes for different stripping thresholds (percentage-based removal of gaps per alignment column). Our summary recommendation for large-scale genomic plastome datasets is to strip alignment columns of all gaps to increase pairwise identity and reduce errant signal from poly A/T bias. To do this we used the "mask alignment" tool in Geneious software. Finally, we determined an overall divergence age for Pooideae of roughly 84.8 Mya, which is in line with, but slightly older than most recent estimates.

Continued Adaptation of C4 Photosynthesis After an Initial Burst of Changes in the Andropogoneae Grasses

C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis is a complex trait that sustains fast growth and high productivity in tropical and subtropical conditions and evolved repeatedly in flowering plants. One of the major C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} lineages is Andropogoneae, a group of \documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$\sim $\end{document}1200 grass species that includes some of the world’s most important crops and species dominating tropical and some temperate grasslands. Previous efforts to understand C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} evolution in the group have compared a few model C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} plants to distantly related C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} species so that changes directly responsible for the transition to C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} could not be distinguished from those that preceded or followed it. In this study, we analyze the genomes of 66 grass species, capturing the earliest diversification within Andropogoneae as well as their C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} relatives. Phylogenomics combined with molecular dating and analyses of protein evolution show that many changes linked to the evolution of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} photosynthesis in Andropogoneae happened in the Early Miocene, between 21 and 18 Ma, after the split from its C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{3}$\end{document} sister lineage, and before the diversification of the group. This initial burst of changes was followed by an extended period of modifications to leaf anatomy and biochemistry during the diversification of Andropogoneae, so that a single C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} origin gave birth to a diversity of C\documentclass[12pt]{minimal} \usepackage{amsmath} \usepackage{wasysym} \usepackage{amsfonts} \usepackage{amssymb} \usepackage{amsbsy} \usepackage{upgreek} \usepackage{mathrsfs} \setlength{\oddsidemargin}{-69pt} \begin{document} }{}$_{4}$\end{document} phenotypes during 18 million years of speciation events and migration across geographic and ecological spaces. Our comprehensive approach and broad sampling of the diversity in the group reveals that one key transition can lead to a plethora of phenotypes following sustained adaptation of the ancestral state. [Adaptive evolution; complex traits; herbarium genomics; Jansenelleae; leaf anatomy; Poaceae; phylogenomics.]

The morphogenesis of fast growth in plants

Growth rate represents a fundamental axis of life history variation. Faster growth associated with C₄ photosynthesis and annual life history has evolved multiple times, and the resulting diversity in growth is typically explained via resource acquisition and allocation. However, the underlying changes in morphogenesis remain unknown. We conducted a phylogenetic comparative experiment with 74 grass species, conceptualising morphogenesis as the branching and growth of repeating modules. We aimed to establish whether faster growth in C₄ and annual grasses, compared with C₃ and perennial grasses, came from the faster growth of individual modules or higher rates of module initiation. Morphogenesis produces fast growth in different ways in grasses using C₄ and C₃ photosynthesis, and in annual compared with perennial species. C₄ grasses grow faster than C₃ species through a greater enlargement of shoot modules and quicker secondary branching of roots. However, leaf initiation is slower and there is no change in shoot branching. Conversely, faster growth in annuals than perennials is achieved through greater branching and enlargement of shoots, and possibly faster root branching. The morphogenesis of fast growth depends on ecological context, with C₄ grasses tending to promote resource capture under competition, and annuals enhancing branching to increase reproductive potential.

The Origin of the Legumes is a Complex Paleopolyploid Phylogenomic Tangle Closely Associated with the Cretaceous-Paleogene (K-Pg) Mass Extinction Event

The consequences of the Cretaceous–Paleogene (K–Pg) boundary (KPB) mass extinction for the evolution of plant diversity remain poorly understood, even though evolutionary turnover of plant lineages at the KPB is central to understanding assembly of the Cenozoic biota. The apparent concentration of whole genome duplication (WGD) events around the KPB may have played a role in survival and subsequent diversification of plant lineages. To gain new insights into the origins of Cenozoic biodiversity, we examine the origin and early evolution of the globally diverse legume family (Leguminosae or Fabaceae). Legumes are ecologically (co-)dominant across many vegetation types, and the fossil record suggests that they rose to such prominence after the KPB in parallel with several well-studied animal clades including Placentalia and Neoaves. Furthermore, multiple WGD events are hypothesized to have occurred early in legume evolution. Using a recently inferred phylogenomic framework, we investigate the placement of WGDs during early legume evolution using gene tree reconciliation methods, gene count data and phylogenetic supernetwork reconstruction. Using 20 fossil calibrations we estimate a revised timeline of legume evolution based on 36 nuclear genes selected as informative and evolving in an approximately clock-like fashion. To establish the timing of WGDs we also date duplication nodes in gene trees. Results suggest either a pan-legume WGD event on the stem lineage of the family, or an allopolyploid event involving (some of) the earliest lineages within the crown group, with additional nested WGDs subtending subfamilies Papilionoideae and Detarioideae. Gene tree reconciliation methods that do not account for allopolyploidy may be misleading in inferring an earlier WGD event at the time of divergence of the two parental lineages of the polyploid, suggesting that the allopolyploid scenario is more likely. We show that the crown age of the legumes dates to the Maastrichtian or early Paleocene and that, apart from the Detarioideae WGD, paleopolyploidy occurred close to the KPB. We conclude that the early evolution of the legumes followed a complex history, in which multiple auto- and/or allopolyploidy events coincided with rapid diversification and in association with the mass extinction event at the KPB, ultimately underpinning the evolutionary success of the Leguminosae in the Cenozoic. [Allopolyploidy; Cretaceous–Paleogene (K–Pg) boundary; Fabaceae, Leguminosae; paleopolyploidy; phylogenomics; whole genome duplication events]

Phylogenomics indicates the "living fossil" Isoetes diversified in the Cenozoic

The fossil record provides an invaluable insight into the temporal origins of extant lineages of organisms. However, establishing the relationships between fossils and extant lineages can be difficult in groups with low rates of morphological change over time. Molecular dating can potentially circumvent this issue by allowing distant fossils to act as calibration points, but rate variation across large evolutionary scales can bias such analyses. In this study, we apply multiple dating methods to genome-wide datasets to infer the origin of extant species of Isoetes, a group of mostly aquatic and semi-aquatic isoetalean lycopsids, which closely resemble fossil forms dating back to the Triassic. Rate variation observed in chloroplast genomes hampers accurate dating, but genome-wide nuclear markers place the origin of extant diversity within this group in the mid-Paleogene, 45-60 million years ago. Our genomic analyses coupled with a careful evaluation of the fossil record indicate that despite resembling forms from the Triassic, extant Isoetes species do not represent the remnants of an ancient and widespread group, but instead have spread around the globe in the relatively recent past.

Fire and grazing determined grasslands of central Madagascar represent ancient assemblages

The ecology of Madagascar's grasslands is under-investigated and the dearth of ecological understanding of how disturbance by fire and grazing shapes these grasslands stems from a perception that disturbance shaped Malagasy grasslands only after human arrival. However, worldwide, fire and grazing shape tropical grasslands over ecological and evolutionary timescales, and it is curious Madagascar should be a global anomaly. We examined the functional and community ecology of Madagascar's grasslands across 71 communities in the Central Highlands. Combining multivariate abundance models of community composition and clustering of grass functional traits, we identified distinct grass assemblages each shaped by fire or grazing. The fire-maintained assemblage is primarily composed of tall caespitose species with narrow leaves and low bulk density. By contrast, the grazer-maintained assemblage is characterized by mat-forming, high bulk density grasses with wide leaves. Within each assemblage, levels of endemism, diversity and grass ages support these as ancient assemblages. Grazer-dependent grasses can only have co-evolved with a now-extinct megafauna. Ironically, the human introduction of cattle probably introduced a megafaunal substitute facilitating modern day persistence of a grazer-maintained grass assemblage in an otherwise defaunated landscape, where these landscapes now support the livelihoods of millions of people.

Implications of plastome evolution in the true lilies (monocot order Liliales)

The families of the monocot order Liliales exhibit highly contrasting characteristic of photosynthetic and mycoheterotrophic life histories. Although previous phylogenetic and morphological studies of Liliales have been conducted, they have not examined molecular evolution associated with this contrasting phenomenon. Here, we conduct the first comparative plastome study of all ten families of Liliales using 29 newly sequenced plastid genomes analyzed together with previously published data. We also present a phylogenetic analysis for Liliales of 78 plastid genes combined with 22 genes from all three genomes (nuclear 18S rDNA and phyC; 17 plastid genes; and mitochondrial matR, atpA, and cob). Within the newly generated phylogenetic tree of Liliales, we evaluate the ancestral state changes of selected morphological traits in the order. There are no significant differences in plastid genome features among species that show divergent characteristics correlated with family circumscriptions. However, the results clearly differentiate between photosynthetic and mycoheterotrophic taxa of Liliales in terms of genome structure, and gene content and order. The newly sequenced plastid genomes and combined three-genome data revealed Smilacaceae as sister to Liliaceae instead of Philesiaceae and Ripogonaceae. Additionally, we propose a revised familial classification system of Liliales that consists of nine families, considering Ripogonaceae a synonym of Philesiaceae. The ancestral state reconstruction indicated synapomorphies for each family of Liliales, except Liliaceae, Melanthiaceae and Colchicaceae. A taxonomic key for all nine families of Liliales is also provided.

Genomic insights from the first chromosome-scale assemblies of oat (Avena spp.) diploid species

BACKGROUND

Cultivated hexaploid oat (Common oat; Avena sativa) has held a significant place within the global crop community for centuries; although its cultivation has decreased over the past century, its nutritional benefits have garnered increased interest for human consumption. We report the development of fully annotated, chromosome-scale assemblies for the extant progenitor species of the As- and Cp-subgenomes, Avena atlantica and Avena eriantha respectively. The diploid Avena species serve as important genetic resources for improving common oat's adaptive and food quality characteristics.

RESULTS

The A. atlantica and A. eriantha genome assemblies span 3.69 and 3.78 Gb with an N50 of 513 and 535 Mb, respectively. Annotation of the genomes, using sequenced transcriptomes, identified ~ 50,000 gene models in each species-including 2965 resistance gene analogs across both species. Analysis of these assemblies classified much of each genome as repetitive sequence (~ 83%), including species-specific, centromeric-specific, and telomeric-specific repeats. LTR retrotransposons make up most of the classified elements. Genome-wide syntenic comparisons with other members of the Pooideae revealed orthologous relationships, while comparisons with genetic maps from common oat clarified subgenome origins for each of the 21 hexaploid linkage groups. The utility of the diploid genomes was demonstrated by identifying putative candidate genes for flowering time (HD3A) and crown rust resistance (Pc91). We also investigate the phylogenetic relationships among other A- and C-genome Avena species.

CONCLUSIONS

The genomes we report here are the first chromosome-scale assemblies for the tribe Poeae, subtribe Aveninae. Our analyses provide important insight into the evolution and complexity of common hexaploid oat, including subgenome origin, homoeologous relationships, and major intra- and intergenomic rearrangements. They also provide the annotation framework needed to accelerate gene discovery and plant breeding.

Phylogenetic Shifts in Gene Body Methylation Correlate with Gene Expression and Reflect Trait Conservation

A subset of genes in plant genomes are labeled with DNA methylation specifically at CG residues. These genes, known as gene-body methylated (gbM), have a number of associated characteristics. They tend to have longer sequences, to be enriched for intermediate expression levels, and to be associated with slower rates of molecular evolution. Most importantly, gbM genes tend to maintain their level of DNA methylation between species, suggesting that this trait is under evolutionary constraint. Given the degree of conservation in gbM, we still know surprisingly little about its function in plant genomes or whether gbM is itself a target of selection. To address these questions, we surveyed DNA methylation across eight grass (Poaceae) species that span a gradient of genome sizes. We first established that genome size correlates with genome-wide DNA methylation levels, but less so for genic levels. We then leveraged genomic data to identify a set of 2,982 putative orthologs among the eight species and examined shifts of methylation status for each ortholog in a phylogenetic context. A total of 55% of orthologs exhibited a shift in gbM, but these shifts occurred predominantly on terminal branches, indicating that shifts in gbM are rarely conveyed over time. Finally, we found that the degree of conservation of gbM across species is associated with increased gene length, reduced rates of molecular evolution, and increased gene expression level, but reduced gene expression variation across species. Overall, these observations suggest a basis for evolutionary pressure to maintain gbM status over evolutionary time.

Molecular Dating of the Emergence of Anaerobic Rumen Fungi and the Impact of Laterally Acquired Genes

The anaerobic gut fungi (AGF), or Neocallimastigomycota, inhabit the rumen and alimentary tract of herbivorous mammals, where they play important roles in the degradation of plant fiber. Comparative genomic and phylogenomic analyses of the AGF have long been hampered by their fastidious growth condition, as well as their large (up to 200 Mb) and AT-biased (78 to 84%) genomes. We sequenced 21 AGF transcriptomes and combined them with 5 available AGF genome sequences to explore their evolutionary relationships, time their divergence, and characterize gene gain/loss patterns associated with their evolution. We estimate that the most recent common ancestor of the AGF diverged 66 (±10) million years ago, a time frame that coincides with the evolution of grasses (Poaceae), as well as the mammalian transition from insectivory to herbivory. The concordance of independent estimations suggests that AGF have been important in shaping the success of mammalian herbivory transition by improving the efficiency of energy acquisition from recalcitrant plant materials. Comparative genomics identified multiple lineage-specific genes in the AGF, two of which were acquired from rumen gut bacteria and animal hosts via horizontal gene transfer (HGT). A third AGF domain, plant-like polysaccharide lyase, represents a novel gene in fungi that potentially aids AGF to degrade pectin. Analysis of genomic and transcriptomic sequences confirmed both the presence and expression of these lineage-specific genes in nearly all AGF clades. These genetic elements may contribute to the exceptional abilities of AGF to degrade plant biomass and enable metabolism of the rumen microbes and animal hosts.IMPORTANCE Anaerobic fungi living in the rumen of herbivorous mammals possess an extraordinary ability to degrade plant biomass. We examined the origin and genomic composition of these poorly characterized anaerobic gut fungi using both transcriptome and genomic data. Phylogenomics and molecular dating analyses found remarkable concurrence of the divergence times of the rumen fungi, the forage grasses, and the dietary shift of ancestral mammals from primarily insectivory to herbivory. Comparative genomics identified unique machinery in these fungi to utilize plant polysaccharides. The rumen fungi were also identified with the ability to code for three protein domains with putative functions in plant pectin degradation and microbial defense, which were absent from all other fungal organisms (examined over 1,000 fungal genomes). Two of these domains were likely acquired from rumen gut bacteria and animal hosts separately via horizontal gene transfer. The third one is a plant-like polysaccharide lyase, representing a unique fungal enzyme with potential pectin breakdown abilities.

Key changes in gene expression identified for different stages of C4 evolution in Alloteropsis semialata

C4 photosynthesis is a complex trait that boosts productivity in tropical conditions. Compared with C3 species, the C4 state seems to require numerous novelties, but species comparisons can be confounded by long divergence times. Here, we exploit the photosynthetic diversity that exists within a single species, the grass Alloteropsis semialata, to detect changes in gene expression associated with different photosynthetic phenotypes. Phylogenetically informed comparative transcriptomics show that intermediates with a weak C4 cycle are separated from the C3 phenotype by increases in the expression of 58 genes (0.22% of genes expressed in the leaves), including those encoding just three core C4 enzymes: aspartate aminotransferase, phosphoenolpyruvate carboxykinase, and phosphoenolpyruvate carboxylase. The subsequent transition to full C4 physiology was accompanied by increases in another 15 genes (0.06%), including only the core C4 enzyme pyruvate orthophosphate dikinase. These changes probably created a rudimentary C4 physiology, and isolated populations subsequently improved this emerging C4 physiology, resulting in a patchwork of expression for some C4 accessory genes. Our work shows how C4 assembly in A. semialata happened in incremental steps, each requiring few alterations over the previous step. These create short bridges across adaptive landscapes that probably facilitated the recurrent origins of C4 photosynthesis through a gradual process of evolution.

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

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

The Estimated Pacemaker for Great Apes Supports the Hominoid Slowdown Hypothesis

The recent surge of genomic data has prompted the investigation of substitution rate variation across the genome, as well as among lineages. Evolutionary trees inferred from distinct genomic regions may display branch lengths that differ between loci by simple proportionality constants, indicating that rate variation follows a pacemaker model, which may be attributed to lineage effects. Analyses of genes from diverse biological clades produced contrasting results, supporting either this model or alternative scenarios where multiple pacemakers exist. So far, an evaluation of the pacemaker hypothesis for all great apes has never been carried out. In this work, we tested whether the evolutionary rates of hominids conform to pacemakers, which were inferred accounting for gene tree/species tree discordance. For higher precision, substitution rates in branches were estimated with a calibration-free approach, the relative rate framework. A predominant evolutionary trend in great apes was evidenced by the recovery of a large pacemaker, encompassing most hominid genomic regions. In addition, the majority of genes followed a pace of evolution that was closely related to the strict molecular clock. However, slight rate decreases were recovered in the internal branches leading to humans, corroborating the hominoid slowdown hypothesis. Our findings suggest that in great apes, life history traits were the major drivers of substitution rate variation across the genome.

Recent range expansion in Australian hummock grasses (Triodia) inferred using genotyping-by-sequencing

The Australian arid zone (AAZ) has undergone aridification and the formation of vast sandy deserts since the mid-Miocene. Studies on AAZ organisms, particularly animals, have shown patterns of mesic ancestry, persistence in rocky refugia and range expansions in arid lineages. There has been limited molecular investigation of plants in the AAZ, particularly of taxa that arrived in Australia after the onset of aridification. Here we investigate populations of the widespread AAZ grass Triodia basedowii to determine whether there is evidence for a recent range expansion, and if so, its source and direction. We also undertake a dating analysis for the species complex to which T. basedowii belongs, in order to place its diversification in relation to changes in AAZ climate and landscapes. We analyse a genomic single nucleotide polymorphism data set from 17 populations of T. basedowii in a recently developed approach for detecting the signal and likely origin of a range expansion. We also use alignments from existing and newly sequenced plastomes from across Poaceae for analysis in BEAST to construct fossil-calibrated phylogenies. Across a range of sampling parameters and outgroups, we detected a consistent signal of westward expansion for T. basedowii, originating in central or eastern Australia. Divergence time estimation indicates that Triodia began to diversify in the late Miocene (crown 7.0-8.8 million years (Ma)), and the T. basedowii complex began to radiate during the Pleistocene (crown 1.4-2.0 Ma). This evidence for range expansion in an arid-adapted plant is consistent with similar patterns in AAZ animals and likely reflects a general response to the opening of new habitat during aridification. Radiation of the T. basedowii complex through the Pleistocene has been associated with preferences for different substrates, providing an explanation why only one lineage is widespread across sandy deserts.

Leaf shape and size track habitat transitions across forest-grassland boundaries in the grass family (Poaceae)

Grass leaf shape is a strong indicator of their habitat with linear leaves predominating in open areas and ovate leaves distinguishing forest-associated grasses. This pattern among extant species suggests that ancestral shifts between forest and open habitats may have coincided with changes in leaf shape or size. We tested relationships between habitat, climate, photosynthetic pathway, and leaf shape and size in a phylogenetic framework to evaluate drivers of leaf shape and size variation over the evolutionary history of the family. We also estimated the ancestral habitat of Poaceae and tested whether forest margins served as transitional zones for shifts between forests and grasslands. We found that grass leaf shape is converging toward different shape optima in the forest understory, forest margins, and open habitats. Leaf size also varies with habitat. Grasses have smaller leaves in open and drier areas, and in areas with high solar irradiance. Direct transitions between linear and ovate leaves are rare as are direct shifts between forest and open habitats. The most likely ancestral habitat of the family was the forest understory and forest margins along with an intermediate leaf shape served as important transitional habitat and morphology, respectively, for subsequent shifts across forest-grassland biome boundaries.

Strong selection of the TLR2 coding region among the Lagomorpha suggests an evolutionary history that differs from other mammals

Toll-like receptors (TLRs) are one of the first lines of defense against pathogens and are crucial for triggering an appropriate immune response. Among TLRs, TLR2 is functional in all vertebrates and has high ability in detecting bacterial and viral pathogen ligands. The mammals' phylogenetic tree of TLR2 showed longer branches for the Lagomorpha clade, raising the hypothesis that lagomorphs experienced an acceleration of the mutation rate. This hypothesis was confirmed by (i) Tajima's test of neutrality that revealed different evolutionary rates between lagomorphs and the remaining mammals with lagomorphs presenting higher nucleotide diversity; (ii) genetic distances were similar among lagomorphs and between lagomorphs and other mammals; and (iii) branch models reinforced the existence of an acceleration of the mutation rate in lagomorphs. These results suggest that the lagomorph TLR2 has been strongly involved in pathogen recognition, which probably caused a host-pathogen arms race that led to the observed acceleration of the mutation rate.

Remembering winter through vernalisation

Vernalisation is the programmed physiological process in which prolonged cold-exposure provides competency to flower in plants; widely found in winter and biennial species, such as Arabidopsis, fruit trees, vegetables and wheat. This phenomenon is regulated by diverse genetic networks, and memory of vernalisation in a life cycle mainly depends on epigenetic mechanisms. However, less is known about how to count winter-dosage for flowering in plants. Here, we compare the vernalisation genetic framework between the dicots Arabidopsis, temperate grasses, wheat, barley and Brachypodium. We discuss vernalisation mechanisms involving crosstalk between phosphorylation and O-GlcNAcylation modification of key proteins, and epigenetic modifications of the key gene VRN1 in wheat. We also highlight the potential evolutionary origins of vernalisation in various species. Current progress toward understanding the regulation of vernalisation requirements provides insight that will inform the design of molecular breeding strategies for winter crops.

Opposite macroevolutionary responses to environmental changes in grasses and insects during the Neogene grassland expansion

The rise of Neogene C4 grasslands is one of the most drastic changes recently experienced by the biosphere. A central - and widely debated - hypothesis posits that Neogene grasslands acted as a major adaptive zone for herbivore lineages. We test this hypothesis with a novel model system, the Sesamiina stemborer moths and their associated host-grasses. Using a comparative phylogenetic framework integrating paleoenvironmental proxies we recover a negative correlation between the evolutionary trajectories of insects and plants. Our results show that paleoenvironmental changes generated opposing macroevolutionary dynamics in this insect-plant system and call into question the role of grasslands as a universal adaptive cradle. This study illustrates the importance of implementing environmental proxies in diversification analyses to disentangle the relative impacts of biotic and abiotic drivers of macroevolutionary dynamics.

Whole-Genome Duplication and Plant Macroevolution

Whole-genome duplication (WGD) is characteristic of almost all fundamental lineages of land plants. Unfortunately, the timings of WGD events are loosely constrained and hypotheses of evolutionary consequence are poorly formulated, making them difficult to test. Using examples from across the plant kingdom, we show that estimates of timing can be improved through the application of molecular clock methodology to multigene datasets. Further, we show that phenotypic change can be quantified in morphospaces and that relative phenotypic disparity can be compared in the light of WGD. Together, these approaches facilitate tests of hypotheses on the role of WGD in plant evolution, underscoring the potential of plants as a model system for investigating the role WGD in macroevolution.

Shifting the limits in wheat research and breeding using a fully annotated reference genome

An annotated reference sequence representing the hexaploid bread wheat genome in 21 pseudomolecules has been analyzed to identify the distribution and genomic context of coding and noncoding elements across the A, B, and D subgenomes. With an estimated coverage of 94% of the genome and containing 107,891 high-confidence gene models, this assembly enabled the discovery of tissue- and developmental stage-related coexpression networks by providing a transcriptome atlas representing major stages of wheat development. Dynamics of complex gene families involved in environmental adaptation and end-use quality were revealed at subgenome resolution and contextualized to known agronomic single-gene or quantitative trait loci. This community resource establishes the foundation for accelerating wheat research and application through improved understanding of wheat biology and genomics-assisted breeding.

Investigation of mitochondrial-derived plastome sequences in the Paspalum lineage (Panicoideae; Poaceae)

BACKGROUND

The grass family (Poaceae), ca. 12,075 species, is a focal point of many recent studies that aim to use complete plastomes to reveal and strengthen relationships within the family. The use of Next Generation Sequencing technology has revealed intricate details in many Poaceae plastomes; specifically the trnI - trnL intergenic spacer region. This study investigates this region and the putative mitochondrial inserts within it in complete plastomes of Paspalum and other Poaceae.

RESULTS

Nine newly sequenced plastomes, seven of which contain an insert within the trnI - trnL intergenic spacer, were combined into plastome phylogenomic and divergence date analyses with 52 other species. A robust Paspalum topology was recovered, originating at 10.6 Ma, with the insert arising at 8.7 Ma. The alignment of the insert across Paspalum reveals 21 subregions with pairwise homology in 19. In an analysis of emergent self-organizing maps of tetranucleotide frequencies, the Paspalum insert grouped with mitochondrial DNA.

CONCLUSIONS

A hypothetical ancestral insert, 17,685 bp in size, was found in the trnI - trnL intergenic spacer for the Paspalum lineage. A different insert, 2808 bp, was found in the same region for Paraneurachne muelleri. Seven different intrastrand deletion events were found within the Paspalum lineage, suggesting selective pressures to remove large portions of noncoding DNA. Finally, a tetranucleotide frequency analysis was used to determine that the origin of the insert in the Paspalum lineage is mitochondrial DNA.

Some like it hot: the physiological ecology of C4 plant evolution

The evolution of C₄ photosynthesis requires an intermediate phase where photorespiratory glycine produced in the mesophyll cells must flow to the vascular sheath cells for metabolism by glycine decarboxylase. This glycine flux concentrates photorespired CO₂ within the sheath cells, allowing it to be efficiently refixed by sheath Rubisco. A modest C₄ biochemical cycle is then upregulated, possibly to support the refixation of photorespired ammonia in sheath cells, with subsequent increases in C₄ metabolism providing incremental benefits until an optimized C₄ pathway is established. 'Why' C₄ photosynthesis evolved is largely explained by ancestral C₃ species exploiting photorespiratory CO₂ to improve carbon gain and thus enhance fitness. While photorespiration depresses C₃ performance, it produces a resource (photorespired CO₂) that can be exploited to build an evolutionary bridge to C₄ photosynthesis. 'Where' C₄ evolved is indicated by the habitat of species branching near C₃-to-C₄ transitions on phylogenetic trees. Consistent with the photorespiratory bridge hypothesis, transitional species show that the large majority of > 60 C₄ lineages arose in hot, dry, and/or saline regions where photorespiratory potential is high. 'When' C₄ evolved has been clarified by molecular clock analyses using phylogenetic data, coupled with isotopic signatures from fossils. Nearly all C₄ lineages arose after 25 Ma when atmospheric CO₂ levels had fallen to near current values. This reduction in CO₂, coupled with persistent high temperature at low-to-mid-latitudes, met a precondition where photorespiration was elevated, thus facilitating the evolutionary selection pressure that led to C₄ photosynthesis.

Comparative plastome genomics and phylogenomics of Brachypodium: flowering time signatures, introgression and recombination in recently diverged ecotypes

Few pan-genomic studies have been conducted in plants, and none of them have focused on the intraspecific diversity and evolution of their plastid genomes. We address this issue in Brachypodium distachyon and its close relatives B. stacei and B. hybridum, for which a large genomic data set has been compiled. We analyze inter- and intraspecific plastid comparative genomics and phylogenomic relationships within a family-wide framework. Major indel differences were detected between Brachypodium plastomes. Within B. distachyon, we detected two main lineages, a mostly Extremely Delayed Flowering (EDF+) clade and a mostly Spanish (S+) - Turkish (T+) clade, plus nine chloroplast capture and two plastid DNA (ptDNA) introgression and micro-recombination events. Early Oligocene (30.9 million yr ago (Ma)) and Late Miocene (10.1 Ma) divergence times were inferred for the respective stem and crown nodes of Brachypodium and a very recent Mid-Pleistocene (0.9 Ma) time for the B. distachyon split. Flowering time variation is a main factor driving rapid intraspecific divergence in B. distachyon, although it is counterbalanced by repeated introgression between previously isolated lineages. Swapping of plastomes between the three different genomic groups, EDF+, T+, S+, probably resulted from random backcrossing followed by stabilization through selection pressure.

CO2 availability influences hydraulic function of C3 and C4 grass leaves

Atmospheric CO2 (ca) has increased since the last glacial period, increasing photosynthetic water use efficiency and improving plant productivity. Evolution of C4 photosynthesis at low ca led to decreased stomatal conductance (gs), which provided an advantage over C3 plants that may be reduced by rising ca. Using controlled environments, we determined how increasing ca affects C4 water use relative to C3 plants. Leaf gas exchange and mass per area (LMA) were measured for four C3 and four C4 annual, crop-related grasses at glacial (200 µmol mol-1), ambient (400 µmol mol-1), and super-ambient (640 µmol mol-1) ca. C4 plants had lower gs, which resulted in a water use efficiency advantage at all ca and was broadly consistent with slower stomatal responses to shade, indicating less pressure on leaf water status. At glacial ca, net CO2 assimilation and LMA were lower for C3 than for C4 leaves, and C3 and C4 grasses decreased leaf hydraulic conductance (Kleaf) similarly, but only C4 leaves decreased osmotic potential at turgor loss. Greater carbon availability in C4 leaves at glacial ca generated a different hydraulic adjustment relative to C3 plants. At current and future ca, C4 grasses have advantages over C3 grasses due to lower gs, lower stomatal sensitivity, and higher absolute water use efficiency.

Enhancing (crop) plant photosynthesis by introducing novel genetic diversity

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

Molecular phylogeny of Panicum s. str. (Poaceae, Panicoideae, Paniceae) and insights into its biogeography and evolution

Panicum sensu stricto is a genus of grasses (Poaceae) with nearly, according to this study, 163 species distributed worldwide. This genus is included in the subtribe Panicinae together with Louisiella, the latter with 2 species. Panicum and subtribe Panicinae are characterized by including annual or perennial taxa with open and lax panicles, and spikelets with the lower glume reduced; all taxa also share a basic chromosome number of x = 9 and a Kranz leaf blade anatomy typical of the NAD-me subtype photosynthetic pathway. Nevertheless, the phylogenetic placements of many Panicum species, and the circumscription of the genus, remained untested. Therefore, phylogenetic analyses were conducted using sequence data from the ndhF plastid region, in an extensive worldwide sampling of Panicum and related genera, in order to infer evolutionary relationships and to provide a phylogenetic framework to review the classification of the genus. Diversification times, historical biogeography and evolutionary patterns of the life history (annual vs. perennial) in the subtribe and Panicum were also studied. Results obtained provide strong support for a monophyletic Panicum including 71 species and 7 sections, of which sections Arthragrostis and Yakirra are new in the genus; 7 new combinations are made here. Furthermore, 32 species traditionally assigned to Panicum were excluded from the genus, and discussed in other subtribes of Paniceae. Our study suggested that early diversification in subtribe Panicinae and Panicum occurred through the Early-Mid Miocene in the Neotropics, while the subsequent diversification of its sections mainly occurred in the Late Miocene-Pleistocene, involving multiple dispersals to all continents. Our analyses also showed that transition rates and changes between annual and perennial life history in Panicum were quite frequent, suggesting considerable lability of this trait. Changes of the life history, together with C4 photosynthesis, and the multiple dispersal events since the Mid Miocene, seem to have facilitated a widespread distribution of the genus. All these findings contribute to a better understanding of the systematics and evolution of Panicum.

The rate and potential relevance of new mutations in a colonizing plant lineage

By following the evolution of populations that are initially genetically homogeneous, much can be learned about core biological principles. For example, it allows for detailed studies of the rate of emergence of de novo mutations and their change in frequency due to drift and selection. Unfortunately, in multicellular organisms with generation times of months or years, it is difficult to set up and carry out such experiments over many generations. An alternative is provided by "natural evolution experiments" that started from colonizations or invasions of new habitats by selfing lineages. With limited or missing gene flow from other lineages, new mutations and their effects can be easily detected. North America has been colonized in historic times by the plant Arabidopsis thaliana, and although multiple intercrossing lineages are found today, many of the individuals belong to a single lineage, HPG1. To determine in this lineage the rate of substitutions-the subset of mutations that survived natural selection and drift-, we have sequenced genomes from plants collected between 1863 and 2006. We identified 73 modern and 27 herbarium specimens that belonged to HPG1. Using the estimated substitution rate, we infer that the last common HPG1 ancestor lived in the early 17th century, when it was most likely introduced by chance from Europe. Mutations in coding regions are depleted in frequency compared to those in other portions of the genome, consistent with purifying selection. Nevertheless, a handful of mutations is found at high frequency in present-day populations. We link these to detectable phenotypic variance in traits of known ecological importance, life history and growth, which could reflect their adaptive value. Our work showcases how, by applying genomics methods to a combination of modern and historic samples from colonizing lineages, we can directly study new mutations and their potential evolutionary relevance.

Global grass (Poaceae) success underpinned by traits facilitating colonization, persistence and habitat transformation

Poaceae (the grasses) is arguably the most successful plant family, in terms of its global occurrence in (almost) all ecosystems with angiosperms, its ecological dominance in many ecosystems, and high species richness. We suggest that the success of grasses is best understood in context of their capacity to colonize, persist, and transform environments (the "Viking syndrome"). This results from combining effective long-distance dispersal, efficacious establishment biology, ecological flexibility, resilience to disturbance and the capacity to modify environments by changing the nature of fire and mammalian herbivory. We identify a diverse set of functional traits linked to dispersal, establishment and competitive abilities. Enhanced long-distance dispersal is determined by anemochory, epizoochory and endozoochory and is facilitated via the spikelet (and especially the awned lemma) which functions as the dispersal unit. Establishment success could be a consequence of the precocious embryo and large starch reserves, which may underpin the extremely short generation times in grasses. Post-establishment genetic bottlenecks may be mitigated by wind pollination and the widespread occurrence of polyploidy, in combination with gametic self-incompatibility. The ecological competitiveness of grasses is corroborated by their dominance across the range of environmental extremes tolerated by angiosperms, facilitated by both C₃ and C₄ photosynthesis, well-developed frost tolerance in several clades, and a sympodial growth form that enabled the evolution of both annual and long-lived life forms. Finally, absence of investment in wood (except in bamboos), and the presence of persistent buds at or below ground level, provides tolerance of repeated defoliation (whether by fire, frost, drought or herbivores). Biotic modification of environments via feedbacks with herbivory or fire reinforce grass dominance leading to open ecosystems. Grasses can be both palatable and productive, fostering high biomass and diversity of mammalian herbivores. Many grasses have a suite of architectural and functional traits that facilitate frequent fire, including a tufted growth form, and tannin-like substances in leaves which slow decomposition. We mapped these traits over the phylogeny of the Poales, spanning the grasses and their relatives, and demonstrated the accumulation of traits since monocots originated in the mid-Cretaceous. Although the sympodial growth form is a monocot trait, tillering resulting in the tufted growth form most likely evolved within the grasses. Similarly, although an ovary apparently constructed of a single carpel evolved in the most recent grass ancestor, spikelets and the awned lemma dispersal units evolved within the grasses. Frost tolerance and C₄ photosynthesis evolved relatively late (late Palaeogene), and the last significant trait to evolve was probably the production of tannins, associated with pyrophytic savannas. This fits palaeobotanical data, suggesting several phases in the grass success story: from a late Cretaceous origin, to occasional tropical grassland patches in the later Palaeogene, to extensive C₃ grassy woodlands in the early-middle Miocene, to the dramatic expansion of the tropical C₄ grass savannas and grasslands in the Pliocene, and the C₃ steppe grasslands during the Pleistocene glacial periods. Modern grasslands depend heavily on strongly seasonal climates, making them sensitive to climate change.

The spatial distribution of African savannah herbivores: species associations and habitat occupancy in a landscape context

Herbivores play an important role in determining the structure and function of tropical savannahs. Here, we (i) outline a framework for how interactions among large mammalian herbivores, carnivores and environmental variation influence herbivore habitat occupancy in tropical savannahs. We then (ii) use a Bayesian hierarchical model to analyse camera trap data to quantify spatial patterns of habitat occupancy for lions and eight common ungulates of varying body size across an approximately 1100 km(2) landscape in the Serengeti ecosystem. Our results reveal strong positive associations among herbivores at the scale of the entire landscape. Lions were positively associated with migratory ungulates but negatively associated with residents. Herbivore habitat occupancy differed with body size and migratory strategy: large-bodied migrants, at less risk of predation and able to tolerate lower quality food, were associated with high NDVI, while smaller residents, constrained to higher quality forage, avoided these areas. Small herbivores were strongly associated with fires, likely due to the subsequent high-quality regrowth, while larger herbivores avoided burned areas. Body mass was strongly related to herbivore habitat use, with larger species more strongly associated with riverine and woodlands than smaller species. Large-bodied migrants displayed diffuse habitat occupancy, whereas smaller species demonstrated fine-scale occupancy reflecting use of smaller patches of high-quality habitat. Our results demonstrate the emergence of strong positive spatial associations among a diverse group of savannah herbivores, while highlighting species-specific habitat selection strongly determined by herbivore body size.This article is part of the themed issue 'Tropical grassy biomes: linking ecology, human use and conservation'.

Polyphyly of Arundinoideae (Poaceae) and evolution of the twisted geniculate lemma awn

BACKGROUND AND AIMS

Subfamily Arundinoideae represents one of the last unsolved taxonomic mysteries in the grass family (Poaceae) due to the narrow and remote distributions of many of its 19 morphologically and ecologically heterogeneous genera. Resolving the phylogenetic relationships of these genera could have substantial implications for understanding character evolution in the grasses, for example the twisted geniculate awn - a hygroscopic awn that has been shown to be important in seed germination for some grass species. In this study, the phylogenetic positions of most arundinoid genera were determined using DNA from herbarium specimens, and their placement affects interpretation of this ecologically important trait.

METHODS

A phylogenetic analysis was conducted on a matrix of full-plastome sequences from 123 species in 107 genera representing all grass subfamilies, with 15 of the 19 genera in subfamily Arundinoideae. Parsimony and maximum likelihood mapping approaches were used to estimate ancestral states for presence of a geniculate lemma awn with a twisted column across Poaceae. Lastly, anatomical characters were examined for former arundinoid taxa using light microscopy and scanning electron microscopy.

KEY RESULTS

Four genera traditionally included in Arundinoideae fell outside the subfamily in the plastome phylogeny, with the remaining 11 genera forming Arundinoideae sensu stricto . The twisted geniculate awn has originated independently at least five times in the PACMAD grasses, in the subfamilies Panicoideae, Danthonioideae/Chloridoideae and Arundinoideae. Morphological and anatomical characters support the new positions of the misplaced arundinoid genera in the phylogeny, but also highlight convergent and parallel evolution in the grasses.

CONCLUSIONS

In placing the majority of arundinoid genera in a phylogenetic framework, our study answers one of the last remaining big questions in grass taxonomy while highlighting examples of convergent evolution in an ecologically important trait, the hygroscopic, twisted geniculate awn.