Patterns and causes of incongruence between plastid and nuclear Senecioneae (Asteraceae) phylogenies

The challenging biogeography of the Juan Fernández Islands and Coast Range of central Chile explained by new models of East Pacific tectonics

Biogeographers have often been puzzled by several unusual features in the Juan Fernández Islands (JFI) biota. These include the very high endemism density, multiple endemics that are older than the current islands, close biogeographic affinities with the central and West Pacific, and affinities with the diverse Coast Range of central Chile. We review aspects of biogeography in the JFI and the Coast Range in light of recent geological studies. These have examined the mantle below the East Pacific and South America, and have produced radical, new ideas on tectonic history. A long-lived, intraoceanic archipelago ~9000 km long is now thought to have existed in the East Pacific (passing between the JFI hotspot and mainland Chile) until the mid-Cretaceous. At this time, South America, which was moving westward with the opening of the Atlantic, collided with the archipelago. The assumption that the JFI biota is no older than its current islands is questionable, as taxa would have survived on prior islands produced at the JFI hotspot. We propose a new interpretation of evolution in the region based on tectonics rather than on island age and incorporating the following factors: the newly described East Pacific Archipelago; a long history for the JFI hotspot; metapopulation dynamics, including metapopulation vicariance; and formation of the Humboldt Current in the Cretaceous. The model accounts for many distinctive features of the JFI and Coast Range biota.

Tracing the evolutionary and genetic footprints of atmospheric tillandsioids transition from land to air

Plant evolution is driven by key innovations of functional traits that enables their survivals in diverse ecological environments. However, plant adaptive evolution from land to atmospheric niches remains poorly understood. In this study, we use the epiphytic Tillandsioideae subfamily of Bromeliaceae as model plants to explore their origin, evolution and diversification. We provide a comprehensive phylogenetic tree based on nuclear transcriptomic sequences, indicating that core tillandsioids originated approximately 11.3 million years ago in the Andes. The geological uplift of the Andes drives the divergence of tillandsioids into tank-forming and atmospheric types. Our genomic and transcriptomic analyses reveal gene variations and losses associated with adaptive traits such as impounding tanks and absorptive trichomes. Furthermore, we uncover specific nitrogen-fixing bacterial communities in the phyllosphere of tillandsioids as potential source of nitrogen acquisition. Collectively, our study provides integrative multi-omics insights into the adaptive evolution of tillandsioids in response to elevated aerial habitats.

Phylogenomics and phylogeographic model testing using convolutional neural networks reveal a history of recent admixture in the Canarian Kleinia neriifolia

Multiple-island endemics (MIE) are considered ideal natural subjects to study patterns of island colonization that involve recent population-level genetic processes. Kleinia neriifolia is a Canarian MIE widespread across the archipelago, which exhibits a close phylogenetic relationship with species in northwest Africa and at the other side of the Sahara Desert. Here, we used target sequencing with plastid skimming (Hyb-Seq), a dense population-level sampling of K. neriifolia, and representatives of its African-southern Arabian relatives to infer phylogenetic relationships and divergence times at the species and population levels. Using population genetic techniques and machine learning (convolutional neural networks [CNNs]), we reconstructed phylogeographic relationships and patterns of genetic admixture based on a multilocus SNP nuclear dataset. Phylogenomic analysis based on the nuclear dataset identifies the northwestern African Kleinia anteuphorbium as the sister species of K. neriifolia, with divergence starting in the early Pliocene. Divergence from its sister clade, comprising species from the Horn of Africa and southern Arabia, is dated to the arid Messinian period, lending support to the climatic vicariance origin of the Rand Flora. Phylogeographic model testing with CNNs supports an initial colonization of the central island of Tenerife followed by eastward and westward migration across the archipelago, which resulted in the observed east/west phylogeographic split. Subsequent population extinctions linked to aridification events, and recolonization from Tenerife, are proposed to explain the patterns of genetic admixture in the eastern Canary Islands. We demonstrate that CNNs based on SNPs can be used to discriminate among complex scenarios of island migration and colonization.

Island biogeography of the megadiverse plant family Asteraceae

The megadiverse plant family Asteraceae forms an iconic component of island floras including many spectacular radiations, but a global picture of its insular diversity is lacking. Here, we uncover the global biogeographical and evolutionary patterns of Asteraceae on islands to reveal the magnitude and potential causes of their evolutionary success. We compile a global checklist of Asteraceae species native and endemic to islands and combine it with macroecological analyses and a phylogenetic review of island radiations. Asteraceae have a global distribution on islands, comprising approximately 6,000 native island species, with 58% endemics. While diversity of the family on islands is lower than expected given its overall diversity, Asteraceae are the most diverse family on oceanic islands, suggesting an exceptional ability to thrive in isolation. In agreement with island biogeography predictions, native Asteraceae diversity increases with area and decreases with isolation, while endemism increases with both. We identify 39 confirmed island radiations and 69 putative radiations, exceeding numbers for other iconic insular groups. Our results reveal Asteraceae offer immense potential for research in ecology and evolution, given their close tracking of island biogeography expectations, large number of both species and radiations, cosmopolitan distribution, and numerous undiscovered radiations.

Comparative analysis of complete chloroplast genomes of Synotis species (Asteraceae, Senecioneae) for identification and phylogenetic analysis

BACKGROUND

The Synotis (C. B. Clarke) C. Jeffrey & Y. L. Chen is an ecologically important genus of the tribe Senecioneae, family Asteraceae. Because most species of the genus bear similar morphology, traditional morphological identification methods are very difficult to discriminate them. Therefore, it is essential to develop a reliable and effective identification method for Synotis species. In this study, the complete chloroplast (cp.) genomes of four Synotis species, S. cavaleriei (H.Lév.) C. Jeffrey & Y.L. Chen, S. duclouxii (Dunn) C. Jeffrey & Y.L. Chen, S. nagensium (C.B. Clarke) C. Jeffrey & Y.L. Chen and S. erythropappa (Bureau & Franch.) C. Jeffrey & Y. L. Chen had been sequenced using next-generation sequencing technology and reported here.

RESULTS

These four cp. genomes exhibited a typical quadripartite structure and contained the large single-copy regions (LSC, 83,288 to 83,399 bp), the small single-copy regions (SSC, 18,262 to 18,287 bp), and the inverted repeat regions (IR, 24,837 to 24,842 bp). Each of the four cp. genomes encoded 134 genes, including 87 protein-coding genes, 37 tRNA genes, 8 rRNA genes, and 2 pseudogenes (ycf1 and rps19). The highly variable regions (trnC-GCA-petN, ccsA-psaC, trnE-UUC-rpoB, ycf1, ccsA and petN) may be used as potential molecular barcodes. The complete cp. genomes sequence of Synotis could be used as the potentially effective super-barcode to accurately identify Synotis species. Phylogenetic analysis demonstrated that the four Synotis species were clustered into a monophyletic group, and they were closed to the Senecio, Crassocephalum and Dendrosenecio in tribe Senecioneae.

CONCLUSIONS

This study will be useful for further species identification, evolution, genetic diversity and phylogenetic studies within this genus Synotis and the tribe Senecioneae.

Insights into comparative genomics, structural features, and phylogenetic relationship of species from Eurasian Aster and its related genera (Asteraceae: Astereae) based on complete chloroplast genome

Aster L. is an economically and phylogenetically important genus in the tribe Astereae. Here, the complete plastomes of the eight Aster species were assembled and characterized using next-generation sequencing datasets. The results indicated the complete plastomes of Aster had a quadripartite structure. These genomes were 152,045-152,729 bp in length and contained 132-133 genes, including 87 protein-coding genes, 37-38 tRNA genes, and eight rRNA genes. Expansion or contraction of inverted repeat regions and forward, palindromic, complement, and reverse repeats were detected in the eight Aster species. Additionally, our analyses showed the richest type of simple sequence repeats was A/T mononucleotides, and 14 highly variable regions were discovered by analyzing the border regions, sequence divergence, and hotspots. Phylogenetic analyses indicated that 27 species in Astereae were clustered into six clades, i.e., A to D, North American, and outgroup clades, and supported that the genera Heteropappus, Kalimeris, and Heteroplexis are nested within Aster. The results indicated the clades B to D might be considered as genera. Divergence time estimate showed the clades A, B, C, and D diverged at 23.15 Mya, 15.13 Mya, 24.29 Mya, and 21.66 Mya, respectively. These results shed light on the phylogenetic relationships of Aster and provided new information on species identification of Aster and its related genera.

Atmospheric CO2 decline and the timing of CAM plant evolution

BACKGROUND AND AIMS

CAM photosynthesis is hypothesized to have evolved in atmospheres of low CO2 concentration in recent geological time because of its ability to concentrate CO2 around Rubisco and boost water use efficiency relative to C3 photosynthesis. We assess this hypothesis by compiling estimates of when CAM clades arose using phylogenetic chronograms for 73 CAM clades. We further consider evidence of how atmospheric CO2 affects CAM relative to C3 photosynthesis.

RESULTS

Where CAM origins can be inferred, strong CAM is estimated to have appeared in the past 30 million years in 46 of 48 examined clades, after atmospheric CO2 had declined from high (near 800 ppm) to lower (<450 ppm) values. In turn, 21 of 25 clades containing CAM species (but where CAM origins are less certain) also arose in the past 30 million years. In these clades, CAM is probably younger than the clade origin. We found evidence for repeated weak CAM evolution during the higher CO2 conditions before 30 million years ago, and possible strong CAM origins in the Crassulaceae during the Cretaceous period prior to atmospheric CO2 decline. Most CAM-specific clades arose in the past 15 million years, in a similar pattern observed for origins of C4 clades.

CONCLUSIONS

The evidence indicates strong CAM repeatedly evolved in reduced CO2 conditions of the past 30 million years. Weaker CAM can pre-date low CO2 and, in the Crassulaceae, strong CAM may also have arisen in water-limited microsites under relatively high CO2. Experimental evidence from extant CAM species demonstrates that elevated CO2 reduces the importance of nocturnal CO2 fixation by increasing the contribution of C3 photosynthesis to daily carbon gain. Thus, the advantage of strong CAM would be reduced in high CO2, such that its evolution appears less likely and restricted to more extreme environments than possible in low CO2.

Synotisjinpingensis (Asteraceae, Senecioneae), a new species with white ray florets from southeastern Yunnan, China

Synotisjinpingensis (Asteraceae, Senecioneae), a new species from Jinping county in southeastern Yunnan province, China, is described and illustrated. This species is distinguished by having white ray florets in the genus Synotis, in which only species with yellow ray florets have been hitherto known. In habit and leaf shape S.jinpingensis is most closely similar to S.duclouxii, a species occurring in southwestern Guizhou, southern Sichuan and northeastern Yunnan, China, but differs, in addition to the color of ray florets, by having fewer lateral veins of leaves, obviously longer bracts of calyculus, and larger phyllaries. The membership of the new species within Synotis is strongly corroborated by evidence from floral micromorphology and phylogenetic analyses based on ITS sequence data. Color photographs of living plants, a distribution map, and provisional IUCN status of S.jinpingensis are provided.

Ancient hybridization in Curcuma (Zingiberaceae)-Accelerator or brake in lineage diversifications?

Hybridization is a widespread phenomenon in the evolution of plants and exploring its role is crucial to understanding diversification processes of many taxonomic groups. Recently, more attention is focused on the role of ancient hybridization that has repeatedly been shown as triggers of evolutionary radiation, although in some cases, it can prevent further diversification. The causes, frequency, and consequences of ancient hybridization remain to be explored. Here, we present an account of several events of ancient hybridization in turmeric, the economically important plant genus Curcuma (Zingiberaceae), which harbors about 130 known species. We analyzed 1094 targeted low-copy genes and plastomes obtained by next-generation sequencing of 37 species of Curcuma, representing the known genetic diversity and spanning the geographical distribution of the genus. Using phylogenetic network analysis, we show that the entire genus Curcuma as well as its most speciose lineage arose via introgression from the genus Pyrgophyllum and one of the extinct lineages, respectively. We also document a single event of ancient hybridization, with C. vamana as a product, that represents an evolutionary dead end. We further discuss distinct circumstances of those hybridization events that deal mainly with (in)congruence in chromosome counts of the parental lineages.

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.

Historical development of karst evergreen broadleaved forests in East Asia has shaped the evolution of a hemiparasitic genus Brandisia (Orobanchaceae)

Brandisia is a shrubby genus of about eight species distributed basically in East Asian evergreen broadleaved forests (EBLFs), with distribution centers in the karst regions of Yunnan, Guizhou, and Guangxi in southwestern China. Based on the hemiparasitic and more or less liana habits of this genus, we hypothesized that its evolution and distribution were shaped by the development of EBLFs there. To test our hypothesis, the most comprehensive phylogenies of Brandisia hitherto were constructed based on plastome and nuclear loci (nrDNA, PHYA and PHYB); then divergence time and ancestral areas were inferred using the combined nuclear loci dataset. Phylogenetic analyses reconfirmed that Brandisia is a member of Orobanchaceae, with unstable placements caused by nuclear-plastid incongruences. Within Brandisia, three major clades were well supported, corresponding to the three subgenera based on morphology. Brandisia was inferred to have originated in the early Oligocene (32.69 Mya) in the Eastern Himalayas-SW China, followed by diversification in the early Miocene (19.45 Mya) in karst EBLFs. The differentiation dates of Brandisia were consistent with the origin of keystone species of EBLFs in this region (e.g., Fagaceae, Lauraceae, Theaceae, and Magnoliaceae) and the colonization of other characteristic groups (e.g., Gesneriaceae and Mahonia). These findings indicate that the distribution and evolution of Brandisia were facilitated by the rise of the karst EBLFs in East Asia. In addition, the woody and parasitic habits, and pollination characteristics of Brandisia may also be the important factors affecting its speciation and dispersal.

Plastome evolution and phylogenomic insights into the evolution of Lysimachia (Primulaceae: Myrsinoideae)

BACKGROUND

Lysimachia L., the second largest genus within the subfamily Myrsinoideae of Primulaceae, comprises approximately 250 species worldwide. China is the species diversity center of Lysimachia, containing approximately 150 species. Despite advances in the backbone phylogeny of Lysimachia, species-level relationships remain poorly understood due to limited genomic information. This study analyzed 50 complete plastomes for 46 Lysimachia species. We aimed to identify the plastome structure features and hypervariable loci of Lysimachia. Additionally, the phylogenetic relationships and phylogenetic conflict signals in Lysimachia were examined.

RESULTS

These fifty plastomes within Lysimachia had the typical quadripartite structure, with lengths varying from 152,691 to 155,784 bp. Plastome size was positively correlated with IR and intron length. Thirteen highly variable regions in Lysimachia plastomes were identified. Additionally, ndhB, petB and ycf2 were found to be under positive selection. Plastid ML trees and species tree strongly supported that L. maritima as sister to subg. Palladia + subg. Lysimachia (Christinae clade), while the nrDNA ML tree clearly placed L. maritima and subg. Palladia as a sister group.

CONCLUSIONS

The structures of these plastomes of Lysimachia were generally conserved, but potential plastid markers and signatures of positive selection were detected. These genomic data provided new insights into the interspecific relationships of Lysimachia, including the cytonuclear discordance of the position of L. maritima, which may be the result of ghost introgression in the past. Our findings have established a basis for further exploration of the taxonomy, phylogeny and evolutionary history within Lysimachia.

Phylogenomics of the Olea europaea complex using 15 whole genomes supports recurrent genetic admixture together with differentiation into seven subspecies

BACKGROUND

The last taxonomic account of Olea recognises six subspecies within Olea europaea L., including the Mediterranean olive tree (subsp. europaea) and five other subspecies (laperrinei, guanchica, maroccana, cerasiformis, and cuspidata) distributed across the Old World, including Macaronesian islands. The evolutionary history of this monophyletic group (O. europaea complex) has revealed a reticulated scenario involving hybridization and polyploidization events, leading to the presence of a polyploid series associated with the subspecies. However, how the polyploids originated, and how the different subspecies contributed to the domestication of the cultivated olive are questions still debated. Tracing the recent evolution and genetic diversification of the species is key for the management and preservation of its genetic resources. To study the recent history of the O. europaea complex, we compared newly sequenced and available genomes for 27 individuals representing the six subspecies.

RESULTS

Our results show discordance between current subspecies distributions and phylogenomic patterns, which support intricate biogeographic patterns. The subspecies guanchica, restricted to the Canary Islands, is closely related to subsp. europaea, and shows a high genetic diversity. The subsp. laperrinei, restricted now to high mountains of the Sahara desert, and the Canarian subsp. guanchica contributed to the formation of the allotetraploid subsp. cerasiformis (Madeira islands) and the allohexaploid subsp. maroccana (western Sahara region). Our phylogenomic data support the recognition of one more taxon (subsp. ferruginea) for the Asian populations, which is clearly segregated from the African subsp. cuspidata.

CONCLUSIONS

In sum, the O. europaea complex underwent several processes of hybridization, polyploidy, and geographical isolation resulting in seven independent lineages with certain morphological traits recognised into subspecies.

Plastome evolution and phylogeny of the tribe Ruteae (Rutaceae)

Rutaceae is a large family, and the genus-level classification in the subfamilies or tribes of this family is not unified based on different taxonomic treatments. Until now, phylogenetic relationships of some genera in traditional tribe Ruteae have not been clearly resolved. In this study, seven new complete plastomes of this tribe were sequenced, and a comparative analysis was performed to investigate their plastome characteristics and evolution. In addition, we inferred the phylogenetic relationships of Ruteae based on complete plastome and nuclear ITS data. All plastomes exhibited a typical quadripartite structure and were relatively conserved in their structure and gene arrangement. Their genome sizes ranged from 154,656 bp to 160,677 bp, and the size variation was found to be associated with differences in IR expansion and gene loss. A total of 112 to 114 genes were identified in the genomes, including 78 to 79 protein-coding genes, 30 tRNA genes, 4 rRNA genes, and 2 pseudogenes. Sequence divergence analysis indicated that non-coding regions exhibited a higher percentage of variable characters, and nine non-coding and six coding regions were identified as divergent hotspots. Phylogenetic results based on different datasets showed that this tribe was divided into three reciprocally exclusive groups. The phylogenetic analyses between plastome and nuclear ITS data were partly incongruent with each other. This study provides new insights into plastome evolution of Ruteae as well as Rutaceae. The availability of these plastomes provides useful genomic resources for molecular DNA barcodes and phylogenetically informative markers and deepens our understanding of the phylogeny in Ruteae.

Molecular and Morphometric Update on Italian Salicornia (Chenopodiaceae), with a Focus on the Species S. procumbens s. l

Salicornia is a highly taxonomically problematic genus due to the reduced morphological observable characters. Ten Eurasian species are currently recognized: S. alpini, S. europaea, S. fruticosa, S. hispanica, S. lagascae, S. perennans, S. perennis, S. persica, S. procumbens, and S. pruinosa. In addition, eleven subspecies are accepted, mainly based on their distribution areas. Along the Venetian coasts and in Sardinia, in the past, an endemic species called S. veneta was recognized, but this name was later synonymized with S. procumbens subsp. procumbens. The aim of the present research is investigating different Italian Salicornia populations by a molecular point of view, using the nuclear ribosomal external transcribed spacer ETS and the plastid psbA-trnH intergenic spacer. A particular focus is on the comparison between Venetian (including those occurring in locus classicus of S. veneta) and Sardinian S. procumbens and other Italian populations of this species. The molecular analyses based on the plastid marker highlight that the Italian S. procumbens populations form two well distinct groups. In particular, some of the Venetian (Locus classicus of S. veneta) and all the Sardinian specimens are genetically distinct (=plastid haplotype 1) from the other investigated populations (=plastid haplotype 2). This indicates that the psbA-trnH haplotype 1 glassworts represent a distinct entity, which we suppose to coincide with the former S. veneta. Therefore, we suggest to recognize this taxonomic entity at the subspecies rank, as S. procumbens subsp. veneta comb. and stat. nov. However, contrary to the results found with the plastid psbA-trnH intergenic spacer, the ETS locus does not show a separation into two distinct clades for S. procumbens, probably due to a different evolution of the two loci. Nevertheless, in the ETS phylogenetic reconstruction, the Sardinian specimens (=ribotypes 2 and 3) are placed, together with a Moroccan sample, in a subclade separated from all the other S. procumbens. These results suggest that the Sardinian populations can represent a subspecies/incipient speciation process, probably due to geographic isolation. In the light of this, morphometric analyses (k-means, MANOVA, PCA, DA, and Box-Plot) have been carried out on the Sardinian and Venetian populations to verify if this distinction is detectable also by a morphological point of view. The morphometric analyses highlight the existence of two groups, concerning both the nuclear and plastid trees. Six characters were found to be diagnostic.

Characterization and phylogenetic analysis of the complete chloroplast genome of Emilia sonchifolia (L.) DC

Emilia sonchifolia is a herb with antioxidant, anti-inflammatory, antitumor, and wound healing properties. The complete chloroplast genome (cp genome) of the genus Emilia was sequenced for the first time. The cp genome of E. sonchifolia is 151,474 bp in length. It contained a large single-copy (LSC) region (84,004 bp), and small single-copy (SSC) region (17,980 bp), and two inverted repeats (IRs, 24,745 bp). Phylogenetic analysis of 24 species was conducted. E. sonchifolia was found to be closely related to Pericallis hybrida and Dendrosenecio spp. The sequenced cp genome would be useful to understand the phylogeny and genomic studies of the genus Emilia.

Evidence from Phylogenomics and Morphology Provide Insights into the Phylogeny, Plastome Evolution, and Taxonomy of Kitagawia

Kitagawia Pimenov is one of the segregate genera of Peucedanum sensu lato within the Apiaceae. The phylogenetic position and morphological delimitation of Kitagawia have been controversial. In this study, we used plastid genome (plastome) and nuclear ribosomal DNA (nrDNA) sequences to reconstruct the phylogeny of Kitagawia, along with comparative plastome and morphological analyses between Kitagawia and related taxa. The phylogenetic results identified that all examined Kitagawia species were divided into Subclade I and Subclade II within the tribe Selineae, and they were all distant from the representative members of Peucedanum sensu stricto. The plastomes of Kitagawia and related taxa showed visible differences in the LSC/IRa junction (JLA) and several hypervariable regions, which separated Subclade I and Subclade II from other taxa. Fruit anatomical and micromorphological characteristics, as well as general morphological characteristics, distinguished the four Kitagawia species within Subclade I from Subclade II and other related genera. This study supported the separation of Kitagawia from Peucedanum sensu lato, confirmed that Kitagawia belongs to Selineae, and two species (K. praeruptora and K. formosana) within Subclade II should be placed in a new genus. We believe that the "core" Kitagawia should be limited to Subclade I, and this genus can be distinguished by the association of a series of morphological characteristics. Overall, our study provides new insights into the phylogeny, plastome evolution, and taxonomy of Kitagawia.

The phylogeny of Seseli (Apiaceae, Apioideae): insights from molecular and morphological data

BACKGROUND

The genus Seseli L., which consists of 125-140 species distributed in the Old World from western Europe and northwestern Africa to China and Japan, is one of the largest and most taxonomically difficult genera of Apiaceae Lindl. Although several previous studies have been conducted on Seseli based on limited morphological characteristics and molecular fragments, a robust and comprehensive phylogeny of Seseli remains elusive. Plastomes provide abundant genetic information and have been widely used in studying plant phylogeny and evolution. Consequently, we newly generated the complete plastomes of eleven Seseli taxa. We combined plastome data and morphological characteristics to investigate the phylogeny of Seseli.

RESULTS

In our study, we observed that the genome length, gene numbers, IR/SC borders, and repeat composition of the eleven Seseli plastomes were variable. Several appropriate mutation hotspot regions may be developed as candidate DNA barcodes for evolution, phylogeny, and species identification of Seseli. The phylogenetic results identified that Seseli was not a monophyletic group. Moreover, the eleven newly sequenced Seseli taxa did not cluster with S. tortuosum (the type species of Seseli, belonging to the tribe Selineae), where S. delavayi clustered with Eriocycla belonging to the tribe Echinophoreae and the other ten belonged to Selineae. The comparative plastome and morphological characteristics analyses confirmed the reliability of the phylogenetic analyses and implied the complex evolution of Seseli.

CONCLUSION

Combining molecular and morphological data is efficient and useful for studying the phylogeny of Seseli. We suggest that "a narrow sense" of Seseli will be meaningful for further study and the current taxonomic system of Seseli needs to be revised. In summary, our study can provide new insights into the phylogenetic relationships and taxonomic framework of Seseli.

A target capture approach for phylogenomic analyses at multiple evolutionary timescales in rosewoods (Dalbergia spp.) and the legume family (Fabaceae)

Understanding the genetic changes associated with the evolution of biological diversity is of fundamental interest to molecular ecologists. The assessment of genetic variation at hundreds or thousands of unlinked genetic loci forms a sound basis to address questions ranging from micro- to macroevolutionary timescales, and is now possible thanks to advances in sequencing technology. Major difficulties are associated with (i) the lack of genomic resources for many taxa, especially from tropical biodiversity hotspots; (ii) scaling the numbers of individuals analysed and loci sequenced; and (iii) building tools for reproducible bioinformatic analyses of such data sets. To address these challenges, we developed target capture probes for genomic studies of the highly diverse, pantropically distributed and economically significant rosewoods (Dalbergia spp.), explored the performance of an overlapping probe set for target capture across the legume family (Fabaceae), and built the general purpose bioinformatic pipeline CaptureAl. Phylogenomic analyses of Malagasy Dalbergia species yielded highly resolved and well supported hypotheses of evolutionary relationships. Population genomic analyses identified differences between closely related species and revealed the existence of a potentially new species, suggesting that the diversity of Malagasy Dalbergia species has been underestimated. Analyses at the family level corroborated previous findings by the recovery of monophyletic subfamilies and many well-known clades, as well as high levels of gene tree discordance, especially near the root of the family. The new genomic and bioinformatic resources, including the Fabaceae1005 and Dalbergia2396 probe sets, will hopefully advance systematics and ecological genetics research in legumes, and promote conservation of the highly diverse and endangered Dalbergia rosewoods.

Multiplex PCR Identification of Aspergillus cristatus and Aspergillus chevalieri in Liupao Tea Based on Orphan Genes

“Golden flower” fungi in dark tea are beneficial to human health. The rapid identification method of “golden flower” fungi can verify the quality of dark tea products and ensure food safety. In this study, 6 strains were isolated from Liupao tea. They were respectively identified as A. cristatus, A. chevalieri, and A. pseudoglaucus. A. pseudoglaucus was reported as Liupao tea “golden flower” fungus for the first time. It was found that the ITS and BenA sequences of A. cristatus and A. chevalieri were highly conserved. It is difficult to clearly distinguish these closely related species by ITS sequencing. To rapidly identify species, multiplex PCR species-specific primers were designed based on orphan genes screened by comparative genomics analysis. Multiplex PCR results showed that orphan genes were specific and effective for the identification of A. cristatus and A. chevalieri isolated from Liupao tea and Fu brick tea. We confirmed that orphan genes can be used for identification of closely related Aspergillus species. Validation showed that the method is convenient, rapid, robust, sequencing-free, and economical. This promising method will be greatly beneficial to the dark tea processing industry and consumers.

The complete chloroplast genome of Gynura japonica and its phylogenetic implications

Gynura japonica (Asteraceae) is a folk herbal medicine with multi-pharmacological functions involving analgesic, hemostatic and antiangiogenic activities. The study was conducted to assemble the complete chloroplast (cp) genome of G. japonica through a genome-skimming approach. The assembled cp genome was 151,023 bp in size, with 62.8% AT content, consisting of a large single copy (LSC) of 83,185 bp, two copies of inverted repeat (IRa and IRb) of 24,847 bp, and a small single copy (SSC) of 18,144 bp. The cp genome of G. japonica contained 133 genes, including eight ribosomal RNA genes (rRNAs), 37 transfer RNA genes (tRNAs), 86 protein-coding genes (PCGs), and two pseudogenes (ψycf1 and ψrps19). Our phylogenomic analysis based on whole plastid genomes strongly supports G. japonica is a sister to the clade including Crassocephalum crepidioides and Jacobaea vulgaris.

Phylogenetic position of Ligusticopsis (Apiaceae, Apioideae): evidence from molecular data and carpological characters

Ligusticopsis (Apiaceae, Apioideae) is now considered to have an East-Asia and Sino-Himalaya distribution. The genus was not recognized as a natural and separate genus and was treated as a synonym of Ligusticum both in Flora Reipublicae Popularis Sinicae and Flora of China since first established, though Pimenov et al. have made many taxonomic revisions to Ligusticopsis, phylogenetic relationships between Ligusticopsis and Ligusticum have been in dispute. Thirty-four plastomes and 35 nrITS from Apioideae were analysed by RAxML and MrBayes to reconstruct the phylogenetic relationships, along with carpology of 10 species and comparative analyses of 17 plastomes to investigate the evidence supporting the independence of Ligusticopsis. As a result, nine species suggested to be Ligusticopsis formed a highly supported monophyletic branch (Subclade A) inside Selineae both in maximum likelihood and Bayesian inference; the results of the comparative analyses further supported the monophyly of Subclade A, mainly in the location of genes at the IRa/LSC boundary, the sequence diversity exhibited by various genes (e.g. trnH-GUG-psbA and ycf2) and same codon biases in terminator TAA (relative synonymous codon usage = 1.75). Species in Subclade A also had shared characters in mericarps, combined with other characters of the plant, 'base clothed in fibrous remnant sheaths, pinnate bracts, pinnate bracteoles longer than rays of umbellule, mericarps strongly compressed dorsally, median and lateral ribs filiform or keeled, marginal ribs winged, and numerous vittae in commissure and each furrow' should be the most important and diagnostic characters of Ligusticopsis. Our phylogenetic trees and other analyses supported the previous taxonomic treatments of Pimenov et al. that Ligusticopsis should be a natural and separate genus rather than a synonym of Ligusticum.

Elsholtziazhongyangii (Lamiaceae), a new species from Sichuan, China

Elsholtziazhongyangii (Lamiaceae), a new species from Sichuan Province, China, is described and illustrated. The new species is morphologically similar to E.feddeif.feddei, but it can be easily distinguished from E.feddeif.feddei by smaller corolla (3.2-3.5 mm vs. 4.5-5.3 mm), bract indumentum (glabrous, except margin ciliate vs. villous, especially on veins abaxially, glabrous adaxially) and bract stalked (ca. 1.2 mm vs. sessile). Phylogenetic analyses, based on two nuclear ribosomal (ETS, ITS) and five plastid (rbcL, matK, trnL-F, ycf1, ycf1-rps15) regions, confirmed that the new species formed a monophyletic clade with robust support. The new species is currently known from western Sichuan.

Phylogenomics and Biogeography of Populus Based on Comprehensive Sampling Reveal Deep-Level Relationships and Multiple Intercontinental Dispersals

Populus not only has significant economic and ecological values, but also serves as a model tree that is widely used in the basic research of tree growth, physiology, and genetics. However, high levels of morphological variation and extensive interspecific hybridization of Populus pose an obstacle for taxonomy, and also to the understanding of phylogenetic interspecific relationships and biogeographical history. In this study, a total of 103 accessions representing almost all the wild species of Populus were collected and whole-genome re-sequenced to examine the phylogenetic relationships and biogeography history. On the basis of 12,916,788 nuclear single nucleotide polymorphisms (SNPs), we reconstructed backbone phylogenies using concatenate and coalescent methods, we highly disentangled the species relationships of Populus, and several problematic taxa were treated as species complexes. Furthermore, the phylogeny of the chloroplast genome showed extensive discordance with the trees from the nuclear genome data, and due to extensive chloroplast capture and hybridization of Populus species, plastomes could not accurately evaluate interspecies relationships. Ancient gene flow between clades and some hybridization events were also identified by ABBA-BABA analysis. The reconstruction of chronogram and ancestral distributions suggested that North America was the original region of this genus, and subsequent long dispersal and migration across land bridges were contributed to the modern range of Populus. The diversification of Populus mainly occurred in East Asia in recent 15 Ma, possibly promoted by the uplift of the Tibetan Plateau. This study provided comprehensive evidence on the phylogeny of Populus and proposed a four-subgeneric classification and a new status, subgenus Abaso. Meanwhile, ancestral distribution reconstruction with nuclear data advanced the understanding of the biogeographic history of Populus.

Phylogenetic incongruence in Cymbidium orchids

Cymbidium, which includes approximately 80 species, is one of the most ornamental and cultivated orchid genera. However, a lack of markers and sparse sampling have posed great challenges to resolving the phylogenetic relationships within the genus. In the present study, we reconstructed the phylogenetic relationships by utilizing one nuclear DNA (nrITS) and seven plastid genes (rbcL, trnS, trnG, matK, trnL, psbA, and atpI) from 70 species (varieties) in Cymbidium. We also examined the occurrence of phylogenetic conflict between nuclear (nrITS) and plastid loci and investigated how phylogenetic conflict bears on taxonomic classification within the genus. We found that phylogenetic conflict and low support values may be explained by hybridization and a lack of informative characteristics. Our results do not support previous classification of the subgenera and sections within Cymbidium. Discordance between gene trees and network analysis indicate that reticulate evolution occurred in the genus Cymbidium. Overall, our study indicates that Cymbidium has undergone a complex evolution.

Chloranthus genome provides insights into the early diversification of angiosperms

Chloranthales remain the last major mesangiosperm lineage without a nuclear genome assembly. We therefore assemble a high-quality chromosome-level genome of Chloranthus spicatus to resolve enigmatic evolutionary relationships, as well as explore patterns of genome evolution among the major lineages of mesangiosperms (eudicots, monocots, magnoliids, Chloranthales, and Ceratophyllales). We find that synteny is highly conserved between genomic regions of Amborella, Vitis, and Chloranthus. We identify an ancient single whole-genome duplication (WGD) (κ) prior to the divergence of extant Chloranthales. Phylogenetic inference shows Chloranthales as sister to magnoliids. Furthermore, our analyses indicate that ancient hybridization may account for the incongruent phylogenetic placement of Chloranthales + magnoliids relative to monocots and eudicots in nuclear and chloroplast trees. Long genes and long introns are found to be prevalent in both Chloranthales and magnoliids compared to other angiosperms. Overall, our findings provide an improved context for understanding mesangiosperm relationships and evolution and contribute a valuable genomic resource for future investigations.

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

PREMISE

The genus Gynoxys and relatives form a species-rich lineage of Andean shrubs and trees with low genetic distances within the sunflower subtribe Tussilaginineae. Previous molecular phylogenetic investigations of the Tussilaginineae have included few, if any, representatives of this Gynoxoid group or reconstructed ambiguous patterns of relationships for it.

METHODS

We sequenced complete plastid genomes of 21 species of the Gynoxoid group and related Tussilaginineae and conducted detailed comparisons of the phylogenetic relationships supported by the gene, intron, and intergenic spacer partitions of these genomes. We also evaluated the impact of manual, motif-based adjustments of automatic DNA sequence alignments on phylogenetic tree inference.

RESULTS

Our results indicate that the inclusion of all plastid genome partitions is needed to infer well-supported phylogenetic trees of the Gynoxoid group. Whole plastome-based tree inference suggests that the genera Gynoxys and Nordenstamia are polyphyletic and form the core clade of the Gynoxoid group. This clade is sister to a clade of Aequatorium and Paragynoxys and also includes some but not all representatives of Paracalia.

CONCLUSIONS

The concatenation and combined analysis of all plastid genome partitions and the construction of manually-curated, motif-based DNA sequence alignments are found to be instrumental in the recovery of well-supported relationships of the Gynoxoid group. We demonstrate that the correct assessment of homology in genome-level plastid sequence data sets is crucial for subsequent phylogeny reconstruction and that the manual post-processing of multiple sequence alignments improves the reliability of such reconstructions amid low genetic distances between taxa.

Complete genus-level plastid phylogenomics of Alismataceae with revisited historical biogeography

Alismataceae, an ancient lineage of monocots, has attracted attention due to its complex evolutionary history, ornamental value, and ecological role. However, the phylogenetic relationships and evolutionary history of the family have not been conclusively resolved. Here, we constructed the first complete genus-level plastid phylogeny of Alismataceae by using 78 genes and updated the historical biogeography based on the phylogenomic tree. Our results divide the Alismataceae into three major clades with robust support values; one clade comprises the former Limnocharitaceae, and the second clade includes the mainly temperate genera Alisma, Baldellia, Damasonium and Luronium, and the monotypic African genus Burnatia as a sister of the temperate genera. The remaining genera are either tropical or have some temperate species in addition to tropical ones, and they constitute the third major clade. Molecular dating and biogeographic analyses suggest that Alismataceae arose in Neotropical, West Palearctic, and Afrotropical regions during the Cretaceous, followed by the split into three main clades due to a combination of vicariance and dispersal events. Unlike earlier studies, we inferred that the mainly temperate clade likely originated from Afrotropical and West Palearctic regions during the Eocene. The most recent common ancestor of the other two clades lived in the Neotropical area during the Late Cretaceous. Long-distance dispersal and vicariance together seem to contribute to the transoceanic distribution of this family.

Species Phylogeny versus Gene Trees: A Case Study of an Incongruent Data Matrix Based on Paphiopedilum Pfitz. (Orchidaceae)

The phylogeny of the genus Paphiopedilum based on the plastome is consistent with morphological analysis. However, to date, none of the analyzed nuclear markers has confirmed this. Topology incongruence among the trees of different nuclear markers concerns entire sections of the subgenus Paphiopedilum. The low-copy nuclear protein-coding gene PHYC was obtained for 22 species representing all sections and subgenera of Paphiopedilum. The nuclear-based phylogeny is supported by morphological characteristics and plastid data analysis. We assumed that an incongruence in nuclear gene trees is caused by ancestral homoploid hybridization. We present a model for inferring the phylogeny of the species despite the incongruence of the different tree topologies. Our analysis, based on six low-copy nuclear genes, is congruent with plastome phylogeny and has been confirmed by phylogenetic network analysis.

Phylotranscriptomic insights into Asteraceae diversity, polyploidy, and morphological innovation

Biodiversity is not evenly distributed among related groups, raising questions about the factors contributing to such disparities. The sunflower family (Asteraceae, >26,000 species) is among the largest and most diverse plant families, but its species diversity is concentrated in a few subfamilies, providing an opportunity to study the factors affecting biodiversity. Phylotranscriptomic analyses here of 244 transcriptomes and genomes produced a phylogeny with strong support for the monophyly of Asteraceae and the monophyly of most subfamilies and tribes. This phylogeny provides a reference for detecting changes in diversification rates and possible factors affecting Asteraceae diversity, which include global climate shifts, whole-genome duplications (WGDs), and morphological evolution. The origin of Asteraceae was estimated at ~83 Mya, with most subfamilies having diverged before the Cretaceous-Paleocene boundary. Phylotranscriptomic analyses supported the existence of 41 WGDs in Asteraceae. Changes to herbaceousness and capitulescence with multiple flower-like capitula, often with distinct florets and scaly pappus/receptacular bracts, are associated with multiple upshifts in diversification rate. WGDs might have contributed to the survival of early Asteraceae by providing new genetic materials to support morphological transitions. The resulting competitive advantage for adapting to different niches would have increased biodiversity in Asteraceae.

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

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

Delimiting the genera of the Ficinia Clade (Cypereae, Cyperaceae) based on molecular phylogenetic data

Generic delimitations in the Ficinia Clade of tribe Cypereae are revisited. In particular, we aim to establish the placement of annual species currently included in Isolepis of which the phylogenetic position is uncertain. Phylogenetic inference is based on two nuclear markers (ETS, ITS) and five plastid markers (the genes matK, ndhF, rbcL and rps16, the trnL intron and trnL-F spacer) data, analyzed using model based methods. Topologies based on nuclear and plastid data show incongruence at the backbone. Therefore, the results are presented separately. The monophyly of the smaller genera (Afroscirpoides, Dracoscirpoides, Erioscirpus, Hellmuthia, Scirpoides) is confirmed. However, Isolepis is paraphyletic as Ficinia is retrieved as one of its clades. Furthermore, Ficinia is paraphyletic if I. marginata and allies are excluded. We take a pragmatic approach based on the nuclear topology, driven by a desire to minimize taxonomic changes, to recircumscribe Ficinia to include the annual Isolepis species characterized by cartilaginous glumes and formally include all the Isolepis species inferred outside the core Isolepis clade. Consequently, the circumscription of Isolepis is narrowed to encompass only those species retrieved as part of the core Isolepis clade. Five new combinations are made (Ficinia neocapensis, Ficinia hemiuncialis, Ficinia incomtula, Ficinia leucoloma, Ficinia minuta). We present nomenclatural summary at genus level, identification keys and diagnostic features.

Phylogenomic analyses based on genome-skimming data reveal cyto-nuclear discordance in the evolutionary history of Cotoneaster (Rosaceae)

As a consequence of hybridization, polyploidization, and apomixis, the genus Cotoneaster (Rosaceae) represents one of the most complicated and controversial lineages in Rosaceae, with ca. 370 species which have been classified into two subgenera and several sections, and is notorious for its taxonomic difficulty. The infrageneric relationships and taxonomy of Cotoneaster have remained poorly understood. Previous studies have focused mainly on natural hybridization involving only several species, and phylogeny based on very limited markers. In the present study, the sequences of complete chloroplast genomes and 204 low-copy nuclear genes of 72 accessions, representing 69 species as ingroups, were used to conduct the most comprehensive phylogenetic analysis so far for Cotoneaster. Based on the sequences of complete chloroplast genomes and many nuclear genes, our analyses yield two robust phylogenetic trees respectively. Chloroplast genome and nuclear data confidently resolved relationships of this genus into two major clades which largely supported current classification based on morphological evidence. However, conflicts between the chloroplast genome and low-copy nuclear phylogenies were observed in both the species level and clade level. Cyto-nuclear discordance in the phylogeny could be caused by frequent hybridization events and incomplete sorting lineage (ILS). In addition, our divergence-time analysis revealed an evolutionary radiation of the genus from late Miocene to date.

Assessing the genetic and chemical diversity of Taraxacum species in the Korean Peninsula

The genetic relationship between Taraxacum species, also known as the dandelion, is complicated because of asexual and mixed sexual apomictic reproduction. The usage of Taraxacum species in traditional medicines make their specialized metabolism important, but interspecific chemical difference has rarely been reported for the genus. In this study, we assembled the chloroplast genome and 45S rDNA of six Taraxacum species that occur in Korea (T. campylodes, T. coreanum, T. erythrospermum, T. mongolicum, T. platycarpum, and T. ussuriense), and performed a comparative analysis, which revealed their phylogenetic relationships and possible natural hybridity. We also performed a liquid chromatography-mass spectrometry-based phytochemical analysis to reveal interspecific chemical diversity. The comparative metabolomics analysis revealed that Taraxacum species could be separated into three chemotypes according to their major defensive specialized metabolites, which were the sesquiterpene lactones, the phenolic inositols, and chlorogenic acid derivatives. The CP DNA- and 45S rDNA-based phylogenetic trees showed a tangled relationship, which supports the notion of ongoing hybridization of wild Taraxacum species. The untargeted LC-MS analysis revealed that each Taraxacum plant exhibits species-specific defensive specialized metabolism. Moreover, 45S rDNA-based phylogenetic tree correlated with the hierarchical cluster relied on metabolite compositions. Given the coincidence between these analyses, we represented that 45S rDNA could well reflect overall nuclear genome variation in Taraxacum species.

Diversification in Qinghai-Tibet Plateau: Orchidinae (Orchidaceae) clades exhibiting pre-adaptations play critical role

We explore the origins of the extraordinary plant diversity in the Qinghai-Tibetan Plateau (QTP) using Orchidinae (Orchidaceae) as a model. Our results indicate that six major clades in Orchidinae exhibited substantial variation in the temporal and spatial sequence of diversification. Our time-calibrated phylogenetic model suggests that the species-richness of Orchidinae arose through a combination of in situ diversification, colonisation, and local recruitment. There are multiple origins of species-richness of Orchidinae in the QTP, and pre-adaptations in clades from North Temperate and alpine regions were crucial for in situ diversification. The geographic analysis identified 29 dispersals from Asia, Africa and Europe into the QTP and 15 dispersals out. Most endemic species of Orchidinae evolved within the past six million years.

Complete Chloroplast Genome Sequencing and Phylogenetic Analysis of Two Dracocephalum Plants

Dracocephalum tanguticum and Dracocephalum moldavica are important herbs from Lamiaceae and have great medicinal value. We used the Illumina sequencing technology to sequence the complete chloroplast genome of D. tanguticum and D. moldavica and then conducted de novo assembly. The two chloroplast genomes have a typical quadripartite structure, with the gene's lengths of 82,221 bp and 81,450 bp, large single-copy region's (LSC) lengths of 82,221 bp and 81,450 bp, and small single-copy region's (SSC) lengths of 17,363 bp and 17,066 bp, inverted repeat region's (IR) lengths of 51,370 bp and 51,352 bp, respectively. The GC content of the two chloroplast genomes was 37.80% and 37.83%, respectively. The chloroplast genomes of the two plants encode 133 and 132 genes, respectively, among which there are 88 and 87 protein-coding genes, respectively, as well as 37 tRNA genes and 8 rRNA genes. Among them, the rps2 gene is unique to D. tanguticum, which is not found in D. moldavica. Through SSR analysis, we also found 6 mutation hotspot regions, which can be used as molecular markers for taxonomic studies. Phylogenetic analysis showed that Dracocephalum was more closely related to Mentha.

Plastomes of eight Ligusticum species: characterization, genome evolution, and phylogenetic relationships

BACKGROUND

The genus Ligusticum consists of approximately 60 species distributed in the Northern Hemisphere. It is one of the most taxonomically difficult taxa within Apiaceae, largely due to the varied morphological characteristics. To investigate the plastome evolution and phylogenetic relationships of Ligusticum, we determined the complete plastome sequences of eight Ligusticum species using a de novo assembly approach.

RESULTS

Through a comprehensive comparative analysis, we found that the eight plastomes were similar in terms of repeat sequence, SSR, codon usage, and RNA editing site. However, compared with the other seven species, L. delavayi exhibited striking differences in genome size, gene number, IR/SC borders, and sequence identity. Most of the genes remained under the purifying selection, whereas four genes showed relaxed selection, namely ccsA, rpoA, ycf1, and ycf2. Non-monophyly of Ligusticum species was inferred from the plastomes and internal transcribed spacer (ITS) sequences phylogenetic analyses.

CONCLUSION

The plastome tree and ITS tree produced incongruent tree topologies, which may be attributed to the hybridization and incomplete lineage sorting. Our study highlighted the advantage of plastome with mass informative sites in resolving phylogenetic relationships. Moreover, combined with the previous studies, we considered that the current taxonomy system of Ligusticum needs to be improved and revised. In summary, our study provides new insights into the plastome evolution, phylogeny, and taxonomy of Ligusticum species.

Genomic evidence for recurrent genetic admixture during the domestication of Mediterranean olive trees (Olea europaea L.)

BACKGROUND

Olive tree (Olea europaea L. subsp. europaea, Oleaceae) has been the most emblematic perennial crop for Mediterranean countries since its domestication around 6000 years ago in the Levant. Two taxonomic varieties are currently recognized: cultivated (var. europaea) and wild (var. sylvestris) trees. However, it remains unclear whether olive cultivars derive from a single initial domestication event followed by secondary diversification, or whether cultivated lineages are the result of more than a single, independent primary domestication event. To shed light into the recent evolution and domestication of the olive tree, here we analyze a group of newly sequenced and available genomes using a phylogenomics and population genomics framework.

RESULTS

We improved the assembly and annotation of the reference genome, newly sequenced the genomes of twelve individuals: ten var. europaea, one var. sylvestris, and one outgroup taxon (subsp. cuspidata)-and assembled a dataset comprising whole genome data from 46 var. europaea and 10 var. sylvestris. Phylogenomic and population structure analyses support a continuous process of olive tree domestication, involving a major domestication event, followed by recurrent independent genetic admixture events with wild populations across the Mediterranean Basin. Cultivated olives exhibit only slightly lower levels of genetic diversity than wild forms, which can be partially explained by the occurrence of a mild population bottleneck 3000-14,000 years ago during the primary domestication period, followed by recurrent introgression from wild populations. Genes associated with stress response and developmental processes were positively selected in cultivars, but we did not find evidence that genes involved in fruit size or oil content were under positive selection. This suggests that complex selective processes other than directional selection of a few genes are in place.

CONCLUSIONS

Altogether, our results suggest that a primary domestication area in the eastern Mediterranean basin was followed by numerous secondary events across most countries of southern Europe and northern Africa, often involving genetic admixture with genetically rich wild populations, particularly from the western Mediterranean Basin.

Diversidad y distribución de los endemismos de Asteraceae (Compositae) en la Flora del Ecuador

La familia Asteraceae está representada en la Flora de Ecuador por un total de 310 táxones endémicos. Están agrupadas en 4 subfamilias (Asteroideae, Barnadesioideae, Cichorioideae y Mutisioideae), 16 tribus y 89 géneros. Doscientas setenta y dos táxones (87,74%) restringen su hábitat a la región andina. Los niveles más altos de endemismos están asociados a los Andes, mayoritariamente desde el bosque andino alto hasta el páramo, y a las Islas Galápagos. Los géneros con mayor riqueza de endemismos son Mikania (26), Pentacalia (23), y Gynoxys (20). Tres géneros monotípicos son endémicos: Cyathomone, Idiopappus y Trigonopterum. Los géneros Darwiniothamnus, Kingianthus, Lecocarpus, y Scalesia también son endémicos. La diversidad de la familia se incrementa desde los 2000 m a 3000 m, alcanzando su mayor riqueza entre los 2900-3000 m, con dominancia de las plantas arbustivas (195 especies, 1 subespecie, 2 variedades) y herbáceas (97 especies). Ciento veinteseis especies tienen categoría de vulnerable, 90 están en peligro de extinción y 24 están en estado crítico de amenaza. Las temperaturas y precipitaciones medias anuales varían significativamente entre los sectores biogeográficos.

Origin and adaptive radiation of the exceptional and threatened bembidiine beetle fauna of St Helena (Coleoptera: Carabidae)

The central peaks of the isolated island of St Helena (south Atlantic Ocean) are home to an extraordinary set of ground beetles of the tribe Bembidiini, which belong to three endemic genus-group taxa. These beetles are strikingly different in overall body form from the many bembidiines found elsewhere in the world. At least some of the St Helena species are likely to be extinct, and all are threatened by habitat destruction and invasive species. Through next-generation sequencing of historical museum specimens, we examine the phylogenetic relationships of the St Helena fauna. We find that, in spite of their morphological disparities, the endemic bembidiines of St Helena form a clade of genetically similar species, with their sister group being Bembidion alsium from the Indian Ocean island of La Réunion, and the sister group of this pair being the African subgenus Omotaphus. We propose that the St Helena Peaks Bembidion are an adaptive radiation that arose from a single dispersal event to St Helena from a now-extinct African lineage (sister to Omotaphus) and that this extinct lineage also served as the ancestral source of B. alsium. Given that the St Helena Peaks Bembidion are deeply nested in the genus Bembidion, we move the three taxa back in that genus as subgenera and provide a new name (Bembidion shepherdae) for the now-homonymous Bembidion wollastoni.

Comparison of chloroplast genomes of Gynura species: sequence variation, genome rearrangement and divergence studies

Background

Some Gynura species have been reported to be natural anti-diabetic plants. Improvement of their traits towards application relies on hybridization. Clearly, phylogenetic relationships could optimize compatible hybridizations. For flowerings plants, chloroplast genomes have been used to solve many phylogenetic relationships. To date, the chloroplast genome sequences of 4 genera of the tribe Senecioneae have been uploaded to GenBank. The internal relationships within the genus Gynura and the relationship of the genus Gynura with other genera in the tribe Senecioneae need further research.

Results

The chloroplast genomes of 4 Gynura species were sequenced, assembled and annotated. In comparison with those of 12 other Senecioneae species, the Gynura chloroplast genome features were analysed in detail. Subsequently, differences in the microsatellite and repeat types in the tribe were found. From the comparison, it was found that IR expansion and contraction are conserved in the genera Gynura, Dendrosenecio and Ligularia. Compared to other regions on the chloroplast genome, the region from 25,000 to 50,000 bp was not conserved. Seven ndh genes in this region are under purifying selection, with small changes in amino acids. The whole chloroplast genome sequences of 16 Senecioneae species were used to build a phylogenetic tree. Based on the oldest Artemisia pollen fossil, the divergence time was estimated.

Conclusions

Sequencing the chloroplast genomes of 4 Gynura species helps us to solve many problems. The phylogenetic relationships and divergence time among 4 Gynura and 16 Senecioneae species were evaluated by comparing their chloroplast genomes. The phylogenetic relationship of the genera Gynura and Ligularia was different from that observed previous work. In a previous phylogenetic tree, the genus Ligularia belonged to the Tussilagininae subtribe, which was in a lineage that diverged earlier than other genera. Further morphology and genome-wide analyses are needed to clarify the genus relationships.

An empirical assessment of a single family-wide hybrid capture locus set at multiple evolutionary timescales in Asteraceae

PREMISE

Hybrid capture with high-throughput sequencing (Hyb-Seq) is a powerful tool for evolutionary studies. The applicability of an Asteraceae family-specific Hyb-Seq probe set and the outcomes of different phylogenetic analyses are investigated here.

METHODS

Hyb-Seq data from 112 Asteraceae samples were organized into groups at different taxonomic levels (tribe, genus, and species). For each group, data sets of non-paralogous loci were built and proportions of parsimony informative characters estimated. The impacts of analyzing alternative data sets, removing long branches, and type of analysis on tree resolution and inferred topologies were investigated in tribe Cichorieae.

RESULTS

Alignments of the Asteraceae family-wide Hyb-Seq locus set were parsimony informative at all taxonomic levels. Levels of resolution and topologies inferred at shallower nodes differed depending on the locus data set and the type of analysis, and were affected by the presence of long branches.

DISCUSSION

The approach used to build a Hyb-Seq locus data set influenced resolution and topologies inferred in phylogenetic analyses. Removal of long branches improved the reliability of topological inferences in maximum likelihood analyses. The Astereaceae Hyb-Seq probe set is applicable at multiple taxonomic depths, which demonstrates that probe sets do not necessarily need to be lineage-specific.

Complete Chloroplast Genome Sequence and Phylogenetic Analysis of Quercus bawanglingensis Huang, Li et Xing, a Vulnerable Oak Tree in China

Quercus bawanglingensis Huang, Li et Xing, an endemic evergreen oak of the genus Quercus (Fagaceae) in China, is currently listed in the Red List of Chinese Plants as a vulnerable (VU) plant. No chloroplast (cp) genome information is currently available for Q. bawanglingensis, which would be essential for the establishment of guidelines for its conservation and breeding. In the present study, the cp genome of Q. bawanglingensis was sequenced and assembled into double-stranded circular DNA with a length of 161,394 bp. Two inverted repeats (IRs) with a total of 51,730 bp were identified, and the rest of the sequence was separated into two single-copy regions, namely, a large single-copy (LSC) region (90,628 bp) and a small single-copy (SSC) region (19,036 bp). The genome of Q. bawanglingensis contains 134 genes (86 protein-coding genes, 40 tRNAs and eight rRNAs). More forward (29) than inverted long repeats (21) are distributed in the cp genome. A simple sequence repeat (SSR) analysis showed that the genome contains 82 SSR loci, involving 84.15% A/T mononucleotides. Sequence comparisons among the nine complete cp genomes, including the genomes of Q. bawanglingensis, Q. tarokoensis Hayata (NC036370), Q. aliena var. acutiserrata Maxim. ex Wenz. (KU240009), Q. baronii Skan (KT963087), Q. aquifolioides Rehd. et Wils. (KX911971), Q. variabilis Bl. (KU240009), Fagus engleriana Seem. (KX852398), Lithocarpus balansae (Drake) A. Camus (KP299291) and Castanea mollissima Bl. (HQ336406), demonstrated that the diversity of SC regions was higher than that of IR regions, which might facilitate identification of the relationships within this extremely complex family. A phylogenetic analysis showed that Fagus engleriana and Trigonobalanus doichangensis form the basis of the produced evolutionary tree. Q. bawanglingensis and Q. tarokoensis, which belong to the group Ilex, share the closest relationship. The analysis of the cp genome of Q. bawanglingensis provides crucial genetic information for further studies of this vulnerable species and the taxonomy, phylogenetics and evolution of Quercus.

Multi-gene incongruence consistent with hybridisation in Cladocopium (Symbiodiniaceae), an ecologically important genus of coral reef symbionts

Coral reefs rely on their intracellular dinoflagellate symbionts (family Symbiodiniaceae) for nutritional provision in nutrient-poor waters, yet this association is threatened by thermally stressful conditions. Despite this, the evolutionary potential of these symbionts remains poorly characterised. In this study, we tested the potential for divergent Symbiodiniaceae types to sexually reproduce (i.e. hybridise) within Cladocopium, the most ecologically prevalent genus in this family. With sequence data from three organelles (cob gene, mitochondrion; psbA^ncr region, chloroplast; and ITS2 region, nucleus), we utilised the Incongruence Length Difference test, Approximately Unbiased test, tree hybridisation analyses and visual inspection of raw data in stepwise fashion to highlight incongruences between organelles, and thus provide evidence of reticulate evolution. Using this approach, we identified three putative hybrid Cladocopium samples among the 158 analysed, at two of the seven sites sampled. These samples were identified as the common Cladocopium types C40 or C1 with respect to the mitochondria and chloroplasts, but the rarer types C3z, C3u and C1# with respect to their nuclear identity. These five Cladocopium types have previously been confirmed as evolutionarily distinct and were also recovered in non-incongruent samples multiple times, which is strongly suggestive that they sexually reproduced to produce the incongruent samples. A concomitant inspection of next generation sequencing data for these samples suggests that other plausible explanations, such as incomplete lineage sorting or the presence of co-dominance, are much less likely. The approach taken in this study allows incongruences between gene regions to be identified with confidence, and brings new light to the evolutionary potential within Symbiodiniaceae.

Unraveling patterns and processes of diversification in the South Andean-Patagonian Nassauvia subgenus Strongyloma (Asteraceae, Nassauvieae)

Congruence among different sources of data is highly desirable in phylogenetic analyses. However, plastid and nuclear DNA may record different evolutionary processes such that incongruence among results from these sources can help unravel complex evolutionary histories. That is the case of Nassauvia subgenus Strongyloma (Asteraceae), a taxon with five putative species distributed in the southern Andes and Patagonian steppe. Morphometric and phylogeographic information cast doubt on the integrity of its species, and previous molecular data even questioned the monophyly of the subgenus. We tested those questions using plastid and nuclear DNA sequences by the application of different methods such as phylogenetic trees, networks, a test of genealogical sorting, an analysis of population structure, calibration of the trees, and hybridization test, assembling non-synchronous incongruent results at subgenus and species levels in a single reconstruction. The integration of our molecular analyses and previous taxonomic, morphological, and molecular studies support subgenus Strongyloma as a monophyletic group. However, the topology of the nuclear trees and the evidence of polyploids within subgenus Nassauvia, suggest a hypothetical origin and initial radiation of Nassauvia related to an ancient hybridization event that occurred around 17-6.3 Myr ago near the Andes in west-central Patagonia. Plastid data suggest a recent diversification within subgenus Strongyloma, at most 9.8 Myr ago, towards the Patagonian steppe east of the Andes. These processes cause phylogenies to deviate from the species tree since each putative species lack exclusive ancestry. The non-monophyly of its species using both plastid and nuclear data is caused mainly by incomplete lineage sorting occurred since the Miocene. The final uplift of the Andes and Pliocene-Pleistocene glacial-interglacial and its consequences on the landscape and climate structured the genetic composition of this group of plants in the Patagonian steppe. The molecular data presented here agree with previous morphological studies, in that the five putative species typically accepted in this subgenus are not independent taxa. This study emphasizes that adding more than one sequence per species, not combining data with dissimilar inheritance patterns without first performed incongruence tests, exploring data through different methodologies, considering the timing of events, and searching for the causes of poorly resolved and/or incongruent phylogenies help to reveal complex biological underlying processes, which might otherwise remain hidden.

Molecular phylogenetics and historical biogeography of the tribe Lilieae (Liliaceae): bi-directional dispersal between biodiversity hotspots in Eurasia

Background and Aims

The role played by the Qinghai-Tibet Plateau (QTP) in the organismal diversification and biogeography of plants in the Northern Hemisphere has attracted much attention from evolutionary biologists. Here we use tribe Lilieae (Liliaceae), including primarily temperate and alpine lineages with disjunct distributions in the North Temperate Zone, as a case study to shed light upon these processes.

Methods

Using 191 taxa (five outgroup taxa) comprising more than 60 % of extant Lilieae species across the entire geographical range, we analyse phylogenetic relationships based on three plastid markers (matK, rbcL, rpl16) and nuclear ITS. Divergence time estimation and ancestral range reconstruction were further inferred.

Key Results

The results support a monophyletic Lilieae divided into four clades. Lilium is nested within Fritillaria, which is paraphyletic and partitioned into two clades, New World and Old World, in the chloroplast DNA (cpDNA) analysis. Incongruences between the ITS and cpDNA trees may be explained by divergent ITS paralogues and hybridization. Lilieae originated around 40-49 (28-67) Mya and probably diversified in the QTP region with four major clades that were established during the Oligocene and the Early Miocene. Uplift of the QTP and climatic changes probably drove early diversification of Lilieae in the QTP region. A rapid radiation occurred during the Late Miocene and the Pleistocene, coinciding temporally with recent orogenic process in the QTP region and climatic oscillations. Several lineages dispersed out of the QTP.

Conclusions

Lineage persistence and explosive radiation were important processes for establishing high species diversity of Lilieae in the QTP region. Both long-distance dispersal and migration across Beringia probably contributed to the modern distribution range of Lilieae. Our study shows that biotic interchanges between the QTP region and Irano-Turanian region and the Mediterranean Basin were bi-directional, suggesting the latter was a secondary centre of diversity.

Revision of Muhlenbergia (Poaceae, Chloridoideae, Cynodonteae, Muhlenbergiinae) in Peru: classification, phylogeny, and a new species, M.romaschenkoi

A taxonomic treatment, phylogeny based on analysis of six DNA sequence markers (ITS, ndhA intron, rpl32-trnL, rps3, rps16 intron and rps16-trnK) and classification of Muhlenbergia for Peru is given. Seventeen species and one presumed hybrid are recognised. Muhlenbergiaromaschenkoi sp. nov. is newly described from the Río Huallaga Valley, northeast of Huánuco. The type of Podosemumangustatum [≡ Muhlenbergiaangustata] clearly aligns with what we had been referring to as the hybrid between this species and M.rigida. Therefore, we adopt the next available heterotypic name, Muhlenbergiacoerulea, for what we had been calling M.angustata and change the hybrid designation to M.coerulea × M.rigida. Lectotypes are designated for Epicampescoerulea Griseb., Muhlenbergiaaffinis Trin., Muhlenbergiaberlandieri Trin., Muhlenbergiabeyrichiana Kunth, Muhlenbergiaelegansvar.atroviolacea Kuntze, Muhlenbergiaelegansvar.subviridis Kuntze and Muhlenbergiaphragmitoides Griseb.