Search

Pogoniopsis likely represents an independent photosynthesis loss in orchids. We use phylogenomic data to better identify the phylogenetic placement of this fully mycoheterotrophic taxon, and investigate its molecular evolution.We performed likelihood analysis of plastid and mitochondrial phylogenomic data to localize the position of Pogoniopsis schenckii in orchid phylogeny, and investigated the evolution of its plastid genome.All analyses place Pogoniopsis in subfamily Epidendroideae, with strongest support from mitochondrial data, which also place it near tribe Sobralieae with moderately strong support. Extreme rate elevation in Pogoniopsis plastid genes broadly depresses branch support; in contrast, mitochondrial genes are only mildly rate elevated and display very modest and localized reductions in bootstrap support. Despite considerable genome reduction, including loss of photosynthesis genes and multiple translation apparatus genes, gene order in Pogoniopsis plastomes is identical to related autotrophs, apart from moderately shifted inverted repeat (IR) boundaries. All cis-spliced introns have been lost in retained genes. Two plastid genes (accD, rpl2) show significant strengthening of purifying selection. A retained plastid tRNA gene (trnE-UUC) of Pogoniopsis lacks an anticodon; we predict that it no longer functions in translation but retains a secondary role in heme biosynthesis.Slowly evolving mitochondrial genes clarify the placement of Pogoniopsis in orchid phylogeny, a strong contrast with analysis of rate-elevated plastome data. We documented the effects of the novel loss of photosynthesis: for example, despite massive gene loss, its plastome is fully colinear with other orchids, and it displays only moderate shifts in selective pressure in retained genes.

SummaryAfrothismiais a genus of non-photosynthetic mycoheterotrophs from the forests of continental tropical Africa. Multiple phylogenetic inferences using molecular data recover the genus as sister to a clade comprising mycoheterotrophic Thismiaceae and the photosynthetic family Taccaceae, contrary to earlier placements ofAfrothismiaand Thismiaceae within Burmanniaceae. Morphological support for separatingAfrothismiafrom the rest of Thismiaceae has depended on the zygomorphic flowers ofAfrothismia(although some South American species ofThismiaare also zygomorphic) and their clusters of root tubers, each with a terminal rootlet. The number of described species ofAfrothismiahas recently increased substantially, from four to 16, which has provided additional morphological characters that support its distinction from Thismiaceae. Most notably, the ovary inAfrothismiahas a single stalked placenta, and circumscissile fruits from which seeds are exserted by placental elevation (in Thismiaceae, in contrast, there are three placentas, a deliquescing fruit lid, and the seeds are not exserted).Afrothismiastamens are inserted in the lower perianth tube where they are attached to the stigma, and individual flowers are subtended by a single large dorsal bract (in Thismiaceae, stamens are inserted at the mouth of the tube, free of and distant from the stigma, and each flower is subtended by a loose whorl of (2-)3(−4) bracts). Here we formally characterise Afrothismiaceae and review what is known of its development, seed germination, interactions with mycorrhizal Glomeromycota, biogeography, phylogeny and pollination biology. All but one (Afrothismia insignis; Vulnerable) of the 13 species assessed on the IUCN Red List of Threatened Species are either Endangered or Critically Endangered; one species (A. pachyanthaSchltr.) is considered to be extinct.

Mycoheterotrophy is an alternative nutritional strategy whereby plants obtain sugars and other nutrients from soil fungi. Mycoheterotrophy and associated loss of photosynthesis have evolved repeatedly in plants, particularly in monocots. Although reductive evolution of plastomes in mycoheterotrophs is well documented, the dynamics of nuclear genome evolution remains largely unknown. Transcriptome datasets were generated from four mycoheterotrophs in three families (Orchidaceae, Burmanniaceae, Triuridaceae) and related green plants and used for phylogenomic analyses to resolve relationships among the mycoheterotrophs, their relatives, and representatives across the monocots. Phylogenetic trees based on 602 genes were mostly congruent with plastome phylogenies, except for an Asparagales + Liliales clade inferred in the nuclear trees. Reduction and loss of chlorophyll synthesis and photosynthetic gene expression and relaxation of purifying selection on retained genes were progressive, with greater loss in older nonphotosynthetic lineages. One hundred seventy-four of 1375 plant benchmark universally conserved orthologous genes were undetected in any mycoheterotroph transcriptome or the genome of the mycoheterotrophic orchid Gastrodia but were expressed in green relatives, providing evidence for massively convergent gene loss in nonphotosynthetic lineages. We designate this set of deleted or undetected genes Missing in Mycoheterotrophs (MIM). MIM genes encode not only mainly photosynthetic or plastid membrane proteins but also a diverse set of plastid processes, genes of unknown function, mitochondrial, and cellular processes. Transcription of a photosystem II gene (psb29) in all lineages implies a nonphotosynthetic function for this and other genes retained in mycoheterotrophs. Nonphotosynthetic plants enable novel insights into gene function as well as gene expression shifts, gene loss, and convergence in nuclear genomes.

Inferring accurate biogeographic history of plant taxa with an East Asia (EA)-North America (NA) is usually hindered by conflicting phylogenies and a poor fossil record. The current distribution of Chamaecyparis (false cypress; Cupressaceae) with four species in EA, and one each in western and eastern NA, and its relatively rich fossil record, make it an excellent model for studying the EA-NA disjunction. Here we reconstruct phylogenomic relationships within Chamaecyparis using > 1400 homologous nuclear and 61 plastid genes. Our phylogenomic analyses using concatenated and coalescent approaches revealed strong cytonuclear discordance and conflicting topologies between nuclear gene trees. Incomplete lineage sorting (ILS) and hybridization are possible explanations of conflict; however, our coalescent analyses and simulations suggest that ILS is the major contributor to the observed phylogenetic discrepancies. Based on a well-resolved species tree and four fossil calibrations, the crown lineage of Chamaecyparis is estimated to have originated in the upper Cretaceous, followed by diversification events in the early and middle Paleogene. Ancestral area reconstructions suggest that Chamaecyparis had an ancestral range spanning both EA and NA. Fossil records further indicate that this genus is a relict of the “boreotropical” flora, and that local extinctions of European species were caused by global cooling. Overall, our results unravel a complex evolutionary history of a Paleogene relict conifer genus, which may have involved ILS, hybridization and the extinction of local species.

Fully mycoheterotrophic plants can be difficult to place in plant phylogeny due to elevated substitution rates associated with photosynthesis loss. This potentially limits the effectiveness of downstream analyses of mycoheterotrophy that depend on accurate phylogenetic inference. Although mitochondrial genomic data sets are rarely used in plant phylogenetics, theory predicts that they should be resilient to long-branch artefacts, thanks to their generally slow evolution, coupled with limited rate elevation in heterotrophs. We examined the utility of mitochondrial genomes for resolving contentious higher-order placements of mycoheterotrophic lineages in two test cases: monocots (focusing on Dioscoreales) and Ericaceae. We find Thismiaceae to be distantly related to Burmanniaceae in the monocot order Dioscoreales, conflicting with current classification schemes based on few gene data sets. We confirm that the unusual Afrothismia is related to Taccaceae-Thismiaceae, with a corresponding independent loss of photosynthesis. In Ericaceae we recovered the first well supported relationships among its five major lineages: mycoheterotrophic Ericaceae are not monophyletic, as pyroloids are inferred to be sister to core Ericaceae, and monotropoids to arbutoids. Genes recovered from mitochondrial genomes collectively resolved previously ambiguous mycoheterotroph higher-order relationships. We propose that mitochondrial genomic data should be considered in standardised gene panels for inferring overall plant phylogeny.

Carnivorous plants consume animals for mineral nutrients that enhance growth and reproduction in nutrient-poor environments. Here, we report that Triantha occidentalis (Tofieldiaceae) represents a previously overlooked carnivorous lineage that captures insects on sticky inflorescences. Field experiments, isotopic data, and mixing models demonstrate significant N transfer from prey to Triantha, with an estimated 64% of leaf N obtained from prey capture in previous years, comparable to levels inferred for the cooccurring round-leaved sundew, a recognized carnivore. N obtained via carnivory is exported from the inflorescence and developing fruits and may ultimately be transferred to next year’s leaves. Glandular hairs on flowering stems secrete phosphatase, as seen in all carnivorous plants that directly digest prey. Triantha is unique among carnivorous plants in capturing prey solely with sticky traps adjacent to its flowers, contrary to theory. However, its glandular hairs capture only small insects, unlike the large bees and butterflies that act as pollinators, which may minimize the conflict between carnivory and pollination.

5 Pág. Centro de Biotecnología y Genómica de Plantas

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.

Background and Aims Fully mycoheterotrophic plants derive carbon and other nutrients from root-associated fungi and have lost the ability to photosynthesize. While mycoheterotroph plastomes are often degraded compared with green plants, the effect of this unusual symbiosis on mitochondrial genome evolution is unknown. By providing the first complete organelle genome data from Polygalaceae, one of only three eudicot families that developed mycoheterotrophy, we explore how both organellar genomes evolved after loss of photosynthesis. Methods We sequenced and assembled four complete plastid genomes and a mitochondrial genome from species of Polygalaceae, focusing on non-photosynthetic Epirixanthes. We compared these genomes with those of other mycoheterotroph and parasitic plant lineages, and assessed whether organelle genes in Epirixanthes experienced relaxed or intensified selection compared with autotrophic relatives. Key Results Plastomes of two species of Epirixanthes have become substantially degraded compared with that of autotrophic Polygala. Although the lack of photosynthesis is presumably homologous in the genus, the surveyed Epirixanthes species have marked differences in terms of plastome size, structural rearrangements, gene content and substitution rates. Remarkably, both apparently replaced a canonical plastid inverted repeat with large directly repeated sequences. The mitogenome of E. elongata incorporated a considerable number of fossilized plastid genes, by intracellular transfer from an ancestor with a less degraded plastome. Both plastid and mitochondrial genes in E. elongata have increased substitution rates, but the plastid genes of E. pallida do not. Despite this, both species have similar selection patterns operating on plastid housekeeping genes. Conclusions Plastome evolution largely fits with patterns of gene degradation seen in other heterotrophic plants, but includes highly unusual directly duplicated regions. The causes of rate elevation in the sequenced Epirixanthes mitogenome and of rate differences in plastomes of related mycoheterotrophic species are not currently understood.

Angiosperm lineages in aquatic environments are characterized by high structural and functional diversity, and wide distributions. A long-standing evolutionary riddle is what processes have caused the relatively low diversity of aquatic angiosperms compared to their terrestrial relatives. We use diversification and ancestral reconstruction models with a comprehensive > 10 000 genus angiosperm phylogeny to elucidate the macroevolutionary dynamics associated with transitions of terrestrial plants to water. Our study reveals that net diversification rates are significantly lower in aquatic than in terrestrial angiosperms due to lower speciation and higher extinction. Shifts from land to water started early in angiosperm evolution, but most events were concentrated during the last c. 25 million years. Reversals to a terrestrial habitat started only 40 million years ago, but occurred at much higher rates. Within aquatic angiosperms, the estimated pattern is one of gradual accumulation of lineages, and relatively low and constant diversification rates throughout the Cenozoic. Low diversification rates, together with infrequent water transitions, account for the low diversity of aquatic angiosperms today. The stressful conditions and small global surface of the aquatic habitat available for angiosperms are hypothesized to explain this pattern.

Organellar genomes serve as useful models for genome evolution and contain some of the most widely used phylogenetic markers, but they are poorly characterized in many lineages. Here, we report 20 novel mitochondrial genomes and 16 novel plastid genomes from the brown algae. We focused our efforts on the orders Chordales and Laminariales but also provide the first plastid genomes (plastomes) from Desmarestiales and Sphacelariales, the first mitochondrial genome (mitome) from Ralfsiales and a nearly complete mitome from Sphacelariales. We then compared gene content, sequence evolution rates, shifts in genome structural arrangements, and intron distributions across lineages. We confirm that gene content is largely conserved in both organellar genomes across the brown algal tree of life, with few cases of gene gain or loss. We further show that substitution rates are generally lower in plastid than mitochondrial genes, but plastomes are more variable in gene arrangement, as mitomes tend to be colinear even among distantly related lineages (with exceptions). Patterns of intron distribution across organellar genomes are complex. In particular, the mitomes of several laminarialean species possess group II introns that have T7-like ORFs, found previously only in mitochondrial genomes of Pylaiella spp. (Ectocarpales). The distribution of these mitochondrial introns is inconsistent with vertical transmission and likely reflects invasion by horizontal gene transfer between lineages. In the most extreme case, the mitome of Hedophyllum nigripes is ∼40% larger than the mitomes of close relatives because of these introns. Our results provide substantial insight into organellar evolution across the brown algae.

We used a bi-organellar phylogenomic approach to address higher-order relationships in Pandanales, including the first molecular phylogenetic study of the panama-hat family, Cyclanthaceae. Our genus-level study of plastid and mitochondrial gene sets includes a comprehensive sampling of photosynthetic lineages across the order, and provides a framework for investigating clade ages, biogeographic hypotheses and organellar molecular evolution. Using multiple inference methods and both organellar genomes, we recovered mostly congruent and strongly supported relationships within and between families, including the placement of fully mycoheterotrophic Triuridaceae. Cyclanthaceae and Pandanaceae plastomes have slow substitution rates, contributing to weakly supported plastid-based relationships in Cyclanthaceae. While generally slowly evolving, mitochondrial genomes exhibit sporadic rate elevation across the order. However, we infer well-supported relationships even for slower evolving mitochondrial lineages in Cyclanthaceae. Clade age estimates across photosynthetic lineages are largely consistent with previous studies, are well correlated between the two organellar genomes (with slightly younger inferences from mitochondrial data), and support several biogeographic hypotheses. We show that rapidly evolving non-photosynthetic lineages may bias age estimates upwards at neighbouring photosynthetic nodes, even using a relaxed clock model. Finally, we uncovered new genome structural variants in photosynthetic taxa at plastid inverted repeat boundaries that show promise as interfamilial phylogenetic markers.

Carbohydrate-active enzymes (CAZymes) are central to the biosynthesis and modification of the plant cell wall. An ancient clade of bifunctional plant endo-glucanases (EG16 members) was recently revealed and proposed to represent a transitional group uniting plant xyloglucan endo-transglycosylase/hydrolase (XTH) gene products and bacterial mixed-linkage endo-glucanases in the phylogeny of glycoside hydrolase family 16 (GH16). To gain broader insights into the distribution and frequency of EG16 and other GH16 members in plants, the PHYTOZOME, PLAZA, NCBI and 1000 PLANTS databases were mined to build a comprehensive census among 1289 species, spanning the broad phylogenetic diversity of multiple algae through recent plant lineages. EG16, newly identified EG16-2 and XTH members appeared first in the green algae. Extant EG16 members represent the early adoption of the β-jellyroll protein scaffold from a bacterial or early-lineage eukaryotic GH16 gene, which is characterized by loop deletion and extension of the N terminus (in EG16-2 members) or C terminus (in XTH members). Maximum-likelihood phylogenetic analysis of EG16 and EG16-2 sequences are directly concordant with contemporary estimates of plant evolution. The lack of expansion of EG16 members into multi-gene families across green plants may point to a core metabolic role under tight control, in contrast to XTH genes that have undergone the extensive duplications typical of cell-wall CAZymes. The present census will underpin future studies to elucidate the physiological role of EG16 members across plant species, and serve as roadmap for delineating the closely related EG16 and XTH gene products in bioinformatic analyses of emerging genomes and transcriptomes.

Premise of the study Phylogenomic studies employing large numbers of genes, including those based on plastid genomes (plastomes), are becoming common. Nonphotosynthetic plants such as mycoheterotrophs (which rely on root-associated fungi for essential nutrients, including carbon) tend to have highly elevated rates of plastome evolution, substantial genome reduction, or both. Mycoheterotroph plastomes therefore provide excellent test cases for investigating how extreme conditions impact phylogenomic inference. Methods We used parsimony and likelihood analysis of protein-coding gene sets from published and newly completed plastomes to infer the phylogenetic placement of taxa from the 10 angiosperm families in which mycoheterotrophy evolved. Key results Despite multiple very long branches that reflect elevated substitution rates, and frequently patchy gene recovery due to genome reduction, inferred phylogenetic placements of most mycoheterotrophic lineages in DNA-based likelihood analyses are both well supported and congruent with other studies. Amino-acid-based likelihood placements are broadly consistent with DNA-based inferences, but extremely rate-elevated taxa can have unexpected placements-albeit with weak support. In contrast, parsimony analysis is strongly misled by long-branch attraction among many distantly related mycoheterotrophic monocots. Conclusions Mycoheterotrophic plastomes provide challenging cases for phylogenomic inference, as substitutional rates can be elevated and genome reduction can lead to sparse gene recovery. Nonetheless, diverse likelihood frameworks provide generally well-supported and mutually concordant phylogenetic placements of mycoheterotrophs, consistent with recent phylogenetic studies and angiosperm-wide classifications. Previous predictions of parallel photosynthesis loss within families are supported for Burmanniaceae, Ericaceae, Gentianaceae, and Orchidaceae. Burmanniaceae and Thismiaceae should not be combined as a single family in Dioscoreales.

The 1,000 Plants (1KP) initiative was the first large-scale effort to collect next-generation sequencing (NGS) data across a phylogenetically representative sampling of species for a major clade of life, in this case the Viridiplantae, or green plants. As an international multidisciplinary consortium, we focused on plant evolution and its practical implications. Among the major outcomes were the inference of a reference species tree for green plants by phylotranscriptomic analysis of low-copy genes, a survey of paleopolyploidy (whole-genome duplications) across the Viridiplantae, the inferred evolutionary histories for many gene families and biological processes, the discovery of novel light-sensitive proteins for optogenetic studies in mammalian neuroscience, and elucidation of the genetic network for a complex trait (C4photosynthesis). Altogether, 1KP demonstrated how value can be extracted from a phylodiverse sequencing data set, providing a template for future projects that aim to generate even more data, including complete de novo genomes, across the tree of life.

In this short paper, considering the significant increase on the number of electric vehicles (EV) and the additional load they bring on the power grid, we propose a static charging scheduling algorithm in order to effectively schedule the charging of EVs during nighttime at private charging stations. Given the vehicles and the volume of their charging demands, we first determine the additional load to the grid assuming all EVs need to charge to the maximum level recommended by the manufacturers within the given time period. Then, the proposed algorithm calculates the energy and the time needed to charge the EV battery to the level requested by the owner. After finding the slack time (the latest start time for charging – the current time), it schedules the charging of EV by giving priority to the one whose slack time is minimum. If it is not possible to charge all vehicles to the requested battery level within the time constraints, the algorithm evenly distributes the available energy to EVs that need more charging.

Premise We developed a target enrichment panel for phylogenomic studies of Dioscorea, an economically important genus with incompletely resolved relationships. Methods Our bait panel comprises 260 low- to single-copy nuclear genes targeted to work in Dioscorea, assessed here using a preliminary taxon sampling that includes both distantly and closely related taxa, including several yam crops and potential crop wild relatives. We applied coalescent-based and maximum likelihood phylogenomic inference approaches to the pilot taxon set, incorporating new and published transcriptome data from additional species. Results The custom panel retrieved ~94% of targets and >80% of full gene length from 88% and 68% of samples, respectively. It has minimal gene overlap with existing panels designed for angiosperm-wide studies and generally recovers longer and more variable targets. Pilot phylogenomic analyses consistently resolve most deep and recent relationships with strong support across analyses and point to revised relationships between the crop species D. alata and candidate crop wild relatives. Discussion Our customized panel reliably retrieves targeted loci from Dioscorea, is informative for resolving relationships in denser samplings, and is suitable for refining our understanding of the independent origins of cultivated yam species; the panel likely has broader promise for phylogenomic studies across Dioscoreales.

Reconstructing phylogenetic topologies and divergence times is essential for inferring the timing of radiations, the appearance of adaptations, and the historical biogeography of key lineages. In temperate marine ecosystems, kelps (Laminariales) drive productivity and form essential habitat but an incomplete understanding of their phylogeny has limited our ability to infer their evolutionary origins and the spatial and temporal patterns of their diversification. Here, we reconstruct the diversification of habitat-forming kelps using a global genus-level phylogeny inferred primarily from organellar genome datasets, and investigate the timing of kelp radiation. We resolve several important phylogenetic features, including relationships among the morphologically simple kelp families and the broader radiation of complex kelps, demonstrating that the initial radiation of the latter resulted from an increase in speciation rate around the Eocene-Oligocene boundary. This burst in speciation rate is consistent with a possible role of recent climatic cooling in triggering the kelp radiation and pre-dates the origin of benthic-foraging carnivores. Historical biogeographical reconstructions point to a northeast Pacific origin of complex kelps, with subsequent colonization of new habitats likely playing an important role in driving their ecological diversification. We infer that complex morphologies associated with modern kelp forests (e.g. branching, pneumatocysts) evolved several times over the past 15–20 MY, highlighting the importance of morphological convergence in establishing modern upright kelp forests. Our phylogenomic findings provide new insights into the geographical and ecological proliferation of kelps and provide a timeline along which feedbacks between kelps and their food-webs could have shaped the structure of temperate ecosystems.

The elongation of very-long-chain fatty acids is a conserved process used for the production of many metabolites, including plant cuticular waxes. The elongation of precursors of the most abundant cuticular wax components of some plants, however, is unique in requiring ECERIFERUM2-LIKE (CER2-LIKE) proteins. CER2-LIKEs are a clade within the BAHD superfamily of acyltransferases. They are known to be required for cuticular wax production in both Arabidopsis and maize based on mutant studies. Heterologous expression of Arabidopsis and rice CER2-LIKEs in Saccharomyces cerevisiae has demonstrated that they modify the chain-length specificity of elongation when paired with particular condensing enzymes. Despite sequence homology, CER2-LIKEs are distinct from the BAHD superfamily in that they do not appear to use acyl transfer activity to fulfill their biological function. Here, we review the discovery and characterization of CER2-LIKEs, propose several models to explain their function, and explore the importance of CER2-LIKE proteins for the evolution of plant cuticles.

Despite the extraordinary significance leaves have for life on Earth, their origin and development remain vigorously debated. More than a century of paleobotanical, morphological, and phylogenetic research has still not resolved fundamental questions about leaves. Developmental genetic data are sparse in ferns, and comparative studies of lycophytes and seed plants have reached opposing conclusions on the conservation of a leaf developmental program. We performed phylogenetic and expression analyses of a leaf developmental regulator (Class III HD-Zip genes; C3HDZs) spanning lycophytes and ferns. We show that a duplication and neofunctionalization of C3HDZs probably occurred in the ancestor of euphyllophytes, and that there is a common leaf developmental mechanism conserved between ferns and seed plants. We show C3HDZ expression in lycophyte and fern sporangia and show that C3HDZs have conserved expression patterns during initiation of lateral primordia (leaves or sporangia). This expression is maintained throughout sporangium development in lycophytes and ferns and indicates an ancestral role of C3HDZs in sporangium development. We hypothesize that there is a deep homology of all leaves and that a sporangium-specific developmental program was coopted independently for the development of lycophyte and euphyllophyte leaves. This provides molecular genetic support for a paradigm shift in theories of lycophyte leaf evolution.

Loss of introns in plant mitochondrial genes is commonly explained by retroprocessing. Under this model, an mRNA is reverse transcribed and integrated back into the genome, simultaneously affecting the contents of introns and edited sites. To evaluate the extent to which retroprocessing explains intron loss, we analyzed patterns of intron content and predicted RNA editing for whole mitochondrial genomes of 30 species in the monocot order Alismatales. In this group, we found an unusually high degree of variation in the intron content, even expanding the hitherto known variation among angiosperms. Some species have lost some two-third of the cis-spliced introns. We found a strong correlation between intron content and editing frequency, and detected 27 events in which intron loss is consistent with the presence of nucleotides in an edited state, supporting retroprocessing. However, we also detected seven cases of intron loss not readily being explained by retroprocession. Our analyses are also not consistent with the entire length of a fully processed cDNA copy being integrated into the genome, but instead indicate that retroprocessing usually occurs for only part of the gene. In some cases, several rounds of retroprocessing may explain intron loss in genes completely devoid of introns. A number of taxa retroprocessing seem to be very common and a possibly ongoing process. It affects the entire mitochondrial genome.

We present the first phylogenomic analysis of relationships among all ten families of Liliales, based on 75 plastid genes from 35 species in 29 genera, and 97 additional plastomes stratified across angiosperm lineages. We used a supermatrix approach to extend our analysis to 58 of 64 genera of Liliales, and calibrated the resulting phylogeny against 17 fossil dates to produce a new timeline for monocot evolution. Liliales diverged from other monocots 124 Mya and began splitting into separate families 113 Mya. Our data support an Australian origin for Liliales, with close relationships between three pairs of lineages (Corsiaceae/Campynemataceae, Philesiaceae/Ripogonaceae, tribes Alstroemerieae/Luzuriageae) in South America and Australia or New Zealand reflecting teleconnections of these areas via Antarctica. Long-distance dispersal (LDD) across the Pacific and Tasman Sea led to re-invasion of New Zealand by two lineages (Luzuriaga, Ripogonum); LDD allowed Campynemanthe to colonize New Caledonia after its submergence until 37 Mya. LDD permitted Colchicaceae to invade East Asia and Africa from Australia, and re-invade Africa from Australia. Periodic desert greening permitted Gloriosa and Iphigenia to colonize Southeast Asia overland from Africa, and Androcymbium-Colchicum to invade the Mediterranean from South Africa. Melanthiaceae and Liliaceae crossed the Bering land-bridge several times from the Miocene to the Pleistocene.

Background and Aims: Rapid evolutionary divergence and reticulate evolution may result in phylogenetic relationships that are difficult to resolve using small nucleotide sequence data sets. Next-generation sequencing methods can generate larger data sets that are better suited to solving these puzzles. One major and long-standing controversy in conifers concerns generic relationships within the subfamily Cupressoideae (105 species, approx. 1/6 of all conifers) of Cupressaceae, in particular the relationship between Juniperus, Cupressus and the Hesperocyparis-Callitropsis-Xanthocyparis (HCX) clade. Here we attempt to resolve this question using transcriptome-derived data.Methods: Transcriptome sequences of 20 species from Cupressoideae were collected. Using MarkerMiner, single-copy nuclear (SCN) genes were extracted. These were applied to estimate phylogenies based on concatenated data, species trees and a phylogenetic network. We further examined the effect of alternative backbone topologies on downstream analyses, including biogeographic inference and dating analysis.Results: Based on the 73 SCN genes (>200 000 bp total alignment length) we considered, all tree-building methods lent strong support for the relationship (HCX, (Juniperus, Cupressus)); however, strongly supported conflicts among individual gene trees were also detected. Molecular dating suggests that these three lineages shared a most recent common ancestor approx. 60 million years ago (Mya), and that Juniperus and Cupressus diverged about 56 Mya. Ancestral area reconstructions (AARs) suggest an Asian origin for the entire clade, with subsequent dispersal to North America, Europe and Africa.Conclusions: Our analysis of SCN genes resolves a controversial phylogenetic relationship in the Cupressoideae, a major clade of conifers, and suggests that rapid evolutionary divergence and incomplete lineage sorting probably acted together as the source for conflicting phylogenetic inferences between gene trees and between our robust results and recently published studies. Our updated backbone topology has not substantially altered molecular dating estimates relative to previous studies; however, application of the latest AAR approaches has yielded a clearer picture of the biogeographic history of Cupressoideae.

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).

et al., 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., The work conducted by the US Department of Energy Joint Genome Institute is supported by the Office of Science of the US Department of Energy under Contract No. DE-AC02-05CH11231. Support to RR and SAR by the German Research Foundation (DFG RE 837/10-2), the Excellence Initiative of the German Federal and State Governments (EXC 294), and by the German Federal Ministry of Education and Research (BMBF FRISYS), is highly appreciated. CoGe is supported by the US National Science Foundation under Award Numbers IOS-339156 and IOS-1444490, CyVerse is supported by the U.S. National Science Foundation under Award Numbers DBI-0735191 and DBI-1265383. YK and ACC are grateful for support from the UK Biological Sciences and Biotechnology Research Council (Grant BB/F001797/1). KV acknowledges the Multidisciplinary Research Partnership ‘Bioinformatics: from nucleotides to networks’ Project (no 01MR0410W) of Ghent University. JC is grateful for support from the Spanish Ministerio de Economía y Competitividad (Grant AGL2013-43244-R). RSQ is grateful to Monsanto (St. Louis, MO, USA) for sequencing genomic DNA of P. patens accession Kaskaskia. The 1000 Plants (1 KP) initiative, led by GKSW, is funded by the Alberta Ministry of Innovation and Advanced Education, Alberta Innovates Technology Futures (AITF), Innovates Centres of Research Excellence (iCORE), Musea Ventures, BGI-Shenzhen and China National Genebank (CNGB). TW was supported by EMBO Long-Term Fellowships (ALTF 1166-2011) and by Marie Curie Actions (European Commission EMBOCOFUND2010, GA-2010-267146). The work conducted at PGSB was supported by the German Research Foundation (SFB924) and German Ministry of Education and Research (BMBF, 031A536/de.NBI).

Understanding plant evolution and diversity in a phylogenomic context is an enormous challenge due, in part, to limited availability of genome-scale data across phylodiverse species. The 10KP (10,000 Plants) Genome Sequencing Project will sequence and characterize representative genomes from every major clade of embryophytes, green algae, and protists (excluding fungi) within the next 5 years. By implementing and continuously improving leading-edge sequencing technologies and bioinformatics tools, 10KP will catalogue the genome content of plant and protist diversity and make these data freely available as an enduring foundation for future scientific discoveries and applications. 10KP is structured as an international consortium, open to the global community, including botanical gardens, plant research institutes, universities, and private industry. Our immediate goal is to establish a policy framework for this endeavor, the principles of which are outlined here.

Mycoheterotrophs, i.e. plants that acquire carbon from root-associated soil fungi, often have highly degraded plastomes, reflecting relaxed selective constraints on plastid genes following the loss of photosynthesis. Geosiris Baill. is the only mycoheterotrophic genus in Iridaceae and comprises two species in Madagascar and nearby islands, and a third recently discovered species in north-eastern Australia. Here, we characterise the plastomes of the Australian and one Madagascan species to compare patterns of plastome degradation in relation to autotrophic and other mycoheterotrophic taxa and investigate the evolutionary and biogeographical history of the genus in Iridaceae. Both examined species have lost approximately half their plastid-encoded genes and a small but significant reduction in purifying selection in retained non-photosynthetic genes was observed. Geosiris is confirmed as monophyletic, with initial divergence of the genus occurring c. 53 million years ago, and subsequent diversification occurring c. 30 million years ago. Africa (including Madagascar) is reconstructed as the most likely ancestral area of the genus, implying a major range-expansion event of one lineage to Australia after its divergence in the Oligocene. Our study has highlighted the dynamic evolutionary history of Geosiris, contributed to the characterisation of mycoheterotrophic plastomes, and furthered our understanding of plastome structure and function.

3 PREMISE OF THE STUDY: Understanding fern (monilophyte) phylogeny and its evolutionary timescale is critical for broad investigations of the evolution of land plants, and for providing the point of comparison necessary for studying the evolution of the fern sister group, seed plants. Molecular phylogenetic investigations have revolutionized our understanding of fern phylogeny, however, to date, these studies have relied almost exclusively on plastid data. METHODS: Here we take a curated phylogenomics approach to infer the fibroad fern phylogeny from multiple nuclear loci, by combining broad taxon sampling (73 ferns and 12 outgroup species) with focused character sampling (25 loci comprising 35 877 bp), along with rigorous alignment, orthology inference and model selection. KEY RESULTS: Our phylogeny corroborates some earlier inferences and provides novel insights; in particular, we fi nd strong support for Equisetales as sister to the rest of ferns, Marattiales as sister to leptosporangiate ferns, and Dennstaedtiaceae as sister to the eupolypods. Our divergence-time analyses reveal that divergences among the extant fern orders all occurred prior to ~200 MYA. Finally, our species-tree inferences are congruent with analyses of concat- enated data, but generally with lower support. Those cases where species-tree support values are higher than expected involve relationships that have been supported by smaller plastid datasets, suggesting that deep coalescence may be reducing support from the concatenated nuclear data. CONCLUSIONS: Our study demonstrates the utility of a curated phylogenomics approach to inferring fern phylogeny, and highlights the need to consider underlying data characteristics, along with data quantity, in phylogenetic studies.

The enigmatic monocot family Triuridaceae provides a potentially useful model system for studying the effects of an ancient loss of photosynthesis on the plant plastid genome, as all of its members are mycoheterotrophic and achlorophyllous. However, few studies have placed the family in a comparative context, and its phylogenetic placement is only partly resolved. It was also unclear whether any taxa in this family have retained a plastid genome. Here, we used genome survey sequencing to retrieve plastid genome data for Sciaphila densiflora (Triuridaceae) and ten autotrophic relatives in the orders Dioscoreales and Pandanales. We recovered a highly reduced plastome for Sciaphila that is nearly colinear with Carludovica palmata, a photosynthetic relative that belongs to its sister group in Pandanales, Cyclanthaceae-Pandanaceae. This phylogenetic placement is well supported and robust to a broad range of analytical assumptions in maximum-likelihood inference, and is congruent with recent findings based on nuclear and mitochondrial evidence. The 28 genes retained in the S. densiflora plastid genome are involved in translation and other nonphotosynthetic functions, and we demonstrate that nearly all of the 18 protein-coding genes are under strong purifying selection. Our study confirms the utility of whole plastid genome data in phylogenetic studies of highly modified heterotrophic plants, even when they have substantially elevated rates of substitution.

A phylogenetic analysis of the early branching lineages of the monocotyledons is performed using data from two plastid genes (rbcL and matK), five mitochondrial genes (atp1, ccmB, cob, mttB and nad5) and morphology. The complete matrix includes 93 terminals representing Acorus, the 14 families currently recognized within Alismatales, and numerous lineages of monocotyledons and other angiosperms. Total evidence analysis results in an almost completely resolved strict consensus tree, but all data partitions, genomic as well as morphological, are incongruent. The effects of RNA editing and potentially processed paralogous sequences are explored and discussed. Despite a decrease in incongruence length differences after exclusion of edited sites, the major data partitions remain significantly incongruent. The 14 families of Alismatales are all found to be monophyletic, but Acorus is found to be included in Alismatales rather than being the sister group to all other monocotyledons. The placement is strongly supported by the mitochondrial data, atp1 in particular, but it cannot be explained as an artifact caused by patterns of editing or by sampling of processed paralogues.

Aim Many plant families have a disjunct distribution across the southern Pacific Ocean, including the mycoheterotrophic family Corsiaceae, which provides a prime example of this biogeographical pattern. A better grasp of the family's evolutionary relationships is needed to understand its historical biogeography. We therefore aimed to (1) test the uncertain monophyly of Corsiaceae, (2) define its phylogenetic position, and (3) estimate divergence times for the family, allowing us to assess whether the distribution of the family is the result of vicariance. Location Southern South America and Australasia. Methods We analysed various combinations of mitochondrial and nuclear data to address the monophyly, phylogenetic position and age of Corsiaceae. To test its monophyly, we used a three-locus data set including most monocot orders, and to infer its exact phylogenetic position, we used a five-locus extended data set. We corroborated these findings using an independent plastome dataset. We then used a two-locus dataset with taxa from all monocot orders, and a three-locus dataset containing only taxa of Liliales, to estimate divergence times using a fossil-calibrated uncorrelated lognormal relaxed-clock approach. Results Corsiaceae is a monophyletic family and the sister group of Campynemataceae. This clade is the sister group of all other Liliales. The crown age of Corsiaceae is estimated to be 53 Ma (95% confidence interval 30–76 Ma). Main conclusions Corsiaceae is an ancient family of mycoheterotrophic plants, whose crown age overlaps with the plate-tectonic split of Gondwana, consistent with a vicariance-based explanation for its current distribution.

Recent re-examinations and new fossil findings have added significantly to the data available for evaluating the evolutionary history of the monocotyledons. Integrating data from the monocot fossil record with molecular dating techniques has the potential to help us to understand better the timing of important evolutionary events and patterns of diversification and extinction in this major and ancient clade of flowering plants. In general, the oldest well-placed fossils are used to constrain the age of nodes in molecular dating analyses. However, substantial error can be introduced if calibration fossils are not carefully evaluated and selected. Here we propose a set of 34 fossils representing 19 families and eight orders for calibrating the ages of major monocot clades. We selected these fossils because they can be placed in particular clades with confidence and they come from well-dated stratigraphic sequences. As more fossils are discovered or re-examined, these criteria can also be applied to expand the list of the fossils that are most suitable for dating the early branches of monocot phylogeny. © 2015 The Linnean Society of London, Botanical Journal of the Linnean Society, 2015, 178, 346–374.

Background Strigolactones (SLs) are a class of plant hormones that control many aspects of plant growth. The SL signalling mechanism is homologous to that of karrikins (KARs), smoke-derived compounds that stimulate seed germination. In angiosperms, the SL receptor is an α/β-hydrolase known as DWARF14 (D14); its close homologue, KARRIKIN INSENSITIVE2 (KAI2), functions as a KAR receptor and likely recognizes an uncharacterized, endogenous signal (‘KL’). Previous phylogenetic analyses have suggested that the KAI2 lineage is ancestral in land plants, and that canonical D14-type SL receptors only arose in seed plants; this is paradoxical, however, as non-vascular plants synthesize and respond to SLs. Results We have used a combination of phylogenetic and structural approaches to re-assess the evolution of the D14/KAI2 family in land plants. We analysed 339 members of the D14/KAI2 family from land plants and charophyte algae. Our phylogenetic analyses show that the divergence between the eu-KAI2 lineage and the DDK (D14/DLK2/KAI2) lineage that includes D14 occurred very early in land plant evolution. We show that eu-KAI2 proteins are highly conserved, and have unique features not found in DDK proteins. Conversely, we show that DDK proteins show considerable sequence and structural variation to each other, and lack clearly definable characteristics. We use homology modelling to show that the earliest members of the DDK lineage structurally resemble KAI2 and that SL receptors in non-seed plants likely do not have D14-like structure. We also show that certain groups of DDK proteins lack the otherwise conserved MORE AXILLARY GROWTH2 (MAX2) interface, and may thus function independently of MAX2, which we show is highly conserved throughout land plant evolution. Conclusions Our results suggest that D14-like structure is not required for SL perception, and that SL perception has relatively relaxed structural requirements compared to KAI2-mediated signalling. We suggest that SL perception gradually evolved by neo-functionalization within the DDK lineage, and that the transition from KAI2-like to D14-like protein may have been driven by interactions with protein partners, rather than being required for SL perception per se. Electronic supplementary material The online version of this article (doi:10.1186/s12915-017-0397-z) contains supplementary material, which is available to authorized users.

The order Alismatales is a hotspot for evolution of plant mitochondrial genomes characterized by remarkable differences in genome size, substitution rates, RNA editing, retrotranscription, gene loss and intron loss. Here we have sequenced the complete mitogenomes of Zostera marina and Stratiotes aloides, which together with previously sequenced mitogenomes from Butomus and Spirodela, provide new evolutionary evidence of genome size reduction, gene loss and transfer to the nucleus. The Zostera mitogenome includes a large portion of DNA transferred from the plastome, yet it is the smallest known mitogenome from a non-parasitic plant. Using a broad sample of the Alismatales, the evolutionary history of ribosomal protein gene loss is analyzed. In Zostera almost all ribosomal protein genes are lost from the mitogenome, but only some can be found in the nucleus.

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.

Strigolactones (SLs) are a class of plant hormones that control many aspects of plant growth. The SL signalling mechanism is homologous to that of karrikins (KARs), smoke-derived compounds that stimulate seed germination. In angiosperms, the SL receptor is an α/β hydrolase known as DWARF14 (D14); its close homologue, KARRIKIN INSENSITIVE2 (KAI2), functions as a KAR receptor, and likely recognizes an uncharacterized, endogenous signal. Previous phylogenetic analyses have suggested that the KAI2 lineage is ancestral in land plants, and that canonical D14-type SL receptors only arose in seed plants; this is paradoxical, however, as non-vascular plants synthesize and respond to SLs. Here, we have used a combination of phylogenetic and structural approaches to re-assess the evolution of the D14/KAI2 family in land plants. We analyzed 339 members of the D14/KAI2 family from land plants and charophyte algae. Our phylogenetic analyses show that the divergence between the eu-KAI2 lineage and the DDK (D14/DLK2/KAI2) lineage that includes D14 occurred very early in land plant evolution. We identify characteristic structural features of D14 and KAI2 proteins, and use homology modelling to show that the earliest members of the DDK lineage structurally resemble KAI2, and not D14 proteins. Furthermore, we show that probable SL receptors in non-seed plants do not have D14-like structure. Our results suggest that SL perception has relatively relaxed structural requirements, and that the evolution from KAI2-like to D14-like protein structure in the DDK lineage may have been driven by interactions with protein partners, rather than being required for SL perception itself.

It is commonly believed that gene duplications provide the raw material for morphological evolution. Both the number of genes and size of gene families have increased during the diversification of land plants. Several small proteins that regulate transcription factors have recently been identified in plants, including the LITTLE ZIPPER (ZPR) proteins. ZPRs are post-translational negative regulators, via heterodimerization, of class III Homeodomain Leucine Zipper (C3HDZ) proteins that play a key role in directing plant form and growth. We show that ZPR genes originated as a duplication of a C3HDZ transcription factor paralog in the common ancestor of euphyllophytes (ferns and seed plants). The ZPRs evolved by degenerative mutations resulting in loss all of the C3HDZ functional domains, except the leucine zipper that modulates dimerization. ZPRs represent a novel regulatory module of the C3HDZ network unique to the euphyllophyte lineage, and their origin correlates to a period of rapid morphological changes and increased complexity in land plants. The origin of the ZPRs illustrates the significance of gene duplications in creating developmental complexity during land plant evolution that likely led to morphological evolution.

Willows (Salix: Salicaceae) form a major ecological component of Holarctic floras and consequently are an obvious target for a DNA-based identification system. We surveyed two to seven plastid genome regions (~3.8 kb; ~3% of the genome) from 71 Salix species across all five subgenera, to assess their performance as DNA barcode markers. Although Salix has a relatively high level of interspecific hybridization, this may not sufficiently explain the near complete failure of barcoding that we observed: only one species had a unique barcode. We recovered 39 unique haplotypes, from more than 500 specimens, that could be partitioned into six major haplotype groups. A unique variant of group I (haplotype 1*) was shared by 53 species in three of five Salix subgenera. This unusual pattern of haplotype sharing across infrageneric taxa is suggestive of either a massive nonrandom coalescence failure (incomplete lineage sorting), or of repeated plastid capture events, possibly including a historical selective sweep of haplotype 1* across taxonomic sections. The former is unlikely as molecular dating indicates that haplotype 1* originated recently and is nested in the oldest major haplotype group in the genus. Further, we detected significant non-neutrality in the frequency spectrum of mutations in group I, but not outside group I, and demonstrated a striking absence of geographical (isolation by distance) effects in the haplotype distributions of this group. The most likely explanation for the patterns we observed involves recent repeated plastid capture events, aided by widespread hybridization and long-range seed dispersal, but primarily propelled by one or more trans-species selective sweeps.