Phylotranscriptomic Analyses of Mycoheterotrophic Monocots Show a Continuum of Convergent Evolutionary Changes in Expressed Nuclear Genes From Three Independent Nonphotosynthetic Lineages.

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.

Phylogenomic resolution of order- and family-level monocot relationships using 602 single-copy nuclear genes and 1375 BUSCO genes.

We assess relationships among 192 species in all 12 monocot orders and 72 of 77 families, using 602 conserved single-copy (CSC) genes and 1375 benchmarking single-copy ortholog (BUSCO) genes extracted from genomic and transcriptomic datasets. Phylogenomic inferences based on these data, using both coalescent-based and supermatrix analyses, are largely congruent with the most comprehensive plastome-based analysis, and nuclear-gene phylogenomic analyses with less comprehensive taxon sampling. The strongest discordance between the plastome and nuclear gene analyses is the monophyly of a clade comprising Asparagales and Liliales in our nuclear gene analyses, versus the placement of Asparagales and Liliales as successive sister clades to the commelinids in the plastome tree. Within orders, around six of 72 families shifted positions relative to the recent plastome analysis, but four of these involve poorly supported inferred relationships in the plastome-based tree. In Poales, the nuclear data place a clade comprising Ecdeiocoleaceae+Joinvilleaceae as sister to the grasses (Poaceae); Typhaceae, (rather than Bromeliaceae) are resolved as sister to all other Poales. In Commelinales, nuclear data place Philydraceae sister to all other families rather than to a clade comprising Haemodoraceae+Pontederiaceae as seen in the plastome tree. In Liliales, nuclear data place Liliaceae sister to Smilacaceae, and Melanthiaceae are placed sister to all other Liliales except Campynemataceae. Finally, in Alismatales, nuclear data strongly place Tofieldiaceae, rather than Araceae, as sister to all the other families, providing an alternative resolution of what has been the most problematic node to resolve using plastid data, outside of those involving achlorophyllous mycoheterotrophs. As seen in numerous prior studies, the placement of orders Acorales and Alismatales as successive sister lineages to all other extant monocots. Only 21.2% of BUSCO genes were demonstrably single-copy, yet phylogenomic inferences based on BUSCO and CSC genes did not differ, and overall functional annotations of the two sets were very similar. Our analyses also reveal significant gene tree-species tree discordance despite high support values, as expected given incomplete lineage sorting (ILS) related to rapid diversification. Our study advances understanding of monocot relationships and the robustness of phylogenetic inferences based on large numbers of nuclear single-copy genes that can be obtained from transcriptomes and genomes.

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

Mitochondrial genome evolution in Alismatales: Size reduction and extensive loss of ribosomal protein genes.

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.

Localized Retroprocessing as a Model of Intron Loss in the Plant Mitochondrial Genome.

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.

Drastic reduction of plastome size in the mycoheterotrophic Thismia tentaculata relative to that of its autotrophic relative Tacca chantrieri.

PREMISE OF THE STUDY: Heterotrophic angiosperms tend to have reduced plastome sizes relative to those of their autotrophic relatives because genes that code for proteins involved in photosynthesis are lost. However, some plastid-encoded proteins may have vital nonphotosynthetic functions, and the plastome therefore may be retained after the loss of photosynthesis. METHODS: We sequenced the plastome of the mycoheterotrophic species Thismia tentaculata and a representative of its sister genus, Tacca chantrieri, using next-generation technology, and we compared sequences and structures of genes and genomes of these species. KEY RESULTS: The plastome of Tacca chantrieri is similar to those of other autotrophic taxa of Dioscoreaceae, except in a few local rearrangements and one gene loss. The plastome of Thismia tentaculata is ca. 16 kbp long with a quadripartite structure and is among the smallest known plastomes. Synteny is minimal between the plastomes of Tacca chantrieri and Thismia tentaculata. The latter includes only 12 candidate genes, with all except accD involved in protein synthesis. Of the 12 genes, trnE, trnfM, and accD are frequently among the few that remain in depauperate plastomes. CONCLUSIONS: The plastome of Thismia tentaculata, like those of most other heterotrophic plants, includes a small number of genes previously suggested to be essential to plastome survival.

A Phylogenomic Assessment of Ancient Polyploidy and Genome Evolution across the Poales.

Comparisons of flowering plant genomes reveal multiple rounds of ancient polyploidy characterized by large intragenomic syntenic blocks. Three such whole-genome duplication (WGD) events, designated as rho (ρ), sigma (σ), and tau (τ), have been identified in the genomes of cereal grasses. Precise dating of these WGD events is necessary to investigate how they have influenced diversification rates, evolutionary innovations, and genomic characteristics such as the GC profile of protein-coding sequences. The timing of these events has remained uncertain due to the paucity of monocot genome sequence data outside the grass family (Poaceae). Phylogenomic analysis of protein-coding genes from sequenced genomes and transcriptome assemblies from 35 species, including representatives of all families within the Poales, has resolved the timing of rho and sigma relative to speciation events and placed tau prior to divergence of Asparagales and the commelinids but after divergence with eudicots. Examination of gene family phylogenies indicates that rho occurred just prior to the diversification of Poaceae and sigma occurred before early diversification of Poales lineages but after the Poales-commelinid split. Additional lineage-specific WGD events were identified on the basis of the transcriptome data. Gene families exhibiting high GC content are underrepresented among those with duplicate genes that persisted following these genome duplications. However, genome duplications had little overall influence on lineage-specific changes in the GC content of coding genes. Improved resolution of the timing of WGD events in monocot history provides evidence for the influence of polyploidization on functional evolution and species diversification.

Phylogeny of the Alismatales (Monocotyledons) and the relationship of Acorus (Acorales?).

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.

Plastid phylogenomics and molecular evolution of Alismatales.

Past phylogenetic studies of the monocot order Alismatales left several higher-order relationships unresolved. We addressed these uncertainties using a nearly complete genus-level sampling of whole plastid genomes (gene sets representing 83 protein-coding and ribosomal genes) from members of the core alismatid families, Tofieldiaceae and additional taxa (Araceae and other angiosperms). Parsimony and likelihood analyses inferred generally highly congruent phylogenetic relationships within the order, and several alternative likelihood partitioning schemes had little impact on patterns of clade support. All families with multiple genera were resolved as monophyletic, and we inferred strong bootstrap support for most inter- and intrafamilial relationships. The precise placement of Tofieldiaceae in the order was not well supported. Although most analyses inferred Tofieldiaceae to be the sister-group of the rest of the order, one likelihood analysis indicated a contrasting Araceae-sister arrangement. Acorus (Acorales) was not supported as a member of the order. We also investigated the molecular evolution of plastid NADH dehydrogenase, a large enzymatic complex that may play a role in photooxidative stress responses. Ancestral-state reconstructions support four convergent losses of a functional NADH dehydrogenase complex in Alismatales, including a single loss in Tofieldiaceae.

Plastid genomes reveal support for deep phylogenetic relationships and extensive rate variation among palms and other commelinid monocots.

Despite progress based on multilocus, phylogenetic studies of the palms (order Arecales, family Arecaceae), uncertainty remains in resolution/support among major clades and for the placement of the palms among the commelinid monocots. Palms and related commelinids represent a classic case of substitution rate heterogeneity that has not been investigated in the genomic era. To address questions of relationships, support and rate variation among palms and commelinid relatives, 39 plastomes representing the palms and related family Dasypogonaceae were generated via genome skimming and integrated within a monocot-wide matrix for phylogenetic and molecular evolutionary analyses. Support was strong for 'deep' relationships among the commelinid orders, among the five palm subfamilies, and among tribes of the subfamily Coryphoideae. Additionally, there was extreme heterogeneity in the plastid substitution rates across the commelinid orders indicated by model based analyses, with c. 22 rate shifts, and significant departure from a global clock. To date, this study represents the most comprehensively sampled matrix of plastomes assembled for monocot angiosperms, providing genome-scale support for phylogenetic relationships of monocot angiosperms, and lays the phylogenetic groundwork for comparative analyses of the drivers and correlates of such drastic differences in substitution rates across a diverse and significant clade.

Resolving relationships within the palm subfamily Arecoideae (Arecaceae) using plastid sequences derived from next-generation sequencing.

PREMISE OF THE STUDY: Several studies have incorporated molecular and morphological data to study the phylogeny of the palms (Arecaceae), but some relationships within the family remain ambiguous-particularly those within Arecoideae, the most diverse subfamily including coconut and oil palm. Here, two next-generation, targeted plastid-enrichment methods were compared and used to elucidate Arecoideae phylogeny. METHODS: Next-generation sequencing techniques were used to generate a plastid genome data set. Long range PCR and hybrid gene capture were used to enrich for chloroplast targets. Ten taxa were enriched using both methods for comparison. Chloroplast sequence data were generated for 31 representatives of the 14 Arecoideae tribes and five outgroup taxa. The phylogeny was reconstructed using maximum likelihood, maximum parsimony, and Bayesian analyses. KEY RESULTS: Long range PCR and hybrid gene capture both enriched the plastid genome and provided similar sequencing coverage. Subfamily Arecoideae was resolved as monophyletic with tribe Chamaedoreeae as the earliest-diverging lineage, implying that the development of flowers in triads defines a synapomorphy for the Arecoideae clade excluding Chamaedoreeae. Three major clades within this group were recovered: Roystoneeae/Reinhardtieae/Cocoseae (RRC), Areceae/Euterpeae/Geonomateae/Leopoldinieae/Manicarieae/Pelagodoxeae (core arecoids), and Podococceae/Oranieae/Sclerospermeae (POS). An Areceae + Euterpeae clade was resolved within the core arecoids. The POS clade was sister to a RRC + core arecoids clade, implying a shared ancestral area in South America for these three clades. CONCLUSIONS: The plastome phylogeny recovered here provides robust resolution of previously ambiguous studies and new insights into palm evolution.

Plastid phylogenomics of the cool-season grass subfamily: clarification of relationships among early-diverging tribes.

Whole plastid genomes are being sequenced rapidly from across the green plant tree of life, and phylogenetic analyses of these are increasing resolution and support for relationships that have varied among or been unresolved in earlier single- and multi-gene studies. Pooideae, the cool-season grass lineage, is the largest of the 12 grass subfamilies and includes important temperate cereals, turf grasses and forage species. Although numerous studies of the phylogeny of the subfamily have been undertaken, relationships among some 'early-diverging' tribes conflict among studies, and some relationships among subtribes of Poeae have not yet been resolved. To address these issues, we newly sequenced 25 whole plastomes, which showed rearrangements typical of Poaceae. These plastomes represent 9 tribes and 11 subtribes of Pooideae, and were analysed with 20 existing plastomes for the subfamily. Maximum likelihood (ML), maximum parsimony (MP) and Bayesian inference (BI) robustly resolve most deep relationships in the subfamily. Complete plastome data provide increased nodal support compared with protein-coding data alone at nodes that are not maximally supported. Following the divergence of Brachyelytrum, Phaenospermateae, Brylkinieae-Meliceae and Ampelodesmeae-Stipeae are the successive sister groups of the rest of the subfamily. Ampelodesmeae are nested within Stipeae in the plastome trees, consistent with its hybrid origin between a phaenospermatoid and a stipoid grass (the maternal parent). The core Pooideae are strongly supported and include Brachypodieae, a Bromeae-Triticeae clade and Poeae. Within Poeae, a novel sister group relationship between Phalaridinae and Torreyochloinae is found, and the relative branching order of this clade and Aveninae, with respect to an Agrostidinae-Brizinae clade, are discordant between MP and ML/BI trees. Maximum likelihood and Bayesian analyses strongly support Airinae and Holcinae as the successive sister groups of a Dactylidinae-Loliinae clade.

Resolving ancient radiations: can complete plastid gene sets elucidate deep relationships among the tropical gingers (Zingiberales)?

BACKGROUND AND AIMS: Zingiberales comprise a clade of eight tropical monocot families including approx. 2500 species and are hypothesized to have undergone an ancient, rapid radiation during the Cretaceous. Zingiberales display substantial variation in floral morphology, and several members are ecologically and economically important. Deep phylogenetic relationships among primary lineages of Zingiberales have proved difficult to resolve in previous studies, representing a key region of uncertainty in the monocot tree of life. METHODS: Next-generation sequencing was used to construct complete plastid gene sets for nine taxa of Zingiberales, which were added to five previously sequenced sets in an attempt to resolve deep relationships among families in the order. Variation in taxon sampling, process partition inclusion and partition model parameters were examined to assess their effects on topology and support. KEY RESULTS: Codon-based likelihood analysis identified a strongly supported clade of ((Cannaceae, Marantaceae), (Costaceae, Zingiberaceae)), sister to (Musaceae, (Lowiaceae, Strelitziaceae)), collectively sister to Heliconiaceae. However, the deepest divergences in this phylogenetic analysis comprised short branches with weak support. Additionally, manipulation of matrices resulted in differing deep topologies in an unpredictable fashion. Alternative topology testing allowed statistical rejection of some of the topologies. Saturation fails to explain observed topological uncertainty and low support at the base of Zingiberales. Evidence for conflict among the plastid data was based on a support metric that accounts for conflicting resampled topologies. CONCLUSIONS: Many relationships were resolved with robust support, but the paucity of character information supporting the deepest nodes and the existence of conflict suggest that plastid coding regions are insufficient to resolve and support the earliest divergences among families of Zingiberales. Whole plastomes will continue to be highly useful in plant phylogenetics, but the current study adds to a growing body of literature suggesting that they may not provide enough character information for resolving ancient, rapid radiations.

Plastid genomes and deep relationships among the commelinid monocot angiosperms.

The commelinid monocots comprise the orders Arecales (A), Commelinales (C), Poales sensu Angiosperm Phylogeny Group III (APGIII) (P), Zingiberales (Z), plus the unplaced family Dasypogonaceae (D), collectively containing numerous economically and ecologically important species and encompassing enormous morphological diversity. Commelinids are supported as monophyletic based on anatomy and molecular data; however, relationships among major commelinid groupings conflict among previous studies, representing a long-standing problem in monocot systematics, with major implications for interpretations of character evolution. In more recent analyses, with whole-plastome sampling largely focused on Poales, areas of conflict remain, suggesting the need for closer investigation of relationships and support. Here, we increased sampling of plastomes among non-Poalean commelinid orders to investigate deep nodal support. Analysis of 83 plastid genes recovered relationships as ((A, D) (ZC, P)) with robust support, regardless of reconstruction method (parsimony/likelihood). However, conflict among genes was evident when grouped by genomic region. Cumulative analyses of genes ranked by decreasing numbers of informative characters indicated continued fluctuation in support, even as small genes were added to a nearly complete matrix, contrary to the expected pattern of stabilization in support. Topology tests among major commelinid groups suggested that the data were not powerful enough to reject all alternatives. This study provides clues to the limits of the plastid genome for resolving deep relationships among the commelinid monocots.

Phylogeny of the Liliales (Monocotyledons) with special emphasis on data partition congruence and RNA editing.

A phylogenetic analysis of the monocot order Liliales was performed using sequence data from three mitochondrial (atp1, cob, nad5) and two plastid genes (rbcL, ndhF). The complete data matrix includes 46 terminals representing all 10 families currently included in Liliales. The two major partitions, mitochondrial and plastid data, were congruent, and parsimony analysis resulted in 50 equally parsimonious trees and a well resolved consensus tree confirming monophyly of all families. Mitochondrial genes are known to include RNA edited sites, and in some cases unprocessed genes are replaced by retro-processed gene copies, that is processed paralogs. To test the effects on phylogeny reconstruction of predicted edited sites and potentially unintentionally sampled processed paralogs, a number of analyses were performed using subsets of the complete data matrix. In general, predicted edited sites were more homoplasious than the other characters and increased incongruence among most data partitions. The predicted edited sites have a non-random phylogenetic signal in conflict with the signal of the non-edited sites. The potentially misleading signal was caused partially by the apparent presence of processed paralogs in Galanthus (Amaryllidaceae), part of the outgroup, but also by a deviating evolutionary pattern of predicted edited sites in Liliaceae compared with the remainder of the Liliales. Despite the problems that processed paralogs may cause, we argue that they should not a priori be excluded from phylogenetic analysis.

The plastid genome of the mycoheterotrophic Corallorhiza striata (Orchidaceae) is in the relatively early stages of degradation.

PREMISE OF THE STUDY: Plastid genomes of nonphotosynthetic, mycoheterotrophic plants represent apt systems in which to study effects of relaxed evolutionary constraints. The few mycoheterotrophic angiosperm plastomes sequenced to date display drastic patterns of degradation/reduction relative to those of photosynthetic relatives. The goal of this study was to focus on a mycoheterotrophic orchid hypothesized to be in the "early" stages of plastome degradation, to provide perspective on this process. METHODS: Short-read sequencing was used to generate a complete plastome sequence for Corallorhiza striata var. vreelandii, a mycoheterotrophic orchid, to investigate the extent of plastome degradation. Patterns of nonsynonymous/synonymous mutations were also assessed, and comparisons were made between Corallorhiza and other heterotrophic plant lineages. KEY RESULTS: Corallorhiza yielded a plastome of 137505 bp, with several photosynthesis-related genes either lost or pseudogenized. Members of all major photosynthesis complexes, except ATP-synthase genes, were affected. "Housekeeping" genes were intact, despite the loss of a single tRNA. Intact photosynthesis genes (excluding atp genes) together displayed elevated nonsynonymous changes, while housekeeping genes did not. CONCLUSIONS: The Corallorhiza plastome is not drastically reduced in overall size (∼6% reduction relative to that of photosynthetic Oncidium), but displays a pattern congruent with a loss of photosynthetic function. Comparing Corallorhiza with other heterotrophs allows some emergent evolutionary patterns to be inferred, but these remain as hypotheses to be tested, especially at lower taxonomic levels, and in lineages illustrating transitions from autotrophy to heterotrophy. The independent, unique processes of plastome modification among mycoheterotrophic lineages illustrate the urgency of their conservation.

ABC1K atypical kinases in plants: filling the organellar kinase void.

Surprisingly few protein kinases have been demonstrated in chloroplasts or mitochondria. Here, we discuss the activity of bc(1) complex kinase (ABC1K) protein family, which we suggest locate in mitochondria and plastids, thus filling the kinase void. The ABC1Ks are atypical protein kinases and their ancestral function is the regulation of quinone synthesis. ABC1Ks have proliferated from one or two members in non-photosynthetic organisms to more than 16 members in algae and higher plants. In this review, we reconstruct the evolutionary history of the ABC1K family, provide a functional domain analysis for angiosperms and a nomenclature for ABC1Ks in Arabidopsis (Arabidopsis thaliana), rice (Oryza sativa) and maize (Zea mays). Finally, we hypothesize that targets of ABC1Ks include enzymes of prenyl-lipid metabolism as well as components of the organellar gene expression machineries.

Comparative floral development in Lithospermum (Boraginaceae) and implications for the evolution and development of heterostyly.

PREMISE OF THE STUDY: The evolution and development of floral developmental patterns were investigated in three heterostylous and three homostylous species of Lithospermum to determine whether species that independently acquired the same floral form follow the same pattern of development or different patterns. METHODS: Using light and scanning electron microscopy, we observed developmental patterns in flowers at different stages of maturity. These patterns were compared within individual species, between heterostylous morphs, and among heterostylous and homostylous species. KEY RESULTS: Although heterostyly has been determined by phylogenetic analysis to have originated independently in each of the heterostylous species, flowers of the long-style morph of each species follow similar patterns of gross development, as do those of the short-style morph. In addition, the flowers of each morph develop in a manner similar to those of their respective homostylous, herkogamous relatives. However, the developmental patterns of the stylar epidermal cells differ among these species and between heterostylous and homostylous species. CONCLUSIONS: Floral developmental patterns in homostylous species provide evidence that modification of specific traits, such as patterns of stylar growth, can lead to the evolution of heterostyly. The developmental changes that affect the positions of the stigmas and anthers in each morph likely involve either temporal or spatial modifications of gene function. The floral developmental patterns described here and the occurrence of multiple types of herkogamy within some species of Lithospermum provide evidence that heterostylous species in the genus have originated via distinct evolutionary developmental pathways.

Phylogeny of the Asparagales based on three plastid and two mitochondrial genes.

PREMISE OF THE STUDY: The Asparagales, with ca. 40% of all monocotyledons, include a host of commercially important ornamentals in families such as Orchidaceae, Alliaceae, and Iridaceae, and several important crop species in genera such as Allium, Aloe, Asparagus, Crocus, and Vanilla. Though the order is well defined, the number of recognized families, their circumscription, and relationships are somewhat controversial. METHODS: Phylogenetic analyses of Asparagales were based on parsimony and maximum likelihood using nucleotide sequence variation in three plastid genes (matK, ndhF, and rbcL) and two mitochondrial genes (atp1 and cob). Branch support was assessed using both jackknife analysis implementing strict-consensus (SC) and bootstrap analysis implementing frequency-within-replicates (FWR). The contribution of edited sites in the mitochondrial genes to topology and branch support was investigated. KEY RESULTS: The topologies recovered largely agree with previous results, though some clades remain poorly resolved (e.g., Ruscaceae). When the edited sites were included in the analysis, the plastid and mitochondrial genes were highly incongruent. However, when the edited sites were removed, the two partitions became congruent. CONCLUSIONS: Some deeper nodes in the Asparagales tree remain poorly resolved or unresolved as do the relationships of certain monogeneric families (e.g., Aphyllanthaceae, Ixioliriaceae, Doryanthaceae), whereas support for many families increases. However, the increased support is dominated by plastid data, and the potential influence of mitochondrial and biparentially inherited single or low-copy nuclear genes should be investigated.

Subunit stoichiometry, evolution, and functional implications of an asymmetric plant plastid ClpP/R protease complex in Arabidopsis.

The caseinolytic protease (Clp) protease system has been expanded in plant plastids compared with its prokaryotic progenitors. The plastid Clp core protease consists of five different proteolytic ClpP proteins and four different noncatalytic ClpR proteins, with each present in one or more copies and organized in two heptameric rings. We determined the exact subunit composition and stoichiometry for the intact core and each ring. The chloroplast ClpP/R protease was affinity purified from clpr4 and clpp3 Arabidopsis thaliana null mutants complemented with C-terminal StrepII-tagged versions of CLPR4 and CLPP3, respectively. The subunit stoichiometry was determined by mass spectrometry-based absolute quantification using stable isotope-labeled proteotypic peptides generated from a synthetic gene. One heptameric ring contained ClpP3,4,5,6 in a 1:2:3:1 ratio. The other ring contained ClpP1 and ClpR1,2,3,4 in a 3:1:1:1:1 ratio, resulting in only three catalytic sites. These ClpP1/R1-4 proteins are most closely related to the two subunits of the cyanobacterial P3/R complex and the identical P:R ratio suggests conserved adaptation. Furthermore, the plant-specific C-terminal extensions of the ClpP/R subunits were not proteolytically removed upon assembly, suggesting a regulatory role in Clp chaperone interaction. These results will now allow testing ClpP/R structure-function relationships using rationale design. The quantification workflow we have designed is applicable to other protein complexes.