Using nuclear gene data for plant phylogenetics: Progress and prospects

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.

Phylogenetic analysis of endogenous viral elements in the rice genome reveals local chromosomal evolution in Oryza AA-genome species

Introduction

Rice genomes contain endogenous viral elements homologous to rice tungro bacilliform virus (RTBV) from the pararetrovirus family Caulimoviridae. These viral elements, known as endogenous RTBV-like sequences (eRTBVLs), comprise five subfamilies, eRTBVL-A, -B, -C, -D, and -X. Four subfamilies (A, B, C, and X) are present to a limited degree in the genomes of the Asian cultivated rice Oryza sativa (spp. japonica and indica) and the closely related wild species Oryza rufipogon.

Methods

The eRTBVL-D sequences are widely distributed within these and other Oryza AA-genome species. Fifteen eRTBVL-D segments identified in the japonica (Nipponbare) genome occur mostly at orthologous chromosomal positions in other AA-genome species. The eRTBVL-D sequences were inserted into the genomes just before speciation of the AA-genome species.

Results and discussion

Ten eRTBVL-D segments are located at six loci, which were used for our evolutionary analyses during the speciation of the AA-genome species. The degree of genetic differentiation varied among the eRTBVL-D segments. Of the six loci, three showed phylogenetic trees consistent with the standard speciation pattern (SSP) of the AA-genome species (Type A), and the other three represented phylogenies different from the SSP (Type B). The atypical phylogenetic trees for the Type B loci revealed chromosome region-specific evolution among the AA-genome species that is associated with phylogenetic incongruences: complex genome rearrangements between eRTBVL-D segments, an introgression between the distant species, and low genetic diversity of a shared eRTBVL-D segment. Using eRTBVL-D as an indicator, this study revealed the phylogenetic incongruence of local chromosomal regions with different topologies that developed during speciation.

Transcriptome-Based Study on the Phylogeny and Hybridization of Marattialean Ferns (Marattiaceae)

Marattiaceae is a phylogenetically isolated family of tropical eusporangiate ferns including six genera with more than one-hundred species. In Marattiaceae, monophyly of genera has been well-supported phylogenetically. However, the phylogenetic relationships among them were elusive and controversial. Here, a dataset of 26 transcriptomes (including 11 newly generated) were used to assess single-copy nuclear genes and to obtain the organelle gene sequences. Through phylotranscriptomic analysis, the phylogeny and hybridization events of Marattiaceae were explored and a robust phylogenomic framework for the evolution of Marattiaceae was provided. Using both concatenation- and coalescent-based phylogenies, the gene-tree discordance, incomplete lineage sorting (ILS) simulations, and network inference were examined. Except the low support with mitochondrial genes of Marattiaceae, nuclear genes and chloroplast genes strongly supported a sister relationship between Marattiaceae and leptosporangiate ferns. At the genus level, all phylogenetic analysis based on nuclear genes datasets recovered five genera in Marattiaceae as monophyletic with strong support. Danaea and Ptisana were the first two diverged clades in turn. Christensenia was a sister clade to the clade Marattia + Angiopteris s.l. In Angiopteris s.l., three clades (Angiopteris s.s., the Archangiopteris group, and An. sparsisora) were well identified with maximum support. The Archangiopteris group was derived from Angiopteris s.s. at ca. 18 Ma. The putative hybrid species An. sparsisora between Angiopteris s.s. and the Archangiopteris group was verified by the species network analyses and the maternal plastid genes. This study will improve our understanding for using the phylotranscriptomic method to explore phylogeny and investigate hybridization events for difficult taxa in ferns.

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.

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

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

Development and phylogenetic utilities of a new set of single-/low-copy nuclear genes in Senecioneae (Asteraceae), with new insights into the tribal position and the relationships within subtribe Tussilagininae

The tribe Senecioneae is one of the largest tribes in Asteraceae, with a nearly cosmopolitan distribution. Despite great efforts devoted to elucidate the evolution of Senecioneae, many questions still remain concerning the systematics of this group, from the tribal circumscription and position to species relationships in many genera. The hybridization-based target enrichment method of next-generation sequencing has been accepted as a promising approach to resolve phylogenetic problems. We herein develop a set of single-/low-copy genes for Senecioneae, and test their phylogenetic utilities. Our results demonstrate that these genes work highly efficiently for Senecioneae, with a high average gene recovery of 98.8% across the tribe and recovering robust phylogenetic hypotheses at different levels. In particular, the delimitation of the Senecioneae has been confirmed to include Abrotanella and exclude Doronicum, with the former sister to core Senecioneae and the latter shown to be more closely related to Calenduleae. Moreover, Doronicum and Calenduleae are inferred to be the closest relatives of Senecioneae, which is a new hypothesis well supported by statistical topology tests, morphological evidence, and the profile of pyrrolizidine alkaloids, a special kind of chemical characters generally used to define Senecioneae. Furthermore, this study suggests a complex reticulation history in the diversification of Senecioneae, accounting for the prevalence of polyploid groups in the tribe. With subtribe Tussilagininae s.str. as a case study showing a more evident pattern of gene duplication, we further explored reconstructing the phylogeny in the groups with high ploidy levels. Our results also demonstrate that tree topologies based on sorted paralogous copies are stable across different methods of phylogenetic inference, and more congruent with the morphological evidence and the results of previous phylogenetic studies.

Resolving Recalcitrant Clades in the Pantropical Ochnaceae: Insights From Comparative Phylogenomics of Plastome and Nuclear Genomic Data Derived From Targeted Sequencing

Plastid DNA sequence data have been traditionally widely used in plant phylogenetics because of the high copy number of plastids, their uniparental inheritance, and the blend of coding and non-coding regions with divergent substitution rates that allow the reconstruction of phylogenetic relationships at different taxonomic ranks. In the present study, we evaluate the utility of the plastome for the reconstruction of phylogenetic relationships in the pantropical plant family Ochnaceae (Malpighiales). We used the off-target sequence read fraction of a targeted sequencing study (targeting nuclear loci only) to recover more than 100 kb of the plastid genome from the majority of the more than 200 species of Ochnaceae and all but two genera using de novo and reference-based assembly strategies. Most of the recalcitrant nodes in the family's backbone were resolved by our plastome-based phylogenetic inference, corroborating the most recent classification system of Ochnaceae and findings from a phylogenomic study based on nuclear loci. Nonetheless, the phylogenetic relationships within the major clades of tribe Ochnineae, which comprise about two thirds of the family's species diversity, received mostly low support. Generally, the phylogenetic resolution was lowest at the infrageneric level. Overall there was little phylogenetic conflict compared to a recent analysis of nuclear loci. Effects of taxon sampling were invoked as the most likely reason for some of the few well-supported discords. Our study demonstrates the utility of the off-target fraction of a target enrichment study for assembling near-complete plastid genomes for a large proportion of samples.

Large-scale phylogenomics reveals ancient introgression in Asian Hepatica and new insights into the origin of the insular endemic Hepatica maxima

Hepatica maxima is native to Ulleungdo, which is one of the oceanic islands in Korea, and it likely originated via anagenetic speciation from the Korean mainland species H. asiatica. However, the relationships among the Asian lineages remain unresolved. Phylogenomics based on plant genomes can provide new insights into the evolutionary history of plants. We first generated plastid, mitochondrial and transcriptome sequences of the insular endemic species H. maxima. Using the genomic data for H. maxima, we obtained a phylogenomic dataset consisting of 76 plastid, 37 mitochondrial and 413 nuclear genes from Asian Hepatica and two outgroups. Coalescent- and concatenation-based methods revealed cytonuclear and organellar discordance in the lineage. The presence of gynodioecy with cytoplasmic male sterility in Asian Hepatica suggests that the discordance is correlated with potential disruption of linkage disequilibrium between the organellar genomes. Species network analyses revealed a deep history of hybridization and introgression in Asian Hepatica. We discovered that ancient and recent introgression events occurred throughout the evolutionary history of the insular endemic species H. maxima. The introgression may serve as an important source of genetic variation to facilitate adaptation to the Ulleungdo environment.

Phylogenomic proof of Recurrent Demipolyploidization and Evolutionary Stalling of the "Triploid Bridge" in Arundo (Poaceae)

Polyploidization is a frequent phenomenon in plants, which entails the increase from one generation to the next by multiples of the haploid number of chromosomes. While tetraploidization is arguably the most common and stable outcome of polyploidization, over evolutionary time triploids often constitute only a transient phase, or a "triploid bridge", between diploid and tetraploid levels. In this study, we reconstructed in a robust phylogenomic and statistical framework the evolutionary history of polyploidization in Arundo, a small genus from the Poaceae family with promising biomass, bioenergy and phytoremediation species. Through the obtainment of 10 novel leaf transcriptomes for Arundo and outgroup species, our results prove that recurrent demiduplication has likely been a major driver of evolution in this species-poor genus. Molecular dating further demonstrates that the species originating by demiduplication stalled in the "triploid bridge" for evolutionary times in the order of millions of years without undergoing tetratploidization. Nevertheless, we found signatures of molecular evolution highlighting some of the processes that accompanied the genus radiation. Our results clarify the complex nature of Arundo evolution and are valuable for future gene functional validation as well as reverse and comparative genomics efforts in the Arundo genus and other Arundinoideae.

A Tale of Two Families: Whole Genome and Segmental Duplications Underlie Glutamine Synthetase and Phosphoenolpyruvate Carboxylase Diversity in Narrow-Leafed Lupin (Lupinus angustifolius L.)

Narrow-leafed lupin (Lupinus angustifolius L.) has recently been supplied with advanced genomic resources and, as such, has become a well-known model for molecular evolutionary studies within the legume family—a group of plants able to fix nitrogen from the atmosphere. The phylogenetic position of lupins in Papilionoideae and their evolutionary distance to other higher plants facilitates the use of this model species to improve our knowledge on genes involved in nitrogen assimilation and primary metabolism, providing novel contributions to our understanding of the evolutionary history of legumes. In this study, we present a complex characterization of two narrow-leafed lupin gene families—glutamine synthetase (GS) and phosphoenolpyruvate carboxylase (PEPC). We combine a comparative analysis of gene structures and a synteny-based approach with phylogenetic reconstruction and reconciliation of the gene family and species history in order to examine events underlying the extant diversity of both families. Employing the available evidence, we show the impact of duplications on the initial complement of the analyzed gene families within the genistoid clade and posit that the function of duplicates has been largely retained. In terms of a broader perspective, our results concerning GS and PEPC gene families corroborate earlier findings pointing to key whole genome duplication/triplication event(s) affecting the genistoid lineage.

Nuclear and Chloroplast Sequences Resolve the Enigmatic Origin of the Concord Grape

Despite the commercial importance of the Concord grape, its origin has remained unresolved for over 150 years without a comprehensive phylogenetic analysis. In this study we aimed to reconstruct the evolutionary history of the Concord grape using sequence data from four nuclear markers (AT103, GAI1, PHYA, and SQD1), six plastid markers (matK, psbA-trnH, petN-trnC, ycf1, trnL-F, and trnS-G), and the plastid genome. We sampled extensively the Vitis species native to northeastern North America as well as representative species from Europe and Asia, including the commercially important Vitis vinifera (wine grape), a native European species with hermaphroditic flowers, and its wild progenitor, V. vinifera subsp. sylvestris. We also sequenced the plastid genome of one accession of the Concord grape and compared the plastid genome data to the recently published data set of Vitis plastomes. Phylogenetic analyses of the plastid and nuclear data using maximum likelihood and Bayesian inference support the hybrid origin of the Concord grape. The results clearly pinpoint the wine grape, V. vinifera, as the maternal donor and the fox grape, Vitis labrusca, which is common in northeastern North America, as the paternal donor. Moreover, we infer that the breeding history of the Concord grape must have involved the backcrossing of the F1 hybrid with the paternal parent V. labrusca. This backcrossing also explains the higher morphological similarity of the Concord grape to V. labrusca than to V. vinifera. This study provides concrete genetic evidence for the hybrid origin of a widespread Vitis cultivar and is, therefore, promising for similar future studies focused on resolving ambiguous origins of major crops or to create successful hybrid fruit crops.

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.

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

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

Phylogenetic Relationships in Orobanchaceae Inferred From Low-Copy Nuclear Genes: Consolidation of Major Clades and Identification of a Novel Position of the Non-photosynthetic Orobanche Clade Sister to All Other Parasitic Orobanchaceae

Molecular phylogenetic analyses have greatly advanced our understanding of phylogenetic relationships in Orobanchaceae, a model system to study parasitism in angiosperms. As members of this group may lack some genes widely used for phylogenetic analysis and exhibit varying degrees of accelerated base substitution in other genes, relationships among major clades identified previously remain contentious. To improve inferences of phylogenetic relationships in Orobanchaceae, we used two pentatricopeptide repeat (PPR) and three low-copy nuclear (LCN) genes, two of which have been developed for this study. Resolving power and level of support strongly differed among markers. Despite considerable incongruence among newly and previously sequenced markers, monophyly of major clades identified in previous studies was confirmed and, especially in analyses of concatenated data, strongly supported after the exclusion of a small group of East Asian genera (Pterygiella and Phtheirospermum) from the Euphrasia-Rhinanthus clade. The position of the Orobanche clade sister to all other parasitic Orobanchaceae may indicate that the shift to holoparasitism occurred early in the evolution of the family. Although well supported in analyses of concatenated data comprising ten loci (five newly and five previously sequenced), relationships among major clades, most prominently the Striga-Alectra clade, the Euphrasia-Rhinanthus clade, and the Castilleja-Pedicularis clade, were uncertain because of strongly supported incongruence also among well-resolving loci. Despite the limitations of using a few selected loci, congruence among markers with respect to circumscription of major clades of Orobanchaceae renders those frameworks for detailed, species-level, phylogenetic studies.

Hybridization and polyploidization within the Chenopodium album aggregate analysed by means of cytological and molecular markers

Hybridization and polyploidization represent an important speciation mechanism in the diploid-polyploid complex of the Chenopodium album aggregate. In the present study we successfully reconstructed the evolutionary histories of the majority of Eurasian representatives of the C. album aggregate, resulting in the most comprehensive phylogenetic analysis of this taxonomically intricate group of species to date. We applied a combination of classical karyology for precise chromosome number determination, genomic in-situ hybridization for the determination of genomic composition, flow cytometry for the estimation of genome size and sequencing of plastid (cpDNA) and nuclear (ribosomal internal transcribed spacer - ITS and the introns of the FLOWERING LOCUS T LIKE genes - FTL) markers for a phylogenetic reconstruction and the identification of parental genomes in polyploid taxa. The FTL markers identified eight well supported evolutionary lineages. Five of them include at least one diploid species, and the remaining three comprise solely the subgenomes of polyploids that probably represent extinct or unknown diploid taxa. The existence of eight basic diploid lineages explains the origin of seven Eurasian polyploid groups and brings evidence of a nearly unlimited number of subgenomic combinations. The supposed promiscuity generated new species wherever different diploid lineages met each other and gave rise to tetraploid species or whenever they met other tetraploid species to produce hexaploid species throughout their evolutionary history. Finally, we unravelled a surprisingly simple scheme of polyploid species formation within the C. album aggregate. We determined seven groups of polyploid species differing in their origin in either Eurasia or Africa and convincingly demonstrated that (1) all Chenopodium polyploid species under study are of allopolyploid origin, (2) there are eight major monophyletic evolutionary lineages represented by extant or extinct/unknown diploid taxa, (3) those monophyletic lineages represent individual subgenomes, (4) hybridization among the lineages created seven subgenomic combinations of polyploid taxa, (5) taxa represented by particular subgenome combinations were further subjected to diversification, and (6) the majority of species are relatively young, not exceeding the age of the Quaternary period.

Comparison of taxon-specific versus general locus sets for targeted sequence capture in plant phylogenomics

PREMISE OF THE STUDY

Targeted sequence capture can be used to efficiently gather sequence data for large numbers of loci, such as single-copy nuclear loci. Most published studies in plants have used taxon-specific locus sets developed individually for a clade using multiple genomic and transcriptomic resources. General locus sets can also be developed from loci that have been identified as single-copy and have orthologs in large clades of plants.

METHODS

We identify and compare a taxon-specific locus set and three general locus sets (conserved ortholog set [COSII], shared single-copy nuclear [APVO SSC] genes, and pentatricopeptide repeat [PPR] genes) for targeted sequence capture in Buddleja (Scrophulariaceae) and outgroups. We evaluate their performance in terms of assembly success, sequence variability, and resolution and support of inferred phylogenetic trees.

RESULTS

The taxon-specific locus set had the most target loci. Assembly success was high for all locus sets in Buddleja samples. For outgroups, general locus sets had greater assembly success. Taxon-specific and PPR loci had the highest average variability. The taxon-specific data set produced the best-supported tree, but all data sets showed improved resolution over previous non-sequence capture data sets.

DISCUSSION

General locus sets can be a useful source of sequence capture targets, especially if multiple genomic resources are not available for a taxon.

Robust Phylogeny of Tetrastigma (Vitaceae) Based on Ten Plastid DNA Regions: Implications for Infrageneric Classification and Seed Character Evolution

Tetrastigma (Miq.) Planch. is one of the most species-rich genera of the economically and agronomically important grape family Vitaceae. It includes ca. 95 species widely distributed in the tropics and subtropics of Asia and Australia. Species of Tetrastigma exhibit great diversity in both vegetative and reproductive characters. Here we inferred a well-supported phylogeny of Tetrastigma based on ten chloroplast DNA regions with an expanded taxon sampling of 72 species and two varieties. Our molecular results support six major clades within Tetrastigma and the relationships among these clades were well-resolved. We also documented seed morphology of 44 species covering the six major clades of the genus. Ancestral states of eight characters (seed shape, seed surface rumination pattern, chalaza length/width ratio, chalaza position, ventral infold position, ventral infold divergence, ventral infold depth in cross section, and endosperm shape) were reconstructed in Mesquite and R with four models. Character optimizations suggest that all character states have evolved multiple times except that the irregular-shaped surface rumination has derived only once in Tetrastigma. We evaluated the taxonomic importance of seed morphology and identified potential morphological evidence to support each major clade. Our comprehensive analyses of Tetrastigma shed insights into the infrageneric classification of this morphologically diverse and ecologically important genus in tropical and subtropical Asia.

Hedysarum L. (Fabaceae: Hedysareae) Is Not Monophyletic - Evidence from Phylogenetic Analyses Based on Five Nuclear and Five Plastid Sequences

The legume family (Fabaceae) exhibits a high level of species diversity and evolutionary success worldwide. Previous phylogenetic studies of the genus Hedysarum L. (Fabaceae: Hedysareae) showed that the nuclear and the plastid topologies might be incongruent, and the systematic position of the Hedysarum sect. Stracheya clade was uncertain. In this study, phylogenetic relationships of Hedysarum were investigated based on the nuclear ITS, ETS, PGDH, SQD1, TRPT and the plastid psbA-trnH, trnC-petN, trnL-trnF, trnS-trnG, petN-psbM sequences. Both nuclear and plastid data support two major lineages in Hedysarum: the Hedysarum s.s. clade and the Sartoria clade. In the nuclear tree, Hedysarum is biphyletic with the Hedysarum s.s. clade sister to the Corethrodendron + Eversmannia + Greuteria + Onobrychis clade (the CEGO clade), whereas the Sartoria clade is sister to the genus Taverniera DC. In the plastid tree, Hedysarum is monophyletic and sister to Taverniera. The incongruent position of the Hedysarum s.s. clade between the nuclear and plastid trees may be best explained by a chloroplast capture hypothesis via introgression. The Hedysarum sect. Stracheya clade is resolved as sister to the H. sect. Hedysarum clade in both nuclear and plastid trees, and our analyses support merging Stracheya into Hedysarum. Based on our new evidence from multiple sequences, Hedysarum is not monophyletic, and its generic delimitation needs to be reconsidered.

Reproductive biology and variation of nuclear ribosomal ITS and ETS sequences in the Calligonum mongolicum complex (Polygonaceae)

To explore the biosystematics of the Calligonum mongolicum complex (Polygonaceae), the flowering phenological period, breeding and pollination characters and seed set of the complex (Calligonum Mongolicum Turze, Calligonum chinense A. Los., Calligonum gobicum A. Los., Calligonum pumilum A. Los. and Calligonum zaidamense A. Los.) were documented in the Turpan Eremophyte Botanical Garden, China. The sequences of the nuclear ribosomal ITS and ETS region were employed to differentiate the Calligonum mongolicum complex and other species in sect. Medusae. The results showed species of the Calligonum mongolicum complex occupied overlapping flowering periods and had consistent pollination agents. Their breeding systems are all self-compatible, tend to be out-crossing and they interbreed amongst each other (out-crossing index, OCI = 4).The crosses within and amongst species had high seed sets (44 - 65%). Phylogenetic analyses of Calligonum sect. Medusae and the network analysis of nrDNA (ITS and ETS) in the complex suggest interbreeding amongst "species" within the complex and provide evidence for taxonomically merging the five species in the complex. The detected hybridisation, occurring within the complex, suggests the need to improve traditional methods of ex situ plant conservation in botanical gardens for maintaining genetic diversity of Calligonum within and amongst species from different geographic areas.

Next-generation polyploid phylogenetics: rapid resolution of hybrid polyploid complexes using PacBio single-molecule sequencing

Difficulties in generating nuclear data for polyploids have impeded phylogenetic study of these groups. We describe a high-throughput protocol and an associated bioinformatics pipeline (Pipeline for Untangling Reticulate Complexes (Purc)) that is able to generate these data quickly and conveniently, and demonstrate its efficacy on accessions from the fern family Cystopteridaceae. We conclude with a demonstration of the downstream utility of these data by inferring a multi-labeled species tree for a subset of our accessions. We amplified four c. 1-kb-long nuclear loci and sequenced them in a parallel-tagged amplicon sequencing approach using the PacBio platform. Purc infers the final sequences from the raw reads via an iterative approach that corrects PCR and sequencing errors and removes PCR-mediated recombinant sequences (chimeras). We generated data for all gene copies (homeologs, paralogs, and segregating alleles) present in each of three sets of 50 mostly polyploid accessions, for four loci, in three PacBio runs (one run per set). From the raw sequencing reads, Purc was able to accurately infer the underlying sequences. This approach makes it easy and economical to study the phylogenetics of polyploids, and, in conjunction with recent analytical advances, facilitates investigation of broad patterns of polyploid evolution.

An expanded nuclear phylogenomic PCR toolkit for Sapindales

PREMISE OF THE STUDY

We tested PCR amplification of 91 low-copy nuclear gene loci in taxa from Sapindales using primers developed for Bursera simaruba (Burseraceae).

METHODS AND RESULTS

Cross-amplification of these markers among 10 taxa tested was related to their phylogenetic distance from B. simaruba. On average, each Sapindalean taxon yielded product for 53 gene regions (range: 16-90). Arabidopsis thaliana (Brassicales), by contrast, yielded product for two. Single representatives of Anacardiaceae and Rutacaeae yielded 34 and 26 products, respectively. Twenty-six primer pairs worked for all Burseraceae species tested if highly divergent Aucoumea klaineana is excluded, and eight of these amplified product in every Sapindalean taxon.

CONCLUSIONS

Our study demonstrates that customized primers for Bursera can amplify product in a range of Sapindalean taxa. This collection of primer pairs, therefore, is a valuable addition to the toolkit for nuclear phylogenomic analyses of Sapindales and warrants further investigation.

Multiple Events of Allopolyploidy in the Evolution of the Racemose Lineages in Prunus (Rosaceae) Based on Integrated Evidence from Nuclear and Plastid Data

Prunus is an economically important genus well-known for cherries, plums, almonds, and peaches. The genus can be divided into three major groups based on inflorescence structure and ploidy levels: (1) the diploid solitary-flower group (subg. Prunus, Amygdalus and Emplectocladus); (2) the diploid corymbose group (subg. Cerasus); and (3) the polyploid racemose group (subg. Padus, subg. Laurocerasus, and the Maddenia group). The plastid phylogeny suggests three major clades within Prunus: Prunus-Amygdalus-Emplectocladus, Cerasus, and Laurocerasus-Padus-Maddenia, while nuclear ITS trees resolve Laurocerasus-Padus-Maddenia as a paraphyletic group. In this study, we employed sequences of the nuclear loci At103, ITS and s6pdh to explore the origins and evolution of the racemose group. Two copies of the At103 gene were identified in Prunus. One copy is found in Prunus species with solitary and corymbose inflorescences as well as those with racemose inflorescences, while the second copy (II) is present only in taxa with racemose inflorescences. The copy I sequences suggest that all racemose species form a paraphyletic group composed of four clades, each of which is definable by morphology and geography. The tree from the combined At103 and ITS sequences and the tree based on the single gene s6pdh had similar general topologies to the tree based on the copy I sequences of At103, with the combined At103-ITS tree showing stronger support in most clades. The nuclear At103, ITS and s6pdh data in conjunction with the plastid data are consistent with the hypothesis that multiple independent allopolyploidy events contributed to the origins of the racemose group. A widespread species or lineage may have served as the maternal parent for multiple hybridizations involving several paternal lineages. This hypothesis of the complex evolutionary history of the racemose group in Prunus reflects a major step forward in our understanding of diversification of the genus and has important implications for the interpretation of its phylogeny, evolution, and classification.

Lineage-specific evolutionary rate in plants: Contributions of a screening for Cereus (Cactaceae)

PREMISE OF THE STUDY

Predictable chloroplast DNA (cpDNA) sequences have been listed for the shallowest taxonomic studies in plants. We investigated whether plastid regions that vary between closely allied species could be applied for intraspecific studies and compared the variation of these plastid segments with two nuclear regions.

METHODS

We screened 16 plastid and two nuclear intronic regions for species of the genus Cereus (Cactaceae) at three hierarchical levels (species from different clades, species of the same clade, and allopatric populations).

RESULTS

Ten plastid regions presented interspecific variation, and six of them showed variation at the intraspecific level. The two nuclear regions showed both inter- and intraspecific variation, and in general they showed higher levels of variability in almost all hierarchical levels than the plastid segments.

DISCUSSION

Our data suggest no correspondence between variation of plastid regions at the interspecific and intraspecific level, probably due to lineage-specific variation in cpDNA, which appears to have less effect in nuclear data. Despite the heterogeneity in evolutionary rates of cpDNA, we highlight three plastid segments that may be considered in initial screenings in plant phylogeographic studies.

Molecular marker development from transcript sequences and germplasm evaluation for cultivated peanut (Arachis hypogaea L.)

Molecular markers are important tools for genotyping in genetic studies and molecular breeding. The SSR and SNP are two commonly used marker systems developed from genomic or transcript sequences. The objectives of this study were to: (1) assemble and annotate the publicly available ESTs in Arachis and the in-house short reads, (2) develop and validate SSR and SNP markers, and (3) investigate the genetic diversity and population structure of the peanut breeding lines and the U.S. peanut mini core collection using developed SSR markers. An NCBI EST dataset with 252,951 sequences and an in-house 454 RNAseq dataset with 288,701 sequences were assembled separately after trimming. Transcript sequence comparison and phylogenetic analysis suggested that peanut is closer to cowpea and scarlet bean than to soybean, common bean and Medicago. From these two datasets, 6455 novel SSRs and 11,902 SNPs were identified. Of the discovered SSRs, 380 representing various SSR types were selected for PCR validation. The amplification rate was 89.2 %. Twenty-two (6.5 %) SSRs were polymorphic between at least one pair of four genotypes. Sanger sequencing of PCR products targeting 110 SNPs revealed 13 true SNPs between tetraploid genotypes and 193 homoeologous SNPs within genotypes. Eight out of the 22 polymorphic SSR markers were selected to evaluate the genetic diversity of Florida peanut breeding lines and the U.S. peanut mini core collection. This marker set demonstrated high discrimination power by displaying an average polymorphism information content value of 0.783, a combined probability of identity of 10(-11), and a combined power of exclusion of 0.99991. The structure analysis revealed four sub-populations among the peanut accessions and lines evaluated. The results of this study enriched the peanut genomic resources, provided over 6000 novel SSR markers and the credentials for true peanut SNP marker development, and demonstrated the power of newly developed SSR markers in genotyping peanut germplasm and breeding materials.

Trends in flower symmetry evolution revealed through phylogenetic and developmental genetic advances

A striking aspect of flowering plant (angiosperm) diversity is variation in flower symmetry. From an ancestral form of radial symmetry (polysymmetry, actinomorphy), multiple evolutionary transitions have contributed to instances of non-radial forms, including bilateral symmetry (monosymmetry, zygomorphy) and asymmetry. Advances in flowering plant molecular phylogenetic research and studies of character evolution as well as detailed flower developmental genetic studies in a few model species (e.g. Antirrhinum majus, snapdragon) have provided a foundation for deep insights into flower symmetry evolution. From phylogenetic studies, we have a better understanding of where during flowering plant diversification transitions from radial to bilateral flower symmetry (and back to radial symmetry) have occurred. From developmental studies, we know that a genetic programme largely dependent on the functional action of the CYCLOIDEA gene is necessary for differentiation along the snapdragon dorsoventral flower axis. Bringing these two lines of inquiry together has provided surprising insights into both the parallel recruitment of a CYC-dependent developmental programme during independent transitions to bilateral flower symmetry, and the modifications to this programme in transitions back to radial flower symmetry, during flowering plant evolution.

Molecular phylogeny and evolutionary history of the Eurasiatic orchid genus Himantoglossum s.l. (Orchidaceae)

BACKGROUND AND AIMS

Lizard orchids of the genus Himantoglossum include many of Eurasia's most spectacular orchids, producing substantial spikes of showy flowers. However, until recently the genus had received only limited, and entirely traditional, systematic study. The aim of the current work was to provide a more robust molecular phylogeny in order to better understand the evolutionary relationships among species of particular conservation concern.

METHODS

All putative species of Himantoglossum s.l. were sampled across its geographical range. A large subsample of the 153 populations studied contributed to an initial survey of nuclear ribosomal internal transcribed spacer (nrITS) ribotypes. Smaller subsets were then sequenced for four plastid regions and the first intron of the low-copy-number nuclear gene LEAFY. Rooted using Steveniella as outgroup, phylogenetic trees were generated using parsimony and Bayesian methods from each of the three datasets, supplemented with a ribotype network.

KEY RESULTS

The resulting trees collectively determined the order of branching of the early divergent taxa as Himantoglossum comperianum > H. robertianum group > H. formosum, events that also involved significant morphological divergence. Relaxed molecular clock dating suggested that these divergences preceded the Pleistocene glaciations (the origin of the H. robertianum group may have coincided with the Messinian salinity crisis) and occurred in Asia Minor and/or the Caucasus. Among more controversial taxa of the H. hircinum-jankae clade, which are only subtly morphologically divergent, topological resolution was poorer and topological incongruence between datasets was consequently greater.

CONCLUSIONS

Plastid sequence divergence is broadly consistent with prior, morphologically circumscribed taxa and indicates a division between H. hircinum-adriaticum to the west of the Carpathians and H. jankae-caprinum (plus local endemics) to the east, a distinction also suggested by nrITS ribotypes. LEAFY phylogenies are less congruent with prior taxonomic arrangements and include one likely example of paralogy. Himantoglossum metlesicsianum fully merits its IUCN Endangered status. Potentially significant genetic variation was detected within Steveniella satyrioides, H. robertianum and H. hircinum. However, confident circumscription of the more derived species of Himantoglossum s.s., including local endemics of hybrid origin, must await future morphometric and population-genetic analyses.

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.

Development of novel, exon-primed intron-crossing (EPIC) markers from EST databases and evaluation of their phylogenetic utility in Commiphora (Burseraceae)

PREMISE OF THE STUDY

Novel nuclear exon-primed intron-crossing (EPIC) markers were developed to increase phylogenetic resolution among recently diverged lineages in the frankincense and myrrh family, Burseraceae, using Citrus, Arabidopsis, and Oryza genome resources. •

METHODS AND RESULTS

Primer pairs for 48 nuclear introns were developed using the genome resource IntrEST and were screened using species of Commiphora and other Burseraceae taxa. Four putative intron regions (RPT6A, BXL2, mtATP Synthase D, and Rab6) sequenced successfully for multiple taxa and recovered phylogenies consistent with those of existing studies. In some cases, these regions yielded informative sequence variation on par with that of the nuclear ribosomal DNA internal transcribed spacer. •

CONCLUSIONS

The combination of freely available genome resources and our design criteria have uncovered four single-copy nuclear intron regions that are useful for phylogenetic reconstruction of Burseraceae taxa. Because our EPIC primers also amplify Arabidopsis, we recommend their trial in other rosid and eudicot lineages.

Molecular phylogenetics of New Caledonian Diospyros (Ebenaceae) using plastid and nuclear markers

To clarify phylogenetic relationships among New Caledonian species of Diospyros, sequences of four plastid markers (atpB, rbcL, trnK-matK and trnS-trnG) and two low-copy nuclear markers (ncpGS and PHYA) were analysed. New Caledonian Diospyros species fall into three clades, two of which have only a few members (1 or 5 species); the third has 21 closely related species for which relationships among species have been mostly unresolved in a previous study. Although species of the third group (NC clade III) are morphologically distinct and largely occupy different habitats, they exhibit little molecular variability. Diospyros vieillardii is sister to the rest of the NC clade III, followed by D. umbrosa and D. flavocarpa, which are sister to the rest of this clade. Species from coastal habitats of western Grande Terre (D. cherrieri and D. veillonii) and some found on coralline substrates (D. calciphila and D. inexplorata) form two well-supported subgroups. The species of NC clade III have significantly larger genomes than found in diploid species of Diospyros from other parts of the world, but they all appear to be diploids. By applying a molecular clock, we infer that the ancestor of the NC clade III arrived in New Caledonia around 9 million years ago. The oldest species are around 7 million years old and the youngest ones probably much less than 1 million years.

Phylogeny of ten species of the genus Hordeum L. as revealed by AFLP markers and seed storage protein electrophoresis

The phylogenetic relationships of 60 accessions representing ten species of the genus Hordeum were investigated based on AFLP markers and seed storage protein SDS-PAGE electrophoresis. A total of 339 AFLP polymorphic markers were scored as a result of fingerprinting the studied taxa using seven AFLP primer combinations, whereas 46 polymorphic protein bands resulted from the water soluble and water non-soluble seed storage protein electrophoresis. The phylogenetic tree deduced from AFLP analysis is concordant in a large extent with that deduced from seed storage protein electrophoresis. The studied taxa were clustered according to their genome type into two main groups representing the Old and New World's species. Inside each group the species were clustered according to their genome type. Highly significant cophenetic correlation coefficient was obtained between both AFLP (0.96) and seed storage protein (0.89) indicating the reliability of the results.

Transcriptome-mining for single-copy nuclear markers in ferns

BACKGROUND

Molecular phylogenetic investigations have revolutionized our understanding of the evolutionary history of ferns-the second-most species-rich major group of vascular plants, and the sister clade to seed plants. The general absence of genomic resources available for this important group of plants, however, has resulted in the strong dependence of these studies on plastid data; nuclear or mitochondrial data have been rarely used. In this study, we utilize transcriptome data to design primers for nuclear markers for use in studies of fern evolutionary biology, and demonstrate the utility of these markers across the largest order of ferns, the Polypodiales.

PRINCIPAL FINDINGS

We present 20 novel single-copy nuclear regions, across 10 distinct protein-coding genes: ApPEFP_C, cryptochrome 2, cryptochrome 4, DET1, gapCpSh, IBR3, pgiC, SQD1, TPLATE, and transducin. These loci, individually and in combination, show strong resolving power across the Polypodiales phylogeny, and are readily amplified and sequenced from our genomic DNA test set (from 15 diploid Polypodiales species). For each region, we also present transcriptome alignments of the focal locus and related paralogs-curated broadly across ferns-that will allow researchers to develop their own primer sets for fern taxa outside of the Polypodiales. Analyses of sequence data generated from our genomic DNA test set reveal strong effects of partitioning schemes on support levels and, to a much lesser extent, on topology. A model partitioned by codon position is strongly favored, and analyses of the combined data yield a Polypodiales phylogeny that is well-supported and consistent with earlier studies of this group.

CONCLUSIONS

The 20 single-copy regions presented here more than triple the single-copy nuclear regions available for use in ferns. They provide a much-needed opportunity to assess plastid-derived hypotheses of relationships within the ferns, and increase our capacity to explore aspects of fern evolution previously unavailable to scientific investigation.

Transcriptome sequences resolve deep relationships of the grape family

Previous phylogenetic studies of the grape family (Vitaceae) yielded poorly resolved deep relationships, thus impeding our understanding of the evolution of the family. Next-generation sequencing now offers access to protein coding sequences very easily, quickly and cost-effectively. To improve upon earlier work, we extracted 417 orthologous single-copy nuclear genes from the transcriptomes of 15 species of the Vitaceae, covering its phylogenetic diversity. The resulting transcriptome phylogeny provides robust support for the deep relationships, showing the phylogenetic utility of transcriptome data for plants over a time scale at least since the mid-Cretaceous. The pros and cons of transcriptome data for phylogenetic inference in plants are also evaluated.

Developing nuclear DNA phylogenetic markers in the angiosperm genus Leucadendron (Proteaceae): a next-generation sequencing transcriptomic approach

Despite the recent advances in generating molecular data, reconstructing species-level phylogenies for non-models groups remains a challenge. The use of a number of independent genes is required to resolve phylogenetic relationships, especially for groups displaying low polymorphism. In such cases, low-copy nuclear exons and non-coding regions, such as 3' untranslated regions (3'-UTRs) or introns, constitute a potentially interesting source of nuclear DNA variation. Here, we present a methodology meant to identify new nuclear orthologous markers using both public-nucleotide databases and transcriptomic data generated for the group of interest by using next generation sequencing technology. To identify PCR primers for a non-model group, the genus Leucadendron (Proteaceae), we adopted a framework aimed at minimizing the probability of paralogy and maximizing polymorphism. We anchored when possible the right-hand primer into the 3'-UTR and the left-hand primer into the coding region. Seven new nuclear markers emerged from this search strategy, three of those included 3'-UTRs. We further compared the phylogenetic potential between our new markers and the ribosomal internal transcribed spacer region (ITS). The sequenced 3'-UTRs yielded higher polymorphism rates than the ITS region did. We did not find strong incongruences with the phylogenetic signal contained in the ITS region and the seven new designed markers but they strongly improved the phylogeny of the genus Leucadendron. Overall, this methodology is efficient in isolating orthologous loci and is valid for any non-model group given the availability of transcriptomic data.

Chromosomal diversification and karyotype evolution of diploids in the cytologically diverse genus Prospero (Hyacinthaceae)

Background

Prospero (Hyacinthaceae) provides a unique system to assess the impact of genome rearrangements on plant diversification and evolution. The genus exhibits remarkable chromosomal variation but very little morphological differentiation. Basic numbers of x = 4, 5, 6 and 7, extensive polyploidy, and numerous polymorphic chromosome variants were described, but only three species are commonly recognized: P. obtusifolium, P. hanburyi, and P. autumnale s.l., the latter comprising four diploid cytotypes. The relationship between evolutionary patterns and chromosomal variation in diploids, the basic modules of the extensive cytological diversity, is presented.

Results

Evolutionary inferences were derived from fluorescence in situ hybridization (FISH) with 5S and 35S rDNA, genome size estimations, and phylogenetic analyses of internal transcribed spacer (ITS) of 35S rDNA of 49 diploids in the three species and all cytotypes of P. autumnale s.l. All species and cytotypes possess a single 35S rDNA locus, interstitial except in P. hanburyi where it is sub-terminal, and one or two 5S rDNA loci (occasionally a third in P. obtusifolium) at fixed locations. The localization of the two rDNA types is unique for each species and cytotype. Phylogenetic data in the P. autumnale complex enable tracing of the evolution of rDNA loci, genome size, and direction of chromosomal fusions: mixed descending dysploidy of x = 7 to x = 6 and independently to x = 5, rather than successive descending dysploidy, is proposed.

Conclusions

All diploid cytotypes are recovered as well-defined evolutionary lineages. The cytogenetic and phylogenetic approaches have provided excellent phylogenetic markers to infer the direction of chromosomal change in Prospero. Evolution in Prospero, especially in the P. autumnale complex, has been driven by differentiation of an ancestral karyotype largely unaccompanied by morphological change. These new results provide a framework for detailed analyses of various types of chromosomal rearrangements and karyotypic variation in polyploids.

Evolutionary consequences, constraints and potential of polyploidy in plants

Polyploidy, the possession of more than 2 complete genomes, is a major force in plant evolution known to affect the genetic and genomic constitution and the phenotype of an organism, which will have consequences for its ecology and geography as well as for lineage diversification and speciation. In this review, we discuss phylogenetic patterns in the incidence of polyploidy including possible underlying causes, the role of polyploidy for diversification, the effects of polyploidy on geographical and ecological patterns, and putative underlying mechanisms as well as chromosome evolution and evolution of repetitive DNA following polyploidization. Spurred by technological advances, a lot has been learned about these aspects both in model and increasingly also in nonmodel species. Despite this enormous progress, long-standing questions about polyploidy still cannot be unambiguously answered, due to frequently idiosyncratic outcomes and insufficient integration of different organizational levels (from genes to ecology), but likely this will change in the near future. See also the sister article focusing on animals by Choleva and Janko in this themed issue.

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.

Gene discordance in phylogenomics of recent plant radiations, an example from Hawaiian Cyrtandra (Gesneriaceae)

Resolving species relationships within recent radiations requires analysis at the interface of phylogenetics and population genetics, where coalescence and hybridization may confound our understanding of relationships. We developed 18 new primer pairs for nuclear loci in Cyrtandra (Gesneriaceae), one of the largest plant radiations in the Pacific Islands, and tested the concordance of 14 loci in establishing the phylogenetic relationships of a small number of Hawaiian species. Four genes yielded tree topologies conflicting with the primary concordance tree, suggesting plastid capture and horizontal transfer via hybridization. Combining all concordant genes yielded a tree with stronger support and a different topology from the total-evidence tree. We conclude that a small number of genes may be insufficient for accurate reconstruction of the phylogenetic relationships among closely related species. Further, the combination of genes for phylogenetic analysis without preliminary concordance tests can yield an erroneous tree topology. It seems that the number of genes needed for phylogenetic analysis of closely related species is significantly greater than the small numbers commonly used, which fail to isolate coalescence, introgression and hybridization.

Phylogenetic relationships between Mediterranean and middle-Asian wild species of the genus Hordeum L. as revealed by biochemical and molecular markers

The phylogenetic relationships of 60 accessions of the genus Hordeum (29 Mediterranean and 20 middle-Asian wild accessions, together with nine American accessions and two of unknown origin), representing together nine species, were investigated by AFLP markers. Three hundred sixty six AFLP fragments were used for studying the molecular genetic diversity among the studied species, 339 out of them were polymorphic. Forty seven protein bands were obtained from the water soluble and the water insoluble seed storage protein by SDS-PAGE electrophoresis of 18 accessions representing nine species (two accessions/species). One band was common to all species and the other 46 bands were polymorphic. The phylogenetic tree deduced from AFLP analysis is concordant to a large extent with that deduced from seed storage protein. Highly significant cophenetic correlation coefficient was obtained between both AFLP (0.96) and seed storage protein (0.89) indicating the reliability of the results. The studied taxa were clustered according to their genome type. All Mediterranean and middle-Asian wild accessions could be integrated into the existing phylogenetic scheme.