Largely flat latitudinal life history clines in the dung fly Sepsis fulgens across Europe (Diptera: Sepsidae)

Geographic Variation in Genomic Signals of Admixture Between Two Closely Related European Sepsid Fly Species

The extent of interspecific gene flow and its consequences for the initiation, maintenance, and breakdown of species barriers in natural systems remain poorly understood. Interspecific gene flow by hybridization may weaken adaptive divergence, but can be overcome by selection against hybrids, which may ultimately promote reinforcement. An informative step towards understanding the role of gene flow during speciation is to describe patterns of past gene flow among extant species. We investigate signals of admixture between allopatric and sympatric populations of the two closely related European dung fly species Sepsis cynipsea and S. neocynipsea (Diptera: Sepsidae). Based on microsatellite genotypes, we first inferred a baseline demographic history using Approximate Bayesian Computation. We then used genomic data from pooled DNA of natural and laboratory populations to test for past interspecific gene flow based on allelic configurations discordant with the inferred population tree (ABBA-BABA test with D-statistic). Comparing the detected signals of gene flow with the contemporary geographic relationship among interspecific pairs of populations (sympatric vs. allopatric), we made two contrasting observations. At one site in the French Cevennes, we detected an excess of past interspecific gene flow, while at two sites in Switzerland we observed lower signals of past microsatellite genotypes gene flow among populations in sympatry compared to allopatric populations. These results suggest that the species boundaries between these two species depend on the past and/or present eco-geographic context in Europe, which indicates that there is no uniform link between contemporary geographic proximity and past interspecific gene flow in natural populations.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11692-023-09612-5.

On the Trail of Spatial Patterns of Genetic Variation

The accurate determination of the spatial trends on the variability of a species’ gene pool is essential to elucidate the underlying demographic-evolutionary events, thus helping to unravel the microevolutionary history of the population under study. Herein we present a new software called GenoCline, mainly addressed to detect genetic clines from allele, haplotype, and genome-wide data. This program package allows identifying the geographic orientation of clinal genetic variation through a system of iterative rotation of a virtual coordinate axis. Besides, GenoCline can perform complementary analyses to explore the potential origin of the genetic clines observed, including spatial autocorrelation, isolation by distance, centroid method, multidimensional scaling and Sammon projection. Among the advantages of this software is the ease in data entry and potential interconnection with other programs. Genetic and geographic data can be entered in spreadsheet table formatting (.xls), whereas genome-wide data can be imported in Eigensoft format. Genetic frequencies can also be exported in a format compatible with other programs dealing with population genetic and evolutionary biology analyses. All illustrations of results are saved in.svg format so that there will be high quality and easily editable vectorial graphs available for the researcher. Being implemented in Java, GenoCline is highly portable, thus working in different operating systems.

Using natural laboratories to study evolution to global warming: contrasting altitudinal, latitudinal, and urbanization gradients

Demonstrating the likelihood of evolution in response to global warming is important, yet challenging. We discuss how three spatial thermal gradients (latitudinal, altitudinal, and urbanization) can be used as natural laboratories to inform about the gradual thermal evolution of populations by applying a space-for-time substitution (SFTS) approach. We compare thermal variables and confounding non-thermal abiotic variables, methodological approaches and evolutionary aspects associated with each type of gradient. On the basis of an overview of recent insect studies, we show that a key assumption of SFTS, local thermal adaptation along these gradients, is often but not always met, requiring explicit validation. To increase realism when applying SFTS, we highlight the importance of integrating daily temperature fluctuations, multiple stressors and multiple interacting species. Finally, comparative studies, especially across gradient types, are important to provide more robust inferences of evolution under gradual global warming. Integrating these research directions will further strengthen the still underused, yet powerful SFTS approach to infer gradual evolution under global warming.

Usefulness and limitations of thermal performance curves in predicting ectotherm development under climatic variability

Thermal performance curves (TPCs) have been estimated in multiple ectotherm species to understand their thermal plasticity and adaptation and to predict the effect of global warming. However, TPCs are typically assessed under constant temperature regimes, so their reliability for predicting thermal responses in the wild where temperature fluctuates diurnally and seasonally remains poorly documented. Here, we use distant latitudinal populations of five species of sepsid flies (Diptera: Sepsidae) from the temperate region (Europe, North Africa, North America) to compare estimates derived from constant TPCs with observed development rate under fluctuating temperatures in laboratory and field conditions. TPCs changed across gradients in that flies originating from higher latitudes showed accelerated development at higher temperatures, an adaptive response. TPCs were then used to predict development rates observed under fluctuating temperatures; these predictions were relatively accurate in the laboratory but not the field. Interestingly, the precision of TPC predictions depended not only on the resolution of temperature data, with daily and overall temperature summing performing better than hourly temperature summing, but also on the frequency of temperatures falling below the estimated critical minimum temperature. Hourly temperature resolution most strongly underestimated actual development rates, because flies apparently either did not stop growing when temperatures dropped below this threshold, or they sped up their growth when the temperature rose again, thus most severely reflecting this error. We conclude that when flies do not encounter cold temperatures, TPC predictions based on constant temperatures can accurately reflect performance under fluctuating temperatures if adequately adjusted for nonlinearities, but when encountering cold temperatures, this method is more error-prone. Our study emphasizes the importance of the resolution of temperature data and cold temperatures in shaping thermal reaction norms.

Comparative reproductive dormancy differentiation in European black scavenger flies (Diptera: Sepsidae)

Seasonality is a key environmental factor that regularly promotes life history adaptation. Insects invading cold-temperate climates need to overwinter in a dormant state. We compared the role of temperature and photoperiod in dormancy induction in the laboratory, as well as winter survival and reproduction in the field and the laboratory, of 5 widespread European dung fly species (Diptera: Sepsidae) to investigate their extent of ecological differentiation and thermal adaptation. Unexpectedly, cold temperature is the primary environmental factor inducing winter dormancy, with short photoperiod playing an additional role mainly in species common at high altitudes and latitudes (Sepsis cynipsea, neocynipsea, fulgens), but not in those species also thriving in southern Europe (thoracica, punctum). All species hibernate as adults rather than juveniles. S. thoracica had very low adult winter survivorship under both (benign) laboratory and (harsh) field conditions, suggesting flexible quiescence rather than genetically fixed winter diapause, restricting their distribution towards the pole. All other species appear well suited for surviving cold, Nordic winters. Females born early in the season reproduce before winter while late-born females reproduce after winter, fulgens transitioning earliest before winter and thoracica and punctum latest; a bet-hedging strategy of reproduction during both seasons occurs rarely but is possible physiologically. Fertility patterns indicate that females can store sperm over winter. Winter dormancy induction mechanisms of European sepsids are congruent with their geographic distribution, co-defining their thermal niches. Flexible adult winter quiescence appears the easiest route for insects spreading towards the poles to evolve the necessary overwinter survival.

Does thermal plasticity align with local adaptation? An interspecific comparison of wing morphology in sepsid flies

Although genetic and plastic responses are sometimes considered as unrelated processes, their phenotypic effects may often align because genetic adaptation is expected to mirror phenotypic plasticity if adaptive, but run counter to it when maladaptive. Because the magnitude and direction of this alignment has further consequences for both the tempo and mode of adaptation, they are relevant for predicting an organisms' reaction to environmental change. To better understand the interplay between phenotypic plasticity and genetic change in mediating adaptive phenotypic variation to climate variability, we here quantified genetic latitudinal variation and thermal plasticity in wing loading and wing shape in two closely related and widespread sepsid flies. Common garden rearing of 16 geographical populations reared across multiple temperatures revealed that wing loading decreases with latitude in both species. This pattern could be driven by selection for increased dispersal capacity in the cold. However, although allometry, sexual dimorphism, thermal plasticity and latitudinal differentiation in wing shape all show similar patterns in the two species, the relationship between the plastic and genetic responses differed between them. Although latitudinal differentiation (south to north) mirrored thermal plasticity (hot to cold) in Sepsis punctum, there was no relationship in Sepsis fulgens. While this suggests that thermal plasticity may have helped to mediate local adaptation in S. punctum, it also demonstrates that genetic wing shape differentiation and its relation to thermal plasticity may be complex and idiosyncratic, even among ecologically similar and closely related species. Hence, genetic responses can, but do not necessarily, align with phenotypic plasticity induced by changing environmental selection pressures.