ANGIOSPERM DIVERGENCE TIMES: THE EFFECT OF GENES, CODON POSITIONS, AND TIME CONSTRAINTS

The dynamic history of gymnosperm plastomes: Insights from structural characterization, comparative analysis, phylogenomics, and time divergence

Gymnosperms are among the most endangered groups of plant species; they include ginkgo, pines (Conifers I), cupressophytes (Conifers II), cycads, and gnetophytes. The relationships among the five extant gymnosperm groups remain equivocal. We analyzed 167 available gymnosperm plastomes and investigated their diversity and phylogeny. We found that plastome size, structure, and gene order were highly variable in the five gymnosperm groups, of which Parasitaxus usta (Vieill.) de Laub. and Macrozamia mountperriensis F.M.Bailey had the smallest and largest plastomes, respectively. The inverted repeats (IRs) of the five groups were shown to have evolved through distinctive evolutionary scenarios. The IRs have been lost in all conifers but retained in cycads and gnetophytes. A positive association between simple sequence repeat (SSR) abundance and plastome size was observed, and the SSRs with the most variation were found in Pinaceae. Furthermore, the number of repeats was negatively correlated with IR length; thus, the highest number of repeats was detected in Conifers I and II, in which the IRs had been lost. We constructed a phylogeny based on 29 shared genes from 167 plastomes. With the plastome tree and 13 calibrations, we estimated the tree height between present-day angiosperms and gymnosperms to be ∼380 million years ago (mya). The placement of Gnetales in the tree agreed with the Gnetales-other gymnosperms hypothesis. The divergence between Ginkgo and cycads was estimated as ∼284 mya; the crown age of the cycads was 251 mya. Our time-calibrated plastid-based phylogenomic tree provides a framework for comparative studies of gymnosperm evolution.

The implications of interrelated assumptions on estimates of divergence times and rates of diversification

Phylogenies are increasingly being used as a basis to provide insight into macroevolutionary history. Here, we use simulation experiments and empirical analyses to evaluate methods that use phylogenies as a basis to make estimates of divergence times and rates of diversification. This is the first study to present a comprehensive assessment of the key variables that underpin analyses in this field - including substitution rates, speciation rates, and extinction, plus character sampling and taxon sampling. We show that in unrealistically simplistic cases (where substitution rates and speciation rates are constant, and where there is no extinction), increased character and taxon sampling lead to more accurate and precise parameter estimates. By contrast, in more complex but realistic cases (where substitution rates, speciation rates, and extinction rates vary), gains in accuracy and precision from increased character and taxon sampling are far more limited. The lack of accuracy and precision even occurs when using methods that are designed to account for more complex cases, such as relaxed clocks, fossil calibrations, and models that allow speciation rates and extinction rates to vary. The problem also persists when analysing genomic scale datasets. These results suggest two interrelated problems that occur when the processes that generated the data are more complex. First, methodological assumptions are more likely to be violated. Second, limitations in the information content of the data become more important.

Evolutionary network genomics of wood formation in a phylogenetic survey of angiosperm forest trees

Wood formation was present in early angiosperms, but has been highly modified through evolution to generate the anatomical diversity seen in extant angiosperm lineages. In this project, we modeled changes in gene coexpression relationships associated with the evolution of wood formation in a phylogenetic survey of 13 angiosperm tree species. Gravitropic stimulation was used as an experimental treatment to alter wood formation and also perturb gene expression. Gene transcript abundances were determined using RNA sequencing of developing wood tissues from upright trees, and from the top (tension wood) and bottom (opposite wood) tissues of gravistimulated trees. A network-based approach was employed to align gene coexpression networks across species based on orthologous relationships. A large-scale, multilayer network was modeled that identified both lineage-specific gene coexpression modules and modules conserved across multiple species. Functional annotation and analysis of modules identified specific regulatory processes associated with conserved modules, including regulation of hormones, protein phosphorylation, meristem development and epigenetic processes. Our results provide novel insights into the evolution and development of wood formation, and demonstrate the ability to identify biological processes and genes important for the evolution of a foundational trait in nonmodel, undomesticated forest trees.

Exceptionally high rates of positive selection on the rbcL gene in the genus Ilex (Aquifoliaceae)

Background

The genus Ilex (Aquifoliaceae) has a near-cosmopolitan distribution in mesic habitats from tropical to temperate lowlands and in alpine forests. It has a high rate of hybridization and plastid capture, and comprises four geographically structured plastid groups. A previous study showed that the plastid rbcL gene, coding for the large subunit of Rubisco, has a particularly high rate of non-synonymous substitutions in Ilex, when compared with other plant lineages. This suggests a strong positive selection on rbcL, involved in yet unknown adaptations. We therefore investigated positive selection on rbcL in 240 Ilex sequences from across the global range.

Results

The rbcL gene shows a much higher rate of positive selection in Ilex than in any other plant lineage studied so far (> 3000 species) by tests in both PAML and SLR. Most positively selected residues are on the surface of the folded large subunit, suggesting interaction with other subunits and associated chaperones, and coevolution between positively selected residues is prevalent, indicating compensatory mutations to recover molecular stability. Coevolution between positively selected sites to restore global stability is common.

Conclusions

This study has confirmed the predicted high incidence of positively selected residues in rbcL in Ilex, and shown that this is higher than in any other plant lineage studied so far. The causes and consequences of this high incidence are unclear, but it is probably associated with the similarly high incidence of hybridization and introgression in Ilex, even between distantly related lineages, resulting in large cytonuclear discordance in the phylogenies.

Evaluating the Impact of Genomic Data and Priors on Bayesian Estimates of the Angiosperm Evolutionary Timescale

The evolutionary timescale of angiosperms has long been a key question in biology. Molecular estimates of this timescale have shown considerable variation, being influenced by differences in taxon sampling, gene sampling, fossil calibrations, evolutionary models, and choices of priors. Here, we analyze a data set comprising 76 protein-coding genes from the chloroplast genomes of 195 taxa spanning 86 families, including novel genome sequences for 11 taxa, to evaluate the impact of models, priors, and gene sampling on Bayesian estimates of the angiosperm evolutionary timescale. Using a Bayesian relaxed molecular-clock method, with a core set of 35 minimum and two maximum fossil constraints, we estimated that crown angiosperms arose 221 (251-192) Ma during the Triassic. Based on a range of additional sensitivity and subsampling analyses, we found that our date estimates were generally robust to large changes in the parameters of the birth-death tree prior and of the model of rate variation across branches. We found an exception to this when we implemented fossil calibrations in the form of highly informative gamma priors rather than as uniform priors on node ages. Under all other calibration schemes, including trials of seven maximum age constraints, we consistently found that the earliest divergences of angiosperm clades substantially predate the oldest fossils that can be assigned unequivocally to their crown group. Overall, our results and experiments with genome-scale data suggest that reliable estimates of the angiosperm crown age will require increased taxon sampling, significant methodological changes, and new information from the fossil record. [Angiospermae, chloroplast, genome, molecular dating, Triassic.].

Extremely thick cell walls and low mesophyll conductance: welcome to the world of ancient living!

Mesophyll conductance is thought to be an important photosynthetic limitation in gymnosperms, but they currently constitute the most understudied plant group in regard to the extent to which photosynthesis and intrinsic water use efficiency are limited by mesophyll conductance. A comprehensive analysis of leaf gas exchange, photosynthetic limitations, mesophyll conductance (calculated by three methods previously used for across-species comparisons), and the underlying ultra-anatomical, morphological and chemical traits in 11 gymnosperm species varying in evolutionary history was performed to gain insight into the evolution of structural and physiological controls on photosynthesis at the lower return end of the leaf economics spectrum. Two primitive herbaceous species were included in order to provide greater evolutionary context. Low mesophyll conductance was the main limiting factor of photosynthesis in the majority of species. The strongest sources of limitation were extremely thick mesophyll cell walls, high chloroplast thickness and variation in chloroplast shape and size, and the low exposed surface area of chloroplasts per unit leaf area. In gymnosperms, the negative relationship between net assimilation per mass and leaf mass per area reflected an increased mesophyll cell wall thickness, whereas the easy-to-measure integrative trait of leaf mass per area failed to predict the underlying ultrastructural traits limiting mesophyll conductance.

A metacalibrated time-tree documents the early rise of flowering plant phylogenetic diversity

The establishment of modern terrestrial life is indissociable from angiosperm evolution. While available molecular clock estimates of angiosperm age range from the Paleozoic to the Late Cretaceous, the fossil record is consistent with angiosperm diversification in the Early Cretaceous. The time-frame of angiosperm evolution is here estimated using a sample representing 87% of families and sequences of five plastid and nuclear markers, implementing penalized likelihood and Bayesian relaxed clocks. A literature-based review of the palaeontological record yielded calibrations for 137 phylogenetic nodes. The angiosperm crown age was bound within a confidence interval calculated with a method that considers the fossil record of the group. An Early Cretaceous crown angiosperm age was estimated with high confidence. Magnoliidae, Monocotyledoneae and Eudicotyledoneae diversified synchronously 135-130 million yr ago (Ma); Pentapetalae is 126-121 Ma; and Rosidae (123-115 Ma) preceded Asteridae (119-110 Ma). Family stem ages are continuously distributed between c. 140 and 20 Ma. This time-frame documents an early phylogenetic proliferation that led to the establishment of major angiosperm lineages, and the origin of over half of extant families, in the Cretaceous. While substantial amounts of angiosperm morphological and functional diversity have deep evolutionary roots, extant species richness was probably acquired later.

Leading dimensions in absorptive root trait variation across 96 subtropical forest species

Absorptive root traits show remarkable cross-species variation, but major root trait dimensions across species have not been defined. We sampled first-order roots and measured 14 root traits for 96 angiosperm woody species from subtropical China, including root diameter, specific root length, stele diameter, cortex thickness, root vessel size and density, mycorrhizal colonization rate, root branching intensity, tissue density, and concentrations of carbon and nitrogen ([N]). Root traits differed in the degree of variation and phylogenetic conservatism, but showed predictable patterns of cross-trait coordination. Root diameter, cortex thickness and stele diameter displayed high variation across species (coefficient of variation (CV)=0.51-0.69), whereas the stele:root diameter ratio and [N] showed low variation (CV<0.32). Root diameter, cortex thickness and stele diameter showed a strong phylogenetic signal across species, whereas root branching traits did not, and these two sets of traits were segregated onto two nearly orthogonal (independent) principal component analysis (PCA) axes. Two major dimensions of root trait variation were found: a diameter-related dimension potentially integrating root construction, maintenance, and persistence with mycorrhizal colonization, and a branching architecture dimension expressing root plastic responses to the environment. These two dimensions may offer a promising path for better understanding root trait economics and root ecological strategies world-wide.

DNA barcode authentication of saw palmetto herbal dietary supplements

Herbal dietary supplements made from saw palmetto (Serenoa repens; Arecaceae) fruit are commonly consumed to ameliorate benign prostate hyperplasia. A novel DNA mini-barcode assay to accurately identify [specificity = 1.00 (95% confidence interval = 0.74-1.00); sensitivity = 1.00 (95% confidence interval = 0.66-1.00); n = 31] saw palmetto dietary supplements was designed from a DNA barcode reference library created for this purpose. The mini-barcodes were used to estimate the frequency of mislabeled saw palmetto herbal dietary supplements on the market in the United States of America. Of the 37 supplements examined, amplifiable DNA could be extracted from 34 (92%). Mini-barcode analysis of these supplements demonstrated that 29 (85%) contain saw palmetto and that 2 (6%) supplements contain related species that cannot be legally sold as herbal dietary supplements in the United States of America. The identity of 3 (9%) supplements could not be conclusively determined.

A new arboreal haramiyid shows the diversity of crown mammals in the Jurassic period

A major unsolved problem in mammalian evolution is the origin of Allotheria, including Multituberculata and Haramiyida. Multituberculates are the most diverse and best known Mesozoic era mammals and ecologically resemble rodents, but haramiyids are known mainly from isolated teeth, hampering our search for their phylogenetic relationships. Here we report a new haramiyid from the Jurassic period of China, which is, to our knowledge the largest reported so far. It has a novel dentition, a mandible resembling advanced multituberculates and postcranial features adapted for arboreal life. Our phylogenetic analysis places Haramiyida within crown Mammalia, suggesting the origin of crown Mammalia in the Late Triassic period and diversification in the Jurassic, which contrasts other estimated divergence times of crown Mammalia. The new haramiyid reveals additional mammalian features of the group, helps to identify other haramiyids represented by isolated teeth, and shows again that, regardless of various phylogenetic scenarios, a complex pattern of evolution involving many convergences and/or reversals existed in Mesozoic mammals.

Establishing a time-scale for plant evolution

• Plants have utterly transformed the planet, but testing hypotheses of causality requires a reliable time-scale for plant evolution. While clock methods have been extensively developed, less attention has been paid to the correct interpretation and appropriate implementation of fossil data. • We constructed 17 calibrations, consisting of minimum constraints and soft maximum constraints, for divergences between model representatives of the major land plant lineages. Using a data set of seven plastid genes, we performed a cross-validation analysis to determine the consistency of the calibrations. Six molecular clock analyses were then conducted, one with the original calibrations, and others exploring the impact on divergence estimates of changing maxima at basal nodes, and prior probability densities within calibrations. • Cross-validation highlighted Tracheophyta and Euphyllophyta calibrations as inconsistent, either because their soft maxima were overly conservative or because of undetected rate variation. Molecular clock analyses yielded estimates ranging from 568-815 million yr before present (Ma) for crown embryophytes and from 175-240 Ma for crown angiosperms. • We reject both a post-Jurassic origin of angiosperms and a post-Cambrian origin of land plants. Our analyses also suggest that the establishment of the major embryophyte lineages occurred at a much slower tempo than suggested in most previous studies. These conclusions are entirely compatible with current palaeobotanical data, although not necessarily with their interpretation by palaeobotanists.

In the Light of Evolution: A Reevaluation of Conservation in the CO-FT Regulon and Its Role in Photoperiodic Regulation of Flowering Time

In order to maximize reproductive success, plants have evolved different strategies to control the critical developmental shift marked by the transition to flowering. As plants have adapted to diverse environments across the globe, these strategies have evolved to recognize and respond to local seasonal cues through the induction of specific downstream genetic pathways, thereby ensuring that the floral transition occurs in favorable conditions. Determining the genetic factors involved in controlling the floral transition in many species is key to understanding how this trait has evolved. Striking genetic discoveries in Arabidopsis thaliana (Arabidopsis) and Oryza sativa (rice) revealed that similar genes in both species control flowering in response to photoperiod, suggesting that this genetic module could be conserved between distantly related angiosperms. However, as we have gained a better understanding of the complex evolution of these genes and their functions in other species, another possibility must be considered: that the genetic module controlling flowering in response to photoperiod is the result of convergence rather than conservation. In this review, we show that while data clearly support a central role of FLOWERING LOCUS T (FT) homologs in floral promotion across a diverse group of angiosperms, there is little evidence for a conserved role of CONSTANS (CO) homologs in the regulation of these loci. In addition, although there is an element of conserved function for FT homologs, even this component has surprising complexity in its regulation and evolution.

Molecular phylogenetics, temporal diversification, and principles of evolution in the mustard family (Brassicaceae)

Brassicaceae is an important family at both the agronomic and scientific level. The family not only includes several model species, but it is also becoming an evolutionary model at the family level. However, resolving the phylogenetic relationships within the family has been problematic, and a large-scale molecular phylogeny in terms of generic sampling and number of genes is still lacking. In particular, the deeper relationships within the family, for example between the three major recognized lineages, prove particularly hard to resolve. Using a slow-evolving mitochondrial marker (nad4 intron 1), we reconstructed a comprehensive phylogeny in generic representation for the family. In addition, and because resolution was very low in previous single marker phylogenies, we adopted a supermatrix approach by concatenating all checked and reliable sequences available on GenBank as well as new sequences for a total 207 currently recognized genera and eight molecular markers representing a comprehensive coverage of all three genomes. The supermatrix was dated under an uncorrelated relaxed molecular clock using a direct fossil calibration approach. Finally, a lineage-through-time-plot and rates of diversification for the family were generated. The resulting tree, the largest in number of genera and markers sampled to date and covering the whole family in a representative way, provides important insights into the evolution of the family on a broad scale. The backbone of the tree remained largely unresolved and is interpreted as the consequence of early rapid radiation within the family. The age of the family was inferred to be 37.6 (24.2-49.4) Ma, which largely agrees with previous studies. The ages of all major lineages and tribes are also reported. Analysis of diversification suggests that Brassicaceae underwent a rapid period of diversification, after the split with the early diverging tribe Aethionemeae. Given the dates found here, the family appears to have originated under a warm and humid climate approximately 37 Ma. We suggest that the rapid radiation detected was caused by a global cooling during the Oligocene coupled with a genome duplication event. This duplication could have allowed the family to rapidly adapt to the changing climate.

Post-translational modification of an R2R3-MYB transcription factor by a MAP Kinase during xylem development

Despite the pivotal role played by R2R3-MYB family members in the regulation of plant gene expression, little is known about post-translational regulation of these proteins. In animals, the MYB family member, c-MYB, is post-translationally modified by a mitogen-activated protein kinase (MAPK), p42(mapk). In order to test the hypothesis that R2R3-MYB proteins may be regulated by MAPK activity, interplay between a R2R3-MYB family member expressed in differentiating pine xylem (Pinus taeda MYB4, PtMYB4) and MAPK proteins expressed in the same tissue was examined. One of the MAPK proteins expressed in pine xylem, PtMAPK6, phosphorylated PtMYB4. Recombinant PtMAPK6 phosphorylated PtMYB4 on serine-236, located in the C-terminal activation domain of this transcription factor in a context that is found in other plant MYB proteins. Modification of the PtMAPK6 target serine in PtMYB4 did not appear to alter DNA binding in vitro but did alter the ability of PtMYB4 to promote transcriptional activation in yeast. PtMAPK6 activity was detected in developing xylem cells that had ceased cell division and formed secondary walls. Together, the data support a role for PtMAPK6 during early xylem development and suggest a function for this kinase in regulating gene expression through phosphorylation of PtMYB4.

The PREP suite: predictive RNA editors for plant mitochondrial genes, chloroplast genes and user-defined alignments

RNA editing alters plant mitochondrial and chloroplast transcripts by converting specific cytidines to uridines, which usually results in a change in the amino acid sequence of the translated protein. Systematic studies have experimentally identified sites of RNA editing in organellar transcriptomes from several species, but these analyses have not kept pace with rate of genome sequencing. The PREP (predictive RNA editors for plants) suite was developed to computationally predict sites of RNA editing based on the well-known principle that editing in plant organelles increases the conservation of proteins across species. The PREP suite provides predictive RNA editors for plant mitochondrial genes (PREP-Mt), for chloroplast genes (PREP-Cp), and for alignments submitted by the user (PREP-Aln). These servers require minimal input, are very fast, and are highly accurate on all seed plants examined to date. PREP-Mt has proved useful in several research studies and the newly developed PREP-Cp and PREP-Aln servers should be of further assistance for analyses that require knowledge of the location of sites of RNA editing. The PREP suite is freely available at http://prep.unl.edu/.

Sequencing and comparative analysis of a conserved syntenic segment in the Solanaceae

Comparative genomics is a powerful tool for gaining insight into genomic function and evolution. However, in plants, sequence data that would enable detailed comparisons of both coding and noncoding regions have been limited in availability. Here we report the generation and analysis of sequences for an unduplicated conserved syntenic segment (CSS) in the genomes of five members of the agriculturally important plant family Solanaceae. This CSS includes a 105-kb region of tomato chromosome 2 and orthologous regions of the potato, eggplant, pepper, and petunia genomes. With a total neutral divergence of 0.73-0.78 substitutions/site, these sequences are similar enough that most noncoding regions can be aligned, yet divergent enough to be informative about evolutionary dynamics and selective pressures. The CSS contains 17 distinct genes with generally conserved order and orientation, but with numerous small-scale differences between species. Our analysis indicates that the last common ancestor of these species lived approximately 27-36 million years ago, that more than one-third of short genomic segments (5-15 bp) are under selection, and that more than two-thirds of selected bases fall in noncoding regions. In addition, we identify genes under positive selection and analyze hundreds of conserved noncoding elements. This analysis provides a window into 30 million years of plant evolution in the absence of polyploidization.

ESTimating plant phylogeny: lessons from partitioning

BACKGROUND

While Expressed Sequence Tags (ESTs) have proven a viable and efficient way to sample genomes, particularly those for which whole-genome sequencing is impractical, phylogenetic analysis using ESTs remains difficult. Sequencing errors and orthology determination are the major problems when using ESTs as a source of characters for systematics. Here we develop methods to incorporate EST sequence information in a simultaneous analysis framework to address controversial phylogenetic questions regarding the relationships among the major groups of seed plants. We use an automated, phylogenetically derived approach to orthology determination called OrthologID generate a phylogeny based on 43 process partitions, many of which are derived from ESTs, and examine several measures of support to assess the utility of EST data for phylogenies.

RESULTS

A maximum parsimony (MP) analysis resulted in a single tree with relatively high support at all nodes in the tree despite rampant conflict among trees generated from the separate analysis of individual partitions. In a comparison of broader-scale groupings based on cellular compartment (ie: chloroplast, mitochondrial or nuclear) or function, only the nuclear partition tree (based largely on EST data) was found to be topologically identical to the tree based on the simultaneous analysis of all data. Despite topological conflict among the broader-scale groupings examined, only the tree based on morphological data showed statistically significant differences.

CONCLUSION

Based on the amount of character support contributed by EST data which make up a majority of the nuclear data set, and the lack of conflict of the nuclear data set with the simultaneous analysis tree, we conclude that the inclusion of EST data does provide a viable and efficient approach to address phylogenetic questions within a parsimony framework on a genomic scale, if problems of orthology determination and potential sequencing errors can be overcome. In addition, approaches that examine conflict and support in a simultaneous analysis framework allow for a more precise understanding of the evolutionary history of individual process partitions and may be a novel way to understand functional aspects of different kinds of cellular classes of gene products.