Cytoplasmic incompatibility (CI) is a mating incompatibility that is induced by maternally inherited endosymbionts in many arthropods. These endosymbionts include, most famously, the alpha-proteobacterium Wolbachia pipientis (Yen & Barr 1971; Werren et al. 2008) but also the Bacteroidetes bacterium Cardinium hertigii (Zchori-Fein et al. 2001), a gamma-proteobacterium of the genus Rickettsiella (Rosenwald et al. 2020) and another, as yet undescribed alpha-proteobacterium (Takano et al. 2017). CI manifests as embryonic mortality in crosses between infected males and females that are uninfected or infected with a different strain, whereas embryos develop normally in all other crosses. This phenotype may enable the endosymbionts to spread rapidly within their host population. Exploiting this, CI-inducing Wolbachia are being harnessed to control insect-borne diseases (e.g., O'Neill 2018). Much progress elucidating the genetic basis and developmental mechanism of CI has been made in recent years, but many open questions remain (Shropshire et al. 2020).
Immediately following the discovery and early study of CI in mosquitoes, Laven (1959, 1967) proposed that CI could be an important driver of speciation. Indeed, bi-directional CI can strongly reduce gene flow between two populations due to the elimination of F1 embryos, so that CI can act as a trigger for genetic differentiation in the host (Telschow et al. 2002, 2005). This idea has received much attention, and a potential role for CI in incipient speciation has been demonstrated in several species (e.g., Bordenstein et al. 2001; Jaenike et al. 2006). However, we still don’t know how commonly CI actually triggers speciation, rather than being merely a minor player or secondary phenomenon. The problem is that in addition to CI, postzygotic reproductive isolation can also be caused by host-induced, nuclear incompatibilities. Determining the relative contributions of these two causes of isolation is difficult and has rarely been done.
The study by Cruz et al. (2020) addresses this problem head-on, using a study system of Tetranychus urticae spider mites. These cosmopolitan mites are infected with different strains of Wolbachia. They come in two different colour forms (red and green) that can co-occur sympatrically on the same host plant but exhibit various degrees of reproductive isolation. A complicating factor in spider mites is that they are haplodiploid: unfertilised eggs develop into haploid males and are therefore not affected by any postzygotic incompatibilities, whereas fertilised eggs normally develop into diploid females. In haplodiploids, Wolbachia-induced CI can either kill diploid embryos (as in diplodiploid species), or turn them into haploid males. In their study, Cruz et al. used three different populations (one of the green and two of the red form) and employed a full factorial experiment involving all possible combinations of crosses of Wolbachia infected or uninfected males and females. For each cross, they measured F1 embryonic and juvenile mortality as well as sex ratio, and they also measured F1 fertility and F2 viability. Their results showed that there is strong reduction in hybrid female production caused by Wolbachia-induced CI. However, independent of this and through a different mechanism, there is an even stronger reduction in hybrid production caused by host-associated incompatibilities. In combination with the also observed near-complete sterility of F1 hybrid females and full F2 hybrid breakdown (neither of which is caused by Wolbachia), the results indicate essentially complete reproductive isolation between the green and red forms of T. urticae.
Overall, this is an elegant study with an admirably clean and comprehensive experimental design. It demonstrates that Wolbachia can contribute to reproductive isolation between populations, but that host-induced mechanisms of reproductive isolation predominate in these spider mite populations. Further studies in this exiting system would be useful that also investigate the contribution of pre-zygotic isolation mechanisms such as assortative mating, ascertain whether the results can be generalised to other populations, and – most challengingly – establish the order in which the different mechanisms of reproductive isolation evolved.

References

Bordenstein, S. R., O'Hara, F. P., and Werren, J. H. (2001). Wolbachia-induced incompatibility precedes other hybrid incompatibilities in Nasonia. Nature, 409(6821), 707-710. doi: https://doi.org/10.1038/35055543
Cruz, M. A., Magalhães, S., Sucena, É., and Zélé, F. (2020) Wolbachia and host intrinsic reproductive barriers contribute additively to post-mating isolation in spider mites. bioRxiv, 2020.06.29.178699, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: https://doi.org/10.1101/2020.06.29.178699
Jaenike, J., Dyer, K. A., Cornish, C., and Minhas, M. S. (2006). Asymmetrical reinforcement and Wolbachia infection in Drosophila. PLoS Biol, 4(10), e325. doi: https://doi.org/10.1371/journal.pbio.0040325
Laven, H. (1959). SPECIATION IN MOSQUITOES Speciation by Cytoplasmic Isolation in the Culex Pipiens-Complex. In Cold Spring Harbor Symposia on Quantitative Biology (Vol. 24, pp. 166-173). Cold Spring Harbor Laboratory Press.
Laven, H. (1967). A possible model for speciation by cytoplasmic isolation in the Culex pipiens complex. Bulletin of the World Health Organization, 37(2), 263-266.
O’Neill S.L. (2018) The Use of Wolbachia by the World Mosquito Program to Interrupt Transmission of Aedes aegypti Transmitted Viruses. In: Hilgenfeld R., Vasudevan S. (eds) Dengue and Zika: Control and Antiviral Treatment Strategies. Advances in Experimental Medicine and Biology, vol 1062. Springer, Singapore. doi: https://doi.org/10.1007/978-981-10-8727-1_24
Rosenwald, L.C., Sitvarin, M.I. and White, J.A. (2020). Endosymbiotic Rickettsiella causes cytoplasmic incompatibility in a spider host. doi: https://doi.org/10.1098/rspb.2020.1107
Shropshire, J. D., Leigh, B., and Bordenstein, S. R. (2020). Symbiont-mediated cytoplasmic incompatibility: what have we learned in 50 years?. Elife, 9, e61989. doi: https://doi.org/10.7554/eLife.61989
Takano et al. (2017). Unique clade of alphaproteobacterial endosymbionts induces complete cytoplasmic incompatibility in the coconut beetle. Proceedings of the National Academy of Sciences, 114(23), 6110-6115. doi: https://doi.org/10.1073/pnas.1618094114
Telschow, A., Hammerstein, P., and Werren, J. H. (2002). The effect of Wolbachia on genetic divergence between populations: models with two-way migration. the american naturalist, 160(S4), S54-S66. doi: https://doi.org/10.1086/342153
Telschow, A., Hammerstein, P., and Werren, J. H. (2005). The effect of Wolbachia versus genetic incompatibilities on reinforcement and speciation. Evolution, 59(8), 1607-1619. doi: https://doi.org/10.1111/j.0014-3820.2005.tb01812.x
Werren, J. H., Baldo, L., and Clark, M. E. (2008). Wolbachia: master manipulators of invertebrate biology. Nature Reviews Microbiology, 6(10), 741-751. doi: https://doi.org/10.1038/nrmicro1969
Yen, J. H., and Barr, A. R. (1971). New hypothesis of the cause of cytoplasmic incompatibility in Culex pipiens L. Nature, 232(5313), 657-658. doi: https://doi.org/10.1038/232657a0
Zchori-Fein, E., Gottlieb, Y., Kelly, S. E., Brown, J. K., Wilson, J. M., Karr, T. L., and Hunter, M. S. (2001). A newly discovered bacterium associated with parthenogenesis and a change in host selection behavior in parasitoid wasps. Proceedings of the National Academy of Sciences, 98(22), 12555-12560. doi: https://doi.org/10.1073/pnas.221467498

The search for physiological markers of health and survival in wild animal populations is attracting a great deal of interest. At present, there is no (and may never be) consensus on such a single, robust marker but of all the proposed physiological markers, telomere length is undoubtedly the most widely studied in the field of evolutionary ecology (Monaghan et al., 2022). 

Broadly speaking, telomeres are non-coding DNA sequences located at the end of chromosomes in eukaryotes, protecting genomic DNA against oxidative stress and various detrimental processes (e.g. DNA end-joining) and thus maintaining genome stability (Blackburn et al., 2015). However, in most somatic cells from the vast majority of the species, telomere sequences are not replicated and telomere length progressively declines with increased age (Remot et al., 2022). This shortening of telomere length upon a critical level is causally linked to cellular senescence and has been invoked as one of the primary causes of the aging process (López-Otín et al., 2023). Studies performed in both captive and wild populations of animals have further demonstrated that short telomeres (or telomere sequences with a fast attrition rate) are to some extent associated with an increased risk of mortality, even if the magnitude of this association largely differs between species and populations (Wilbourn et al., 2018).

The repeated observations of associations between telomere length and mortality risk have called for studies seeking to identify the ecological and biological factors that – beyond chronological age – shape the between-individual variability in telomere length. A wide spectrum of environmental stressors such as the level of exposure to pathogens or the degree of human disturbances has been proposed as possible modulators of telomere dynamics (see Chatelain et al., 2019). However, within species, the relative contribution of various ecological and biological factors on telomere length has been rarely quantified. In that context, the study of Criscuolo and colleagues (2023) constitutes a timely attempt to decipher the relative contribution of environmental and biological factors driving telomere length in nestlings (i.e. when individuals are between 15 and 35 days of age) from a wild population of little owls, Athene noctua.

In addition to chronological age, Criscuolo and colleagues (2023) analysed the effects of two environmental variables (i.e. cohort and habitat quality) as well as three life history traits (i.e. hatching rank, sex and body condition). Among these traits, sex was found to impact nestling’s telomere length with females carrying longer telomeres than males. Traditionally, the among-individuals variability in telomere length during the juvenile period is interpreted as a direct consequence of differences in growth allocation. Fast-growing individuals are typically supposed to undergo more cell divisions and a higher exposure to oxidative stress, which ultimately shortens telomeres (Monaghan & Ozanne, 2018). Whether - despite a slightly female-biased sexual size dimorphism - male little owls display a condensed period of fast growth that could explain their shorter telomere is yet to be determined. Future studies should also explore the consequences of these sex differences in telomere length in terms of mortality risk. In birds, it has been observed that telomere length during early life can predict lifespan (see Heidinger et al., 2012 in zebra finches, Taeniopygia guttata), suggesting that females little owls might live longer than their conspecific males. Yet, adult mortality is generally female-biased in birds (Liker & Székely, 2005) and whether little owls constitute an exception to this rule - possibly mediated by sex-specific telomere dynamics - remains to be explored.   

Quite surprisingly, the present study in little owls did not evidence any clear effect of environmental conditions on nestling’s telomere length, at both temporal and special scales. While a trend for a temporal effect was detected with telomere length being slightly shorter for nestling born the last year of the study (out of 4 years analysed), habitat quality (measured by the proportion of meadow and orchards in the nest environment) had absolutely no impact on nestling telomere length. Recently published studies in wild populations of vertebrates have highlighted the detrimental effects of harsh environmental conditions on telomere length (e.g. Dupoué et al., 2022 in common lizards, Zootoca vivipara), arguing for a key role of telomere dynamics in the emerging field of conservation physiology. While we can recognize the relevance of such an integrative approach, especially in the current context of climate change, the study by Criscuolo and colleagues (2023) reminds us that the relationships between environmental conditions and telomere dynamics are far from straightforward. Depending on the species and its life history, telomere length in early life could indeed capture very different environmental signals.

References

Blackburn, E. H., Epel, E. S., & Lin, J. (2015). Human telomere biology: A contributory and interactive factor in aging, disease risks, and protection. Science, 350(6265), 1193-1198. https://doi.org/10.1126/science.aab3389
 
Chatelain, M., Drobniak, S. M., & Szulkin, M. (2019). The association between stressors and telomeres in non-human vertebrates: A meta-analysis. Ecology Letters, 23, 381-398. https://doi.org/10.1111/ele.13426
 
Criscuolo, F., Fache, I., Scaar, B., Zahn, S. & Bleu, J. (2023). Telomere length vary with sex, hatching rank and year of birth in little owls, Athene noctua. EcoEvoRxiv, ver.4, peer-reviewed and recommended by PCI Evol Biol. https://doi.org/10.32942/X2BS3S
 
Dupoué, A., Blaimont, P., Angelier, F., Ribout, C., Rozen-Rechels, D., Richard, M., & Le Galliard, J. F. (2022). Lizards from warm and declining populations are born with extremely short telomeres. Proceedings of the National Academy of Sciences, 119(33), 2201371119. https://doi.org/10.1073/pnas.2201371119
 
Heidinger, B. J., Blount, J. D., Boner, W., Griffiths, K., Metcalfe, N. B., & Monaghan, P. (2012). Telomere length in early life predicts lifespan. Proceedings of the National Academy of Sciences, 109(5), 1743-1748. https://doi.org/10.1073/pnas.1113306109
 
Liker, A., & Székely, T. (2005). Mortality costs of sexual selection and parental care in natural populations of birds. Evolution, 59(4), 890-897. https://doi.org/10.1111/j.0014-3820.2005.tb01762.x
 
López-Otín, C., Blasco, M. A., Partridge, L., Serrano, M., & Kroemer, G. (2023). Hallmarks of aging: An expanding universe. Cell, 186(2), 243-278. https://doi.org/10.1016/j.cell.2022.11.001
 
Monaghan, P., Olsson, M., Richardson, D. S., Verhulst, S., & Rogers, S. M. (2022). Integrating telomere biology into the ecology and evolution of natural populations: Progress and prospects. Molecular Ecology, 31(23), 5909-5916. https://doi.org/10.1111/mec.16768
 
Monaghan, P., & Ozanne, S. E. (2018). Somatic growth and telomere dynamics in vertebrates: Relationships, mechanisms and consequences. Phil. Trans. R. Soc. B, 373(1741), 20160446. https://doi.org/10.1098/rstb.2016.0446
 
Remot, F., Ronget, V., Froy, H., Rey, B., Gaillard, J., Nussey, D. H., & Lemaitre, J. (2022). Decline in telomere length with increasing age across nonhuman vertebrates: A meta‐analysis. Molecular Ecology, 31(23), 5917-5932. https://doi.org/10.1111/mec.16145
 
Wilbourn, R. V., Moatt, J. P., Froy, H., Walling, C. A., Nussey, D. H., & Boonekamp, J. J. (2018). The relationship between telomere length and mortality risk in non-model vertebrate systems: A meta-analysis. Phil. Trans. R. Soc. B, 373(1741), 20160447. https://doi.org/10.1098/rstb.2016.0447

The Latitudinal Diversity Gradient (LDG) has fascinated natural historians, ecologists and evolutionary biologists ever since [1] described it about 200 years ago [2]. Despite such interest, agreement on the origin and nature of this gradient has been elusive. Several tens of hypotheses and models have been put forward as explanations for the LDG [2-3], that can be grouped in ecological, evolutionary and historical explanations [4] (see also [5]). These explanations can be reduced to no less than 26 hypotheses, which account for variations in ecological limits for the establishment of progressively larger assemblages, diversification rates, and time for species accumulation [5]. Besides that, although in general the tropics hold more species, different taxa show different shapes and rates of spatial variation [6], and a considerable number of groups show reverse patterns, with richer assemblages in cold temperate regions (see e.g. [7-9]).
Understanding such complexity needs integrating ecological and evolutionary research into the wide temporal and spatial perspectives provided by the burgeoning field of biogeography. This integrative discipline ¬–that traces back to Humboldt himself (e.g. [10])– seeks to put together historical and functional explanations to explain the complex dynamics of Earth’s biodiversity. Different to quantum physicists, biogeographers cannot pursue the ultimate principle behind the patterns we observe in nature due to the interplay of causes and effects, which in fact tell us that there is not such a single principle. Rather, they need to identify an array of basic principles coming from different perspectives, to then integrate them into models that provide realistic –but never simple– explanations to biodiversity gradients such as LDG (see, e.g., [5; 11]). That is, rather than searching for a sole explanation, research on the LDG must aim to identify as many signals hidden in the pattern as possible, and provide hypotheses or models that account for these signals. To later integrate them and, whenever possible, to validate them with empirical data on the organisms’ distribution, ecology and traits, phylogenies, fossils, etc.
Within this context, Meseguer & Condamine [12] provide a novel perspective to LDG research using phylogenetic and fossil evidence on the origin and extinction of taxa within the turtle, crocodile and lizard (i.e. squamate) lineages. By digging into deep time down to the Triassic (about 250 Myr ago) they are able to identify several episodes of flattening and steepening of the LDG for these three clades. Strikingly, their results show similar diversification rates in the northern hemisphere and in the equator during the over 100 Myr long global greenhouse period that extends from the late Jurassic to the Cretaceous and early Neogene. During this period, the LDG for these three groups would have appeared quite even across a mainly tropical Globe, although the equatorial regions were apparently much more evolutionarily dynamic. The equator shows much higher rates of origination and extinction of branches throughout the Cretaceous, but they counteract each other so net diversification is similar to that of the northern hemisphere in all three groups. The transition to a progressively colder Earth in the Paleogene (starting around 50 Myr ago) provokes a mass extinction in the three clades, which is compensated in the equator by the dispersal of many taxa from the areas that currently pertain to the Holarctic biogeographical realm. Finally, during the coldhouse Earth’s climatic conditions of the Neogene only squamates show significant positive diversification rates in extratropical areas, while the diversity of testudines remains, and crocodiles continue declining progressively towards oblivion in the whole world.
Meseguer & Condamine [12] attribute these temporal patterns to the so-called asymmetric gradient of extinction and dispersal (AGED) framework. Here, the dynamics of extinction-at and dispersal-from high latitudes during colder periods increase the steepness of the LDG. Whereas the gradient flattens when Earth warms up as a result of dispersal from the equator followed by increased diversification in extratropical regions. This idea in itself is not new, for the influence of climatic oscillations on diversification rates is well known, at least for the Pleistocene Ice Ages [13], as is the effect of niche conservatism on the LDG [14]. Nevertheless, Meseguer & Condamine’s AGED provides a synthetic verbal model that could allow integrating the three main types of processes behind the LDG into a single framework. To do this it would be necessary to combine AGED’s cycles of dispersal and diversification with realistic models of: (1) the ecological limits to host rich assemblages in the colder and less productive temperate climatic domains; (2) the variations in diversification rates with shifts in temperature and/or energy regimes; and (3) the geographical patterns of climatic oscillation through time that determine the time for species accumulation in each region.
Integrating these models may allow transposing Meseguer & Condamine’s [12] framework into the more mechanistic macroecological models advocated by Pontarp et al. [5]. This type of mechanistic models has been already used to understand the development of biodiversity gradients through the climatic oscillations of the Pleistocene and the Quaternary (e.g. [11]). So the challenge in this case would be to generate a realistic scenario of geographical dynamics that accounts for plate tectonics and long-term climatic oscillations. This is still a major gap and we would benefit from the integrated work by historical geologists and climatologists here. For instance, there is little doubt about the progressive cooling through the Cenozoic based in isotope recording in sea floor sediments [15]. Meseguer & Condamine [12] use this evidence for separating greenhouse, transition and coldhouse world scenarios, which should not be a problem for these rough classes. However, a detailed study of the evolutionary correlation of true climate variables across the tree of life is still pending, as temperature is inferred only for sea water in an ice-free ocean, say earlier half of the Cenozoic [15]. Precipitation regime is even less known. Such scenario would provide a scaffold upon which the temporal dynamics of several aspects of the generation and loss of biodiversity can be modelled. Additionally, one of the great advantages of selecting key clades to study the LDG would be to determine the functional basis of diversification. There are species traits that are well known to affect speciation and extinction probabilities, such as reproductive strategies or life histories (e.g. [16]). Whereas these traits might also be a somewhat redundant effect of climatic causes, they might foster (i.e. “extended reinforcement”, [17]) or slow diversification. Even so, it is unlikely that such a model would account for all the latitudinal variation in species richness. But it will at least provide a baseline for the main latitudinal variations in the diversity of the regional communities (sensu [18]) worldwide. Within this context the effects of recent ecological, evolutionary and historical processes, such as environmental heterogeneity, current diversification rates or glacial cycles, will only modify the general LDG pattern resulting from the main processes contained in Meseguer & Condamine’s AGED, thereby providing a more comprehensive understanding of the geographical gradients of diversity.

References
[1] Humboldt, A. v. (1808). Ansichten der Natur, mit wissenschaftlichen Erläuterungen. J. G. Cotta, Tübingen.
[2] Hawkins, B. A. (2001). Ecology's oldest pattern? Trends in Ecology & Evolution, 16, 470. doi: 10.1016/S0169-5347(01)02197-8
[3] Lomolino, M. V., Riddle, B. R. & Whittaker, R. J. (2017). Biogeography. Fifth Edition. Sinauer Associates, Inc., Sunderland, Massachussets.
[4] Mittelbach, G. G., Schemske, D. W., Cornell, H. V., Allen, A. P., Brown, J. M., Bush, M. B., Harrison, S. P., Hurlbert, A. H., Knowlton, N., Lessios, H. A., McCain, C. M., McCune, A. R., McDade, L. A., McPeek, M. A., Near, T. J., Price, T. D., Ricklefs, R. E., Roy, K., Sax, D. F., Schluter, D., Sobel, J. M. & Turelli, M. (2007). Evolution and the latitudinal diversity gradient: speciation, extinction and biogeography. Ecology Letters, 10, 315-331. doi: 10.1111/j.1461-0248.2007.01020.x
[5] Pontarp, M., Bunnefeld L., Cabral, J. S., Etienne, R. S., Fritz, S. A., Gillespie, R. Graham, C. H., Hagen, O., Hartig, F., Huang, S., Jansson, R., Maliet, O., Münkemüller, T., Pellissier, L., Rangel, T. F., Storch, D., Wiegand, T. & Hurlbert, A. H. (2019). The latitudinal diversity gradient: novel understanding through mechanistic eco-evolutionary models. Trends in ecology & evolution, 34, 211-223. doi: 10.1016/j.tree.2018.11.009
[6] Hillebrand, H. (2004). On the generality of the latitudinal diversity gradient. The American Naturalist, 163, 192-211. doi: 10.1086/381004
[7] Santos, A. M. C. & Quicke, D. L. J. (2011). Large-scale diversity patterns of parasitoid insects. Entomological Science, 14, 371-382. doi: 10.1111/j.1479-8298.2011.00481.x
[8] Morinière, J., Van Dam, M. H., Hawlitschek, O., Bergsten, J., Michat, M. C., Hendrich, L., Ribera, I., Toussaint, E. F. A. & Balke, M. (2016). Phylogenetic niche conservatism explains an inverse latitudinal diversity gradient in freshwater arthropods. Scientific Reports, 6, 26340. doi: 10.1038/srep26340
[9] Weiser, M. D., Swenson, N. G., Enquist, B. J., Michaletz, S. T., Waide, R. B., Zhou, J. & Kaspari, M. (2018). Taxonomic decomposition of the latitudinal gradient in species diversity of North American floras. Journal of Biogeography, 45, 418-428. doi: 10.1111/jbi.13131
[10] Humboldt, A. v. (1805). Essai sur la geographie des plantes; accompagné d'un tableau physique des régions equinoxiales. Levrault, Paris.
[11] Rangel, T. F., Edwards, N. R., Holden, P. B., Diniz-Filho, J. A. F., Gosling, W. D., Coelho, M. T. P., Cassemiro, F. A. S., Rahbek, C. & Colwell, R. K. (2018). Modeling the ecology and evolution of biodiversity: Biogeographical cradles, museums, and graves. Science, 361, eaar5452. doi: 10.1126/science.aar5452
[12] Meseguer, A. S. & Condamine, F. L. (2019). Ancient tropical extinctions contributed to the latitudinal diversity gradient. bioRxiv, 236646, ver. 4 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/236646
[13] Jansson, R., & Dynesius, M. (2002). The fate of clades in a world of recurrent climatic change: Milankovitch oscillations and evolution. Annual review of ecology and systematics, 33(1), 741-777. doi: 10.1146/annurev.ecolsys.33.010802.150520
[14] Wiens, J. J., & Donoghue, M. J. (2004). Historical biogeography, ecology and species richness. Trends in ecology & evolution, 19, 639-644. doi: 10.1016/j.tree.2004.09.011
[15] Zachos, J. C., Dickens, G. R., & Zeebe, R. E. (2008). An early Cenozoic perspective on greenhouse warming and carbon-cycle dynamics. Nature, 451, 279-283. doi: 10.1038/nature06588
[16] Zúñiga-Vega, J. J., Fuentes-G, J. A., Ossip-Drahos, A. G., & Martins, E. P. (2016). Repeated evolution of viviparity in phrynosomatid lizards constrained interspecific diversification in some life-history traits. Biology letters, 12, 20160653. doi: 10.1098/rsbl.2016.0653
[17] Butlin, R. K., & Smadja, C. M. (2018). Coupling, reinforcement, and speciation. The American Naturalist, 191, 155-172. doi: 10.1086/695136
[18] Ricklefs, R. E. (2015). Intrinsic dynamics of the regional community. Ecology letters, 18, 497-503. doi: 10.1111/ele.12431

Plants display an amazing diversity of reproductive strategies with and without sex. This diversity is particularly remarkable in flowering plants, as highlighted by Charles Darwin, who wrote several botanical books scrutinizing plant reproduction. One particularly influential work concerned floral variation [1]. Darwin recognized that flowers may present different forms within a single population, with or without sex specialization. The number of species concerned is small, but they display recurrent patterns, which made it possible for Darwin to invoke natural and sexual selection to explain them. Most of early evolutionary theory on the evolution of reproductive strategies was developed in the first half of the 20th century and was based on animals. However, the pioneering work by David Lloyd from the 1970s onwards excited interest in the diversity of plant sexual strategies as models for testing adaptive hypotheses and predicting reproductive outcomes [2]. The sex specialization of individual flowers and plants has since become one of the favorite topics of evolutionary biologists. However, attention has focused mostly on cases related to sex differentiation (dioecy and associated conditions [3]). Separate unisexual flower types on the same plant (monoecy and related cases, rendering the plant functionally hermaphroditic) have been much less studied, apart from their possible role in the evolution of dioecy [4] or their association with particular modes of pollination [5].
Two specific non-mutually exclusive hypotheses on the evolution of separate sexes in flowers (dicliny) have been proposed, both anchored in Lloyd’s views and Darwin’s legacy, with selfing avoidance and optimal limited resource allocation. Intermediate sex separation, in which sex morphs have different combinations of unisexual and hermaphrodite flowers, has been crucial for testing these hypotheses through comparative analyses of optimal conditions in suggested transitions. Again, cases in which floral unisexuality does not lead to sex separation have been studied much less than dioecious plants, at both the microevolutionary and macroevolutionary levels. It is surprising that the increasing availability of plant phylogenies and powerful methods for testing evolutionary transitions and correlations have not led to more studies, even though the frequency of monoecy is probably highest among diclinous species (those with unisexual flowers in any distribution among plants within a population [6]).
The study by Torices et al. [7] aims to fill this gap, offering a different perspective to that provided by Diggle & Miller [8] on the evolution of monoecious conditions. The authors use heads of a number of species of the sunflower family (Asteraceae) to test specifically the effect of resource limitation on the expression of sexual morphs within the head. They make use of the very particular and constant architecture of inflorescences in these species (the flower head or “capitulum”) and the diversity of sexual conditions (hermaphrodite, gynomonoecious, monoecious) and their spatial pattern within the flower head in this plant family to develop an elegant means of testing this hypothesis. Their results are consistent with their expectations on the effect of resource limitation on the head, as determined by patterns of fruit size within the head, assuming that female fecundity is more strongly limited by resource availability than male function.
The authors took on a huge challenge in choosing to study the largest plant family (about 25 thousand species). Their sample was limited to only about a hundred species, but species selection was very careful, to ensure that the range of sex conditions and the available phylogenetic information were adequately represented. The analytical methods are robust and cast no doubt on the reported results. However, I can’t help but wonder what would happen if the antiselfing hypothesis was tested simultaneously. This would require self-incompatibility (SI) data for the species sample, as the presence of SI is usually invoked as a powerful antiselfing mechanism, rendering the unisexuality of flowers unnecessary. However, SI is variable and frequently lost in the sunflower family [9]. I also wonder to what extent the very specific architecture of flower heads imposes an idiosyncratic resource distribution that may have fixed these sexual systems in species and lineages of the family. Although not approached in this study, intraspecific variation seems to be low. It would be very interesting to use similar approaches in other plant groups in which inflorescence architecture is lax and resource distribution may differ. A whole-plant approach might be required, rather than investigations of single inflorescences as in this study. This study has no flaws, but instead paves the way for further testing of a long-standing dual hypothesis, probably with different outcomes in different ecological and evolutionary settings. In the end, sex is not only in the head.

References

[1] Darwin, C. (1877). The different forms of flowers on plants of the same species. John Murray.
[2] Barrett, S. C. H., and Harder, L. D. (2006). David G. Lloyd and the evolution of floral biology: from natural history to strategic analysis. In L.D. Harder, L. D., and Barrett, S. C. H. (eds) Ecology and Evolution of Flowers. OUP, Oxford. Pp 1-21.
[3] Geber, M. A., Dawson, T. E., and Delph, L. F. (eds) (1999). Gender and sexual dimorphism in flowering plants. Springer, Berlin.
[4] Charlesworth, D. (1999). Theories of the evolution of dioecy. In Geber, M. A., Dawson T. E. and Delph L. F. (eds) (1999). Gender and sexual dimorphism in flowering plants. Springer, Berlin. Pp. 33-60.
[5] Friedman, J., and Barrett, S. C. (2008). A phylogenetic analysis of the evolution of wind pollination in the angiosperms. International Journal of Plant Sciences, 169(1), 49-58. doi: 10.1086/523365
[6] Renner, S. S. (2014). The relative and absolute frequencies of angiosperm sexual systems: dioecy, monoecy, gynodioecy, and an updated online database. American Journal of botany, 101(10), 1588-1596. doi: 10.3732/ajb.1400196
[7] Torices, R., Afonso, A., Anderberg, A. A., Gómez, J. M., and Méndez, M. (2019). Architectural traits constrain the evolution of unisexual flowers and sexual segregation within inflorescences: an interspecific approach. bioRxiv, 356147, ver. 3 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/356147
[8] Diggle, P. K., and Miller, J. S. (2013). Developmental plasticity, genetic assimilation, and the evolutionary diversification of sexual expression in Solanum. American journal of botany, 100(6), 1050-1060. doi: 10.3732/ajb.1200647
[9] Ferrer, M. M., and Good‐Avila, S. V. (2007). Macrophylogenetic analyses of the gain and loss of self‐incompatibility in the Asteraceae. New Phytologist, 173(2), 401-414. doi: 10.1111/j.1469-8137.2006.01905.x

Modelling the evolution of complete genome sequences in populations requires accounting for the recombination process, as a single tree can no longer describe the underlying genealogy. The sequentially Markov coalescent (SMC, McVean and Cardin 2005; Marjoram and Wall 2006) approximates the standard coalescent with recombination process and permits estimating population genetic parameters (e.g., population sizes, recombination rates) using population genomic datasets. As such datasets become available for an increasing number of species, more fine-tuned models are needed to encompass the diversity of life cycles of organisms beyond the model species on which most methods have been benchmarked.

The work by Korfmann et al. (Korfmann et al. 2024) represents a significant step forward as it accounts for multiple mergers in SMC models. Multiple merger models account for simultaneous coalescence events so that more than two lineages find a common ancestor in a given generation. This feature is not allowed in standard coalescent models and may result from selection or skewed offspring distributions, conditions likely met by a broad range of species, particularly microbial.

Yet, this work goes beyond extending the SMC, as it introduces several methodological innovations. The "classical" SMC-based inference approaches rely on hidden Markov models to compute the likelihood of the data while efficiently integrating over the possible ancestral recombination graphs (ARG). Following other recent works (e.g. Gattepaille et al. 2016), Korfmann et al. propose to separate the ARG inference from model parameter estimation under maximum likelihood (ML). They introduce a procedure where the ARG is first reconstructed from the data and then taken as input in the model fitting step. While this approach does not permit accounting for the uncertainty in the ARG reconstruction (which is typically large), it potentially allows for the extraction of more information from the ARG, such as the occurrence of multiple merging events. Going away from maximum likelihood inference, the authors trained a graph neural network (GNN) on simulated ARGs, introducing a new, flexible way to estimate population genomic parameters.

The authors used simulations under a beta-coalescent model with diverse demographic scenarios and showed that the ML and GNN approaches introduced can reliably recover the simulated parameter values. They further show that when the true ARG is given as input, the GNN outperforms the ML approach, demonstrating its promising power as ARG reconstruction methods improve. In particular, they showed that trained GNNs can disentangle the effects of selective sweeps and skewed offspring distributions while inferring past population size changes.

This work paves the way for new, exciting applications, though many questions must be answered. How frequent are multiple mergers? As the authors showed that these events "erase" the record of past demographic events, how many genomes are needed to conduct reliable inference, and can the methods computationally cope with the resulting (potentially large) amounts of required data? This is particularly intriguing as micro-organisms, prone to strong selection and skewed offspring distributions, also tend to carry smaller genomes.

References

Gattepaille L, Günther T, Jakobsson M. 2016. Inferring Past Effective Population Size from Distributions of Coalescent Times. Genetics 204:1191-1206.
https://doi.org/10.1534/genetics.115.185058
 
Korfmann K, Sellinger T, Freund F, Fumagalli M, Tellier A. 2024. Simultaneous Inference of Past Demography and Selection from the Ancestral Recombination Graph under the Beta Coalescent. bioRxiv, 2022.09.28.508873. ver. 5 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.09.28.508873
 
Marjoram P, Wall JD. 2006. Fast "coalescent" simulation. BMC Genet. 7:16.
https://doi.org/10.1186/1471-2156-7-16
 
McVean GAT, Cardin NJ. 2005. Approximating the coalescent with recombination. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 360:1387-1393.
https://doi.org/10.1098/rstb.2005.1673