Samuel F. Brockington

Angiosperm phylogeny: 17 genes, 640 taxa

PREMISE OF THE STUDY Recent analyses employing up to five genes have provided numerous insights into angiosperm phylogeny, but many relationships have remained unresolved or poorly supported. In the hope of improving our understanding of angiosperm phylogeny, we expanded sampling of taxa and genes beyond previous analyses. METHODS We conducted two primary analyses based on 640 species representing 330 families. The first included 25260 aligned base pairs (bp) from 17 genes (representing all three plant genomes, i.e., nucleus, plastid, and mitochondrion). The second included 19846 aligned bp from 13 genes (representing only the nucleus and plastid). KEY RESULTS Many important questions of deep-level relationships in the nonmonocot angiosperms have now been resolved with strong support. Amborellaceae, Nymphaeales, and Austrobaileyales are successive sisters to the remaining angiosperms (Mesangiospermae), which are resolved into Chloranthales + Magnoliidae as sister to Monocotyledoneae + [Ceratophyllaceae + Eudicotyledoneae]. Eudicotyledoneae contains a basal grade subtending Gunneridae. Within Gunneridae, Gunnerales are sister to the remainder (Pentapetalae), which comprises (1) Superrosidae, consisting of Rosidae (including Vitaceae) and Saxifragales; and (2) Superasteridae, comprising Berberidopsidales, Santalales, Caryophyllales, Asteridae, and, based on this study, Dilleniaceae (although other recent analyses disagree with this placement). Within the major subclades of Pentapetalae, most deep-level relationships are resolved with strong support. CONCLUSIONS Our analyses confirm that with large amounts of sequence data, most deep-level relationships within the angiosperms can be resolved. We anticipate that this well-resolved angiosperm tree will be of broad utility for many areas of biology, including physiology, ecology, paleobiology, and genomics.

Rosid radiation and the rapid rise of angiosperm-dominated forests

The rosid clade (70,000 species) contains more than one-fourth of all angiosperm species and includes most lineages of extant temperate and tropical forest trees. Despite progress in elucidating relationships within the angiosperms, rosids remain the largest poorly resolved major clade; deep relationships within the rosids are particularly enigmatic. Based on parsimony and maximum likelihood (ML) analyses of separate and combined 12-gene (10 plastid genes, 2 nuclear; >18,000 bp) and plastid inverted repeat (IR; 24 genes and intervening spacers; >25,000 bp) datasets for >100 rosid species, we provide a greatly improved understanding of rosid phylogeny. Vitaceae are sister to all other rosids, which in turn form 2 large clades, each with a ML bootstrap value of 100%: (i) eurosids I (Fabidae) include the nitrogen-fixing clade, Celastrales, Huaceae, Zygophyllales, Malpighiales, and Oxalidales; and (ii) eurosids II (Malvidae) include Tapisciaceae, Brassicales, Malvales, Sapindales, Geraniales, Myrtales, Crossosomatales, and Picramniaceae. The rosid clade diversified rapidly into these major lineages, possibly over a period of <15 million years, and perhaps in as little as 4 to 5 million years. The timing of the inferred rapid radiation of rosids [108 to 91 million years ago (Mya) and 107–83 Mya for Fabidae and Malvidae, respectively] corresponds with the rapid rise of angiosperm-dominated forests and the concomitant diversification of other clades that inhabit these forests, including amphibians, ants, placental mammals, and ferns.

Phylogeny of the Caryophyllales sensu lato: revisiting hypotheses on pollination biology and perianth differentiation in the core caryophyllales.

Molecular phylogenetics has revolutionized our understanding of the Caryophyllales, and yet many relationships have remained uncertain, particularly at deeper levels. We have performed parsimony and maximum likelihood analyses on separate and combined data sets comprising nine plastid genes (∼12,000 bp), two nuclear genes (∼5000 bp), and the plastid inverted repeat (∼24,000 bp), giving a combined analyzed length of 42,006 bp for 36 species of Caryophyllales and four outgroups. We have recovered strong support for deep‐level relationships across the order. Two major subclades are well supported, the noncore and core Caryophyllales; Rhabdodendron followed by Simmondsia are sisters to the core Caryophyllales, Limeum and Stegnosperma are successive sisters to the “globular inclusion” clade, Gisekia is a distinct lineage well separated from Rivina within the “raphide” clade, and Rivina and Phytolaccaceae are disparate lineages, with Rivina sister to Nyctaginaceae. The placement of Sarcobatus and relationships within the portulacaceous cohort remain problematic. Within the latter, Halophytum is sister to Basellaceae and Didiereaceae, and the clade comprising Portulaca, Talinum, and Cactaceae is well supported. Classical hypotheses argued that the early Caryophyllales had evolved in open, dry, marginal environments at a time when pollinators were scarce, and, as such, the ancestral caryophyllid flower was wind pollinated with an undifferentiated perianth. We reevaluated these hypotheses in light of our phylogeny and find little support for anemophily as the ancestral condition; however, the early caryophyllid flower is suggested to have possessed an undifferentiated perianth. A subsequent minimum of nine origins of differentiated perianth is inferred. We discuss the evidence for independent origins of differentiated perianth and highlight the research opportunities that this pattern offers to the field of evolutionary developmental genetics.

Dissecting molecular evolution in the highly diverse plant clade Caryophyllales using transcriptome sequencing

Many phylogenomic studies based on transcriptomes have been limited to “single-copy” genes due to methodological challenges in homology and orthology inferences. Only a relatively small number of studies have explored analyses beyond reconstructing species relationships. We sampled 69 transcriptomes in the hyperdiverse plant clade Caryophyllales and 27 outgroups from annotated genomes across eudicots. Using a combined similarity- and phylogenetic tree-based approach, we recovered 10,960 homolog groups, where each was represented by at least eight ingroup taxa. By decomposing these homolog trees, and taking gene duplications into account, we obtained 17,273 ortholog groups, where each was represented by at least ten ingroup taxa. We reconstructed the species phylogeny using a 1,122-gene data set with a gene occupancy of 92.1%. From the homolog trees, we found that both synonymous and nonsynonymous substitution rates in herbaceous lineages are up to three times as fast as in their woody relatives. This is the first time such a pattern has been shown across thousands of nuclear genes with dense taxon sampling. We also pinpointed regions of the Caryophyllales tree that were characterized by relatively high frequencies of gene duplication, including three previously unrecognized whole-genome duplications. By further combining information from homolog tree topology and synonymous distance between paralog pairs, phylogenetic locations for 13 putative genome duplication events were identified. Genes that experienced the greatest gene family expansion were concentrated among those involved in signal transduction and oxidoreduction, including a cytochrome P450 gene that encodes a key enzyme in the betalain synthesis pathway. Our approach demonstrates a new approach for functional phylogenomic analysis in nonmodel species that is based on homolog groups in addition to inferred ortholog groups.

A Targeted Enrichment Strategy for Massively Parallel Sequencing of Angiosperm Plastid Genomes

Premise of the study: We explored a targeted enrichment strategy to facilitate rapid and low-cost next-generation sequencing (NGS) of numerous complete plastid genomes from across the phylogenetic breadth of angiosperms. Methods and Results: A custom RNA probe set including the complete sequences of 22 previously sequenced eudicot plastomes was designed to facilitate hybridization-based targeted enrichment of eudicot plastid genomes. Using this probe set and an Agilent SureSelect targeted enrichment kit, we conducted an enrichment experiment including 24 angiosperms (22 eudicots, two monocots), which were subsequently sequenced on a single lane of the Illumina GAIIx with single-end, 100-bp reads. This approach yielded nearly complete to complete plastid genomes with exceptionally high coverage (mean coverage: 717×), even for the two monocots. Conclusions: Our enrichment experiment was highly successful even though many aspects of the capture process employed were suboptimal. Hence, significant improvements to this methodology are feasible. With this general approach and probe set, it should be possible to sequence more than 300 essentially complete plastid genomes in a single Illumina GAIIx lane (achieving ∼50× mean coverage). However, given the complications of pooling numerous samples for multiplex sequencing and the limited number of barcodes (e.g., 96) available in commercial kits, we recommend 96 samples as a current practical maximum for multiplex plastome sequencing. This high-throughput approach should facilitate large-scale plastid genome sequencing at any level of phylogenetic diversity in angiosperms.

Phylogenetic Analysis of the Plastid Inverted Repeat for 244 Species: Insights into Deeper-Level Angiosperm Relationships from a Long, Slowly Evolving Sequence Region

Recent plastid phylogenomic studies have helped clarify the backbone phylogeny of angiosperms. However, the relatively limited taxon sampling in these studies has precluded strongly supported resolution of some regions of angiosperm phylogeny. Other recent work has suggested that the 25,000-bp plastid inverted repeat (IR) region may be a valuable source of characters for resolving these remaining problematic nodes. Consequently, we aligned all available angiosperm IR sequences to produce a matrix of 24,702 aligned bases for 246 accessions, including 36 new accessions. Maximum likelihood analyses of the complete data set yielded a generally well-supported topology that is highly congruent with those of recent plastid phylogenomic analyses. However, reducing taxon sampling to match a recent 83-gene plastid analysis resulted in significant changes in bootstrap support at some nodes. Notably, IR analyses resolved Pentapetalae into three well-supported clades: (1) superasterids (comprising Santalales, Caryophyllales, Berberidopsidales, and Asteridae), (2) superrosids (comprising Vitaceae, Saxifragales, and Rosidae), and (3) Dilleniaceae. These results provide important new evidence for a stable, well-supported phylogenetic framework for angiosperms and demonstrate the utility of IR data for resolving the deeper levels of angiosperm phylogeny. They also reiterate the importance of carefully considering taxon sampling in phylogenomic studies.

Paralogous Radiations of PIN Proteins with Multiple Origins of Noncanonical PIN Structure

The plant hormone auxin is a conserved regulator of development which has been implicated in the generation of morphological novelty. PIN-FORMED1 (PIN) auxin efflux carriers are central to auxin function by regulating its distribution. PIN family members have divergent structures and cellular localizations, but the origin and evolutionary significance of this variation is unresolved. To characterize PIN family evolution, we have undertaken phylogenetic and structural analyses with a massive increase in taxon sampling over previous studies. Our phylogeny shows that following the divergence of the bryophyte and lycophyte lineages, two deep duplication events gave rise to three distinct lineages of PIN proteins in euphyllophytes. Subsequent independent radiations within each of these lineages were taxonomically asymmetric, giving rise to at least 21 clades of PIN proteins, of which 15 are revealed here for the first time. Although most PIN protein clades share a conserved canonical structure with a modular central loop domain, a small number of noncanonical clades dispersed across the phylogeny have highly divergent protein structure. We propose that PIN proteins underwent sub- and neofunctionalization with substantial modification to protein structure throughout plant evolution. Our results have important implications for plant evolution as they suggest that structurally divergent PIN proteins that arose in paralogous radiations contributed to the convergent evolution of organ systems in different land plant lineages.

Lineage-specific gene radiations underlie the evolution of novel betalain pigmentation in Caryophyllales.

Betalain pigments are unique to the Caryophyllales and structurally and biosynthetically distinct from anthocyanins. Two key enzymes within the betalain synthesis pathway have been identified: 4,5-dioxygenase (DODA) that catalyzes the formation of betalamic acid and CYP76AD1, a cytochrome P450 gene that catalyzes the formation of cyclo-DOPA. We performed phylogenetic analyses to reveal the evolutionary history of the DODA and CYP76AD1 lineages and in the context of an ancestral reconstruction of pigment states we explored the evolution of these genes in relation to the complex evolution of pigments in Caryophylalles. Duplications within the CYP76AD1 and DODA lineages arose just before the origin of betalain pigmentation in the core Caryophyllales. The duplications gave rise to DODA-α and CYP76AD1-α isoforms that appear specific to betalain synthesis. Both betalain-specific isoforms were then lost or downregulated in the anthocyanic Molluginaceae and Caryophyllaceae. Our findings suggest a single origin of the betalain synthesis pathway, with neofunctionalization following gene duplications in the CYP76AD1 and DODA lineages. Loss of DODA-α and CYP76AD1-α in anthocyanic taxa suggests that betalain pigmentation has been lost twice in Caryophyllales, and exclusion of betalain pigments from anthocyanic taxa is mediated through gene loss or downregulation. [Correction added after online publication 13 May 2015: in the last two paragraphs of the Summary the gene name CYP761A was changed to CYP76AD1.]

Conservation and canalization of gene expression during angiosperm diversification accompany the origin and evolution of the flower.

The origin and rapid diversification of the angiosperms (Darwins “Abominable Mystery”) has engaged generations of researchers. Here, we examine the floral genetic programs of phylogenetically pivotal angiosperms (water lily, avocado, California poppy, and Arabidopsis) and a nonflowering seed plant (a cycad) to obtain insight into the origin and subsequent evolution of the flower. Transcriptional cascades with broadly overlapping spatial domains, resembling the hypothesized ancestral gymnosperm program, are deployed across morphologically intergrading organs in water lily and avocado flowers. In contrast, spatially discrete transcriptional programs in distinct floral organs characterize the more recently derived angiosperm lineages represented by California poppy and Arabidopsis. Deep evolutionary conservation in the genetic programs of putatively homologous floral organs traces to those operating in gymnosperm reproductive cones. Female gymnosperm cones and angiosperm carpels share conserved genetic features, which may be associated with the ovule developmental program common to both organs. However, male gymnosperm cones share genetic features with both perianth (sterile attractive and protective) organs and stamens, supporting the evolutionary origin of the floral perianth from the male genetic program of seed plants.

Evolutionary Analysis of the MIXTA Gene Family Highlights Potential Targets for the Study of Cellular Differentiation

Differentiated epidermal cells such as trichomes and conical cells perform numerous essential functions in plant biology and are important for our understanding of developmental patterning and cell shape regulation. Many are also commercially significant, such as cotton fibers and trichomes that secrete pharmaceutically useful or herbivore-deterring compounds. Here, we focus on the phylogeny and evolution of the subgroup 9 R2R3 MYB gene transcription factors, which include the MIXTA gene, and that are important for the specification and regulation of plant cellular differentiation. We have sequenced 49 subgroup 9 R2R3 MYB genes from key experimental taxa and combined these sequences with those identified by an exhaustive bioinformatic search, to compile a data set of 223 subgroup 9 R2R3 MYB genes. Our phylogenetic analyses demonstrate, for the first time, the complex evolutionary history of the subgroup 9 R2R3 MYB genes. A duplication event is inferred before the origin of seed plants giving rise to two major gene lineages, here termed SBG9-A and SBG9-B. The evolutionary conservation of the SBG9-B gene lineage has not been previously recognized and its role in cellular differentiation is unknown, thus an entire clade of potential candidate genes for epidermal cell regulation remains to be explored. Using a heterologous transformation bioassay, we provide functional data that implicate members of the SBG9-B lineage in the specification of epidermal projections. Furthermore, we reveal numerous putative duplication events in both SBG9-A and SBG9-B lineages, resolving uncertainty about orthology and paralogy among the subgroup 9 R2R3 MYB genes. Finally, we provide a robust framework over which to interpret existing functional data and to direct ongoing comparative genetic research into the evolution of plant cellular diversity.