Plastid phylogenomics of the Gynoxoid group (Senecioneae, Asteraceae) highlights the importance of motif-based sequence alignment amid low genetic distances

On the importance of sequence alignment inspections in plastid phylogenomics - an example from revisiting the relationships of the water-lilies

The water-lily clade represents the second earliest-diverging branch of angiosperms. Most of its species belong to Nymphaeaceae, of which the "core Nymphaeaceae"-comprising the genera Euryale, Nymphaea and Victoria-is the most diverse clade. Despite previous molecular phylogenetic studies on the core Nymphaeaceae, various aspects of their evolutionary relationships have remained unresolved. The length-variable introns and intergenic spacers are known to contain most of the sequence variability within the water-lily plastomes. Despite the challenges with multiple sequence alignment, any new molecular phylogenetic investigation on the core Nymphaeaceae should focus on these noncoding plastome regions. For example, a new plastid phylogenomic study on the core Nymphaeaceae should generate DNA sequence alignments of all plastid introns and intergenic spacers based on the principle of conserved sequence motifs. In this investigation, we revisit the phylogenetic history of the core Nymphaeaceae by employing such an approach. Specifically, we use a plastid phylogenomic analysis strategy in which all coding and noncoding partitions are separated and then undergo software-driven DNA sequence alignment, followed by a motif-based alignment inspection and adjustment. This approach allows us to increase the reliability of the character base compared to the default practice of aligning complete plastomes through software algorithms alone. Our approach produces significantly different phylogenetic tree reconstructions for several of the plastome regions under study. The results of these reconstructions underscore that Nymphaea is paraphyletic in its current circumscription, that each of the five subgenera of Nymphaea is monophyletic, and that the subgenus Nymphaea is sister to all other subgenera of Nymphaea. Our results also clarify many evolutionary relationships within the Nymphaea subgenera Brachyceras, Hydrocallis and Nymphaea. In closing, we discuss whether the phylogenetic reconstructions obtained through our motif-based alignment adjustments are in line with morphological evidence on water-lily evolution.

Generic concepts and species diversity within the Gynoxyoid clade (Senecioneae, Compositae)

The Gynoxyoid clade of the Senecioneae (Asteraceae) until now included the five genera Aequatorium, Gynoxys, Nordenstamia, Paracalia and Paragynoxys as diagnosed using selected morphological characters. In their pre-phylogenetic circumscription, the genera Aequatorium and Paragynoxys were considered to inhabit the northern Andes in contrast to Nordenstamia and Paracalia that occur in the central Andes. The most species-rich genus, Gynoxys, was believed to be distributed throughout the Andes. We use a recently established plastid phylogenomic framework that rendered Gynoxys paraphyletic to further evaluate the delimitation of genera in the Gynoxyoid clade. We examine the morphological variation of all members of the Gynoxyoid to identify characters potentially informative at genus level. This results in a matrix of eleven, mostly multistate characters, including those originally used to diagnose these genera. The ancestral character state inference displays a high level of homoplasy, but nevertheless supports the recognition of four genera. Aequatorium is characterised by white radiate capitula. Paracalia and Paragynoxys share white flowers and floral characteristics, such as flower opening and length of disc flowers lobes, as plesiomorphic states, but differ in habit (scandent shrubs vs. trees). Paracalia also retained white flowers, but its two species are characterised by the absence of outer phyllaries. The genera Gynoxys and Nordenstamia comprise species with yellow capitula which appear to be a derived feature in the Gynoxyoids. The genus Nordenstamia, with eight species, is synonymised under Gynoxys since molecular evidence shows its species nested within various parts of the Gynoxys subclade and the morphological variation of Nordenstamia falls well within that of Gynoxys. With the goal to assign all species to four genera (Aequatorium, Gynoxys, Paracalia and Paragynoxys), we assess the states for the eleven characters for all members of the Gynoxyoids and generate new ETS and ITS sequences for 171 specimens belonging to 49 species to further support their generic placement. We provide a taxonomic treatment for the four genera recognised here including amended diagnoses and morphological descriptions. Furthermore, a species-level taxonomic backbone is elaborated for all genera using electronic tools that list 158 currently accepted names and synonyms (209 names in total) with the respective protologue and type information, as well as notes on the current understanding of species limits. Eleven names are newly synonymised, two are lectotypified and eight are newly transferred to other genera.

Plastid phylogenomics and species discrimination in the "Chinese" clade of Roscoea (Zingiberaceae)

Roscoea is an alpine or subalpine genus from the pan-tropical family Zingiberaceae, which consists of two disjunct groups in geography, namely the "Chinese" clade and the "Himalayan" clade. Despite extensive research on the genus, Roscoea species remain poorly defined and relationships between these species are not well resolved. In this study, we used plastid genomes of nine species and one variety to resolve phylogenetic relationships within the "Chinese" clade of Roscoea and as DNA super barcodes for species discrimination. We found that Roscoea plastid genomes ranged in length from 163,063 to 163,796 bp, and encoded 113 genes, including 79 protein-coding genes, 30 tRNA genes, four rRNA genes. In addition, expansion and contraction of the IR regions showed obvious infraspecific conservatism and interspecific differentiation. Plastid phylogenomics revealed that species belonging to the "Chinese" clade of Roscoea can be divided into four distinct subclades. Furthermore, our analysis supported the independence of R. cautleoides var. pubescens, the recovery of Roscoea pubescens Z.Y. Zhu, and a close relationship between R. humeana and R. cautloides. When we used the plastid genome as a super barcode, we found that it possessed strong discriminatory power (90%) with high support values. Intergenic regions provided similar resolution, which was much better than that of protein-coding regions, hypervariable regions, and DNA universal barcodes. However, plastid genomes could not completely resolve Roscoea phylogeny or definitively discriminate species. These limitations are likely related to the complex history of Roscoea speciation, poorly defined species within the genus, and the maternal inheritance of plastid genomes.

Software Choice and Sequencing Coverage Can Impact Plastid Genome Assembly-A Case Study in the Narrow Endemic Calligonum bakuense

Most plastid genome sequences are assembled from short-read whole-genome sequencing data, yet the impact that sequencing coverage and the choice of assembly software can have on the accuracy of the resulting assemblies is poorly understood. In this study, we test the impact of both factors on plastid genome assembly in the threatened and rare endemic shrub Calligonum bakuense. We aim to characterize the differences across plastid genome assemblies generated by different assembly software tools and levels of sequencing coverage and to determine if these differences are large enough to affect the phylogenetic position inferred for C. bakuense compared to congeners. Four assembly software tools (FastPlast, GetOrganelle, IOGA, and NOVOPlasty) and seven levels of sequencing coverage across the plastid genome (original sequencing depth, 2,000x, 1,000x, 500x, 250x, 100x, and 50x) are compared in our analyses. The resulting assemblies are evaluated with regard to reproducibility, contig number, gene complement, inverted repeat length, and computation time; the impact of sequence differences on phylogenetic reconstruction is assessed. Our results show that software choice can have a considerable impact on the accuracy and reproducibility of plastid genome assembly and that GetOrganelle produces the most consistent assemblies for C. bakuense. Moreover, we demonstrate that a sequencing coverage between 500x and 100x can reduce both the sequence variability across assembly contigs and computation time. When comparing the most reliable plastid genome assemblies of C. bakuense, a sequence difference in only three nucleotide positions is detected, which is less than the difference potentially introduced through software choice.

New Insights Into the Relationships Within Subtribe Scorzonerinae (Cichorieae, Asteraceae) Using Hybrid Capture Phylogenomics (Hyb-Seq)

Subtribe Scorzonerinae (Cichorieae, Asteraceae) contains 12 main lineages and approximately 300 species. Relationships within the subtribe, either at inter- or intrageneric levels, were largely unresolved in phylogenetic studies to date, due to the lack of phylogenetic signal provided by traditional Sanger sequencing markers. In this study, we employed a phylogenomics approach (Hyb-Seq) that targets 1,061 nuclear-conserved ortholog loci designed for Asteraceae and obtained chloroplast coding regions as a by-product of off-target reads. Our objectives were to evaluate the potential of the Hyb-Seq approach in resolving the phylogenetic relationships across the subtribe at deep and shallow nodes, investigate the relationships of major lineages at inter- and intrageneric levels, and examine the impact of the different datasets and approaches on the robustness of phylogenetic inferences. We analyzed three nuclear datasets: exon only, excluding all potentially paralogous loci; exon only, including loci that were only potentially paralogous in 1-3 samples; exon plus intron regions (supercontigs); and the plastome CDS region. Phylogenetic relationships were reconstructed using both multispecies coalescent and concatenation (Maximum Likelihood and Bayesian analyses) approaches. Overall, our phylogenetic reconstructions recovered the same monophyletic major lineages found in previous studies and were successful in fully resolving the backbone phylogeny of the subtribe, while the internal resolution of the lineages was comparatively poor. The backbone topologies were largely congruent among all inferences, but some incongruent relationships were recovered between nuclear and plastome datasets, which are discussed and assumed to represent cases of cytonuclear discordance. Considering the newly resolved phylogenies, a new infrageneric classification of Scorzonera in its revised circumscription is proposed.

Plastid Phylogenomic Analysis of Tordylieae Tribe (Apiaceae, Apioideae)

Based on the nrDNA ITS sequence data, the Tordylieae tribe is recognized as monophyletic with three major lineages: the subtribe Tordyliinae, the Cymbocarpum clade, and the Lefebvrea clade. Recent phylogenomic investigations showed incongruence between the nuclear and plastid genome evolution in the tribe. To assess phylogenetic relations and structure evolution of plastomes in Tordylieae, we generated eleven complete plastome sequences using the genome skimming approach and compared them with the available data from this tribe and close relatives. Newly assembled plastomes had lengths ranging from 141,148 to 150,103 base pairs and contained 122–127 genes, including 79–82 protein-coding genes, 35–37 tRNAs, and 8 rRNAs. We observed substantial differences in the inverted repeat length and gene content, accompanied by a complex picture of multiple J_LA and J_LB shifts. In concatenated phylogenetic analyses, Tordylieae plastomes formed at least three not closely related lineages with plastomes of the Lefebvrea clade as a sister group to plastomes from the Selineae tribe. The newly obtained data have increased our knowledge on the range of plastome variability in Apiaceae.