Dating nodes in a phylogeny is an important problem in evolution and is typically performed by using molecular clocks and fossil age estimates [1]. The manuscript by Chauve et al. [2] reports a novel method, which uses lateral gene transfers to help ordering nodes in a species tree. The idea is that a lateral gene transfer can only occur between two species living at the same time, which indirectly informs on node relative ages in a phylogeny: the donor species cannot be more recent than the recipient species. Horizontal gene transfers are increasingly recognized as frequent, even in eukaryotes, and especially in micro-organisms that have little fossil records [3-7]. Yet, such an important source of information has been very rarely used so far for inferring relative node ages in phylogenies. In this context, the method by Chauve et al. [2] represents an innovative and original approach to a difficult problem. An obvious limitation of the approach is that it relies on inferences of horizontal transfers, which detection is in itself a difficult problem. Incomplete taxon sampling, or the extinction of the true donor lineage may render patterns difficult to interpret in a temporary fashion. Yet, for clades with no fossils this may be the only piece of information we have at hand, and the growing amount of sequence data is likely to minimize issues derived from incomplete sampling.

The developed method, MaxTiC (for Maximal Time Consistency) [2], represents a very nice application of theoretical developments on the well-known « Feedback Arc Set » computer science problem to the evolutionary question of ordering nodes in a phylogeny. MaxTiC uses as input a species tree and a set of time constraints based on lateral gene transfers inferred using other softwares, and minimizes conflicts between node ordering and these time constraints. The application of MaxTiC on simulated datasets indicated that node ordering was fairly accurate [2]. MaxTiC is implemented in a freely available software, which represents original and relevant contribution to the field of evolutionary biology.

References

[1] Donoghue P and Smith M, editors. 2003. Telling the evolutionary time. CRC press.

[2] Chauve C, Rafiey A, Davin AA, Scornavacca C, Veber P, Boussau B, Szöllősi GJ, Daubin V and Tannier E. 2017. MaxTiC: Fast ranking of a phylogenetic tree by Maximum Time Consistency with lateral gene transfers. bioRxiv 127548, ver. 6 of 6th November 2017. doi: 10.1101/127548

[3] Ropars J, Rodríguez de la Vega RC, Lopez-Villavicencio M, Gouzy J, Sallet E, Debuchy R, Dupont J, Branca A and Giraud T. 2015. Adaptive horizontal gene transfers between multiple cheese-associated fungi. Current Biology 19, 2562–2569. doi: 10.1016/j.cub.2015.08.025

[4] Novo M, Bigey F, Beyne E, Galeote V, Gavory F, Mallet S, Cambon B, Legras JL, Wincker P, Casaregola S and Dequin S. 2009. Eukaryote-to-eukaryote gene transfer events revealed by the genome sequence of the wine yeast Saccharomyces cerevisiae EC1118. Proceeding of the National Academy of Science USA, 106, 16333–16338. doi: 10.1073/pnas.0904673106

[5] Naranjo-Ortíz MA, Brock M, Brunke S, Hube B, Marcet-Houben M, Gabaldón T. 2016. Widespread inter- and intra-domain horizontal gene transfer of d-amino acid metabolism enzymes in Eukaryotes. Frontiers in Microbiology 7, 2001. doi: 10.3389/fmicb.2016.02001

[6] Alexander WG, Wisecaver JH, Rokas A, Hittinger CT. 2016. Horizontally acquired genes in early-diverging pathogenic fungi enable the use of host nucleosides and nucleotides. Proceeding of the National Academy of Science USA. 113, 4116–4121. doi: 10.1073/pnas.1517242113

[7] Marcet-Houben M, Gabaldón T. 2010. Acquisition of prokaryotic genes by fungal genomes. Trends in Genetics. 26, 5–8. doi: 10.1016/j.tig.2009.11.007

Domesticated organisms are excellent models for understanding ecology and evolution and they are important for our food production and safety. While less studied than plants and animals, micro-organisms have also been domesticated, in particular for food fermentation [1]. The most studied domesticated micro-organism is the yeast used to make wine, beer and bread, Saccharomyces cerevisiae [2, 3, 4].

Filamentous fungi used for cheese-making have recently gained interest, for example Penicillium roqueforti used to make blue cheeses and P. camemberti to make soft cheeses [5, 6, 7, 8]. As for plants and animals, domestication has led to beneficial traits for food production in fermenting fungi, but also to bottlenecks and degeneration [6, 7, 9]; P. camemberti for example does not produce enough spores any more for optimal culture and inoculation and P. roqueforti has lost sexual fertility [9]. The loss of genetic diversity and of species diversity in our food production system is concerning for multiple reasons : i) it jeopardizes future improvement in the face of global changes ; ii) it causes the loss of evolved diversity during centuries under human selection, and therefore of beneficial characteristics and specificities that we may never be able to recover ; iii) it leads to degeneration in the few cultivated strains; iv) it impoverishes the diversity of our food products and local adaptation of production practices. 

The study of domesticated fungi used for food fermentation has focused so far on the evolution of lineages and on their metabolic specificities. Microbiological assemblages and species diversity have been much less studied, while they likely also have a strong impact on the quality and safety of final products. This study by Elisa Michel and colleagues [10] addresses this question, using an interdisciplinary participatory research approach including bakers, psycho-sociologists and microbiologists to analyse bread-making practices and their impact on microbial communities in sourdough. Elisa Michel and colleagues [10] identified two distinct groups of bread-making practices based on interviews and surveys, with farmer-like practices (low bread production, use of ancient wheat populations, manual kneading, working at ambient temperature, long fermentation periods and no use of commercial baker’s yeast) versus more intensive, artisanal-like practices. Metabarcoding and microbial culture-based analyses showed that the well-known baker’s yeast, Saccharomyces cerevisiae, was, surprisingly, not the most common species in French organic sourdoughs. Kazachstania was the most represented yeast genus over all sourdoughs, both in terms of read abundance and of species diversity. Kazachstania species were also often dominant in individual sourdoughs, but Saccharomyces uvarum or Torulaspora delbrueckii could also be the dominant yeast species.

Metabarcoding analyses further revealed that the composition of the fungal communities differed between the farmer-like and more intensive practices, representing the first evidence of the influence of artisanal practices on microbial communities. The fungal communities were impacted by a combination of bread-making variables including the type of wheat varieties, the length of fermentation, the quantity of bread made per week and the use of commercial yeast. Maintaining on farm less intensive bread-making practices, may allow the preservation of typical species and phenotypic diversity in microbial communities in sourdough. Farmer-like practices did not lead to higher diversity within sourdoughs but, overall, the diversity of bread-making practices allow maintaining a larger diversity in sourdoughs. For example, different Kazachstania species were most abundant in sourdoughs from artisanal-like and farmer-like practices. Interviews with the bakers suggested the role of dispersal of Kazachstania species in shaping sourdough microbial communities, dispersal occurring by seed exchanges, sourdough mixing or gifts, bread-making training in common or working in one another’s bakery. Nikolai Vavilov [11] had already highlighted for crops the importance of isolated cultures and selection in different farms for generating and preserving crop diversity, but also the importance of seed exchange for fostering adaptation. 

Furthermore, one of the yeast frequently found in artisanal sourdoughs, Kazachstania humilis, displayed phenotypic differences between sourdough and non-sourdough strains, suggesting domestication. The sourdough strains exhibited significantly higher CO2 production rate and a lower fermentation latency-phase time. 

The study by Elisa Michel and colleagues [10] is thus novel and inspiring in showing the importance of interdisciplinary studies, combining metabarcoding, microbiology and interviews for assessing the composition and diversity of microbial communities in human-made food, and in revealing the impact of artisanal-like bread-making practices in preserving microbial community diversity.

Interdisciplinary studies are still rare but have already shown the importance of combining ethno-ecology, biology and evolution to decipher the role of human practices on genetic diversity in crops, animals and food microorganisms and to help preserving genetic resources [12]. For example, in the case of the bread wheat Triticum aestivum, such interdisciplinary studies have shown that genetic diversity has been shaped by farmers’ seed diffusion and farming practices [13]. We need more of such interdisciplinary studies on the impact of farmer versus industrial agricultural and food-making practices as we urgently need to preserve the diversity of micro-organisms used in food production that we are losing at a rapid pace [6, 7, 14]. 

References

[1] Dupont J, Dequin S, Giraud T, Le Tacon F, Marsit S, Ropars J, Richard F, Selosse M-A (2017) Fungi as a Source of Food. Microbiology Spectrum, 5, 5.3.09. https://doi.org/10.1128/microbiolspec.FUNK-0030-2016

[2] Legras J-L, Galeote V, Bigey F, Camarasa C, Marsit S, Nidelet T, Sanchez I, Couloux A, Guy J, Franco-Duarte R, Marcet-Houben M, Gabaldon T, Schuller D, Sampaio JP, Dequin S (2018) Adaptation of S. cerevisiae to Fermented Food Environments Reveals Remarkable Genome Plasticity and the Footprints of Domestication. Molecular Biology and Evolution, 35, 1712–1727. https://doi.org/10.1093/molbev/msy066

[3] Bai F-Y, Han D-Y, Duan S-F, Wang Q-M (2022) The Ecology and Evolution of the Baker’s Yeast Saccharomyces cerevisiae. Genes, 13, 230. https://doi.org/10.3390/genes13020230

[4] Fay JC, Benavides JA (2005) Evidence for Domesticated and Wild Populations of Saccharomyces cerevisiae. PLOS Genetics, 1, e5. https://doi.org/10.1371/journal.pgen.0010005

[5] Ropars J, Rodríguez de la Vega RC, López-Villavicencio M, Gouzy J, Sallet E, Dumas É, Lacoste S, Debuchy R, Dupont J, Branca A, Giraud T (2015) Adaptive Horizontal Gene Transfers between Multiple Cheese-Associated Fungi. Current Biology, 25, 2562–2569. https://doi.org/10.1016/j.cub.2015.08.025

[6] Dumas E, Feurtey A, Rodríguez de la Vega RC, Le Prieur S, Snirc A, Coton M, Thierry A, Coton E, Le Piver M, Roueyre D, Ropars J, Branca A, Giraud T (2020) Independent domestication events in the blue-cheese fungus Penicillium roqueforti. Molecular Ecology, 29, 2639–2660. https://doi.org/10.1111/mec.15359

[7] Ropars J, Didiot E, Rodríguez de la Vega RC, Bennetot B, Coton M, Poirier E, Coton E, Snirc A, Le Prieur S, Giraud T (2020) Domestication of the Emblematic White Cheese-Making Fungus Penicillium camemberti and Its Diversification into Two Varieties. Current Biology, 30, 4441-4453.e4. https://doi.org/10.1016/j.cub.2020.08.082

[8] Caron T, Piver ML, Péron A-C, Lieben P, Lavigne R, Brunel S, Roueyre D, Place M, Bonnarme P, Giraud T, Branca A, Landaud S, Chassard C (2021) Strong effect of Penicillium roqueforti populations on volatile and metabolic compounds responsible for aromas, flavor and texture in blue cheeses. International Journal of Food Microbiology, 354, 109174. https://doi.org/10.1016/j.ijfoodmicro.2021.109174

[9] Ropars J, Lo Y-C, Dumas E, Snirc A, Begerow D, Rollnik T, Lacoste S, Dupont J, Giraud T, López-Villavicencio M (2016) Fertility depression among cheese-making Penicillium roqueforti strains suggests degeneration during domestication. Evolution, 70, 2099–2109. https://doi.org/10.1111/evo.13015

[10] Michel E, Masson E, Bubbendorf S, Lapicque L, Nidelet T, Segond D, Guézenec S, Marlin T, Devillers H, Rué O, Onno B, Legrand J, Sicard D, Bakers TP (2022) Artisanal and farmer bread making practices differently shape fungal species community composition in French sourdoughs. bioRxiv, 679472, ver. 6 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/679472

[11] Vavilov NI, Vavylov MI, Dorofeev VF (1992) Origin and Geography of Cultivated Plants. Cambridge University Press.

[12] Saslis-Lagoudakis CH, Clarke AC (2013) Ethnobiology: the missing link in ecology and evolution. Trends in Ecology & Evolution, 28, 67–68. https://doi.org/10.1016/j.tree.2012.10.017

[13] Thomas M, Demeulenaere E, Dawson JC, Khan AR, Galic N, Jouanne-Pin S, Remoue C, Bonneuil C, Goldringer I (2012) On-farm dynamic management of genetic diversity: the impact of seed diffusions and seed saving practices on a population-variety of bread wheat. Evolutionary Applications, 5, 779–795. https://doi.org/10.1111/j.1752-4571.2012.00257.x

[14] Demeulenaere É, Lagrola M (2021) Des indicateurs pour accompagner “ les éleveurs de microbes” : Une communauté épistémique face au problème des laits “ paucimicrobiens ” dans la production fromagère au lait cru (1995-2015). Revue d’anthropologie des connaissances, 15. http://journals.openedition.org/rac/24953

Artificial selection experiments are key experiments in evolutionary biology. The demonstration that application of selective pressure across multiple generations results in heritable phenotypic changes is a tangible and reproducible proof of the evolution by natural selection.
Artificial selection experiments are used to evaluate the joint effects of selection on multiple traits, their genetic covariances and differences in responses in different environments. Most studies on artificial selection experiments report and base their analyses on phenotypic changes [1]. More recently, changes in allele frequency and other patterns of molecular genetic diversity have been used to identify genomic locations where selection has had an effect. However, so far the changes in gene expression have not been in the focus of artificial selection experiment studies (see [2] for an example though).
In plants, one of the most famous artificial selection experiments is the Illinois Corn Experiment where maize (Zea mays) is selected for oil and protein content [3], but in addition, similar experiments have been conducted also for other traits in maize. In Saclay divergent selection experiment [4] two maize inbred lines (F252 and MBS847) have been selected for early and late flowering for 13 generations, resulting in two week difference in flowering time.
In ”Transcriptomic response to divergent selection for flowering time in maize reveals convergence and key players of the underlying gene regulatory network ” [5] Maud Tenaillon and her coworkers study the gene expression differences among these two independently selected maize populations. Their experiments cover two years in field conditions and they use samples of shoot apical meristem at three different developmental stages: vegetative, transitioning and reproductive. They use RNA-seq transcriptome level differences and qRT-PCR for gene expression pattern investigation. The work is continuation to earlier genetic and phenotypic studies on the same material [4, 6].
The reviewers and I agree that dataset is unique and its major benefit is that it has been obtained from field conditions similar to those that species may face under natural setting during selection. Their tissue sampling is supported by flowering time phenotypic observations and covers the developmental transition stage, making a good effort to identify key transcriptional and phenotypic changes and their timing affected by selection.
Tenaillon et al. [5] identify more than 2000 genes that are differentially expressed among early and late flowering populations. Expectedly, they are enriched for known flowering time genes. As they point out, differential expression of thousands of genes does not mean that they all were independently affected by selection, but rather that the whole transcriptional network has shifted, possibly due to just few upstream or hub-genes. Also, the year-to-year variation had smaller effect in gene expression compared to developmental stage or genetic background, possibly indicating selection for stability across environmental fluctuation for such an important phenotype as flowering time.
Another noteworthy observation is that they find convergent patterns of transcriptional changes among the two selected lines. 115 genes expression patterns are shifted due to selection in both genetic backgrounds. This convergent pattern can be a result of either selection on standing variation or de novo mutations. The data does not allow testing which process is underlying the observed convergence. However, their results show that this is an interesting future question that can be addressed using genotype and gene expression data from the same ancestral and derived material and possibly their hybrids.

References

[1] Hill, W. G., & Caballero, A. (1992). Artificial selection experiments. Annual Review of Ecology and Systematics, 23(1), 287-310. doi: 10.1146/annurev.es.23.110192.001443
[2] Konczal, M., Babik, W., Radwan, J., Sadowska, E. T., & Koteja, P. (2015). Initial molecular-level response to artificial selection for increased aerobic metabolism occurs primarily through changes in gene expression. Molecular biology and evolution, 32(6), 1461-1473. doi: 10.1093/molbev/msv038
[3] Moose, S. P., Dudley, J. W., & Rocheford, T. R. (2004). Maize selection passes the century mark: a unique resource for 21st century genomics. Trends in plant science, 9(7), 358-364. doi: 10.1016/j.tplants.2004.05.005
[4] Durand, E., Tenaillon, M. I., Ridel, C., Coubriche, D., Jamin, P., Jouanne, S., Ressayre, A., Charcosset, A. and Dillmann, C. (2010). Standing variation and new mutations both contribute to a fast response to selection for flowering time in maize inbreds. BMC evolutionary biology, 10(1), 2. doi: 10.1186/1471-2148-10-2
[5] Tenaillon, M. I., Seddiki, K., Mollion, M., Le Guilloux, M., Marchadier, E., Ressayre, A. and Dillmann C. (2019). Transcriptomic response to divergent selection for flowering time in maize reveals convergence and key players of the underlying gene regulatory network. BioRxiv, 461947 ver. 5 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/461947
[6] Durand, E., Tenaillon, M. I., Raffoux, X., Thépot, S., Falque, M., Jamin, P., Bourgais A., Ressayre, A. and Dillmann, C. (2015). Dearth of polymorphism associated with a sustained response to selection for flowering time in maize. BMC evolutionary biology, 15(1), 103. doi: 10.1186/s12862-015-0382-5

Most mutations affecting fitness are deleterious and they have many evolutionary consequences. The dynamics and consequences of deleterious mutations are a long-standing question in evolutionary biology and a strong theoretical background has already been developed, for example, to predict the mutation load, inbreeding depression or background selection. One of the classical results is that inbreeding helps purge partially recessive deleterious mutations by exposing them to selection in homozygotes. However, this mainly results from single-locus considerations. When interactions among several, more or less linked, deleterious mutations are taken into account, peculiar dynamics can emerge. One of them, called pseudo-overdominance (POD), corresponds to the maintenance in a population of two (or more) haplotype blocks composed of several recessive deleterious mutations in repulsion that mimics overdominance. Indeed, homozygote individuals for one of the haplotype blocks expose many deleterious mutations to selection whereas they are reciprocally masked in heterozygotes, leading to higher fitness of heterozygotes compared to both homozygotes. A related process, called associative overdominance (AOD) is the effect of such deleterious alleles in repulsion on the linked neutral variation that can be increased by AOD. Although this possibility has been recognized for a long time (Otha and Kimura 1969), it has been mainly considered an anecdotal process. Recently, both theoretical (Zhao and Charlesworth 2016) and genomic analyses (Gilbert et al. 2020) have renewed interest in such a process, suggesting that it could be important in weakly recombining regions of a genome. Donald Waller (2021) - one of the co-authors of the current work - also recently proposed that POD could be quantitatively important with broad implications, and could resolve some unexplained observations such as the maintenance of inbreeding depression in highly selfing species. Yet, a proper theoretical framework analysing the effect of inbreeding on POD was lacking.

In this theoretical work, Diala Abu Awad and Donald Waller (2022) addressed this question through an elegant combination of analytical predictions and intensive multilocus simulations. They determined the conditions under which POD can be maintained and how long it could resist erosion by recombination, which removes the negative association between deleterious alleles (repulsion) at the core of the mechanism. They showed that under tight linkage, POD regions can persist for a long time and generate substantial segregating load and inbreeding depression, even under inbreeding, so opposing (for a while) to the purging effect. They also showed that background selection can affect the genomic structure of POD regions by rapidly erasing weak POD regions but maintaining strong POD regions (i.e with many tightly linked deleterious alleles).

These results have several implications. They can explain the maintenance of inbreeding depression despite inbreeding (as anticipated by Waller 2021), which has implications for the evolution of mating systems. If POD can hardly emerge under high selfing, it can persist from an outcrossing ancestor long after the transition towards a higher selfing rate and could explain the maintenance of mixed mating systems(which is possible with true overdominance, see Uyenoyama and Waller 1991). The results also have implications for genomic analyses, pointing to regions of low or no recombination where POD could be maintained, generating both higher diversity and heterozygosity than expected and variance in fitness. As structural variations are likely widespread in genomes with possible effects on suppressing recombination (Mérot et al. 2020), POD regions should be checked more carefully in genomic analyses (see also Gilbert et al. 2020).

Overall, this work should stimulate new theoretical and empirical studies, especially to assess how quantitatively strong and widespread POD can be. It also stresses the importance of properly considering genetic linkage genome-wide, and so the role of recombination landscapes in determining patterns of diversity and fitness effects.

Awad DA, Waller D (2022) Conditions for maintaining and eroding pseudo-overdominance and its contribution to inbreeding depression. bioRxiv, 2021.12.16.473022, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.12.16.473022

Gilbert KJ, Pouyet F, Excoffier L, Peischl S (2020) Transition from Background Selection to Associative Overdominance Promotes Diversity in Regions of Low Recombination. Current Biology, 30, 101-107.e3. https://doi.org/10.1016/j.cub.2019.11.063

Mérot C, Oomen RA, Tigano A, Wellenreuther M (2020) A Roadmap for Understanding the Evolutionary Significance of Structural Genomic Variation. Trends in Ecology & Evolution, 35, 561–572. https://doi.org/10.1016/j.tree.2020.03.002

Ohta T, Kimura M (1969) Linkage disequilibrium at steady state determined by random genetic drift and recurrent mutation. Genetics, 63, 229–238. https://doi.org/10.1093/genetics/63.1.229

Uyenoyama MK, Waller DM (1991) Coevolution of self-fertilization and inbreeding depression II. Symmetric overdominance in viability. Theoretical Population Biology, 40, 47–77. https://doi.org/10.1016/0040-5809(91)90046-I

Waller DM (2021) Addressing Darwin’s dilemma: Can pseudo-overdominance explain persistent inbreeding depression and load? Evolution, 75, 779–793. https://doi.org/10.1111/evo.14189

Zhao L, Charlesworth B (2016) Resolving the Conflict Between Associative Overdominance and Background Selection. Genetics, 203, 1315–1334. https://doi.org/10.1534/genetics.116.188912

The evolution of selfing results from a balance between multiple evolutionary forces. Selfing provides an "automatic advantage" due to the higher efficiency of selfers to transmit their genes via selfed and outcrossed offspring. Selfed offspring, however, may suffer from inbreeding depression. In principle the ultimate evolutionary outcome is easy to predict from the relative magnitude of these two evolutionary forces [1,2]. Yet, several studies explicitly taking into account the genetic architecture of inbreeding depression noted that these predictions are often too restrictive because selfing can evolve in a broader range of conditions [3,4].
The present work by Abu Awad and Roze [5] provides an analytic understanding of these results. Abu Awad and Roze analyse the evolution of selfing in a multilocus model where some loci are coding for selfing while others are under direct selection. The evolution of selfing depends on (i) the classical benefit of selfing (automatic advantage), (ii) the cost of selfing due to inbreeding depression, (iii) the association between the loci coding for selfing and the loci under direct selection (likely to be positive because selfing is expected to be found in better purged genetic backgrounds) and (iv) the association between the loci coding for selfing and the linkage between loci under selection (this final term depends on the magnitude and the type of epistasis). Because these last two terms depend on genetic associations they are expected to play in when selection is strong and recombination is small. These last two terms explain why selfing is evolving under a range of conditions which is broader than predicted by earlier theoretical models. The match between the approximations for the different terms acting on the evolution of selfing and individual based simulations (for different fitness landscapes) is very convincing. In particular, this analysis also yields new results on the effect of different types of epistasis on inbreeding depression.
Another remarkable and important feature of this work is its readability. The analysis of multilocus models rely on several steps and approximations that often result in overwhelmingly complex papers. Abu Awad and Roze’s paper [5] is dense but it provides a very clear and comprehensive presentation of the interplay between multiple evolutionary forces acting on the evolution of selfing.

References

[1] Holsinger, K. E., Feldman, M. W., and Christiansen, F. B. (1984). The evolution of self-fertilization in plants: a population genetic model. The American Naturalist, 124(3), 446-453. doi: 10.1086/284287
[2] Lande, R., and Schemske, D. W. (1985). The evolution of self‐fertilization and inbreeding depression in plants. I. Genetic models. Evolution, 39(1), 24-40. doi: 10.1111/j.1558-5646.1985.tb04077.x
[3] Charlesworth, D., Morgan, M. T., and Charlesworth, B. (1990). Inbreeding depression, genetic load, and the evolution of outcrossing rates in a multilocus system with no linkage. Evolution, 44(6), 1469-1489. doi: 10.1111/j.1558-5646.1990.tb03839.x
[4] Uyenoyama, M. K., and Waller, D. M. (1991). Coevolution of self-fertilization and inbreeding depression I. Mutation-selection balance at one and two loci. Theoretical population biology, 40(1), 14-46. doi: 10.1016/0040-5809(91)90045-H
[5] Abu Awad, D. and Roze, D. (2020). Epistasis, inbreeding depression and the evolution of self-fertilization. bioRxiv, 809814, ver. 4 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/809814