Additivity of the genetic load revealed to be a product of both synergistic and antagonistic epistasis in Drosophila.

Additivity of the genetic load revealed to be a product of both synergistic and antagonistic epistasis in Drosophila. | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Article Additivity of the genetic load revealed to be a product of both synergistic and antagonistic epistasis in Drosophila . Imran Sayyed, Halle McNamara, Troy Day, Adam Chippindale This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-7910627/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Epistasis has been theoretically implicated in several major evolutionary processes, including the evolution of sex and speciation, but empirical evidence for its impact is still lacking. We tested for epistatic interactions using a full-sib inbreeding design in Drosophila melanogaster . We examined the effects of deleterious mutations on two sequential components of life-history: pupal productivity and pupa-to-adult viability. We observed an accelerating decline in pupal productivity with increasing inbreeding coefficient, consistent with synergistic epistasis, whereas pupa-to-adult viability exhibited a decelerating decline, consistent with antagonistic epistasis. This decoupling led to an overall linear decline in the production of adult offspring with increasing inbreeding. Our findings suggest that epistatic interactions can be highly trait-specific, underscoring the importance of analyzing individual, physiologically interconnected, fitness traits to understand the role of epistasis in evolutionary processes. Biological sciences/Evolution/Evolutionary theory Biological sciences/Evolution/Experimental evolution Figures Figure 1 Figure 2 Figure 3 Figure 4 Introduction Epistasis, defined as a non-additive interaction between two or more genes, is hypothesized to play a significant role in regulating long-range linkage disequilibrium, thereby shaping pleiotropic associations and the organization of conserved genes (Singhal et al., 2023). Such genetic mechanisms are purported to shape the evolution of several fundamental biological processes, including speciation, sexual reproduction and perhaps, aging, suggesting epistasis may be a major driving force in these phenomena (Charlesworth, 1998; de Visser & Elena, 2007; Kondrashov, 1988; Wachter et al., 2013; Wagner et al., 1994; Weinreich et al., 2005; Zhang et al., 2024). Epistatic interactions can be categorized as either synergistic or antagonistic. Synergistic epistasis occurs when the combined deleterious effect of interacting alleles on fitness is greater than the sum of their individual effects, relative to an additive or multiplicative model. Conversely, antagonistic epistasis describes a scenario where the combined deleterious effect is less than that predicted from the individual effects (Weinreich et al., 2005; Zhang et al., 2024). Such non-additivity can be observed relative to the additive model of fitness or multiplicative model when fitness is analyzed on a logarithmic scale. However, empirical investigations of epistasis have proven challenging and often produce equivocal results, as explained below. Two approaches have been employed to investigate patterns of fitness decline potentially associated with epistatic interactions among deleterious mutations. The first approach involves mutation accumulation (MA) experiments, which allow for the investigation of fitness decline patterns due to the increasing effects of spontaneously arising deleterious mutations, either under relaxed selection or using mutagens (Dickinson, 2008; Elena & Lenski, 1997; Kondrashov, 1994). These newly accumulated mutations may interact, causing deviations from additive fitness decline patterns. However, the evidence from mutation accumulation experiments regarding fitness decline potentially related to epistasis among deleterious alleles remains varied in quality and equivocal in outcomes (Rivero et al., 2003; Visser et al., 1997; West et al., 1998; Whitlock & Bourguet, 2000). Another approach to studying fitness decline patterns involves inbreeding depression, which can be explained by either the partial dominance hypothesis or the overdominance hypothesis. Under the partial dominance hypothesis, fitness decline occurs due to the increased homozygosity of pre-existing recessive deleterious mutations maintained within a population. In contrast, the overdominance hypothesis posits that fitness reduction arises from a loss of heterozygote advantage. Given that heterozygote advantage is typically observed in a small percentage of loci in an equilibrium population, it is generally assumed that the majority of fitness decline due to inbreeding occurs under the partial dominance hypothesis (Charlesworth et al., 1991; Crow, 1999; Hedrick, 2012; Roff, 2002). Studies of inbreeding depression have reported conflicting evidence, generally showing either no epistasis or antagonistic epistasis for fitness (Ávila et al., 2006; Rand et al., 2006; Sharp & Agrawal, 2016; Sohail et al., 2017; Szafraniec et al., 2003). Some of these discrepancies may be attributable to experimental biases, such as the repeated analysis of the same lines throughout the inbreeding process, sampling bias from the loss of lines due to lethal recessive deleterious mutations, and the effects of purging. However, a recent study by Domínguez-García et al., (2019) provided empirical evidence of synergistic epistasis for productivity using a full-sib inbred design in both laboratory and natural populations of Drosophila melanogaster . While Domínguez-García et al. (2019) addressed some of the extant challenges, certain experimental limitations persisted in their experiments, such as only analysing surviving lines which potentially creates a bias against synergistic epistasis among strong deleterious mutations. Another frequently overlooked factor in studies of inbreeding and the pattern of fitness decline is the direct influence of the inbreeding coefficient on fitness components separately. A critical assumption often made is that the accumulation of deleterious mutations, and consequently their homozygous expression, increases at a consistent linear rate across a variety of fitness components affected by these mutations. Under such an assumption, the inbreeding coefficient is expected to affect various life history traits equally leading to a linear decline in the overall fitness. While a linear relationship between increasing homozygosity of recessive deleterious mutations and their phenotypic impact is often assumed across all fitness components, it is plausible that different life history traits may exhibit varying rates or non-linear patterns of decline in response to increasing inbreeding coefficients (Charlesworth, 1998; Falconer, 1996; Willis, 1993). Investigating these potential differential responses in fitness components is crucial for a comprehensive understanding of genetic load and epistatic interactions. This conflicted but developing literature inspired our study, in which we investigated the effect of increasing inbreeding coefficient on several composite traits which determine the fitness of fruit flies. We used a full-sib inbreeding design but replicated parental lines to create ten copies of each initiating genotype. We investigated the patterns of decline in pupal productivity, pupa-to-adult viability, and overall offspring productivity with increasing inbreeding coefficients. Our results revealed an accelerated decline in pupal productivity but a decelerated decline in pupa-to-adult viability. This culminated in an overall linear decline in offspring productivity with increasing inbreeding coefficients. These findings highlight the complexity of studying epistasis among deleterious mutations affecting fitness, especially given that fitness is a highly composite trait. Furthermore, our results underscore the fact that distinct epistatic interactions can influence specific life-history traits without necessarily altering the expected additive effects on overall fitness. Materials and Methods Study system: Our experimental design utilized a laboratory population of Drosophila melanogaster designated Control (C 1 ) that is descended from the Ives (IV) base population. IV was established from a wild caught sample of 400 flies (1:1 sex ratio) in Amherst, Massachusetts in 1975 (Rose & Charlesworth, 1981a, 1981b). The C 1 population has a varied demographic history but is outbred and had been maintained under standardized environmental conditions since 1980. These conditions include discrete generations, a population size of approximately 2000 individuals per generation, cultivation on a banana-agar-killed yeast food medium, a moderate larval density (80-100 larvae per vial), a constant temperature of 25°C, and a 12-hour light:12-hour dark photoperiod. This population had been maintained for over 150 generations on a two-week life cycle. With more than 400 total generations of controlled laboratory maintenance, the C 1 population is expected to be at, or close to, an evolutionary equilibrium for standing genetic variation. This characteristic makes it particularly well-suited for estimating mutational load at an equilibrium and measuring fitness under conditions of evolutionary history by assessing various life-history traits. Full-sib Inbreeding: To establish the inbred lines, 30 unmated male and female pairs were isolated from the C 1 population, with each pair initiating a distinct line within a separate vial. In the second generation, 10 unmated male and female pairs were isolated from the offspring of each established inbred line to create 10 replicate 'sub-lines' (300 lines in total). For the subsequent five generations, each sub-line was maintained by randomly collecting a single pair of unmated male and female offspring; a sub-line was terminated if no viable adult pairs were produced in a generation. This ‘sub-line’ design helps us to account for the possible extinction of lines in the analysis due to strong deleterious effect arising from synergistic epistasis leading to a bias against accelerated decline in surviving lines. In total, these experimental lines underwent six generations of inbreeding, during which fitness decline was assessed. Assessing fitness: The pattern of fitness decline due to increasing inbreeding was assessed by measuring several key life-history traits. These included pupal productivity as the number of pupae produced by an inbred pair, the pupae to adult viability during pupal stage, and offspring productivity measured as the total number of adults eclosed for each surviving sub-line. Concurrently, these same traits were measured in the outbred control population (C 1 ) by sampling 50 male and female pairs each generation, synchronized with the measurements from the inbred lines. For any sub-line that did not survive in a given generation, the zero values were still treated as data in analyses. Analysis: All analyses were performed using R version 4.4.2 (R Core Team, 2024). Since pupal productivity and offspring productivity are composite traits that account for the reproductive fitness of parents, number of eggs laid, hatched and pupated successfully, the expected inbreeding coefficient ( ICf ) for pupal productivity and offspring productivity was calculated and assigned for each generation of inbreeding by averaging the inbreeding coefficients of the mother and the offspring based on the previous work by Domínguez-García et al. (2019). This resulted in expected ICf values of 0, 0.125, 0.313, 0.438, 0.547, 0.633 for the generations 1 to 6, respectively, for the entire genome in a full-sib line. For pupa-to-adult viability and sub-line survival, the expected inbreeding coefficient was calculated based solely on the offspring's inbreeding coefficient ( IC ), yielding values of 0, 0.25, 0.375, 0.5, 0.594, 0.671. Using the lme4 package in R, we fitted linear and quadratic hurdle models to the measured traits against the expected inbreeding coefficient. In our regression models, inbred lines were included as a random effect to account for stochasticity inherent to each homozygous line. Model comparisons between linear and quadratic fits for each trait were performed using Akaike Information Criterion corrected for small sample sizes (AICc), Bayesian Information Criterion (BIC), and log-likelihood comparisons to determine whether the observed decline was linear or non-linear with increasing inbreeding coefficient. Results 1. Survival of sublines: The number of surviving sub-lines per inbred line exhibited an accelerating decline of 58% over 6 generations with increasing inbreeding coefficient ( ICf ), as inferred from comparisons between quadratic and linear model fits (Figure 1, Table 1). This non-linear decline in sub-line survival suggests a significant effect of synergistic epistasis among lethal mutations or synergistic interactions between recessive deleterious mutations affecting sterility or pre-copulatory mortality. 2. Pupal productivity: Overall pupal productivity declined by 56.02% in a non-linear fashion with increasing inbreeding coefficient (Figure 2, Table 1). Among the 30 inbred lines, 24 lines showed a significant accelerated decline in pupal productivity, providing strong evidence for synergistic epistasis involving mutations that affect egg-to-pupae viability and development. 3. Pupa-to-adult viability: Overall pupa-to-adult viability, representing survival during the pupal stage of pre-adult development, declined by 18.65% in a non-linear fashion with increasing inbreeding coefficient ( IC ) (Figure 3, Table 1). Out of the 30 inbred lines, 22 lines exhibited a significant deceleration in pupa-to-adult viability, suggesting strong evidence for antagonistic epistasis among mutations affecting this developmental stage. 4. Number of adults eclosed : The overall number of adults eclosing, representing survival through the entire pre-adult developmental stage, declined by 64.37% in a linear fashion with increasing inbreeding coefficient (Figure 4, Table 1). Among the 30 inbred lines, 27 lines showed a significant linear decline in the number of adults eclosed, which appears counter-intuitive given the results for pupal productivity and pupa-to-adult viability. Instead, this linear pattern can be interpreted as an outcome of the combined effects of accelerated decline in pupal productivity and de-accelerated decline in pupa-to-adult viability, resulting in an overall linear reduction in the number of adults eclosed with increasing ICf . Discussion Experiments aimed at deducing the effects of increasing mutation load have yielded a range of outcomes and are equivocal, as a literature, on the question of epistasis for fitness effects. These studies have varied in scale and design, with limits to inference of any one approach. For example, mutation accumulation (MA) experiments have not provided consistent evidence for epistasis among deleterious mutations (Ávila et al., 2006; Dickinson, 2008; Elena & Lenski, 1997; Kibota & Lynch, 1996; Visser et al., 1997). These investigations typically focus on newly arising spontaneous mutations that would be subject to natural selection and subsequent purging from the population; such spontaneous mutations may not be the ones segregating naturally and having the greatest impact in populations. Consequently, an inbreeding approach offers a potentially valuable alternative approach for examining epistatic interactions among pre-existing recessive deleterious mutations that have accumulated in a population under mutation-selection balance in natural populations (Charlesworth & Hughes, 1996; Domínguez-García et al., 2019; Swindell & Bouzat, 2006). Our study population has been maintained for over 400 generations in the lab at large population size, a duration that is expected to result in something approaching an equilibrium state, minimizing the impact of purging and balancing selection on the analysis of epistatic interactions. Few inbreeding studies have been directed at the question of epistasis, however several do suggest the presence of epistatic interactions, both synergistic and antagonistic, affecting various fitness-related traits (Domínguez-García et al., 2019; Rosa et al., 2005; Willis, 1993). Given that fitness is a complex, polygenic trait, and researchers often use proxies for total fitness, it is important to dissect epistatic interactions within its constituent traits individually, where feasible. Our current investigation was inspired by Domínguez-García et al. (2019) who suggested the existence of synergistic epistasis for fitness, evidenced by an accelerated decline in pupal productivity in D. melanogaster laboratory and natural populations. Consistent with their findings, our results also indicated an accelerating decline in a similar metric of pupal productivity. However, we also observed a deceleration in the effects of inbreeding on pupa-to-adult viability, which ultimately led to a linear decline in the total number of eclosed adults (total fitness) with increasing inbreeding. Why then is the impact of inbreeding different for pupal productivity and pupal viability? Pupal productivity is a highly composite character, encompassing almost all aspects of total fitness: adult fertility in the parental generation and larval viability in the subsequent generation. Potentially the sheer complexity and number of different ways that productivity could fail through infertility or poor juvenile competition and development, with performance deficits in one phase potentially compounding in the next, creates the conditions for synergism across mutations. In contrast, pupal survival is likely to be an intensely canalized trait, with the metamorphic animal executing a highly choreographed developmental phase highly dependent on the interactions of rare alleles with environmental stress (Bendixsen et al., 2017; Borne et al., 2021; Collet et al., 2023). In summary, although we are unable to pinpoint the driving elements from the current analysis, we speculate that the expression of deleterious mutations in one life stage could have a potential cascading or ‘synergistic’ effect in later stages leading to larger than expected fitness decline, without directly interacting with other deleterious mutations. This hypothesis involves condition-dependence, in which erosion of condition below thresholds results in catastrophic failures and morbidity. This mechanism may be complimentary to a study from Sohail et al. (2017), who detected rare loss-of-function alleles exhibiting synergistic epistasis under purifying selection in both human and Drosophila populations. On the other hand, it explains why purifying selection during mutation-selection balance could promote the accumulation of recessive deleterious mutations exhibiting antagonistic epistasis during the pupal stage in D. melanogaster . This hypothesis is consistent with evidence of a higher rate of non-synonymous substitutions relative to mutation rate (ω) in coding regions during the pupal stage of D. melanogaster compared to ω in the coding regions involved in pupal productivity (Coronado-Zamora et al., 2019). Our findings demonstrate that epistatic interactions among deleterious mutations may be strongly trait-specific, to the point of shifting from synergistic to antagonistic depending upon the particular fitness components being assessed. This means that, at least in our study population, even measurements of fertility and survival that capture nearly all aspects of total fitness may misrepresent that parameter, which is of greatest importance from the theoretical standpoint. While synergistic epistasis for total fitness has been implicated in broadly important evolutionary phenomena, such as the evolution of sexual reproduction, our results do not support such a role. However, our study does raise intriguing questions about why different life stages and combinations of fitness characters would yield opposing results regarding epistasis. Hypotheses regarding the relationship between purifying selection, epistasis, and the ratio of non-synonymous to synonymous substitutions (ω) can serve as a starting point to investigate such trait-specific epistatic interactions, but their combination in sequence poses special challenges, particularly for highly composite traits. Overall, understanding the genetic and physiological underpinnings of non-additivity will allow a greater insight into the evolutionary architecture of fitness and into the genetic load of populations. Declarations Conflict of interest: The authors declare they have no conflict of interest. Data archiving: To be uploaded on https://datadryad.org/ Acknowledgments: Funding for the work was from NSERC Discovery grants to AKC & TD. We thank all the undergraduate students involved in the population maintenance. We thank Harshavardhan Thygarajan, Mindy Baroody, Avery Want, Josh Kowal for their efforts in the lab. Author contribution statement: IS was responsible for designing the protocol, analysing data, creating figures and tables and interpreting results. IS and HM contributed equally to conducting the research and collecting the data. IS wrote the first draft of the paper. IS and AKC contributed to the updating of reference list and writing the report. IS, TD and AKC contributed to the final draft and revisions. References Ávila, V., Chavarrías, D., Sánchez, E., Manrique, A., López-Fanjul, C., & García-Dorado, A. (2006). Increase of the spontaneous mutation rate in a long-term experiment with Drosophila melanogaster . Genetics , 173 (1), 267–277. https://doi.org/10.1534/genetics.106.056200 Bendixsen, D. P., Østman, B., & Hayden, E. J. (2017). Negative epistasis in experimental RNA fitness landscapes. Journal of Molecular Evolution , 85 (5), 159–168. https://doi.org/10.1007/s00239-017-9817-5 Borne, F., Kulathinal, R. J., & Courtier-Orgogozo, V. (2021). Glue genes are subjected to diverse selective forces during Drosophila development. Genome Biology and Evolution , 13 (12), evab248. https://doi.org/10.1093/gbe/evab248 Charlesworth, B. (1998). The effect of synergistic epistasis on the inbreeding load. Genetics Research , 71 (1), 85–89. https://doi.org/10.1017/S0016672398003140 Charlesworth, B., & Hughes, K. A. (1996). Age-specific inbreeding depression and components of genetic variance in relation to the evolution of senescence. Proceedings of the National Academy of Sciences , 93 (12), 6140–6145. https://doi.org/10.1073/pnas.93.12.6140 Charlesworth, B., Morgan, M. T., & Charlesworth, D. (1991). Multilocus models of inbreeding depression with synergistic selection and partial self-fertilization. Genetics Research , 57 (2), 177–194. https://doi.org/10.1017/S0016672300029256 Collet, J. M., Nidelet, S., & Fellous, S. (2023). Genetic independence between traits separated by metamorphosis is widespread but varies with biological function. Proceedings of the Royal Society B: Biological Sciences , 290 (2010), 20231784. https://doi.org/10.1098/rspb.2023.1784 Coronado-Zamora, M., Salvador-Martínez, I., Castellano, D., Barbadilla, A., & Salazar-Ciudad, I. (2019). Adaptation and conservation throughout the Drosophila melanogaster life-cycle. Genome Biology and Evolution , 11 (5), 1463–1482. https://doi.org/10.1093/gbe/evz086 Crow, J. F. (1999). Dominance and overdominance. In Genetics and Exploitation of Heterosis in Crops (pp. 49–58). https://doi.org/10.2134/1999.geneticsandexploitation.c5 de Visser, J. A. G. M., & Elena, S. F. (2007). The evolution of sex: Empirical insights into the roles of epistasis and drift. Nature Reviews Genetics , 8 (2), Article 2. https://doi.org/10.1038/nrg1985 Dickinson, W. J. (2008). Synergistic fitness interactions and a high frequency of beneficial changes among mutations accumulated under relaxed selection in Saccharomyces cerevisiae . Genetics , 178 (3), 1571–1578. https://doi.org/10.1534/genetics.107.080853 Domínguez-García, S., García, C., Quesada, H., & Caballero, A. (2019). Accelerated inbreeding depression suggests synergistic epistasis for deleterious mutations in Drosophila melanogaster . Heredity , 123 (6), 709–722. https://doi.org/10.1038/s41437-019-0263-6 Elena, S. F., & Lenski, R. E. (1997). Test of synergistic interactions among deleterious mutations in bacteria. Nature , 390 (6658), 395–398. https://doi.org/10.1038/37108 Falconer, D. S. (1996). Introduction to quantitative genetics 4th edition. Pearson Education , ISBN 8131727408, 9788131727409. Hedrick, P. W. (2012). What is the evidence for heterozygote advantage selection? Trends in Ecology & Evolution , 27 (12), 698–704. https://doi.org/10.1016/j.tree.2012.08.012 Kibota, T. T., & Lynch, M. (1996). Estimate of the genomic mutation rate deleterious to overall fitness in E. coli . Nature , 381 (6584), 694–696. https://doi.org/10.1038/381694a0 Kondrashov, A. S. (1988). Deleterious mutations and the evolution of sexual reproduction. Nature , 336 (6198), Article 6198. https://doi.org/10.1038/336435a0 Kondrashov, A. S. (1994). Muller’s ratchet under epistatic selection. Genetics , 136 (4), 1469–1473. https://doi.org/10.1093/genetics/136.4.1469 R Core Team. (2024). R: The R Project for Statistical Computing . https://www.r-project.org/ Rand, D. M., Fry, A., & Sheldahl, L. (2006). Nuclear–mitochondrial epistasis and Drosophila aging: introgression of Drosophila simulans mtdna modifies longevity in D. melanogaster nuclear backgrounds. Genetics , 172 (1), 329–341. https://doi.org/10.1534/genetics.105.046698 Rivero, A., Balloux, F., & West, S. A. (2003). Testing for epistasis between deleterious mutations in a parasitoid wasp. Evolution , 57 (7), 1698–1703. https://doi.org/10.1111/j.0014-3820.2003.tb00375.x Roff, D. A. (2002). Inbreeding depression: Tests of the overdominance and partial dominance hypotheses. Evolution , 56 (4), 768–775. https://doi.org/10.1111/j.0014-3820.2002.tb01387.x Rosa, J. M., Camacho, S., & García‐dorado, A. (2005). A measure of the within‐chromosome synergistic epistasis for Drosophila viability. Journal of Evolutionary Biology , 18 (4), 1130–1137. https://doi.org/10.1111/j.1420-9101.2005.00892.x Rose, M. R., & Charlesworth, B. (1981a). Genetics of life history in Drosophila melanogaster . I. Sib analysis of adult females. Genetics , 97 (1), 173–186. https://doi.org/10.1093/genetics/97.1.173 Rose, M. R., & Charlesworth, B. (1981b). Genetics of life history in Drosophila melanogaster . II. Exploratory selection experiments. Genetics , 97 (1), 187–196. https://doi.org/10.1093/genetics/97.1.187 Sharp, N. P., & Agrawal, A. F. (2016). The decline in fitness with inbreeding: Evidence for negative dominance‐by‐dominance epistasis in Drosophila melanogaster . Journal of Evolutionary Biology , 29 (4), 857–864. https://doi.org/10.1111/jeb.12815 Singhal, P., Veturi, Y., Dudek, S. M., Lucas, A., Frase, A., van Steen, K., Schrodi, S. J., Fasel, et al. (2023). Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets. American Journal of Human Genetics , 110 (4), 575–591. https://doi.org/10.1016/j.ajhg.2023.03.007 Sohail, M., Vakhrusheva, O. A., Sul, J. H., Pulit, S. L., Francioli, L. C., Genome of the Netherlands Consortium, Alzheimer’s Disease Neuroimaging Initiative, van den Berg, L. H., Veldink, J. H., de Bakker, P. I. W., Bazykin, G. A., Kondrashov, A. S., & Sunyaev, S. R. (2017). Negative selection in humans and fruit flies involves synergistic epistasis. Science , 356 (6337), 539–542. https://doi.org/10.1126/science.aah5238 Swindell, W. R., & Bouzat, J. L. (2006). Inbreeding Depression and Male Survivorship in Drosophila : Implications for Senescence Theory. Genetics , 172 (1), 317–327. https://doi.org/10.1534/genetics.105.045740 Szafraniec, K., Wloch, D. M., Sliwa, P., Borts, R. H., & Korona, R. (2003). Small fitness effects and weak genetic interactions between deleterious mutations in heterozygous loci of the yeast Saccharomyces cerevisiae . Genetics Research , 82 (1), 19–31. https://doi.org/10.1017/S001667230300630X Visser, J. A. G. M. de, Hoekstra, R. F., & Ende, H. van den. (1997). Test of interaction between genetic markers that affect fitness in Aspergillus niger . Evolution , 51 (5), 1499–1505. https://doi.org/10.1111/j.1558-5646.1997.tb01473.x Wachter, K. W., Evans, S. N., & Steinsaltz, D. (2013). The age-specific force of natural selection and biodemographic walls of death. Proceedings of the National Academy of Sciences of the United States of America , 110 (25), 10141–10146. https://doi.org/10.1073/pnas.1306656110 Wagner, A., Wagner, G. P., & Similion, P. (1994). Epistasis can facilitate the evolution of reproductive isolation by peak shifts: A two-locus two-allele model. Genetics , 138 (2), 533–545. https://doi.org/10.1093/genetics/138.2.533 Weinreich, D. M., Watson, R. A., & Chao, L. (2005). Perspective: Sign epistasis and genetic constraint on evolutionary trajectories. Evolution; International Journal of Organic Evolution , 59 (6), 1165–1174. PMID: 16050094. West, S. A., Peters, A. D., & Barton, N. H. (1998). Testing for epistasis between deleterious mutations. Genetics , 149 (1), 435–444. https://doi.org/10.1093/genetics/149.1.435 Whitlock, M. C., & Bourguet, D. (2000). Factors affecting the genetic load in Drosophila : synergistic epistasis and correlations among fitness components. Evolution , 54 (5), 1654–1660. https://doi.org/10.1111/j.0014-3820.2000.tb00709.x Willis, J. H. (1993). Effects of different levels of inbreeding on fitness components in mimulus guttatus . Evolution , 47 (3), 864–876. https://doi.org/10.1111/j.1558-5646.1993.tb01240.x Zhang, J., Chen, F., & Li, X. (2024). Mechanistic causes of sign epistasis and its applications. Frontiers in Genetics , 15 , 1366917. https://doi.org/10.3389/fgene.2024.1366917 Table Table 1 is available in the Supplementary Files section. Additional Declarations There is no duality of interest Supplementary Files Tables.xlsx Table 1 Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-7910627","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Article","associatedPublications":[],"authors":[{"id":539981803,"identity":"93f2a9c2-20de-4def-b86e-359d7f046578","order_by":0,"name":"Imran Sayyed","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA9klEQVRIiWNgGAWjYBACNmSOxAeStUjOINlKaR5iVPGJnX34uYChTt7gRvLB2zZlh+X4pdufSTDU2OF2mHS6sfQMhsOGG26kJVvnnDtsLDnnjJkEw7FkPFrSQO45wLjhRo6ZdG7b4UQgg02CsYEZnxbm3zwMdfYbbuR/k7ZsO1y//0b6M6CWenxa2IC2MIMNl2ZsO5xgIJFgBtRyGK8Wax6Dw8kzzzwztuw5l24440aOsUXCseM4tcjPTmO+zVNRZ9t3PPnhjR9l1vL8M9If3vhQU41TCwQYMDAoXEhgQMRsAgENEOv6DzCgJYZRMApGwSgYBRAAAHJ3TSL7tzltAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0009-0004-5762-3142","institution":"Queen's University","correspondingAuthor":true,"prefix":"","firstName":"Imran","middleName":"","lastName":"Sayyed","suffix":""},{"id":539981804,"identity":"beff1e14-4343-4cfd-9652-2b8c2a79eed0","order_by":1,"name":"Halle McNamara","email":"","orcid":"","institution":"Queen's University","correspondingAuthor":false,"prefix":"","firstName":"Halle","middleName":"","lastName":"McNamara","suffix":""},{"id":539981805,"identity":"bc369e48-2826-4c5a-9312-15142c56f293","order_by":2,"name":"Troy Day","email":"","orcid":"","institution":"Queen's University, Ontario","correspondingAuthor":false,"prefix":"","firstName":"Troy","middleName":"","lastName":"Day","suffix":""},{"id":539981806,"identity":"892964e9-c07f-44a1-be1d-40246e3322a9","order_by":3,"name":"Adam Chippindale","email":"","orcid":"https://orcid.org/0000-0003-3384-9556","institution":"Queens University","correspondingAuthor":false,"prefix":"","firstName":"Adam","middleName":"","lastName":"Chippindale","suffix":""}],"badges":[],"createdAt":"2025-10-21 05:10:29","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-7910627/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-7910627/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":96244840,"identity":"fab8701d-0966-4ef8-88cb-ca81127e3a3d","added_by":"auto","created_at":"2025-11-19 07:19:23","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":24918,"visible":true,"origin":"","legend":"\u003cp\u003eBoxplot of number of sublines surviving at a given IC. The black solid line represents the best linear fit, while dotted non-linear line represents quadratic fit for the decline in subline survival with IC.\u003c/p\u003e","description":"","filename":"Binder11.png","url":"https://assets-eu.researchsquare.com/files/rs-7910627/v1/b1ec52906578547198746655.png"},{"id":95949994,"identity":"e2a26cdb-a2c5-4bf2-aa93-016b27d61ce4","added_by":"auto","created_at":"2025-11-14 19:02:34","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":76481,"visible":true,"origin":"","legend":"\u003cp\u003eQuadratic regression for pupal productivity (on log-scale) with increasing average inbreeding coefficient (ICf). The relative scale is with respect to the outbred lines. Each point indicates pupal productivity for each inbred sub-line at a given ICf.\u003c/p\u003e","description":"","filename":"Binder12.png","url":"https://assets-eu.researchsquare.com/files/rs-7910627/v1/ad9080861a4618dac9bd7662.png"},{"id":96244142,"identity":"a7db4765-cf9f-471e-b39e-8af1ff80a6f6","added_by":"auto","created_at":"2025-11-19 07:17:52","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":56371,"visible":true,"origin":"","legend":"\u003cp\u003eQuadratic regression for pupae-to-adult viability (on log-scale) with increasing inbreeding coefficient (IC). The relative scale is with respect to the outbred lines. Each point indicates pupae-to-adult viability for each surviving sub-line at a given IC, as this viability ratio can not be calculated for extinct lines.\u003c/p\u003e","description":"","filename":"Binder13.png","url":"https://assets-eu.researchsquare.com/files/rs-7910627/v1/9b6524e580c9ea227b3098b4.png"},{"id":95949997,"identity":"e918c59f-53c3-417e-8452-20fe7120074b","added_by":"auto","created_at":"2025-11-14 19:02:34","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":76449,"visible":true,"origin":"","legend":"\u003cp\u003eLinear regression for number of adults eclosed (on log-scale) with increasing average inbreeding coefficient (ICf). The relative scale is with respect to the outbred lines. Each point indicates number of adults eclosed for each inbred sub-line at a given ICf.\u003c/p\u003e","description":"","filename":"Binder14.png","url":"https://assets-eu.researchsquare.com/files/rs-7910627/v1/f10cf34e321aa1a8aafd45b1.png"},{"id":97895695,"identity":"b18a6bda-aefc-4dd3-9f9d-3df1ec6c49b7","added_by":"auto","created_at":"2025-12-10 15:34:43","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":628714,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-7910627/v1/eb6c5eb5-a7cd-4fb9-abb1-b3fb8a6e9965.pdf"},{"id":96244491,"identity":"200bf85a-6c88-422a-a97c-3bec7b5f3659","added_by":"auto","created_at":"2025-11-19 07:18:35","extension":"xlsx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":10615,"visible":true,"origin":"","legend":"Table 1","description":"","filename":"Tables.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-7910627/v1/b5d8584e96e2253a64a6e407.xlsx"}],"financialInterests":"There is no duality of interest","formattedTitle":"Additivity of the genetic load revealed to be a product of both synergistic and antagonistic epistasis in \u003ci\u003eDrosophila\u003c/i\u003e.","fulltext":[{"header":"Introduction","content":"\u003cp\u003eEpistasis, defined as a non-additive interaction between two or more genes, is hypothesized to play a significant role in regulating long-range linkage disequilibrium, thereby shaping pleiotropic associations and the organization of conserved genes (Singhal et al., 2023). Such genetic mechanisms are purported to shape the evolution of several fundamental biological processes, including speciation, sexual reproduction and perhaps, aging, suggesting epistasis may be a major driving force in these phenomena (Charlesworth, 1998; de Visser \u0026amp; Elena, 2007; Kondrashov, 1988; Wachter et al., 2013; Wagner et al., 1994; Weinreich et al., 2005; Zhang et al., 2024). Epistatic interactions can be categorized as either synergistic or antagonistic. Synergistic epistasis occurs when the combined deleterious effect of interacting alleles on fitness is greater than the sum of their individual effects, relative to an additive or multiplicative model. Conversely, antagonistic epistasis describes a scenario where the combined deleterious effect is less than that predicted from the individual effects (Weinreich et al., 2005; Zhang et al., 2024). Such non-additivity can be observed relative to the additive model of fitness or multiplicative model when fitness is analyzed on a logarithmic scale. However, empirical investigations of epistasis have proven challenging and often produce equivocal results, as explained below.\u003c/p\u003e\n\u003cp\u003eTwo approaches have been employed to investigate patterns of fitness decline potentially associated with epistatic interactions among deleterious mutations. The first approach involves mutation accumulation (MA) experiments, which allow for the investigation of fitness decline patterns due to the increasing effects of spontaneously arising deleterious mutations, either under relaxed selection or using mutagens (Dickinson, 2008; Elena \u0026amp; Lenski, 1997; Kondrashov, 1994). These newly accumulated mutations may interact, causing deviations from additive fitness decline patterns. However, the evidence from mutation accumulation experiments regarding fitness decline potentially related to epistasis among deleterious alleles remains varied in quality and equivocal in outcomes (Rivero et al., 2003; Visser et al., 1997; West et al., 1998; Whitlock \u0026amp; Bourguet, 2000).\u003c/p\u003e\n\u003cp\u003eAnother approach to studying fitness decline patterns involves inbreeding depression, which can be explained by either the partial dominance hypothesis or the overdominance hypothesis. Under the partial dominance hypothesis, fitness decline occurs due to the increased homozygosity of pre-existing recessive deleterious mutations maintained within a population. In contrast, the overdominance hypothesis posits that fitness reduction arises from a loss of heterozygote advantage. Given that heterozygote advantage is typically observed in a small percentage of loci in an equilibrium population, it is generally assumed that the majority of fitness decline due to inbreeding occurs under the partial dominance hypothesis (Charlesworth et al., 1991; Crow, 1999; Hedrick, 2012; Roff, 2002).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eStudies of inbreeding depression have reported conflicting evidence, generally showing either no epistasis or antagonistic epistasis for fitness (\u0026Aacute;vila et al., 2006; Rand et al., 2006; Sharp \u0026amp; Agrawal, 2016; Sohail et al., 2017; Szafraniec et al., 2003). Some of these discrepancies may be attributable to experimental biases, such as the repeated analysis of the same lines throughout the inbreeding process, sampling bias from the loss of lines due to lethal recessive deleterious mutations, and the effects of purging. However, \u0026nbsp;a recent study by Dom\u0026iacute;nguez-Garc\u0026iacute;a et al., (2019) provided empirical evidence of synergistic epistasis for productivity using a full-sib inbred design in both laboratory and natural populations of \u003cem\u003eDrosophila melanogaster\u003c/em\u003e. While Dom\u0026iacute;nguez-Garc\u0026iacute;a et al. (2019) addressed some of the extant challenges, certain experimental limitations persisted in their experiments, such as only analysing surviving lines which potentially creates a bias against synergistic epistasis among strong deleterious mutations. Another frequently overlooked factor in studies of inbreeding and the pattern of fitness decline is the direct influence of the inbreeding coefficient on fitness components separately. A critical assumption often made is that the accumulation of deleterious mutations, and consequently their homozygous expression, increases at a consistent linear rate across a variety of fitness components affected by these mutations. Under such an assumption, the inbreeding coefficient is expected to affect various life history traits equally leading to a linear decline in the overall fitness.\u003c/p\u003e\n\u003cp\u003eWhile a linear relationship between increasing homozygosity of recessive deleterious mutations and their phenotypic impact is often assumed across all fitness components, it is plausible that different life history traits may exhibit varying rates or non-linear patterns of decline in response to increasing inbreeding coefficients (Charlesworth, 1998; Falconer, 1996; Willis, 1993). Investigating these potential differential responses in fitness components is crucial for a comprehensive understanding of genetic load and epistatic interactions.\u003c/p\u003e\n\u003cp\u003eThis conflicted but developing literature inspired our study, in which we investigated the effect of increasing inbreeding coefficient on several composite traits which determine the fitness of fruit flies. We used a full-sib inbreeding design but replicated parental lines to create ten copies of each initiating genotype. We investigated the patterns of decline in pupal productivity, pupa-to-adult viability, and overall offspring productivity with increasing inbreeding coefficients. Our results revealed an accelerated decline in pupal productivity but a decelerated decline in pupa-to-adult viability. This culminated in an overall linear decline in offspring productivity with increasing inbreeding coefficients. These findings highlight the complexity of studying epistasis among deleterious mutations affecting fitness, especially given that fitness is a highly composite trait. Furthermore, our results underscore the fact that distinct epistatic interactions can influence specific life-history traits without necessarily altering the expected additive effects on overall fitness.\u003c/p\u003e"},{"header":"Materials and Methods","content":"\u003cp\u003eStudy system:\u003c/p\u003e\n\u003cp\u003eOur experimental design utilized a laboratory population of \u003cem\u003eDrosophila melanogaster\u003c/em\u003e designated \u003cem\u003eControl\u003c/em\u003e (C\u003csub\u003e1\u003c/sub\u003e) that is descended from the \u003cem\u003eIves\u0026nbsp;\u003c/em\u003e(IV) base population. IV was established from a wild caught sample of 400 flies (1:1 sex ratio) in Amherst, Massachusetts in 1975 (Rose \u0026amp; Charlesworth, 1981a, 1981b). The C\u003csub\u003e1\u003c/sub\u003e population has a varied demographic history but is outbred and had been maintained under standardized environmental conditions since 1980. These conditions include discrete generations, a population size of approximately 2000 individuals per generation, cultivation on a banana-agar-killed yeast food medium, a moderate larval density (80-100 larvae per vial), a constant temperature of 25\u0026deg;C, and a 12-hour light:12-hour dark photoperiod. This population had been maintained for over 150 generations on a two-week life cycle. With more than 400 total generations of controlled laboratory maintenance, the C\u003csub\u003e1\u003c/sub\u003e population is expected to be at, or close to, an evolutionary equilibrium for standing genetic variation. This characteristic makes it particularly well-suited for estimating mutational load at an equilibrium and measuring fitness under conditions of evolutionary history by assessing various life-history traits.\u003c/p\u003e\n\u003cp\u003eFull-sib Inbreeding:\u003c/p\u003e\n\u003cp\u003eTo establish the inbred lines, 30 unmated male and female pairs were isolated from the C\u003csub\u003e1\u003c/sub\u003e population, with each pair initiating a distinct line within a separate vial. In the second generation, 10 unmated male and female pairs were isolated from the offspring of each established inbred line to create 10 replicate \u0026apos;sub-lines\u0026apos; (300 lines in total). For the subsequent five generations, each sub-line was maintained by randomly collecting a single pair of unmated male and female offspring; a sub-line was terminated if no viable adult pairs were produced in a generation. This \u0026lsquo;sub-line\u0026rsquo; design helps us to account for the possible extinction of lines in the analysis due to strong deleterious effect arising from synergistic epistasis leading to a bias against accelerated decline in surviving lines. In total, these experimental lines underwent six generations of inbreeding, during which fitness decline was assessed.\u003c/p\u003e\n\u003cp\u003eAssessing fitness:\u003c/p\u003e\n\u003cp\u003eThe pattern of fitness decline due to increasing inbreeding was assessed by measuring several key life-history traits.\u0026nbsp;These included pupal productivity as the number of pupae produced by an inbred pair, the pupae to adult viability during pupal stage, and offspring productivity measured as the total number of adults eclosed for each surviving sub-line. Concurrently, these same traits were measured in the outbred control population (C\u003csub\u003e1\u003c/sub\u003e) by sampling 50 male and female pairs each generation, synchronized with the measurements from the inbred lines. For any sub-line that did not survive in a given generation, the zero values were still treated as data in analyses.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAnalysis:\u003c/p\u003e\n\u003cp\u003eAll analyses were performed using R version 4.4.2 (R Core Team, 2024).\u003c/p\u003e\n\u003cp\u003eSince pupal productivity and offspring productivity are composite traits that account for the reproductive fitness of parents, number of eggs laid, hatched and pupated successfully, the expected inbreeding coefficient (\u003cem\u003eICf\u003c/em\u003e) for pupal productivity and offspring productivity was calculated and assigned for each generation of inbreeding by averaging the inbreeding coefficients of the mother and the offspring based on the previous work by\u0026nbsp;Dom\u0026iacute;nguez-Garc\u0026iacute;a et al. (2019).\u0026nbsp;This resulted in expected \u003cem\u003eICf\u003c/em\u003e values of 0, 0.125, 0.313, 0.438, 0.547, 0.633 for the generations 1 to 6, respectively, for the entire genome in a full-sib line. For pupa-to-adult viability and sub-line survival, the expected inbreeding coefficient was calculated based solely on the offspring\u0026apos;s inbreeding coefficient (\u003cem\u003eIC\u003c/em\u003e), yielding values of 0, 0.25, 0.375, 0.5, 0.594, 0.671. Using the lme4 package in R, we fitted linear and quadratic hurdle models to the measured traits against the expected inbreeding coefficient. In our regression models, inbred lines were included as a random effect to account for stochasticity inherent to each homozygous line. Model comparisons between linear and quadratic fits for each trait were performed using Akaike Information Criterion corrected for small sample sizes (AICc), Bayesian Information Criterion (BIC), and log-likelihood comparisons to determine whether the observed decline was linear or non-linear with increasing inbreeding coefficient.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e1. Survival of sublines: The number of surviving sub-lines per inbred line exhibited an accelerating decline of 58% over 6 generations with increasing inbreeding coefficient (\u003cem\u003eICf\u003c/em\u003e), as inferred from comparisons between quadratic and linear model fits (Figure 1, Table 1). This non-linear decline in sub-line survival suggests a significant effect of synergistic epistasis among lethal mutations or synergistic interactions between recessive deleterious mutations affecting sterility or pre-copulatory mortality.\u003c/p\u003e\n\u003cp\u003e2. Pupal productivity: Overall pupal productivity declined by 56.02% in a non-linear fashion with increasing inbreeding coefficient (Figure 2, Table 1). Among the 30 inbred lines, 24 lines showed a significant accelerated decline in pupal productivity, providing strong evidence for synergistic epistasis involving mutations that affect egg-to-pupae viability and development.\u003c/p\u003e\n\u003cp\u003e3. Pupa-to-adult viability: Overall pupa-to-adult viability, representing survival during the pupal stage of pre-adult development, declined by 18.65% in a non-linear fashion with increasing inbreeding coefficient (\u003cem\u003eIC\u003c/em\u003e) (Figure 3, Table 1). Out of the 30 inbred lines, 22 lines exhibited a significant deceleration in pupa-to-adult viability, suggesting strong evidence for antagonistic epistasis among mutations affecting this developmental stage.\u003c/p\u003e\n\u003cp\u003e4. Number of adults eclosed\u003cstrong\u003e:\u003c/strong\u003e The overall number of adults eclosing, representing survival through the entire pre-adult developmental stage, declined by 64.37% in a linear fashion with increasing inbreeding coefficient (Figure 4, Table 1). Among the 30 inbred lines, 27 lines showed a significant linear decline in the number of adults eclosed, which appears counter-intuitive given the results for pupal productivity and pupa-to-adult viability. Instead, this linear pattern can be interpreted as an outcome of the combined effects of accelerated decline in pupal productivity and de-accelerated decline in pupa-to-adult viability, resulting in an overall linear reduction in the number of adults eclosed with increasing \u003cem\u003eICf\u003c/em\u003e.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003eExperiments aimed at deducing the effects of increasing mutation load have yielded a range of outcomes and are equivocal, as a literature, on the question of epistasis for fitness effects. These studies have varied in scale and design, with limits to inference of any one approach. For example, mutation accumulation (MA) experiments have not provided consistent evidence for epistasis among deleterious mutations (\u0026Aacute;vila et al., 2006; Dickinson, 2008; Elena \u0026amp; Lenski, 1997; Kibota \u0026amp; Lynch, 1996; Visser et al., 1997). These investigations typically focus on newly arising spontaneous mutations that would be subject to natural selection and subsequent purging from the population; such spontaneous mutations may not be the ones segregating naturally and having the greatest impact in populations. Consequently, an inbreeding approach offers a potentially valuable alternative approach for examining epistatic interactions among pre-existing recessive deleterious mutations that have accumulated in a population under mutation-selection balance in natural populations (Charlesworth \u0026amp; Hughes, 1996; Dom\u0026iacute;nguez-Garc\u0026iacute;a et al., 2019; Swindell \u0026amp; Bouzat, 2006). Our study population has been maintained for over 400 generations in the lab at large population size, a duration that is expected to result in something approaching an equilibrium state, minimizing the impact of purging and balancing selection on the analysis of epistatic interactions.\u003c/p\u003e\n\u003cp\u003eFew inbreeding studies have been directed at the question of epistasis, however several do suggest the presence of epistatic interactions, both synergistic and antagonistic, affecting various fitness-related traits (Dom\u0026iacute;nguez-Garc\u0026iacute;a et al., 2019; Rosa et al., 2005; Willis, 1993). Given that fitness is a complex, polygenic trait, and researchers often use proxies for total fitness, it is important to dissect epistatic interactions within its constituent traits individually, where feasible. Our current investigation was inspired by Dom\u0026iacute;nguez-Garc\u0026iacute;a et al. (2019) who suggested the existence of synergistic epistasis for fitness, evidenced by an accelerated decline in pupal productivity in \u003cem\u003eD. melanogaster\u0026nbsp;\u003c/em\u003elaboratory and natural populations. Consistent with their findings, our results also indicated an accelerating decline in a similar metric of pupal productivity. However, we also observed a deceleration in the effects of inbreeding on pupa-to-adult viability, which ultimately led to a linear decline in the total number of eclosed adults (total fitness) with increasing inbreeding.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eWhy then is the impact of inbreeding different for pupal productivity and pupal viability? Pupal productivity is a highly composite character, encompassing almost all aspects of total fitness: adult fertility in the parental generation and larval viability in the subsequent generation. \u0026nbsp;Potentially the sheer complexity and number of different ways that productivity could fail through infertility or poor juvenile competition and development, with performance deficits in one phase potentially compounding in the next, creates the conditions for synergism across mutations. In contrast, pupal survival is likely to be an intensely canalized trait, with the metamorphic animal executing a highly choreographed developmental phase highly dependent on the interactions of rare alleles with environmental stress (Bendixsen et al., 2017; Borne et al., 2021; Collet et al., 2023).\u003c/p\u003e\n\u003cp\u003eIn summary, although we are unable to pinpoint the driving elements from the current analysis, we speculate that the expression of deleterious mutations in one life stage could have a potential cascading or \u0026lsquo;synergistic\u0026rsquo; effect in later stages leading to larger than expected fitness decline, without directly interacting with other deleterious mutations. This hypothesis involves condition-dependence, in which erosion of condition below thresholds results in catastrophic failures and morbidity. This mechanism may be complimentary to a study from Sohail et al. (2017), who detected rare loss-of-function alleles exhibiting synergistic epistasis under purifying selection in both human and \u003cem\u003eDrosophila\u003c/em\u003e populations. \u0026nbsp;On the other hand, it explains why purifying selection during mutation-selection balance could promote the accumulation of recessive deleterious mutations exhibiting antagonistic epistasis during the pupal stage in \u003cem\u003eD. melanogaster\u003c/em\u003e. This hypothesis is consistent with evidence of a higher rate of non-synonymous substitutions relative to mutation rate (\u0026omega;) in coding regions during the pupal stage of \u003cem\u003eD. melanogaster\u003c/em\u003e compared to \u0026omega; in the coding regions involved in pupal productivity (Coronado-Zamora et al., 2019).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eOur findings demonstrate that epistatic interactions among deleterious mutations may be strongly trait-specific, to the point of shifting from synergistic to antagonistic depending upon the particular fitness components being assessed. This means that, at least in our study population, even measurements of fertility and survival that capture nearly all aspects of total fitness may misrepresent that parameter, which is of greatest importance from the theoretical standpoint. While synergistic epistasis for total fitness has been implicated in broadly important evolutionary phenomena, such as the evolution of sexual reproduction, our results do not support such a role. However, our study does raise intriguing questions about why different life stages and combinations of fitness characters would yield opposing results regarding epistasis. Hypotheses regarding the relationship between purifying selection, epistasis, and the ratio of non-synonymous to synonymous substitutions (\u0026omega;) can serve as a starting point to investigate such trait-specific epistatic interactions, but their combination in sequence poses special challenges, particularly for highly composite traits. Overall, understanding the genetic and physiological underpinnings of non-additivity will allow a greater insight into the evolutionary architecture of fitness and into the genetic load of populations.\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003e\u003cstrong\u003eConflict of interest:\u0026nbsp;\u003c/strong\u003eThe authors declare they have no conflict of interest.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eData archiving:\u0026nbsp;\u003c/strong\u003eTo be uploaded on\u003cstrong\u003e\u0026nbsp;\u003c/strong\u003ehttps://datadryad.org/\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAcknowledgments:\u003c/strong\u003e\u0026nbsp;\u003cbr\u003e\u0026nbsp;Funding for the work was from NSERC Discovery grants\u0026nbsp;to AKC \u0026amp; TD. We thank all the undergraduate students involved in the population maintenance. We thank Harshavardhan Thygarajan, Mindy Baroody, Avery Want, Josh Kowal for their efforts in the lab.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eAuthor contribution statement:\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIS was responsible for designing the protocol, analysing data, creating figures and tables and interpreting results. IS and HM contributed equally to conducting the research and collecting the data. IS wrote the first draft of the paper. IS and AKC contributed to the updating of reference list and writing the report. IS, TD and AKC contributed to the final draft and revisions.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\n\u003cli\u003e\u0026Aacute;vila, V., Chavarr\u0026iacute;as, D., S\u0026aacute;nchez, E., Manrique, A., L\u0026oacute;pez-Fanjul, C., \u0026amp; Garc\u0026iacute;a-Dorado, A. (2006). Increase of the spontaneous mutation rate in a long-term experiment with \u003cem\u003eDrosophila\u003c/em\u003e \u003cem\u003emelanogaster\u003c/em\u003e. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e173\u003c/em\u003e(1), 267\u0026ndash;277. https://doi.org/10.1534/genetics.106.056200\u003c/li\u003e\n\u003cli\u003eBendixsen, D. P., \u0026Oslash;stman, B., \u0026amp; Hayden, E. J. (2017). Negative epistasis in experimental RNA fitness landscapes. \u003cem\u003eJournal of Molecular Evolution\u003c/em\u003e, \u003cem\u003e85\u003c/em\u003e(5), 159\u0026ndash;168. https://doi.org/10.1007/s00239-017-9817-5\u003c/li\u003e\n\u003cli\u003eBorne, F., Kulathinal, R. J., \u0026amp; Courtier-Orgogozo, V. (2021). Glue genes are subjected to diverse selective forces during \u003cem\u003eDrosophila\u003c/em\u003e development. \u003cem\u003eGenome Biology and Evolution\u003c/em\u003e, \u003cem\u003e13\u003c/em\u003e(12), evab248. \u003cbr\u003e https://doi.org/10.1093/gbe/evab248\u003c/li\u003e\n\u003cli\u003eCharlesworth, B. (1998). The effect of synergistic epistasis on the inbreeding load. \u003cem\u003eGenetics Research\u003c/em\u003e, \u003cem\u003e71\u003c/em\u003e(1), 85\u0026ndash;89. \u003cbr\u003e https://doi.org/10.1017/S0016672398003140\u003c/li\u003e\n\u003cli\u003eCharlesworth, B., \u0026amp; Hughes, K. A. (1996). Age-specific inbreeding depression and components of genetic variance in relation to the evolution of senescence. \u003cem\u003eProceedings of the National Academy of Sciences\u003c/em\u003e, \u003cem\u003e93\u003c/em\u003e(12), 6140\u0026ndash;6145. https://doi.org/10.1073/pnas.93.12.6140\u003c/li\u003e\n\u003cli\u003eCharlesworth, B., Morgan, M. T., \u0026amp; Charlesworth, D. (1991). Multilocus models of inbreeding depression with synergistic selection and partial self-fertilization. \u003cem\u003eGenetics Research\u003c/em\u003e, \u003cem\u003e57\u003c/em\u003e(2), 177\u0026ndash;194. \u003cbr\u003e https://doi.org/10.1017/S0016672300029256\u003c/li\u003e\n\u003cli\u003eCollet, J. M., Nidelet, S., \u0026amp; Fellous, S. (2023). Genetic independence between traits separated by metamorphosis is widespread but varies with biological function. \u003cem\u003eProceedings of the Royal Society B: Biological Sciences\u003c/em\u003e, \u003cem\u003e290\u003c/em\u003e(2010), 20231784. https://doi.org/10.1098/rspb.2023.1784\u003c/li\u003e\n\u003cli\u003eCoronado-Zamora, M., Salvador-Mart\u0026iacute;nez, I., Castellano, D., Barbadilla, A., \u0026amp; Salazar-Ciudad, I. (2019). Adaptation and conservation throughout the \u003cem\u003eDrosophila\u003c/em\u003e \u003cem\u003emelanogaster\u003c/em\u003e life-cycle. \u003cem\u003eGenome Biology and Evolution\u003c/em\u003e, \u003cem\u003e11\u003c/em\u003e(5), 1463\u0026ndash;1482. https://doi.org/10.1093/gbe/evz086\u003c/li\u003e\n\u003cli\u003eCrow, J. F. (1999). Dominance and overdominance. In \u003cem\u003eGenetics and Exploitation of Heterosis in Crops\u003c/em\u003e (pp. 49\u0026ndash;58). https://doi.org/10.2134/1999.geneticsandexploitation.c5\u003c/li\u003e\n\u003cli\u003ede Visser, J. A. G. M., \u0026amp; Elena, S. F. (2007). The evolution of sex: Empirical insights into the roles of epistasis and drift. \u003cem\u003eNature Reviews Genetics\u003c/em\u003e, \u003cem\u003e8\u003c/em\u003e(2), Article 2. https://doi.org/10.1038/nrg1985\u003c/li\u003e\n\u003cli\u003eDickinson, W. J. (2008). Synergistic fitness interactions and a high frequency of beneficial changes among mutations accumulated under relaxed selection in \u003cem\u003eSaccharomyces cerevisiae\u003c/em\u003e. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e178\u003c/em\u003e(3), 1571\u0026ndash;1578. https://doi.org/10.1534/genetics.107.080853\u003c/li\u003e\n\u003cli\u003eDom\u0026iacute;nguez-Garc\u0026iacute;a, S., Garc\u0026iacute;a, C., Quesada, H., \u0026amp; Caballero, A. (2019). Accelerated inbreeding depression suggests synergistic epistasis for deleterious mutations in \u003cem\u003eDrosophila melanogaster\u003c/em\u003e. \u003cem\u003eHeredity\u003c/em\u003e, \u003cem\u003e123\u003c/em\u003e(6), 709\u0026ndash;722.\u003cbr\u003e https://doi.org/10.1038/s41437-019-0263-6\u003c/li\u003e\n\u003cli\u003eElena, S. F., \u0026amp; Lenski, R. E. (1997). Test of synergistic interactions among deleterious mutations in bacteria. \u003cem\u003eNature\u003c/em\u003e, \u003cem\u003e390\u003c/em\u003e(6658), 395\u0026ndash;398.\u003cbr\u003e https://doi.org/10.1038/37108\u003c/li\u003e\n\u003cli\u003eFalconer, D. S. (1996). Introduction to quantitative genetics 4th edition. \u003cem\u003ePearson Education\u003c/em\u003e, ISBN 8131727408, 9788131727409.\u003c/li\u003e\n\u003cli\u003eHedrick, P. W. (2012). What is the evidence for heterozygote advantage selection? \u003cem\u003eTrends in Ecology \u0026amp; Evolution\u003c/em\u003e, \u003cem\u003e27\u003c/em\u003e(12), 698\u0026ndash;704. https://doi.org/10.1016/j.tree.2012.08.012\u003c/li\u003e\n\u003cli\u003eKibota, T. T., \u0026amp; Lynch, M. (1996). Estimate of the genomic mutation rate deleterious to overall fitness in \u003cem\u003eE. coli\u003c/em\u003e. \u003cem\u003eNature\u003c/em\u003e, \u003cem\u003e381\u003c/em\u003e(6584), 694\u0026ndash;696.\u003cbr\u003e https://doi.org/10.1038/381694a0\u003c/li\u003e\n\u003cli\u003eKondrashov, A. S. (1988). Deleterious mutations and the evolution of sexual reproduction. \u003cem\u003eNature\u003c/em\u003e, \u003cem\u003e336\u003c/em\u003e(6198), Article 6198.\u003cbr\u003e https://doi.org/10.1038/336435a0\u003c/li\u003e\n\u003cli\u003eKondrashov, A. S. (1994). Muller\u0026rsquo;s ratchet under epistatic selection. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e136\u003c/em\u003e(4), 1469\u0026ndash;1473.\u003cbr\u003e https://doi.org/10.1093/genetics/136.4.1469\u003c/li\u003e\n\u003cli\u003eR Core Team. (2024). \u003cem\u003eR: The R Project for Statistical Computing\u003c/em\u003e.\u003cbr\u003e https://www.r-project.org/\u003c/li\u003e\n\u003cli\u003eRand, D. M., Fry, A., \u0026amp; Sheldahl, L. (2006). Nuclear\u0026ndash;mitochondrial epistasis and \u003cem\u003eDrosophila\u003c/em\u003e aging: introgression of\u003cem\u003e Drosophila simulans \u003c/em\u003emtdna modifies longevity in \u003cem\u003eD. melanogaster\u003c/em\u003e nuclear backgrounds. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e172\u003c/em\u003e(1), 329\u0026ndash;341. https://doi.org/10.1534/genetics.105.046698\u003c/li\u003e\n\u003cli\u003eRivero, A., Balloux, F., \u0026amp; West, S. A. (2003). Testing for epistasis between deleterious mutations in a parasitoid wasp. \u003cem\u003eEvolution\u003c/em\u003e, \u003cem\u003e57\u003c/em\u003e(7), 1698\u0026ndash;1703. https://doi.org/10.1111/j.0014-3820.2003.tb00375.x\u003c/li\u003e\n\u003cli\u003eRoff, D. A. (2002). Inbreeding depression: Tests of the overdominance and partial dominance hypotheses. \u003cem\u003eEvolution\u003c/em\u003e, \u003cem\u003e56\u003c/em\u003e(4), 768\u0026ndash;775.\u003cbr\u003e https://doi.org/10.1111/j.0014-3820.2002.tb01387.x\u003c/li\u003e\n\u003cli\u003eRosa, J. M., Camacho, S., \u0026amp; Garc\u0026iacute;a‐dorado, A. (2005). A measure of the within‐chromosome synergistic epistasis for \u003cem\u003eDrosophila\u003c/em\u003e viability. \u003cem\u003eJournal of Evolutionary Biology\u003c/em\u003e, \u003cem\u003e18\u003c/em\u003e(4), 1130\u0026ndash;1137.\u003cbr\u003e https://doi.org/10.1111/j.1420-9101.2005.00892.x\u003c/li\u003e\n\u003cli\u003eRose, M. R., \u0026amp; Charlesworth, B. (1981a). Genetics of life history in \u003cem\u003eDrosophila\u003c/em\u003e \u003cem\u003emelanogaster\u003c/em\u003e. I. Sib analysis of adult females. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e97\u003c/em\u003e(1), 173\u0026ndash;186. https://doi.org/10.1093/genetics/97.1.173\u003c/li\u003e\n\u003cli\u003eRose, M. R., \u0026amp; Charlesworth, B. (1981b). Genetics of life history in \u003cem\u003eDrosophila\u003c/em\u003e \u003cem\u003emelanogaster\u003c/em\u003e. II. Exploratory selection experiments. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e97\u003c/em\u003e(1), 187\u0026ndash;196. https://doi.org/10.1093/genetics/97.1.187\u003c/li\u003e\n\u003cli\u003eSharp, N. P., \u0026amp; Agrawal, A. F. (2016). The decline in fitness with inbreeding: Evidence for negative dominance‐by‐dominance epistasis in \u003cem\u003eDrosophila\u003c/em\u003e \u003cem\u003emelanogaster\u003c/em\u003e. \u003cem\u003eJournal of Evolutionary Biology\u003c/em\u003e, \u003cem\u003e29\u003c/em\u003e(4), 857\u0026ndash;864.\u003cbr\u003e https://doi.org/10.1111/jeb.12815\u003c/li\u003e\n\u003cli\u003eSinghal, P., Veturi, Y., Dudek, S. M., Lucas, A., Frase, A., van Steen, K., Schrodi, S. J., Fasel, et al. (2023). Evidence of epistasis in regions of long-range linkage disequilibrium across five complex diseases in the UK Biobank and eMERGE datasets. \u003cem\u003eAmerican Journal of Human Genetics\u003c/em\u003e, \u003cem\u003e110\u003c/em\u003e(4), 575\u0026ndash;591.\u003cbr\u003e https://doi.org/10.1016/j.ajhg.2023.03.007\u003c/li\u003e\n\u003cli\u003eSohail, M., Vakhrusheva, O. A., Sul, J. H., Pulit, S. L., Francioli, L. C., Genome of the Netherlands Consortium, Alzheimer\u0026rsquo;s Disease Neuroimaging Initiative, van den Berg, L. H., Veldink, J. H., de Bakker, P. I. W., Bazykin, G. A., Kondrashov, A. S., \u0026amp; Sunyaev, S. R. (2017). Negative selection in humans and fruit flies involves synergistic epistasis. \u003cem\u003eScience\u003c/em\u003e, \u003cem\u003e356\u003c/em\u003e(6337), 539\u0026ndash;542. \u003cbr\u003e https://doi.org/10.1126/science.aah5238\u003c/li\u003e\n\u003cli\u003eSwindell, W. R., \u0026amp; Bouzat, J. L. (2006). Inbreeding Depression and Male Survivorship in \u003cem\u003eDrosophila\u003c/em\u003e: Implications for Senescence Theory. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e172\u003c/em\u003e(1), 317\u0026ndash;327. https://doi.org/10.1534/genetics.105.045740\u003c/li\u003e\n\u003cli\u003eSzafraniec, K., Wloch, D. M., Sliwa, P., Borts, R. H., \u0026amp; Korona, R. (2003). Small fitness effects and weak genetic interactions between deleterious mutations in heterozygous loci of the yeast \u003cem\u003eSaccharomyces\u003c/em\u003e \u003cem\u003ecerevisiae\u003c/em\u003e. \u003cem\u003eGenetics Research\u003c/em\u003e, \u003cem\u003e82\u003c/em\u003e(1), 19\u0026ndash;31. https://doi.org/10.1017/S001667230300630X\u003c/li\u003e\n\u003cli\u003eVisser, J. A. G. M. de, Hoekstra, R. F., \u0026amp; Ende, H. van den. (1997). Test of interaction between genetic markers that affect fitness in \u003cem\u003eAspergillus niger\u003c/em\u003e. \u003cem\u003eEvolution\u003c/em\u003e, \u003cem\u003e51\u003c/em\u003e(5), 1499\u0026ndash;1505.\u003cbr\u003e https://doi.org/10.1111/j.1558-5646.1997.tb01473.x\u003c/li\u003e\n\u003cli\u003eWachter, K. W., Evans, S. N., \u0026amp; Steinsaltz, D. (2013). The age-specific force of natural selection and biodemographic walls of death. \u003cem\u003eProceedings of the National Academy of Sciences of the United States of America\u003c/em\u003e, \u003cem\u003e110\u003c/em\u003e(25), 10141\u0026ndash;10146. https://doi.org/10.1073/pnas.1306656110\u003c/li\u003e\n\u003cli\u003eWagner, A., Wagner, G. P., \u0026amp; Similion, P. (1994). Epistasis can facilitate the evolution of reproductive isolation by peak shifts: A two-locus two-allele model. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e138\u003c/em\u003e(2), 533\u0026ndash;545. \u003cbr\u003e https://doi.org/10.1093/genetics/138.2.533\u003c/li\u003e\n\u003cli\u003eWeinreich, D. M., Watson, R. A., \u0026amp; Chao, L. (2005). Perspective: Sign epistasis and genetic constraint on evolutionary trajectories. \u003cem\u003eEvolution; International Journal of Organic Evolution\u003c/em\u003e, \u003cem\u003e59\u003c/em\u003e(6), 1165\u0026ndash;1174. PMID: 16050094.\u003c/li\u003e\n\u003cli\u003eWest, S. A., Peters, A. D., \u0026amp; Barton, N. H. (1998). Testing for epistasis between deleterious mutations. \u003cem\u003eGenetics\u003c/em\u003e, \u003cem\u003e149\u003c/em\u003e(1), 435\u0026ndash;444. https://doi.org/10.1093/genetics/149.1.435\u003c/li\u003e\n\u003cli\u003eWhitlock, M. C., \u0026amp; Bourguet, D. (2000). Factors affecting the genetic load in \u003cem\u003eDrosophila\u003c/em\u003e: synergistic epistasis and correlations among fitness components. \u003cem\u003eEvolution\u003c/em\u003e, \u003cem\u003e54\u003c/em\u003e(5), 1654\u0026ndash;1660.\u003cbr\u003e https://doi.org/10.1111/j.0014-3820.2000.tb00709.x\u003c/li\u003e\n\u003cli\u003eWillis, J. H. (1993). Effects of different levels of inbreeding on fitness components in \u003cem\u003emimulus\u003c/em\u003e \u003cem\u003eguttatus\u003c/em\u003e. \u003cem\u003eEvolution\u003c/em\u003e, \u003cem\u003e47\u003c/em\u003e(3), 864\u0026ndash;876.\u003cbr\u003e https://doi.org/10.1111/j.1558-5646.1993.tb01240.x\u003c/li\u003e\n\u003cli\u003eZhang, J., Chen, F., \u0026amp; Li, X. (2024). Mechanistic causes of sign epistasis and its applications. \u003cem\u003eFrontiers in Genetics\u003c/em\u003e, \u003cem\u003e15\u003c/em\u003e, 1366917. https://doi.org/10.3389/fgene.2024.1366917\u003cstrong\u003e\u003cbr\u003e \u003c/strong\u003e\u003c/li\u003e\n\u003c/ol\u003e"},{"header":"Table","content":"\u003cp\u003eTable 1 is available in the Supplementary Files section.\u003c/p\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"","lastPublishedDoi":"10.21203/rs.3.rs-7910627/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-7910627/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"Epistasis has been theoretically implicated in several major evolutionary processes, including the evolution of sex and speciation, but empirical evidence for its impact is still lacking. We tested for epistatic interactions using a full-sib inbreeding design in \u003ci\u003eDrosophila melanogaster\u003c/i\u003e. We examined the effects of deleterious mutations on two sequential components of life-history: pupal productivity and pupa-to-adult viability. We observed an accelerating decline in pupal productivity with increasing inbreeding coefficient, consistent with synergistic epistasis, whereas pupa-to-adult viability exhibited a decelerating decline, consistent with antagonistic epistasis. This decoupling led to an overall linear decline in the production of adult offspring with increasing inbreeding. Our findings suggest that epistatic interactions can be highly trait-specific, underscoring the importance of analyzing individual, physiologically interconnected, fitness traits to understand the role of epistasis in evolutionary processes.","manuscriptTitle":"Additivity of the genetic load revealed to be a product of both synergistic and antagonistic epistasis in Drosophila.","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2025-11-14 19:02:29","doi":"10.21203/rs.3.rs-7910627/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"4a443445-7aaa-4143-8164-abdc4fa42478","owner":[],"postedDate":"November 14th, 2025","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":57438787,"name":"Biological sciences/Evolution/Evolutionary theory"},{"id":57438788,"name":"Biological sciences/Evolution/Experimental evolution"}],"tags":[],"updatedAt":"2025-12-08T16:57:19+00:00","versionOfRecord":[],"versionCreatedAt":"2025-11-14 19:02:29","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-7910627","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-7910627","identity":"rs-7910627","version":["v1"]},"buildId":"8U1c8b4HqxoKbykW_rLl7","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}