Diversification Models Conflate Likelihood and Prior, and Cannot be Compared Using Conventional Model-Comparison Tools.

Time-calibrated phylogenetic trees are a tremendously powerful tool for studying evolutionary, ecological, and epidemiological phenomena. Such trees are predominantly inferred in a Bayesian framework, with the phylogeny itself treated as a parameter with a prior distribution (a "tree prior"). However, we show that the tree "parameter" consists, in part, of data, in the form of taxon samples. Treating the tree as a parameter fails to account for these data and compromises our ability to compare among models using standard techniques (e.g., marginal likelihoods estimated using path-sampling and stepping-stone sampling algorithms). Since accuracy of the inferred phylogeny strongly depends on how well the tree prior approximates the true diversification process that gave rise to the tree, the inability to accurately compare competing tree priors has broad implications for applications based on time-calibrated trees. We outline potential remedies to this problem, and provide guidance for researchers interested in assessing the fit of tree models. [Bayes factors; Bayesian model comparison; birth-death models; divergence-time estimation; lineage diversification].

Unearthing Modes of Climatic Adaptation in Underground Storage Organs Across Liliales.

Testing adaptive hypotheses about how continuous traits evolve in association with developmentally structured discrete traits, while accounting for the confounding influence of other, hidden, evolutionary forces, remains a challenge in evolutionary biology. For example, geophytes are herbaceous plants-with underground buds-that use underground storage organs (USOs) to survive extended periods of unfavorable conditions. Such plants have evolved multiple times independently across all major vascular plant lineages. Even within closely related lineages, however, geophytes show impressive variation in the morphological modifications and structures (i.e.,"types" of USOs) that allow them to survive underground. Despite the developmental and structural complexity of USOs, the prevailing hypothesis is that they represent convergent evolutionary "solutions" to a common ecological problem, though some recent research has drawn this conclusion into question. We extend existing phylogenetic comparative methods to test for links between the hierarchical discrete morphological traits associated with USOs and adaptation to environmental variables, using a phylogeny of 621 species in Liliales. We found that plants with different USO types do not differ in climatic niche more than expected by chance, with the exception of root morphology, where modified roots are associated with lower temperature seasonality. These findings suggest that root tubers may reflect adaptations to different climatic conditions than those represented by other types of USOs. Thus, the tissue type and developmental origin of the USO structure may influence the way it mediates ecological relationships, which draws into question the appropriateness of ascribing broad ecological patterns uniformly across geophytic taxa. This work provides a new framework for testing adaptive hypotheses and for linking ecological patterns across morphologically varying taxa while accounting for developmental (non-independent) relationships in morphological data. [Climatic niche evolution; geophytes; imperfect correspondence; macroevolution.].

Polyploid goldback and silverback ferns (Pentagramma) occupy a wider, colder, and wetter bioclimatic niche than diploid counterparts.

PREMISE: The western North American fern genus Pentagramma (Pteridaceae) is characterized by complex patterns of ploidy variation, an understanding of which is critical to comprehending both the evolutionary processes within the genus and its current diversity. METHODS: We undertook a cytogeographic study across the range of the genus, using a combination of chromosome counts and flow cytometry to infer ploidy level. Bioclimatic variables and elevation were used to compare niches. RESULTS: We found that diploids and tetraploids are common and widespread, and triploids are rare and sporadic; in contrast with genome size inferences in earlier studies, no hexaploids were found. Diploids and tetraploids show different geographic ranges: only tetraploids were found in the northernmost portion of the range (Washington, Oregon, and British Columbia) and only diploids were found in the Sierra Nevada of California. Diploid, triploid, and tetraploid cytotypes were found to co-occur in relatively few localities: in the southern (San Diego County, California) and desert Southwest (Arizona) parts of the range, and along the Pacific Coast of California. CONCLUSIONS: Tetraploids occupy a wider bioclimatic niche than diploids both within P. triangularis and at the genus-wide scale. It is unknown whether the wider niche of tetraploids is due to their expansion upon the diploid niche, if diploids have contracted their niche due to competition or changing abiotic conditions, or if this wider niche occupancy is due to multiple origins of tetraploids.

Inferring the Total-Evidence Timescale of Marattialean Fern Evolution in the Face of Model Sensitivity.

Phylogenetic divergence-time estimation has been revolutionized by two recent developments: 1) total-evidence dating (or "tip-dating") approaches that allow for the incorporation of fossils as tips in the analysis, with their phylogenetic and temporal relationships to the extant taxa inferred from the data and 2) the fossilized birth-death (FBD) class of tree models that capture the processes that produce the tree (speciation, extinction, and fossilization) and thus provide a coherent and biologically interpretable tree prior. To explore the behavior of these methods, we apply them to marattialean ferns, a group that was dominant in Carboniferous landscapes prior to declining to its modest extant diversity of slightly over 100 species. We show that tree models have a dramatic influence on estimates of both divergence times and topological relationships. This influence is driven by the strong, counter-intuitive informativeness of the uniform tree prior, and the inherent nonidentifiability of divergence-time models. In contrast to the strong influence of the tree models, we find minor effects of differing the morphological transition model or the morphological clock model. We compare the performance of a large pool of candidate models using a combination of posterior-predictive simulation and Bayes factors. Notably, an FBD model with epoch-specific speciation and extinction rates was strongly favored by Bayes factors. Our best-fitting model infers stem and crown divergences for the Marattiales in the mid-Devonian and Late Cretaceous, respectively, with elevated speciation rates in the Mississippian and elevated extinction rates in the Cisuralian leading to a peak diversity of ${\sim}$2800 species at the end of the Carboniferous, representing the heyday of the Psaroniaceae. This peak is followed by the rapid decline and ultimate extinction of the Psaroniaceae, with their descendants, the Marattiaceae, persisting at approximately stable levels of diversity until the present. This general diversification pattern appears to be insensitive to potential biases in the fossil record; despite the preponderance of available fossils being from Pennsylvanian coal balls, incorporating fossilization-rate variation does not improve model fit. In addition, by incorporating temporal data directly within the model and allowing for the inference of the phylogenetic position of the fossils, our study makes the surprising inference that the clade of extant Marattiales is relatively young, younger than any of the fossils historically thought to be congeneric with extant species. This result is a dramatic demonstration of the dangers of node-based approaches to divergence-time estimation, where the assignment of fossils to particular clades is made a priori (earlier node-based studies that constrained the minimum ages of extant genera based on these fossils resulted in much older age estimates than in our study) and of the utility of explicit models of morphological evolution and lineage diversification. [Bayesian model comparison; Carboniferous; divergence-time estimation; fossil record; fossilized birth-death; lineage diversification; Marattiales; models of morphological evolution; Psaronius; RevBayes.].

Comparative transcriptomics of a monocotyledonous geophyte reveals shared molecular mechanisms of underground storage organ formation.

Many species from across the vascular plant tree-of-life have modified standard plant tissues into tubers, bulbs, corms, and other underground storage organs (USOs), unique innovations which allow these plants to retreat underground. Our ability to understand the developmental and evolutionary forces that shape these morphologies is limited by a lack of studies on certain USOs and plant clades. We take a comparative transcriptomics approach to characterizing the molecular mechanisms of tuberous root formation in Bomarea multiflora (Alstroemeriaceae) and compare these mechanisms to those identified in other USOs across diverse plant lineages; B. multiflora fills a key gap in our understanding of USO molecular development as the first monocot with tuberous roots to be the focus of this kind of research. We sequenced transcriptomes from the growing tip of four tissue types (aerial shoot, rhizome, fibrous root, and root tuber) of three individuals of B. multiflora. We identified differentially expressed isoforms between tuberous and non-tuberous roots and tested the expression of a priori candidate genes implicated in underground storage in other taxa. We identify 271 genes that are differentially expressed in root tubers versus non-tuberous roots, including genes implicated in cell wall modification, defense response, and starch biosynthesis. We also identify a phosphatidylethanolamine-binding protein, which has been implicated in tuberization signalling in other taxa and, through gene-tree analysis, place this copy in a phylogenetic context. These findings suggest that some similar molecular processes underlie the formation of USOs across flowering plants despite the long evolutionary distances among taxa and non-homologous morphologies (e.g., bulbs vs. tubers). (Plant development, tuberous roots, comparative transcriptomics, geophytes).

Polyploid phylogenetics.

Polyploidy is a dominant feature of extant plant diversity. However, major research questions, including whether polyploidy is important to long-term evolution or is just 'evolutionary noise', remain unresolved due to difficulties associated with the generation and analysis of data from polyploid lineages. Many of these difficulties have been recently overcome, such that it is now often relatively straightforward to infer the full and often reticulate phylogenetic history of groups with recently formed polyploids. This nascent field of 'polyploid phylogenetics' allows researchers to tackle long-standing questions of polyploid macroevolution, supplies the foundation for mechanistic models of ploidy change, and provides the opportunity to include a more complete and representative sample of plant taxa in our analyses in general.

A global plastid phylogeny of the fern genus Asplenium (Aspleniaceae).

The infrageneric relationships and taxonomy of the largest fern genus, Asplenium (Aspleniaceae), have remained poorly understood. Previous studies have focused mainly on specific species complexes involving a few or dozens of species only, or have achieved a large taxon sampling but only one plastid marker was used. In the present study, DNA sequences from six plastid markers (atpB, rbcL, rps4, rps4-trnS, trnL and trnL-F) of 1030 accessions (616 of them newly sequenced here) representing c. 420 species of Asplenium (60% of estimated species diversity), 16 species of Hymenasplenium, three Diplaziopsidaceae, and four Rhachidosoraceae were used to produce the largest genus-level phylogeny yet for ferns. Our major results include: (i) Asplenium as broadly circumscribed is monophyletic based on our inclusion of representatives of 32 of 38 named segregate genera; (ii) 11 major clades in Asplenium are identified, and their relationships are mostly well-resolved and strongly supported; (iii) numerous species, unsampled in previous studies, suggest new relationships and numerous cryptic species and species complexes in Asplenium; and (iv) the accrued molecular evidence provides an essential foundation for further investigations of complex patterns of geographical diversification, speciation and reticulate evolution in this family.

Infraspecific diversification of the star cloak fern (Notholaena standleyi) in the deserts of the United States and Mexico.

PREMISE: Not all ferns grow in moist and shaded habitats. One well-known example is Notholaena standleyi, a species that thrives in deserts of the southwestern United States and Mexico. This species exhibits several "chemotypes" that differ in farina (flavonoid exudates) color and chemistry. By integrating data from molecular phylogenetics, cytology, biochemistry, and biogeography, we circumscribed the major evolutionary lineages within N. standleyi and reconstructed their diversification histories. METHODS: Forty-eight samples were selected from across the geographic distribution of N. standleyi. Phylogenetic relationships were inferred using four plastid and five nuclear markers. Ploidy levels were inferred using spore sizes calibrated by chromosome counts, and farina chemistry was compared using thin-layer chromatography. RESULTS: Four clades are recognized, three of which roughly correspond to previously recognized chemotypes. The diploid clades G and Y are found in the Sonoran and Chihuahuan deserts, respectively; they are estimated to have diverged in the Pleistocene, congruent with the postulated timing of climatological events separating these two deserts. Clade P/YG is tetraploid and partially overlaps the distribution of clade Y in the eastern Chihuahuan Desert. It is apparently confined to limestone, a geologic substrate rarely occupied by members of the other clades. The cryptic (C) clade, a diploid group known only from southern Mexico and highly disjunct from the other three clades, is newly recognized here. CONCLUSIONS: Our results reveal a complex intraspecific diversification history of N. standleyi, traceable to a variety of evolutionary drivers including classic allopatry, parapatry with or without changes in geologic substrate, and sympatric divergence through polyploidization.

Organellomic data sets confirm a cryptic consensus on (unrooted) land-plant relationships and provide new insights into bryophyte molecular evolution.

PREMISE: Phylogenetic trees of bryophytes provide important evolutionary context for land plants. However, published inferences of overall embryophyte relationships vary considerably. We performed phylogenomic analyses of bryophytes and relatives using both mitochondrial and plastid gene sets, and investigated bryophyte plastome evolution. METHODS: We employed diverse likelihood-based analyses to infer large-scale bryophyte phylogeny for mitochondrial and plastid data sets. We tested for changes in purifying selection in plastid genes of a mycoheterotrophic liverwort (Aneura mirabilis) and a putatively mycoheterotrophic moss (Buxbaumia), and compared 15 bryophyte plastomes for major structural rearrangements. RESULTS: Overall land-plant relationships conflict across analyses, generally weakly. However, an underlying (unrooted) four-taxon tree is consistent across most analyses and published studies. Despite gene coverage patchiness, relationships within mosses, liverworts, and hornworts are largely congruent with previous studies, with plastid results generally better supported. Exclusion of RNA edit sites restores cases of unexpected non-monophyly to monophyly for Takakia and two hornwort genera. Relaxed purifying selection affects multiple plastid genes in mycoheterotrophic Aneura but not Buxbaumia. Plastid genome structure is nearly invariant across bryophytes, but the tufA locus, presumed lost in embryophytes, is unexpectedly retained in several mosses. CONCLUSIONS: A common unrooted tree underlies embryophyte phylogeny, [(liverworts, mosses), (hornworts, vascular plants)]; rooting inconsistency across studies likely reflects substantial distance to algal outgroups. Analyses combining genomic and transcriptomic data may be misled locally for heavily RNA-edited taxa. The Buxbaumia plastome lacks hallmarks of relaxed selection found in mycoheterotrophic Aneura. Autotrophic bryophyte plastomes, including Buxbaumia, hardly vary in overall structure.

The Physcomitrella patens chromosome-scale assembly reveals moss genome structure and evolution.

The draft genome of the moss model, Physcomitrella patens, comprised approximately 2000 unordered scaffolds. In order to enable analyses of genome structure and evolution we generated a chromosome-scale genome assembly using genetic linkage as well as (end) sequencing of long DNA fragments. We find that 57% of the genome comprises transposable elements (TEs), some of which may be actively transposing during the life cycle. Unlike in flowering plant genomes, gene- and TE-rich regions show an overall even distribution along the chromosomes. However, the chromosomes are mono-centric with peaks of a class of Copia elements potentially coinciding with centromeres. Gene body methylation is evident in 5.7% of the protein-coding genes, typically coinciding with low GC and low expression. Some giant virus insertions are transcriptionally active and might protect gametes from viral infection via siRNA mediated silencing. Structure-based detection methods show that the genome evolved via two rounds of whole genome duplications (WGDs), apparently common in mosses but not in liverworts and hornworts. Several hundred genes are present in colinear regions conserved since the last common ancestor of plants. These syntenic regions are enriched for functions related to plant-specific cell growth and tissue organization. The P. patens genome lacks the TE-rich pericentromeric and gene-rich distal regions typical for most flowering plant genomes. More non-seed plant genomes are needed to unravel how plant genomes evolve, and to understand whether the P. patens genome structure is typical for mosses or bryophytes.

Phylogeny of Paullinia L. (Paullinieae: Sapindaceae), a diverse genus of lianas with dynamic fruit evolution.

Paullinia L. is a genus of c. 220 mostly Neotropical forest-dwelling lianas that display a wide diversity of fruit morphologies. Paullinia resembles other members of the Paullinieae tribe in being a climber with stipulate compound leaves and paired inflorescence tendrils. However, it is distinct in having capsular fruits with woody, coriaceous, or crustaceous pericarps. While consistent in this basic plan, the pericarps of Paullinia fruits are otherwise highly variable-in some species they are winged, whereas in others they are without wings or covered with spines. With the exception of the water-dispersed indehiscent spiny fruits of some members of Paullinia sect. Castanella, all species are dehiscent, opening their capsules while they are still attached to the branch, to reveal arillate animal-dispersed seeds. Here we present a molecular phylogeny of Paullinia derived from 11 molecular markers, including nine newly-developed single-copy nuclear markers amplified by microfluidics PCR. This is the first broadly sampled molecular phylogeny for the genus. Paullinia is supported as monophyletic and is sister to Cardiospermum L., which together are sister to Serjania Mill + Urvillea Kunth. We apply this novel phylogenetic hypothesis to test previous infrageneric classifications and to infer that unwinged fruits represent the ancestral condition, from which there were repeated evolutionary transitions and reversals. However, because the seeds of both winged and unwinged fruits are dispersed by animals, we conclude that the repeated transitions in fruit morphology may relate to visual display strategies to attract animal dispersers, and do not represent transitions to wind dispersal.

Low-copy nuclear sequence data confirm complex patterns of farina evolution in notholaenid ferns (Pteridaceae).

Notholaenids are an unusual group of ferns that have adapted to, and diversified within, the deserts of Mexico and the southwestern United States. With approximately 40 species, this group is noted for being desiccation-tolerant and having "farina"-powdery exudates of lipophilic flavonoid aglycones-that occur on both the gametophytic and sporophytic phases of their life cycle. The most recent circumscription of notholaenids based on plastid markers surprisingly suggests that several morphological characters, including the expression of farina, are homoplasious. In a striking case of convergence, Notholaena standleyi appears to be distantly related to core Notholaena, with several taxa not before associated with Notholaena nested between them. Such conflicts can be due to morphological homoplasy resulting from adaptive convergence or, alternatively, the plastid phylogeny itself might be misleading, diverging from the true species tree due to incomplete lineage sorting, hybridization, or other factors. In this study, we present a species phylogeny for notholaenid ferns, using four low-copy nuclear loci and concatenated data from three plastid loci. A total of 61 individuals (49 notholaenids and 12 outgroup taxa) were sampled, including 31 out of 37 recognized notholaenid species. The homeologous/allelic nuclear sequences were retrieved using PacBio sequencing and the PURC bioinformatics pipeline. Each dataset was first analyzed individually using maximum likelihood and Bayesian inference, and the species phylogeny was inferred using *BEAST. Although we observed several incongruences between the nuclear and plastid phylogenies, our principal results are broadly congruent with previous inferences based on plastid data. By mapping the presence of farina and their biochemical constitutions on our consensus phylogenetic tree, we confirmed that the characters are indeed homoplastic and have complex evolutionary histories. Hybridization among recognized species of the notholaenid clade appears to be relatively rare compared to that observed in other well-studied fern genera.

Genetic Analysis of Physcomitrella patens Identifies ABSCISIC ACID NON-RESPONSIVE, a Regulator of ABA Responses Unique to Basal Land Plants and Required for Desiccation Tolerance.

The anatomically simple plants that first colonized land must have acquired molecular and biochemical adaptations to drought stress. Abscisic acid (ABA) coordinates responses leading to desiccation tolerance in all land plants. We identified ABA nonresponsive mutants in the model bryophyte Physcomitrella patens and genotyped a segregating population to map and identify the ABA NON-RESPONSIVE (ANR) gene encoding a modular protein kinase comprising an N-terminal PAS domain, a central EDR domain, and a C-terminal MAPKKK-like domain. anr mutants fail to accumulate dehydration tolerance-associated gene products in response to drought, ABA, or osmotic stress and do not acquire ABA-dependent desiccation tolerance. The crystal structure of the PAS domain, determined to 1.7-Å resolution, shows a conserved PAS-fold that dimerizes through a weak dimerization interface. Targeted mutagenesis of a conserved tryptophan residue within the PAS domain generates plants with ABA nonresponsive growth and strongly attenuated ABA-responsive gene expression, whereas deleting this domain retains a fully ABA-responsive phenotype. ANR orthologs are found in early-diverging land plant lineages and aquatic algae but are absent from more recently diverged vascular plants. We propose that ANR genes represent an ancestral adaptation that enabled drought stress survival of the first terrestrial colonizers but were lost during land plant evolution.

Widespread co-occurrence of multiple ploidy levels in fragile ferns (Cystopteris fragilis complex; Cystopteridaceae) probably stems from similar ecology of cytotypes, their efficient dispersal and inter-ploidy hybridization.

BACKGROUND AND AIMS: Polyploidy has played an important role in the evolution of ferns. However, the dearth of data on cytotype diversity, cytotype distribution patterns and ecology in ferns is striking in comparison with angiosperms and prevents an assessment of whether cytotype coexistence and its mechanisms show similar patterns in both plant groups. Here, an attempt to fill this gap was made using the ploidy-variable and widely distributed Cystopteris fragilis complex. METHODS: Flow cytometry was used to assess DNA ploidy level and monoploid genome size (Cx value) of 5518 C. fragilis individuals from 449 populations collected over most of the species' global distributional range, supplemented with data from 405 individuals representing other related species from the complex. Ecological preferences of C. fragilis tetraploids and hexaploids were compared using field-recorded parameters and database-extracted climate data. KEY RESULTS: Altogether, five different ploidy levels (2x, 4x, 5x, 6x, 8x) were detected and three species exhibited intraspecific ploidy-level variation: C. fragilis, C. alpina and C. diaphana. Two predominant C. fragilis cytotypes, tetraploids and hexaploids, co-occur over most of Europe in a diffuse, mosaic-like pattern. Within this contact zone, 40 % of populations were mixed-ploidy and most also contained pentaploid hybrids. Environmental conditions had only a limited effect on the distribution of cytotypes. Differences were found in the Cx value of tetraploids and hexaploids: between-cytotype divergence was higher in uniform-ploidy than in mixed-ploidy populations. CONCLUSIONS: High ploidy-level diversity and widespread cytotype coexistence in the C. fragilis complex match the well-documented patterns in some angiosperms. While ploidy coexistence in C. fragilis is not driven by environmental factors, it could be facilitated by the perennial life-form of the species, its reproductive modes and efficient wind dispersal of spores. Independent origins of hexaploids and/or inter-ploidy gene flow may be expected in mixed-ploidy populations according to Cx value comparisons.

Super-resolution ribosome profiling reveals unannotated translation events in Arabidopsis.

Deep sequencing of ribosome footprints (ribosome profiling) maps and quantifies mRNA translation. Because ribosomes decode mRNA every 3 nt, the periodic property of ribosome footprints could be used to identify novel translated ORFs. However, due to the limited resolution of existing methods, the 3-nt periodicity is observed mostly in a global analysis, but not in individual transcripts. Here, we report a protocol applied to Arabidopsis that maps over 90% of the footprints to the main reading frame and thus offers super-resolution profiles for individual transcripts to precisely define translated regions. The resulting data not only support many annotated and predicted noncanonical translation events but also uncover small ORFs in annotated noncoding RNAs and pseudogenes. A substantial number of these unannotated ORFs are evolutionarily conserved, and some produce stable proteins. Thus, our study provides a valuable resource for plant genomics and an efficient optimization strategy for ribosome profiling in other organisms.

Microbial-type terpene synthase genes occur widely in nonseed land plants, but not in seed plants.

The vast abundance of terpene natural products in nature is due to enzymes known as terpene synthases (TPSs) that convert acyclic prenyl diphosphate precursors into a multitude of cyclic and acyclic carbon skeletons. Yet the evolution of TPSs is not well understood at higher levels of classification. Microbial TPSs from bacteria and fungi are only distantly related to typical plant TPSs, whereas genes similar to microbial TPS genes have been recently identified in the lycophyte Selaginella moellendorffii The goal of this study was to investigate the distribution, evolution, and biochemical functions of microbial terpene synthase-like (MTPSL) genes in other plants. By analyzing the transcriptomes of 1,103 plant species ranging from green algae to flowering plants, putative MTPSL genes were identified predominantly from nonseed plants, including liverworts, mosses, hornworts, lycophytes, and monilophytes. Directed searching for MTPSL genes in the sequenced genomes of a wide range of seed plants confirmed their general absence in this group. Among themselves, MTPSL proteins from nonseed plants form four major groups, with two of these more closely related to bacterial TPSs and the other two to fungal TPSs. Two of the four groups contain a canonical aspartate-rich "DDxxD" motif. The third group has a "DDxxxD" motif, and the fourth group has only the first two "DD" conserved in this motif. Upon heterologous expression, representative members from each of the four groups displayed diverse catalytic functions as monoterpene and sesquiterpene synthases, suggesting these are important for terpene formation in nonseed plants.

Insights into the Evolution of Hydroxyproline-Rich Glycoproteins from 1000 Plant Transcriptomes.

The carbohydrate-rich cell walls of land plants and algae have been the focus of much interest given the value of cell wall-based products to our current and future economies. Hydroxyproline-rich glycoproteins (HRGPs), a major group of wall glycoproteins, play important roles in plant growth and development, yet little is known about how they have evolved in parallel with the polysaccharide components of walls. We investigate the origins and evolution of the HRGP superfamily, which is commonly divided into three major multigene families: the arabinogalactan proteins (AGPs), extensins (EXTs), and proline-rich proteins. Using motif and amino acid bias, a newly developed bioinformatics pipeline, we identified HRGPs in sequences from the 1000 Plants transcriptome project (www.onekp.com). Our analyses provide new insights into the evolution of HRGPs across major evolutionary milestones, including the transition to land and the early radiation of angiosperms. Significantly, data mining reveals the origin of glycosylphosphatidylinositol (GPI)-anchored AGPs in green algae and a 3- to 4-fold increase in GPI-AGPs in liverworts and mosses. The first detection of cross-linking (CL)-EXTs is observed in bryophytes, which suggests that CL-EXTs arose though the juxtaposition of preexisting SPn EXT glycomotifs with refined Y-based motifs. We also detected the loss of CL-EXT in a few lineages, including the grass family (Poaceae), that have a cell wall composition distinct from other monocots and eudicots. A key challenge in HRGP research is tracking individual HRGPs throughout evolution. Using the 1000 Plants output, we were able to find putative orthologs of Arabidopsis pollen-specific GPI-AGPs in basal eudicots.

Rapid allopolyploid radiation of moonwort ferns (Botrychium; Ophioglossaceae) revealed by PacBio sequencing of homologous and homeologous nuclear regions.

Polyploidy is a major speciation process in vascular plants, and is postulated to be particularly important in shaping the diversity of extant ferns. However, limitations in the availability of bi-parental markers for ferns have greatly limited phylogenetic investigation of polyploidy in this group. With a large number of allopolyploid species, the genus Botrychium is a classic example in ferns where recurrent polyploidy is postulated to have driven frequent speciation events. Here, we use PacBio sequencing and the PURC bioinformatics pipeline to capture all homeologous or allelic copies of four long (∼1 kb) low-copy nuclear regions from a sample of 45 specimens (25 diploids and 20 polyploids) representing 37 Botrychium taxa, and three outgroups. This sample includes most currently recognized Botrychium species in Europe and North America, and the majority of our specimens were genotyped with co-dominant nuclear allozymes to ensure species identification. We analyzed the sequence data using maximum likelihood (ML) and Bayesian inference (BI) concatenated-data ("gene tree") approaches to explore the relationships among Botrychium species. Finally, we estimated divergence times among Botrychium lineages and inferred the multi-labeled polyploid species tree showing the origins of the polyploid taxa, and their relationships to each other and to their diploid progenitors. We found strong support for the monophyly of the major lineages within Botrychium and identified most of the parental donors of the polyploids; these results largely corroborate earlier morphological and allozyme-based investigations. Each polyploid had at least two distinct homeologs, indicating that all sampled polyploids are likely allopolyploids (rather than autopolyploids). Our divergence-time analyses revealed that these allopolyploid lineages originated recently-within the last two million years-and thus that the genus has undergone a recent radiation, correlated with multiple independent allopolyploidizations across the phylogeny. Also, we found strong parental biases in the formation of allopolyploids, with individual diploid species participating multiple times as either the maternal or paternal donor (but not both). Finally, we discuss the role of polyploidy in the evolutionary history of Botrychium and the interspecific reproductive barriers possibly involved in these parental biases.

Inferring the evolutionary reduction of corm lobation in Isoëtes using Bayesian model-averaged ancestral state reconstruction.

PREMISE OF THE STUDY: Inferring the evolution of characters in Isoëtes has been problematic, as these plants are morphologically conservative and yet highly variable and homoplasious within that conserved base morphology. However, molecular phylogenies have given us a valuable tool for testing hypotheses of character evolution within the genus, such as the hypothesis of ongoing morphological reductions. METHODS: We examined the reduction in lobe number on the underground trunk, or corm, by combining the most recent molecular phylogeny with morphological descriptions gathered from the literature and observations of living specimens. Ancestral character states were inferred using nonstationary evolutionary models, reversible-jump MCMC, and Bayesian model averaging. KEY RESULTS: Our results support the hypothesis of a directional reduction in lobe number in Isoëtes, with the best-supported model of character evolution being one of irreversible reduction. Furthermore, the most probable ancestral corm lobe number of extant Isoëtes is three, and a reduction to two lobes has occurred at least six times. CONCLUSIONS: From our results, we can infer that corm lobation, like many other traits in Isoëtes, shows a degree of homoplasy, and yet also shows ongoing evolutionary reduction.

A novel chloroplast gene reported for flagellate plants.

PREMISE OF THE STUDY: Gene space in plant plastid genomes is well characterized and annotated, yet we discovered an unrecognized open reading frame (ORF) in the fern lineage that is conserved across flagellate plants. METHODS: We initially detected a putative uncharacterized ORF by the existence of a highly conserved region between rps16 and matK in a series of matK alignments of leptosporangiate ferns. We mined available plastid genomes for this ORF, which we now refer to as ycf94, to infer evolutionary selection pressures and assist in functional prediction. To further examine the transcription of ycf94, we assembled the plastid genome and sequenced the transcriptome of the leptosporangiate fern Adiantum shastense Huiet & A.R. Sm. KEY RESULTS: The ycf94 predicted protein has a distinct transmembrane domain but with no sequence homology to other proteins with known function. The nonsynonymous/synonymous substitution rate ratio of ycf94 is on par with other fern plastid protein-encoding genes, and additional homologs can be found in a few lycophyte, moss, hornwort, and liverwort plastid genomes. Homologs of ycf94 were not found in seed plants. In addition, we report a high level of RNA editing for ycf94 transcripts-a hallmark of protein-coding genes in fern plastomes. CONCLUSIONS: The degree of sequence conservation, together with the presence of a distinct transmembrane domain and RNA-editing sites, suggests that ycf94 is a protein-coding gene of functional significance in ferns and, potentially, bryophytes and lycophytes. However, the origin and exact function of this gene require further investigation.