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

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

Evolutionary dynamics of Euphorbia carniolica suggest a complex Plio-Pleistocene history of understorey species of deciduous forest in southeastern Europe.

Deciduous forests form the dominant natural vegetation of Europe today, but were restricted to small refugia during Pleistocene cold stages, implying an evolutionary past shaped by recurrent range contractions and expansions. Cold-stage forest refugia were probably widespread in southern and central Europe, with the northwestern Balkan Peninsula being of particular importance. However, the actual number and location of deciduous forest refugia, as well as the connections between them, remain disputed. Here, we address the evolutionary dynamics of the deciduous forest understorey species Euphorbia carniolica as a proxy for past forest dynamics. To do so, we obtained genomic and morphometric data from populations representing the species' entire range, investigated phylogenetic position and intraspecific genetic variation, tested explicit demographic scenarios and applied species distribution models. Our data support two disjoint groups linked to separate refugia on the northwestern and central Balkan Peninsula. We find that genetic differentiation between groups started in the early Pleistocene via vicariance, suggesting a larger distribution in the past. Both refugia acted as sources for founder events to the southeastern Alps and the Carpathians; the latter were likely colonised before the last cold stage. In line with traditional views on the pre-Pleistocene origin of many southeastern European deciduous forest species, the origin of E. carniolica was dated to the late Pliocene. The fact that E. carniolica evolved at a time when a period of continuous forestation was ending in much of Eurasia provides an interesting biogeographical perspective on the past links between Eurasian deciduous forests and their biota.

Meeting in Liguria: Hybridisation between Apennine endemic Euphorbia barrelieri and western Mediterranean E. nicaeensis led to the allopolyploid origin of E. ligustica.

The Mediterranean Basin is renowned for its extremely rich biota and is considered as one of the 25 Global Biodiversity Hotspots, but its diversity is not homogeneously distributed. Outstanding in the number of (endemic) species are the Ligurian Alps (Italy). At the foot of the Ligurian Alps, little above the Mediterranean Sea, a disjunct occurrence of Italian endemic Euphorbia barrelieri was reported. Using an array of integrative methods ranging from cytogenetic (chromosome number and relative genome size estimation), over phylogenetic approaches (plastid, ITS and RAD sequencing) to multivariate morphometrics we disentangled the origin of these populations that were shown to be tetraploid. We performed phylogenetic analyses of the nuclear ITS and plastid regions of a broad taxonomic sampling of Euphorbia sect. Pithyusa to identify possible species involved in the origin of the tetraploid populations and then applied various analyses of RADseq data to identify the putative parental species. Our results have shown that the Ligurian populations of E. barrelieri are of allotetraploid origin that involved E. barrelieri and western Mediterranean E. nicaeensis as parental species. We thus describe a new species, E. ligustica, and hypothesise that its adaptation to similar environments in which E. barrelieri occurs, triggered development of similar morphology, whereas its genetic composition appears to be closer to that of E. nicaeensis. Our study emphasises the importance of polyploidisation for plant diversification, highlights the value of the Ligurian Alps as a hotspot of biodiversity and endemism and underlines the importance of integrative taxonomic approaches in uncovering cryptic diversity.

Phylogenomics and genome size evolution in Amomum s. s. (Zingiberaceae): Comparison of traditional and modern sequencing methods.

BACKGROUND AND AIMS: A targeted enrichment NGS approach was used to construct the phylogeny of Amomum Roxb. (Zingiberaceae). Phylogenies based on hundreds of nuclear genes, the whole plastome and the rDNA cistron were compared with an ITS-based phylogeny. Trends in genome size (GS) evolution were examined, chromosomes were counted and the geographical distribution of phylogenetic lineages was evaluated. METHODS: In total, 92 accessions of 54 species were analysed. ITS was obtained for 79 accessions, 37 accessions were processed with Hyb-Seq and sequences from 449 nuclear genes, the whole cpDNA, and the rDNA cistron were analysed using concatenation, coalescence and supertree approaches. The evolution of absolute GS was analysed in a phylogenetic and geographical context. The chromosome numbers of 12 accessions were counted. KEY RESULTS: Four groups were recognised in all datasets though their mutual relationships differ among datasets. While group A (A. subulatum and A. petaloideum) is basal to the remaining groups in the nuclear gene phylogeny, in the cpDNA topology it is sister to group B (A. repoeense and related species) and, in the ITS topology, it is sister to group D (the Elettariopsis lineage). The former Elettariopsis makes a monophyletic group. There is an increasing trend in GS during evolution. The largest GS values were found in group D in two tetraploid taxa, A. cinnamomeum and A. aff. biphyllum (both 2n = 96 chromosomes). The rest varied in GS (2C = 3.54-8.78 pg) with a constant chromosome number 2n = 48. There is a weak connection between phylogeny, GS and geography in Amomum. CONCLUSIONS: Amomum consists of four groups, and the former Elettariopsis is monophyletic. Species in this group have the largest GS. Two polyploids were found and GS greatly varied in the rest of Amomum.

Post-glacial range formation of temperate forest understorey herbs - Insights from a spatio-temporally explicit modelling approach.

Aim: Our knowledge of Pleistocene refugia and post-glacial recolonization routes of forest understorey plants is still very limited. The geographical ranges of these species are often rather narrow and show highly idiosyncratic, often fragmented patterns indicating either narrow and species-specific ecological tolerances or strong dispersal limitations. However, the relative roles of these factors are inherently difficult to disentangle. Location: Central and south-eastern Europe. Time period: 17,100 BP - present. Major taxa studied: Five understorey herbs of European beech forests: Aposeris foetida, Cardamine trifolia, Euphorbia carniolica, Hacquetia epipactis and Helleborus niger. Methods: We used spatio-temporally explicit modelling to reconstruct the post-glacial range dynamics of the five forest understorey herbs. We varied niche requirements, demographic rates and dispersal abilities across plausible ranges and simulated the spread of species from potential Pleistocene refugia identified by phylogeographical analyses. Then we identified the parameter settings allowing for the most accurate reconstruction of their current geographical ranges. Results: We found a largely homogenous pattern of optimal parameter settings among species. Broad ecological niches had to be combined with very low but non-zero rates of long-distance dispersal via chance events and low rates of seed dispersal over moderate distances by standard dispersal vectors. However, long-distance dispersal events, although rare, led to high variation among replicated simulation runs. Main conclusions: Small and fragmented ranges of many forest understorey species are best explained by a combination of broad ecological niches and rare medium- and long-distance dispersal events. Stochasticity is thus an important determinant of current species ranges, explaining the idiosyncratic distribution patterns of the study species despite strong similarities in refugia, ecological tolerances and dispersal abilities.

Repeat proliferation and partial endoreplication jointly shape the patterns of genome size evolution in orchids.

Although the evolutionary drivers of genome size change are known, the general patterns and mechanisms of plant genome size evolution are yet to be established. Here we aim to assess the relative importance of proliferation of repetitive DNA, chromosomal variation (including polyploidy), and the type of endoreplication for genome size evolution of the Pleurothallidinae, the most species-rich orchid lineage. Phylogenetic relationships between 341 Pleurothallidinae representatives were refined using a target enrichment hybrid capture combined with high-throughput sequencing approach. Genome size and the type of endoreplication were assessed using flow cytometry supplemented with karyological analysis and low-coverage Illumina sequencing for repeatome analysis on a subset of samples. Data were analyzed using phylogeny-based models. Genome size diversity (0.2-5.1 Gbp) was mostly independent of profound chromosome count variation (2n = 12-90) but tightly linked with the overall content of repetitive DNA elements. Species with partial endoreplication (PE) had significantly greater genome sizes, and genomic repeat content was tightly correlated with the size of the non-endoreplicated part of the genome. In PE species, repetitive DNA is preferentially accumulated in the non-endoreplicated parts of their genomes. Our results demonstrate that proliferation of repetitive DNA elements and PE together shape the patterns of genome size diversity in orchids.

Multiple auto- and allopolyploidisations marked the Pleistocene history of the widespread Eurasian steppe plant Astragalus onobrychis (Fabaceae).

The Eurasian steppes occupy a significant portion of the worldwide land surface and their biota have been affected by specific past range dynamics driven by ice ages-related climatic fluctuations. The dynamic alterations in conditions during the Pleistocene often triggered reticulate evolution and whole genome duplication events. Employing genomic, genetic and cytogenetic tools as well as morphometry we investigate the intricate evolution of Astragalus onobrychis, a widespread Eurasian steppe plant with diploid, tetraploid and octoploid cytotypes. To analyse the heteroploid RADseq dataset we employ both genotype-based and genotype-free methods that result in highly consistent results, and complement our inference with information from the plastid ycf1 region. We uncover a complex and reticulate evolutionary history, including at least one auto-tetraploidization event and two allo-octoploidization events; one of them involved also genetic contributions from other species, most likely A. goktschaicus. The present genetic structure points to the existence of four main clades within A. onobrychis, which only partly correspond to different ploidies. Time-calibrated diffusion models suggest that diversification within A. onobrychis was associated with ice age-related climatic fluctuations during the last million years. We finally argue for the usefulness of uniparentally inherited plastid markers, even in the genomic era, especially when investigating heteroploid systems.

Integrating phylogenomics, phylogenetics, morphometrics, relative genome size and ecological niche modelling disentangles the diversification of Eurasian Euphorbia seguieriana s. l. (Euphorbiaceae).

Next generation sequencing has revolutionised biology. Restriction-associated DNA sequencing (RADseq) has primarily been used to study infraspecific relationships but has also been applied in multi-species phylogenomic analyses. In this study, we used a combination of phylogenomic (with RADseq data) and phylogenetic (with sequences of the nuclear internal transcribed spacer, ITS) methods to explore relationships within the taxonomically intricate Euphorbia seguieriana s. l., one of the most widespread Euphorbia taxa inhabiting zonal and extrazonal steppes from Iberia to Central Asia. In the inferred phylogenies the southeastern Balkan and Anatolian populations were clearly separated, supporting the distinction of E. niciciana from E. seguieriana at the species level. Within E. seguieriana, the populations from the Caucasus, Iran, and easternmost Anatolia were sister to all other populations based on RADseq, making necessary the description of a new, morphologically divergent subspecies, E. seguieriana subsp. armeniaca. Conversely, additional studies are needed to understand the status of E. seguieriana subsp. hohenackeri, which is sympatric with E. seguieriana subsp. armeniaca. Niche analyses indicated that differences in the climatic niche between E. niciciana and E. seguieriana are relatively small compared with the climatic differences between the regions over which they are distributed. Contrary to previous believes, E. niciciana and E. seguieriana are allopatric and have likely diverged during the Pleistocene in two different glacial refugia as suggested by distribution modelling. Euphorbia niciciana nowadays has a submediterranean distribution, occupying habitats that are slightly warmer, moister, and less seasonal in temperature but more seasonal in precipitation than E. seguieriana, a characteristic species of continental steppes. Using flow cytometry, we demonstrate that the relative genome sizes of E. niciciana and E. seguieriana differ significantly. Additionally, multivariate morphometric analyses of 56 morphological characters indicated clear morphological divergence of the two species. Importantly, we also provide a revised taxonomic treatment including formal nomenclatural changes, an identification key and species descriptions. Our study demonstrates that an integrative approach, combining modern phylogenomic methods with traditional phylogenetic, cytogenetic, environmental and morphological analyses can result in satisfactorily resolved relationships in intricate groups of closely related species. Finally, phylogenetic inference using ITS sequences is still a useful tool for resolving relationships among the taxa at the species level, but the phylogenomic approach based on RADseq data certainly provides better resolution both among and within species.

Diversification dynamics and transoceanic Eurasian-Australian disjunction in the genus Picris (Compositae) induced by the interplay of shifts in intrinsic/extrinsic traits and paleoclimatic oscillations.

Understanding transcontinental biogeographic patterns has been one of the main foci of the field of biogeography. While multiple explanations for transcontinental disjunctions have been proposed, little is still known about the relative importance of intrinsic and extrinsic traits for the diversification dynamics of disjunct taxa. Here, we study the evolutionary history of the genus Picris L. (Compositae), a great model for investigating the diversification dynamics of transoceanic bipolar disjunct organisms. Ancestral state reconstructions indicate that the most recent common ancestor (MRCA) of Picris was a semelparous and heterocarpic herb that lived in unpredictable environments of North Africa and West Asia. Diversification analyses suggest a significant shift in speciation ca. 1 million years ago, likely associated with the onset of the mid-Pleistocene revolution. Longevity characters are correlated with the evolution of particular fruit types and with environmental conditions. Heterocarpic species are mostly semelparous herbs strongly linked with unpredictable habitats, while homocarpic taxa are mostly iteroparous plants occurring in predictable environments. Binary-state speciation and extinction analyses suggest that homocarpy, iteroparity, and habitats predictability accelerate diversification. Although the combination of homocarpy and iteroparity evolved in several lineages, only members of the P. hieracioides group were able to colonise Eurasia and expand to Australia by transoceanic dispersal. Those findings indicate that large-scale colonisation events depend on a complex interplay of intrinsic and extrinsic factors.

The Mediterranean: the cradle of Anthoxanthum (Poaceae) diploid diversity.

Background and Aims: Knowledge of diploid phylogeny and ecogeography provide a foundation for understanding plant evolutionary history, diversification patterns and taxonomy. The genus Anthoxanthum (vernal grasses, Poaceae) represents a taxonomically intricate polyploid complex with large phenotypic variation and poorly resolved evolutionary relationships. The aims of the study were to reveal: (1) evolutionary lineages of the diploid taxa and their genetic differentiation; (2) the past distribution of the rediscovered 'Mediterranean diploid'; and (3) possible migration routes of diploids in the Mediterranean. Methods: A combined approach involving sequencing of two plastid regions ( trnL-trnF and rpl32-trnL ), nrDNA ITS, rDNA FISH analyses, climatic niche characterization and spatio-temporal modelling was used. Key Results: Among the examined diploid species, only two well-differentiated evolutionary lineages were recognized: Anthoxanthum gracile and A. alpinum . The other taxa - A. aristatum, A. ovatum, A. maderense and the 'Mediterranean diploid' - form a rather intermixed group based on the examined molecular data. In situ rDNA localization enabled identification of the ancestral Anthoxanthum karyotype, shared by A. gracile and two taxa from the crown group. For the studied taxa, ancestral location probabilities for six discrete geographical regions in the Mediterranean were proposed and likely scenarios of gradual expansion from them were suggested. Modelling past and present distributions shows that the 'Mediterranean diploid' has already been occurring in the same localities for 120 000 years. Conclusions: Highly congruent results were obtained and dated the origin and first diversification of Anthoxanthum to the Miocene. The later divergence probably took place in the Pleistocene and started polyploid evolution within the genus. The most recent diversification event is still occurring, and incomplete lineage sorting prevents full diversification of taxa at the molecular level, despite clear separation based on climatic niches. The 'Mediterranean diploid' is hypothesized to be a possible relic of the most recent common ancestor of Anthoxanthum due to their sharing of ancestral features.

Northern glacial refugia and altitudinal niche divergence shape genome-wide differentiation in the emerging plant model Arabidopsis arenosa.

Quaternary climatic oscillations profoundly impacted temperate biodiversity. For many diverse yet undersampled areas, however, the consequences of this impact are still poorly known. In Europe, particular uncertainty surrounds the role of Balkans, a major hotspot of European diversity, in postglacial recolonization of more northerly areas, and the Carpathians, a debatable candidate for a northern 'cryptic' glacial refugium. Using genome-wide SNPs and microsatellites, we examined how the interplay of historical processes and niche shifts structured genetic diversity of diploid Arabidopsis arenosa, a little-known member of the plant model genus that occupies a wide niche range from sea level to alpine peaks across eastern temperate Europe. While the northern Balkans hosted one isolated endemic lineage, most of the genetic diversity was concentrated further north in the Pannonian Basin and the Carpathians, where it likely survived the last glaciation in northern refugia. Finally, a distinct postglacial environment in northern Europe was colonized by populations of admixed origin from the two Carpathian lineages. Niche differentiation along altitude-related bioclimatic gradients was the main trend in the phylogeny of A. arenosa. The most prominent niche shifts, however, characterized genetically only slightly divergent populations that expanded into narrowly defined alpine and northern coastal postglacial environments. Our study highlights the role of eastern central European mountains not only as refugia for unique temperate diversity but also sources for postglacial expansion into novel high-altitude and high-latitude niches. Knowledge of distinct genetic substructure of diploid A. arenosa also opens new opportunities for follow-up studies of this emerging model of evolutionary biology.

Hybridization among distantly related species: Examples from the polyploid genus Curcuma (Zingiberaceae).

Discerning relationships among species evolved by reticulate and/or polyploid evolution is not an easy task, although it is widely discussed. The economically important genus Curcuma (ca. 120 spp.; Zingiberaceae), broadly distributed in tropical SE Asia, is a particularly interesting example of a group of palaeopolyploid origin whose evolution is driven mainly by hybridization and polyploidization. Although a phylogeny and a new infrageneric classification of Curcuma, based on commonly used molecular markers (ITS and cpDNA), have recently been proposed, significant evolutionary questions remain unresolved. We applied a multilocus approach and a combination of modern analytical methods to this genus to distinguish causes of gene tree incongruence and to identify hybrids and their parental species. Five independent regions of nuclear DNA (DCS, GAPDH, GLOBOSA3, LEAFY, ITS) and four non-coding cpDNA regions (trnL-trnF, trnT-trnL, psbA-trnH and matK), analysed as a single locus, were employed to construct a species tree and hybrid species trees using (*)BEAST and STEM-hy. Detection of hybridogenous species in the dataset was also conducted using the posterior predictive checking approach as implemented in JML. The resulting species tree outlines the relationships among major evolutionary lineages within Curcuma, which were previously unresolved or which conflicted depending upon whether they were based on ITS or cpDNA markers. Moreover, by using the additional markers in tests of plausible topologies of hybrid species trees for C. vamana, C. candida, C. roscoeana and C. myanmarensis suggested by previous molecular and morphological evidence, we found strong evidence that all the species except C. candida are of subgeneric hybrid origin.

The origin of unique diversity in deglaciated areas: traces of Pleistocene processes in north-European endemics from the Galium pusillum polyploid complex (Rubiaceae).

The role of glacial oscillations in shaping plant diversity has been only rarely addressed in endemics of formerly glaciated areas. The Galium pusillum group represents a rare example of an ecologically diverse and ploidy-variable species complex that exhibits substantial diversity in deglaciated northern Europe. Using AFLP and plastid and nuclear DNA sequences of 67 populations from northern, central, and western Europe with known ecological preferences, we elucidate the evolutionary history of lineages restricted to deglaciated areas and identify the eco-geographic partitioning of their genetic variation. We reveal three distinct endemic northern lineages: (i) diploids from southern Sweden + the British Isles, (ii) tetraploids from southern Scandinavia and the British Isles that show signs of ancient hybridization between the first lineage and populations from unglaciated central Europe, and (iii) tetraploids from Iceland + central Norway. Available evidence supports a stepwise differentiation of these three lineages that started at least before the last glacial maximum by processes of genome duplication, interlineage hybridization and/or allopatric evolution in distinct periglacial refugia. We reject the hypothesis of more recent postglacial speciation. Ecological characteristics of the populations under study only partly reflect genetic variation and suggest broad niches of postglacial colonizers. Despite their largely allopatric modern distributions, the north-European lineages of the G. pusillum group do not show signs of rapid postglacial divergence, in contrast to most other northern endemics. Our study suggests that plants inhabiting deglaciated areas outside the Arctic may exhibit very different evolutionary histories compared with their more thoroughly investigated high-arctic counterparts.

Ancient hybridization and repetitive element proliferation in the evolutionary history of the monocot genus Amomum (Zingiberaceae).

Genome size variation is a crucial aspect of plant evolution, influenced by a complex interplay of factors. Repetitive elements, which are fundamental components of genomic architecture, often play a role in genome expansion by selectively amplifying specific repeat motifs. This study focuses on Amomum, a genus in the ginger family (Zingiberaceae), known for its 4.4-fold variation in genome size. Using a robust methodology involving PhyloNet reconstruction, RepeatExplorer clustering, and repeat similarity-based phylogenetic network construction, we investigated the repeatome composition, analyzed repeat dynamics, and identified potential hybridization events within the genus. Our analysis confirmed the presence of four major infrageneric clades (A-D) within Amomum, with clades A-C exclusively comprising diploid species (2n = 48) and clade D encompassing both diploid and tetraploid species (2n = 48 and 96). We observed an increase in the repeat content within the genus, ranging from 84% to 89%, compared to outgroup species with 75% of the repeatome. The SIRE lineage of the Ty1-Copia repeat superfamily was prevalent in most analyzed ingroup genomes. We identified significant difference in repeatome structure between the basal Amomum clades (A, B, C) and the most diverged clade D. Our investigation revealed evidence of ancient hybridization events within Amomum, coinciding with a substantial proliferation of multiple repeat groups. This finding supports the hypothesis that ancient hybridization is a driving force in the genomic evolution of Amomum. Furthermore, we contextualize our findings within the broader context of genome size variations and repeatome dynamics observed across major monocot lineages. This study enhances our understanding of evolutionary processes within monocots by highlighting the crucial roles of repetitive elements in shaping genome size and suggesting the mechanisms that drive these changes.

From Western Asia to the Mediterranean Basin: Diversification of the Widespread Euphorbia nicaeensis Alliance (Euphorbiaceae).

The Mediterranean Basin is an important biodiversity hotspot and one of the richest areas in the world in terms of plant diversity. Its flora parallels in several aspects that of the Eurasian steppes and the adjacent Irano-Turanian floristic region. The Euphorbia nicaeensis alliance spans this immense area from the western Mediterranean to Central Asia. Using an array of complementary methods, ranging from phylogenomic and phylogenetic data through relative genome size (RGS) estimation to morphometry, we explored relationships and biogeographic connections among taxa of this group. We identified the main evolutionary lineages, which mostly correspond to described taxa. However, despite the use of highly resolving Restriction Site Associated DNA (RAD) sequencing data, relationships among the main lineages remain ambiguous. This is likely due to hybridisation, lineage sorting triggered by rapid range expansion, and polyploidisation. The phylogenomic data identified cryptic diversity in the Mediterranean, which is also correlated with RGS and, partly, also, morphological divergence, rendering the description of a new species necessary. Biogeographic analyses suggest that Western Asia is the source area for the colonisation of the Mediterranean by this plant group and highlight the important contribution of the Irano-Turanian region to the high diversity in the Mediterranean Basin. The diversification of the E. nicaeensis alliance in the Mediterranean was triggered by vicariance in isolated Pleistocene refugia, morphological adaptation to divergent ecological conditions, and, to a lesser extent, by polyploidisation.

Stability in the South, Turbulence Toward the North: Evolutionary History of Aurinia saxatilis (Brassicaceae) Revealed by Phylogenomic and Climatic Modelling Data.

The Balkan Peninsula played an important role in the evolution of many Mediterranean plants and served as a major source for post-Pleistocene colonisation of central and northern Europe. Its complex geo-climatic history and environmental heterogeneity significantly influenced spatiotemporal diversification and resulted in intricate phylogeographic patterns. To explore the evolutionary dynamics and phylogeographic patterns within the widespread eastern Mediterranean and central European species Aurinia saxatilis, we used a combination of phylogenomic (restriction-site associated DNA sequencing, RADseq) and phylogenetic (sequences of the plastid marker ndhF) data as well as species distribution models generated for the present and the Last Glacial Maximum (LGM). The inferred phylogenies retrieved three main geographically distinct lineages. The southern lineage is restricted to the eastern Mediterranean, where it is distributed throughout the Aegean area, the southern Balkan Peninsula, and the southern Apennine Peninsula, and corresponds to the species main distribution area during the LGM. The eastern lineage extends from the eastern Balkan Peninsula over the Carpathians to central Europe, while the central lineage occupies the central Balkan Peninsula. Molecular dating places the divergence among all the three lineages to the early to middle Pleistocene, indicating their long-term independent evolutionary trajectories. Our data revealed an early divergence and stable in situ persistence of the southernmost, eastern Mediterranean lineage, whereas the mainland, south-east European lineages experienced more complex and turbulent evolutionary dynamics triggered by Pleistocene climatic oscillations. Our data also support the existence of multiple glacial refugia in southeast Europe and highlight the central Balkan Peninsula not only as a cradle of lineage diversifications but also as a source of lineage dispersal. Finally, the extant genetic variation within A. saxatilis is congruent with the taxonomic separation of peripatric A. saxatilis subsp. saxatilis and A. saxatilis subsp. orientalis, whereas the taxonomic status of A. saxatilis subsp. megalocarpa remains doubtful.

Congruent evolutionary responses of European steppe biota to late Quaternary climate change.

Quaternary climatic oscillations had a large impact on European biogeography. Alternation of cold and warm stages caused recurrent glaciations, massive vegetation shifts, and large-scale range alterations in many species. The Eurasian steppe biome and its grasslands are a noteworthy example; they underwent climate-driven, large-scale contractions during warm stages and expansions during cold stages. Here, we evaluate the impact of these range alterations on the late Quaternary demography of several phylogenetically distant plant and insect species, typical of the Eurasian steppes. We compare three explicit demographic hypotheses by applying an approach combining convolutional neural networks with approximate Bayesian computation. We identified congruent demographic responses of cold stage expansion and warm stage contraction across all species, but also species-specific effects. The demographic history of the Eurasian steppe biota reflects major paleoecological turning points in the late Quaternary and emphasizes the role of climate as a driving force underlying patterns of genetic variance on the biome level.

Integrative Study of Genotypic and Phenotypic Diversity in the Eurasian Orchid Genus Neotinea.

Knowledge of population variation across species' ranges is a prerequisite for correctly assessing the overall variability of any group of organisms and provides an invaluable basis for unraveling evolutionary history, optimizing taxonomy and devising effective conservation strategies. Here, we examine the genus Neotinea, which represents a relatively recently delimited monophyletic genus of orchids, for which a detailed study of its overall variability was lacking. We applied a suite of biosystematic methods, consisting of flow cytometry, multivariate and geometric morphometrics, and analysis of genomic SNP data, to identify phylogenetic lineages within the genus, to delineate phenotypic variation relevant to these lineages, and to identify potential cryptic taxa within lineages. We found clear differentiation into four major lineages corresponding to the groups usually recognized within the genus: Neotinea maculata as a distinct and separate taxon, the Neotinea lactea group comprising two Mediterranean taxa N. lactea and Neotinea conica, the Neotinea ustulata group comprising two phenologically distinct varieties, and the rather complex Neotinea tridentata group comprising two major lineages and various minor lineages of unclear taxonomic value. N. conica constitutes both a monophyletic group within N. lactea and a distinct phenotype within the genus and merits its proposed subspecies-level recognition. By contrast, the spring and summer flowering forms of N. ustulata (var. ustulata and var. aestivalis) were confirmed to be distinct only morphologically, not phylogenetically. The most complex pattern emerged in the N. tridentata group, which splits into two main clades, one containing lineages from the Balkans and eastern Mediterranean and the other consisting of plants from Central Europe and the central Mediterranean. These individual lineages differ in genome size and show moderate degrees of morphological divergence. The tetraploid Neotinea commutata is closely related to the N. tridentata group, but our evidence points to an auto- rather than an allopolyploid origin. Our broad methodological approach proved effective in recognizing cryptic lineages among the orchids, and we propose the joint analysis of flow cytometric data on genome size and endopolyploidy as a useful and beneficial marker for delineating orchid species with partial endoreplication.

Performance comparison of two reduced-representation based genome-wide marker-discovery strategies in a multi-taxon phylogeographic framework.

Multi-locus genetic data are pivotal in phylogenetics. Today, high-throughput sequencing (HTS) allows scientists to generate an unprecedented amount of such data from any organism. However, HTS is resource intense and may not be accessible to wide parts of the scientific community. In phylogeography, the use of HTS has concentrated on a few taxonomic groups, and the amount of data used to resolve a phylogeographic pattern often seems arbitrary. We explore the performance of two genetic marker sampling strategies and the effect of marker quantity in a comparative phylogeographic framework focusing on six species (arthropods and plants). The same analyses were applied to data inferred from amplified fragment length polymorphism fingerprinting (AFLP), a cheap, non-HTS based technique that is able to straightforwardly produce several hundred markers, and from restriction site associated DNA sequencing (RADseq), a more expensive, HTS-based technique that produces thousands of single nucleotide polymorphisms. We show that in four of six study species, AFLP leads to results comparable with those of RADseq. While we do not aim to contest the advantages of HTS techniques, we also show that AFLP is a robust technique to delimit evolutionary entities in both plants and animals. The demonstrated similarity of results from the two techniques also strengthens biological conclusions that were based on AFLP data in the past, an important finding given the wide utilization of AFLP over the last decades. We emphasize that whenever the delimitation of evolutionary entities is the central goal, as it is in many fields of biodiversity research, AFLP is still an adequate technique.

Evidence for Glacial Refugia of the Forest Understorey Species Helleborus niger (Ranunculaceae) in the Southern as Well as in the Northern Limestone Alps.

Glacial refugia of alpine and subnival biota have been intensively studied in the European Alps but the fate of forests and their understory species in that area remains largely unclear. In order to fill this gap, we aimed at disentangling the spatiotemporal diversification of disjunctly distributed black hellebore Helleborus niger (Ranunculaceae). We applied a set of phylogeographic analyses based on restriction-site associated DNA sequencing (RADseq) data and plastid DNA sequences to a range-wide sampling of populations. These analyses were supplemented with species distribution models generated for the present and the Last Glacial Maximum (LGM). We used exploratory analyses to delimit genomically coherent groups and then employed demographic modeling to reconstruct the history of these groups. We uncovered a deep split between two major genetic groups with western and eastern distribution within the Southern Limestone Alps, likely reflecting divergent evolution since the mid-Pleistocene in two glacial refugia situated along the unglaciated southern margin of the Alps. Long-term presence in the Southern Limestone Alps is also supported by high numbers of private alleles, elevated levels of nucleotide diversity and the species' modeled distribution at the LGM. The deep genetic divergence, however, is not reflected in leaf shape variation, suggesting that the morphological discrimination of genetically divergent entities within H. niger is questionable. At a shallower level, populations from the Northern Limestone Alps are differentiated from those in the Southern Limestone Alps in both RADseq and plastid DNA data sets, reflecting the North-South disjunction within the Eastern Alps. The underlying split was dated to ca. 0.1 mya, which is well before the LGM. In the same line, explicit tests of demographic models consistently rejected the hypothesis that the partial distribution area in the Northern Limestone Alps is the result of postglacial colonization. Taken together, our results strongly support that forest understory species such as H. niger have survived the LGM in refugia situated along the southern, but also along the northern or northeastern periphery of the Alps. Being a slow migrator, the species has likely survived repeated glacial-interglacial circles in distributional stasis while the composition of the tree canopy changed in the meanwhile.