Reprint of: using nuclear gene data for plant phylogenetics: progress and prospects

A digital PCR approach to assess the purity of oregano

Herbs and spices are food categories known to be at high risk of adulteration. Presence of undeclared foreign plant species has often been reported in oregano and may have a direct impact on its organoleptic quality and potentially the safety of this aromatic herb. A droplet digital PCR approach was developed to assess the purity of oregano by quantifying the DNA copies of oregano versus the total plant DNA copies. Nuclear single-copy genes were selected by targeting the terpene synthase 5 gene from oregano and the plant phosphoenolpyruvate carboxylase 2 gene. The reactions were specific to the Origanum genus and plant materials respectively, whereas trueness and precision data confirmed the reliability of the method to quantify oregano. The applicability of the method was further verified on proficiency test samples before being applied on commercial oregano samples.

Granule-bound starch synthase in plants: Towards an understanding of their evolution, regulatory mechanisms, applications, and perspectives

Amylose content (AC) is a significant quality trait in starchy crops, affecting their processing and application by the food and non-food industries. Therefore, fine-tuning AC in these crops has become a focus for breeders. Granule-bound starch synthase (GBSS) is the core enzyme that directly determines the AC levels. Several excellent reviews have summarized key progress in various aspects of GBSS research in recent years, but they mostly focus on cereals. Herein, we provide an in-depth review of GBSS research in monocots and dicots, focusing on the molecular characteristics, evolutionary relationships, expression patterns, molecular regulation mechanisms, and applications. We also discuss future challenges and directions for controlling AC in starchy crops, and found simultaneously increasing both the PTST and GBSS gene expression levels may be an effective strategy to increase amylose content.

Allopolyploid origin and niche expansion of Rhodiola integrifolia (Crassulaceae)

Polyploidy after hybridization between species can lead to immediate post-zygotic isolation, causing saltatory origin of new species. Although the incidence of polyploidization in plants is high, it is thought that a new polyploid lineage can succeed only if it establishes a new ecological niche divergent from its progenitor lineages. We tested the hypothesis that Rhodiola integrifolia from North America is an allopolyploid produced by R. rhodantha and R. rosea and determined whether its survival can be explained by the niche divergence hypothesis. To this end, we sequenced two low-copy nuclear genes (ncpGS and rpb2) in a phylogenetic analysis of 42 Rhodiola species and tested for niche equivalency and similarity using Schoener's D as the index of niche overlap. Our phylogeny-based approach showed that R. integrifolia possesses alleles from both R. rhodantha and R. rosea. Dating analysis showed that the hybridization event that led to R. integrifolia occurred ca. 1.67 Mya and niche modeling analysis showed that at this time, both R. rosea and R. rhodantha may have been present in Beringia, providing the opportunity for the hybridization event. We also found that the niche of R. integrifolia differs from that of its progenitors in both niche breadth and optimum. Taken together, these results confirm the hybrid origin of R. integrifolia and support the niche divergence hypothesis for this tetraploid species. Our results underscore the fact that lineages with no current overlapping distribution could produce hybrid descendants in the past, when climate oscillations made their distributions overlap.

A pipeline for assembling low copy nuclear markers from plant genome skimming data for phylogenetic use

Background

Genome skimming is a popular method in plant phylogenomics that do not include a biased enrichment step, relying on random shallow sequencing of total genomic DNA. From these data the plastome is usually readily assembled and constitutes the bulk of phylogenetic information generated in these studies. Despite a few attempts to use genome skims to recover low copy nuclear loci for direct phylogenetic use, such endeavor remains neglected. Causes might include the trade-off between libraries with few reads and species with large genomes (i.e., missing data caused by low coverage), but also might relate to the lack of pipelines for data assembling.

Methods

A pipeline and its companion R package designed to automate the recovery of low copy nuclear markers from genome skimming libraries are presented. Additionally, a series of analyses aiming to evaluate the impact of key assembling parameters, reference selection and missing data are presented.

Results

A substantial amount of putative low copy nuclear loci was assembled and proved useful to base phylogenetic inference across the libraries tested (4 to 11 times more data than previously assembled plastomes from the same libraries).

Discussion

Critical aspects of assembling low copy nuclear markers from genome skims include the minimum coverage and depth of a sequence to be used. More stringent values of these parameters reduces the amount of assembled data and increases the relative amount of missing data, which can compromise phylogenetic inference, in turn relaxing the same parameters might increase sequence error. These issues are discussed in the text, and parameter tuning through multiple comparisons tracking their effects on support and congruence is highly recommended when using this pipeline. The skimmingLoci pipeline (https://github.com/mreginato/skimmingLoci) might stimulate the use of genome skims to recover nuclear loci for direct phylogenetic use, increasing the power of genome skimming data to resolve phylogenetic relationships, while reducing the amount of sequenced DNA that is commonly wasted.

Evolutionary patterns of plastome resolve multiple origins of the Ns-containing polyploid species in Triticeae

Tracing evolutionary history proves challenging for polyploid groups that have evolved rapidly, especially if an ancestor of a polyploid is extinct. The Ns-containing polyploids are recognized as the NsXm and StHNsXm genomic constitutions in Triticeae. The Ns originated from Psathyrostachys, while the Xm represented a genome of unknown origin. Here, we use genetic information in plastome to trace the complex lineage history of the Ns-containing polyploid species by sampling 26 polyploids and 90 diploid taxa representing 23 basic genomes in Triticeae. Phylogenetic reconstruction, cluster plot of genetic distance matrix, and migration event demonstrated that (1) the Ns plastome originated from different Psathyrostachys species, and the Xm plastome may originate from an ancestral lineage of Henrardia, Agropyron, and Eremopyrum; (2) the Ns, Xm, and St genome donors separately served as the maternal parents during the speciation of the Ns-containing polyploid species, resulting in a maternal haplotype polymorphism; (3) North AmericanLeymusspecies might originate from colonization during late Miocene via the Bering land bridge and were the paternal donor of the StHNsXm genome Pascopyrum species. Our results shed new light on our understanding of the rich diversity and ecological adaptation of the Ns-containing polyploid species.

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.

Subspecies divergence and pronounced phylogenetic incongruence in the East-Asia-endemic shrub Magnolia sieboldii

BACKGROUND AND AIMS

The biogeographic patterns of the East-Asia-endemic shrub Magnolia sieboldii, in which the range of the subsp. sieboldii is interposed with the disjunct distribution of subsp. japonica, implies a complex evolutionary history, involving rapid speciation and hybridization. Here, we aim to reveal the evolutionary and phylogeographic histories of the species with a particular focus on the time of subspecies divergence, the hypothesis of secondary hybridization and the Pleistocene survival of each subspecies, using a combination of genetic analyses and ecological niche modelling.

METHODS

Genetic variation, genetic structures and phylogenetic relationships were elucidated based on nuclear low-copy genes, chloroplast DNA, and nuclear simple sequence repeats (SSRs). A scenario selection analysis and divergence time estimation were performed using coalescent simulation in DIYABC and *BEAST. Ecological niche modelling and a test of niche differentiation were performed using Maxent and ENMTools.

KEY RESULTS

All marker types showed deep, but pronouncedly incongruent, west-east genetic divergences, with the subspecies being delineated only by the nuclear low-copy genes. Phylogenetic tree topologies suggested that ancient hybridization and introgression were likely to have occurred; however, this scenario did not receive significant support in the DIYABC analysis. The subspecies differentiated their niches, but both showed a dependence on high humidity and were predicted to have persisted during the last glacial cycle by maintaining a stable latitudinal distribution via migration to lower altitudes.

CONCLUSIONS

We found a deep genetic divergence and a pronounced phylogenetic incongruence among the two subspecies of M. sieboldii, which may have been driven by major paleogeographic and paleoclimatic events that have occurred since the Neogene in East Asia, including global cooling, climate oscillations and the formation of land bridges. Both subspecies were, however, considered to persist in situ in stable climatic conditions during the late Pleistocene.

Guidelines for the Choice of Sequences for Molecular Plant Taxonomy

This chapter presents an overview of the major plant DNA sequences and molecular methods available for plant taxonomy. Guidelines are provided for the choice of sequences and methods to be used, based on the DNA compartment (nuclear, chloroplastic, mitochondrial), evolutionary mechanisms, and the level of taxonomic differentiation of the plants under survey.

Nuclear Ribosomal RNA Genes: ITS Region

Despite possible drawbacks (intraspecific polymorphisms and possible fungal contamination), sequencing of the ITS region of the ribosomal RNA genes remains one of the most popular nuclear sequences used for plant taxonomy and phylogeny. A protocol for PCR amplification and sequencing of this region using universal plant primers is provided.

Genetic Characterization of an Endangered Chilean Endemic Species, Prosopis burkartii Muñoz, Reveals its Hybrids Parentage

The hybridization of Prosopis burkartii, a critically endangered endemic species, and the identification of its paternal species has not been genetically studied before. In this study we aimed to genetically confirm the origin of this species. To resolve the parental status of P. burkartii, inter-simple sequence repeat (ISSR), simple sequence repeats (SSR) and intron trnL molecular markers were used, and compared with Chilean species from the Algarobia and Strombocarpa sections. Out of seven ISSRs, a total of 70 polymorphic bands were produced in four species of the Strombocarpa section. An Multi-dimensional scaling (MDS) and Bayasian (STRUCTURE) analysis showed signs of introgression of genetic material in P. burkartii. Unweighted pair group method with arithmetic average (UPGMA) cluster analysis showed three clusters, and placed the P. burkartii cluster nested within the P. tamarugo group. Sequencing of the trnL intron showed a fragment of 535 bp and 529 bp in the species of the Algarobia and Strombocarpa sections, respectively. Using maximum parsimony (MP) and maximum likelihood (ML) trees with the trnL intron, revealed four clusters. A species-specific diagnostic method was performed, using the trnL intron Single Nucleotide Polymorphism (SNP). This method identified if individuals of P. burkartii inherited their maternal DNA from P. tamarugo or from P. strombulifera. We deduced that P. tamarugo and P. strombulifera are involved in the formation of P. burkartii.

Allopolyploid origin in Rubus (Rosaceae) inferred from nuclear granule-bound starch synthase I (GBSSI) sequences

BACKGROUND

Polyploidy and hybridization are ubiquitous in Rubus L., a large and taxonomically challenging genus. Chinese Rubus are mainly concentrated into two major sections, the diploid Idaeobatus and the polyploid Malachobatus. However, it remains unclear to be auto- or allo- polyploid origin of polyploids in Rubus. We investigated the homoeologs and the structure of the GBSSI-1 (granule-bound starch synthase I) gene in 140 Rubus individuals representing 102 taxa in 17 (out of the total 24) subsections of 7 (total of 12) sections at different ploidy levels.

RESULTS

Based on the gene structure and sequence divergence, we defined three gene variants, GBSSI-1a, GBSSI-1b, and GBSSI-1c. When compared with GBSSI-1a, both GBSSI-1b and GBSSI-1c have a shorter fourth intron, and GBSSI-1c had an additional deletion in the fifth intron. For diploids, either GBSSI-1a or GBSSI-1b was detected in 56 taxa consisting of 82 individuals from sect. Idaeobatus, while both alleles existed in R. pentagonus and R. peltatus. Both homoeologs GBSSI-1a and GBSSI-1b were identified in 39 taxa (48 individuals) of Malachobatus polyploids. They were also observed in two sect. Dalibardastrum taxa, in one sect. Chamaebatus taxon, and in three taxa from sect. Cylactis. Interestingly, all three homoeologs were observed in the three tetraploid taxa. Phylogenetic trees and networks suggested two clades (I and II), corresponding to GBSSI-1a, and GBSSI-1b/1c sequences, respectively. GBSSI-1 homoeologs from the same polyploid individual were resolved in different well-supported clades, and some of these homoelogs were more closely related to homoelogs in other species than they were to each other. This implied that the homoeologs of these polyploids were donated by different ancestral taxa, indicating their allopolyploid origin. Two kinds of diploids hybridized to form most allotetraploid species. The early-divergent diploid species with GBSSI-1a or -1b emerged before polyploid formation in the evolutionary history of Rubus.

CONCLUSION

This study provided new insights into allopolyploid origin and evolution from diploid to polyploid within the genus Rubus at the molecular phylogenetic level, consistent with the taxonomic treatment by Yü et al. and Lu.

Potato virus Y; the Andean connection

Potato virus Y (PVY) causes disease in potatoes and other solanaceous crops. The appearance of its necrogenic strains in the 1980s made it the most economically important virus of potatoes. We report the isolation and genomic sequences of 32 Peruvian isolates of PVY which, together with 428 published PVY genomic sequences, gave an alignment of 460 sequences. Of these 190 (41%) were non-recombinant, and 162 of these provided a dated phylogeny, that corresponds well with the likely history of PVY, and show that PVY originated in South America which is where potatoes were first domesticated. The most basal divergences of the PVY population produced the N and C: O phylogroups; the origin of the N phylogroup is clearly Andean, but that of the O and C phylogroups is unknown, although they may have been first to establish in European crops. The current PVY population originated around 156 CE. PVY was probably first taken from South America to Europe in the 16th century in tubers. Most of the present PVY diversity emerged in the second half of the 19th century, after the Phytophthora infestans epidemics of the mid-19th century destroyed the European crop and stimulated potato breeding. Imported breeding lines were shared, and there was no quarantine. The early O population was joined later by N phylogroup isolates and their recombinants generated the R1 and R2 populations of damaging necrogenic strains. Our dating study has confirmed that human activity has dominated the phylodynamics of PVY for the last two millennia.

A New Perspective on Polyploid Fragaria (Strawberry) Genome Composition Based on Large-Scale, Multi-Locus Phylogenetic Analysis

The subgenomic compositions of the octoploid (2n = 8× = 56) strawberry (Fragaria) species, including the economically important cultivated species Fragaria x ananassa, have been a topic of long-standing interest. Phylogenomic approaches utilizing next-generation sequencing technologies offer a new window into species relationships and the subgenomic compositions of polyploids. We have conducted a large-scale phylogenetic analysis of Fragaria (strawberry) species using the Fluidigm Access Array system and 454 sequencing platform. About 24 single-copy or low-copy nuclear genes distributed across the genome were amplified and sequenced from 96 genomic DNA samples representing 16 Fragaria species from diploid (2×) to decaploid (10×), including the most extensive sampling of octoploid taxa yet reported. Individual gene trees were constructed by different tree-building methods. Mosaic genomic structures of diploid Fragaria species consisting of sequences at different phylogenetic positions were observed. Our findings support the presence in octoploid species of genetic signatures from at least five diploid ancestors (F. vesca, F. iinumae, F. bucharica, F. viridis, and at least one additional allele contributor of unknown identity), and questions the extent to which distinct subgenomes are preserved over evolutionary time in the allopolyploid Fragaria species. In addition, our data support divergence between the two wild octoploid species, F. virginiana and F. chiloensis.

Plant phylogenomics based on genome-partitioning strategies: Progress and prospects

The rapid expansion of next-generation sequencing (NGS) has generated a powerful array of approaches to address fundamental questions in biology. Several genome-partitioning strategies to sequence selected subsets of the genome have emerged in the fields of phylogenomics and evolutionary genomics. In this review, we summarize the applications, advantages and limitations of four NGS-based genome-partitioning approaches in plant phylogenomics: genome skimming, transcriptome sequencing (RNA-seq), restriction site associated DNA sequencing (RAD-Seq), and targeted capture (Hyb-seq). Of these four genome-partitioning approaches, targeted capture (especially Hyb-seq) shows the greatest promise for plant phylogenetics over the next few years. This review will aid researchers in their selection of appropriate genome-partitioning approaches to address questions of evolutionary scale, where we anticipate continued development and expansion of whole-genome sequencing strategies in the fields of plant phylogenomics and evolutionary biology research.

Evolution of Rosaceae Fruit Types Based on Nuclear Phylogeny in the Context of Geological Times and Genome Duplication

Fruits are the defining feature of angiosperms, likely have contributed to angiosperm successes by protecting and dispersing seeds, and provide foods to humans and other animals, with many morphological types and important ecological and agricultural implications. Rosaceae is a family with ∼3000 species and an extraordinary spectrum of distinct fruits, including fleshy peach, apple, and strawberry prized by their consumers, as well as dry achenetum and follicetum with features facilitating seed dispersal, excellent for studying fruit evolution. To address Rosaceae fruit evolution and other questions, we generated 125 new transcriptomic and genomic datasets and identified hundreds of nuclear genes to reconstruct a well-resolved Rosaceae phylogeny with highly supported monophyly of all subfamilies and tribes. Molecular clock analysis revealed an estimated age of ∼101.6 Ma for crown Rosaceae and divergence times of tribes and genera, providing a geological and climate context for fruit evolution. Phylogenomic analysis yielded strong evidence for numerous whole genome duplications (WGDs), supporting the hypothesis that the apple tribe had a WGD and revealing another one shared by fleshy fruit-bearing members of this tribe, with moderate support for WGDs in the peach tribe and other groups. Ancestral character reconstruction for fruit types supports independent origins of fleshy fruits from dry-fruit ancestors, including the evolution of drupes (e.g., peach) and pomes (e.g., apple) from follicetum, and drupetum (raspberry and blackberry) from achenetum. We propose that WGDs and environmental factors, including animals, contributed to the evolution of the many fruits in Rosaceae, which provide a foundation for understanding fruit evolution.

Contrasting evolutionary patterns of multiple loci uncover new aspects in the genome origin and evolutionary history of Leymus (Triticeae; Poaceae)

Leymus Hochst. (Triticeae: Poaceae), a group of allopolyploid species with the NsXm genomes, is a perennial genus with diversity in morphology, cytology, ecology, and distribution in the Triticeae. To investigate the genome origin and evolutionary history of Leymus, three unlinked low-copy nuclear genes (Acc1, Pgk1, and GBSSI) and three chloroplast regions (trnL-F, matK, and rbcL) of 32 Leymus species were analyzed with those of 36 diploid species representing 18 basic genomes in the Triticeae. The phylogenetic relationships were reconstructed using Bayesian inference, Maximum parsimony, and NeighborNet methods. A time-calibrated phylogeny was generated to estimate the evolutionary history of Leymus. The results suggest that reticulate evolution has occurred in Leymus species, with several distinct progenitors contributing to the Leymus. The molecular data in resolution of the Xm-genome lineage resulted in two apparently contradictory results, with one placing the Xm-genome lineage as closely related to the P/F genome and the other splitting the Xm-genome lineage as sister to the Ns-genome donor. Our results suggested that (1) the Ns genome of Leymus was donated by Psathyrostachys, and additional Ns-containing alleles may be introgressed into some Leymus polyploids by recurrent hybridization; (2) The phylogenetic incongruence regarding the resolution of the Xm-genome lineage suggested that the Xm genome of Leymus was closely related to the P genome of Agropyron; (3) Both Ns- and Xm-genome lineages served as the maternal donor during the speciation of Leymus species; (4) The Pseudoroegneria, Lophopyrum and Australopyrum genomes contributed to some Leymus species.

Phylogenetic Insights into Chinese Rubus (Rosaceae) from Multiple Chloroplast and Nuclear DNAs

Rubus L. is a large and taxonomically complex genus, species of which exhibit apomixis, polyploidy, and frequent hybridization. Most of Chinese Rubus are assigned in two major sections, Idaeobatus and Malachobatus. To explore the phylogenetic relationships within Chinese Rubus, inferences upon three chloroplast DNA (rbcL, rpl20-rps12, and trnG-trnS), nuclear ribosomal ITS, and two low-copy nuclear markers (GBSSI-2 and PEPC) were deduced in 142 Rubus taxa from 17 subsections in 6 sections. nrITS and GBSSI-2 were the most informative among the six DNA regions assessed. Phylogenetic relationships within Rubus were well-resolved by combined nuclear datasets rather than chloroplast markers. The phylogenetic inferences strongly supported that section Idaeobatus was a polyphyletic group with four distant clades. All samples of sect. Malachobatus formed a monophyletic clade, in which R. tsangorum and R. amphidasys of sect. Dalibardastrum, and R. peltatus from subsection Peltati of sect. Idaeobatus were always nested. Rubus pentagonus (2n = 2x = 14) from subsect. Alpestres of sect. Idaeobatus was a sister group to the polyploid sect. Malachobatus, as well as the polytomy of three sect. Cyalctis members. This suggests that some polyploids of Malachobatus might originate from common ancestors, via polyploidization of hybrids between R. pentagonus and sect. Cylactis species. They had experienced species explosion in a short time. Section Dalibardastrum species have potential parental lineages from subsects. Moluccani and Stipulosi of sect. Malachobatus. Based on molecular phylogenies, we also provided recommendations for the taxonomic treatments of four taxa. In addition, our results showed certain incongruence between chloroplast and nuclear markers, which might be due to hybridization and introgression.

The impacts of polyploidy, geographic and ecological isolations on the diversification of Panax (Araliaceae)

BACKGROUND

Panax L. is a medicinally important genus within family Araliaceae, where almost all species are of cultural significance for traditional Chinese medicine. Previous studies suggested two independent origins of the East Asia and North America disjunct distribution of this genus and multiple rounds of whole genome duplications (WGDs) might have occurred during the evolutionary process.

RESULTS

We employed multiple chloroplast and nuclear markers to investigate the evolution and diversification of Panax. Our phylogenetic analyses confirmed previous observations of the independent origins of disjunct distribution and both ancient and recent WGDs have occurred within Panax. The estimations of divergence time implied that the ancient WGD might have occurred before the establishment of Panax. Thereafter, at least two independent recent WGD events have occurred within Panax, one of which has led to the formation of three geographically isolated tetraploid species P. ginseng, P. japonicus and P. quinquefolius. Population genetic analyses showed that the diploid species P. notoginseng harbored significantly lower nucleotide diversity than those of the two tetraploid species P. ginseng and P. quinquefolius and the three species showed distinct nucleotide variation patterns at exon regions.

CONCLUSION

Our findings based on the phylogenetic and population genetic analyses, coupled with the species distribution patterns of Panax, suggested that the two rounds of WGD along with the geographic and ecological isolations might have together contributed to the evolution and diversification of this genus.

Tobamoviruses have probably co-diverged with their eudicotyledonous hosts for at least 110 million years

A phylogeny has been calculated by maximum likelihood comparisons of the concatenated consensus protein sequences of 29 tobamoviruses shown to be non-recombinant. This phylogeny has statistically significant support throughout, including its basal branches. The viruses form eight lineages that are congruent with the taxonomy of the hosts from which each was first isolated and, with the exception of three of the twenty-nine species, all fall into three clusters that have either asterid or rosid or caryophyllid hosts (i.e. the major subdivisions of eudicotyledonous plants). A modified Mantel permutation test showed that the patristic distances of virus and host phylogenies are significantly correlated, especially when the three anomalously placed viruses are removed. When the internal branches of the virus phylogeny were collapsed the congruence decreased. The simplest explanation of this congruence of the virus and host phylogenies is that most tobamovirus lineages have co-diverged with their primary plant hosts for more than 110 million years, and only the brassica-infecting lineage originated from a major host switch from asterids to rosids. Their co-divergence seems to have been ‘fuzzy’ rather than ‘strict’, permitting viruses to switch hosts within major host clades. Our conclusions support those of a coalesence analysis of tobamovirus sequences, that used proxy node dating, but not a similar analysis of nucleotide sequences from dated samples, which concluded that the tobamoviruses originated only 100 thousand years ago.

Congruent Deep Relationships in the Grape Family (Vitaceae) Based on Sequences of Chloroplast Genomes and Mitochondrial Genes via Genome Skimming

Vitaceae is well-known for having one of the most economically important fruits, i.e., the grape (Vitis vinifera). The deep phylogeny of the grape family was not resolved until a recent phylogenomic analysis of 417 nuclear genes from transcriptome data. However, it has been reported extensively that topologies based on nuclear and organellar genes may be incongruent due to differences in their evolutionary histories. Therefore, it is important to reconstruct a backbone phylogeny of the grape family using plastomes and mitochondrial genes. In this study,next-generation sequencing data sets of 27 species were obtained using genome skimming with total DNAs from silica-gel preserved tissue samples on an Illumina NextSeq 500 instrument [corrected]. Plastomes were assembled using the combination of de novo and reference genome (of V. vinifera) methods. Sixteen mitochondrial genes were also obtained via genome skimming using the reference genome of V. vinifera. Extensive phylogenetic analyses were performed using maximum likelihood and Bayesian methods. The topology based on either plastome data or mitochondrial genes is congruent with the one using hundreds of nuclear genes, indicating that the grape family did not exhibit significant reticulation at the deep level. The results showcase the power of genome skimming in capturing extensive phylogenetic data: especially from chloroplast and mitochondrial DNAs.

Resolution of Brassicaceae Phylogeny Using Nuclear Genes Uncovers Nested Radiations and Supports Convergent Morphological Evolution

Brassicaceae is one of the most diverse and economically valuable angiosperm families with widely cultivated vegetable crops and scientifically important model plants, such as Arabidopsis thaliana. The evolutionary history, ecological, morphological, and genetic diversity, and abundant resources and knowledge of Brassicaceae make it an excellent model family for evolutionary studies. Recent phylogenetic analyses of the family revealed three major lineages (I, II, and III), but relationships among and within these lineages remain largely unclear. Here, we present a highly supported phylogeny with six major clades using nuclear markers from newly sequenced transcriptomes of 32 Brassicaceae species and large data sets from additional taxa for a total of 55 species spanning 29 out of 51 tribes. Clade A consisting of Lineage I and Macropodium nivale is sister to combined Clade B (with Lineage II and others) and a new Clade C. The ABC clade is sister to Clade D with species previously weakly associated with Lineage II and Clade E (Lineage III) is sister to the ABCD clade. Clade F (the tribe Aethionemeae) is sister to the remainder of the entire family. Molecular clock estimation reveals an early radiation of major clades near or shortly after the Eocene–Oligocene boundary and subsequent nested divergences of several tribes of the previously polytomous Expanded Lineage II. Reconstruction of ancestral morphological states during the Brassicaceae evolution indicates prevalent parallel (convergent) evolution of several traits over deep times across the entire family. These results form a foundation for future evolutionary analyses of structures and functions across Brassicaceae.

Thousands of RAD-seq Loci Fully Resolve the Phylogeny of the Highly Disjunct Arctic-Alpine Genus Diapensia (Diapensiaceae)

Restriction-site associated DNA sequencing (RAD-seq) has recently become an important method to generate genome-wide molecular data for species delimitation, phylogeography, and population genetic studies. However, very few empirical studies have so far tested its applicability in phylogenetic reconstruction. The alpine-arctic genus Diapensia was selected to study the origin of the disjunction between the Arctic and the Himalayan-Hengduan Mountains (HHM). However, a previous phylogenetic analysis based on one nuclear and four plastid DNA regions failed to resolve the oldest divergences in Diapensia as well as the relationship between the two HHM species. Here we reconstruct a fully resolved phylogeny of Diapensia and address the conflict between the currently accepted taxonomy and the gene trees in the HHM species using RAD-seq. Based on a data set containing 2,650 loci selected to maximize the number of parsimony informative sites and allowing for a high level of missing data (51%), the phylogeny of Diapensia was fully resolved and each of the four species was reciprocally monophyletic. Whereas the arctic D. lapponica was inferred as sister to the HHM clade in the previous study, the RAD-seq data resolved the two arctic species as sisters to the HHM clade. Similar relationships were inferred from a differently filtered data set with far fewer loci (114) and less missing data (21%), but with lower support and with one of the two HHM species as non-monophyletic. Bayesian concordance analysis and Patterson’s D-statistic tests suggested that admixture has occurred between the two HHM species.

A genome-scale mining strategy for recovering novel rapidly-evolving nuclear single-copy genes for addressing shallow-scale phylogenetics in Hydrangea

Background

Identifying orthologous molecular markers that potentially resolve relationships at and below species level has been a major challenge in molecular phylogenetics over the past decade. Non-coding regions of nuclear low- or single-copy markers are a vast and promising source of data providing information for shallow-scale phylogenetics. Taking advantage of public transcriptome data from the One Thousand Plant Project (1KP), we developed a genome-scale mining strategy for recovering potentially orthologous single-copy markers to address low-scale phylogenetics. Our marker design targeted the amplification of intron-rich nuclear single-copy regions from genomic DNA. As a case study we used Hydrangea section Cornidia, one of the most recently diverged lineages within Hydrangeaceae (Cornales), for comparing the performance of three of these nuclear markers to other “fast” evolving plastid markers.

Results

Our data mining and filtering process retrieved 73 putative nuclear single-copy genes which are potentially useful for resolving phylogenetic relationships at a range of divergence depths within Cornales. The three assessed nuclear markers showed considerably more phylogenetic signal for shallow evolutionary depths than conventional plastid markers. Phylogenetic signal in plastid markers increased less markedly towards deeper evolutionary divergences. Potential phylogenetic noise introduced by nuclear markers was lower than their respective phylogenetic signal across all evolutionary depths. In contrast, plastid markers showed higher probabilities for introducing phylogenetic noise than signal at the deepest evolutionary divergences within the tribe Hydrangeeae (Hydrangeaceae).

Conclusions

While nuclear single-copy markers are highly informative for shallow evolutionary depths without introducing phylogenetic noise, plastid markers might be more appropriate for resolving deeper-level divergences such as the backbone relationships of the Hydrangeaceae family and deeper, at which non-coding parts of nuclear markers could potentially introduce noise due to elevated rates of evolution. The herein developed and demonstrated transcriptome based mining strategy has a great potential for the design of novel and highly informative nuclear markers for a range of plant groups and evolutionary scales.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0416-z) contains supplementary material, which is available to authorized users.

Phylogeography of Thlaspi arvense (Brassicaceae) in China Inferred from Chloroplast and Nuclear DNA Sequences and Ecological Niche Modeling

Thlaspi arvense is a well-known annual farmland weed with worldwide distribution, which can be found from sea level to above 4000 m high on the Qinghai-Tibetan Plateau (QTP). In this paper, a phylogeographic history of T. arvense including 19 populations from China was inferred by using three chloroplast (cp) DNA segments (trnL-trnF, rpl32-trnL and rps16) and one nuclear (n) DNA segment (Fe-regulated transporter-like protein, ZIP). A total of 11 chloroplast haplotypes and six nuclear alleles were identified, and haplotypes unique to the QTP were recognized (C4, C5, C7 and N4). On the basis of molecular dating, haplotypes C4, C5 and C7 have separated from others around 1.58 Ma for cpDNA, which corresponds to the QTP uplift. In addition, this article suggests that the T. arvense populations in China are a mixture of diverged subpopulations as inferred by hT/vT test (hT ≤ vT, cpDNA) and positive Tajima's D values (1.87, 0.05 < p < 0.10 for cpDNA and 3.37, p < 0.01 for nDNA). Multimodality mismatch distribution curves and a relatively large shared area of suitable environmental conditions between the Last Glacial Maximum (LGM) as well as the present time recognized by MaxEnt software reject the sudden expansion population model.

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

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

Selection and Validation of Housekeeping Genes as Reference for Gene Expression Studies in Pigeonpea (Cajanus cajan) under Heat and Salt Stress Conditions

To identify stable housekeeping genes as a reference for expression analysis under heat and salt stress conditions in pigeonpea, the relative expression variation for 10 commonly used housekeeping genes (EF1α, UBQ10, GAPDH, 18Sr RNA, 25Sr RNA, TUB6, ACT1, IF4α, UBC, and HSP90) was studied in root, stem, and leaves tissues of Asha (ICPL 87119), a leading pigeonpea variety. Three statistical algorithms geNorm, NormFinder, and BestKeeper were used to define the stability of candidate genes. Under heat stress, UBC, HSP90, and GAPDH were found to be the most stable reference genes. In the case of salinity stress, GAPDH followed by UBC and HSP90 were identified to be the most stable reference genes. Subsequently, the above identified genes were validated using qRT-PCR based gene expression analysis of two universal stress-resposive genes namely uspA and uspB. The relative quantification of these two genes varied according to the internal controls (most stable, least stable, and combination of most stable and least stable housekeeping genes) and thus confirmed the choice as well as validation of internal controls in such experiments. The identified and validated housekeeping genes will facilitate gene expression studies under heat and salt stress conditions in pigeonpea.

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

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

Hyb-Seq: Combining target enrichment and genome skimming for plant phylogenomics

PREMISE OF THE STUDY

Hyb-Seq, the combination of target enrichment and genome skimming, allows simultaneous data collection for low-copy nuclear genes and high-copy genomic targets for plant systematics and evolution studies. •

METHODS AND RESULTS

Genome and transcriptome assemblies for milkweed (Asclepias syriaca) were used to design enrichment probes for 3385 exons from 768 genes (>1.6 Mbp) followed by Illumina sequencing of enriched libraries. Hyb-Seq of 12 individuals (10 Asclepias species and two related genera) resulted in at least partial assembly of 92.6% of exons and 99.7% of genes and an average assembly length >2 Mbp. Importantly, complete plastomes and nuclear ribosomal DNA cistrons were assembled using off-target reads. Phylogenomic analyses demonstrated signal conflict between genomes. •

CONCLUSIONS

The Hyb-Seq approach enables targeted sequencing of thousands of low-copy nuclear exons and flanking regions, as well as genome skimming of high-copy repeats and organellar genomes, to efficiently produce genome-scale data sets for phylogenomics.

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

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

Nuclear ribosomal RNA genes: ITS region

Despite possible drawbacks (intraspecific polymorphisms and possible fungal contamination), sequencing of the ribosomal RNA gene ITS region remains one of the most popular nuclear sequences used for plant taxonomy and phylogeny. A protocol for PCR amplification and sequencing of this region using universal plant primers is provided.

Guidelines for the choice of sequences for molecular plant taxonomy

This chapter presents an overview of the major plant DNA sequences and molecular methods available for plant taxonomy. Guidelines are provided for the choice of sequences and methods to be used, based on the DNA compartment (nuclear, chloroplastic, mitochondrial), evolutionary mechanisms, and the level of taxonomic differentiation of the plants under survey.

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

BACKGROUND AND AIMS

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

METHODS

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

KEY RESULTS

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

CONCLUSIONS

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