Extraordinary preservation of gene collinearity over three hundred million years revealed in homosporous lycophytes.

Homosporous lycophytes (Lycopodiaceae) are a deeply diverged lineage in the plant tree of life, having split from heterosporous lycophytes (Selaginella and Isoetes) ~400 Mya. Compared to the heterosporous lineage, Lycopodiaceae has markedly larger genome sizes and remains the last major plant clade for which no chromosome-level assembly has been available. Here, we present chromosomal genome assemblies for two homosporous lycophyte species, the allotetraploid Huperzia asiatica and the diploid Diphasiastrum complanatum. Remarkably, despite that the two species diverged ~350 Mya, around 30% of the genes are still in syntenic blocks. Furthermore, both genomes had undergone independent whole genome duplications, and the resulting intragenomic syntenies have likewise been preserved relatively well. Such slow genome evolution over deep time is in stark contrast to heterosporous lycophytes and is correlated with a decelerated rate of nucleotide substitution. Together, the genomes of H. asiatica and D. complanatum not only fill a crucial gap in the plant genomic landscape but also highlight a potentially meaningful genomic contrast between homosporous and heterosporous species.

Green plant genomes: What we know in an era of rapidly expanding opportunities.

Green plants play a fundamental role in ecosystems, human health, and agriculture. As de novo genomes are being generated for all known eukaryotic species as advocated by the Earth BioGenome Project, increasing genomic information on green land plants is essential. However, setting standards for the generation and storage of the complex set of genomes that characterize the green lineage of life is a major challenge for plant scientists. Such standards will need to accommodate the immense variation in green plant genome size, transposable element content, and structural complexity while enabling research into the molecular and evolutionary processes that have resulted in this enormous genomic variation. Here we provide an overview and assessment of the current state of knowledge of green plant genomes. To date fewer than 300 complete chromosome-scale genome assemblies representing fewer than 900 species have been generated across the estimated 450,000 to 500,000 species in the green plant clade. These genomes range in size from 12 Mb to 27.6 Gb and are biased toward agricultural crops with large branches of the green tree of life untouched by genomic-scale sequencing. Locating suitable tissue samples of most species of plants, especially those taxa from extreme environments, remains one of the biggest hurdles to increasing our genomic inventory. Furthermore, the annotation of plant genomes is at present undergoing intensive improvement. It is our hope that this fresh overview will help in the development of genomic quality standards for a cohesive and meaningful synthesis of green plant genomes as we scale up for the future.

150 Years of Coevolution Research: Evolution and Ecology of Yucca Moths (Prodoxidae) and Their Hosts.

Yucca moths (Tegeticula and Parategeticula) are specialized pollinators of yucca plants, possessing unique, tentacle-like mouthparts used to actively collect pollen and deposit it onto the flowers of their hosts. The moths' larvae feed on the developing seeds and fruit tissue. First described in 1873, the yucca-yucca moth pollination system is now considered the archetypical example of a coevolved intimate mutualism. Research conducted over the past three decades has transformed our understanding of yucca moth diversity and host plant interactions. We summarize the current understanding of the diversity, ecology, and evolution of this group, review evidence for coevolution of the insects and their hosts, and describe how the nature of the interaction varies across evolutionary time and ecological contexts. Finally, we identify unresolved questions and areas for future research.

Interchromosomal linkage disequilibrium and linked fitness cost loci associated with selection for herbicide resistance.

The adaptation of weeds to herbicide is both a significant problem in agriculture and a model of rapid adaptation. However, significant gaps remain in our knowledge of resistance controlled by many loci and the evolutionary factors that influence the maintenance of resistance. Here, using herbicide-resistant populations of the common morning glory (Ipomoea purpurea), we perform a multilevel analysis of the genome and transcriptome to uncover putative loci involved in nontarget-site herbicide resistance (NTSR) and to examine evolutionary forces underlying the maintenance of resistance in natural populations. We found loci involved in herbicide detoxification and stress sensing to be under selection and confirmed that detoxification is responsible for glyphosate (RoundUp) resistance using a functional assay. We identified interchromosomal linkage disequilibrium (ILD) among loci under selection reflecting either historical processes or additive effects leading to the resistance phenotype. We further identified potential fitness cost loci that were strongly linked to resistance alleles, indicating the role of genetic hitchhiking in maintaining the cost. Overall, our work suggests that NTSR glyphosate resistance in I. purpurea is conferred by multiple genes which are potentially maintained through generations via ILD, and that the fitness cost associated with resistance in this species is likely a by-product of genetic hitchhiking.

A derived ZW chromosome system in Amborella trichopoda, representing the sister lineage to all other extant flowering plants.

The genetic basis and evolution of sex determination in dioecious plants is emerging as an active area of research with exciting advances in genome sequencing and analysis technologies. As the sole species within the sister lineage to all other extant flowering plants, Amborella trichopoda is an important model for understanding the evolution and development of flowers. Plants typically produce only male or female flowers, but sex determination mechanisms are unknown for the species. Sequence data derived from plants of natural origin and an F1 mapping population were used to identify sex-linked genes and the nonrecombining region. Amborella trichopoda has a ZW sex determination system. Analysis of genes in a 4 Mb nonrecombining sex-determination region reveals recent divergence of Z and W gametologs, and few Z- and W-specific genes. The sex chromosomes of A. trichopoda evolved less than 16.5 Myr ago, long after the divergence of the extant angiosperms.

Evolution of chloroplast retrograde signaling facilitates green plant adaptation to land.

Chloroplast retrograde signaling networks are vital for chloroplast biogenesis, operation, and signaling, including excess light and drought stress signaling. To date, retrograde signaling has been considered in the context of land plant adaptation, but not regarding the origin and evolution of signaling cascades linking chloroplast function to stomatal regulation. We show that key elements of the chloroplast retrograde signaling process, the nucleotide phosphatase (SAL1) and 3'-phosphoadenosine-5'-phosphate (PAP) metabolism, evolved in streptophyte algae-the algal ancestors of land plants. We discover an early evolution of SAL1-PAP chloroplast retrograde signaling in stomatal regulation based on conserved gene and protein structure, function, and enzyme activity and transit peptides of SAL1s in species including flowering plants, the fern Ceratopteris richardii, and the moss Physcomitrella patens Moreover, we demonstrate that PAP regulates stomatal closure via secondary messengers and ion transport in guard cells of these diverse lineages. The origin of stomata facilitated gas exchange in the earliest land plants. Our findings suggest that the conquest of land by plants was enabled by rapid response to drought stress through the deployment of an ancestral SAL1-PAP signaling pathway, intersecting with the core abscisic acid signaling in stomatal guard cells.

Plant 24-nt reproductive phasiRNAs from intramolecular duplex mRNAs in diverse monocots.

In grasses, two pathways that generate diverse and numerous 21-nt (premeiotic) and 24-nt (meiotic) phased siRNAs are highly enriched in anthers, the male reproductive organs. These "phasiRNAs" are analogous to mammalian piRNAs, yet their functions and evolutionary origins remain largely unknown. The 24-nt meiotic phasiRNAs have only been described in grasses, wherein their biogenesis is dependent on a specialized Dicer (DCL5). To assess how evolution gave rise to this pathway, we examined reproductive phasiRNA pathways in nongrass monocots: garden asparagus, daylily, and lily. The common ancestors of these species diverged approximately 115-117 million years ago (MYA). We found that premeiotic 21-nt and meiotic 24-nt phasiRNAs were abundant in all three species and displayed spatial localization and temporal dynamics similar to grasses. The miR2275-triggered pathway was also present, yielding 24-nt reproductive phasiRNAs, and thus originated more than 117 MYA. In asparagus, unlike in grasses, these siRNAs are largely derived from inverted repeats (IRs); analyses in lily identified thousands of precursor loci, and many were also predicted to form foldback substrates for Dicer processing. Additionally, reproductive phasiRNAs were present in female reproductive organs and thus may function in both male and female germinal development. These data describe several distinct mechanisms of production for 24-nt meiotic phasiRNAs and provide new insights into the evolution of reproductive phasiRNA pathways in monocots.

Evolution of rapid blue-light response linked to explosive diversification of ferns in angiosperm forests.

Ferns appear in the fossil record some 200 Myr before angiosperms. However, as angiosperm-dominated forest canopies emerged in the Cretaceous period there was an explosive diversification of modern (leptosporangiate) ferns, which thrived in low, blue-enhanced light beneath angiosperm canopies. A mechanistic explanation for this transformative event in the diversification of ferns has remained elusive. We used physiological assays, transcriptome analysis and evolutionary bioinformatics to investigate a potential connection between the evolution of enhanced stomatal sensitivity to blue light in modern ferns and the rise of angiosperm-dominated forests in the geological record. We demonstrate that members of the largest subclade of leptosporangiate ferns, Polypodiales, have significantly faster stomatal response to blue light than more ancient fern lineages and a representative angiosperm. We link this higher sensitivity to levels of differentially expressed genes in blue-light signaling, particularly in the cryptochrome (CRY) signaling pathway. Moreover, CRYs of the Polypodiales examined show gene duplication events between 212.9-196.9 and 164.4-151.8 Ma, when angiosperms were emerging, which are lacking in other major clades of extant land plants. These findings suggest that evolution of stomatal blue-light sensitivity helped modern ferns exploit the shady habitat beneath angiosperm forest canopies, fueling their Cretaceous hyperdiversification.

Comprehending Cornales: phylogenetic reconstruction of the order using the Angiosperms353 probe set.

PREMISE: Cornales is an order of flowering plants containing ecologically and horticulturally important families, including Cornaceae (dogwoods) and Hydrangeaceae (hydrangeas), among others. While many relationships in Cornales are strongly supported by previous studies, some uncertainty remains with regards to the placement of Hydrostachyaceae and to relationships among families in Cornales and within Cornaceae. Here we analyzed hundreds of nuclear loci to test published phylogenetic hypotheses and estimated a robust species tree for Cornales. METHODS: Using the Angiosperms353 probe set and existing data sets, we generated phylogenomic data for 158 samples, representing all families in the Cornales, with intensive sampling in the Cornaceae. RESULTS: We curated an average of 312 genes per sample, constructed maximum likelihood gene trees, and inferred a species tree using the summary approach implemented in ASTRAL-III, a method statistically consistent with the multispecies coalescent model. CONCLUSIONS: The species tree we constructed generally shows high support values and a high degree of concordance among individual nuclear gene trees. Relationships among families are largely congruent with previous molecular studies, except for the placement of the nyssoids and the Grubbiaceae-Curtisiaceae clades. Furthermore, we were able to place Hydrostachyaceae within Cornales, and within Cornaceae, the monophyly of known morphogroups was well supported. However, patterns of gene tree discordance suggest potential ancient reticulation, gene flow, and/or ILS in the Hydrostachyaceae lineage and the early diversification of Cornus. Our findings reveal new insights into the diversification process across Cornales and demonstrate the utility of the Angiosperms353 probe set.

A Universal Probe Set for Targeted Sequencing of 353 Nuclear Genes from Any Flowering Plant Designed Using k-Medoids Clustering.

Sequencing of target-enriched libraries is an efficient and cost-effective method for obtaining DNA sequence data from hundreds of nuclear loci for phylogeny reconstruction. Much of the cost of developing targeted sequencing approaches is associated with the generation of preliminary data needed for the identification of orthologous loci for probe design. In plants, identifying orthologous loci has proven difficult due to a large number of whole-genome duplication events, especially in the angiosperms (flowering plants). We used multiple sequence alignments from over 600 angiosperms for 353 putatively single-copy protein-coding genes identified by the One Thousand Plant Transcriptomes Initiative to design a set of targeted sequencing probes for phylogenetic studies of any angiosperm group. To maximize the phylogenetic potential of the probes, while minimizing the cost of production, we introduce a k-medoids clustering approach to identify the minimum number of sequences necessary to represent each coding sequence in the final probe set. Using this method, 5-15 representative sequences were selected per orthologous locus, representing the sequence diversity of angiosperms more efficiently than if probes were designed using available sequenced genomes alone. To test our approximately 80,000 probes, we hybridized libraries from 42 species spanning all higher-order groups of angiosperms, with a focus on taxa not present in the sequence alignments used to design the probes. Out of a possible 353 coding sequences, we recovered an average of 283 per species and at least 100 in all species. Differences among taxa in sequence recovery could not be explained by relatedness to the representative taxa selected for probe design, suggesting that there is no phylogenetic bias in the probe set. Our probe set, which targeted 260 kbp of coding sequence, achieved a median recovery of 137 kbp per taxon in coding regions, a maximum recovery of 250 kbp, and an additional median of 212 kbp per taxon in flanking non-coding regions across all species. These results suggest that the Angiosperms353 probe set described here is effective for any group of flowering plants and would be useful for phylogenetic studies from the species level to higher-order groups, including the entire angiosperm clade itself.

Horizontal Transfer of Non-LTR Retrotransposons from Arthropods to Flowering Plants.

Even though lateral movements of transposons across families and even phyla within multicellular eukaryotic kingdoms have been found, little is known about transposon transfer between the kingdoms Animalia and Plantae. We discovered a novel non-LTR retrotransposon, AdLINE3, in a wild peanut species. Sequence comparisons and phylogenetic analyses indicated that AdLINE3 is a member of the RTE clade, originally identified in a nematode and rarely reported in plants. We identified RTE elements in 82 plants, spanning angiosperms to algae, including recently active elements in some flowering plants. RTE elements in flowering plants were likely derived from a single family we refer to as An-RTE. Interestingly, An-RTEs show significant DNA sequence identity with non-LTR retroelements from 42 animals belonging to four phyla. Moreover, the sequence identity of RTEs between two arthropods and two plants was higher than that of homologous genes. Phylogenetic and evolutionary analyses of RTEs from both animals and plants suggest that the An-RTE family was likely transferred horizontally into angiosperms from an ancient aphid(s) or ancestral arthropod(s). Notably, some An-RTEs were recruited as coding sequences of functional genes participating in metabolic or other biochemical processes in plants. This is the first potential example of horizontal transfer of transposons between animals and flowering plants. Our findings help to understand exchanges of genetic material between the kingdom Animalia and Plantae and suggest arthropods likely impacted on plant genome evolution.

Concomitant loss of NDH complex-related genes within chloroplast and nuclear genomes in some orchids.

The chloroplast NAD(P)H dehydrogenase-like (NDH) complex consists of about 30 subunits from both the nuclear and chloroplast genomes and is ubiquitous across most land plants. In some orchids, such as Phalaenopsis equestris, Dendrobium officinale and Dendrobium catenatum, most of the 11 chloroplast genome-encoded ndh genes (cp-ndh) have been lost. Here we investigated whether functional cp-ndh genes have been completely lost in these orchids or whether they have been transferred and retained in the nuclear genome. Further, we assessed whether both cp-ndh genes and nucleus-encoded NDH-related genes can be lost, resulting in the absence of the NDH complex. Comparative analyses of the genome of Apostasia odorata, an orchid species with a complete complement of cp-ndh genes which represents the sister lineage to all other orchids, and three published orchid genome sequences for P. equestris, D. officinale and D. catenatum, which are all missing cp-ndh genes, indicated that copies of cp-ndh genes are not present in any of these four nuclear genomes. This observation suggests that the NDH complex is not necessary for some plants. Comparative genomic/transcriptomic analyses of currently available plastid genome sequences and nuclear transcriptome data showed that 47 out of 660 photoautotrophic plants and all the heterotrophic plants are missing plastid-encoded cp-ndh genes and exhibit no evidence for maintenance of a functional NDH complex. Our data indicate that the NDH complex can be lost in photoautotrophic plant species. Further, the loss of the NDH complex may increase the probability of transition from a photoautotrophic to a heterotrophic life history.

Parallel evolution of storage roots in morning glories (Convolvulaceae).

BACKGROUND: Storage roots are an ecologically and agriculturally important plant trait that have evolved numerous times in angiosperms. Storage roots primarily function to store carbohydrates underground as reserves for perennial species. In morning glories, storage roots are well characterized in the crop species sweetpotato, where starch accumulates in storage roots. This starch-storage tissue proliferates, and roots thicken to accommodate the additional tissue. In morning glories, storage roots have evolved numerous times. The primary goal of this study is to understand whether this was through parallel evolution, where species use a common genetic mechanism to achieve storage root formation, or through convergent evolution, where storage roots in distantly related species are formed using a different set of genes. Pairs of species where one forms storage roots and the other does not were sampled from two tribes in the morning glory family, the Ipomoeeae and Merremieae. Root anatomy in storage roots and fine roots was examined. Furthermore, we sequenced total mRNA from storage roots and fine roots in these species and analyzed differential gene expression. RESULTS: Anatomical results reveal that storage roots of species in the Ipomoeeae tribe, such as sweetpotato, accumulate starch similar to species in the Merremieae tribe but differ in vascular tissue organization. In both storage root forming species, more genes were found to be upregulated in storage roots compared to fine roots. Further, we find that fifty-seven orthologous genes were differentially expressed between storage roots and fine roots in both storage root forming species. These genes are primarily involved in starch biosynthesis, regulation of starch biosynthesis, and transcription factor activity. CONCLUSIONS: Taken together, these results demonstrate that storage roots of species from both morning glory tribes are anatomically different but utilize a common core set of genes in storage root formation. This is consistent with a pattern of parallel evolution, thus highlighting the importance of examining anatomy together with gene expression to understand the evolutionary origins of ecologically and economically important plant traits.

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.

Shifts in gene expression profiles are associated with weak and strong Crassulacean acid metabolism.

PREMISE OF THE STUDY: The relative ease of high throughput sequencing is facilitating comprehensive phylogenomic and gene expression studies, even for nonmodel groups. To date, however, these two approaches have not been merged; while phylogenomic methods might use transcriptome sequences to resolve relationships, assessment of gene expression patterns in a phylogenetic context is less common. Here we analyzed both carbon assimilation and gene expression patterns of closely related species within the Agavoideae (Asparagaceae) to elucidate changes in gene expression across weak and strong phenotypes for Crassulacean acid metabolism (CAM). METHODS: Gene expression patterns were compared across four genera: Agave (CAM), which is paraphyletic with Polianthes (weak CAM) and Manfreda (CAM), and Beschorneria (weak CAM). RNA-sequencing was paired with measures of gas exchange and titratable acidity. Climate niche space was compared across the four lineages to examine abiotic factors and their correlation to CAM. KEY RESULTS: Expression of homologous genes showed both shared and variable patterns in weak and strong CAM species. Network analysis highlights that despite shared expression patterns, highly connected genes differ between weak and strong CAM, implicating shifts in regulatory gene function as key for the evolution of CAM. Variation in carbohydrate metabolism between weak and strong CAM supports the importance of sugar turnovers for CAM physiology. CONCLUSIONS: Integration of phylogenetics and RNA-sequencing provides a powerful tool to study the evolution of CAM photosynthesis across closely related but photosynthetically variable species. Our findings regarding shared or shifted gene expression and regulation of CAM via carbohydrate metabolism have important implications for efforts to engineer the CAM pathway into C3 food and biofuel crops.

Phylotranscriptomic analysis and genome evolution of the Cypripedioideae (Orchidaceae).

PREMISE OF THE STUDY: The slipper orchids (Cypripedioideae) are a morphologically distinct subfamily of Orchidaceae. They also have some of the largest genomes in the orchids, which may be due to polyploidy or some other mechanism of genome evolution. We generated 10 transcriptomes and incorporated existing RNA-seq data to infer a multilocus nuclear phylogeny of the Cypripedioideae and to determine whether a whole-genome duplication event (WGD) correlated with the large genome size of this subfamily. Knowing more about timing of ancient polyploidy events can help us understand the evolution of one of the most species-rich plant families. METHODS: Transcriptome data were used to identify low-copy orthologous genes to infer a phylogeny of Orchidaceae and to identify paralogs to place any WGD events on the species tree. KEY RESULTS: Our transcriptome phylogeny confirmed relationships published in previous studies that used fewer markers but incorporated more taxa. We did not find a WGD event at the base of the slipper orchids; however, we did identify one on the Orchidaceae stem lineage. We also confirmed the presence of a previously identified WGD event deeper in the monocot phylogeny. CONCLUSIONS: Although WGD has played a role in the evolution of Orchidaceae, polyploidy does not appear to be responsible for the large genome size of slipper orchids. The conserved set of 775 largely single-copy nuclear genes identified in this study should prove useful in future studies of orchid evolution.

Brassicales phylogeny inferred from 72 plastid genes: A reanalysis of the phylogenetic localization of two paleopolyploid events and origin of novel chemical defenses.

PREMISE OF THE STUDY: Previous phylogenetic studies employing molecular markers have yielded various insights into the evolutionary history across Brassicales, but many relationships between families remain poorly supported or unresolved. A recent phylotranscriptomic approach utilizing 1155 nuclear markers obtained robust estimates for relationships among 14 of 17 families. Here we report a complete family-level phylogeny estimated using the plastid genome. METHODS: We conducted phylogenetic analyses on a concatenated data set comprising 44,926 bp from 72 plastid genes for species distributed across all 17 families. Our analysis includes three additional families, Tovariaceae, Salvadoraceae, and Setchellanthaceae, that were omitted in the previous phylotranscriptomic study. KEY RESULTS: Our phylogenetic analyses obtained fully resolved and strongly supported estimates for all nodes across Brassicales. Importantly, these findings are congruent with the topology reported in the phylotranscriptomic study. This consistency suggests that future studies could utilize plastid genomes as markers for resolving relationships within some notoriously difficult clades across Brassicales. We used this new phylogenetic framework to verify the placement of the At-α event near the origin of Brassicaceae, with median date estimates of 31.8 to 42.8 million years ago and restrict the At-ß event to one of two nodes with median date estimates between 85 to 92.2 million years ago. These events ultimately gave rise to novel chemical defenses and are associated with subsequent shifts in net diversification rates. CONCLUSIONS: We anticipate that these findings will aid future comparative evolutionary studies across Brassicales, including selecting candidates for whole-genome sequencing projects.

Monocot plastid phylogenomics, timeline, net rates of species diversification, the power of multi-gene analyses, and a functional model for the origin of monocots.

PREMISE OF THE STUDY: We present the first plastome phylogeny encompassing all 77 monocot families, estimate branch support, and infer monocot-wide divergence times and rates of species diversification. METHODS: We conducted maximum likelihood analyses of phylogeny and BAMM studies of diversification rates based on 77 plastid genes across 545 monocots and 22 outgroups. We quantified how branch support and ascertainment vary with gene number, branch length, and branch depth. KEY RESULTS: Phylogenomic analyses shift the placement of 16 families in relation to earlier studies based on four plastid genes, add seven families, date the divergence between monocots and eudicots+Ceratophyllum at 136 Mya, successfully place all mycoheterotrophic taxa examined, and support recognizing Taccaceae and Thismiaceae as separate families and Arecales and Dasypogonales as separate orders. Only 45% of interfamilial divergences occurred after the Cretaceous. Net species diversification underwent four large-scale accelerations in PACMAD-BOP Poaceae, Asparagales sister to Doryanthaceae, Orchidoideae-Epidendroideae, and Araceae sister to Lemnoideae, each associated with specific ecological/morphological shifts. Branch ascertainment and support across monocots increase with gene number and branch length, and decrease with relative branch depth. Analysis of entire plastomes in Zingiberales quantifies the importance of non-coding regions in identifying and supporting short, deep branches. CONCLUSIONS: We provide the first resolved, well-supported monocot phylogeny and timeline spanning all families, and quantify the significant contribution of plastome-scale data to resolving short, deep branches. We outline a new functional model for the evolution of monocots and their diagnostic morphological traits from submersed aquatic ancestors, supported by convergent evolution of many of these traits in aquatic Hydatellaceae (Nymphaeales).

Phylotranscriptomic analysis of the origin and early diversification of land plants.

Reconstructing the origin and evolution of land plants and their algal relatives is a fundamental problem in plant phylogenetics, and is essential for understanding how critical adaptations arose, including the embryo, vascular tissue, seeds, and flowers. Despite advances in molecular systematics, some hypotheses of relationships remain weakly resolved. Inferring deep phylogenies with bouts of rapid diversification can be problematic; however, genome-scale data should significantly increase the number of informative characters for analyses. Recent phylogenomic reconstructions focused on the major divergences of plants have resulted in promising but inconsistent results. One limitation is sparse taxon sampling, likely resulting from the difficulty and cost of data generation. To address this limitation, transcriptome data for 92 streptophyte taxa were generated and analyzed along with 11 published plant genome sequences. Phylogenetic reconstructions were conducted using up to 852 nuclear genes and 1,701,170 aligned sites. Sixty-nine analyses were performed to test the robustness of phylogenetic inferences to permutations of the data matrix or to phylogenetic method, including supermatrix, supertree, and coalescent-based approaches, maximum-likelihood and Bayesian methods, partitioned and unpartitioned analyses, and amino acid versus DNA alignments. Among other results, we find robust support for a sister-group relationship between land plants and one group of streptophyte green algae, the Zygnematophyceae. Strong and robust support for a clade comprising liverworts and mosses is inconsistent with a widely accepted view of early land plant evolution, and suggests that phylogenetic hypotheses used to understand the evolution of fundamental plant traits should be reevaluated.