Against all odds: reconstructing the evolutionary history of Scrophularia (Scrophulariaceae) despite high levels of incongruence and reticulate evolution

Using RAD-seq to identify and differentiate the medicinal herb Scrophularia ningpoensis and its adulterants

BACKGROUND

Scrophulariae Radix, the dried root of Scrophularia ningpoensis Hemsl., belongs to the Scrophulariaceae family and is known for its diverse pharmacological effects, such as antihypertensive, anti-inflammatory, and antitumor properties. It is common that the Scrophulariae Radix is often contaminated by various adulterants like S. buergeriana, S. kakudensis, and S. yoshimurae in its medicinal usage. However, traditional methods for identifying Chinese medicinal materials have been found to be insufficient in accurately distinguishing S. ningpoensis from the three adulterant species.

RESULTS

In this study, RAD-seq technology was employed to sequence 27 individuals from the four species, yielding 55,250 high-quality SNP sites. Genetic differentiation analysis using these sites revealed significant genetic divergence between S. ningpoensis and the three adulterants. Through assessment genetic structure, principal component, and phylogenetic analysis, these four species were effectively distinguished with a high degree of confidence. The phylogenetic results indicated that S. ningpoensis is more closely related to S. yoshimurae, while S. buergeriana shows a closer relationship with S. kakudensis, providing new reference information for understanding interspecific relationships among the four Scrophularia species.

CONCLUSIONS

This study successfully identified S. ningpoensis and the three adulterants, enriched the genomic information of Scrophularia species, and contributed to a more profound understanding of the evolutionary dynamics within this genus. Furthermore, it offers novel insights and methodologies for the identification of medicinal plants.

Comparative plastome analyses and evolutionary relationships of 25 East Asian species within the medicinal plant genus Scrophularia (Scrophulariaceae)

Backgroud

Scrophularia L., a genus of the Scrophulariaceae, is a group of important medicinal plants used for eliminating heat and detoxifying. East Asia has an abundance of potentially medicinal Scrophularia species, and it serves as a secondary diversity center of the genus. However, the genomic resources available for germplasm identification and pharmaceutical exploration of East Asian Scrophularia are insufficient, hindering its commercial and industrial development. Additionally, the interspecific relationships of most East Asian Scrophularia species remain unclear.

Methods

In this study, we sequenced the leaves of 25 East Asian species of the genus Scrophularia, assembled and annotated the complete chloroplast genomes, and subsequently performed comparative and phylogenetic analyses on these genomes.

Results and discussion

The conserved plastome length of these 25 species ranged from 151,582 bp to 153,239 bp, containing a total of 132 coding genes, including 18 duplicated genes and 114 unique genes. Through genome alignment of these 25 species, 38-53 repeated sequences and 7 shared SSRs were identified, along with regions with high nucleotide polymorphism (Pi), which could potentially serve as molecular markers for species identification. The genome structure, gene content, and arrangement showed conservation, while variations were observed in the IR boundary regions and IGS. Phylogenetic inferences based on whole plastomes or on coding sequences (CDS) only yielded congruent results. We categorized the 25 East Asian Scrophularia species into six distinct clades and further explored their interspecies relationships using morphological characteristics, such as flower color, the relative position of stamens and corolla, and plant height. This could lay a genetic basis for future resource development of Scrophularia in East Asia.

Scalable neighbour search and alignment with uvaia

Despite millions of SARS-CoV-2 genomes being sequenced and shared globally, manipulating such data sets is still challenging, especially selecting sequences for focused phylogenetic analysis. We present a novel method, uvaia, which is based on partial and exact sequence similarity for quickly extracting database sequences similar to query sequences of interest. Many SARS-CoV-2 phylogenetic analyses rely on very low numbers of ambiguous sites as a measure of quality since ambiguous sites do not contribute to single nucleotide polymorphism (SNP) differences. Uvaia overcomes this limitation by using measures of sequence similarity which consider partially ambiguous sites, allowing for more ambiguous sequences to be included in the analysis if needed. Such fine-grained definition of similarity allows not only for better phylogenetic analyses, but could also lead to improved classification and biogeographical inferences. Uvaia works natively with compressed files, can use multiple cores and efficiently utilises memory, being able to analyse large data sets on a standard desktop.

Phylogenomic analyses reveal reticulate evolution between Neomicrocalamus and Temochloa (Poaceae: Bambusoideae)

Neomicrocalamus and Temochloa are closely related to bamboo genera. However, when considered with newly discovered and morphologically similar material from China and Vietnam, the phylogenetic relationship among these three groups was ambiguous in the analyses based on DNA regions. Here, as a means of investigating the relationships among the three bamboo groups and exploring potential sources of genomic conflicts, we present a phylogenomic examination based on the whole plastome, single-nucleotide polymorphism (SNP), and single-copy nuclear (SCN) gene datasets. Three different phylogenetic hypotheses were found. The inconsistency is attributed to the combination of incomplete lineage sorting and introgression. The origin of newly discovered bamboos is from introgressive hybridization between Temochloa liliana (which contributed 80.7% of the genome) and Neomicrocalamus prainii (19.3%), indicating that the newly discovered bamboos are closer to T. liliana in genetics. The more similar morphology and closer distribution elevation also imply a closer relationship between Temochloa and newly discovered bamboos.

Phylogenomics sheds new light on the drivers behind a long-lasting systematic riddle: the figwort family Scrophulariaceae

The figwort family, Scrophulariaceae, comprises c. 2000 species whose evolutionary relationships at the tribal level have proven difficult to resolve, hindering our ability to understand their origin and diversification. We designed a specific probe kit for Scrophulariaceae, targeting 849 nuclear loci and obtaining plastid regions as by-products. We sampled c. 87% of the genera described in the family and use the nuclear dataset to estimate evolutionary relationships, timing of diversification, and biogeographic patterns. Ten tribes, including two new tribes, Androyeae and Camptolomeae, are supported, and the phylogenetic positions of Androya, Camptoloma, and Phygelius are unveiled. Our study reveals a major diversification at c. 60 million yr ago in some Gondwanan landmasses, where two different lineages diversified, one of which gave rise to nearly 81% of extant species. A Southern African origin is estimated for most modern-day tribes, with two exceptions, the American Leucophylleae, and the mainly Australian Myoporeae. The rapid mid-Eocene diversification is aligned with geographic expansion within southern Africa in most tribes, followed by range expansion to tropical Africa and multiple dispersals out of Africa. Our robust phylogeny provides a framework for future studies aimed at understanding the role of macroevolutionary patterns and processes that generated Scrophulariaceae diversity.

Characterization and Comparative Analysis of Chloroplast Genomes of Medicinal Herb Scrophularia ningpoensis and Its Common Adulterants (Scrophulariaceae)

Scrophularia ningpoensis, a perennial medicinal plant from the Scrophulariaceae family, is the original species of Scrophulariae Radix (SR) in the Chinese Pharmacopoeia. This medicine is usually deliberately substituted or accidentally contaminated with other closely related species including S. kakudensis, S. buergeriana, and S. yoshimurae. Given the ambiguous identification of germplasm and complex evolutionary relationships within the genus, the complete chloroplast genomes of the four mentioned Scrophularia species were sequenced and characterized. Comparative genomic studies revealed a high degree of conservation in genomic structure, gene arrangement, and content within the species, with the entire chloroplast genome spanning 153,016-153,631 bp in full length, encoding 132 genes, including 80 protein-coding genes, 4 rRNA genes, 30 tRNA genes, and 18 duplicated genes. We identified 8 highly variable plastid regions and 39-44 SSRs as potential molecular markers for further species identification in the genus. The consistent and robust phylogenetic relationships of S. ningpoensis and its common adulterants were firstly established using a total of 28 plastid genomes from the Scrophulariaceae family. In the monophyletic group, S. kakudensis was determined to be the earliest diverging species, succeeded by S. ningpoensis. Meanwhile, S. yoshimurae and S. buergeriana were clustered together as sister clades. Our research manifestly illustrates the efficacy of plastid genomes in identifying S. ningpoensis and its counterfeits and will also contribute to a deeper understanding of the evolutionary processes within Scrophularia.

Comparative and Phylogenetic Analyses of Complete Chloroplast Genomes of Scrophularia incisa Complex (Scrophulariaceae)

The Scrophularia incisa complex is a group of closely related desert and steppe subshrubs that includes S. incisa, S. kiriloviana and S. dentata, which are the only S. sect. Caninae components found in Northwest China. Based on earlier molecular evidence, the species boundaries and phylogenetic relationships within this complex remain poorly resolved. Here, we characterized seven complete chloroplast genomes encompassing the representatives of the three taxa in the complex and one closely related species, S. integrifolia, as well as three other species of Scrophularia. Comparative genomic analyses indicated that the genomic structure, gene order and content were highly conserved among these eleven plastomes. Highly variable plastid regions and simple sequence repeats (SSRs) were identified. The robust and consistent phylogenetic relationships of the S. incisa complex were firstly constructed based on a total of 26 plastid genomes from Scrophulariaceae. Within the monophyletic complex, a S. kiriloviana individual from Pamirs Plateau was identified as the earliest diverging clade, followed by S. dentata from Tibet, while the remaining individuals of S. kiriloviana from the Tianshan Mountains and S. incisa from Qinghai-Gansu were clustered into sister clades. Our results evidently demonstrate the capability of plastid genomes to improve phylogenetic resolution and species delimitation, particularly among closely related species, and will promote the understanding of plastome evolution in Scrophularia.

Phylogeny and Historical Biogeography of Veronica Subgenus Pentasepalae (Plantaginaceae): Evidence for Its Origin and Subsequent Dispersal

Simple Summary

The Irano-Turanian phytogeographical region is considered a biodiversity reservoir for adjacent regions. The present phylogeographic study suggests that Veronica subgenus Pentasepalae originated in the Iranian plateau and was dispersed via a North African route to the Mediterranean and the Euro-Siberian regions. These findings highlight the importance of the Iranian plateau as a center of origin for many temperate plant species. Our results also resolve several taxonomic and phylogenetic issues surrounding the Southwest Asian species of this subgenus.

Abstract

Veronica subgenus Pentasepalae is the largest subgenus of Veronica in the Northern Hemisphere with approximately 80 species mainly from Southwest Asia. In order to reconstruct the phylogenetic relationships among the members of V. subgenus Pentasepalae and to test the “out of the Iranian plateau” hypothesis, we applied thorough taxonomic sampling, employing nuclear DNA (ITS) sequence data complimented with morphological studies and chromosome number counts. Several high or moderately supported clades are reconstructed, but the backbone of the phylogenetic tree is generally unresolved, and many Southwest Asian species are scattered along a large polytomy. It is proposed that rapid diversification of the Irano-Turanian species in allopatric glacial refugia and a relatively high rate of extinction during interglacial periods resulted in such phylogenetic topology. The highly variable Asian V. orientalis–V. multifida complex formed a highly polyphyletic assemblage, emphasizing the idea of cryptic speciation within this group. The phylogenetic results allow the re-assignment of two species into this subgenus. In addition, V. bombycina subsp. bolkardaghensis, V. macrostachya subsp. schizostegia and V. fuhsii var. linearis are raised to species rank and the new name V. parsana is proposed for the latter. Molecular dating and ancestral area reconstructions indicate a divergence age of about 9 million years ago and a place of origin on the Iranian Plateau. Migration to the Western Mediterranean region has likely taken place through a North African route during early quaternary glacial times. This study supports the assumption of the Irano-Turanian region as a source of taxa for neighboring regions, particularly in the alpine flora.

Plio-Pleistocene climatic change drives allopatric speciation and population divergence within the Scrophularia incisa complex (Scrophulariaceae) of desert and steppe subshrubs in Northwest China

Numerous temperate plants and animals on the Qinghai-Tibet Plateau (QTP) are hypothesized to have differentiated due to vicariant allopatric speciation associated with the geologic uplifts. However, this hypothesis has rarely been tested through a phylogeographic study of relative species in a broader geographic context, including the QTP, Tianshan Mountains, Mongolian Plateau, and surrounding regions. To understand the speciation and diversification process of plants across this wide area, phylogeographic analysis were examined from Scrophularia incisa and two other closely relative species comprising S. kiriloviana and S. dentata. Thirty-two populations of the three close relatives were genotyped using chloroplast DNA fragments and nuclear microsatellite loci to assess population structure and diversity, supplemented by phylogenetic dating, ancestral area reconstructions and species distribution modelings, as well as niche identity tests. Our chloroplast DNA (cpDNA) phylogeny showed that this monophyletic group of desert and steppe semi-shrub is derived from a Middle Pliocene ancestor of the Central Asia. Lineages in Central Asia vs. China diverged through climate/tectonic-induced vicariance during Middle Pliocene. Genetic and ENM data in conjunction with niche differentiation analyses support that the divergence of S. incisa, S. dentata and S. kiriloviana in China lineage proceeded through allopatric speciation, might triggered by early Pleistocene climate change of increase of aridification and enlargement of deserts, while subsequent climate-induced cycles of range contractions/expansions enhanced the geographical isolation and habit fragmentation of these taxa. These findings highlight the importance of the Plio-Pleistocene climate change in shaping genetic diversity and driving speciation in temperate steppes and deserts of Northwestern China.

Origin of the critically endangered endemic species Scrophularia takesimensis (Scrophulariaceae) on Ulleung Island, Korea: implications for conservation

Scrophularia takesimensis is a critically endangered endemic species of Ulleung Island, Korea. A previous molecular phylogenetic study based on nuclear ribosomal DNA (nrDNA) internal transcribed spacer (ITS) sequences with very limited sampling suggested that it is most closely related to the clade comprising S. alata and S. grayanoides. To determine the origin of S. takesimensis, we sampled a total of 171 accessions including S. takesimensis (9 populations and 63 individuals) and two closely related species, S. alata (11 populations and 68 individuals) and S. grayanoides (5 populations and 40 individuals) from eastern Asia and sequenced ITS and two chloroplast DNA (cpDNA) non-coding regions. Previously sequenced representative species of Scrophularia (109 taxa for ITS and 80 taxa for cpDNA) were combined with our data set and analyzed. While the global scale ITS phylogenetic tree suggests monophyly for each of the three eastern Asian species, S. takesimensis appears to be more closely related (albeit weakly) to a clade containing eastern North American/Caribbean species than to either S. alata or S. grayanoides. By contrast, the global scale cpDNA phylogenetic tree demonstrates that the eastern North America/Caribbean clade is sister to a clade comprising the three eastern Asian species. In addition, the monophyletic S. takesimensis is deeply embedded within paraphyletic S. alata, sharing its most recent common ancestor with populations from Japan/Sakhalin. Two divergent, geographically structured cp haplotype groups within S. takesimensis suggest at least two independent introductions from different source areas. A new and accurate chromosome number of S. takesimensis (2n = 94) is reported and some conservation strategies are discussed.

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

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