Mathematical models play an essential role in current evolutionary biology, and evolutionary epidemiology is not an exception [1]. While the issues of virulence evolution and drug-resistance evolution resonate in the literature for quite some time [2, 3], the study by Alizon [4] is one of a few that consider co-evolution of both these traits [5]. The idea behind this study is the following: treating individuals with more severe symptoms at a higher rate (which appears to be quite natural) leads to an appearance of virulent drug-resistant strains, via treatment failure. The author then shows that virulence in drug-resistant strains may face different selective pressures than in drug-sensitive strains and hence proceed at different rates. Hence, treatment itself modulates evolution of virulence. As one of the reviewers emphasizes, the present manuscript offers a mathematical view on why the resistant and more virulent strains can be selected in epidemics. Also, we both find important that the author highlights that the topic and results of this study can be attributed to public health policies and development of optimal treatment protocols [6].
Mathematical models are simplified representations of reality, created with a particular purpose. It can be simple as well as complex, but even simple models can produce relatively complex and knitted results. The art of modelling thus lies not only in developing a model, but also in interpreting and unknitting the results. And this is what Alizon [4] indeed does carefully and exhaustively. Using two contrasting theoretical approaches to study co-evolution, the Price equation approach to study short-term evolution and the adaptive dynamics approach to study long-term evolution, Alizon [4] shows that a positive correlation between the rate of treatment and infection severity causes virulence in drug-sensitive strains to decrease. Clearly, no single model can describe and explain an examined system in its entirety, and even this aspect of the work is taken seriously. Many possible extensions of the study are laid out, providing a wide opportunity to pursue this topic even further. Personally, I have had an opportunity to read many Alizon’s papers and use, teach or discuss many of his models and results. All, including the current one, keep high standard and pursue the field of theoretical (evolutionary) epidemiology.

References

[1] Gandon S, Day T, Metcalf JE, Grenfell BT (2016) Forecasting epidemiological and evolutionary dynamics of infectious diseases. Trends Ecol Evol 31: 776-788. doi: https://doi.org/10.1016/j.tree.2016.07.010
[2] Berngruber TW, Froissart R, Choisy M, Gandon S (2013) Evolution of virulence in emerging epidemics. PLoS Pathog 9(3): e1003209. doi: https://doi.org/10.1371/journal.ppat.1003209
[3] Spicknall IH, Foxman B, Marrs CF, Eisenberg JNS (2013) A modeling framework for the evolution and spread of antibiotic resistance: literature review and model categorization. Am J Epidemiol 178: 508-520. doi: https://doi.org/10.1093/aje/kwt017
[4] Alizon S (2020) Treating symptomatic infections and the co-evolution of virulence and drug resistance. bioRxiv, 2020.02.29.970905, ver. 3 peer-reviewed and recommended by PCI Evol Biol. doi: https://doi.org/10.1101/2020.02.29.970905
[5] Carval D, Ferriere R (2010) A unified model for the coevolution of resistance, tolerance, and virulence. Evolution 64: 2988–3009. doi: https://doi.org/10.1111/j.1558-5646.2010.01035.x
[6] Read AF, T Day, and S Huijben (2011). The evolution of drug resistance and the curious orthodoxy of aggressive chemotherapy. Proc Natl Acad Sci USA 108 Suppl 2, 10871–7. doi: https://doi.org/10.1073/pnas.1100299108

The genome of eukaryotic species is a complex structure that experience many different interactions within itself and with the surrounding environment. The genetic architecture of a phenotype (that is, the set of genetic elements affecting a trait of the organism) plays a fundamental role in understanding the adaptation process of a species to, for example, different climate environments, or to its interaction with other species. Thus, it is fundamental to study the different aspects of the genetic architecture of the species and its relationship with its surronding environment. Aspects such as modularity (the number of genetic units and the degree to which each unit is affecting a trait of the organism), pleiotropy (the number of different effects that a genetic unit can have on an organism) or linkage (the degree of association between the different genetic units) are essential to understand the genetic architecture and to interpret the effects of selection on the genome. Indeed, the knowledge of the different aspects of the genetic architecture could clarify whether genes are affected by multiple aspects of the environment or, on the contrary, are affected by only specific aspects [1,2].

The work performed by Lotterhos et al. [3] sought to understand the genetic architecture of the adaptation to different environments in lodgepole pine (Pinus contorta), considering as candidate SNPs those previously detected as a result of its extreme association patterns to different environmental variables or to extreme population differentiation. This consideration is very important because the study is only relevant if the studied markers are under the effect of selection. Otherwise, the genetic architecture of the adaptation to different environments would be masked by other (neutral) kind of associations that would be difficult to interpret [4,5]. In order to understand the relationship between genetic architecture and adaptation, it is relevant to detect the association networks of the candidate SNPs with climate variables (a way to measure modularity) and if these SNPs (and loci) are affected by single or multiple environments (a way to measure pleiotropy).

The authors used co-association networks, an innovative approach in this field, to analyse the interaction between the environmental information and the genetic polymorphism of each individual. This methodology is more appropriate than other multivariate methods - such as analysis based on principal components - because it is possible to cluster SNPs based on associations with similar environmental variables. In this sense, the co-association networks allowed to both study the genetic and physical linkage between different co-associations modules but also to compare two different models of evolution: a Modular environmental response architecture (specific genes are affected by specific aspects of the environment) or a Universal pleiotropic environmental response architecture (all genes are affected by all aspects of the environment). The representation of different correlations between allelic frequency and environmental factors (named galaxy biplots) are especially informative to understand the effect of the different clusters on specific aspects of the environment (for example, the co-association network ‘Aridity’ shows strong associations with hot/wet versus cold/dry environments).

The analysis performed by Lotterhos et al. [3], although it has some unavoidable limitations (e.g., only extreme candidate SNPs are selected, limiting the results to the stronger effects; the genetic and physical map is incomplete in this species), includes relevant results and also implements new methodologies in the field. To highlight some of them: the preponderance of a Modular environmental response architecture (evolution in separated modules), the detection of physical linkage among SNPs that are co-associated with different aspects of the environment (which was unexpected a priori), the implementation of co-association networks and galaxy biplots to see the effect of modularity and pleiotropy on different aspects of environment. Finally, this work contains remarkable introductory Figures and Tables explaining unambiguously the main concepts [6] included in this study. This work can be treated as a starting point for many other future studies in the field.

References

[1] Hancock AM, Brachi B, Faure N, Horton MW, Jarymowycz LB, Sperone FG, Toomajian C, Roux F & Bergelson J. 2011. Adaptation to climate across the Arabidopsis thaliana genome. Science 334: 83–86. doi: 10.1126/science.1209244
[2] Wagner GP & Zhang J. The pleiotropic structure of the genotype­phenotype map: the evolvability of complex organisms. Nature Review Genetics 12: 204–213. doi: 10.1038/nrg2949
[3] Lotterhos KE, Yeaman S, Degner J, Aitken S, Hodgins K. 2018. Modularity of genes involved in local adaptation to climate despite physical linkage. bioRxiv 202481, ver. 4 peer-reviewed by Peer Community In Evolutionary Biology. doi: 10.1101/202481
[4] Lotterhos KE & Whitlock MC. 2014. Evaluation of demographic history and neutral parameterization on the performance of FST outlier tests. Molecular Ecology 23: 2178–2192. doi: 10.1111/mec.12725
[5] Lotterhos KE & Whitlock MC. 2015. The relative power of genome scans to detect local adaptation depends on sampling design and statistical method. Molecular Ecology 24: 1031–1046. doi: 10.1111/mec.13100
[6] Paaby AB & Rockman MV. 2013. The many faces of pleiotropy. Trends in Genetics 29: 66-73. doi: 10.1016/j.tig.2012.10.010

The introduction of Ludwigia species as ornamental plants in both North America and Europe dates back almost two centuries, during which time they expanded as naturalized and later invasive species in these territories (Dandelot et al. 2005, Okada et al. 2009). Repeated deliberate or non-deliberate introductions over time of this species complex that can hybridize has given rise to an evolutionarily complex scenario, which is compounded by the difficulty in delimiting some of these species and by the diversity of their modes of reproduction.

Dandelot (2004) and Dandelot et al. (2005) determined the presence of two Ludwigia taxa in France, L. peploides subsp. montevidensis (Spreng.) P.H.Raven (here after Lpm), and L. grandiflora subsp. hexapetala (Hook. & Arn.) G.L.Nesom & Kartesz (here after Lgh), based on their cytotypes (2n = 16, and 2n = 80, respectively) and without evidence of hybridization between them. Furthermore, despite a predominantly vegetative reproduction observed for both species, they differed in their breeding systems. While Lpm is self-compatible and produce a high number of viable seeds in all populations, Lgh is self-incompatible and its populations may drastically differ in seed viability (Dandelot 2004). Several years later, Portillo-Lemus et al. (2021) determined that the differences in seed production between some populations of Lgh are due to the existence of a heteromorphic reproductive system in this taxon, involving a self-incompatible morph (long-style morph; hereafter L-morph), and a self-compatible morph (short-style morph: hereafter S-morph). Moreover, Portillo-Lemus et al. (2022) observed that self-pollen in the L-morph flowers stop growing lately (i.e., in the ovarian area) without fertilizing the ovules, concluding that a late-acting self-incompatible system (hereafter, LSI) is present in this morph.

At this point, it is relevant to understand the possible interactions between populations of different morphs in Lgh, and the implications that they may have on their expansive success in non-native areas in order to develop more effective management plans. To achieve this goal, Stoeckel et al. (2024) analyzed the population genetics in 53 Lgh populations in western Europe (without finding any Lpm population in the sampling area), 40 of which exclusively presented the L-morph and 13 the S-morph. This fact offered the opportunity to compare and interpret the differences between populations of different morphs in Lgh. Other previous works on genetic diversity of Lgh in peripatric or non-native areas pointed to a high clonality and an extremely low genetic diversity (Okada et al. 2009, Armitage et al. 2013), concluding in a monoclonal or few ancestral original clones for these invasive populations.

However, the investigations of Stoeckel et al. (2024) found a high genetic diversity in all populations of Lgh studied despite their predominant clonal reproduction. Interestingly, they found that sexual reproduction is also present, not only in the S-morph by selfing, but also in the L-morph, although limited and preferably by allogamy. They discuss the advantages and drawbacks of the different modes of reproduction observed in Lgh populations, the interactions among them, and the implications that both, the scarcely documented LSI (Gibbs 2014) and selfing, have in the reproductive success and in the maintenance of the high genetic diversity observed in Lgh in western Europe.

The contrasting results with the previous ones (Okada et al. 2009, Armitage et al. 2013) stress the relevance of using appropriate markers and analyses to assess the genetic diversity in autoployploid species, as well as the necessity of knowing the modes of reproduction in the populations studied for an optimal interpretation of the genetic metrics. The approach of the study by Stoeckel et al. (2024) had the challenge of having found suitable markers to deal with a taxon of complex origin such as Lgh, whose genome is compound by a set of autotetraploid chromosomes shared with Lpm and traces of ancient hybridizations of other diploid lineages (Barloy et al. 2024). Using RAD-Seq, Stoeckel et al. (2024) generated an original set of 36 polymorphic SNPs shared between Lgh and Lpm ensuring that these SNPs belong to the tetraploid part of the Lgh genome derived from Lpm. Another interesting contribution of this work is the exhaustive analysis of several genetic descriptors (indexes) and the interpretative guide they provide for each of them in relation to the different modes of reproduction of the study system. Finally, they propose a pair of very useful synthetic indices (i.e., clonality index and selfing index), since they allow to classify populations according to their levels of clonality and selfing.

Stoeckel et al. (2024) conclude the relevance that selfing and LSI populations, and the hybridization between them may have on the expansion and success of invasive plant species, and the necessity to know the modes of reproduction of these populations jointly with their genetic diversity in order to develop appropriate management plans. This study raises new questions such as the modes of reproduction and genetic diversity and structure have other Lgh populations, both invasive and native, and the dynamics of these populations under different future scenarios.

References

Armitage, J. D., Könyves, K., Bailey, J. P., David, J. C., & Culham, A. (2013). A molecular, morphological and cytological investigation of the identity of non-native Ludwigia (Onagraceae) populations in Britain. New Journal of Botany, 3(2), 88–95. https://doi.org/10.1179/2042349713Y.0000000023

Barloy, D., Lemus, L. P.-, Krueger-Hadfield, S. A., Huteau, V., & Coriton, O. (2024). Genomic relationships among diploid and polyploid species of the genus Ludwigia L. section Jussiaea using a combination of cytogenetic, morphological, and crossing investigations. ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2023.01.02.522458

Dandelot, S. (2004). Les Ludwigia spp. invasives du Sud de la France: Historique, Biosystématique, Biologie et Ecologie [PhD thesis, Aix-Marseille 3]. https://www.theses.fr/2004AIX30052

Dandelot, S., Verlaque, R., Dutartre, A., & Cazaubon, A. (2005). Ecological, dynamic and taxonomic problems due to Ludwigia (Onagraceae) in France. Hydrobiologia, 551(1), 131–136. https://doi.org/10.1007/s10750-005-4455-0

Gibbs, P. E. (2014). Late-acting self-incompatibility – the pariah breeding system in flowering plants. New Phytologist, 203(3), 717–734. https://doi.org/10.1111/nph.12874

Okada, M., Grewell, B. J., & Jasieniuk, M. (2009). Clonal spread of invasive Ludwigia hexapetala and L. grandiflora in freshwater wetlands of California. Aquatic Botany, 91(3), 123–129. https://doi.org/10.1016/j.aquabot.2009.03.006

Portillo Lemus, L. O., Bozec, M., Harang, M., Coudreuse, J., Haury, J., Stoeckel, S., & Barloy, D. (2021). Self-incompatibility limits sexual reproduction rather than environmental conditions in an invasive water primrose. Plant-Environment Interactions, 2(2), 74–86. https://doi.org/10.1002/pei3.10042

Portillo Lemus, L. O., Harang, M., Bozec, M., Haury, J., Stoeckel, S., & Barloy, D. (2022). Late-acting self-incompatible system, preferential allogamy and delayed selfing in the heteromorphic invasive populations of Ludwigia grandiflora subsp. hexapetala. Peer Community Journal, 2. https://doi.org/10.24072/pcjournal.108

Stoeckel, S., Becheler, R., Portillo-Lemus, L., Harang, M., Besnard, A.-L., Lassalle, G., Causse-Védrines, R., Michon-Coudouel, S., Park D. J., Pope, B. J., Petit, E. J. & Barloy, D. (2024) Reproductive modes in populations of late-acting self-incompatible and self-compatible polyploid Ludwigia grandiflora subsp. hexapetala in western Europe. biorxiv, ver.4 peer-reviewed and recommended by PCI Evol Biol https://doi.org/10.1101/2024.03.21.586104

Genetic markers are used for in modern population genetics/genomics to uncover the past neutral and selective history of population and species. Besides Single Nucleotide Polymorphisms (SNPs) obtained from whole genome data, microsatellites (or Short Tandem Repeats, SSR) have been common markers of choice in numerous population genetics studies of non-model species with large sample sizes [1]. Microsatellites can be used to uncover and draw inference of the past population demography (e.g. expansion, decline, bottlenecks…), population split, population structure and gene flow, but also life history traits and modes of reproduction (e.g. [2,3]). These markers are widely used in conservation genetics [4] or to study parasites or disease vectors [5]. Microsatellites do show higher mutation rate than SNPs increasing, on the one hand, the statistical power to infer recent events (for example crop domestication, [2,3]), while, on the other hand, decreasing their statistical power over longer time scales due to homoplasy [6].
To perform such analyses, however, an excellent and reliable quality of data is required. As emphasized in the article by De Meeûs et al. [7] three main issues do bias the observed heterozygosity at microsatellites: null alleles, short allele dominance (SAD) and stuttering. These originates from poor PCR amplification. As a result, an excess of homozygosity is observed at the microsatellite loci leading to overestimation of the variation statistics FIS and FST as well as increased linage disequilibrium (LD). For null alleles, several methods and software do help to reduce the bias, and in the present study, De Meeûs et al. [7] propose a way to tackle issues with SAD and stuttering.
The authors study a dataset consisting of 387 samples from 61 subsamples genotyped at nine loci of the species Ixodes scapularis, i.e. ticks transmitting the Lyme disease. Based on correlation methods and FST, FIS they can uncover null alleles and SAD. Stuttering is detected by evaluating the heterozygote deficit between alleles displaying a single repeat difference. Without correction, six loci are affected by one of these amplification problems generating a large deficit of heterozygotes (measured by significant FIS and FST) remaining so after correction for the false discovery rate (FDR). These results would be classically interpreted as a strong Wahlund effect and/or selection at several loci.
After correcting for null alleles, the authors apply two novel corrections: 1) a re-examination of the chromatograms reveals previously disregarded larger alleles thus decreasing SAD, and 2) pooling alleles close in size decreasing stuttering. The corrected dataset shows then a significant excess of heterozygotes as could be expected in a dioecious species with strong population structure. The FDR correction removes then the significant excess of homozygotes and LD between pairs of loci. FST on the cured dataset is used to demonstrate the strong population structure and small effective subpopulation sizes. This is confirmed by a clustering analysis using discriminant analysis of principal components (DAPC).
While based on a specific dataset of ticks from different populations sampled across the USA, the generality of the authors’ approach is presented in Figure 6 in which they provide a step by step flowchart to cure microsatellite datasets from null alleles, SAD and stuttering. Several criteria based on FIS, FST and LD between loci are used as decision keys in the flowchart. An excel file is also provided as help for the curation steps. This study and the proposed methodology are thus extremely useful for all population geneticists working on non-model species with large number of samples genotyped at microsatellite markers. The method not only allows more accurate estimates of heterozygosity but also prevents the thinning of datasets due to the removal of problematic loci. As a follow-up and extension of this work, an exhaustive simulation study could investigate the influence of these data quality issues on past demographic and population structure inference under a wide range of scenarios. This would allow to quantify the current biases in the literature and the robustness of the methodology devised by De Meeûs et al. [7].

References

[1] Jarne, P., and Lagoda, P. J. (1996). Microsatellites, from molecules to populations and back. Trends in ecology & evolution, 11(10), 424-429. doi: 10.1016/0169-5347(96)10049-5
[2] Cornille, A., Giraud, T., Bellard, C., Tellier, A., Le Cam, B., Smulders, M. J. M., Kleinschmit, J., Roldan-Ruiz, I. and Gladieux, P. (2013). Postglacial recolonization history of the E uropean crabapple (Malus sylvestris M ill.), a wild contributor to the domesticated apple. Molecular Ecology, 22(8), 2249-2263. doi: 10.1111/mec.12231
[3] Parat, F., Schwertfirm, G., Rudolph, U., Miedaner, T., Korzun, V., Bauer, E., Schön C.-C. and Tellier, A. (2016). Geography and end use drive the diversification of worldwide winter rye populations. Molecular ecology, 25(2), 500-514. doi: 10.1111/mec.13495
[4] Broquet, T., Ménard, N., & Petit, E. (2007). Noninvasive population genetics: a review of sample source, diet, fragment length and microsatellite motif effects on amplification success and genotyping error rates. Conservation Genetics, 8(1), 249-260. doi: 10.1007/s10592-006-9146-5
[5] Koffi, M., De Meeûs, T., Séré, M., Bucheton, B., Simo, G., Njiokou, F., Salim, B., Kaboré, J., MacLeod, A., Camara, M., Solano, P., Belem, A. M. G. and Jamonneau, V. (2015). Population genetics and reproductive strategies of African trypanosomes: revisiting available published data. PLoS neglected tropical diseases, 9(10), e0003985. doi: 10.1371/journal.pntd.0003985
[6] Estoup, A., Jarne, P., & Cornuet, J. M. (2002). Homoplasy and mutation model at microsatellite loci and their consequences for population genetics analysis. Molecular ecology, 11(9), 1591-1604. doi: 10.1046/j.1365-294X.2002.01576.x
[7] De Meeûs, T., Chan, C. T., Ludwig, J. M., Tsao, J. I., Patel, J., Bhagatwala, J., and Beati, L. (2019). Deceptive combined effects of short allele dominance and stuttering: an example with Ixodes scapularis, the main vector of Lyme disease in the USA. bioRxiv, 622373, ver. 4 peer-reviewed and recommended by Peer Community In Evolutionary Biology. doi: 10.1101/622373

Despite advances in genomic research, many views of genome evolution are still based on what we know from a handful of species, such as humans. This also applies to our knowledge of sex chromosomes. We've apparently been too much used to the situation in which a highly degenerate Y chromosome coexists with an almost normal X chromosome to be able to fully grasp all the questions implied by this situation. Lately, many more sex chromosomes have been studied in other organisms, such as in plants, and the view is changing radically: there is a large diversity of situations, ranging from young highly divergent sex chromosomes to old ones that are so similar that they're hard to detect. Undoubtedly inspired by these recent findings, a few theoretical studies have been published around 2 years ago that put an entirely new light on the evolution of sex chromosomes. The differences between these models have however remained somewhat difficult to appreciate by non-specialists. 

In particular, the models by Lenormand & Roze (2022) and by Jay et al. (2022) seemed quite similar. Indeed, both rely on the same mechanism for initial recombination suppression: a ``lucky'' inversion, i.e. one with less deleterious mutations than the population average, encompassing the sex-determination locus, is initially selected. However, as it doesn't recombine, it will quickly accumulate deleterious mutations lowering its fitness. And it's at this point the models diverge: according to Lenormand & Roze (2022), nascent dosage compensation not only limits the deleterious effects on fitness by the ongoing degeneration, but it actually opposes recombination restoration as this would lead gene expression away from the optimum that has been reached. On the other hand, in the model by Jay et al. (2022), no additional ingredient is required: they argue that once an inversion had been fixed, reversions that restore recombination are extremely unlikely.

This is what Lenormand & Roze (2024) now call a ``constraint'': in Jay et al.'s model, recombination restoration is impossible for mechanistic reasons. Lenormand & Roze (2024) argue such constraints cannot explain long-term recombination suppression. Instead, a mechanism should evolve to limit the negative fitness effects of recombination arrest, otherwise recombination is either restored, or the population goes extinct due to a dramatic drop in the fitness of the heterogametic sex. These two arguments work together: given the huge fitness cost of the lack of ongoing degeneration of the non-recombining Y, in the absence of compensatory mechanisms, there is a very strong selection for the restoration of recombination, so that even when restoration a priori is orders of magnitude less likely than inversion (leading to recombination suppression), it will eventually happen. 

One way the negative fitness effects of recombination suppression can be limited, is the way the authors propose in their own model: dosage compensation evolves through regulatory evolution right at the start of recombination suppression. This changes our classical, simplistic view that dosage compensation evolves in response to degeneration: rather, Lenormand & Roze (2024) argue, that degeneration can only happen when dosage compensation is effective.

The reasoning is convincing and exposes the difference between the models to readers without a firm background in mathematical modelling. Although Lenormand & Roze (2024) target the "constraint theory", it seems likely that other theories for the maintenance of recombination suppression that don't imply the compensation of early degeneration are subject to the same criticism. Indeed, they mention the widely-cited "sexual antagonism" theory, in which mutations with a positive effect in males but a negative in females will select for recombination suppression that will link them to the sex-determining gene on the Y. However, once degeneration starts, the sexually-antagonistic benefits should be huge to overcome the negative effects of degeneration, and it's unlikely they'll be large enough.

A convincing argument by Lenormand & Roze (2024) is that there are many ways recombination could be restored, allowing to circumvent the possible constraints that might be associated with reverting an inversion. First, reversions don't have to be exact to restore recombination. Second, the sex-determining locus can be transposed to another chromosome pair, or an entirely new sex-determining locus might evolve, leading to sex-chromosome turnover which has effectively been observed in several groups.

These modelling studies raise important questions that need to be addressed with both theoretical and empirical work. First, is the regulatory hypothesis proposed by Lenormand & Roze (2022) the only plausible mechanism for the maintenance of long-term recombination suppression? The female- and male-specific trans regulators of gene expression that are required for this model, are they readily available or do they need to evolve first? Both theoretical work and empirical studies of nascent sex chromosomes will help to answer these questions. However, nascent sex chromosomes are difficult to detect and dosage compensation is difficult to reveal.

Second, how many species today actually have "stable" recombination suppression? Maybe many species are in a transient phase, with different populations having different inversions that are either on their way to being fixed or starting to get counterselected. The models have now shown us some possibilities qualitatively but can they actually be quantified to be able to fit the data and to predict whether an observed case of recombination suppression is transient or stable? 

The debate will continue, and we need the active contribution of theoretical biologists to help clarify the underlying hypotheses of the proposed mechanisms. 

Conflict of interest statement: I did co-author a manuscript with D. Roze in 2023, but do not consider this a conflict of interest. The manuscript is the product of discussions that have taken place in a large consortium mainly in 2019. It furthermore deals with an entirely different topic of evolutionary biology.

References

Jay P, Tezenas E, Véber A, and Giraud T. (2022) Sheltering of deleterious mutations explains the stepwise extension of recombination suppression on sex chromosomes and other supergenes. PLoS Biol.;20:e3001698. https://doi.org/10.1371/journal.pbio.3001698
 
Lenormand T and Roze D. (2022) Y recombination arrest and degeneration in the absence of sexual dimorphism. Science;375:663-6. https://doi.org/10.1126/science.abj1813
 
Lenormand T and Roze D. (2024) Can mechanistic constraints on recombination reestablishment explain the long-term maintenance of degenerate sex chromosomes? bioRxiv, ver. 5 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2023.02.17.528909