Understanding the spectacular failure of DNA barcoding in willows (Salix): Does this result from a trans-specific selective sweep?

Disentangling serial chloroplast captures in willows

PREMISE

Chloroplast capture is a process through which the chloroplast of a focal species is replaced by the chloroplast from another species during repeated backcrossing of an initial hybrid. Here we investigated serial chloroplast capture from Salix nigra in willows during sequential hybridization events that led to the capture of the same chloroplast lineage across multiple Salix species.

METHODS

Previously generated sequences of nuclear and chloroplast regions from several Salix species were used to identify cases of cytonuclear phylogenetic discordance, a pattern indicating chloroplast capture. Serial chloroplast captures were identified by comparing phylogenetic topologies of the chloroplast trees to discriminate among (1) a single chloroplast capture and subsequent speciation of the lineage with the captured chloroplast, (2) multiple chloroplast captures from the same parent species, and (3) serial chloroplast captures. We also looked for hybridization in genes involved in cytonuclear interactions and in photosynthesis.

RESULTS

We identified cases of serial chloroplast capture and speciation after chloroplast capture in Salix. Although these chloroplast capture events were accompanied by signals of hybridization in the nuclear genomes, nuclear genes that functionally interact with chloroplast genes and nuclear genes involved in photosynthesis were no more likely to introgress in species with chloroplast captures than in species without chloroplast captures.

CONCLUSIONS

This study illuminates the complex evolution of the chloroplast genomes in Salix and the potential for hybridization and introgression to influence genomic evolution.

Backbone phylogeny of Salix based on genome skimming data

The genus Salix is a common component of the Northern Hemisphere dendroflora with important ecological and economic value. However, taxonomy and systematics of Salix is extremely difficult and relationships between main lineages, especially deep phylogenies, remain largely unresolved. In this study, we used genome-skimming, plastome assembly, and single-copy orthologs (SCOs) from 66 Salix accessions, along with publicly available plastome and sequence read archive (SRA) datasets to obtain a robust backbone phylogeny of Salix, clarify relationships between its main lineages, and gain a more precise understanding of the origin and diversification of this species-rich genus. The plastome and SCO datasets resolved Salix into two robust clades, with plastome-based phylogenies lacking inner resolution and SCO offering fully resolved phylogenies. Our results support the classification of Salix into five subgenera: Salix, Urbaniana, Triandrae, Longifoliae and Vetrix. We observed a significant acceleration in the diversification rate within the Chamaetia-Vetrix clade, while Salix exhibited increased rates of diversification spanning from the early Oligocene to the late Miocene. These changes coincided with contemporaneous tectonic and climate change events. Our results provide a foundation for future systematic and evolutionary studies of Salix. Additionally, we showed that genome skimming data is an efficient, rapid, and reliable approach for obtaining extensive genomic data for phylogenomic studies, enabling the comprehensive elucidation of Salix relationships.

Impacts of Climate Change Conditions on the Potential Distribution of Anoplophora glabripennis and Its Host Plants, Salix babylonica and Salix matsudana, in China

The Anoplophora glabripennis (Motschulsky) is a phytophagous pest that is seriously endangering Salix babylonica Linn. and S. matsudana Koidz. Poor control can damage local ecosystems, resulting in economic losses and management risks. In the context of climate change, the climatic ecological niche of organisms is no longer compatible with the surrounding environment. To mitigate the effects of climate change, some organisms respond adaptively to climate change through different mechanisms and in different ways. In this study, an optimized MaxEnt model was used to explore the potential distribution areas of A. glabripennis and its host plants, S. babylonica and S. matsudana, in response to current and future climate and to determine their movement routes and relative dynamics. The results show that the optimized model exhibits the lowest complexity and excellent prediction accuracy. It is important to note that both temperature and precipitation are the main environmental factors affecting the distribution of suitable habitats for A. glabripennis and its host plants. This is evidenced by the mean temperature of the warmest quarter and precipitation of the wettest month being the main environmental factors affecting the distribution of suitable habitats for the host plants. Similarly, the minimum temperature of the coldest month and precipitation seasonality are the primary bioclimatic variables constraining the dispersal of A. glabripennis. Under climate change, the suitable areas of both S. babylonica and S. matsudana are declining, while the suitable areas of A. glabripennis are expanding in future climates. Furthermore, three species exhibited a proclivity for migration to higher latitudes in response to climate change. In conclusion, this study contributes to our understanding of the biogeographic characteristics of these A. glabripennis, S. babylonica, and S. matsudana and provides a basis for the formulation of timely conservation strategies to reduce the potential impacts of climate change. This is of great significance for the rational management, utilization, and protection of forest ecosystems in China.

Characterization of the complete chloroplast genome and development of molecular markers of Salix

Salix, an economically and ecologically multifunctional tree species widely distributed in China, encompasses five ornamental species sequenced in this study, which are highly beneficial for plant phytoremediation due to their ability to absorb heavy metals. This research utilized high-throughput sequencing to acquire chloroplast genome sequences of Salix, analyzing their gene composition and structural characteristics, identifying potential molecular markers, and laying a foundation for Salix identification and resource classification. Chloroplast DNA was extracted from the leaves of Salix argyracea, Salix dasyclados, Salix eriocephala, Salix integra 'Hakuro Nishiki', and Salix suchowensis using an optimized CTAB method. Sequencing was conducted on the Illumina NovaSeq PE150 platform, and bioinformatics tools were employed to compare the structural features and variations within the chloroplast genomes of the Salix. Analysis revealed high similarity among the chloroplast genome sequences of the five Salix species, with a subsequent examination identifying 276, 269, 270, 273, and 273 SSR loci, respectively, along with unique simple repeat sequences in each variety. Comparison of chloroplast genomes across 22 Salix highlighted variations in regions such as matK-trnQ, ndhC-trnV, psbE-petL, rpl36-rps8, and ndhB-rps7, which may serve as valuable molecular markers for willow resource classification studies. In this study, chloroplast genome sequencing and structural analysis of Salix not only enhances the genetic resources of Salix but also forms a critical basis for the development of molecular markers and the exploration of interspecific phylogeny in the genus.

Skmer approach improves species discrimination in taxonomically problematic genus Schima (Theaceae)

Genome skimming has dramatically extended DNA barcoding from short DNA fragments to next generation barcodes in plants. However, conserved DNA barcoding markers, including complete plastid genome and nuclear ribosomal DNA (nrDNA) sequences, are inadequate for accurate species identification. Skmer, a recently proposed approach that estimates genetic distances among species based on unassembled genome skims, has been proposed to effectively improve species discrimination rate. In this study, we used Skmer to identify species based on genomic skims of 47 individuals representing 10 out of 13 species of Schima (Theaceae) from China. The unassembled reads identified six species, with a species identification rate of 60%, twice as high as previous efforts that used plastid genomes (27.27%). In addition, Skmer was able to identify Schima species with only 0.5× sequencing depth, as six species were well-supported with unassembled data sizes as small as 0.5 Gb. These findings demonstrate the potential for Skmer approach in species identification, where nuclear genomic data plays a crucial role. For taxonomically difficult taxa such as Schima, which have diverged recently and have low levels of genetic variation, Skmer is a promising alternative to next generation barcodes.

Phylogenetic insights into the Salicaceae: The evolution of willows and beyond

The Salicaceae includes approximately 54 genera and over 1,400 species with a cosmopolitan distribution. Members of the family are well-known for their diverse secondary plant metabolites, and they play crucial roles in tropical and temperate forest ecosystems. Phylogenetic reconstruction of the Salicaceae has been historically challenging due to the limitations of molecular markers and the extensive history of hybridization and polyploidy within the family. Our study employs whole-genome sequencing of 74 species to generate an extensive phylogeny of the Salicaceae. We generated two RAD-Seq enriched whole-genome sequence datasets and extracted two additional gene sets corresponding to the universal Angiosperms353 and Salicaceae-specific targeted-capture arrays. We reconstructed maximum likelihood-based molecular phylogenies using supermatrix and coalescent-based supertree approaches. Our fossil-calibrated phylogeny estimates that the Salicaceae originated around 128 million years ago and unravels the complex taxonomic relationships within the family. Our findings confirm the non-monophyly of the subgenus Salix s.l. and further support the merging of subgenera Chamaetia and Vetrix, both of which exhibit intricate patterns within and among different sections. Overall, our study not only enhances our understanding of the evolution of the Salicaceae, but also provides valuable insights into the complex relationships within the family.

Sex chromosome turnover plays an important role in the maintenance of barriers to post-speciation introgression in willows

Almost all species in the genus Salix (willow) are dioecious and willows have variable sex-determining systems, the role of this variation in maintaining species barriers is relatively untested. We first analyzed the sex determination systems (SDS) of two species, Salix cardiophylla and Salix interior, whose positions in the Salix phylogeny make them important for understanding a sex chromosome turnover that has been detected in their relatives, and that changed the system from male (XX/XY) to female (ZW/ZZ) heterogamety. We show that both species have male heterogamety, with sex-linked regions (SLRs) on chromosome 15 (termed a 15XY system). The SLRs occupy 21.3% and 22.8% of the entire reference chromosome, respectively. By constructing phylogenetic trees, we determined the phylogenetic positions of all the species with known SDSs. Reconstruction of ancestral SDS character states revealed that the 15XY system is likely the ancestral state in willows. Turnovers of 15XY to 15ZW and 15XY to 7XY likely contributed to early speciation in Salix and gave rise to major groups of the Vetrix and Salix clades. Finally, we tested introgression among species in the phylogenetic trees based on both autosomes and SLRs separately. Frequent introgression was observed among species with 15XY, 15ZW, and 7XY on autosomes, in contrast to the SLR datasets, which showed less introgression, and in particular no gene flow between 15ZW and 7XY species. We argue that, although SDS turnovers in willow speciation may not create complete reproductive barriers, the evolution of SLRs plays important roles in preventing introgression and maintaining species boundaries.

Challenge accepted: Evolutionary lineages versus taxonomic classification of North American shrub willows (Salix)

PREMISE

The huge diversity of Salix subgenus Chamaetia/Vetrix clade in North America and the lack of phylogenetic resolution within this clade has presented a difficult but fascinating challenge for taxonomists to resolve. Here we tested the existing taxonomic classification with molecular tools.

METHODS

In this study, 132 samples representing 46 species from 22 described sections of shrub willows from the United States and Canada were analyzed and combined with 67 samples from Eurasia. The ploidy levels of the samples were determined using flow cytometry and nQuire. Sequences were produced using a RAD sequencing approach and subsequently analyzed with ipyrad, then used for phylogenetic reconstructions (RAxML, SplitsTree), dating analyses (BEAST, SNAPPER), and character evolution analyses of 14 selected morphological traits (Mesquite).

RESULTS

The RAD sequencing approach allowed the production of a well-resolved phylogeny of shrub willows. The resulting tree showed an exclusively North American (NA) clade in sister position to a Eurasian clade, which included some North American endemics. The NA clade began to diversify in the Miocene. Polyploid species appeared in each observed clade. Character evolution analyses revealed that adaptive traits such as habit and adaxial nectaries evolved multiple times independently.

CONCLUSIONS

The diversity in shrub willows was shaped by an evolutionary radiation in North America. Most species were monophyletic, but the existing sectional classification could not be supported by molecular data. Nevertheless, monophyletic lineages share several morphological characters, which might be useful in the revision of the taxonomic classification of shrub willows.

A cryptic syngameon within Betula shrubs revealed: Implications for conservation in changing subarctic environments

Arctic and subarctic ecosystems are rapidly transforming due to global warming, emphasizing the need to understand the genetic diversity and adaptive strategies of northern plant species for effective conservation. This study focuses on Betula glandulosa, a native North American tundra shrub known as dwarf birch, which demonstrates an apparent capacity to adapt to changing climate conditions. To address the taxonomic challenges associated with shrub birches and logistical complexities of sampling in the northernmost areas where species' ranges overlap, we adopted a multicriteria approach. Incorporating molecular data, ploidy level assessment and leaf morphology, we aimed to distinguish B. glandulosa individuals from other shrub birch species sampled. Our results revealed three distinct species and their hybrids within the 537 collected samples, suggesting the existence of a shrub birch syngameon, a reproductive network of interconnected species. Additionally, we identified two discrete genetic clusters within the core species, B. glandulosa, that likely correspond to two different glacial lineages. A comparison between the nuclear and chloroplast SNP data emphasizes a long history of gene exchange between different birch species and genetic clusters. Furthermore, our results highlight the significance of incorporating interfertile congeneric species in conservation strategies and underscores the need for a holistic approach to conservation in the context of climate change, considering the complex dynamics of species interactions. While further research will be needed to describe this shrub birches syngameon and its constituents, this study is a first step in recognizing its existence and disseminating awareness among ecologists and conservation practitioners. This biological phenomenon, which offers evolutionary flexibility and resilience beyond what its constituent species can achieve individually, may have significant ecological implications.

A Geometric Morphometrics Approach to the Study of Natural Variations and Hybrid Detection in Populations of Alnus incana (L.) Moench and Alnus rohlenae Vít, Douda and Mandák

Landmark-based geometric morphometrics (GM) was used to examine, for the first time, spontaneous hybridization between Alnus incana (L.) Moench and Alnus rohlenae Vít, Douda and Mandák, and to assess inter- and intrapopulation variability in leaf shape, leaf size and venation in natural populations in Serbia (Western Balkans). Two geographically distant (30 km) and two close (1.2 km) populations were selected to examine hybridization. The variability in leaf shapes was assessed by canonical variate analysis and linear discriminant analysis performed on the symmetric component of variation. Covariation between the symmetric component of shape variation and the number of pairs of secondary leaf veins was investigated with partial least squares analysis. Static allometry was examined for the first time in the genus Alnus Mill. A higher proportion of A. incana leaves was classified as A. rohlenae in geographically close populations, which is in accordance with the hypothesis about spontaneous hybridization. No single leaf of A. rohlenae was classified as A. incana, indicating that putative hybrids can only be found in grey alder populations. This study demonstrates that GM is a powerful tool for species delimitation and hybrid detection in the genus Alnus and it can be used for preliminary screening in hybrid zones.

Phylogenomics reveals patterns of ancient hybridization and differential diversification that contribute to phylogenetic conflict in willows, poplars, and close relatives

Despite the economic, ecological, and scientific importance of the genera Salix L. (willows) and Populus L. (poplars, cottonwoods, and aspens) Salicaceae, we know little about the sources of differences in species diversity between the genera and of the phylogenetic conflict that often confounds estimating phylogenetic trees. Salix subgenera and sections, in particular, have been difficult to classify, with one recent attempt termed a "spectacular failure" due to a speculated radiation of the subgenera Vetrix and Chamaetia. Here, we use targeted sequence capture to understand the evolutionary history of this portion of the Salicaceae plant family. Our phylogenetic hypothesis was based on 787 gene regions and identified extensive phylogenetic conflict among genes. Our analysis supported some previously described subgeneric relationships and confirmed the polyphyly of others. Using an fbranch analysis, we identified several cases of hybridization in deep branches of the phylogeny, which likely contributed to discordance among gene trees. In addition, we identified a rapid increase in diversification rate near the origination of the Vetrix-Chamaetia clade in Salix. This region of the tree coincided with several nodes that lacked strong statistical support, indicating a possible increase in incomplete lineage sorting due to rapid diversification. The extraordinary level of both recent and ancient hybridization in both Salix and Populus have played important roles in the diversification and diversity in these two genera.

Plastome phylogenomics and morphological traits analyses provide new insights into the phylogenetic position, species delimitation and speciation of Triplostegia (Caprifoliaceae)

BACKGROUND

The genus Triplostegia contains two recognized species, T. glandulifera and T. grandiflora, but its phylogenetic position and species delimitation remain controversial. In this study, we assembled plastid genomes and nuclear ribosomal DNA (nrDNA) cistrons sampled from 22 wild Triplostegia individuals, each from a separate population, and examined these with 11 recently published Triplostegia plastomes. Morphological traits were measured from herbarium specimens and wild material, and ecological niche models were constructed.

RESULTS

Triplostegia is a monophyletic genus within the subfamily Dipsacoideae comprising three monophyletic species, T. glandulifera, T. grandiflora, and an unrecognized species Triplostegia sp. A, which occupies much higher altitude than the other two. The new species had previously been misidentified as T. glandulifera, but differs in taproot, leaf, and other characters. Triplotegia is an old genus, with stem age 39.96 Ma, and within it T. glandulifera diverged 7.94 Ma. Triplostegia grandiflora and sp. A diverged 1.05 Ma, perhaps in response to Quaternary climate fluctuations. Niche overlap between Triplostegia species was positively correlated with their phylogenetic relatedness.

CONCLUSIONS

Our results provide new insights into the species delimitation of Triplostegia, and indicate that a taxonomic revision of Triplostegia is needed. We also identified that either rpoB-trnC or ycf1 could serve as a DNA barcode for Triplostegia.

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.

Genomic evidence for contrasting patterns of host-associated genetic differentiation across shared host-plant species in leaf- and bud-galling sawflies

Resource specialization and ecological speciation arising through host-associated genetic differentiation (HAD) are frequently invoked as an explanation for the high diversity of plant-feeding insects and other organisms with a parasitic lifestyle. While genetic studies have demonstrated numerous examples of HAD in insect herbivores, the rarity of comparative studies means that we still lack an understanding of how deterministic HAD is, and whether patterns of host shifts can be predicted over evolutionary timescales. We applied genome-wide single nucleotide polymorphism and mitochondrial DNA sequence data obtained through genome resequencing to define species limits and to compare host-plant use in population samples of leaf- and bud-galling sawflies (Hymenoptera: Tenthredinidae: Nematinae) collected from seven shared willow (Salicaceae: Salix) host species. To infer the repeatability of long-term cophylogenetic patterns, we also contrasted the phylogenies of the two galler groups with each other as well as with the phylogeny of their Salix hosts estimated based on RADseq data. We found clear evidence for host specialization and HAD in both of the focal galler groups, but also that leaf gallers are more specialized to single host species compared with most bud gallers. In contrast to bud gallers, leaf gallers also exhibited statistically significant cophylogenetic signal with their Salix hosts. The observed discordant patterns of resource specialization and host shifts in two related galler groups that have radiated in parallel across a shared resource base indicate a lack of evolutionary repeatability in the focal system, and suggest that short- and long-term host use and ecological diversification in plant-feeding insects are dominated by stochasticity and/or lineage-specific effects.

Evolution of a hybrid zone of two willow species (Salix L.) in the European Alps analyzed by RAD-seq and morphometrics

Natural hybridization of plants can result in many outcomes with several evolutionary consequences, such as hybrid speciation and introgression. Natural hybrid zones can arise in mountain systems as a result of fluctuating climate during the exchange of glacial and interglacial periods, where species retract and expand their territories, resulting in secondary contacts. Willows are a large genus of woody plants with an immense capability of interspecific crossing. In this study, the sympatric area of two diploid sister species, S. foetida and S. waldsteiniana in the eastern European Alps, was investigated to study the genomic structure of populations within and outside their contact zone and to analyze congruence of morphological phenotypes with genetic data. Eleven populations of the two species were sampled across the Alps and examined using phylogenetic network and population genetic structure analyses of RAD Seq data and morphometric analyses of leaves. The results showed that a homoploid hybrid zone between the two species was established within their sympatric area. Patterns of genetic admixture in homoploid hybrids indicated introgression with asymmetric backcrossing to not only one of the parental species but also one hybrid population forming a separate lineage. The lack of F1 hybrids indicated a long-term persistence of the hybrid populations. Insignificant isolation by distance suggests that gene flow can act over large geographical scales. Morphometric characteristics of hybrids supported the molecular data and clearly separated populations of the parental species, but showed intermediacy in the hybrid zone populations with a bias toward S. waldsteiniana. The homoploid hybrid zone might have been established via secondary contact hybridization, and its establishment was fostered by the low genetic divergence of parental species and a lack of strong intrinsic crossing barriers. Incomplete ecological separation and the ability of long-distance dispersal of willows could have contributed to the spatial expansion of the hybrid zone.

Super-Mitobarcoding in Plant Species Identification? It Can Work! The Case of Leafy Liverworts Belonging to the Genus Calypogeia

Molecular identification of species is especially important where traditional taxonomic methods fail. The genus Calypogeia belongs to one of the tricky taxons. The simple morphology of these species and a tendency towards environmental plasticity make them complicated in identification. The finding of the universal single-locus DNA barcode in plants seems to be 'the Holy Grail'; therefore, researchers are increasingly looking for multiloci DNA barcodes or super-barcoding. Since the mitochondrial genome has low sequence variation in plants, species delimitation is usually based on the chloroplast genome. Unexpectedly, our research shows that super-mitobarcoding can also work! However, our outcomes showed that a single method of molecular species delimitation should be avoided. Moreover, it is recommended to interpret the results of molecular species delimitation alongside other types of evidence, such as ecology, population genetics or comparative morphology. Here, we also presented genetic data supporting the view that C. suecica is not a homogeneous species.

Hybridisation and chloroplast capture between distinct Themeda triandra lineages in Australia

Ecotypes are distinct populations within a species that are adapted to specific environmental conditions. Understanding how these ecotypes become established, and how they interact when reunited, is fundamental to elucidating how ecological adaptations are maintained. This study focuses on Themeda triandra, a dominant grassland species across Asia, Africa and Australia. It is the most widespread plant in Australia, where it has distinct ecotypes that are usually restricted to either wetter and cooler coastal regions or the drier and hotter interior. We generate a reference genome for T. triandra and use whole genome sequencing for over 80 Themeda accessions to reconstruct the evolutionary history of T. triandra and related taxa. Organelle phylogenies confirm that Australia was colonized by T. triandra twice, with the division between ecotypes predating their arrival in Australia. The nuclear genome provides evidence of differences in the dominant ploidal level and gene-flow among the ecotypes. In northern Queensland there appears to be a hybrid zone between ecotypes with admixed nuclear genomes and shared chloroplast haplotypes. Conversely, in the cracking claypans of Western Australia, there is cytonuclear discordance with individuals possessing the coastal chloroplast and interior clade nuclear genome. This chloroplast capture is potentially a result of adaptive introgression, with selection detected in the rpoC2 gene which is associated with water use efficiency. The reason that T. triandra is the most widespread plant in Australia appears to be a result of distinct ecotypic genetic variation and genome duplication, with the importance of each depending on the geographic scale considered.

Testing the effectiveness of rbcLa DNA-barcoding for species discrimination in tropical montane cloud forest vascular plants (Oaxaca, Mexico) using BLAST, genetic distance, and tree-based methods

DNA-barcoding is a species identification tool that uses a short section of the genome that provides a genetic signature of the species. The main advantage of this novel technique is that it requires a small sample of tissue from the tested organism. In most animal groups, this technique is very effective. However, in plants, the recommended standard markers, such as rbcLa, may not always work, and their efficacy remains to be tested in many plant groups, particularly from the Neotropical region. We examined the discriminating power of rbcLa in 55 tropical cloud forest vascular plant species from 38 families (Oaxaca, Mexico). We followed the CBOL criteria using BLASTn, genetic distance, and monophyly tree-based analyses (neighbor-joining, NJ, maximum likelihood, ML, and Bayesian inference, BI). rbcLa universal primers amplified 69.0% of the samples and yielded 91.3% bi-directional sequences. Sixty-three new rbcLa sequences were established. BLAST discriminates 80.8% of the genus but only 15.4% of the species. There was nil minimum interspecific genetic distances in Quercus, Oreopanax, and Daphnopsis. Contrastingly, Ericaceae (5.6%), Euphorbiaceae (4.6%), and Asteraceae (3.3%) species displayed the highest within-family genetic distances. According to the most recent angiosperm classification, NJ and ML trees successfully resolved (100%) monophyletic species. ML trees showed the highest mean branch support value (87.3%). Only NJ and ML trees could successfully discriminate Quercus species belonging to different subsections: Quercus martinezii (white oaks) from Q. callophylla and Q. laurina (red oaks). The ML topology could distinguish species in the Solanaceae clade with similar BLAST matches. Also, the BI topology showed a polytomy in this clade, and the NJ tree displayed low-support values. We do not recommend genetic-distance approaches for species discrimination. Severe shortages of rbcLa sequences in public databases of neotropical species hindered effective BLAST comparisons. Instead, ML tree-based analysis displays the highest species discrimination among the tree-based analyses. With the ML topology in selected genera, rbcLa helped distinguish infra-generic taxonomic categories, such as subsections, grouping affine species within the same genus, and discriminating species. Since the ML phylogenetic tree could discriminate 48 species out of our 55 studied species, we recommend this approach to resolve tropical montane cloud forest species using rbcLa, as an initial step and improve DNA amplification methods.

Measuring the Invisible: The Sequences Causal of Genome Size Differences in Eyebrights (Euphrasia) Revealed by k-mers

Genome size variation within plant taxa is due to presence/absence variation, which may affect low-copy sequences or genomic repeats of various frequency classes. However, identifying the sequences underpinning genome size variation is challenging because genome assemblies commonly contain collapsed representations of repetitive sequences and because genome skimming studies by design miss low-copy number sequences. Here, we take a novel approach based on k-mers, short sub-sequences of equal length k, generated from whole-genome sequencing data of diploid eyebrights (Euphrasia), a group of plants that have considerable genome size variation within a ploidy level. We compare k-mer inventories within and between closely related species, and quantify the contribution of different copy number classes to genome size differences. We further match high-copy number k-mers to specific repeat types as retrieved from the RepeatExplorer2 pipeline. We find genome size differences of up to 230Mbp, equivalent to more than 20% genome size variation. The largest contributions to these differences come from rDNA sequences, a 145-nt genomic satellite and a repeat associated with an Angela transposable element. We also find size differences in the low-copy number class (copy number ≤ 10×) of up to 27 Mbp, possibly indicating differences in gene space between our samples. We demonstrate that it is possible to pinpoint the sequences causing genome size variation within species without the use of a reference genome. Such sequences can serve as targets for future cytogenetic studies. We also show that studies of genome size variation should go beyond repeats if they aim to characterise the full range of genomic variants. To allow future work with other taxonomic groups, we share our k-mer analysis pipeline, which is straightforward to run, relying largely on standard GNU command line tools.

Species discrimination in Schima (Theaceae): Next-generation super-barcodes meet evolutionary complexity

Plastid genome and nrDNA arrays, proposed recently as "super barcodes", might provide additional discriminatory power and overcome the limitations of traditional barcoding loci, yet super barcodes need to be tested for their effectiveness in more plant groups. Morphological homoplasy among Schima species makes the genus a model for testing the efficacy of super barcodes. In this study, we generated multiple datasets comprising standard DNA barcodes (matK, rbcL, trnH-psbA, nrITS) and super-barcodes (plastid genome, nrDNA arrays) across 58 individuals from 12 out of 13 species of Schima from China. No samples were correctly assigned to species using standard DNA barcodes and nrDNA arrays, while only 27.27% of species with multiple accessions were distinguished using the plastid genome and its partitioned datasets-the lowest estimated rate of super barcode success in the literature so far. For Schima and other taxa with similarly recently divergence and low levels of genetic variation, incomplete lineage sorting, hybridization, or taxonomic oversplitting are all possible causes of the failure. Taken together, our study suggests that by no means are super barcodes immune to the challenges imposed by evolutionary complexity. We therefore call for developing multi-locus nuclear markers for species discrimination in plant groups.

Multilocus marker-based delimitation of Salicornia persica and its population discrimination assisted by supervised machine learning approach

The Salicornia L. has been considered one of the most taxonomically challenging genera due to high morphological plasticity, intergradation between related species, and lack of diagnostic features in preserved herbarium specimens. In the United Arab Emirates (UAE), only one species of this genus, Salicornia europaea, has been reported, though investigating its identity at the molecular level has not yet been undertaken. Moreover, based on growth form and morphology variation between the Ras-Al-Khaimah (RAK) population and the Umm-Al-Quwain (UAQ) population, we suspect the presence of different species or morphotypes. The present study aimed to initially perform species identification using multilocus DNA barcode markers from chloroplast DNA (cpDNA) and nuclear ribosomal DNA (nrDNA), followed by the genetic divergence between two populations (RAK and UAQ) belonging to two different coastal localities in the UAE. The analysis resulted in high-quality multilocus barcode sequences subjected to species discrimination through the unsupervised OTU picking and supervised learning methods. The ETS sequence data from our study sites had high identity with the previously reported sequences of Salicornia persica using NCBI blast and was further confirmed using OTU picking methods viz., TaxonDNAs Species identifier and Assemble Species by Automatic Partitioning (ASAP). Moreover, matK sequence data showed a non-monophyletic relationship, and significant discrimination between the two populations through alignment-based unsupervised OTU picking, alignment-free Co-Phylog, and alignment & alignment-free supervised learning approaches. Other markers viz., rbcL, trnH-psbA, ITS2, and ETS could not distinguish the two populations individually, though their combination with matK (cpDNA & cpDNA+nrDNA) showed enough population discrimination. However, the ITS2+ETS (nrDNA) exhibited much higher genetic divergence, further splitting both the populations into four haplotypes. Based on the observed morphology, genetic divergence, and the number of haplotypes predicted using the matK marker, it can be suggested that two distinct populations (RAK and UAQ) do exist. Further extensive morpho-taxonomic studies are required to determine the inter-population variability of Salicornia in the UAE. Altogether, our results suggest that S. persica is the species that grow in the present study area in UAE, and do not support previous treatments as S. europaea.

Testing the Complete Plastome for Species Discrimination, Cryptic Species Discovery and Phylogenetic Resolution in Cephalotaxus (Cephalotaxaceae)

Species of Cephalotaxus have great economic and ecological values. However, the taxonomy and interspecific phylogenetic relationships within the genus have been controversial and remained not fully resolved until now. To date, no study examined the efficiency of the complete plastome as super-barcode across Cephalotaxus species with multiple samples per taxon. In this study, we have evaluated the complete plastome in species discrimination and phylogenetic resolution in Cephalotaxus by including 32 individuals of all eight recognized species and five varieties following Farjon's classification (2010) with multiple samples per taxon. Our results indicated that not all species recognized in recent taxonomic revisions of Cephalotaxus could be distinguished and not all were monophyletic. Based on the plastome phylogeny, a new taxonomic classification for the genus comprising nine species and two varieties, including a cryptic species, was proposed. The phylogeny also resolved all interspecific relationships. Compared to the plastome based classification, standard DNA barcodes, alone or in combination, only recognized a maximum of seven out of the nine species. Moreover, two highly variable single loci, ycf1 and rps16, each alone achieved full species discrimination. With the moderate length of 1079 bp, rps16 is proposed as a specific barcode to discriminate Cephalotaxus species. The super-barcodes and specific barcode candidates will aid in the identification of endangered Cephalotaxus species, and to help focus conservation measures.

Highly Diverse Shrub Willows (Salix L.) Share Highly Similar Plastomes

Plastome phylogenomics is used in a broad range of studies where single markers do not bear enough information. Phylogenetic reconstruction in the genus Salix is difficult due to the lack of informative characters and reticulate evolution. Here, we use a genome skimming approach to reconstruct 41 complete plastomes of 32 Eurasian and North American Salix species representing different lineages, different ploidy levels, and separate geographic regions. We combined our plastomes with published data from Genbank to build a comprehensive phylogeny of 61 samples (50 species) using RAxML (Randomized Axelerated Maximum Likelihood). Additionally, haplotype networks for two observed subclades were calculated, and 72 genes were tested to be under selection. The results revealed a highly conserved structure of the observed plastomes. Within the genus, we observed a variation of 1.68%, most of which separated subg. Salix from the subgeneric Chamaetia/Vetrix clade. Our data generally confirm previous plastid phylogenies, however, within Chamaetia/Vetrix phylogenetic results represented neither taxonomical classifications nor geographical regions. Non-coding DNA regions were responsible for most of the observed variation within subclades and 5.6% of the analyzed genes showed signals of diversifying selection. A comparison of nuclear restriction site associated DNA (RAD) sequencing and plastome data on a subset of 10 species showed discrepancies in topology and resolution. We assume that a combination of (i) a very low mutation rate due to efficient mechanisms preventing mutagenesis, (ii) reticulate evolution, including ancient and ongoing hybridization, and (iii) homoplasy has shaped plastome evolution in willows.

Testing genome skimming for species discrimination in the large and taxonomically difficult genus Rhododendron

Standard plant DNA barcodes based on 2-3 plastid regions, and nrDNA ITS show variable levels of resolution, and fail to discriminate among species in many plant groups. Genome skimming to recover complete plastid genome sequences and nrDNA arrays has been proposed as a solution to address these resolution limitations. However, few studies have empirically tested what gains are achieved in practice. Of particular interest is whether adding substantially more plastid and nrDNA characters will lead to an increase in discriminatory power, or whether the resolution limitations of standard plant barcodes are fundamentally due to plastid genomes and nrDNA not tracking species boundaries. To address this, we used genome skimming to recover near-complete plastid genomes and nuclear ribosomal DNA from Rhododendron species and compared discrimination success with standard plant barcodes. We sampled 218 individuals representing 145 species of this species-rich and taxonomically difficult genus, focusing on the global biodiversity hotspots of the Himalaya-Hengduan Mountains. Only 33% of species were distinguished using ITS+matK+rbcL+trnH-psbA. In contrast, 55% of species were distinguished using plastid genome and nrDNA sequences. The vast majority of this increase is due to the additional plastid characters. Thus, despite previous studies showing an asymptote in discrimination success beyond 3-4 plastid regions, these results show that a demonstrable increase in discriminatory power is possible with extensive plastid genome data. However, despite these gains, many species remain unresolved, and these results also reinforce the need to access multiple unlinked nuclear loci to obtain transformative gains in species discrimination in plants.

Plant sex chromosomes defy evolutionary models of expanding recombination suppression and genetic degeneration

Hundreds of land plant lineages have independently evolved separate sexes in either gametophytes (dioicy) or sporophytes (dioecy), but 43% of all dioecious angiosperms are found in just 34 entirely dioecious clades, suggesting that their mode of sex determination evolved a long time ago. Here, we review recent insights on the molecular mechanisms that underlie the evolutionary change from individuals that each produce male and female gametes to individuals specializing in the production of just one type of gamete. The canonical model of sex chromosome evolution in plants predicts that two sex-determining genes will become linked in a sex-determining region (SDR), followed by expanding recombination suppression, chromosome differentiation and, ultimately, degeneration. Experimental work, however, is showing that single genes function as master regulators in model systems, such as the liverwort Marchantia and the angiosperms Diospyros and Populus. In Populus, this type of regulatory function has been demonstrated by genome editing. In other systems, including Actinidia, Asparagus and Vitis, two coinherited factors appear to independently regulate female and male function, yet sex chromosome differentiation has remained low. We discuss the best-understood systems and evolutionary pathways to dioecy, and present a meta-analysis of the sizes and ages of SDRs. We propose that limited sexual conflict explains why most SDRs are small and sex chromosomes remain homomorphic. It appears that models of increasing recombination suppression with age do not apply because selection favours mechanisms in which sex determination depends on minimal differences, keeping it surgically precise.

HPLC-UV, Metabarcoding and Genome Skims of Botanical Dietary Supplements: A Case Study in Echinacea

The use of DNA-based methods to authenticate botanical dietary supplements has been vigorously debated for a variety of reasons. More comparisons of DNA-based and chemical methods are needed, and concordant evaluation of orthogonal approaches on the same products will provide data to better understand the strengths and weaknesses of both approaches. The overall application of DNA-based methods is already firmly integrated into a wide array of continually modernizing stand alone and complementary authentication protocols. Recently, the use of full-length chloroplast genome sequences provided enhanced discriminatory capacity for closely related species of Echinacea compared to traditional DNA barcoding approaches (matK and rbcL). Here, two next-generation sequencing approaches were used: (1) genome skimming and (2) PCR amplicon (metabarcoding). The two genetic approaches were then combined with HPLC-UV to evaluate 20 commercially available dietary supplements of Echinacea representing "finished" products. The trade-offs involved in different DNA approaches were discussed, with a focus on how DNA methods support existing, accepted chemical methods. In most of the products (19/20), HPLC-UV suggested the presence of Echinacea spp. While metabarcoding was not useful with this genus and instead only resolved 7 products to the family level, genome skimming was able to resolve to species (9) or genus (1) with the 10/20 products where it was successful. Additional ingredients that HPLC-UV was unable to identify were also found in four products along with the relative sequence proportion of the constituents. Additionally, genome skimming was able to identify one product that was a different Echinacea species entirely.

The Evolutionary History, Diversity, and Ecology of Willows (Salix L.) in the European Alps

The genus Salix (willows), with 33 species, represents the most diverse genus of woody plants in the European Alps. Many species dominate subalpine and alpine types of vegetation. Despite a long history of research on willows, the evolutionary and ecological factors for this species richness are poorly known. Here we will review recent progress in research on phylogenetic relationships, evolution, ecology, and speciation in alpine willows. Phylogenomic reconstructions suggest multiple colonization of the Alps, probably from the late Miocene onward, and reject hypotheses of a single radiation. Relatives occur in the Arctic and in temperate Eurasia. Most species are widespread in the European mountain systems or in the European lowlands. Within the Alps, species differ ecologically according to different elevational zones and habitat preferences. Homoploid hybridization is a frequent process in willows and happens mostly after climatic fluctuations and secondary contact. Breakdown of the ecological crossing barriers of species is followed by introgressive hybridization. Polyploidy is an important speciation mechanism, as 40% of species are polyploid, including the four endemic species of the Alps. Phylogenomic data suggest an allopolyploid origin for all taxa analyzed so far. Further studies are needed to specifically analyze biogeographical history, character evolution, and genome evolution of polyploids.

High Sequence Divergence but Limited Architectural Rearrangements in Organelle Genomes of Cyanophora (Glaucophyta) Species

Cyanophora is the glaucophyte model taxon. Following the sequencing of the nuclear genome of C. paradoxa, studies based on single organelle and nuclear molecular markers revealed previously unrecognized species diversity within this glaucophyte genus. Here, we present the complete plastid (ptDNA) and mitochondrial (mtDNA) genomes of C. kugrensii, C. sudae, and C. biloba. The respective sizes and coding capacities of both ptDNAs and mtDNAs are conserved among Cyanophora species with only minor differences due to specific gene duplications. Organelle phylogenomic analyses consistently recover the species C. kugrensii and C. paradoxa as a clade and C. sudae and C. biloba as a separate group. The phylogenetic affiliations of the four Cyanophora species are consistent with architectural similarities shared at the organelle genomic level. Genetic distance estimations from both organelle sequences are also consistent with phylogenetic and architecture evidence. Comparative analyses confirm that the Cyanophora mitochondrial genes accumulate substitutions at 3-fold higher rates than plastid counterparts, suggesting that mtDNA markers are more appropriate to investigate glaucophyte diversity and evolutionary events that occur at a population level. The study of complete organelle genomes is becoming the standard for species delimitation and is particularly relevant to study cryptic diversity in microbial groups.

A targeted sequence capture array for phylogenetics and population genomics in the Salicaceae

PREMISE

The family Salicaceae has proved taxonomically challenging, especially in the genus Salix, which is speciose and features frequent hybridization and polyploidy. Past efforts to reconstruct the phylogeny with molecular barcodes have failed to resolve the species relationships of many sections of the genus.

METHODS

We used the wealth of sequence data in the family to design sequence capture probes to target regions of 300-1200 bp of exonic regions of 972 genes.

RESULTS

We recovered sequence data for nearly all of the targeted genes in three species of Populus and three species of Salix. We present a species tree, discuss concordance among gene trees, and present population genomic summary statistics for these loci.

CONCLUSIONS

Our sequence capture array has extremely high capture efficiency within the genera Populus and Salix, resulting in abundant phylogenetic information. Additionally, these loci show promise for population genomic studies.

The low copy nuclear region, RPB2 as a novel DNA barcode region for species identification in the rattan genus Calamus (Arecaceae)

Taxonomic complexities, like environmental plasticity and homoplasy, make precise identification challenging in Calamus, the genus of spiny climbing palms of the subfamily Calamoideae (Arecaceae). In the present study, the species discriminatory power of twelve potential DNA barcode regions (rbcL, matK, psbA-trnH, rpoC, rpoB, psbK-psbI, atpF-atpH, psbZ-trnfM, ITS1, ITS2, PRK, and RPB2) were evaluated in 21 species of Calamus from the Western Ghats region of India, using distance, tree, and similarity based statistical methods. Except for the low copy nuclear region, RPB2, none of the tested plastid loci or nuclear loci ITS, either singly or in combinations, could discriminate all the species of Calamus due to low substitution rate of plastid regions and multiple copies of ITS respectively. The RPB2 locus showed highest species resolution with 96% accuracy in similarity based analysis, indicating its potential and efficiency as a barcode locus for the genus. The putative "Calamus gamblei complex" based on overlapping morphology was successfully resolved as six distinct, though closely related, species. The analysis also indicates that C. delessertianus is a morphological variant of C. dransfieldii. In spite of being a low copy nuclear gene region, RPB2 provided an efficient barcode to delineate Calamus species and has the potential to further extend its use as a prospective barcode to other Palm genera.

Phylogenomic Relationships and Evolution of Polyploid Salix Species Revealed by RAD Sequencing Data

Polyploidy is common in the genus Salix. However, little is known about the origin, parentage and genomic composition of polyploid species because of a lack of suitable molecular markers and analysis tools. We established a phylogenomic framework including species of all described sections of Eurasian shrub willows. We analyzed the genomic composition of seven polyploid willow species in comparison to putative diploid parental species to draw conclusions on their origin and the effects of backcrossing and post-origin evolution. We applied recently developed programs like SNAPP, HyDe, and SNiPloid to establish a bioinformatic pipeline for unravelling the complexity of polyploid genomes. RAD sequencing revealed 23,393 loci and 320,010 high quality SNPs for the analysis of relationships of 35 species of Eurasian shrub willows (Salix subg. Chamaetia/Vetrix). Polyploid willow species appear to be predominantly of allopolyploid origin. More ancient allopolyploidization events were observed for two hexaploid and one octoploid species, while our data suggested a more recent allopolyploid origin for the included tetraploids and identified putative parental taxa. SNiPloid analyses disentangled the different genomic signatures resulting from hybrid origin, backcrossing, and secondary post-origin evolution in the polyploid species. Our RAD sequencing data demonstrate that willow genomes are shaped by ancient and recent reticulate evolution, polyploidization, and post-origin divergence of species.

Complete plastome sequencing resolves taxonomic relationships among species of Calligonum L. (Polygonaceae) in China

BACKGROUND

Calligonum (Polygonaceae) is distributed from southern Europe through northern Africa to central Asia, and is typically found in arid, desert regions. Previous studies have revealed that standard DNA barcodes fail to discriminate Calligonum species. In this study, the complete plastid genomes (plastome) for 32 accessions of 21 Calligonum species is sequenced to not only generate the first complete plastome sequence for the genus Calligonum but to also 1) Assess the ability of the complete plastome sequence to discern species within the group, and 2) screen the plastome sequence for a cost-effective DNA barcode that can be used in future studies to resolve taxonomic relationships within the group.

RESULTS

The whole plastomes of Calligonum species possess a typical quadripartite structure. The size of the Calligonum plastome is approximately 161 kilobase pairs (kbp), and encodes 113 genes, including 79 protein-coding genes, 30 tRNA genes, and four rRNA genes. Based on ML phylogenetic tree analyses, the complete plastome has higher species identification (78%) than combinations of standard DNA barcodes (rbcL + matK + nrITS, 56%). Five newly screened gene regions (ndhF, trnS-G, trnC-petN, ndhF-rpl32, rpl32-trnL) had high species resolution, where ndhF and trnS-G were able to distinguish the highest proportion of Calligonum species (56%).

CONCLUSIONS

The entire plastid genome was the most effective barcode for the genus Calligonum, although other gene regions showed great potential as taxon-specific barcodes for species identification in Calligonum.

Microfluidic Enrichment Barcoding (MEBarcoding): a new method for high throughput plant DNA barcoding

DNA barcoding is a valuable tool to support species identification with broad applications from traditional taxonomy, ecology, forensics, food analysis, and environmental science. We introduce Microfluidic Enrichment Barcoding (MEBarcoding) for plant DNA Barcoding, a cost-effective method for high-throughput DNA barcoding. MEBarcoding uses the Fluidigm Access Array to simultaneously amplify targeted regions for 48 DNA samples and hundreds of PCR primer pairs (producing up to 23,040 PCR products) during a single thermal cycling protocol. As a proof of concept, we developed a microfluidic PCR workflow using the Fluidigm Access Array and Illumina MiSeq. We tested 96 samples for each of the four primary DNA barcode loci in plants: rbcL, matK, trnH-psbA, and ITS. This workflow was used to build a reference library for 78 families and 96 genera from all major plant lineages - many currently lacking in public databases. Our results show that this technique is an efficient alternative to traditional PCR and Sanger sequencing to generate large amounts of plant DNA barcodes and build more comprehensive barcode databases.

Salix transect of Europe: patterns in the distribution of willow-feeding psyllids (Hemiptera: Psylloidea) from Greece to arctic Norway

Background Psyllids are oligophagous phytophagous insects with many specialist willow (Salix spp.) feeding species in two genera (Cacopsylla and Bactericera). We examine the patterns of distribution and co-occurrence of willow-feeding species at 42 willow sites across Europe forming a transect from Greece (lat. 38.8 °N) to arctic Norway (lat. 70.6 °N). The transect and sites have been described in previous papers. New information A total of 1245 individual psyllids were examined from 23 species of willow over the transect, representing 17 willow-feeding species (11 Cacopsylla and 6 Bactericera). Numerous species were very widely distributed, with two species, Bactericera albiventris (Foerster, 1848) and Cacopsylla pulchra (Zetterstedt, 1840), occurring from Greece to Finland. Other widespread species (Romania to Finland) were Cacopsylla ambigua (Foerster, 1848) and Bactericera curvatinervis (Foerster, 1848). The mean number of psyllid species per site was 2.4 (1.3 Cacopsylla, 1.1 Bactericera).

Genomic biosurveillance of forest invasive alien enemies: A story written in code

The world's forests face unprecedented threats from invasive insects and pathogens that can cause large irreversible damage to the ecosystems. This threatens the world's capacity to provide long-term fiber supply and ecosystem services that range from carbon storage, nutrient cycling, and water and air purification, to soil preservation and maintenance of wildlife habitat. Reducing the threat of forest invasive alien species requires vigilant biosurveillance, the process of gathering, integrating, interpreting, and communicating essential information about pest and pathogen threats to achieve early detection and warning and to enable better decision-making. This process is challenging due to the diversity of invasive pests and pathogens that need to be identified, the diverse pathways of introduction, and the difficulty in assessing the risk of establishment. Genomics can provide powerful new solutions to biosurveillance. The process of invasion is a story written in four chapters: transport, introduction, establishment, and spread. The series of processes that lead to a successful invasion can leave behind a DNA signature that tells the story of an invasion. This signature can help us understand the dynamic, multistep process of invasion and inform management of current and future introductions. This review describes current and future application of genomic tools and pipelines that will provide accurate identification of pests and pathogens, assign outbreak or survey samples to putative sources to identify pathways of spread, and assess risk based on traits that impact the outbreak outcome.

Comparison of different annotation tools for characterization of the complete chloroplast genome of Corylus avellana cv Tombul

BACKGROUND

Several bioinformatics tools have been designed for assembly and annotation of chloroplast (cp) genomes, making it difficult to decide which is most useful and applicable to a specific case. The increasing number of plant genomes provide an opportunity to accurately obtain cp genomes from whole genome shotgun (WGS) sequences. Due to the limited genetic information available for European hazelnut (Corylus avellana L.) and as part of a genome sequencing project, we analyzed the complete chloroplast genome of the cultivar 'Tombul' with multiple annotation tools.

RESULTS

Three different annotation strategies were tested, and the complete cp genome of C. avellana cv Tombul was constructed, which was 161,667 bp in length, and had a typical quadripartite structure. A large single copy (LSC) region of 90,198 bp and a small single copy (SSC) region of 18,733 bp were separated by a pair of inverted repeat (IR) regions of 26,368 bp. In total, 125 predicted functional genes were annotated, including 76 protein-coding, 25 tRNA, and 4 rRNA unique genes. Comparative genomics indicated that the cp genome sequences were relatively highly conserved in species belonging to the same order. However, there were still some variations, especially in intergenic regions, that could be used as molecular markers for analyses of phylogeny and plant identification. Simple sequence repeat (SSR) analysis showed that there were 83 SSRs in the cp genome of cv Tombul. Phylogenetic analysis suggested that C. avellana cv Tombul had a close affinity to the sister group of C. fargesii and C. chinensis, and then a closer evolutionary relationship with Betulaceae family than other species of Fagales.

CONCLUSION

In this study, the complete cp genome of Corylus avellana cv Tombul, the most widely cultivated variety in Turkey, was obtained and annotated, and additionally phylogenetic relationships were predicted among Fagales species. Our results suggest a very accurate assembly of chloroplast genome from next generation whole genome shotgun (WGS) sequences. Enhancement of taxon sampling in Corylus species provide genomic insights into phylogenetic analyses. The nucleotide sequences of cv Tombul cp genomes can provide comprehensive genetic insight into the evolution of genus Corylus.

Ancient events and climate adaptive capacity shaped distinct chloroplast genetic structure in the oak lineages

BACKGROUND

Understanding the origin of genetic variation is the key to predict how species will respond to future climate change. The genus Quercus is a species-rich and ecologically diverse woody genus that dominates a wide range of forests and woodland communities of the Northern Hemisphere. Quercus thus offers a unique opportunity to investigate how adaptation to environmental changes has shaped the spatial genetic structure of closely related lineages. Furthermore, Quercus provides a deep insight into how tree species will respond to future climate change. This study investigated whether closely related Quercus lineages have similar spatial genetic structures and moreover, what roles have their geographic distribution, ecological tolerance, and historical environmental changes played in the similar or distinct genetic structures.

RESULTS

Despite their close relationships, the three main oak lineages (Quercus sections Cyclobalanopsis, Ilex, and Quercus) have different spatial genetic patterns and occupy different climatic niches. The lowest level and most homogeneous pattern of genetic diversity was found in section Cyclobalanopsis, which is restricted to warm and humid climates. The highest genetic diversity and strongest geographic genetic structure were found in section Ilex, which is due to their long-term isolation and strong local adaptation. The widespread section Quercus is distributed across the most heterogeneous range of environments; however, it exhibited moderate haplotype diversity. This is likely due to regional extinction during Quaternary climatic fluctuation in Europe and North America.

CONCLUSIONS

Genetic variations of sections Ilex and Quercus were significantly predicted by geographic and climate variations, while those of section Cyclobalanopsis were poorly predictable by geographic or climatic diversity. Apart from the different historical environmental changes experienced by different sections, variation of their ecological or climatic tolerances and physiological traits induced varying responses to similar environment changes, resulting in distinct spatial genetic patterns.

Phylogenetic analysis of Fritillaria cirrhosa D. Don and its closely related species based on complete chloroplast genomes

Fritillaria cirrhosa D. Don, whose bulb is used in a well-known traditional Chinese medicine to relieve cough and eliminate phlegm, is one of the most important medicinal plants of Fritillaria L. The species is widely distributed among the alpine regions in southwestern China and possesses complex morphological variations in different distributions. A series of newly related species were reported, based on obscure morphological differences. As a result, F. cirrhosa and its closely related species constitute a taxonomically complex group. However, it is difficult to accurately identify these species and reveal their phylogenetic relationships using traditional taxonomy. Molecular markers and gene fragments have been adopted but they are not able to afford sufficient phylogenetic resolution in the genus. Here, we report the complete chloroplast genome sequences of F. cirrhosa and its closely related species using next generation sequencing (NGS) technology. Eight plastid genomes ranged from 151,058 bp to 152,064 bp in length and consisted of 115 genes. Gene content, gene order, GC content, and IR/SC boundary structures were highly similar among these genomes. SSRs and five large repeat sequences were identified and the total number of them ranged from 73 to 79 and 63 to 75, respectively. Six highly divergent regions were successfully identified that could be used as potential genetic markers of Fritillaria. Phylogenetic analyses revealed that eight Fritillaria species were clustered into three clades with strong supports and F. cirrhosa was closely related to F. przewalskii and F. sinica. Overall, this study indicated that the complete chloroplast genome sequence was an efficient tool for identifying species in taxonomically complex groups and exploring their phylogenetic relationships.

Exploring DNA quantity and quality from raw materials to botanical extracts

OBJECTIVES

The aim of this study was to explore the variability in DNA quality and quantity along a gradient of industrial processing of botanical ingredients from raw materials to extracts.

METHODS

A data matrix was assembled for 1242 botanical ingredient samples along a gradient of industrial processing commonly used in the Natural Health Product (NHP) industry. Multivariate statistics was used to explore dependant variables for quality and quantity. The success of attaining a positive DNA test result along a gradient of industrial processing was compared among four biotechnologies: DNA barcoding, NGS, Sanger sequencing and qPCR.

RESULTS

There was considerable variance in DNA quality and quantity among the samples, which could be interpreted along a gradient from raw materials with greater quantities (50-120 ng/μL) of DNA and longer DNA (400-500bp) sequences to extracts, which were characterized by lower quantities (0.1-10.0 ng/μL) and short fragments (50-150bp).

CONCLUSIONS

Targeted molecular diagnostic tests for species identity can be used in the NHP industry for raw and processed samples. Non-targeted tests or the use of NGS for any identity test needs considerable research and development and must be validated before it can be used in commercial operations as these methods are subject to considerable risk of false negative and positive results. Proper use of these tools can be used to ensure ingredient authenticity, and to avert adulteration, and contamination with plants that are a health concern. Lastly these tools can be used to prevent the exploitation of rare herbal species and the harvesting of native biodiversity for commercial purposes.

Phylogenomic Approaches to DNA Barcoding of Herbal Medicines: Developing Clade-Specific Diagnostic Characters for Berberis

DNA barcoding of herbal medicines has been mainly concerned with authentication of products in trade and has raised awareness of species substitution and adulteration. More recently DNA barcodes have been included in pharmacopoeias, providing tools for regulatory purposes. The commonly used DNA barcoding regions in plants often fail to resolve identification to species level. This can be especially challenging in evolutionarily complex groups where incipient or reticulate speciation is ongoing. In this study, we take a phylogenomic approach, analyzing whole plastid sequences from the evolutionarily complex genus Berberis in order to develop DNA barcodes for the medicinally important species Berberis aristata. The phylogeny reconstructed from an alignment of ∼160 kbp of chloroplast DNA for 57 species reveals that the pharmacopoeial species in question is polyphyletic, complicating development of a species-specific DNA barcode. Instead we propose a DNA barcode that is clade specific, using our phylogeny to define Operational Phylogenetic Units (OPUs). The plastid alignment is then reduced to small, informative DNA regions including nucleotides diagnostic for these OPUs. These DNA barcodes were tested on commercial samples, and shown to discriminate plants in trade and therefore to meet the requirement of a pharmacopoeial standard. The proposed method provides an innovative approach for inferring DNA barcodes for evolutionarily complex groups for regulatory purposes and quality control.

Prevalence of isomeric plastomes and effectiveness of plastome super-barcodes in yews (Taxus) worldwide

Taxus (yew) is both the most species-rich and taxonomically difficult genus in Taxaceae. To date, no study has elucidated the complexities of the plastid genome (plastome) or examined the possibility of whole plastomes as super-barcodes across yew species worldwide. In this study, we sequenced plastomes from two to three individuals for each of the 16 recognized yew species (including three potential cryptics) and Pseudotaxus chienii. Our comparative analyses uncovered several gene loss events that independently occurred in yews, resulting in a lower plastid gene number than other Taxaceous genera. In Pseudotaxus and Taxus, we found two isomeric arrangements that differ by the orientation of a 35 kb fragment flanked by "trnQ-IRs". These two arrangements exist in different ratios within each sampled individual, and intraspecific shifts in major isomeric arrangements are first reported here in Taxus. Moreover, we demonstrate that entire plastomes can be used to successfully discriminate all Taxus species with 100% support, suggesting that they are useful as super-barcodes for species identification. We also propose that accD and rrn16-rrn23 are promising special barcodes to discriminate yew species. Our newly developed Taxus plastomic sequences provide a resource for super-barcodes and conservation genetics of several endangered yews and serve as comprehensive data to improve models of plastome complexity in Taxaceae as a whole and authenticate Taxus species.

Species delimitation and interspecific relationships of the endangered herb genus Notopterygium inferred from multilocus variations

Species identification and discrimination is the basis of biodiversity research. In general, it is considered that numerous nucleotide variations (e.g., whole chloroplast genomes) can identify species with higher resolution than a few loci, e.g., partial chloroplast or nuclear gene fragments. In this study, we tested this hypothesis by sampling population genetics samples of the endangered herb genus Notopterygium. We sequenced the complete plastomes, five nuclear gene regions, three chloroplast DNA fragments, and a nuclear internal transcribed spacer (nrITS) region for 18 populations sampled throughout most of the geographic ranges of all six Notopterygium species. Species identification analysis showed that four DNA barcodes (matK, rbcL, trnS-trnG, and nrITS) and/or combinations of these markers achieved Notopterygium species discrimination at higher resolution than the general plastomes and nuclear gene sequences. In particular, nrITS had the highest discriminatory power among all of the individual markers. Molecular data sets and morphological evidence indicated that all six Notopterygium species could be reclassified unambiguously to four putative species clades. N. oviforme and N. franchetii had the closest relationship. Molecular dating showed that the origin and divergence of Notopterygium species was significantly associated with geological and climatic fluctuations during the middle of the Pliocene. In conclusion, our results suggest that a few nucleotide variations can achieve species discrimination with higher resolution than numerous plastomes and general nuclear gene fragments when discerning related Notopterygium species.

RAD sequencing resolved phylogenetic relationships in European shrub willows (Salix L. subg. Chamaetia and subg. Vetrix) and revealed multiple evolution of dwarf shrubs

The large and diverse genus Salix L. is of particular interest for decades of biological research. However, despite the morphological plasticity, the reconstruction of phylogenetic relationships was so far hampered by the lack of informative molecular markers. Infrageneric classification based on morphology separates dwarf shrubs (subg. Chamaetia) and taller shrubs (subg. Vetrix), while previous phylogenetic studies placed species of these two subgenera just in one largely unresolved clade. Here we want to test the utility of genomic RAD sequencing markers for resolving relationships at different levels of divergence in Salix. Based on a sampling of 15 European species representing 13 sections of the two subgenera, we used five different RAD sequencing datasets generated by ipyrad to conduct phylogenetic analyses. Additionally we reconstructed the evolution of growth form and analyzed the genetic composition of the whole clade. The results showed fully resolved trees in both ML and BI analysis with high statistical support. The two subgenera Chamaetia and Vetrix were recognized as nonmonophyletic, which suggests that they should be merged. Within the Vetrix/Chamaetia clade, a division into three major subclades could be observed. All species were confirmed to be monophyletic. Based on our data, arctic-alpine dwarf shrubs evolved four times independently. The structure analysis showed five mainly uniform genetic clusters which are congruent in sister relationships observed in the phylogenies. Our study confirmed RAD sequencing as a useful genomic tool for the reconstruction of relationships on different taxonomic levels in the genus Salix.

DNA barcoding herbaceous and woody plant species at a subalpine forest dynamics plot in Southwest China

Although DNA barcoding has been widely used to identify plant species composition in temperate and tropical ecosystems, relatively few studies have used DNA barcodes to document both herbaceous and woody components of forest plot. A total of 201 species (72 woody species and 129 herbaceous species) representing 135 genera distributed across 64 families of seed plants were collected in a 25 ha CForBio subalpine forest dynamics plot. In total, 491 specimens were screened for three DNA regions of the chloroplast genome (rbcL, matK, and trnH-psbA) as well as the internal transcribed spacers (ITS) of nuclear ribosomal DNA. We quantified species resolution for each barcode separately or in combination using a ML tree-based method. Amplification and sequencing success were highest for rbcL, followed by trnH-psbA, which performed better than ITS and matK. The rbcL + ITS barcode had slightly higher species resolution rates (88.60%) compared with rbcL + matK (86.60%) and rbcL + trnH-psbA (86.01%). The addition of trnH-psbA or ITS to the rbcL + matK barcode only marginally increased species resolution rates, although in combination the four barcodes had the highest discriminatory power (90.21%). The situations where DNA barcodes did not discriminate among species were typically associated with higher numbers of co-occurring con-generic species. In addition, herbaceous species were much better resolved than woody species. Our study represents one of the first applications of DNA barcodes in a subalpine forest dynamics plot and contributes to our understanding of patterns of genetic divergence among woody and herbaceous plant species.

DNA barcoding a taxonomically complex hemiparasitic genus reveals deep divergence between ploidy levels but lack of species-level resolution

DNA barcoding is emerging as a useful tool not only for species identification but also for studying evolutionary and ecological processes. Although plant DNA barcodes do not always provide species-level resolution, the generation of large DNA barcode data sets can provide insights into the mechanisms underlying the generation of species diversity. Here, we study evolutionary processes in taxonomically complex British Euphrasia (Orobanchaceae), a group with multiple ploidy levels, frequent self-fertilization, young species divergence and widespread hybridization. We use a phylogenetic approach to investigate the colonization history of British Euphrasia, followed by a DNA barcoding survey and population genetic analyses to reveal the causes of shared sequence variation. Phylogenetic analysis shows Euphrasia have colonized Britain from mainland Europe on multiple occasions. DNA barcoding reveals that no British Euphrasia species has a consistent diagnostic sequence profile, and instead, plastid haplotypes are either widespread across species, or are population specific. The partitioning of nuclear genetic variation suggests differences in ploidy act as a barrier to gene exchange, while the divergence between diploid and tetraploid ITS sequences supports the polyploids being allotetraploid in origin. Overall, these results show that even when lacking species-level resolution, analyses of DNA barcoding data can reveal evolutionary patterns in taxonomically complex genera.

A nuclear DNA barcode for eastern North American oaks and application to a study of hybridization in an Arboretum setting

DNA barcoding has proved difficult in a number of woody plant genera, including the ecologically important oak genus Quercus. In this study, we utilized restrictionsite-associated DNA sequencing (RAD-seq) to develop an economical single nucleotide polymorphism (SNP) DNA barcoding system that suffices to distinguish eight common, sympatric eastern North American white oak species. Two de novo clustering pipelines, PyRAD and Stacks, were used in combination with postclustering bioinformatic tools to generate a list of 291 potential SNPs, 80 of which were included in a barcoding toolkit that is easily implemented using MassARRAY mass spectrometry technology. As a proof-of-concept, we used the genotyping toolkit to infer potential hybridization between North American white oaks transplanted outside of their native range (Q. michauxii, Q. montana, Q muehlenbergii/Q. prinoides, and Q. stellata) into a horticultural collection surrounded by natural forests of locally native trees (Q. alba and Q. macrocarpa) in the living collection at The Morton Arboretum (Lisle, IL, USA). Phylogenetic and clustering analyses suggested low rates of hybridization between cultivated and native species, with the exception of one Q. michauxii mother tree, the acorns of which exhibited high admixture from either Q. alba or Q. stellata and Q. macrocarpa, and a hybrid between Q. stellata that appears to have backcrossed almost exclusively to Q. alba. Together, RAD-seq and MassARRAY technologies allow for efficient development and implementation of a multispecies barcode for one of the more challenging forest tree genera.

The Application and Limitation of Universal Chloroplast Markers in Discriminating East Asian Evergreen Oaks

The East Asian subtropics mostly occupied by evergreen broad-leaved forests (EBLFs), is one of the global diversity centers for evergreen oaks. Evergreen oaks are keystone canopy trees in EBLFs with important ecosystem function and crucial significance for regional biodiversity conservation. However, the species composition and diversity of Asian evergreen oaks are poorly understood. Here, we test whether the four chloroplast markers atpI-atpH, matK, psbA-trnH, and ycf1, can discriminate the two evergreen oak sections in Asia - Cyclobalanopsis and Ilex. Two hundred and seventy-two individuals representing 57 species were scanned and 17 species from other oaks sections were included for phylogenetic reconstruction. The genetic diversity of the Quercus sections was also compared. Overall, we found that universal chloroplast DNA (cpDNA) barcoding markers could resolve two clades in Quercus, i.e., subgenus Cerris (Old World Clade) and subgenus Quercus (New World Clade). The chloroplast markers distinguished the main sections, with few exceptions. Each cpDNA region showed no barcoding gap and none of them provided good resolution at the species level. The best species resolution (27.78%) was obtained when three or four markers were combined and analyzed using BLAST. The high conservation of the cpDNA and complicated evolutionary patterns, due to incomplete lineage sorting, interspecific hybridization and introgressions may hinder the ability of cpDNA markers to discriminate different species. When comparing diversification pattern across Quercus sections (Cyclobalanopsis, Ilex, Cerris, Quercus, and Protobalanus), we found that section Ilex was the most genetically diverse, and section Cyclobalanopsis was lower genetically diverse. This diversification pattern may have resulted from the interplay of the Eurasia Cenozoic tectonic movements, climate changes and different niches of their ancestral lineages.