Heredity

Heat Shock Protein 90 (HSP90) functions as an evolutionary capacitor, allowing populations to store cryptic genetic variation that can be released under stress. While former studies have described the release of morphological variation, its behavioural consequences remain unexplored. In the red flour beetle, Tribolium castaneum, HSP90 inhibition released a phenotype with much smaller, less defined eyes that confers fitness benefits in continuous light and was subsequently assimilated. We hypothesized that altered eye morphology affects light perception and thereby changes light-dependent behaviours. To test whether phenotypes released via evolutionary capacitance can beneficially alter behaviour, we examined locomotor activity rhythm entrainment to light-dark cycles as well as individual and group light choice behaviour. Males of the reduced-eye phenotype exhibited a diminished startle response to sudden light exposure in locomotor activity assays. We also found reduced negative phototaxis in groups of beetles with reduced eyes. This modified behaviour, indicating reduced light sensitivity, may stem from impaired light perception caused by altered eye morphology. Lower light sensitivity could be beneficial under stressful environmental conditions by promoting the exploration of alternative niches. Therefore, this study provides the first evidence for potentially beneficial behavioural changes in a HSP90-released phenotype, reinforcing HSP90s role as an evolutionary capacitor.

Alleles with opposing effects on fitness characters are said to exhibit selectional antagonistic pleiotropy (broadly construed so that effects are not necessarily confined to the same individual). A number of theoretical investigations considered the case where a pair of alleles at a locus influences two fitness components and derived the conditions giving rise to stable polymorphism under various assumptions about the mode of trait-interaction. Strikingly, many of these analyses concluded that the potential for maintaining polymorphism is strongly constrained by the joint influence of two factors: (1) the prevalence of weak selection coefficients over coefficients of large magnitude, and (2) the absence of beneficial dominance reversals (where the deleterious effects of each allele are partially or completely masked in the heterozygous genotype). Consequently, the conclusion that selective polymorphism is unlikely to be maintained by intralocus mechanisms of antagonistic pleiotropy has achieved widespread acceptance. Here we argue that such conclusions do not apply to any of the following models of antagonism: (i) additive trait-interaction, (ii) multiplicative trait-interaction, (iii) bivoltine selection, (iv) soft selection, (v) hard selection, and (vi) sexual antagonism. We demonstrate that the parameter space giving rise to stable allelic variation is quite large throughout, and moreover, the plenitude of suitable parameters neither depends on the strength of selection nor requires dominance reversal. Dominance coefficients associated with stringent conditions for stable polymorphism are shown to be atypical as compared to all feasible parameters, and best regarded as an outcome of adherence to a special relation: dominance with a constant magnitude and direction, which includes the case of additive allelic effects at a locus. Properties of single-locus equilibria (heterozygosity, allele frequency differentiation) are investigated, as well as the contribution of dominance schemes to the genetic variance in fitness characters in populations at multilocus linkage equilibrium. Author summaryAllelic variants at a locus with opposing effects on multiple fitness components (antagonistic fitness pleiotropy) have long been appreciated as a possible source of balancing selection. The prevalence of polymorphism owing to this form of natural selection, however, has been doubted on theoretical grounds due to the fact that standard assumptions of genetic models (namely, constant magnitudes for the dominance coefficients) are hardly conducive to the maintenance of polymorphism. The major exception to this conclusion lies with schemes that exhibit dominance reversal (where the direction of dominance for antagonistic alleles flips across fitness components). Here we conduct a geometric analysis of the space of polymorphism-promoting dominance parameters and conclude that the conditions for maintaining balanced alleles is unrestrictive, with non-reversals playing an underappreciated role.

The extent to which phenotypic evolution can be constrained by genetic correlations is an important question in evolutionary biology. To address this question, biological invasions are opportune models where derived, invasive populations can be compared to their extant ancestors, allowing to track the evolution of genetic correlations from the ancestor, throughout the invasion process. In this paper, we focused on the worldwide invasion of Drosophila suzukii (Matsumara, 1931), and investigated the evolution of the genetic covariance matrix G of wing shape between ancestral native, and derived invasive populations. Leveraging demographic history resolved by population genetics approaches, we tested whether G remained stable during the invasion. Using a multivariate QST-FST approach, we further tested whether or not the observed phenotypic divergence in wing shape aligned with a neutral scenario of evolution. Our results show moderate yet significant quantitative genetic differentiation of wing shape among D. suzukii populations and a relative stability in the structure of G, presenting a roughly spherical shape but slightly different volumes. These characteristics likely reflect the demographic history of populations and suggest a low level of genetic constraint on wing shape evolution. The divergence between populations was greater than expected under a purely neutral model of evolution, compatible with an effect of divergent selection among them. Overall, our study suggests that selection and drift, but not ancestral genetic constraints, affected the early stages of wing shape evolution during D. suzukii invasion.

Estimation of the effective size (Ne) of large populations with a continuous distribution across wide geographic areas and limited dispersal of individuals has been elusive so far. Estimates of the contemporary Ne from genetic markers for such large, structured populations, typically of plant and marine species, tend to be strongly biased downwards, which has led to question their relevance. Here we show that a recently proposed estimation method of Ne from linkage disequilibrium between markers, which accounts for population structure, yields estimates of metapopulation Ne when the sampling area is sufficiently large. The method is applied to empirical data of maritime pine (Pinus pinaster Aiton). While previous estimates of Ne in pine populations were of the order of a few hundred individuals, we show that estimates of the metapopulation Ne can reach values of the order of tens of thousands of individuals. This result is especially relevant from a conservation point of view, as populations with Ne lower than 500 individuals are considered to be under the risk of extinction.

Deciphering the genetic architecture of complex quantitative phenotypes remains challenging in quantitative genetics. These traits not only depend of multiple genetic factors but are also established over time and environments. Although quantitative genetics has investigated the genetic determinism of phenotypic plasticity in contrasted environmental conditions, the time related phenotypic plasticity has received less attention. Here we proposed a multivariate Bayesian framework, the Bayesian Varying Coefficient Model, designed for analysing the genetic architecture of the time related phenotypic plasticity by a multilocus approach. We applied the BVCM to time series phenotypes measured at various time scales (daily, monthly, yearly) across a diverse set of biological species. We included in this study: yeast (Saccharomyces cerevisiae), fungi (Fusarium graminearum), eucalyptus (Eucalyptus urophylla x E. grandis), and sweet cherry tree (Prunus avium). The BVCM results were compared with those obtained with a known genome-wide association method carried out time by time. For all species and traits, the BVCM was able to detect the major QTL identified by marker-trait association methods and revealed additional genetic regions of weak effect. It also increased the phenotypic variance explained for most of the phenotypes considered. It revealed dynamic QTLs with transitory, increasing or decreasing effects over time. By considering both the temporal and genetic multivariate structures in a single statistical model, we increased our understanding of the genetic architecture of complex traits notably by reducing the issue of missing heritability. More broadly, this work raises the foundation for extended applications in functional genomics, evolutionary ecology, and crop breeding programs, in which time-related phenotypic plasticity remains crucial for predicting and selecting key quantitative complex traits. Key messageBy capturing the genetic factors influencing the time related phenotypic plasticity, our approach contributes to a deeper understanding of the dynamic nature of genotype-phenotype relationships.

The ectoparasitic mite Varroa destructor is the most significant parasite of the Western honey bee (Apis mellifera) and a major driver of colony losses worldwide. Although extensive research has focused on behavioral and physiological mechanisms of host resistance, comparatively little attention has been paid to potential phenotypic responses of the parasite itself. Here we investigated body size variation in Varroa destructor associated with varroa-resistant and non-resistant honey bee colonies across four European countries. We quantified the dorsal shield area of adult female mites from multiple colonies differing in the honey bee colonies resistance status, using standardized digital image analysis. Across geographically distant non-resistant populations, mite body size was remarkably consistent, with a median dorsal shield area of 1.47 mm2. In contrast, mites originating from varroa-resistant colonies were consistently smaller, with a median dorsal shield area of 1.37 mm2, representing an approximately 6.8% reduction in body size. This pattern was reproducible across different geographical areas, honey bee genetic backgrounds and beekeeping practices. The striking stability of mite body size in non-resistant populations contrasted with the consistent reduction observed in mites associated with resistant hosts, suggesting a host associated shift in parasite phenotype. Because body size in arthropods integrates developmental conditions, nutritional availability and resource allocation, the observed pattern may reflect altered developmental environments and selective pressures imposed by resistant hosts. Our results show a consistent morphological shift in this globally important parasite associated with resistant hosts and suggest that dorsal shield size in Varroa could serve as a new selection marker for varroa-resistant honey bee colonies.

Trees accumulate somatic mutations throughout their long lifespan, resulting in genetic mosaicism among branches. While recent genomic studies quantified these mutations, they were largely limited to describing static patterns of variation. In this study, we developed a mathematical model to infer the dynamic processes of somatic mutation accumulation from snapshot genomic data obtained from four tropical trees (Dipterocarpaceae), which dominate tropical rain forests in Southeast Asia. Our model focus on genetic differences between shoot apical meristems (SAMs) at branch tips and explicitly incorporate stem cell dynamics within SAMs during shoot elongation and branching, enabling us to quantify somatic genetic drift arising from stem cell lineage replacement. By comparing model predictions with empirical data from Dipterocarpaceae trees, we estimated key parameters governing stem cell dynamics and somatic mutation rates. Our results indicate that both shoot elongation and branching involve replacement of stem cell lineages, leading to a moderate degree of somatic genetic drift. Accounting for stem cell dynamics resulted in slightly lower mutation rate estimates than previous approaches that ignored these processes. Using the estimated parameters, we further performed stochastic simulations to predict patterns of somatic mutations, including features not directly observed in the sampled trees, such as occasional deviations of somatic mutation phylogenies from physical architecture. Together, our modeling framework provides insights into how genetic mosaicism is shaped within tropical trees and reveals the stem cell dynamics underlying their long-term growth and accumulation of somatic mutations. (236 words) Highlights- We built mathematical models to predict the genetic differences between branch tips by somatic mutations. - The model considers the varying dynamics of stem cells in shoot meristem during shoot elongation and branching. - We compared the model prediction with empirical data from tropical trees, Dipterocarpaceae, and estimated the dynamics of stem cells and mutation rate. - Somatic mutation dynamics were shaped by somatic genetic drift arising from stem cell lineage replacement during shoot elongation and branching. - Accounting for stem cell dynamics led to slightly smaller estimates of mutation rates compared with previous estimates that ignored the dynamics. - Our models offer insights into how genetic variability is shaped in the tropical trees and the stem cell dynamics underlying their long-term growth.

Population genetic theory predicts that natural selection will be more efficient in large than small populations because genetic drift reduces the efficiency of selection in small populations. Small populations adapting to new environments can also be expected to evolve higher recombination rates to facilitate adaptation as well as to dissociate and purge harmful mutations. We tested these hypotheses (1) by investigating differences in the strength of association between nucleotide diversity ({pi}) and recombination rate across the genomes of nine-spined sticklebacks (Pungitius pungitius) from four small freshwater (mean Ne {approx} 2 578) and four large marine (mean Ne = 86 742) populations, as well as (2) by comparing recombination rates between small and large populations using population specific linkage maps. We found the predicted positive correlation of{pi} with recombination rate from all but the smallest freshwater populations, suggesting prevalent linked selection even after accounting for variation in GC/CpG content, and gene density. Mean recombination rates did not differ between freshwater and marine populations, except that the smallest Ne freshwater population exhibited significantly elevated recombination rate. GWAS analyses suggested a polygenic basis for recombination rates. These results suggest an important role for linked selection in reducing{pi} in low recombination regions especially in large populations. Moreover, as predicted by theory, at least one of the small freshwater populations appears to have evolved a higher recombination rate than its marine ancestors.

Genotype-by-environment interaction (GEI) has been studied to identify environment-stable/favorable genotypes. The GEI simulation could help refine the inference by incorporating tangible factors such as genomic and environmental information. The Bayesian additive main effect and multiplicative interaction (Bayesian AMMI) model captures the genotype-specific responses across environments, reflecting directional relationships between genotypes and environments. Thus, we propose a Bayesian AMMI-based GEI simulation framework that utilizes high-throughput environmental covariance matrices to generate GEI effects with interpretable directional structure. To demonstrate the proposed approach, two simulated phenotypes were assessed under four levels of GEI variance. In the first simulation (Sim1), GEI effects were sampled from a multivariate normal distribution defined by the GEI matrix. In the second simulation (Sim2), GEI effects were generated by extending Sim1 with the Bayesian AMMI model. In both simulations, increasing GEI variance resulted in lower correlations of phenotypes across environments and stronger genotype-specific sensitivity to environmental variation. Across five cross-validation designs, models accounting for GEI consistently outperformed one that did not, with prediction accuracy generally decreasing as GEI variance increased. Clear distinctions between the two simulated phenotypes were evident from biplot analyses: Sim2 successfully captured environmental relatedness and genotype-specific responses, whereas such structure was absent in Sim1. These results demonstrate that the proposed Bayesian AMMI-based GEI simulation framework enables interpretable visualization of GEI and supports genomic selection strategies under complex environmental conditions.

Climate warming alters the thermal environment experienced by ectotherms, whose physiological performance and fitness are constrained by temperature. Early life stages are often the temperature-sensitive phases of the life cycle, with potential consequences for population persistence, particularly in freshwater stenotherms such as the Arctic charr (Salvelinus alpinus). The persistence of populations will partly depend on the adaptive potential of critical life stages to environmental changes. In this study, we used a common garden approach to compare the response and phenotypic plasticity of four charr populations to warmer conditions. These populations inhabit thermally contrasted lakes and differ in origin (native/introduced) and management history. We reared embryos at either an optimal (5{degrees}C) temperature for larval development or a warmer but realistic (8.5 {degrees}C) temperature. We tested adaptive divergence among populations in four traits (survival, incubation duration, body length and yolk sac volume), using Qst - Fst comparisons. We report negative effects of temperature on body size, survival and earlier hatching. Thermal reaction norms differed among populations, indicating adaptive divergence. Contrary to expectations, populations originating from warmer environments did not consistently exhibit higher trait values under elevated temperatures. In contrast, the unmanaged and colder high-altitude population exhibited higher survival rates and lower yolk reserves for a given size under heat stress than the other populations. Our results suggested that evolutionary trajectories specific to each population are shaped by factors related to the populations history, including introductions, demographic fluctuations and long-term repopulation practices, which can jointly influence the potential for adaptation to heat stress.

Many simultaneous hermaphrodites use selfing for reproductive assurance only when outcrossing opportunities are limited, owing to inbreeding depression in selfed progeny. However, scenarios that enforce substantial selfing (such as during recolonisation) can rapidly select for a high selfing propensity, a shift in mating system that is expected to eliminate both inbreeding depression and the delayed reproductive onset under selfing that is typically associated with it. We tested these predictions in the flatworm Macrostomum hystrix, using a line derived from an outcrossing population that had been subjected to enforced selfing for multiple generations followed by several years of relaxed selection. As predicted, isolated individual forced to self and individuals with constant partner access (i.e. outcrossing opportunities) did not differ in reproductive onset nor in inbreeding depression estimated through offspring survival. However, a third treatment group that provided intermittent partner access (to allow outcrossing but minimise potential competition effects) showed a different pattern: no inbreeding depression in offspring but a substantially accelerated reproductive onset. Whilst our results thus support the effective purging of inbreeding depression and increased selfing propensity under enforced selfing, we suggest that cues of an unstable social or physical environment nevertheless exert a major influence on reproductive timing.

Epimutations are changes in chromatin modifications, such as DNA methylation or histone modifications. Some of these epigenetic changes can be inherited for several generations, and they potentially contribute to evolutionary processes. Estimates of epimutation rates now exists in a few species, but the presence and function of epigenetic marks are not conserved across different species. To understand the properties of epimutations in fungi, we performed a mutation accumulation experiment with the filamentous fungus Neurospora crassa and investigated spontaneous changes in DNA methylation and trimethylation of lysine 9 on histone H3 (H3K9me3) in the mutation accumulation lines. We observed that centromeric regions are hotspots of spontaneous DNA methylation changes in N. crassa. In these hotspot regions, DNA methylation changes were transmitted across mitoses, but changes occurring in euchromatin were not maintained. The rate of DNA methylation changes was around 30 000 fold faster than the genetic mutation rate. We did not observe spontaneous changes in H3K9me3 that were transmitted across mitoses. Our results show that while spontaneous epimutations occur in this species, they occur predominantly in gene poor heterochromatic regions, so their impact for evolutionary adaptation may be limited.

Premise of studyUnder increasingly frequent pollinator-limited environments, plants need to rely on modes of reproductive assurance such as selfing and cloning. However, few studies investigate the interplay between selfing and cloning in plants that can do both. Here, we explore mechanisms determining the relative expression of selfing and cloning throughout the European distribution of the wild woodland strawberry (Fragaria vesca) under a pollinator-free environment. MethodsWe established an outdoor common garden with 121 woodland strawberry genotypes from across Europe and excluded them from pollinators. For each genotype, we recorded reproductive traits and performed hand-pollination treatments. Key ResultsWe found a weak trade-off between cloning and selfing, driven by increased seed and fruit provisioning rather than flower production. The capacity to autonomously self-fertilize was determined by the lateral proximity of the anthers to the pistils (lateral herkogamy), but not by early inbreeding depression. Genotypes sampled at lower latitudes and altitudes were better at self-fertilizing and had smaller petals. The propensity to clone increased towards the east, where genotypes also had smaller petals, particularly at higher latitudes. ConclusionAt the species level, we detected a trade-off between the propensity for clonal reproduction and the capacity for self-fertilization. At a continental scale, the capacity to self-fertilize varied along a north-south gradient, whereas clonal propensity varied along an east-west gradient. Our results suggest that these two modes of reproductive assurance may compensate for reduced pollinator attractiveness (smaller petals) in regions where each mode is most strongly expressed.

Linkage disequilibrium (LD) and haplotype block structure govern the resolution and utility of genomic selection, marker-assisted selection, and genome-wide association studies (GWAS) in livestock. We performed a comprehensive genome-wide characterization of LD decay, haplotype block architecture, and population diversity across all 24 autosomes in Nili-Ravi buffalo (Bubalus bubalis; n = 85), using 43,543 post-quality-control SNPs. Mean genome-wide r2 was 0.124 (median 0.074) and mean D was 0.540 (median 0.481), with LD half-decay at {approx}70 kb. A total of 133 haplotype blocks encompassing 721 SNPs were identified (Gabriel et al., 2002). Haploview analysis of nine chromosomes harbouring bTB resistance candidate genes revealed contrasting selection signatures: directional selection at innate immune loci (IFNG, TLR1; H < 0.55) versus balancing selection at adaptive immune loci (BoLA-DRB3, SP110; H > 1.0). Critically, BBU15 Block 3 (28.6 kb; OR52E5/NCR1 locus, 47.16 Mb) showed a genome-wide significant integrated haplotype score (iHS; -log1 0 p = 5.408), directly co-localising with the published bTB susceptibility QTL (Bermingham et al., 2014). The TAA haplotype (frequency 53.3%) at this block represents a candidate resistance-associated haplotype for marker-assisted selection. These findings provide essential parameters for SNP panel design and bTB resistance breeding in South Asian buffalo.

Life history traits are often correlated, creating trade-offs that may impede the response to natural selection and be responsible for the evolution of senescence. These trade-offs may arise through pleiotropic effects, which can affect the response to selection in ways that resemble intra-locus sexual antagonism. Despite these hypothesized relationships, we lack clear connections between pleiotropy, sexual antagonism, and the evolution of life histories. Empirical tests for inter-sexual differences in life-history traits, including sex-specific aging, can be used to evaluate hypotheses about how pleiotropy and sexual conflict affect evolutionary trade-offs. To those ends, we measured lifespan, development time, and body size in Drosophila pseudoobscura males and females, each of which carried one of six third chromosome inversion genotypes. Temperature affected lifespan and development more than any other factor; higher temperatures increased mortality rate, decreased lifespan, and accelerated development. However, we also observed sex differences in mortality rates and development times that depended on genotype and temperature. Notably, temperature elevated the initial mortality rate across all flies, yet increasing temperatures reduced the rate of aging in some genotype-sex combinations. Similarly, direct effects of genotype on mortality rate and development time depended greatly on sex and temperature, but there was no genotype effect on body size. Despite these context-dependent genotype effects on life history traits, we failed to identify any correlations that would serve as clear evidence for sexual conflict or trade-offs. Our results therefore suggest that either historical conflicts have been resolved or any conflicts that may exist do not result in the correlations predicted by existing models.

Honeybee breeding plans are relatively recent in most countries. In France, diverse small-scale breeding groups are emerging. Beekeepers are highly diverse in their motivations, farm productions and services, practices and management techniques. Yet, little is known about what beekeepers would consider as relevant breeding goals in the design of breeding plans. We therefore conducted an online survey answered by about 250 French beekeepers, mostly professionals, to assess their perceived importance of including 20 pre-defined traits in breeding goals and to identify how beekeeping profiles might influence these priorities. Respondents rated each trait as essential, useful, or useless, and indicated if they wished useful or essential traits to be genetically improved or merely maintained at their current level. Results indicated a strong preference for multi-trait selection, with a median of 13 traits considered useful or essential. Honey yield, disease resistance, swarming tendency, gentleness, and summer feed autonomy, emerged as the main traits of interest with about 90% of beekeepers finding them at least useful. About 40% or more only wished to maintain these traits at their current level rather than to directionally improve them. A major exception to this was disease resistance, that 75% wanted to improve. Bees genetic background influenced the most the importance attributed to breeding goal traits, while other beekeeping profile characteristics only had a marginal effect on breeding goal trait priorities. Some poorly studied traits, such as summer and winter feed autonomy, winter diapause, and longevity, were considered at least useful in a breeding goal by over 70% of beekeepers. Future research is needed to explore possible selection criteria for these traits and estimate the potential for their genetic improvement. ImplicationsOur survey shows that French beekeepers wish to improve or maintain through selective breeding usual colony production and behavioral traits, but also colony resilience, especially disease resistance and feed autonomy. However, trait priorities differ depending on the genetic background of the bees used. This knowledge is essential for designing breeding programs that truly match beekeeper needs and for identifying which traits deserve research attention. In France, beekeepers are increasingly starting breeding efforts to adapt their bees to current conditions, facing growing pressures from climate change, diseases, invasive species, and pesticides. Well-designed breeding programs can support sustainable beekeeping and essential pollination services.

Reproductive Character Displacement (RCD) often occurs when species with mating-related polymorphism come into secondary contact, leading to divergence in reproductive traits. Ischnura elegans and Ischnura graellsii have formed two independent hybrid zones in Spain where reinforcement has strengthened a mechanical barrier, and RCD has shaped mating-related structures, although reinforcement is asymmetric only in gynochrome females. This study examines the link between asymmetric reinforcement and asymmetric RCD. Using geometric morphometrics, we analyze prothorax shape and size in both female morphs and males, and male caudal appendages, to assess morphological divergence, determine whether gynochrome females show stronger divergence, and test for morphological covariation between male traits involved in the tandem position. Our results reveal consistent patterns of size and shape variation across species and zones: in I. elegans, androchromes are larger and resemble males in size, with clear shape differentiation between female morphs that diminishes in hybrid zones. In contrast, I. graellsii shows less consistent size differences between males and morphs, and weaker shape differentiation. Our results confirm RCD in prothorax shape in I. elegans females from both hybrid zones, but reveal that RCD in prothorax size is asymmetric, occurring only in gynochrome females from the NC hybrid zone. We also detected RCD in the prothorax shape of I. elegans males from the NC hybrid zone, extending previous evidence of RCD in male caudal appendages, while morphological covariation between male cerci and the prothorax was limited to size in I. elegans. Together, these findings illustrate how hybridization may generate morph-specific patterns of reproductive divergence.

Understanding relatedness in sharks is challenging due to uncertainty in distributions, low population densities and difficulties in sampling across life stages. In Fiji, bull sharks (Carcharhinus leucas), with an effective population size estimate of [~]258, aggregate at the Shark Reef Marine Reserve (SRMR), but gravid females disperse at the end of the year to give birth in adjacent rivers. Questions remain regarding reproductive connectivity, female returns across years, and kinship structure. Using population genomics on 296 bull sharks across age classes (neonates, young-of-the-year, juveniles, and adults) collected over a decade at the SRMR and in three adjacent rivers, we assessed familial connections. Direct genetic links, including first- and second-degree relationships, connected SRMR adults with young age classes in the Navua and Rewa rivers, providing evidence of reproductive connectivity. Within rivers, genetic similarities across cohorts revealed reproductive philopatry. Remarkably, several individuals sampled years apart were assigned to the same sire-dam pairs, indicating repeated pairings across breeding seasons. However, the few related links detected between the SRMR and the rivers may reflect incomplete sampling. Altogether, bull shark reproduction in Fiji seems influenced by reproductive philopatry and repeated pairings, suggesting added complexity in their reproductive behaviour.

Genetic drift in natural populations reduces the efficacy of selection, promoting the fixation of deleterious recessive alleles with consequences for maladaptation and population persistence. Heterosis, or increased F1 fitness relative to the parental mean, has been proposed as a tool for investigating the role of drift on genetic variation in fitness, but its genetic basis and environmental dependence remain unclear in natural populations. We used heterozygous near-isogenic lines (NILs) derived from a cross between locally adapted Arabidopsis thaliana ecotypes to assess how specific genomic regions influence heterosis. Cumulative fitness, estimated as fruits per seedling, was evaluated in a greenhouse and two simulated native environments. F1s showed strong heterosis in the greenhouse and one simulated environment. Non-additive effects in heterozygous NILs were highly environment- and background-dependent, varying in magnitude and sign, and no NIL had consistently effects across environments. The relative fitness of NILs was not correlated with gene number or genomic load in the introgressed regions. Small heterozygous regions often had large effects, indicating that complementation of mildly deleterious alleles alone does not fully explain heterosis and suggesting that overdominance or pseudo-overdominance may play a role. Evidence of epistasis was also observed, including outbreeding depression in some NILs, likely due to negative additive-by-dominance interactions. Summed effects of NILs often exceeded the fitness increase of the F1 suggesting dominance-by-dominance epistasis, but the direction of these epistatic effects depended on both genetic background and environment. Our results demonstrate that F1 fitness reflects both positive dominance and different epistatic interactions that are environment- and background-dependent.

Genetic evidence for fluctuating selection has begun to accumulate for different species over the past few decades, especially for the Drosophila genus where studies have reported hundreds of loci undergoing putatively adaptive oscillations across successive seasons. However, most theoretical and simulation studies of fluctuating selection have relied on abstract or weakly parameterized models, making it difficult to assess their relevance for natural populations. In this study, we simulate multilocus seasonally fluctuating selection under a recently developed model and examine its effect on the variance effective population size (Ne) at a genome-wide scale. By recapitulating genomic, demographic, and evolutionary parameters from natural Drosophila populations in our simulations, we were able to reproduce allele frequency oscillations reported in recent studies and show that these lead to [~]50% genome-wide reductions in Ne. We also demonstrate that Ne reductions are well predicted by the maximum frequency amplitude among all adaptively fluctuating loci, and that the frequency amplitudes are largely determined by the number of adaptively fluctuating loci and the strength of their epistatic interactions. Our results demonstrate that fluctuating selection can substantially reduce effective population size and underscore the importance of temporally variable selection in shaping genome-wide patterns of variation beyond classical models. Article SummaryGenetic studies of fluctuating selection in natural populations have grown steadily over the past decade, with reports suggesting that hundreds of loci undergo adaptive oscillations over seasonal timescales in cosmopolitan Drosophila populations. By simulating seasonally fluctuating selection under a recently developed model and ecological scenarios informed by published studies, the authors show that this mode of selection can reduce effective population size by [~]50%, with the magnitude of the reduction correlated with the locus exhibiting the largest allele frequency fluctuations. These findings highlight fluctuating selection as an important factor shaping genome-wide patterns of genetic variation and effective population size.