Molecular biogeography of Europe: Pleistocene cycles and postglacial trends

Active around the year: Butterflies and moths adapt their life cycles to a warming world

Living in a warming world requires adaptations to altered annual temperature regimes. In Europe, spring is starting earlier, and the vegetation period is ending later in the year. These climatic changes are leading not only to shifts in distribution ranges of flora and fauna, but also to phenological shifts. Using long-term observation data of butterflies and moths collected during the past decades across northern Austria, we test for phenological shifts over time and changes in the number of generations. On average, Lepidoptera adults emerged earlier in the year and tended to extend their flight periods in autumn. Many species increased the annual number of generations. These changes were more pronounced at lower altitudes than at higher altitudes, leading to an altered phenological zonation. Our findings indicate that climate change does not only affect community composition but also the life history of insects. Increased activity and reproductive periods might alter Lepidoptera-host plant associations and food webs.

Patterns of speciation in a parapatric pair of Saturnia moths as revealed by target capture

The focus of this study has been to understand the evolutionary relationships and taxonomy of a widely distributed parapatric species pair of wild silk moths in Europe: Saturnia pavonia and Saturnia pavoniella (Lepidoptera: Saturniidae). To address species delimitation in these parapatric taxa, target enrichment and mtDNA sequencing was employed alongside phylogenetic, admixture, introgression, and species delimitation analyses. The dataset included individuals from both species close to and farther away from the contact zone as well as two hybrids generated in the lab. Nuclear markers strongly supported both S. pavonia and S. pavoniella as two distinct species, with hybrids forming a sister group to S. pavoniella. However, the Maximum Likelihood (ML) tree generated from mtDNA sequencing data presented a different picture, showing both taxa to be phylogenetically intermixed. This inconsistency is likely attributable to mitonuclear discordance, which can arise from biological factors (e.g., introgressive hybridization and/or incomplete lineage sorting). Our analyses indicate that past introgressions have taken place, but that there is no evidence to suggest an ongoing admixture between the two species, demonstrating that the taxa have reached full postzygotic reproductive isolation and hence represent two distinct biological species. Finally, we discuss our results from an evolutionary point of view taking into consideration the past climatic oscillations that have likely shaped the present dynamics between the two species. Overall, our study demonstrates the effectiveness of the target enrichment approach in resolving shallow phylogenetic relationships under complex evolutionary circumstances and that this approach is useful in establishing robust and well-informed taxonomic delimitations involving parapatric taxa.

Distribution patterns of Quercus ilex from the last interglacial period to the future by ecological niche modeling

The plants' geographic distribution is affected by natural or human-induced climate change. Numerous studies at both the global and regional levels currently focus on the potential changes in plant distribution areas. Ecological niche modeling can help predict the likely distribution of species according to environmental variables under different climate scenarios. In this study, we predicted the potential geographic distributions of Quercus ilex L. (holm oak), a keystone species of the Mediterranean ecosystem, for the Last Interglacial period (LIG: ~130 Ka), the Last Glacial Maximum (LGM: ~22 Ka), mid-Holocene (MH: ~6 Ka), and future climate scenarios (Representative Concentration Pathway (RCP) 4.5 and 8.5 scenarios) for 2050-2070 obtained from CCSM4 and MIROC-ESM global climate scenarios respectively. The models were produced with algorithms from the R-package "biomod2" and assessed by AUC of the receiver operating characteristic plot and true skill statistics. Aside from BIOCLIM (SRE), all model algorithms performed similarly and produced projections that are supported by good evaluation scores, although random forest (RF) slightly outperformed all the others. Additionally, distribution maps generated for the past period were validated through a comparison with pollen data acquired from the Neotoma Pollen Database. The results revealed that southern areas of the Mediterranean Basin, particularly coastal regions, served as long-term refugia for Q. ilex, which was supported by fossil pollen data. Furthermore, the models suggest long-term refugia role for Anatolia and we argue that Anatolia may have served as a founding population for the species. Future climate scenarios indicated that Q. ilex distribution varied by region, with some areas experiencing range contractions and others range expands. This study provides significant insights into the vulnerability of the Q. ilex to future climate change in the Mediterranean ecosystem and highlights the crucial role of Anatolia in the species' historical distribution.

Genomic analysis reveals complex population structure within the smooth newt, Lissotriton vulgaris, in Central Europe

Species with wide-range distributions usually display high genetic variation. This variation can be partly explained by historical lineages that were temporally isolated from each other and are back into secondary reproductive contact, and partly by local adaptations. The smooth newt (Lissotriton vulgaris) is one of the most widely distributed amphibians species across Eurasia and forms a species complex with a partially overlapping distribution and morphology. In the present study, we explored the population genomic structure of smooth newt lineages in the Carpathian Basin (CB) relying on single-nucleotide polymorphisms. Our dataset included new and previously published data to study the secondary contact zone between lineages in the CB and also tested for the barrier effect of rivers to gene flow between these lineages. We confirmed the presence of the South L. v. vulgaris Lineage distributed in Transdanubia and we provided new distribution records of L. v. ampelensis inhabiting the eastern territories of the CB. High genetic diversity of smooth newts was observed, especially in the North Hungarian Mountains and at the interfluves of the main rivers in the South with four distinct lineages of L. v. vulgaris and one lineage of L. v. ampelensis showing a low level of admixture with the spatially closest L. v. vulgaris lineage. Moreover, admixture detected at the interfluve of the main rivers (i.e. Danube and Tisza) suggested a secondary contact zone in the area. Finally, we found that the river Danube has a very weak effect on population divergence, while the river Tisza is a geographical barrier limiting gene flow between smooth newt lineages. As the range boundaries of L. v. vulgaris and L. v. ampelensis in the CB coincide with the river Tisza, our study underpins the influence of rivers in lineage diversification.

Range‐wide population genomics of the spongy moth, Lymantria dispar (Erebidae): Implications for biosurveillance, subspecies classification and phylogeography of a destructive moth

The spongy moth, Lymantria dispar, is an irruptive forest pest native to Eurasia where its range extends from coast to coast and overspills into northern Africa. Accidentally introduced from Europe in Massachusetts in 1868-1869, it is now established in North America where it is considered a highly destructive invasive pest. A fine-scale characterization of its population genetic structure would facilitate identification of source populations for specimens intercepted during ship inspections in North America and would enable mapping of introduction pathways to help prevent future incursions into novel environments. In addition, detailed knowledge of L. dispar's global population structure would provide new insight into the adequacy of its current subspecies classification system and its phylogeographic history. To address these issues, we generated >2000 genotyping-by-sequencing-derived SNPs from 1445 contemporary specimens sampled at 65 locations in 25 countries/3 continents. Using multiple analytical approaches, we identified eight subpopulations that could be further partitioned into 28 groups, achieving unprecedented resolution for this species' population structure. Although reconciliation between these groupings and the three currently recognized subspecies proved to be challenging, our genetic data confirmed circumscription of the japonica subspecies to Japan. However, the genetic cline observed across continental Eurasia, from L. dispar asiatica in East Asia to L. d. dispar in Western Europe, points to the absence of a sharp geographical boundary (e.g., the Ural Mountains) between these two subspecies, as suggested earlier. Importantly, moths from North America and the Caucasus/Middle East displayed high enough genetic distances from other populations to warrant their consideration as separate subspecies of L. dispar. Finally, in contrast with earlier mtDNA-based investigations that identified the Caucasus as L. dispar's place of origin, our analyses suggest continental East Asia as its evolutionary cradle, from where it spread to Central Asia and Europe, and to Japan through Korea.

Evidence of an additional centre of apple domestication in Iran, with contributions from the Caucasian crab apple Malus orientalis Uglitzk. to the cultivated apple gene pool

Divergence processes in crop-wild fruit tree complexes in pivotal regions for plant domestication such as the Caucasus and Iran remain little studied. We investigated anthropogenic and natural divergence processes in apples in these regions using 26 microsatellite markers amplified in 550 wild and cultivated samples. We found two genetically distinct cultivated populations in Iran that are differentiated from Malus domestica, the standard cultivated apple worldwide. Coalescent-based inferences showed that these two cultivated populations originated from specific domestication events of Malus orientalis in Iran. We found evidence of substantial wild-crop and crop-crop gene flow in the Caucasus and Iran, as has been described in apple in Europe. In addition, we identified seven genetically differentiated populations of wild apple (M. orientalis), not introgressed by the cultivated apple. Niche modelling combined with genetic diversity estimates indicated that these wild populations likely resulted from range changes during past glaciations. This study identifies Iran as a key region in the domestication of apple and M. orientalis as an additional contributor to the cultivated apple gene pool. Domestication of the apple tree therefore involved multiple origins of domestication in different geographic locations and substantial crop-wild hybridization, as found in other fruit trees. This study also highlights the impact of climate change on the natural divergence of a wild fruit tree and provides a starting point for apple conservation and breeding programmes in the Caucasus and Iran.

The evolution of the Aristolochia pallida complex (Aristolochiaceae) challenges traditional taxonomy and reflects large‐scale glacial refugia in the Mediterranean

The taxonomy of the Mediterranean Aristolochia pallida complex has been under debate since several decades with the following species currently recognized: A. pallida, A. lutea, A. nardiana, A. microstoma, A. merxmuelleri, A. croatica, and A. castellana. These taxa are distributed from Iberia to Turkey. To reconstruct phylogenetic and biogeographic patterns, we employed cpDNA sequence variation using both noncoding (intron and spacer) and protein‐coding regions (i.e., trnK intron, matK gene, and trnK‐psbA spacer). Our results show that the morphology‐based traditional taxonomy was not corroborated by our phylogenetic analyses. Aristolochia pallida, A. lutea, A. nardiana, and A. microstoma were not monophyletic. Instead, strong geographic signals were detected. Two major clades, one exclusively occurring in Greece and a second one of pan‐Mediterranean distribution, were found. Several subclades distributed in Greece, NW Turkey, Italy, as well as amphi‐Adriatic subclades, and a subgroup of southern France and Spain, were revealed. The distribution areas of these groups are in close vicinity to hypothesized glacial refugia areas in the Mediterranean. According to molecular clock analyses the diversification of this complex started around 3–3.3 my, before the onset of glaciation cycles, and the further evolution of and within major lineages falls into the Pleistocene. Based on these data, we conclude that the Aristolochia pallida alliance survived in different Mediterranean refugia rarely with low, but often with a high potential for range extension, and a high degree of morphological diversity.