Capturing single-copy nuclear genes,organellar genomes,and nuclear ribosomal DNA from deep genome skimming data for plant phylogenetics: A case study in Vitaceae

With the decreasing cost and availability of many newly developed bioinformatics pipelines, next-generation sequencing (NGS) has revolutionized plant systematics in recent years. Genome skimming has been widely used to obtain high-copy fractions of the genomes, including plastomes, mitochondrial DNA (mtDNA), and nuclear ribosomal DNA (nrDNA). In this study, through simulations, we evaluated the optimal (minimum) sequencing depth and performance for recovering single-copy nuclear genes (SCNs) from genome skimming data, by subsampling genome resequencing data and generating 10 data sets with different sequencing coverage in silico. We tested the performance of four data sets (plastome, nrDNA, mtDNA, and SCNs) obtained from genome skimming based on phylogenetic analyses of the Vitis clade at the genus level and Vitaceae at the family level, respectively. Our results showed that optimal minimum sequencing depth for high-quality SCNs assembly via genome skimming was about 10× coverage. Without the steps of synthesizing baits and enrichment experiments, coupled with incredibly low sequencing costs, we showcase that deep genome skimming (DGS) is as effective for capturing large data sets of SCNs as the widely used Hyb-Seq approach, in addition to capturing plastomes, mtDNA, and entire nrDNA repeats. DGS may serve as an efficient and economical alternative and may be superior to the popular target enrichment/Hyb-Seq approach.     

Nuclear and plastid phylogenomic analyses provide insights into the reticulate evolution,species delimitation,and biogeography of the Sino-Japanese disjunctive Diabelia (Caprifoliaceae)

Understanding biological diversity and the mechanisms of the Sino-Japanese disjunctions are major challenges in eastern Asia biogeography. The Sino-Japanese flora has been broadly studied as an ideal model for plant phylogeography. Diabelia Landrein (Caprifoliaceae) is an East Asian genus, with a disjunctive distribution across the Sino-Japanese region. However, relationships within Diabelia remain elusive. In this study, we reconstructed the phylogeny of Diabelia and inferred historical biogeography and evolutionary patterns based on nuclear and plastid sequences from target enrichment and genome skimming approaches, respectively. We found that the main clades within Diabelia were discordant between nuclear and plastid trees. Both nuclear and plastid phylogenetic analyses supported five main clades: Diabelia serrata (Siebold & Zucc.) Landrein, Diabelia tetrasepala (Koidz.) Landrein, Diabelia sanguinea (Makino) Landrein, Diabelia stenophylla (Honda) Landrein, and Diabelia spathulata (Siebold & Zucc.) Landrein. Species network analyses revealed that Diabelia tetrasepala is likely the result of a hybridization event. Divergence time estimation and ancestral area reconstructions showed that Diabelia originated in Japan during the early Miocene, with subsequent vicariance and dispersal events between Japan and Korea, and between Japan and China. Overall, our results support the division of Diabelia into five main clades and the recognition of five species in the genus. This research provides new insights into the species delimitation and speciation processes of taxonomically complex lineages such as Diabelia.     

Herbarium phylogenomics: Resolving the generic status of the enigmatic Pseudobartsia (Orobanchaceae)

The millions of herbarium specimens in collections around the world provide historical resources for phylogenomics and evolutionary studies. Many rare and endangered species exist only as historical specimens. Here, we report a case study of the monotypic Pseudobartsia yunnanensis D. Y. Hong (=Pseudobartsia glandulosa[Bentham] W. B. Yu & D. Z. Li: Orobanchaceae) known from a single Chinese collection taken in 1940. We obtained genomic data of Pseudobartsia glandulosa using high-throughput short-read sequencing, and then assembled a complete chloroplast genome and nuclear ribosome DNA region in this study. We found that the newly assembled three plastid DNA regions (atpB-rbcL, rpl16, and trnS-G) and nuclear ribosomal internal transcribed spacer (nrITS) of Pseudobartsia glandulosa were more than 99.98% similar to published sequences obtained by target sequencing. Phylogenies of Orobanchaceae using 30 plastomes (including 10 new plastomes), using both supermatrix and multispecies coalescent approaches following a novel plastid phylogenomic workflow, recovered seven recognized tribes and two unranked groups, both of which were proposed as new tribes, that is, Brandisieae and Pterygielleae. Within Pterygielleae, all analyses strongly supported Xizangia D. Y. Hong as the first diverging genus, with Pseudobartsia D. Y. Hong as sister to Pterygiella Oliver + Phtheirospermum Bunge (excluding Phtheirospermum japonicum [Thunberg] Kanitz); this supports reinstatement of Pseudobartsia and Xizangia. Although elements of Buchnereae-Cymbarieae-Orobancheae and Brandisieae-Pterygielleae-Rhinantheae showed incongruence among gene trees, the topology of the supermatrix tree was congruent with the majority of gene trees and functional-group trees. Therefore, most plastid genes are evolving as a linkage group, allowing the supermatrix tree approach to yield internally consistent phylogenies for Orobanchaceae.     

Phylogenomic relationships and species identification of the olive genus Olea (Oleaceae)

The olive genus Olea includes c. 30–40 taxa in three subgenera (Olea, Tetrapilus, and Paniculatae) within the family Oleaceae. Historically, the Olea genus was classified into four groups that were overall well supported by reconstructed phylogenies, despite incomplete sampling of subgenus Tetrapilus and poor resolution within clades. These analyses also showed that the genus was not monophyletic. Reliable identification of Olea species is important for both their conservation and utilization of this economically important genus. In this study, we used phylogenomic data from genome skimming to resolve relationships within Olea and to identify molecular markers for species identification. We assembled the complete plastomes, and nrDNA of 26 individuals representing 13 species using next-generation sequencing and added 18 publicly available accessions of Olea. We also developed nuclear SNPs using the genome skimming data to infer the phylogenetic relationships of Olea. Large-scale phylogenomic analyses of 138 samples of tribe Oleeae supported the polyphyly of Olea, with Olea caudatilimba and Olea subgenus Tetrapilus not sharing their most recent common ancestor with the main Olea clade (subgenus Paniculatae and subgenus Olea). The interspecific phylogenetic resolution was poor owing to a possible rapid radiation. By comparing with the plastome data, we identified the markers ycf1b and psbE-petL as the best Olea-specific chloroplast DNA barcodes. Compared with universal barcodes, specific DNA barcodes and super-barcode exhibited higher discriminatory power. Our results demonstrated the power of phylogenomics to improve phylogenetic relationships of intricate groups and provided new insights into barcodes that allow for accurate identification of Olea species.     

A decade of progress in plant molecular phylogenetics

Over the past decade, botanists have produced several thousand phylogenetic analyses based on molecular data, with particular emphasis on sequencing rbcL, the plastid gene encoding the large subunit of Rubisco (ribulose bisphosphate carboxylase). Because phylogenetic trees retrieved from the three plant genomes (plastid, nuclear and mitochondrial) have been highly congruent, the ‘Angiosperm Phylogeny Group’ has used these DNA-based phylogenetic trees to reclassify all families of flowering plants. However, in addition to taxonomy, these major phylogenetic efforts have also helped to define strategies to reconstruct the ‘tree of life’, and have revealed the size of the ancestral plant genome, uncovered potential candidates for the ancestral flower, identified molecular living fossils, and linked the rate of neutral substitutions with species diversity. With an increased interest in DNA sequencing programmes in non-model organisms, the next decade will hopefully see these phylogenetic findings integrated into new genetic syntheses, from genomes to taxa.     

Phylogenetic utility of the external transcribed spacer (ETS) of 18S-26S rDNA: congruence of ETS and ITS trees of Calycadenia (Compositae)

The 3' region of the external transcribed spacer (ETS) of 18S-26S nuclear ribosomal DNA was sequenced in 19 representatives of Calycadenia/Osmadenia and two outgroup species (Compositae) to assess its utility for phylogeny reconstruction compared to rDNA internal transcribed spacer (ITS) data. Universal primers based on plant, fungal, and animal sequences were designed to amplify the intergenic spacer (IGS) and an angiosperm primer was constructed to sequence the 3' end of the ETS in members of tribe Heliantheae. Based on these sequences, an internal ETS primer useful across Heliantheae sensu lato was designed to amplify and sequence directly the 3' ETS region in the study taxa, which were the subjects of an earlier phylogenetic investigation based on ITS sequences. Size variation in the amplified ETS region varied across taxa of Heliantheae sensu lato from approximately 350 to 700 bp, in part attributable to an approximately 200-bp tandem duplication in a common ancestor of Calycadenia/Osmadenia. Phylogenetic analysis of the 200-bp subrepeats and examination of apomorphic changes in the duplicated region demonstrate that the subrepeats in Calycadenia/Osmadenia have evolved divergently. Phylogenetic analyses of the entire amplified ETS region yielded a highly resolved strict consensus tree that is nearly identical in topology to the ITS tree, with strong bootstrap and decay support on most branches. Parsimony analyses of combined ETS and ITS data yielded a strict consensus tree that is better resolved and generally better supported than trees based on either data set analyzed separately. We calculated an approximately 1.3- to 2.4-fold higher rate of sequence evolution by nucleotide substitution in the ETS region studied than in ITS-1 + ITS-2. A similar disparity in the proportion of variable (1.3 ETS:1 ITS) and potentially informative (1.5 ETS:1 ITS) sites was observed for the ingroup. Levels of homoplasy are similar in the ETS and ITS data. We conclude that the ETS holds great promise for augmenting ITS data for phylogenetic studies of young lineages.     

DNA barcodes and microsatellites: How they complement for species identification in the complex genus Tamarix (Tamaricaceae)

DNA barcoding allows the identification of an organism by comparing the sequence of selected DNA regions (barcodes) with a previously compiled database, and it can be useful for taxonomic identification of species in complex genera, such as Tamarix. Many species of this genus show convergent morphology, which leads to frequent errors in their identification. Highly variable genetic markers, such as microsatellites or short sequence repeats (SSR), could be used to differentiate species where DNA barcodes fail. Here, we tested the ability of both, 5 different marker regions (rbcL, matK, ITS, trnH-psbA, and ycf1), and 14 microsatellites, to properly identify Tamarix species, especially those from the Mediterranean Basin, and compared the pros and cons of the different analytical methods for species identification. DNA barcoding allows the genetic identification of certain species in Tamarix. The two-locus barcodes matK + ITS and ITS + ycf1 were the best-performing combinations, allowing up to 69% and 70%, respectively, correct identification. However, DNA barcoding failed in phylogenetically close groups, such as many Mediterranean species. The use of SSR can aid the identification of species, and the combination of both types of data (DNA barcoding and SSR) improved the success. The combination of data was especially relevant in detecting the presence of hybridization processes, which are common in the genus. However, caution must be exercised when choosing the clustering methods for the SSR datasince different methods can lead to very different results.     

Resolving robust phylogenetic relationships of core Brassicaceae using genome skimming data

The Brassicaceae is an economically and scientifically important family distributed globally, including oilseed rape and the model plant, Arabidopsis thaliana. Although growing molecular data have been used in phylogenetic studies, the relationships among major clades and tribes of Brassicaceae are still controversial. Here, we investigated the core Brassicaceae phylogenetics using 222 plastomes and 235 nrDNA cistrons, including 106 plastomes and 112 nrDNA cistrons assembled from newly sequenced genome skimming data of 112 taxa. The sampling covered 73 genera from 61.5% tribes and four unassigned genera and species. Three well supported lineages LI, LII, and LIII were revealed in our plastomic analyses, with LI sister to LII + LIII. In addition, the monophyly of the newly delimitated LII was strongly supported by three different partition strategies, concatenated methods under Bayesian and Maximum Likelihood analyses. LII comprised 13 tribes, including four tribes previously unassigned to any lineage, that is Biscutelleae as the earliest diverging clade and Cochlearieae as the sister to Megacarpaeeae + Anastaticeae. Within LII, the intertribal relationships were also well resolved, except that a conflicting position of Orychophragmus was detected among different datasets. In LIII, Shehbazia was resolved as a member of Chorisproreae, but Chorisproreae, Dontostemoneae, and Euclidieae were all resolved as paraphyletic, which was also confirmed by nrDNA analyses. Moreover, the loss of the rps16 gene was detected as likely to be a synapomorphy of the tribes Arabideae and Alysseae. Overall, using genome skimming data, we resolved robust phylogenetic relationships of core Brassicaceae and shed new light on the complex evolutionary history of this family.     

Phylogenetic analyses suggest a hybrid origin of the figs (Moraceae: Ficus) that are endemic to the Ogasawara (Bonin) Islands, Japan

The Ogasawara Islands are oceanic islands and harbor a unique endemic flora. There are three fig species (Ficus boninsimae, F. nishimurae and F. iidaiana) endemic to the Ogasawara Islands, and these species have been considered to be closely related to Ficus erecta, and to have diverged within the islands. However, this hypothesis remains uncertain. To investigate this issue, we assessed the phylogenetic relationships of the Ogasawara figs and their close relatives occurring in Japan, Taiwan and South China based on six plastid genome regions, nuclear ITS region and two nuclear genes. The plastid genome-based tree indicated a close relationship between the Ogasawara figs and F. erecta, whereas some of the nuclear gene-based trees suggested this relationship was not so close. In addition, the phylogenetic analyses of the pollinating wasps associated with these fig species based on the nuclear 28S rRNA and mitochondrial cytB genes suggested that the fig-pollinating wasps of F. erecta are not sister to those of the Ogasawara figs These results suggest the occurrence of an early hybridization event(s) in the lineage leading to the Ogasawara figs.     

Out of the Middle East: New phylogenetic insights in the genus Tamarix (Tamaricaceae)

Tamarix is one of the taxonomically most complex genera among the angiosperms, and there is little consensus regarding its infrageneric classification. Here we present the most complete phylogenetic reconstruction of the genus to date. This includes a DNA phylogenetic tree based on nuclear ribosomal ITS, and a plastid DNA phylogeny based on three intergenic spacers (trnS‐trnG, ndhF‐rpl32, and trnQ‐rps16). In total, both nuclear and plastid phylogenetic analyses include more than 70 samples of 39 species from 27 countries, which represent close to 60% of the diversity of the genus. Two complementary trees, based only on one plastid marker, are also included. The first, based on trnS‐trnG, is used to increase the number of species related to T. amplexicaulis. The second, based on ndhF‐rpl32, is used to investigate the separation between T. tetrandra and T. parviflora. The incongruence between the available infrageneric classifications and the molecular results is confirmed. A reticulate evolution is inferred from the trees, showing characters such as vaginate leaves appearing at different stages along the evolutionary history of the genus. The presence of T. canariensis outside the Canary Islands is cast into doubt, and all such records from NW Africa and Europe are here considered to belong to T. gallica. The results also suggest independence of T. karelinii from T. hispida, and T. parviflora from T. tetrandra. Relationships between a number of species are still not resolved, and additional studies will be needed to further refine the complex taxonomy of Tamarix.     

Phylogenomic conflict analyses in the apple genus Malus s.l. reveal widespread hybridization and allopolyploidy driving diversification,with insights into the complex biogeographic history in the Northern Hemisphere

Phylogenomic evidence from an increasing number of studies has demonstrated that different data sets and analytical approaches often reconstruct strongly supported but conflicting relationships. In this study, 785 single-copy nuclear genes and 75 complete plastomes were used to infer the phylogenetic relationships and estimate the historical biogeography of the apple genus Malus sensu lato, an economically important lineage disjunctly distributed in the Northern Hemisphere and involved in known and suspected hybridization and allopolyploidy events. The nuclear phylogeny recovered the monophyly of Malus s.l. (including Docynia); however, the genus was supported to be biphyletic in the plastid phylogeny. An ancient chloroplast capture event in the Eocene in western North America best explains the cytonuclear discordance. Our conflict analysis demonstrated that ILS, hybridization, and allopolyploidy could explain the widespread nuclear gene tree discordance. One deep hybridization event (Malus doumeri) and one recent event (Malus coronaria) were detected in Malus s.l. Furthermore, our historical biogeographic analysis integrating living and fossil data supported a widespread East Asian-western North American origin of Malus s.l. in the Eocene, followed by several extinction and dispersal events in the Northern Hemisphere. We also propose a general workflow for assessing phylogenomic discordance and biogeographic analysis using deep genome skimming data sets.     

Phylogenetic information from three mitochondrial genomes of Terebelliformia (Annelida) worms and duplication of the methionine tRNA

Mitochondrial genomes have been useful for inferring animal phylogeny across a wide range of clades, however they are still poorly sampled in some animal taxa, limiting our knowledge of mtDNA evolution. For example, despite being one of the most diverse animal phyla, only 5 complete annelid mitochrondial genomes have been published. To address this paucity of information, we obtained complete mitochondrial genomic sequences from Pista cristata (Terebellidae) and Terebellides stroemi (Trichobranchidae) as well as one nearly complete mitochondrial genome from Eclysippe vanelli (Ampharetidae). These taxa are within Terebelliformia (Annelida), which include spaghetti worms, icecream cone worms and their relatives. In contrast to the 37 genes found in most bilaterian metazoans, we recover 38 genes in the mitochondrial genomes of T. stroemi and P. cristata due to the presence of a second methionine tRNA (trnM). Interestingly, the two trnMs are located next to each other and are possibly a synapomorphy of these two taxa. The E. vanelli partial mitochondrial genome lacks this additional trnM at the same position, but it may be present in the region not sampled. Compared to other annelids, gene orders of these three mitochondrial genomes are generally conserved except for the atp6-mSSU region. Phylogenetic analyses reveal that mtDNA data strongly supports a Trichobranchidae/Terebellidae clade.     

Hybrid genera in Liatrinae (Asteraceae: Eupatorieae)

Liatrinae is a small subtribe of Eupatorieae that occurs in North America with a center of generic-level diversity in the southeastern United States. Molecular phylogenetic data were sought to assess whether two monotypic genera, Garberia and Hartwrightia, are accurately placed in the subtribe, and to resolve questions of the generic-level classification of Carphephorus. Phylogenetic analyses of nuclear ITS/ETS and plastid DNA data indicated that Garberia is the basalmost diverging lineage, and that Hartwrightia is phylogenetically embedded in the subtribe. There was significant incongruence between the ITS/ETS and plastid DNA datasets in the placement of Hartwrightia and another monotypic genus, Litrisa, suggesting that both are of original hybrid origin. The results also showed that Carphephorus s.l. is not monophyletic, and even after removal of the two species of Trilisa, it is still paraphyletic to Liatris. The apparent hybrid origin of Hartwrightia, which is morphologically transgressive relative to its inferred parental lineages, suggests that reticulation between phylogenetically distinct lineages may be a recurrent problem for phylogenetic estimation in Asteraceae.     

鲁桑叶绿体基因组序列及特征分析

鲁桑(Morus multicaulis)是亚洲地区栽培的重要经济作物。以鲁桑品种日本胡橙为实验材料, 利用高通量测序技术对鲁桑叶绿体基因组进行测序, 获得NCBI登录号(KU355297), 并研究鲁桑的叶绿体基因组结构。结合前人对蒙桑(M. mongolica)、印度桑(M. indica)和川桑(M. notabilis)的研究结果, 对鲁桑的系统进化关系进行了探讨。研究结果表明: 鲁桑叶绿体基因组是一个典型的四部分结构, 全长159 154 bp, 共注释130个基因, 包含85个蛋白质编码基因(18个基因在反向重复区重复)、37个转运RNA (tRNA)基因和8个核糖体RNA (rRNA)基因。生物信息学分析表明, 在鲁桑中共搜索到82个SSR位点, 单核苷酸、二核苷酸、三核苷酸、四核苷酸和五核苷酸重复基序个数分别为63、7、2、9和1个, 并没有发现六核苷酸; 其中单核苷酸重复在鲁桑的叶绿体基因组SSR中占76.8%。采用MEGA 6.0软件, 通过最大似然法和近邻结合法对包括4个桑属物种在内的15个物种的叶绿体基因组序列进行聚类分析, 2种方法得到的聚类结果均为鲁桑和蒙桑聚在一起。研究结果对叶绿体基因组工程研究及桑属种间的分子标记开发和优良品种培育具有一定的参考价值。     

Reconciling gene and genome duplication events: using multiple nuclear gene families to infer the phylogeny of the aquatic plant family Pontederiaceae

Most plant phylogenetic inference has used DNA sequence data from the plastid genome. This genome represents a single genealogical sample with no recombination among genes, potentially limiting the resolution of evolutionary relationships in some contexts. In contrast, nuclear DNA is inherently more difficult to employ for phylogeny reconstruction because major mutational events in the genome, including polyploidization, gene duplication, and gene extinction can result in homologous gene copies that are difficult to identify as orthologs or paralogs. Gene tree parsimony (GTP) can be used to infer the rooted species tree by fitting gene genealogies to species trees while simultaneously minimizing the estimated number of duplications needed to reconcile conflicts among them. Here, we use GTP for five nuclear gene families and a previously published plastid data set to reconstruct the phylogenetic backbone of the aquatic plant family Pontederiaceae. Plastid-based phylogenetic studies strongly supported extensive paraphyly of Eichhornia (one of the four major genera) but also depicted considerable ambiguity concerning the true root placement for the family. Our results indicate that species trees inferred from the nuclear genes (alone and in combination with the plastid data) are highly congruent with gene trees inferred from plastid data alone. Consideration of optimal and suboptimal gene tree reconciliations place the root of the family at (or near) a branch leading to the rare and locally restricted E. meyeri. We also explore methods to incorporate uncertainty in individual gene trees during reconciliation by considering their individual bootstrap profiles and relate inferred excesses of gene duplication events on individual branches to whole-genome duplication events inferred for the same branches. Our study improves understanding of the phylogenetic history of Pontederiaceae and also demonstrates the utility of GTP for phylogenetic analysis.     

Complex evolution in Arundinarieae (Poaceae: Bambusoideae): incongruence between plastid and nuclear GBSSI gene phylogenies

The monophyly of tribe Arundinarieae (the temperate woody bamboos) has been unequivocally recovered in previous molecular phylogenetic studies. In a recent phylogenetic study, 10 major lineages in Arundinarieae were resolved based on eight non-coding plastid regions, which conflicted significantly with morphological classifications both at the subtribal and generic levels. Nevertheless, relationships among and within the 10 lineages remain unclear. In order to further unravel the evolutionary history of Arundinarieae, we used the nuclear GBSSI gene sequences along with those of eight plastid regions for phylogenetic reconstruction, with an emphasis on Chinese species. The results of the plastid analyses agreed with previous studies, whereas 13 primary clades revealed in the GBSSI phylogeny were better resolved at the generic level than the plastid phylogeny. Our analyses also revealed many inconsistencies between the plastid DNA and the nuclear GBSSI trees. These results implied that the nuclear genome and the plastid genome had different evolutionary trajectories. The patterns of incongruence suggested that lack of informative characters, incomplete lineage sorting, and/or hybridization (introgression) could be the causes. Seven putative hybrid species were hypothesized, four of which are discussed in detail on the basis of topological incongruence, chromosome numbers, morphology, and distribution patterns, and those taxa probably resulted from homoploid hybrid speciation. Overall, our study indicates that the tribe Arundinarieae has undergone a complex evolution.