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.

Species delimitation 4.0: integrative taxonomy meets artificial intelligence.

Although species are central units for biological research, recent findings in genomics are raising awareness that what we call species can be ill-founded entities due to solely morphology-based, regional species descriptions. This particularly applies to groups characterized by intricate evolutionary processes such as hybridization, polyploidy, or asexuality. Here, challenges of current integrative taxonomy (genetics/genomics + morphology + ecology, etc.) become apparent: different favored species concepts, lack of universal characters/markers, missing appropriate analytical tools for intricate evolutionary processes, and highly subjective ranking and fusion of datasets. Now, integrative taxonomy combined with artificial intelligence under a unified species concept can enable automated feature learning and data integration, and thus reduce subjectivity in species delimitation. This approach will likely accelerate revising and unraveling eukaryotic biodiversity.

Contrasting levels of ß-diversity and underlying phylogenetic trends indicate different paths to chemical diversity in highland and lowland willow species.

Diverse specialised metabolites contributed to the success of vascular plants in colonising most terrestrial habitats. Understanding how distinct aspects of chemical diversity arise through heterogeneous environmental pressures can help us understand the effects of abiotic and biotic stress on plant evolution and community assembly. We examined highland and lowland willow species within a phylogenetic framework to test for trends in their chemical α-diversity (richness) and ß-diversity (variation among species sympatric in elevation). We show that differences in chemistry among willows growing at different elevations occur mainly through shifts in chemical ß-diversity and due to convergence or divergence among species sharing their elevation level. We also detect contrasting phylogenetic trends in concentration and α-diversity of metabolites in highland and lowland willow species. The resulting elevational patterns contribute to the chemical diversity of willows and suggest that variable selective pressure across ecological gradients may, more generally, underpin complex changes in plant chemistry.

Insights into the Taxonomically Challenging Hexaploid Alpine Shrub Willows of Salix Sections Phylicifoliae and Nigricantes (Salicaceae).

The complex genomic composition of allopolyploid plants leads to morphologically diverse species. The traditional taxonomical treatment of the medium-sized, hexaploid shrub willows distributed in the Alps is difficult based on their variable morphological characters. In this study, RAD sequencing data, infrared-spectroscopy, and morphometric data are used to analyze the phylogenetic relationships of the hexaploid species of the sections Nigricantes and Phylicifoliae in a phylogenetic framework of 45 Eurasian Salix species. Both sections comprise local endemics as well as widespread species. Based on the molecular data, the described morphological species appeared as monophyletic lineages (except for S. phylicifolia s.str. and S. bicolor, which are intermingled). Both sections Phylicifoliae and Nigricantes are polyphyletic. Infrared-spectroscopy mostly confirmed the differentiation of hexaploid alpine species. The morphometric data confirmed the molecular results and supported the inclusion of S. bicolor into S. phylicifolia s.l., whereas the alpine endemic S. hegetschweileri is distinct and closely related to species of the section Nigricantes. The genomic structure and co-ancestry analyses of the hexaploid species revealed a geographical pattern for widespread S. myrsinifolia, separating the Scandinavian from the alpine populations. The newly described S. kaptarae is tetraploid and is grouped within S. cinerea. Our data reveal that both sections Phylicifoliae and Nigricantes need to be redefined.

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.

Ancient DNA extraction methods for herbarium specimens: When is it worth the effort?

Premise: Herbaria harbor a tremendous number of plant specimens that are rarely used for molecular systematic studies, largely due to the difficulty in extracting sufficient amounts of high-quality DNA from the preserved plant material. Methods: We compared the standard Qiagen DNeasy Plant Mini Kit and a specific protocol for extracting ancient DNA (aDNA) (the N-phenacylthiazolium bromide and dithiothreitol [PTB-DTT] extraction method) from two different plant genera (Xanthium and Salix). The included herbarium materials covered about two centuries of plant collections. To analyze the success of DNA extraction using each method, a subset of samples was subjected to a standard library preparation as well as target-enrichment approaches. Results: The PTB-DTT method produced a higher DNA yield of better quality than the Qiagen kit; however, extracts from the Qiagen kit over a certain DNA yield and quality threshold produced comparable sequencing results. The sequencing resulted in high proportions of endogenous reads. We were able to successfully sequence 200-year-old samples. Discussion: This method comparison revealed that, for younger specimens, DNA extraction using a standard kit might be sufficient. For old and precious herbarium specimens, aDNA extraction methods are better suited to meet the requirements for next-generation sequencing.

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.

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.

Ecological range shift in the polyploid members of the South American genus Fosterella (Bromeliaceae).

Background and Aims: The distribution of polyploidy along a relatively steep Andean elevation and climatic gradient is studied using the genus Fosterella L.B. Sm. (Bromeliaceae) as a model system. Ecological differentiation of cytotypes and the link of polyploidy with historical biogeographic processes such as dispersal events and range shift are assessed. Methods: 4',6-Diamidino-2-phenylindole (DAPI) staining of nuclei and flow cytometry were used to estimate the ploidy levels of 159 plants from 22 species sampled throughout the distribution range of the genus. Ecological differentiation among ploidy levels was tested by comparing the sets of climatic variables. Ancestral chromosome number reconstruction was carried out on the basis of a previously generated phylogeographic framework. Key Results: This study represents the first assessment of intrageneric, intraspecific and partially intrapopulational cytotype diversity in a genus of the Bromeliaceae family. In Fosterella , the occurrence of polyploidy was limited to the phylogenetically isolated penduliflora and rusbyi groups. Cytotypes were found to be ecologically differentiated, showing that polyploids preferentially occupy colder habitats with high annual temperature variability (seasonality). The combined effects of biogeographic history and adaptive processes are presumed to have shaped the current cytotype distribution in the genus. Conclusions: The results provide indirect evidence for both adaptive ecological and non-adaptive historical processes that jointly influenced the cytotype distribution in the predominantly Andean genus Fosterella (Bromeliaceae). The results also exemplify the role of polyploidy as an important driver of speciation in a topographically highly structured and thus climatically diverse landscape.