Speciation via hybridization and polyploidization is a major evolutionary force in plant evolution but is still poorly understood for neopolyploid groups. Challenges are attributed to high heterozygosity, low genetic divergence, and missing information on progenitors, ploidy, and reproduction. We study the large Eurasian Ranunculus auricomus species complex and use a comprehensive workflow integrating reduced-representation sequencing (RRS) genomic data to unravel reticulate evolution, genome diversity and composition of polyploids.We rely on 97 312 restriction site-associated DNA sequencing (RAD-Seq) loci, 576 targeted nuclear genes (48 phased), and 71 plastid regions derived from 78 polyploid apomictic taxa and four diploid and one tetraploid putative sexual progenitor species. We applied (phylo)genomic structure, network, and single nucleotide polymorphism (SNP)-origin analyses.Results consistently showed only 3-5 supported and geographically structured polyploid genetic groups, each containing extant sexual and one unknown progenitor species. Combined analyses demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual progenitor species. Young allotetraploids were characterized by subgenome dominance and nonhybrid SNPs, suggesting substantial post-origin but little lineage-specific evolution.The biodiversity of neopolyploid complexes can result from multiple hybrid origins involving different progenitors and substantial post-origin evolution (e.g. homoeologous exchanges, hybrid segregation, gene flow). Reduced-representation sequencing genomic data including multi-approach information is efficient to delimit shallow reticulate relationships.

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.

Herbaria harbor a tremendous amount of plant specimens that are rarely used for plant systematic studies. The main reason is the difficulty to extract a decent quantity of good quality DNA from the preserved plant material. While the extraction of ancient DNA in animals is well established, studies including old plant material are still underrepresented. In our study we compared the standard Qiagen DNeasy Plant Mini Kit and a specific PTB-DTT protocol on to two different plant genera (Xanthium L. and Salix L.).  The included herbarium material covered about two centuries of plant collections. A selected subset of samples was used for a standard library preparation as well as a target enrichment approach. The results revealed that PTB-PTT resulted in higher quantity and quality regarding DNA yield. Despite the lower overall yield of DNA, the Qiagen Kit resulted in better sequencing results regarding the number of filtered and mapped reads. We were able to successfully sequence a sample from 1820 and conclude that it is possible to include old herbarium specimens in NGS approaches. This opens a treasure box in phylogenomic research.

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.

Complex genome evolution of young polyploid complexes is poorly understood. Besides challenges caused by hybridization, polyploidization, and incomplete lineage sorting, bioinformatic analyses are often exacerbated by missing information on progenitors, ploidy, and reproduction modes. By using a comprehensive, self-developed bioinformatic pipeline covering tree, structure, network, and SNP-origin analyses, we for the first time unraveled polyploid phylogenetic relationships and genome evolution within the large Eurasian Ranunculus auricomus species complex comprising more than 840 taxa. Our results rely on 97,312 genomic RADseq loci, target enrichment of 576 nuclear genes (48 phased), and 71 plastid regions (Hybseq; OMICS-data) derived from the 75 most widespread polyploid apomictic taxa and four di- and one tetraploid potential sexual progenitor species. Phylogenetic tree and structure analyses consistently showed 3-5 supported polyploid groups, each containing sexual progenitor species. In total, analyses revealed four diploid sexual progenitors and a one unknown, probably extinct progenitor, contributing to the genome composition of R. auricomus polyploids. Phylogenetic network, structure, and SNP-origin analyses based on RADseq loci and phased nuclear genes completed by plastid data demonstrated predominantly allopolyploid origins, each involving 2-3 different diploid sexual subgenomes. Allotetraploid genomes were characterized by subgenome dominance and large proportions of interspecific, non-hybrid SNPs, indicating an enormous degree of post-origin evolution (i.e., Mendelian segregation of the diploid hybrid generations, back-crossings, and gene flow due to facultative sexuality of apomicts), but only low proportions of lineage-specific SNPs. The R. auricomus model system is the first large European polyploid species complex studied with reduced representation OMICS data. Our bioinformatic pipeline underlines the importance of combining different approaches and datasets to successfully unveil how reticulate evolution and post-origin processes shape the diversity of polyploid plant complexes.