Highly conserved low-copy nuclear genes as effective markers for phylogenetic analyses in angiosperms

Comparative transcriptome resources of two Dysosma species (Berberidaceae) and molecular evolution of the CYP719A gene in Podophylloideae

Dysosma species (Berberidaceae, Podophylloideae) are of great medicinal pharmacogenetic importance and used as model systems to study the drivers and mechanisms of species diversification of temperate plants in East Asia. Recently, we have sequenced the transcriptome of the low-elevation D. versipellis. In this study, we sequenced the transcriptome of the high-elevation D. aurantiocaulis and used comparative genomic approaches to investigate the transcriptome evolution of the two species. We retrieved 53,929 unigenes from D. aurantiocaulis by de novo transcriptome assemblies using the Illumina HiSeq 2000 platform. Comparing the transcriptomes of both species, we identified 4593 orthologs. Estimation of Ka/Ks ratios for 3126 orthologs revealed that none had a Ka/Ks significantly greater than 1, whereas 1273 (Ka/Ks < 0.5, P < 0.05) were inferred to be under purifying selection. A total of 51 primer pairs were successfully designed from 461 EST-SSRs contained in 4593 orthologs. Marker validation assay revealed that 26 (51%) and 41 (80.4%) produced clear fragments with the expected sizes in all Podophylloideae species. Specifically, 19 different sequences of CYP719A were identified from PCR-amplified genomic DNA of all 12 species of Podophylloideae using primers designed from the assembled transcripts. The data further indicated that CYP719A was likely subject to strong selective constraints maintaining only one copy per genome. In Dysosma, there was relaxed purifying selection or more positive selection for high-elevation species. Overall, this study has generated a wealth of molecular resources potentially useful for pharmacogenetic and evolutionary studies in Dysosma and allied taxa.

Assessing the impact of phylogenetic incongruence on taxonomy, floral evolution, biogeographical history, and phylogenetic diversity

PREMISE OF THE STUDY

Phylogenetic incongruence between "gene trees" and "species trees" has been widely acknowledged in phylogenetic research. Conflicts may emerge from several processes including paralogy, hybridization, and incomplete lineage sorting. Although phylogenetic incongruence appears common, its impact on many phylogeny-based analyses remains poorly understood.

METHODS

We examined the occurrence of phylogenetic conflict between nuclear (ribosome ITS) and plastid (rbcL, trnL-F, rpl20-rps12, and rps16 intron) loci in the ancient angiosperm family Chloranthaceae. Then we investigated how phylogenetic conflict bears on taxonomic classification within the family as well as on inferences on biogeographical history, floral evolution, and measures of phylogenetic diversity (PD).

KEY RESULTS

We found evidence for significant phylogenetic incongruence between plastid and nuclear data in the genus Hedyosmum. Within Hedyosmum, our results did not support previous subgeneric classification of the genus. Division of sections within subgenus Tafalla was supported by the ITS data but not by the plastid data set. As a consequence, we showed that inferring the evolution of key floral characters and geographical history within Hedyosmum depends on the phylogenetic data used. Both data sets yielded similar PD measures across genera, but we found contrasting PD measures in Hedyosmum, even after correcting for rate heterogeneity.

CONCLUSIONS

Our study demonstrated that phylogenetic conflict not only affects the inference of organismal relationships but also impacts our understanding of biogeographical history, morphological evolution, and phylogenetic diversity.

Genome-wide analysis reveals diverged patterns of codon bias, gene expression, and rates of sequence evolution in picea gene families

The recent sequencing of several gymnosperm genomes has greatly facilitated studying the evolution of their genes and gene families. In this study, we examine the evidence for expression-mediated selection in the first two fully sequenced representatives of the gymnosperm plant clade (Picea abies and Picea glauca). We use genome-wide estimates of gene expression (>50,000 expressed genes) to study the relationship between gene expression, codon bias, rates of sequence divergence, protein length, and gene duplication. We found that gene expression is correlated with rates of sequence divergence and codon bias, suggesting that natural selection is acting on Picea protein-coding genes for translational efficiency. Gene expression, rates of sequence divergence, and codon bias are correlated with the size of gene families, with large multicopy gene families having, on average, a lower expression level and breadth, lower codon bias, and higher rates of sequence divergence than single-copy gene families. Tissue-specific patterns of gene expression were more common in large gene families with large gene expression divergence than in single-copy families. Recent family expansions combined with large gene expression variation in paralogs and increased rates of sequence evolution suggest that some Picea gene families are rapidly evolving to cope with biotic and abiotic stress. Our study highlights the importance of gene expression and natural selection in shaping the evolution of protein-coding genes in Picea species, and sets the ground for further studies investigating the evolution of individual gene families in gymnosperms.

Polyploidy-associated genome modifications during land plant evolution

The occurrence of polyploidy in land plant evolution has led to an acceleration of genome modifications relative to other crown eukaryotes and is correlated with key innovations in plant evolution. Extensive genome resources provide for relating genomic changes to the origins of novel morphological and physiological features of plants. Ancestral gene contents for key nodes of the plant family tree are inferred. Pervasive polyploidy in angiosperms appears likely to be the major factor generating novel angiosperm genes and expanding some gene families. However, most gene families lose most duplicated copies in a quasi-neutral process, and a few families are actively selected for single-copy status. One of the great challenges of evolutionary genomics is to link genome modifications to speciation, diversification and the morphological and/or physiological innovations that collectively compose biodiversity. Rapid accumulation of genomic data and its ongoing investigation may greatly improve the resolution at which evolutionary approaches can contribute to the identification of specific genes responsible for particular innovations. The resulting, more 'particulate' understanding of plant evolution, may elevate to a new level fundamental knowledge of botanical diversity, including economically important traits in the crop plants that sustain humanity.

Intercontinental long-distance dispersal of Canellaceae from the New to the Old World revealed by a nuclear single copy gene and chloroplast loci

Canellales, a clade consisting of Winteraceae and Canellaceae, represent the smallest order of magnoliid angiosperms. The clade shows a broad distribution throughout the Southern Hemisphere, across a diverse range of dry to wet tropical forests. In contrast to their sister-group, Winteraceae, the phylogenetic relations and biogeography within Canellaceae remain poorly studied. Here we present the phylogenetic relationships of all currently recognized genera of Canellales with a special focus on the Old World Canellaceae using a combined dataset consisting of the chloroplast trnK-matK-trnK-psbA and the nuclear single copy gene mag1 (Maigo 1). Within Canellaceae we found high statistical support for the monophyly of Warburgia and Cinnamosma. However, we also found relationships that differ from previous studies. Cinnamodendron splitted into two clades, a South American clade and a second clade confined to the Antilles and adjacent areas. Cinnamodendron from the Antilles, as well as Capsicodendron, South American Cinnamodendron and Pleodendron were not monophyletic. Consequently, Capsicodendron should be included in the South American Cinnamodendron clade and the genus Pleodendron merged with the Cinnamodendron clade from the Antilles. We also found that Warburgia (restricted to mainland eastern Africa) together with the South American Cinnamodendron and Capsicodendron are sister to the Malagasy genus Cinnamosma. In addition to the unexpected geographical relationships, both biogeographic and molecular clock analyses suggest vicariance, extinction, and at least one intercontinental long-distance-dispersal event. Our dating result contrasts previous work on Winteraceae. Diversification of Winteraceae took place in the Paleocene, predating the Canellaceae diversification by 13 MA in the Eocene. The phylogenetic relationships for Canellaceae supported here offer a solid framework for a future taxonomic revision of the Canellaceae.

Genomic repeat abundances contain phylogenetic signal

A large proportion of genomic information, particularly repetitive elements, is usually ignored when researchers are using next-generation sequencing. Here we demonstrate the usefulness of this repetitive fraction in phylogenetic analyses, utilizing comparative graph-based clustering of next-generation sequence reads, which results in abundance estimates of different classes of genomic repeats. Phylogenetic trees are then inferred based on the genome-wide abundance of different repeat types treated as continuously varying characters; such repeats are scattered across chromosomes and in angiosperms can constitute a majority of nuclear genomic DNA. In six diverse examples, five angiosperms and one insect, this method provides generally well-supported relationships at interspecific and intergeneric levels that agree with results from more standard phylogenetic analyses of commonly used markers. We propose that this methodology may prove especially useful in groups where there is little genetic differentiation in standard phylogenetic markers. At the same time as providing data for phylogenetic inference, this method additionally yields a wealth of data for comparative studies of genome evolution.

Comparative transcriptome resources of eleven Primulina species, a group of 'stone plants' from a biodiversity hot spot

The genus Primulina is an emerging model system in studying the drivers and mechanisms of species diversification, for its high species richness and endemism, together with high degree of habitat specialization. In this study, we sequenced transcriptomes for eleven Primulina species across the phylogeny of the genus using the Illumina HiSeq 2000 platform. A total of 336 million clean reads were processed into 355 573 unigenes with a mean length of 1336 bp and an N50 value of 2191 bp after pooling and reassembling twelve individual pre-assembled unigene sets. Of these unigenes, 249 973 (70%) were successfully annotated and 256 601 (72%) were identified as coding sequences (CDSs). We identified a total of 38 279 simple sequence repeats (SSRs) and 367 123 single nucleotide polymorphisms (SNPs). Marker validation assay revealed that 354 (27.3%) of the 1296 SSR and 795 (39.6%) of the 2008 SNP loci showed successful genotyping performance and exhibited expected polymorphism profiles. We screened 834 putative single-copy nuclear genes and proved their high effectiveness in phylogeny construction and estimation of ancestral population parameters. We identified a total of 85 candidate orthologs under positive selection for 46 of the 66 species pairs. This study provided an efficient application of RNA-seq in development of genomic resources for a group of 'stone plants' from south China Karst regions, a biodiversity hot spot of the World. The assembled unigenes with annotations and the massive gene-associated molecular markers would help guide further molecular systematic, population genetic and ecological genomics studies in Primulina and its relatives.

Resolution of deep angiosperm phylogeny using conserved nuclear genes and estimates of early divergence times

Angiosperms are the most successful plants and support human livelihood and ecosystems. Angiosperm phylogeny is the foundation of studies of gene function and phenotypic evolution, divergence time estimation and biogeography. The relationship of the five divergent groups of the Mesangiospermae (~99.95% of extant angiosperms) remains uncertain, with multiple hypotheses reported in the literature. Here transcriptome data sets are obtained from 26 species lacking sequenced genomes, representing each of the five groups: eudicots, monocots, magnoliids, Chloranthaceae and Ceratophyllaceae. Phylogenetic analyses using 59 carefully selected low-copy nuclear genes resulted in highly supported relationships: sisterhood of eudicots and a clade containing Chloranthaceae and Ceratophyllaceae, with magnoliids being the next sister group, followed by monocots. Our topology allows a re-examination of the evolutionary patterns of 110 morphological characters. The molecular clock estimates of Mesangiospermae diversification during the late to middle Jurassic correspond well to the origins of some insects, which may have been a factor facilitating early angiosperm radiation.

Phylogeny and divergence times of gymnosperms inferred from single-copy nuclear genes

Phylogenetic reconstruction is fundamental to study evolutionary biology and historical biogeography. However, there was not a molecular phylogeny of gymnosperms represented by extensive sampling at the genus level, and most published phylogenies of this group were constructed based on cytoplasmic DNA markers and/or the multi-copy nuclear ribosomal DNA. In this study, we use LFY and NLY, two single-copy nuclear genes that originated from an ancient gene duplication in the ancestor of seed plants, to reconstruct the phylogeny and estimate divergence times of gymnosperms based on a complete sampling of extant genera. The results indicate that the combined LFY and NLY coding sequences can resolve interfamilial relationships of gymnosperms and intergeneric relationships of most families. Moreover, the addition of intron sequences can improve the resolution in Podocarpaceae but not in cycads, although divergence times of the cycad genera are similar to or longer than those of the Podocarpaceae genera. Our study strongly supports cycads as the basal-most lineage of gymnosperms rather than sister to Ginkgoaceae, and a sister relationship between Podocarpaceae and Araucariaceae and between Cephalotaxaceae-Taxaceae and Cupressaceae. In addition, intergeneric relationships of some families that were controversial, and the relationships between Taxaceae and Cephalotaxaceae and between conifers and Gnetales are discussed based on the nuclear gene evidence. The molecular dating analysis suggests that drastic extinctions occurred in the early evolution of gymnosperms, and extant coniferous genera in the Northern Hemisphere are older than those in the Southern Hemisphere on average. This study provides an evolutionary framework for future studies on gymnosperms.

Infrageneric phylogeny and temporal divergence of Sorghum (Andropogoneae, Poaceae) based on low-copy nuclear and plastid sequences

The infrageneric phylogeny and temporal divergence of Sorghum were explored in the present study. Sequence data of two low-copy nuclear (LCN) genes, phosphoenolpyruvate carboxylase 4 (Pepc4) and granule-bound starch synthase I (GBSSI), from 79 accessions of Sorghum plus Cleistachne sorghoides together with those from outgroups were used for maximum likelihood (ML) and Bayesian inference (BI) analyses. Bayesian dating based on three plastid DNA markers (ndhA intron, rpl32-trnL, and rps16 intron) was used to estimate the ages of major diversification events in Sorghum. The monophyly of Sorghum plus Cleistachne sorghoides (with the latter nested within Sorghum) was strongly supported by the Pepc4 data using BI analysis, and the monophyly of Sorghum was strongly supported by GBSSI data using both ML and BI analyses. Sorghum was divided into three clades in the Pepc4, GBSSI, and plastid phylograms: the subg. Sorghum lineage; the subg. Parasorghum and Stiposorghum lineage; and the subg. Chaetosorghum and Heterosorghum lineage. Two LCN homoeologous loci of Cleistachne sorghoides were first discovered in the same accession. Sorghum arundinaceum, S. bicolor, S. x drummondii, S. propinquum, and S. virgatum were closely related to S. x almum in the Pepc4, GBSSI, and plastid phylograms, suggesting that they may be potential genome donors to S. almum. Multiple LCN and plastid allelic variants have been identified in S. halepense of subg. Sorghum. The crown ages of Sorghum plus Cleistachne sorghoides and subg. Sorghum are estimated to be 12.7 million years ago (Mya) and 8.6 Mya, respectively. Molecular results support the recognition of three distinct subgenera in Sorghum: subg. Chaetosorghum with two sections, each with a single species, subg. Parasorghum with 17 species, and subg. Sorghum with nine species and we also provide a new nomenclatural combination, Sorghum sorghoides.

Identification of putative orthologous genes for the phylogenetic reconstruction of temperate woody bamboos (Poaceae: Bambusoideae)

The temperate woody bamboos (Arundinarieae) are highly diverse in morphology but lack a substantial amount of genetic variation. The taxonomy of this lineage is intractable, and the relationships within the tribe have not been well resolved. Recent studies indicated that this tribe could have a complex evolutionary history. Although phylogenetic studies of the tribe have been carried out, most of these phylogenetic reconstructions were based on plastid data, which provide lower phylogenetic resolution compared with nuclear data. In this study, we intended to identify a set of desirable nuclear genes for resolving the phylogeny of the temperate woody bamboos. Using two different methodologies, we identified 209 and 916 genes, respectively, as putative single copy orthologous genes. A total of 112 genes was successfully amplified and sequenced by next-generation sequencing technologies in five species sampled from the tribe. As most of the genes exhibited intra-individual allele heterozygotes, we investigated phylogenetic utility by reconstructing the phylogeny based on individual genes. Discordance among gene trees was observed and, to resolve the conflict, we performed a range of analyses using BUCKy and HybTree. While caution should be taken when inferring a phylogeny from multiple conflicting genes, our analysis indicated that 74 of the 112 investigated genes are potential markers for resolving the phylogeny of the temperate woody bamboos.

Conserved genetic regions across angiosperms as tools to develop single-copy nuclear markers in gymnosperms: an example using cycads

Several individuals of the Caribbean Zamia clade and other cycad genera were used to identify single-copy nuclear genes for phylogeographic and phylogenetic studies in Cycadales. Two strategies were employed to select target loci: (i) a tblastX search of Arabidopsis conserved ortholog sequence (COS) set and (ii) a tblastX search of Arabidopsis-Populus-Vitis-Oryza Shared Single-Copy genes (APVO SSC) against the EST Zamia databases in GenBank. From the first strategy, 30 loci were selected, and from the second, 16 loci. In both cases, the matching GenBank accessions of Zamia were used as a query for retrieving highly similar sequences from Cycas, Picea, Pinus species or Ginkgo biloba. After retrieving and aligning all the sequences in each locus, intron predictions were completed to assist in primer design. PCR was carried out in three rounds to detect paralogous loci. A total of 29 loci were successfully amplified as a single band of which 20 were likely single-copy loci. These loci showed different diversity and divergence levels. A preliminary screening allowed us to select 8 promising loci (40S, ATG2, BG, GroES, GTP, LiSH, PEX4 and TR) for the Zamia pumila complex and 4 loci (COS26, GroES, GTP and HTS) for all other cycad genera.

Calcium-dependent protein kinases in plants: evolution, expression and function

Calcium-dependent protein kinases (CPKs) are plant proteins that directly bind calcium ions before phosphorylating substrates involved in metabolism, osmosis, hormone response and stress signaling pathways. CPKs are a large multigene family of proteins that are present in all plants studied to date, as well as in protists, oomycetes and green algae, but are not found in animals and fungi. Despite the increasing evidence of the importance of CPKs in developmental and stress responses from various plants, a comprehensive genome-wide analysis of CPKs from algae to higher plants has not been undertaken. This paper describes the evolution of CPKs from green algae to plants using a broadly sampled phylogenetic analysis and demonstrates the functional diversification of CPKs based on expression and functional studies in different plant species. Our findings reveal that CPK sequence diversification into four major groups occurred in parallel with the terrestrial transition of plants. Despite significant expansion of the CPK gene family during evolution from green algae to higher plants, there is a high level of sequence conservation among CPKs in all plant species. This sequence conservation results in very little correlation between CPK evolutionary groupings and functional diversity, making the search for CPK functional orthologs a challenge.

Widespread and frequent horizontal transfers of transposable elements in plants

Vertical, transgenerational transmission of genetic material occurs through reproduction of living organisms. In addition to vertical inheritance, horizontal gene transfer between reproductively isolated species has recently been shown to be an important, if not dominant, mechanism in the evolution of prokaryotic genomes. In contrast, only a few horizontal transfer (HT) events have been characterized so far in eukaryotes and mainly concern transposable elements (TEs). Whether these are frequent and have a significant impact on genome evolution remains largely unknown. We performed a computational search for highly conserved LTR retrotransposons among 40 sequenced eukaryotic genomes representing the major plant families. We found that 26 genomes (65%) harbor at least one case of horizontal TE transfer (HTT). These transfers concern species as distantly related as palm and grapevine, tomato and bean, or poplar and peach. In total, we identified 32 cases of HTTs, which could translate into more than 2 million among the 13,551 monocot and dicot genera. Moreover, we show that these TEs have remained functional after their transfer, occasionally causing a transpositional burst. This suggests that plants can frequently exchange genetic material through horizontal transfers and that this mechanism may be important in TE-driven genome evolution.

The Amborella genome and the evolution of flowering plants

Amborella trichopoda is strongly supported as the single living species of the sister lineage to all other extant flowering plants, providing a unique reference for inferring the genome content and structure of the most recent common ancestor (MRCA) of living angiosperms. Sequencing the Amborella genome, we identified an ancient genome duplication predating angiosperm diversification, without evidence of subsequent, lineage-specific genome duplications. Comparisons between Amborella and other angiosperms facilitated reconstruction of the ancestral angiosperm gene content and gene order in the MRCA of core eudicots. We identify new gene families, gene duplications, and floral protein-protein interactions that first appeared in the ancestral angiosperm. Transposable elements in Amborella are ancient and highly divergent, with no recent transposon radiations. Population genomic analysis across Amborella's native range in New Caledonia reveals a recent genetic bottleneck and geographic structure with conservation implications.

Single-copy nuclear genes place haustorial Hydnoraceae within piperales and reveal a cretaceous origin of multiple parasitic angiosperm lineages

Extreme haustorial parasites have long captured the interest of naturalists and scientists with their greatly reduced and highly specialized morphology. Along with the reduction or loss of photosynthesis, the plastid genome often decays as photosynthetic genes are released from selective constraint. This makes it challenging to use traditional plastid genes for parasitic plant phylogenetics, and has driven the search for alternative phylogenetic and molecular evolutionary markers. Thus, evolutionary studies, such as molecular clock-based age estimates, are not yet available for all parasitic lineages. In the present study, we extracted 14 nuclear single copy genes (nSCG) from Illumina transcriptome data from one of the “strangest plants in the world”, Hydnora visseri (Hydnoraceae). A ∼15,000 character molecular dataset, based on all three genomic compartments, shows the utility of nSCG for reconstructing phylogenetic relationships in parasitic lineages. A relaxed molecular clock approach with the same multi-locus dataset, revealed an ancient age of ∼91 MYA for Hydnoraceae. We then estimated the stem ages of all independently originated parasitic angiosperm lineages using a published dataset, which also revealed a Cretaceous origin for Balanophoraceae, Cynomoriaceae and Apodanthaceae. With the exception of Santalales, older parasite lineages tend to be more specialized with respect to trophic level and have lower species diversity. We thus propose the “temporal specialization hypothesis” (TSH) implementing multiple independent specialization processes over time during parasitic angiosperm evolution.

Phylogeny of the cycads based on multiple single-copy nuclear genes: congruence of concatenated parsimony, likelihood and species tree inference methods

BACKGROUND AND AIMS

Despite a recent new classification, a stable phylogeny for the cycads has been elusive, particularly regarding resolution of Bowenia, Stangeria and Dioon. In this study, five single-copy nuclear genes (SCNGs) are applied to the phylogeny of the order Cycadales. The specific aim is to evaluate several gene tree-species tree reconciliation approaches for developing an accurate phylogeny of the order, to contrast them with concatenated parsimony analysis and to resolve the erstwhile problematic phylogenetic position of these three genera.

METHODS

DNA sequences of five SCNGs were obtained for 20 cycad species representing all ten genera of Cycadales. These were analysed with parsimony, maximum likelihood (ML) and three Bayesian methods of gene tree-species tree reconciliation, using Cycas as the outgroup. A calibrated date estimation was developed with Bayesian methods, and biogeographic analysis was also conducted.

KEY RESULTS

Concatenated parsimony, ML and three species tree inference methods resolve exactly the same tree topology with high support at most nodes. Dioon and Bowenia are the first and second branches of Cycadales after Cycas, respectively, followed by an encephalartoid clade (Macrozamia-Lepidozamia-Encephalartos), which is sister to a zamioid clade, of which Ceratozamia is the first branch, and in which Stangeria is sister to Microcycas and Zamia.

CONCLUSIONS

A single, well-supported phylogenetic hypothesis of the generic relationships of the Cycadales is presented. However, massive extinction events inferred from the fossil record that eliminated broader ancestral distributions within Zamiaceae compromise accurate optimization of ancestral biogeographical areas for that hypothesis. While major lineages of Cycadales are ancient, crown ages of all modern genera are no older than 12 million years, supporting a recent hypothesis of mostly Miocene radiations. This phylogeny can contribute to an accurate infrafamilial classification of Zamiaceae.

Phylogenomic analyses of nuclear genes reveal the evolutionary relationships within the BEP clade and the evidence of positive selection in Poaceae

BEP clade of the grass family (Poaceae) is composed of three subfamilies, i.e. Bambusoideae, Ehrhartoideae, and Pooideae. Controversies on the phylogenetic relationships among three subfamilies still persist in spite of great efforts. However, previous evidence was mainly provided from plastid genes with only a few nuclear genes utilized. Given different evolutionary histories recorded by plastid and nuclear genes, it is indispensable to uncover their relationships based on nuclear genes. Here, eleven species with whole-sequenced genome and six species with transcriptomic data were included in this study. A total of 121 one-to-one orthologous groups (OGs) were identified and phylogenetic trees were reconstructed by different tree-building methods. Genes which might have undergone positive selection and played important roles in adaptive evolution were also investigated from 314 and 173 one-to-one OGs in two bamboo species and 14 grass species, respectively. Our results support the ((B, P) E) topology with high supporting values. Besides, our findings also indicate that 24 and nine orthologs with statistically significant evidence of positive selection are mainly involved in abiotic and biotic stress response, reproduction and development, plant metabolism and enzyme etc. from two bamboo species and 14 grass species, respectively. In summary, this study demonstrates the power of phylogenomic approach to shed lights on the evolutionary relationships within the BEP clade, and offers valuable insights into adaptive evolution of the grass family.

Advances in the floral structural characterization of the major subclades of Malpighiales, one of the largest orders of flowering plants

BACKGROUND AND AIMS

Malpighiales are one of the largest angiosperm orders and have undergone radical systematic restructuring based on molecular phylogenetic studies. The clade has been recalcitrant to molecular phylogenetic reconstruction, but has become much more resolved at the suprafamilial level. It now contains so many newly identified clades that there is an urgent need for comparative studies to understand their structure, biology and evolution. This is especially true because the order contains a disproportionally large diversity of rain forest species and includes numerous agriculturally important plants. This study is a first broad systematic step in this endeavour. It focuses on a comparative structural overview of the flowers across all recently identified suprafamilial clades of Malpighiales, and points towards areas that desperately need attention.

METHODS

The phylogenetic comparative analysis of floral structure for the order is based on our previously published studies on four suprafamilial clades of Malpighiales, including also four related rosid orders (Celastrales, Crossosomatales, Cucurbitales, Oxalidales). In addition, the results are compiled from a survey of over 3000 publications on macrosystematics, floral structure and embryology across all orders of the core eudicots.

KEY RESULTS

Most new suprafamilial clades within Malpighiales are well supported by floral structural features. Inner morphological structures of the gynoecium (i.e. stigmatic lobes, inner shape of the locules, placentation, presence of obturators) and ovules (i.e. structure of the nucellus, thickness of the integuments, presence of vascular bundles in the integuments, presence of an endothelium in the inner integument) appear to be especially suitable for characterizing suprafamilial clades within Malpighiales.

CONCLUSIONS

Although the current phylogenetic reconstruction of Malpighiales is much improved compared with earlier versions, it is incomplete, and further focused phylogenetic and morphological studies are needed. Once all major subclades of Malpighiales are elucidated, more in-depth studies on promising structural features can be conducted. In addition, once the phylogenetic tree of Malpighiales, including closely related orders, is more fully resolved, character optimization studies will be possible to reconstruct evolution of structural and biological features within the order.