Revisiting the Zingiberales: using multiplexed exon capture to resolve ancient and recent phylogenetic splits in a charismatic plant lineage

Nuclear phylogenomics of angiosperms and insights into their relationships and evolution

Angiosperms (flowering plants) are by far the most diverse land plant group with over 300,000 species. The sudden appearance of diverse angiosperms in the fossil record was referred to by Darwin as the "abominable mystery," hence contributing to the heightened interest in angiosperm evolution. Angiosperms display wide ranges of morphological, physiological, and ecological characters, some of which have probably influenced their species richness. The evolutionary analyses of these characteristics help to address questions of angiosperm diversification and require well resolved phylogeny. Following the great successes of phylogenetic analyses using plastid sequences, dozens to thousands of nuclear genes from next-generation sequencing have been used in angiosperm phylogenomic analyses, providing well resolved phylogenies and new insights into the evolution of angiosperms. In this review we focus on recent nuclear phylogenomic analyses of large angiosperm clades, orders, families, and subdivisions of some families and provide a summarized Nuclear Phylogenetic Tree of Angiosperm Families. The newly established nuclear phylogenetic relationships are highlighted and compared with previous phylogenetic results. The sequenced genomes of Amborella, Nymphaea, Chloranthus, Ceratophyllum, and species of monocots, Magnoliids, and basal eudicots, have facilitated the phylogenomics of relationships among five major angiosperms clades. All but one of the 64 angiosperm orders were included in nuclear phylogenomics with well resolved relationships except the placements of several orders. Most families have been included with robust and highly supported placements, especially for relationships within several large and important orders and families. Additionally, we examine the divergence time estimation and biogeographic analyses of angiosperm on the basis of the nuclear phylogenomic frameworks and discuss the differences compared with previous analyses. Furthermore, we discuss the implications of nuclear phylogenomic analyses on ancestral reconstruction of morphological, physiological, and ecological characters of angiosperm groups, limitations of current nuclear phylogenomic studies, and the taxa that require future attention.

Comparative Chloroplast Genomics of 21 Species in Zingiberales with Implications for Their Phylogenetic Relationships and Molecular Dating

Zingiberales includes eight families and more than 2600 species, with many species having important economic and ecological value. However, the backbone phylogenetic relationships of Zingiberales still remain controversial, as demonstrated in previous studies, and molecular dating based on chloroplast genomes has not been comprehensively studied for the whole order. Herein, 22 complete chloroplast genomes from 21 species in Zingiberales were sequenced, assembled, and analyzed. These 22 genomes displayed typical quadripartite structures, which ranged from 161,303 bp to 163,979 bp in length and contained 111-112 different genes. The genome structures, gene contents, simple sequence repeats, long repeats, and codon usage were highly conserved, with slight differences among these genomes. Further comparative analysis of the 111 complete chloroplast genomes of Zingiberales, including 22 newly sequenced ones and the remaining ones from the national center for biotechnology information (NCBI) database, identified three highly divergent regions comprising ccsA, psaC, and psaC-ndhE. Maximum likelihood and Bayesian inference phylogenetic analyses based on chloroplast genome sequences found identical topological structures and identified a strongly supported backbone of phylogenetic relationships. Cannaceae was sister to Marantaceae, forming a clade that was collectively sister to the clade of (Costaceae, Zingiberaceae) with strong support (bootstrap (BS) = 100%, and posterior probability (PP) = 0.99-1.0); Heliconiaceae was sister to the clade of (Lowiaceae, Strelitziaceae), then collectively sister to Musaceae with strong support (BS = 94-100%, and PP = 0.93-1.0); the clade of ((Cannaceae, Marantaceae), (Costaceae, Zingiberaceae)) was sister to the clade of (Musaceae, (Heliconiaceae, (Lowiaceae, Strelitziaceae))) with robust support (BS = 100%, and PP = 1.0). The results of divergence time estimation of Zingiberales indicated that the crown node of Zingiberales occurred approximately 85.0 Mya (95% highest posterior density (HPD) = 81.6-89.3 million years ago (Mya)), with major family-level lineages becoming from 46.8 to 80.5 Mya. These findings proved that chloroplast genomes could contribute to the study of phylogenetic relationships and molecular dating in Zingiberales, as well as provide potential molecular markers for further taxonomic and phylogenetic studies of Zingiberales.

A near-complete genome assembly of Thalia dealbata Fraser (Marantaceae)

This study presents a chromosome-level, near-complete genome assembly of Thalia dealbata (Marantaceae), a typical emergent wetland plant with high ornamental and environmental value. Based on 36.99 Gb PacBio HiFi reads and 39.44 Gb Hi-C reads, we obtained a 255.05 Mb assembly, of which 251.92 Mb (98.77%) were anchored into eight pseudo-chromosomes. Five pseudo-chromosomes were completely assembled, and the other three had one to two gaps. The final assembly had a high contig N50 value (29.80 Mb) and benchmarking universal single-copy orthologs (BUSCO) recovery score (97.52%). The T. dealbata genome had 100.35 Mb repeat sequences, 24,780 protein-coding genes, and 13,679 non-coding RNAs. Phylogenetic analysis revealed that T. dealbata was closest to Zingiber officinale, whose divergence time was approximately 55.41 million years ago. In addition, 48 and 52 significantly expanded and contracted gene families were identified within the T. dealbata genome. Moreover, 309 gene families were specific to T. dealbata, and 1,017 genes were positively selected. The T. dealbata genome reported in this study provides a valuable genomic resource for further research on wetland plant adaptation and the genome evolution dynamics. This genome is also beneficial for the comparative genomics of Zingiberales species and flowering plants.

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.

The banana genome hub: a community database for genomics in the Musaceae

The Banana Genome Hub provides centralized access for genome assemblies, annotations, and the extensive related omics resources available for bananas and banana relatives. A series of tools and unique interfaces are implemented to harness the potential of genomics in bananas, leveraging the power of comparative analysis, while recognizing the differences between datasets. Besides effective genomic tools like BLAST and the JBrowse genome browser, additional interfaces enable advanced gene search and gene family analyses including multiple alignments and phylogenies. A synteny viewer enables the comparison of genome structures between chromosome-scale assemblies. Interfaces for differential expression analyses, metabolic pathways and GO enrichment were also added. A catalogue of variants spanning the banana diversity is made available for exploration, filtering, and export to a wide variety of software. Furthermore, we implemented new ways to graphically explore gene presence-absence in pangenomes as well as genome ancestry mosaics for cultivated bananas. Besides, to guide the community in future sequencing efforts, we provide recommendations for nomenclature of locus tags and a curated list of public genomic resources (assemblies, resequencing, high density genotyping) and upcoming resources-planned, ongoing or not yet public. The Banana Genome Hub aims at supporting the banana scientific community for basic, translational, and applied research and can be accessed at https://banana-genome-hub.southgreen.fr.

Correlation between Inflorescence Architecture and Floral Asymmetry-Evidence from Aberrant Flowers in Canna L. (Cannaceae)

Floral symmetry studies often focus on the development of monosymmetric and polysymmetric flowers, whereas asymmetric flowers and their position and function within the inflorescence structure are largely neglected. Cannaceae is one of the few families that possesses truly asymmetric flowers, serving as a model to study the characters and mechanisms involved in the development of floral asymmetry and its context within the developing and mature inflorescence. In this study, inflorescence structure and floral morphology of normal asymmetric flowers and 16 aberrant flower collections from Canna indica L. and C. glauca L. were photographed, analyzed, and compared with attention to stamen petaloidy, floral symmetry, and inflorescence branching patterns anterior and posterior to the aberrant flower. In comparison with normal flowers, the aberrant flowers are arranged into abnormal partial florescences, and vary in floral symmetry, orientation, and degree of androecial petaloidy. The appendage of the fertile stamen is universally located distal from the higher order bract, indicating an underlying influence of inflorescence architecture. A synthetic model is proposed to explain the relationship between floral symmetry and inflorescence structure. Data from the observation of aberrant phenotypes strongly support the hypothesis that irregular petaloidy of the stamens is correlated with an asymmetric morphogenetic field within the inflorescence that contributes to the overall floral asymmetry in Canna flowers.

The genetic mechanisms underlying the convergent evolution of pollination syndromes in the Neotropical radiation of Costus L

Selection together with variation in floral traits can act to mold floral form, often driven by a plant's predominant or most effective pollinators. To investigate the evolution of traits associated with pollination, we developed a phylogenetic framework for evaluating tempo and mode of pollination shifts across the genus Costus L., known for its evolutionary toggle between traits related to bee and bird pollination. Using a target enrichment approach, we obtained 957 loci for 171 accessions to expand the phylogenetic sampling of Neotropical Costus. In addition, we performed whole genome resequencing for a subset of 20 closely related species with contrasting pollination syndromes. For each of these 20 genomes, a high-quality assembled transcriptome was used as reference for consensus calling of candidate loci hypothesized to be associated with pollination-related traits of interest. To test for the role these candidate genes may play in evolutionary shifts in pollinators, signatures of selection were estimated as dN/dS across the identified candidate loci. We obtained a well-resolved phylogeny for Neotropical Costus despite conflict among gene trees that provide evidence of incomplete lineage sorting and/or reticulation. The overall topology and the network of genome-wide single nucleotide polymorphisms (SNPs) indicate that multiple shifts in pollination strategy have occurred across Costus, while also suggesting the presence of previously undetected signatures of hybridization between distantly related taxa. Traits related to pollination syndromes are strongly correlated and have been gained and lost in concert several times throughout the evolution of the genus. The presence of bract appendages is correlated with two traits associated with defenses against herbivory. Although labellum shape is strongly correlated with overall pollination syndrome, we found no significant impact of labellum shape on diversification rates. Evidence suggests an interplay of pollination success with other selective pressures shaping the evolution of the Costus inflorescence. Although most of the loci used for phylogenetic inference appear to be under purifying selection, many candidate genes associated with functional traits show evidence of being under positive selection. Together these results indicate an interplay of phylogenetic history with adaptive evolution leading to the diversification of pollination-associated traits in Neotropical Costus.

distAngsd: Fast and accurate inference of genetic distances for Next Generation Sequencing data

Commonly used methods for inferring phylogenies were designed before the emergence of high throughput sequencing and can generally not accommodate the challenges associated with noisy, diploid sequencing data. In many applications, diploid genomes are still treated as haploid through the use of ambiguity characters; while the uncertainty in genotype calling - arising as a consequence of the sequencing technology - is ignored. In order to address this problem we describe two new probabilistic approaches for estimating genetic distances: distAngsd-geno and distAngsd-nuc, both implemented in a software suite named distAngsd. These methods are specifically designed for next generation sequencing data, utilize the full information from the data, and take uncertainty in genotype calling into account. Through extensive simulations, we show that these new methods are markedly more accurate and have more stable statistical behaviors than other currently available methods for estimating genetic distances - even for very low depth data with high error rates.

Comparative genome sequence and phylogenetic analysis of chloroplast for evolutionary relationship among Pinus species

Genus Pinus is a widely dispersed genus of conifer plants in the Northern Hemisphere. However, the inadequate accessibility of genomic knowledge limits our understanding of molecular phylogeny and evolution of Pinus species. In this study, the evolutionary features of complete plastid genome and the phylogeny of the Pinus genus were studied. A total of thirteen divergent hotspot regions (trnk-UUU, matK, trnQ-UUG, atpF, atpH, rpoC1, rpoC2, rpoB, ycf2, ycf1, trnD-GUC, trnY-GUA, and trnH-GUG) were identified that would be utilized as possible genetic markers for determination of phylogeny and population genetics analysis of Pinus species. Furthermore, seven genes (petD, psaI, psaM, matK, rps18, ycf1, and ycf2) with positive selection site in Pinus species were identified. Based on the whole genome this phylogenetic study showed that twenty-four Pinus species form a significant genealogical clade. Divergence time showed that the Pinus species originated about 100 million years ago (MYA) (95% HPD, 101.76.35–109.79 MYA), in lateral stages of Cretaceous. Moreover, two of the subgenera are consequently originated in 85.05 MYA (95% HPD, 81.04–88.02 MYA). This study provides a phylogenetic relationship and a chronological framework for the future study of the molecular evolution of the Pinus species.

Whence Came These Plants Most Foul? Phylogenomics and Biogeography of Lowiaceae (Zingiberales)

Lowiaceae (order Zingiberales) is a small family of forest herbs in Southeast Asia. All species belong to the genus Orchidantha. They are known for possessing orchid-like flowers that are smelly, apparently mimicking dead animals, feces, or mushrooms. Little is known of the biogeographic patterns or character evolution of the family. We sampled the family extensively, including many recently discovered species, and reconstructed the phylogeny of the family using HybSeq with Lowiaceae-specific RNA baits. Our phylogenetic reconstructions confirm that the family is most closely related to Strelitziaceae, and that species with dark, foul-smelling flowers form a grade in which a clade of species with paler flowers are embedded. The pale-flowered species produce a distinct odor, resembling edible mushrooms. Apart from a single species, the species from Borneo form a clade, and the same is true for Indochinese species. The remaining species form a more widespread clade. A biogeographic analysis shows that the distribution of Lowiaceae can explained by vicariance and gradual dispersal from a shared ancestral range of Borneo and Indochina. There is no evidence of long-distance dispersal, only a later extension in distribution to Peninsular Malaysia which coincides with the presence of a land bridge. Different directions of spread are possible, but none require long-distance dispersal. The results are consistent with the geological history of Southeast Asia. In particular, the relatively early isolation between Indochina and Borneo could be explained by the presence of a sea barrier that developed 10–15 MYA, and the continuous movement of plant species between Borneo and Peninsular Malaysia could be explained by a land bridge that existed until c. 5 MYA. The lack of an extensive land bridge with a suitable habitat may explain the absence of this genus from Sumatra and other Indonesian islands aside from Borneo. The strict reliance on a continuous habitat for the range expansion of Lowiaceae can be explained by their fruits and seeds, which lack obvious adaptations for long-distance dispersal. The inability to disperse to new areas may also explain why the extant species have very restricted distributions.

Aiming off the target: recycling target capture sequencing reads for investigating repetitive DNA

BACKGROUND AND AIMS

With the advance of high-throughput sequencing, reduced-representation methods such as target capture sequencing (TCS) emerged as cost-efficient ways of gathering genomic information, particularly from coding regions. As the off-target reads from such sequencing are expected to be similar to genome skimming (GS), we assessed the quality of repeat characterization in plant genomes using these data.

METHODS

Repeat composition obtained from TCS datasets of five Rhynchospora (Cyperaceae) species were compared with GS data from the same taxa. In addition, a FISH probe was designed based on the most abundant satellite found in the TCS dataset of Rhynchospora cephalotes. Finally, repeat-based phylogenies of the five Rhynchospora species were constructed based on the GS and TCS datasets and the topologies were compared with a gene-alignment-based phylogenetic tree.

KEY RESULTS

All the major repetitive DNA families were identified in TCS, including repeats that showed abundances as low as 0.01 % in the GS data. Rank correlations between GS and TCS repeat abundances were moderately high (r = 0.58-0.85), increasing after filtering out the targeted loci from the raw TCS reads (r = 0.66-0.92). Repeat data obtained by TCS were also reliable in developing a cytogenetic probe of a new variant of the holocentromeric satellite Tyba. Repeat-based phylogenies from TCS data were congruent with those obtained from GS data and the gene-alignment tree.

CONCLUSIONS

Our results show that off-target TCS reads can be recycled to identify repeats for cyto- and phylogenomic investigations. Given the growing availability of TCS reads, driven by global phylogenomic projects, our strategy represents a way to recycle genomic data and contribute to a better characterization of plant biodiversity.

Taxon-specific or universal? Using target capture to study the evolutionary history of rapid radiations

Target capture has emerged as an important tool for phylogenetics and population genetics in nonmodel taxa. Whereas developing taxon‐specific capture probes requires sustained efforts, available universal kits may have a lower power to reconstruct relationships at shallow phylogenetic scales and within rapidly radiating clades. We present here a newly developed target capture set for Bromeliaceae, a large and ecologically diverse plant family with highly variable diversification rates. The set targets 1776 coding regions, including genes putatively involved in key innovations, with the aim to empower testing of a wide range of evolutionary hypotheses. We compare the relative power of this taxon‐specific set, Bromeliad1776, to the universal Angiosperms353 kit. The taxon‐specific set results in higher enrichment success across the entire family; however, the overall performance of both kits to reconstruct phylogenetic trees is relatively comparable, highlighting the vast potential of universal kits for resolving evolutionary relationships. For more detailed phylogenetic or population genetic analyses, for example the exploration of gene tree concordance, nucleotide diversity or population structure, the taxon‐specific capture set presents clear benefits. We discuss the potential lessons that this comparative study provides for future phylogenetic and population genetic investigations, in particular for the study of evolutionary radiations.

New genetic markers for Sapotaceae phylogenomics: More than 600 nuclear genes applicable from family to population levels

Some tropical plant families, such as the Sapotaceae, have a complex taxonomy, which can be resolved using Next Generation Sequencing (NGS). For most groups however, methodological protocols are still missing. Here we identified 531 monocopy genes and 227 Short Tandem Repeats (STR) markers and tested them on Sapotaceae using target capture and NGS. The probes were designed using two genome skimming samples from Capurodendron delphinense and Bemangidia lowryi, both from the Tseboneae tribe, as well as the published Manilkara zapota transcriptome from the Sapotoideae tribe. We combined our probes with 261 additional ones previously published and designed for the entire angiosperm group. On a total of 792 low-copy genes, 638 showed no signs of paralogy and were used to build a phylogeny of the family with 231 individuals from all main lineages. A highly supported topology was obtained at high taxonomic ranks but also at the species level. This phylogeny revealed the existence of more than 20 putative new species. Single nucleotide polymorphisms (SNPs) extracted from the 638 genes were able to distinguish lineages within a species complex and to highlight geographical structuration. STR were recovered efficiently for the species used as reference (C. delphinense) but the recovery rate decreased dramatically with the phylogenetic distance to the focal species. Altogether, the new loci will help reaching a sound taxonomic understanding of the family Sapotaceae for which many circumscriptions and relationships are still debated, at the species, genus and tribe levels.

Hybrid capture of 964 nuclear genes resolves evolutionary relationships in the mimosoid legumes and reveals the polytomous origins of a large pantropical radiation

PREMISE

Targeted enrichment methods facilitate sequencing of hundreds of nuclear loci to enhance phylogenetic resolution and elucidate why some parts of the "tree of life" are difficult (if not impossible) to resolve. The mimosoid legumes are a prominent pantropical clade of ~3300 species of woody angiosperms for which previous phylogenies have shown extensive lack of resolution, especially among the species-rich and taxonomically challenging ingoids.

METHODS

We generated transcriptomes to select low-copy nuclear genes, enrich these via hybrid capture for representative species of most mimosoid genera, and analyze the resulting data using de novo assembly and various phylogenomic tools for species tree inference. We also evaluate gene tree support and conflict for key internodes and use phylogenetic network analysis to investigate phylogenetic signal across the ingoids.

RESULTS

Our selection of 964 nuclear genes greatly improves phylogenetic resolution across the mimosoid phylogeny and shows that the ingoid clade can be resolved into several well-supported clades. However, nearly all loci show lack of phylogenetic signal for some of the deeper internodes within the ingoids.

CONCLUSIONS

Lack of resolution in the ingoid clade is most likely the result of hyperfast diversification, potentially causing a hard polytomy of six or seven lineages. The gene set for targeted sequencing presented here offers great potential to further enhance the phylogeny of mimosoids and the wider Caesalpinioideae with denser taxon sampling, to provide a framework for taxonomic reclassification, and to study the ingoid radiation.

Unraveling the Spiraling Radiation: A Phylogenomic Analysis of Neotropical Costus L

The family of pantropical spiral gingers (Costaceae Nakai; c. 125 spp.) can be used as a model to enhance our understanding of the mechanisms underlying Neotropical diversity. Costaceae has higher taxonomic diversity in South and Central America (c. 72 Neotropical species, c. 30 African, c. 23 Southeast Asian), particularly due to a radiation of Neotropical species of the genus Costus L. (c. 57 spp.). However, a well-supported phylogeny of the Neotropical spiral gingers including thorough sampling of proposed species encompassing their full morphologic and geographic variation is lacking, partly due to poor resolution recovered in previous analyses using a small sampling of loci. Here we use a phylogenomic approach to estimate the phylogeny of a sample of Neotropical Costus species using a targeted enrichment approach. Baits were designed to capture conserved elements' variable at the species level using available genomic sequences of Costus species and relatives. We obtained 832 loci (generating 791,954 aligned base pairs and 31,142 parsimony informative sites) for samples that encompassed the geographical and/or morphological diversity of some recognized species. Higher support values that improve the results of previous studies were obtained when including all the available loci, even those producing unresolved gene trees and having a low proportion of variable sites. Concatenation and coalescent-based species trees methods converge in almost the same topology suggesting a robust estimation of the relationships, even under the high levels of gene tree conflict presented here. The bait set design here presented made inferring a robust phylogeny to test taxonomic hypotheses possible and will improve our understanding of the origins of the charismatic diversity of the Neotropical spiral gingers.

Can the anatomy of abnormal flowers elucidate relationships of the androecial members in the ginger (Zingiberaceae)?

BACKGROUND

Interpretation of the floral structure of Zingiberaceae has long concentrated on the relationships of the androecial members. It suggested that labellum is composed of two structures rather than three or five, and glands are interpreted either as gynoecial part or as androecial members.

METHODS

Serial sections were used to observe the vasculature of normal and two-staminate flowers in Alpinia intermedia 'shengzhen'. Floral diagrams were drawn to interpret the morphological category of the floral organs and the relationships of the androecial members. Androecial vascular bundles were associated with carpellary dorsal bundles (CDBs) and parietal bundles (PBs) in a Zingiberales phylogeny setting using ancestral state reconstruction.

RESULTS

Anatomical observations demonstrate that the fertile stamen(s) incorporate parietal bundles both in normal and two-staminate flowers. The three appendages represent the three members of the outer whorl of the androecium, while the labellum represents the inner whorl of the androecium in the two-staminate flower. Reconstruction of the origin of the vascular system in the androecium suggests that the outer whorl of androecium receives its vascular supply from the CDBs, and the inner whorl of androecium receives from the PBs in both the basal banana group and the more derived ginger clade.

CONCLUSIONS

The present study adds to a growing body of literature suggesting that the anatomy of abnormal flowers may not provide enough evidence for elucidating the relationships of the androecial members, and help us to better understand how the vascular system is constructed during the androecial petaloidy evolution.

Target Nuclear and Off-Target Plastid Hybrid Enrichment Data Inform a Range of Evolutionary Depths in the Orchid Genus Epidendrum

Universal angiosperm enrichment probe sets designed to enrich hundreds of putatively orthologous nuclear single-copy loci are increasingly being applied to infer phylogenetic relationships of different lineages of angiosperms at a range of evolutionary depths. Studies applying such probe sets have focused on testing the universality and performance of the target nuclear loci, but they have not taken advantage of off-target data from other genome compartments generated alongside the nuclear loci. Here we do so to infer phylogenetic relationships in the orchid genus Epidendrum and closely related genera of subtribe Laeliinae. Our aims are to: 1) test the technical viability of applying the plant anchored hybrid enrichment (AHE) method (Angiosperm v.1 probe kit) to our focal group, 2) mine plastid protein coding genes from off-target reads; and 3) evaluate the performance of the target nuclear and off-target plastid loci in resolving and supporting phylogenetic relationships along a range of taxonomical depths. Phylogenetic relationships were inferred from the nuclear data set through coalescent summary and site-based methods, whereas plastid loci were analyzed in a concatenated partitioned matrix under maximum likelihood. The usefulness of target and flanking non-target nuclear regions and plastid loci was assessed through the estimation of their phylogenetic informativeness. Our study successfully applied the plant AHE probe kit to Epidendrum, supporting the universality of this kit in angiosperms. Moreover, it demonstrated the feasibility of mining plastome loci from off-target reads generated with the Angiosperm v.1 probe kit to obtain additional, uniparentally inherited sequence data at no extra sequencing cost. Our analyses detected some strongly supported incongruences between nuclear and plastid data sets at shallow divergences, an indication of potential lineage sorting, hybridization, or introgression events in the group. Lastly, we found that the per site phylogenetic informativeness of the ycf1 plastid gene surpasses that of all other plastid genes and several nuclear loci, making it an excellent candidate for assessing phylogenetic relationships at medium to low taxonomic levels in orchids.

Sinopse de Marantaceae no estado de Pernambuco, Brasil

Resumo O território de Pernambuco engloba partes da Caatinga e da Mata Atlântica, abrigando áreas de extrema importância biológica, ricas em espécies de Marantaceae. Este estudo apresenta uma sinopse taxonômica de Marantaceae ocorrentes no estado, sendo baseada na análise morfológica de materiais herborizados depositados em herbários nacionais e de espécimes provenientes de excursões de campo realizadas no período de 2013 a 2018. Foram encontradas 26 espécies pertencentes a 10 gêneros, sendo Maranta (sete spp.) e Goeppertia (cinco spp.) os mais representativos. Dentre as espécies, Ctenanthe casupoides e Hylaeanthe hexantha são novos registros para o estado, e Goeppertia yoshida-arnsiae, G. widgrenii, e Maranta gigantea são espécies que se encontram sob algum grau de ameaça. Dez espécies são endêmicas da Mata Atlântica e quatro restritas a região Nordeste. Com relação à distribuição local, 22 spp. (85%) ocorrem na Floresta Estacional Semidecidual de Terras Baixas. São apresentadas descrições, chaves de identificação, comentários sobre distribuição geográfica e habitats, e mapas de distribuição local, juntamente com prancha de fotos das espécies.

The complete chloroplast genome of Curcuma longa L. (Zingiberaceae)

Curcuma longa L. is a widely distributed medicinal plant in Southeast Asia, especially in China. In this study, we reported and characterized the complete chloroplast genome sequence of C. longa. The chloroplast genome was determined to be 162,255 bp in length. It contained large single-copy (LSC) and small single-copy (SSC) regions of 86,975 bp and 15,776 bp, respectively, which were separated by a pair of 29,752 bp inverted repeat (IR) regions. The genome is predicted to contain 133 genes, including 87 protein-coding genes, 38 tRNA genes, and eight rRNA genes. The overall GC content of the genome is 36.2%. A phylogenetic tree reconstructed by 23 chloroplast genomes reveals that C. longa is mostly related to Curcuma roscoeana, with bootstrap support values of 100%.

Targeted Capture of Hundreds of Nuclear Genes Unravels Phylogenetic Relationships of the Diverse Neotropical Palm Tribe Geonomateae

The tribe Geonomateae is a widely distributed group of 103 species of Neotropical palms which contains six ecologically important understory or subcanopy genera. Although it has been the focus of many studies, our understanding of the evolutionary history of this group, and in particular of the taxonomically complex genus Geonoma, is far from complete due to a lack of molecular data. Specifically, the previous Sanger sequencing-based studies used a few informative characters and partial sampling. To overcome these limitations, we used a recently developed Arecaceae-specific target capture bait set to undertake a phylogenomic analysis of the tribe Geonomateae. We sequenced 3,988 genomic regions for 85% of the species of the tribe, including 84% of the species of the largest genus, Geonoma. Phylogenetic relationships were inferred using both concatenation and coalescent methods. Overall, our phylogenetic tree is highly supported and congruent with taxonomic delimitations although several morphological taxa were revealed to be non-monophyletic. It is the first time that such a large genomic dataset is provided for an entire tribe within the Arecaceae. Our study lays the groundwork not only for detailed macro- and micro-evolutionary studies within the group, but also sets a workflow for understanding other species complexes across the tree of life.

A customized nuclear target enrichment approach for developing a phylogenomic baseline for Dioscorea yams (Dioscoreaceae)

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.

Leaf Transcriptome Assembly of Protium copal (Burseraceae) and Annotation of Terpene Biosynthetic Genes

Plants in the Burseraceae are globally recognized for producing resins and essential oils with medicinal properties and have economic value. In addition, most of the aromatic and non-aromatic components of Burseraceae resins are derived from a variety of terpene and terpenoid chemicals. Although terpene genes have been identified in model plant crops (e.g., Citrus, Arabidopsis), very few genomic resources are available for non-model groups, including the highly diverse Burseraceae family. Here we report the assembly of a leaf transcriptome of Protium copal, an aromatic tree that has a large distribution in Central America, describe the functional annotation of putative terpene biosynthetic genes and compare terpene biosynthetic genes found in P. copal with those identified in other Burseraceae taxa. The genomic resources of Protium copal can be used to generate novel sequencing markers for population genetics and comparative phylogenetic studies, and to investigate the diversity and evolution of terpene genes in the Burseraceae.

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

Comparative analysis of whole flower transcriptomes in the Zingiberales

The advancement of next generation sequencing technologies (NGS) has revolutionized our ability to generate large quantities of data at a genomic scale. Despite great challenges, these new sequencing technologies have empowered scientists to explore various relevant biological questions on non-model organisms, even in the absence of a complete sequenced reference genome. Here, we analyzed whole flower transcriptome libraries from exemplar species across the monocot order Zingiberales, using a comparative approach in order to gain insight into the evolution of the molecular mechanisms underlying flower development in the group. We identified 4,153 coding genes shared by all floral transcriptomes analyzed, and 1,748 genes that are only retrieved in the Zingiberales. We also identified 666 genes that are unique to the ginger lineage, and 2,001 that are only found in the banana group, while in the outgroup species Dichorisandra thyrsiflora J.C. Mikan (Commelinaceae) we retrieved 2,686 unique genes. It is possible that some of these genes underlie lineage-specific molecular mechanisms of floral diversification. We further discuss the nature of these lineage-specific datasets, emphasizing conserved and unique molecular processes with special emphasis in the Zingiberales. We also briefly discuss the strengths and shortcomings of de novo assembly for the study of developmental processes across divergent taxa from a particular order. Although this comparison is based exclusively on coding genes, with particular emphasis in transcription factors, we believe that the careful study of other regulatory mechanisms, such as non-coding RNAs, might reveal new levels of complexity, which were not explored in this work.

Building the monocot tree of death: Progress and challenges emerging from the macrofossil-rich Zingiberales

PREMISE OF THE STUDY

Inclusion of fossils in phylogenetic analyses is necessary in order to construct a comprehensive "tree of death" and elucidate evolutionary history of taxa; however, such incorporation of fossils in phylogenetic reconstruction is dependent on the availability and interpretation of extensive morphological data. Here, the Zingiberales, whose familial relationships have been difficult to resolve with high support, are used as a case study to illustrate the importance of including fossil taxa in systematic studies.

METHODS

Eight fossil taxa and 43 extant Zingiberales were coded for 39 morphological seed characters, and these data were concatenated with previously published molecular sequence data for analysis in the program MrBayes.

KEY RESULTS

Ensete oregonense is confirmed to be part of Musaceae, and the other seven fossils group with Zingiberaceae. There is strong support for Spirematospermum friedrichii, Spirematospermum sp. 'Goth', S. wetzleri, and Striatornata sanantoniensis in crown Zingiberaceae while "Musa" cardiosperma, Spirematospermum chandlerae, and Tricostatocarpon silvapinedae are best considered stem Zingiberaceae. Inclusion of fossils explains how different topologies from morphological and molecular data sets is due to shared plesiomorphic characters shared by Musaceae, Zingiberaceae, and Costaceae, and most of the fossils.

CONCLUSIONS

Inclusion of eight fossil taxa expands the Zingiberales tree and helps explain the difficulty in resolving relationships. Inclusion of fossils was possible in part due to a large morphological data set built using nondestructive microcomputed tomography data. Collaboration between paleo- and neobotanists and technology such as microcomputed tomography will help to build the tree of death and ultimately improve our understanding of the evolutionary history of monocots.

Exploring data processing strategies in NGS target enrichment to disentangle radiations in the tribe Cardueae (Compositae)

Target enrichment is a cost-effective sequencing technique that holds promise for elucidating evolutionary relationships in fast-evolving lineages. However, potential biases and impact of bioinformatic sequence treatments in phylogenetic inference have not been thoroughly explored yet. Here, we investigate this issue with an ultimate goal to shed light into a highly diversified group of Compositae (Asteraceae) constituted by four main genera: Arctium, Cousinia, Saussurea, and Jurinea. Specifically, we compared sequence data extraction methods implemented in two easy-to-use workflows, PHYLUCE and HybPiper, and assessed the impact of two filtering practices intended to reduce phylogenetic noise. In addition, we compared two phylogenetic inference methods: (1) the concatenation approach, in which all loci were concatenated in a supermatrix; and (2) the coalescence approach, in which gene trees were produced independently and then used to construct a species tree under coalescence assumptions. Here we confirm the usefulness of the set of 1061 COS targets (a nuclear conserved orthology loci set developed for the Compositae) across a variety of taxonomic levels. Intergeneric relationships were completely resolved: there are two sister groups, Arctium-Cousinia and Saussurea-Jurinea, which are in agreement with a morphological hypothesis. Intrageneric relationships among species of Arctium, Cousinia, and Saussurea are also well defined. Conversely, conflicting species relationships remain for Jurinea. Methodological choices significantly affected phylogenies in terms of topology, branch length, and support. Across all analyses, the phylogeny obtained using HybPiper and the strictest scheme of removing fast-evolving sites was estimated as the optimal. Regarding methodological choices, we conclude that: (1) trees obtained under the coalescence approach are topologically more congruent between them than those inferred using the concatenation approach; (2) refining treatments only improved support values under the concatenation approach; and (3) branch support values are maximized when fast-evolving sites are removed in the concatenation approach, and when a higher number of loci is analyzed in the coalescence approach.

Genome sequence data from 17 accessions of Ensete ventricosum, a staple food crop for millions in Ethiopia

We present raw sequence reads and genome assemblies derived from 17 accessions of the Ethiopian orphan crop plant enset (Ensete ventricosum (Welw.) Cheesman) using the Illumina HiSeq and MiSeq platforms. Also presented is a catalogue of single-nucleotide polymorphisms inferred from the sequence data at an average density of approximately one per kilobase of genomic DNA.

Phylogenomic analyses reveal a deep history of hybridization and polyploidy in the Neotropical genus Lachemilla (Rosaceae)

Hybridization, incomplete lineage sorting, and phylogenetic error produce similar incongruence patterns, representing a great challenge for phylogenetic reconstruction. Here, we use sequence capture data and multiple species tree and species network approaches to resolve the backbone phylogeny of the Neotropical genus Lachemilla, while distinguishing among sources of incongruence. We used 396 nuclear loci and nearly complete plastome sequences from 27 species to clarify the relationships among the major groups of Lachemilla, and explored multiple sources of conflict between gene trees and species trees inferred with a plurality of approaches. All phylogenetic methods recovered the four major groups previously proposed for Lachemilla, but species tree methods recovered different topologies for relationships between these four clades. Species network analyses revealed that one major clade, Orbiculate, is likely of ancient hybrid origin, representing one of the main sources of incongruence among the species trees. Additionally, we found evidence for a potential whole genome duplication event shared by Lachemilla and allied genera. Lachemilla shows clear evidence of ancient and recent hybridization throughout the evolutionary history of the group. Also, we show the necessity to use phylogenetic network approaches that can simultaneously accommodate incomplete lineage sorting and gene flow when studying groups that show patterns of reticulation.

Dense infraspecific sampling reveals rapid and independent trajectories of plastome degradation in a heterotrophic orchid complex

Heterotrophic plants provide excellent opportunities to study the effects of altered selective regimes on genome evolution. Plastid genome (plastome) studies in heterotrophic plants are often based on one or a few highly divergent species or sequences as representatives of an entire lineage, thus missing important evolutionary-transitory events. Here, we present the first infraspecific analysis of plastome evolution in any heterotrophic plant. By combining genome skimming and targeted sequence capture, we address hypotheses on the degree and rate of plastome degradation in a complex of leafless orchids (Corallorhiza striata) across its geographic range. Plastomes provide strong support for relationships and evidence of reciprocal monophyly between C. involuta and the endangered C. bentleyi. Plastome degradation is extensive, occurring rapidly over a few million years, with evidence of differing rates of genomic change among the two principal clades of the complex. Genome skimming and targeted sequence capture differ widely in coverage depth overall, with depth in targeted sequence capture datasets varying immensely across the plastome as a function of GC content. These findings will help to fill a knowledge gap in models of heterotrophic plastid genome evolution, and have implications for future studies in heterotrophs.

Flexible colour patterns obscure identification and mimicry in Indo-Pacific Chromodoris nudibranchs (Gastropoda: Chromodorididae)

Chromodoris is a genus of colourful nudibranchs that feed on sponges and is found across the Indo-Pacific. While this was once the most diverse chromodorid genus, recent work has shown that the genus should be restricted to a monophyletic lineage that contains only 22 species, all of which exhibit black pigmentation and planar spawning behaviour. Earlier phylogenies of this group are poorly resolved and thus additional work is needed to clarify species boundaries within Chromodoris. This study presents a maximum-likelihood phylogeny based on mitochondrial loci (COI, 16S) for 345 Chromodoris specimens, including data from 323 new specimens and 22 from GenBank, from across the Indo-Pacific. Species hypotheses and phylogenetic analysis uncovered 39 taxa in total containing 18 undescribed species, with only five of 39 taxa showing stable colour patterns and distinct morphotypes. This study also presents the first evidence for regional mimicry in this genus, with C. colemani and C. joshi displaying geographically-based variation in colour patterns which appear to match locally abundant congenerics, highlighting the flexibility of these colour patterns in Chromodoris nudibranchs. The current phylogeny contains short branch lengths, polytomies and poor support at interior nodes, which is indicative of a recent radiation. As such, future work will employ a transcriptome-based exon capture approach for resolving the phylogeny of this group. In all, this study included 21 of the 22 described species in the Chromodoris sensu stricto group with broad sampling coverage from across the Indo-Pacific, constituting the most comprehensive sampling of this group to date. This work highlights several cases of undocumented diversity, ultimately expanding our knowledge of species boundaries in this group, while also demonstrating the limitations of colour patterns for species identification in this genus.

Allele phasing is critical to revealing a shared allopolyploid origin of Medicago arborea and M. strasseri (Fabaceae)

BACKGROUND

Whole genome duplication plays a central role in plant evolution. There are two main classes of polyploid formation: autopolyploids which arise within one species by doubling of similar homologous genomes; in contrast, allopolyploidy (hybrid polyploidy) arise via hybridization and subsequent doubling of nonhomologous (homoeologous) genomes. The distinction between polyploid origins can be made using gene phylogenies, if alleles from each genome can be correctly retrieved. We examined whether two closely related tetraploid Mediterranean shrubs (Medicago arborea and M. strasseri) have an allopolyploid origin - a question that has remained unsolved despite substantial previous research. We sequenced and analyzed ten low-copy nuclear genes from these and related species, phasing all alleles. To test the efficacy of allele phasing on the ability to recover the evolutionary origin of polyploids, we compared these results to analyses using unphased sequences.

RESULTS

In eight of the gene trees the alleles inferred from the tetraploids formed two clades, in a non-sister relationship. Each of these clades was more closely related to alleles sampled from other species of Medicago, a pattern typical of allopolyploids. However, we also observed that alleles from one of the remaining genes formed two clades that were sister to one another, as is expected for autopolyploids. Trees inferred from unphased sequences were very different, with the tetraploids often placed in poorly supported and different positions compared to results obtained using phased alleles.

CONCLUSIONS

The complex phylogenetic history of M. arborea and M. strasseri is explained predominantly by shared allotetraploidy. We also observed that an increase in woodiness is correlated with polyploidy in this group of species and present a new possibility that woodiness could be a transgressive phenotype. Correctly phased homoeologues are likely to be critical for inferring the hybrid origin of allopolyploid species, when most genes retain more than one homoeologue. Ignoring homoeologous variation by merging the homoeologues can obscure the signal of hybrid polyploid origins and produce inaccurate results.

Development of single-copy nuclear intron markers for species-level phylogenetics: Case study with Paullinieae (Sapindaceae)

PREMISE OF THE STUDY

We developed a bioinformatic pipeline that leverages a publicly available genome and published transcriptomes to design primers in conserved coding sequences flanking targeted introns of single-copy nuclear loci. Paullinieae (Sapindaceae) is used to demonstrate the pipeline.

METHODS AND RESULTS

Transcriptome reads phylogenetically closer to the lineage of interest are aligned to the closest genome. Single-nucleotide polymorphisms are called, generating a "pseudoreference" closer to the lineage of interest. Several filters are applied to meet the criteria of single-copy nuclear loci with introns of a desired size. Primers are designed in conserved coding sequences flanking introns. Using this pipeline, we developed nine single-copy nuclear intron markers for Paullinieae.

CONCLUSIONS

This pipeline is highly flexible and can be used for any group with available genomic and transcriptomic resources. This pipeline led to the development of nine variable markers for phylogenetic study without generating sequence data de novo.

Change of Fate and Staminodial Laminarity as Potential Agents of Floral Diversification in the Zingiberales

The evolution of floral morphology in the monocot order Zingiberales shows a trend in which androecial whorl organs are progressively modified into variously conspicuous "petaloid" structures with differing degrees of fertility. Petaloidy of androecial members results from extensive laminarization of an otherwise radially symmetric structure. The genetic basis of the laminarization of androecial members has been addressed through recent candidate gene studies focused on understanding the spatiotemporal expression patterns of genes known to be necessary to floral organ formation. Here, we explore the correlation between gene duplication events and floral and inflorescence morphological diversification across the Zingiberales by inferring ancestral character states and gene copy number using the most widely accepted phylogenetic hypotheses. Our results suggest that the duplication and differential loss of GLOBOSA (GLO) copies is correlated with a change in the degree of the laminarization of androecial members. We also find an association with increased diversification in most families. We hypothesize that retention of paralogs in flower development genes could have led to a developmental shift affecting androecial organs with potential adaptive consequences, thus favoring diversification in some lineages but not others.

The phylogeny of Heliconia (Heliconiaceae) and the evolution of floral presentation

Heliconia (Heliconiaceae, order Zingiberales) is among the showiest plants of the Neotropical rainforest and represent a spectacular co-evolutionary radiation with hummingbirds. Despite the attractiveness and ecological importance of many Heliconia, the genus has been the subject of limited molecular phylogenetic studies. We sample seven markers from the plastid and nuclear genomes for 202 samples of Heliconia. This represents ca. 75% of accepted species and includes coverage of all taxonomic subgenera and sections. We date this phylogeny using fossils associated with other families in the Zingiberales; in particular we review and evaluate the Eocene fossil Ensete oregonense. We use this dated phylogenetic framework to evaluate the evolution of two components of flower orientation that are hypothesized to be important for modulating pollinator discrimination and pollen placement: resupination and erect versus pendant inflorescence habit. Our phylogenetic results suggest that the monophyletic Melanesian subgenus Heliconiopsis and a small clade of Ecuadorian species are together the sister group to the rest of Heliconia. Extant diversity of Heliconia originated in the Late Eocene (39Ma) with rapid diversification through the Early Miocene, making it the oldest known clade of hummingbird-pollinated plants. Most described subgenera and sections are not monophyletic, though closely related groups of species, often defined by shared geography, mirror earlier morphological cladistic analyses. Evaluation of changes in resupination and inflorescence habit suggests that these characters are more homoplasious than expected, and this largely explains the non-monophyly of previously circumscribed subgenera, which were based on these characters. We also find strong evidence for the correlated evolution of resupination and inflorescence habit. The correlated model suggests that the most recent common ancestor of all extant Heliconia had resupinate flowers and erect inflorescences. Finally, we note our nearly complete species sampling and dated phylogeny allow for an assessment of taxonomic history in terms of phylogenetic diversity. We find approximately half of the currently recognized species, corresponding to half of the phylogenetic diversity, have been described since 1975, highlighting the continued importance of basic taxonomic research and conservation initiatives to preserve both described and undiscovered species of Heliconia.

Species diversity driven by morphological and ecological disparity: a case study of comparative seed morphology and anatomy across a large monocot order

Phenotypic variation can be attributed to genetic heritability as well as biotic and abiotic factors. Across Zingiberales, there is a high variation in the number of species per clade and in phenotypic diversity. Factors contributing to this phenotypic variation have never been studied in a phylogenetic or ecological context. Seeds of 166 species from all eight families in Zingiberales were analyzed for 51 characters using synchrotron based 3D X-ray tomographic microscopy to determine phylogenetically informative characters and to understand the distribution of morphological disparity within the order. All families are distinguishable based on seed characters. Non-metric multidimensional scaling analyses show Zingiberaceae occupy the largest seed morphospace relative to the other families, and environmental analyses demonstrate that Zingiberaceae inhabit both temperate and tropical regions, while other Zingiberales are almost exclusively tropical. Temperate species do not cluster in morphospace nor do they share a common suite of character states. This suggests that the diversity seen is not driven by adaptation to temperate niches; rather, the morphological disparity seen likely reflects an underlying genetic plasticity that allowed Zingiberaceae to repeatedly colonize temperate environments. The notable morphoanatomical variety in Zingiberaceae seeds may account for their extraordinary ecological success and high species diversity as compared to other Zingiberales.