In most animal species, males and females display distinct survival prospect, a phenomenon known as sex gap in longevity (SGL, Marais et al. 2018). The study of SGLs is crucial not only for having a full picture of the causes underlying organisms’ health, aging and death but also to initiate the development of sex-specific anti-aging interventions in humans (Austad and Bartke 2015). Three non-mutually evolutionary causes have been proposed to underlie SGLs (Marais et al. 2018). First, SGLs could be the consequences of sex-differences in life history strategies. For example, evolving dimorphic traits (e.g. body size, ornaments or armaments) may imply unequal physiological costs (e.g. developmental, maintenance) between the sexes and this may result in differences in longevity and aging. Second, mitochondria are usually transmitted by the mother and thus selection is blind to mitochondrial deleterious mutations affecting only males. Such mutations can freely accumulate in the mitochondrial genome and may reduce male longevity, a phenomenon called the mother’s curse (Frank and Hurst 1996). Third, in species with sex chromosomes, all recessive deleterious mutations will be expressed on the single X chromosome in XY males and may reduce their longevity (the unguarded X effect). In addition, the numerous transposable elements (TEs) on the Y chromosome may affect aging. TE activity is normally repressed by epigenetic regulation (DNA methylation, histone modifications and small RNAs). However, it is known that this regulation is disrupted with increasing age. Because of the TE-rich Y chromosome, more TEs may become active in old males than in old females, generating more somatic mutations, accelerating aging and reducing longevity in males (the toxic Y effect, Marais et al. 2018).
The relative contributions of these different effects to SGLs remain unknown. Sex-differences in life history strategies have been considered as the most important cause of SGLs for long (Tidière et al. 2015) but this effect remain equivocal (Lemaître et al. 2020) and cannot explain alone the diversity of patterns observed across species (Marais et al. 2018). Similarly, while studies in Drosophila and humans have shown that the mother’s curse contributes to SGLs in those organisms (e.g. Milot et al. 2007), its contribution may not be strong. Recently, two large-scale comparative analyses have shown that in species with XY chromosomes males show a shorter lifespan compared to females, while in species with ZW chromosomes (a system in which the female are the heterogametic sex and are ZW, and the males ZZ) the opposite pattern is observed (Pipoly et al. 2015; Xirocostas et al. 2020). Apart from these correlational studies, very little experimental tests of the effect of sex chromosomes on longevity have been conducted. In Drosophila, the evidence suggests that the unguarded X effect does not contribute to SGLs (Brengdahl et al. 2018). Whether a toxic Y effect exists in this species was unknown.
In a very elegant study, Brown et al. (2020) provided strong evidence for such a toxic Y effect in Drosophila melanogaster. First, they checked that in the D. melanogaster strain that they were studying (Canton-S), males were indeed dying younger than females. They also confirmed that in this strain, as in others, the male genomes include more repeats and heterochromatin than the female ones using cytometry. A careful analysis of the heterochromatin (using H3K9me2, a repressive histone modification typical of heterochromatin, as a proxy) in old flies revealed that heterochromatin loss was much more important in males than in females, in particular on the Y chromosome (but also to a lesser extent at the pericentromric regions of the autosomes). This change in heterochromatin had two outputs, they found. First, the expression of the genes in those regions was affected. They highlighted that many of such genes are involved in immunity and regulation with a potential impact on longevity. Second, they found a striking TE reactivation. These two effects were stronger in males. While females showed clear reactivation of 6 TEs, with the total fraction of repeats in the transcriptome going from 2% (young females) to 4.6% (old females), males experienced the reactivation of 32 TEs, with the total fraction of repeats in the transcriptome going from 1.6% (young males) to 5.8% (old males). It appeared that most of these TEs are Y-linked. And when focusing on Y-linked repeats, they found that 32 Y-linked TEs became upregulated during male aging and the fraction of Y-linked TEs in the transcriptome increased ninefold.
All these observations clearly suggested that male longevity was decreased because of a toxic Y effect. To really uncover a causal relationship between having a Y chromosome and shorter longevity, Brown et al. (2020) artificially produced flies with atypical karyotypes: X0 males, XXY females and XYY males. This is very interesting as they could uncouple the effect of the phenotypical sex (being male or female) and having a Y chromosome or not, as in fruit flies sex is determined not by the Y chromosome but by the X/autosome ratio. Their results are striking. They found that longevity of the X0 males was the highest (higher than XX females in fact), and that of the XYY males the lowest. Females XXY had intermediate longevities. Importantly, this was found to be robust to genomic background as results were the same using crosses from different strains. When analysing TEs of these flies, they found a particularly strong expression of the Y-linked TEs in old XXY and XYY flies. Interestingly, in young XXY and XYY flies Y-linked TEs expression was also strong, suggesting the chromatin regulation of the Y chromosome is disrupted in these flies.
This work points to the idea that SGLs in D. melanogaster are mainly explained by the toxic Y effect. The molecular details however remain to be elucidated. The effect of the Y chromosome on aging might be more complex than envisioned in the toxic Y model presented above. Brown et al. (2020) indeed found that heterochromatin loss was globally faster in males, both at the Y chromosome and the autosomes. The organisation of the nucleus, in particular of the nucleolus, which is involved in heterochromatin maintenance, involves the sex chromosomes in D. melanogaster as discussed in the paper, and may explain this observation. The epigenetic status of the Y chromosome is known to affect that of all the autosomes in Drosophila (Lemos et al. 2008). Also, in Brown et al. (2020) most of the work (in particular the genomic part) has been done on Canton-S. Only D. melanogaster was studied but limited data suggest different Drosophila species may have different SGLs. The TE analysis is known to be tricky, different tools to analyse TE expression exist (e.g. Lerat et al. 2017; Lanciano and Cristofari 2020). Future work should focus on testing the toxic Y effect on other D. melanogaster strains and other Drosophila species, using different tools to study TE expression, and on dissecting the molecular details of the toxic Y effect.

References

Austad, S. N., and Bartke, A. (2015). Sex differences in longevity and in responses to anti-aging interventions: A Mini-Review. Gerontology, 62(1), 40–46. 10.1159/000381472
Brengdahl, M., Kimber, C. M., Maguire-Baxter, J., and Friberg, U. (2018). Sex differences in life span: Females homozygous for the X chromosome do not suffer the shorter life span predicted by the unguarded X hypothesis. Evolution; international journal of organic evolution, 72(3), 568–577. 10.1111/evo.13434
Brown, E. J., Nguyen, A. H., and Bachtrog, D. (2020). The Y chromosome may contribute to sex-specific ageing in Drosophila. Nature ecology and evolution, 4(6), 853–862. 10.1038/s41559-020-1179-5 or preprint link on bioRxiv
Frank, S. A., and Hurst, L. D. (1996). Mitochondria and male disease. Nature, 383(6597), 224. 10.1038/383224a0
Lanciano, S., and Cristofari, G. (2020). Measuring and interpreting transposable element expression. Nature reviews. Genetics, 10.1038/s41576-020-0251-y. Advance online publication. 10.1038/s41576-020-0251-y
Lemaître, J. F., Ronget, V., Tidière, M., Allainé, D., Berger, V., Cohas, A., Colchero, F., Conde, D. A., Garratt, M., Liker, A., Marais, G., Scheuerlein, A., Székely, T., and Gaillard, J. M. (2020). Sex differences in adult lifespan and aging rates of mortality across wild mammals. Proceedings of the National Academy of Sciences of the United States of America, 117(15), 8546–8553. 10.1073/pnas.1911999117
Lemos, B., Araripe, L. O., and Hartl, D. L. (2008). Polymorphic Y chromosomes harbor cryptic variation with manifold functional consequences. Science (New York, N.Y.), 319(5859), 91–93. 10.1126/science.1148861
Lerat, E., Fablet, M., Modolo, L., Lopez-Maestre, H., and Vieira, C. (2017). TEtools facilitates big data expression analysis of transposable elements and reveals an antagonism between their activity and that of piRNA genes. Nucleic acids research, 45(4), e17. 10.1093/nar/gkw953
Marais, G., Gaillard, J. M., Vieira, C., Plotton, I., Sanlaville, D., Gueyffier, F., and Lemaitre, J. F. (2018). Sex gap in aging and longevity: can sex chromosomes play a role?. Biology of sex differences, 9(1), 33. 10.1186/s13293-018-0181-y
Milot, E., Moreau, C., Gagnon, A., Cohen, A. A., Brais, B., and Labuda, D. (2017). Mother's curse neutralizes natural selection against a human genetic disease over three centuries. Nature ecology and evolution, 1(9), 1400–1406. 10.1038/s41559-017-0276-6
Pipoly, I., Bókony, V., Kirkpatrick, M., Donald, P. F., Székely, T., and Liker, A. (2015). The genetic sex-determination system predicts adult sex ratios in tetrapods. Nature, 527(7576), 91–94. 10.1038/nature15380
Tidière, M., Gaillard, J. M., Müller, D. W., Lackey, L. B., Gimenez, O., Clauss, M., and Lemaître, J. F. (2015). Does sexual selection shape sex differences in longevity and senescence patterns across vertebrates? A review and new insights from captive ruminants. Evolution; international journal of organic evolution, 69(12), 3123–3140. 10.1111/evo.12801
Xirocostas, Z. A., Everingham, S. E., and Moles, A. T. (2020). The sex with the reduced sex chromosome dies earlier: a comparison across the tree of life. Biology letters, 16(3), 20190867. 10.1098/rsbl.2019.0867

The advent of genome-wide association studies (GWAS) has been a great promise for our understanding of the connection between genotype and phenotype. Today, the NHGRI-EBI GWAS catalog contains 251,401 associations from 4,961 studies (1). This wealth of studies has also generated interest to use the summary statistics beyond the few top hits in order to make predictions for individuals without known phenotype, e.g. to predict polygenic risk scores or to study polygenic selection by comparing different groups. For instance, polygenic selection acting on the most studied polygenic trait, height, has been subject to multiple studies during the past decade (e.g. 2–6). They detected north-south gradients in Europe which were consistent with expectations. However, their GWAS summary statistics were based on the GIANT consortium data set, a meta-analysis of GWAS conducted in different European cohorts (7,8). The availability of large data sets with less stratification such as the UK Biobank (9) has led to a re-evaluation of those results. The nature of the GIANT consortium data set was realized to represent a potential problem for studies of polygenic adaptation which led several of the authors of the original articles to caution against the interpretations of polygenic selection on height (10,11). This was a great example on how the scientific community assessed their own earlier results in a critical way as more data became available. At the same time it left the question whether there is detectable polygenic selection separating populations more open than ever.

Generally, recent years have seen several articles critically assessing the portability of GWAS results and risk score predictions to other populations (12–14). Refoyo-Martínez et al. (15) are now presenting a systematic assessment on the robustness of cross-population signatures of polygenic adaptation in humans. They compiled GWAS results for complex traits which have been studied in more than one cohort and then use allele frequencies from the 1000 Genomes Project data (16) set to detect signals of polygenic score overdispersion. As the source for the allele frequencies is kept the same across all tests, differences between the signals must be caused by the underlying GWAS. The results are concerning as the level of overdispersion largely depends on the choice of GWAS cohort. Cohorts with homogenous ancestries show little to no overdispersion compared to cohorts of mixed ancestries such as meta-analyses. It appears that the meta-analyses fail to fully account for stratification in their data sets.

The authors based most of their analyses on the heavily studied trait height. Additionally, they use educational attainment (measured as the number of school years of an individual) as an example. This choice was due to the potential over- or misinterpretation of results by the media, the general public and by far right hate groups. Such traits are potentially confounded by unaccounted cultural and socio-economic factors. Showing that previous results about polygenic selection on educational attainment are not robust is an important result that needs to be communicated well. This forms a great example for everyone working in human genomics. We need to be aware that our results can sometimes be misinterpreted. And we need to make an effort to write our papers and communicate our results in a way that is honest about the limitations of our research and that prevents the misuse of our results by hate groups.

This article represents an important contribution to the field. It is cruicial to be aware of potential methodological biases and technical artifacts. Future studies of polygenic adaptation need to be cautious with their interpretations of polygenic score overdispersion. A recommendation would be to use GWAS results obtained in homogenous cohorts. But even if different biobank-scale cohorts of homogeneous ancestry are employed, there will always be some remaining risk of unaccounted stratification. These conclusions may seem sobering but they are part of the scientific process. We need additional controls and new, different methods than polygenic score overdispersion for assessing polygenic selection. Last year also saw the presentation of a novel approach using sequence data and GWAS summary statistics to detect directional selection on a polygenic trait (17). This new method appears to be robust to bias stemming from stratification in the GWAS cohort as well as other confounding factors. Such new developments show light at the end of the tunnel for the use of GWAS summary statistics in the study of polygenic adaptation.

References

1. Buniello A, MacArthur JAL, Cerezo M, Harris LW, Hayhurst J, Malangone C, et al. The NHGRI-EBI GWAS Catalog of published genome-wide association studies, targeted arrays and summary statistics 2019. Nucleic Acids Research. 2019 Jan 8;47(D1):D1005–12. doi: https://doi.org/10.1093/nar/gky1120

2. Turchin MC, Chiang CW, Palmer CD, Sankararaman S, Reich D, Hirschhorn JN. Evidence of widespread selection on standing variation in Europe at height-associated SNPs. Nature Genetics. 2012 Sep;44(9):1015–9. doi: https://doi.org/10.1038/ng.2368

3. Berg JJ, Coop G. A Population Genetic Signal of Polygenic Adaptation. PLOS Genetics. 2014 Aug 7;10(8):e1004412. doi: https://doi.org/10.1371/journal.pgen.1004412

4. Robinson MR, Hemani G, Medina-Gomez C, Mezzavilla M, Esko T, Shakhbazov K, et al. Population genetic differentiation of height and body mass index across Europe. Nature Genetics. 2015 Nov;47(11):1357–62. doi: https://doi.org/10.1038/ng.3401

5. Mathieson I, Lazaridis I, Rohland N, Mallick S, Patterson N, Roodenberg SA, et al. Genome-wide patterns of selection in 230 ancient Eurasians. Nature. 2015 Dec;528(7583):499–503. doi: https://doi.org/10.1038/nature16152

6. Racimo F, Berg JJ, Pickrell JK. Detecting polygenic adaptation in admixture graphs. Genetics. 2018. Arp;208(4):1565–1584. doi: https://doi.org/10.1534/genetics.117.300489

7. Lango Allen H, Estrada K, Lettre G, Berndt SI, Weedon MN, Rivadeneira F, et al. Hundreds of variants clustered in genomic loci and biological pathways affect human height. Nature. 2010 Oct;467(7317):832–8. doi: https://doi.org/10.1038/nature09410

8. Wood AR, Esko T, Yang J, Vedantam S, Pers TH, Gustafsson S, et al. Defining the role of common variation in the genomic and biological architecture of adult human height. Nat Genet. 2014 Nov;46(11):1173–86. doi: https://doi.org/10.1038/ng.3097

9. Bycroft C, Freeman C, Petkova D, Band G, Elliott LT, Sharp K, et al. The UK Biobank resource with deep phenotyping and genomic data. Nature. 2018 Oct;562(7726):203–9. doi: https://doi.org/10.1038/s41586-018-0579-z

10. Berg JJ, Harpak A, Sinnott-Armstrong N, Joergensen AM, Mostafavi H, Field Y, et al. Reduced signal for polygenic adaptation of height in UK Biobank. eLife. 2019 Mar 21;8:e39725. doi: https://doi.org/10.7554/eLife.39725

11. Sohail M, Maier RM, Ganna A, Bloemendal A, Martin AR, Turchin MC, et al. Polygenic adaptation on height is overestimated due to uncorrected stratification in genome-wide association studies. eLife. 2019 Mar 21;8:e39702. doi: https://doi.org/10.7554/eLife.39702

12. Martin AR, Kanai M, Kamatani Y, Okada Y, Neale BM, Daly MJ. Clinical use of current polygenic risk scores may exacerbate health disparities. Nature Genetics. 2019 Apr;51(4):584–91. doi: https://doi.org/10.1038/s41588-019-0379-x

13. Bitarello BD, Mathieson I. Polygenic Scores for Height in Admixed Populations. G3: Genes, Genomes, Genetics. 2020 Nov 1;10(11):4027–36. doi: https://doi.org/10.1534/g3.120.401658

14. Uricchio LH, Kitano HC, Gusev A, Zaitlen NA. An evolutionary compass for detecting signals of polygenic selection and mutational bias. Evolution Letters. 2019;3(1):69–79. doi: https://doi.org/10.1002/evl3.97

15. Refoyo-Martínez A, Liu S, Jørgensen AM, Jin X, Albrechtsen A, Martin AR, Racimo F. How robust are cross-population signatures of polygenic adaptation in humans? bioRxiv, 2021, 2020.07.13.200030, version 5 peer-reviewed and recommended by Peer community in Evolutionary Biology. doi: https://doi.org/10.1101/2020.07.13.200030

16. Auton A, Abecasis GR, Altshuler DM, Durbin RM, Abecasis GR, Bentley DR, et al. A global reference for human genetic variation. Nature. 2015 Sep 30;526(7571):68–74. doi: https://doi.org/10.1038/nature15393

17. Stern AJ, Speidel L, Zaitlen NA, Nielsen R. Disentangling selection on genetically correlated polygenic traits using whole-genome genealogies. bioRxiv. 2020 May 8;2020.05.07.083402. doi: https://doi.org/10.1101/2020.05.07.083402

The social amoeba Dictyostelium discoideum is an important model system for the study of cooperation and multicellularity as is has both unicellular and aggregative life phases. In the aggregative phase, which typically occurs when nutrients are limiting, individual cells eventually gather together to form a fruiting bodies whose spores may be dispersed to another, better, location and whose stalk cells, which support the spores, die. This extreme form of cooperation has been the focus of numerous studies that have revealed the importance genetic relatedness and kin selection (Hamilton 1964; Lehmann and Rousset 2014) in explaining the maintenance of this cooperative collective behavior (Strassmann et al. 2000; Kuzdzal-Fick et al. 2011; Strassmann and Queller 2011). However, much remains unknown with respect to how the interactions between individual cells, their neighbors, and their environment produce cooperative behavior at the scale of whole groups or collectives. In this preprint, Forget et al. (2021) describe how the D. discoideum system is crucial in this respect because it allows these cellular-level interactions to be studied in a systematic and tractable manner.
Spore bias, which is the tendency of a particular genotype or strain to disproportionately migrate to the spore instead of the stalk, is often used to define which strains are "cheaters" (positive spore bias) and which are "cooperative" (negative spore bias). Forget et al. (2021) note that spore bias depends on a number of stochastic factors including external drivers such as variation in environmental (or nutrient) quality and internal drivers like cell-cycle phase at the time of starvation. Spore bias is also affected by the social environment where the fraction of cheater strains in a spore may be limited by the ability of the remaining stalk cells to support the spore. The social environment can also affect cells through their differential responsiveness to the chemical factors that induce differentiation into stalk cells; responsiveness is partly a function of nutrient quality (Thompson and Kay 2000), which in turn can be a function of cell density. Thus, Forget et al. (2021) highlight a number of mechanisms that could generate frequency-dependent selection that would lead to the stable maintenance of multiple strains with different spore biases; in other words, both cheater and cooperative strains might stably coexist due to these cellular-level interactions.
The cellular-level interactions that Forget et al. (2021) highlight are particularly important because they pose a challenge evolutionary theory: some evolutionary models of social and collective behavior neglect or simplify these interactions. For example, Forget et al. (2021) note that the developmental, behavior, and environmental timescales relevant for Dictyostelium fruiting body formation all overlap. Evolutionary analyses often assume some of these timescales, for example developmental and behavior, are separate in order to simplify the analysis of any interactions. Thus, new theoretical work that allows these timescales to overlap may shed light on how cellular-level interactions can produce environmental, physiological, and behavioral feedbacks that drive the evolution of cooperation and other collective behaviors.

References

Forget, M., Adiba, S. and De Monte, S.(2021) Social conflicts in *Dictyostelium discoideum *: a matter of scales. HAL, hal-03088868, ver. 2 peer-reviewed and recommended by PCI Evolutionary Biology. https://hal.archives-ouvertes.fr/hal-03088868/

Hamilton, W. D. (1964). The genetical evolution of social behaviour. II. Journal of theoretical biology, 7(1), 17-52. doi: https://doi.org/10.1016/0022-5193(64)90039-6

Kuzdzal-Fick, J. J., Fox, S. A., Strassmann, J. E., and Queller, D. C. (2011). High relatedness is necessary and sufficient to maintain multicellularity in Dictyostelium. Science, 334(6062), 1548-1551. doi: https://doi.org/10.1126/science.1213272

Lehmann, L., and Rousset, F. (2014). The genetical theory of social behaviour. Philosophical Transactions of the Royal Society B: Biological Sciences, 369(1642), 20130357. doi: https://doi.org/10.1098/rstb.2013.0357

Strassmann, J. E., and Queller, D. C. (2011). Evolution of cooperation and control of cheating in a social microbe. Proceedings of the National Academy of Sciences, 108(Supplement 2), 10855-10862. doi: https://doi.org/10.1073/pnas.1102451108

Strassmann, J. E., Zhu, Y., & Queller, D. C. (2000). Altruism and social cheating in the social amoeba Dictyostelium discoideum. Nature, 408(6815), 965-967. doi: https://doi.org/10.1038/35050087

Thompson, C. R., & Kay, R. R. (2000). Cell-fate choice in Dictyostelium: intrinsic biases modulate sensitivity to DIF signaling. Developmental biology, 227(1), 56-64. doi: https://doi.org/10.1006/dbio.2000.9877

The role of microbial symbionts in governing social traits of their hosts is an exciting and developing research area. Just like symbionts influence host reproductive behaviour and can cause mating incompatibilities to promote symbiont transmission through host populations (Engelstadter and Hurst 2009; Correa and Ballard 2016; Johnson and Foster 2018) (see also discussion on conflict resolution in Johnsen and Foster 2018), microbial symbionts could enhance transmission by promoting the social behaviour of their hosts (Ezenwa et al. 2012; Lewin-Epstein et al. 2017; Gurevich et al. 2020). Here I apply the term ‘symbiosis’ in the broad sense, following De Bary 1879 as “the living together of two differently named organisms“ independent of effects on the organisms involved (De Bary 1879), i.e. the biological interaction between the host and its symbionts may include mutualism, parasitism and commensalism.
So far, we have relative few studies that explore the role of symbionts in promoting social behaviours such as parental care. Clearly, disentangling cause and effect when assessing the functional significance of symbiotic relationships in general is extremely challenging, and perhaps even more caution is needed when assessing the role of symbionts in the evolution of parental care, given the high fitness benefits to the offspring of receiving care. An interesting study on the symbiotic relationship between termites and their eukaryotic gut symbionts proposes a role of gut flagellates in the origin of subsocial behaviour (extended offspring care) in the termites through proctodeal trophallaxis (i.e. anus-to-mouth feeding), driven by mutualistic beneficial interactions (Nalepa 2020). Van Meyel et al. (2021) hypothesized a role of gut symbionts in promoting maternal care behaviour in the European earwig, and set out to test this idea in a carefully executed experimental study. They used a broad-spectrum antibiotic treatment to alter gut microbiota in mothers and examined its effect on maternal care provisioning. While the antibiotic treatment altered the gut microbiome, no effect on pre- or post-hatching maternal care was detected. The authors also investigated a broad range of physiological and reproductive traits measured over a major part of a female’s lifetime, and detected no effect of microbiome alteration on these traits. The study therefore does not provide evidence for a direct role of the gut microbiome in shaping offspring care in this population of European earwigs.
Within populations, earwigs show inter-individual variation in the expression of maternal care (Meunier et al. 2012; Ratz et al. 2016), and there is evidence that genetic and environmental factors contribute to this this variation (Meunier and Kolliker 2012; Kramer et al. 2017). The study by Van Meyel et al. (2021) is the first to analyse microbiome composition of the European earwig, and they study host-symbiont associations in a single population. A next step could be to explore among population variation in the gut microbiome, to achieve a better understanding on host-microbiome variation and dynamics in wild populations. Depending on the nature of host-symbiont associations across populations, new perspectives on their functional significance may arise (Hird 2017; Johnson and Foster 2018). It is therefore too early to conclusively confirm or reject the role of microbial symbionts in the expression of parental care in this system.

References

Correa, C. C., and Ballard, J. W. O. (2016). Wolbachia associations with insects: winning or losing against a master manipulator. Frontiers in Ecology and Evolution, 3, 153. doi: https://doi.org/10.3389/fevo.2015.00153

De Bary, A. (1879). Die Erscheinung der Symbiose. Verlag von Karl J. Trubner, Strassburg.

Engelstädter, J., and Hurst, G. D. (2009). The ecology and evolution of microbes that manipulate host reproduction. Annual Review of Ecology, Evolution, and Systematics, 40, 127-149. doi: https://doi.org/10.1146/annurev.ecolsys.110308.120206

Ezenwa, V. O., Gerardo, N. M., Inouye, D. W., Medina, M., and Xavier, J. B. (2012). Animal behavior and the microbiome. Science, 338(6104), 198-199. doi: https://doi.org/10.1126/science.1227412

Gurevich, Y., Lewin-Epstein, O., and Hadany, L. (2020). The evolution of paternal care: a role for microbes?. Philosophical Transactions of the Royal Society B, 375(1808), 20190599. doi: https://doi.org/10.1098/rstb.2019.0599

Hird, S. M. (2017). Evolutionary biology needs wild microbiomes. Frontiers in microbiology, 8, 725. doi: https://doi.org/10.3389/fmicb.2017.00725

Johnson, K. V. A., and Foster, K. R. (2018). Why does the microbiome affect behaviour?. Nature reviews microbiology, 16(10), 647-655. doi: https://doi.org/10.1038/s41579-018-0014-3

Kramer et al. (2017). When earwig mothers do not care to share: parent–offspring competition and the evolution of family life. Functional Ecology, 31(11), 2098-2107. doi: https://doi.org/10.1111/1365-2435.12915

Lewin-Epstein, O., Aharonov, R., and Hadany, L. (2017). Microbes can help explain the evolution of host altruism. Nature communications, 8(1), 1-7. doi: https://doi.org/10.1038/ncomms14040

Meunier, J., and Kölliker, M. (2012). Parental antagonism and parent–offspring co-adaptation interact to shape family life. Proceedings of the Royal Society B: Biological Sciences, 279(1744), 3981-3988. doi: https://doi.org/10.1098/rspb.2012.1416

Meunier, J., Wong, J. W., Gómez, Y., Kuttler, S., Röllin, L., Stucki, D., and Kölliker, M. (2012). One clutch or two clutches? Fitness correlates of coexisting alternative female life-histories in the European earwig. Evolutionary Ecology, 26(3), 669-682. doi: https://doi.org/10.1007/s10682-011-9510-x

Nalepa, C. A. (2020). Origin of mutualism between termites and flagellated gut protists: transition from horizontal to vertical transmission. Frontiers in Ecology and Evolution, 8, 14. doi: https://doi.org/10.3389/fevo.2020.00014

Ratz, T., Kramer, J., Veuille, M., and Meunier, J. (2016). The population determines whether and how life-history traits vary between reproductive events in an insect with maternal care. Oecologia, 182(2), 443-452. doi: https://doi.org/10.1007/s00442-016-3685-3

Van Meyel, S., Devers, S., Dupont, S., Dedeine, F. and Meunier, J. (2021) Alteration of gut microbiota with a broad-spectrum antibiotic does not impair maternal care in the European earwig. bioRxiv, 2020.10.08.331363. ver. 5 peer-reviewed and recommended by PCI Evol Biol. https://doi.org/10.1101/2020.10.08.331363

Historically neglected or considered a rare aberration in animals under the biological species concept, interspecific hybridisation has by now been recognised to be taxonomically widespread, particularly in rapidly diversifying groups (Dagilis et al. 2021; Edelman & Mallet 2021; Mallet et al. 2016; Seehausen 2004). Yet the prevalence of introgressive hybridizations, its evolutionary significance, and its impact on species diversification remain a hot topic of research in evolutionary biology. The rapid increase in genomic resources now available for non-model species has significantly contributed to the detection of introgressive hybridization across taxa showing that reticulated evolution is far more common in the animal kingdom than historically considered. Yet, detecting it, quantifying its magnitude, and assessing its evolutionary significance remains a challenging endeavour with constantly evolving methodologies to better explore and exploit genomic data (Blair & Ané 2020; Degnan & Rosenberg 2009; Edelman & Mallet 2021; Hibbins & Hahn 2022).

In the marine realm, the dearth of geographic barriers and the large dispersal abilities of pelagic species like cetaceans have raised the questions of how populations and species can diverge and adapt to distinct ecological conditions in face of potentially large gene-flow, the so-called marine speciation paradox (Bierne et al. 2003). Contemporaneous hybridization among cetacean species has been widely documented in nature despite large phenotypic differences (Crossman et al. 2016). The historical prevalence of reticulated evolution, its evolutionary significance, and how it might have impacted the evolutionary history and diversification of the cetaceans have however remained elusive so far. Recent phylogenomic studies suggested that introgression has been prevalent in cetacean evolutionary history with instances reported among baleen whales (mysticetes) (Árnason et al. 2018) and among toothed whales (odontocetes), especially in the rapidly diversifying dolphins family of the Delphininae (Guo et al. 2021; Moura et al. 2020).

Analysing publicly available whole-genome data from nine cetacean species across three families of Delphinoidae – dolphins, porpoises, and monondontidae – using phylogenomics and demo-genetics approaches, Westbury and colleagues (2022) take a step further in documenting that evolution among these species has been far from a simple bifurcating tree. Instead, their study describes widespread occurrences of introgression among Delphinoidae, drawing a complex picture of reticulated evolutionary history. After describing major topology discordance in phylogenetic gene trees along the genome, the authors use a panel of approaches to disentangle introgression from incomplete lineage sorting (ILS), the two most common causes of tree topology discordances (Hibbins & Hahn 2022). Applying popular tests that separate introgression from ILS, such as the Patterson’s D (a.k.a. ABBA-BABA) test (Durand et al. 2011; Green et al. 2010), QuIBL (Edelman et al. 2019), and D-FOIL (Pease & Hahn 2015), the authors report that signals of introgression are present in the genomes of most (if not all) the cetacean species included in their study. However, this picture needs to be nuanced. Most introgression signals seem to echo old introgression events that occurred primarily among ancestors. This is suggested by the differential signals of topology discordance along the genome when considering sliding windows along the genome of varying sizes (50kb, 100kb, and 1Mb), and by patterns of excess derived allele sharing along branches of the species tree, as captured by the f-branch test (Malinsky et al. 2021; Malinsky et al. 2018). The authors further investigated the dynamic of cessation of gene flow (and/or ILS) between species using the F1 hybrid PSMC (or hPSMC) approach (Cahill et al. 2016). By estimating the cross-coalescent rates (CRR) between species pairs with time in light of previously estimated species divergence times (McGowen et al. 2020), the authors report that gene flow (and/or ILS) most likely has stopped by now among most species, but it may have lasted for more than half of the time since the species split from each other. According to the author, this result may reflect the slow process by which reproductive isolation would have evolved between cetacean lineages, and that species isolation was marked by significant introgression events.

Now, while the present study provides a good overview of how complex is the reticulated evolutionary history of the Delphinoidae, getting a complete picture will require overcoming a few important limitations. The first ones are methodological and related to the phylogenomic analyses. Given the sampling design with one diploid genome per species, the authors could not phase the data into the parental haplotypes, but instead relied on a majority consensus creating mosaic haploidized genomes representing a mixture between the two parental copies. Moreover, by using large genomic windows (≥50kb) that likely span multiple independent loci, phylogenetic analyses in windows encompassed distinct phylogenetic signals, potentially leading to bias and inaccuracy in the inferences. Thawornwattana et al (2018) previously showed that this “concatenation approach”  could significantly impact phylogenetic inferences. They proposed instead to use loci small enough to minimise the risk of intra-locus recombination and to consider them in blocks of non-recombining loci along the genome in order to conduct phylogenetic analysed, ideally under the multi-species coalescent (MSC) that can account for ILS (e.g. BPP; Flouri et al. 2018; Jiao et al. 2020; Yang 2015). Such an approach applied to the diversification of the Delphinidae may reveal substantial changes compared to the currently admitted species tree.

Inaccuracy in the species tree estimation may have a major impact on the introgression analyses conducted in this study since the species tree and branching order must be supplied in the introgression analyses to properly disentangle introgression from ILS. Here, the authors rely on the tree topology that was previously estimated in McGowen et al. (2020), which they also recovered using their consensus estimation from ASTRAL-III (Zhang et al. 2018). While the methodologies accounted to a certain extent for ILS, albeit with potential bias induced by the concatenation approach, they ignore the presumably large amount of introgression among species during the diversification process. Estimating species branching order while ignoring introgression can lead to major bias in the phylogenetic inference and can lead to incorrect topologies. Even if these MSC-based methods account for ILS, inferences can become very inaccurate or even break down as gene flow increases (see for ex. Jiao et al. 2020; Müller et al. in press; Solís-Lemus et al. 2016). Dedicated approaches have been developed to model explicitly introgression together with ILS to estimate phylogenetic networks (Blair & Ané 2020; Rabier et al. 2021) or in isolation-with-migration model (Müller et al. in press; Wang et al. 2020). Future studies revisiting the reticulated evolutionary history of the Delphinoidae with such approaches may not only precise the species branching order, but also deliver a finer view of the magnitude and prevalence of introgression during the evolutionary history of these species.

A final part of Westbury et al. (2022)'s study set out to test whether historical periods of low abundance could have facilitated hybridization among Delphinoidae species. During these periods of low abundance, species may encounter only a limited number of conspecifics and may consider individuals from other species as suitable mating partners, leading to hybridisation (Crossman et al. 2016; Edwards et al. 2011; Westbury et al. 2019). The authors tested this hypothesis by considering genome-wide genetic diversity (or heterozygosity) as a proxy of historical effective population size (Ne), itself as a proxy of the evolution of census size with time. They also try to link historical Ne variation (from PSMC, Li & Durbin 2011) with their estimated time to cessation of gene flow or ILS (from the CRR of hPSMC). However, no straightforward relationship was found between the genetic diversity and the propensity of species to hybridize, nor was there any clear link between Ne variation through time and the cessation of gene flow or ILS. Such a lack of relationship may not come as a surprise, since the determinants of genome-wide genetic diversity and its variation through evolutionary time-scale are far more diverse and complex than just a direct link with hybridization, introgression, or even with the census population size. In fact, genetic diversity results from the balance between all the evolutionary processes at play in the species' evolutionary history (see the review of Ellegren & Galtier 2016). Other important factors can strongly impact genetic diversity, including demography and structure, but also linked selection (Boitard et al. 2022; Buffalo 2021; Ellegren & Galtier 2016). 

All in all, Westbury and coll. (2022) present here an interesting study providing an additional step towards resolving and understanding the complex evolutionary history of the Delphinoidae, and shedding light on the importance of introgression during the diversification of these cetacean species. Prospective work improving upon the taxonomic sampling, with additional genomic data for each species, considered with dedicated approaches tailored at estimating species tree or network while accounting for ILS and introgression will be key for refining the picture depicted in this study. Down the road, altogether these studies will contribute to assessing the evolutionary significance of introgression on the diversification of Delphinoides, and more generally on the diversification of cetacean species, which still remain an open and exciting perspective. 

References

Árnason Ú, Lammers F, Kumar V, Nilsson MA, Janke A (2018) Whole-genome sequencing of the blue whale and other rorquals finds signatures for introgressive gene flow. Science Advances 4, eaap9873. https://doi.org/10.1126/sciadv.aap9873

Bierne N, Bonhomme F, David P (2003) Habitat preference and the marine-speciation paradox. Proceedings of the Royal Society of London. Series B: Biological Sciences 270, 1399-1406. https://doi.org/10.1098/rspb.2003.2404

Blair C, Ané C (2020) Phylogenetic Trees and Networks Can Serve as Powerful and Complementary Approaches for Analysis of Genomic Data. Systematic Biology 69, 593-601. https://doi.org/10.1093/sysbio/syz056

Boitard S, Arredondo A, Chikhi L, Mazet O (2022) Heterogeneity in effective size across the genome: effects on the inverse instantaneous coalescence rate (IICR) and implications for demographic inference under linked selection. Genetics 220, iyac008. https://doi.org/10.1093/genetics/iyac008

Buffalo V (2021) Quantifying the relationship between genetic diversity and population size suggests natural selection cannot explain Lewontin's Paradox. e-Life 10, e67509. https://doi.org/10.7554/eLife.67509

Cahill JA, Soares AE, Green RE, Shapiro B (2016) Inferring species divergence times using pairwise sequential Markovian coalescent modelling and low-coverage genomic data. Philos Trans R Soc Lond B Biol Sci 371, 20150138. https://doi.org/10.1098/rstb.2015.0138

Crossman CA, Taylor EB, Barrett‐Lennard LG (2016) Hybridization in the Cetacea: widespread occurrence and associated morphological, behavioral, and ecological factors. Ecology and Evolution 6, 1293-1303. https://doi.org/10.1002/ece3.1913

Dagilis AJ, Peede D, Coughlan JM, Jofre GI, D’Agostino ERR, Mavengere H, Tate AD, Matute DR (2021) 15 years of introgression studies: quantifying gene flow across Eukaryotes. biorXiv, 2021.1106.1115.448399. https://doi.org/10.1101/2021.06.15.448399

Degnan JH, Rosenberg NA (2009) Gene tree discordance, phylogenetic inference and the multispecies coalescent. Trends Ecol Evol 24, 332-340. https://doi.org/10.1016/j.tree.2009.01.009

Durand EY, Patterson N, Reich D, Slatkin M (2011) Testing for ancient admixture between closely related populations. Mol Biol Evol 28, 2239-2252. https://doi.org/10.1093/molbev/msr048

Edelman NB, Frandsen PB, Miyagi M, Clavijo B, Davey J, Dikow RB, Garcia-Accinelli G, Van Belleghem SM, Patterson N, Neafsey DE, Challis R, Kumar S, Moreira GRP, Salazar C, Chouteau M, Counterman BA, Papa R, Blaxter M, Reed RD, Dasmahapatra KK, Kronforst M, Joron M, Jiggins CD, McMillan WO, Di Palma F, Blumberg AJ, Wakeley J, Jaffe D, Mallet J (2019) Genomic architecture and introgression shape a butterfly radiation. Science 366, 594-599. https://doi.org/10.1126/science.aaw2090

Edelman NB, Mallet J (2021) Prevalence and Adaptive Impact of Introgression. Annual Review of Genetics 55, 265-283. https://doi.org/10.1146/annurev-genet-021821-020805

Edwards CJ, Suchard MA, Lemey P, Welch JJ, Barnes I, Fulton TL, Barnett R, O'Connell TC, Coxon P, Monaghan N, Valdiosera CE, Lorenzen ED, Willerslev E, Baryshnikov GF, Rambaut A, Thomas MG, Bradley DG, Shapiro B (2011) Ancient hybridization and an Irish origin for the modern polar bear matriline. Curr Biol 21, 1251-1258. https://doi.org/10.1016/j.cub.2011.05.058

Ellegren H, Galtier N (2016) Determinants of genetic diversity. Nat Rev Genet 17, 422-433. https://doi.org/10.1038/nrg.2016.58

Flouri T, Jiao X, Rannala B, Yang Z (2018) Species Tree Inference with BPP Using Genomic Sequences and the Multispecies Coalescent. Mol Biol Evol 35, 2585-2593. https://doi.org/10.1093/molbev/msy147

Green RE, Krause J, Briggs AW, Maricic T, Stenzel U, Kircher M, Patterson N, Li H, Zhai W, Fritz MH, Hansen NF, Durand EY, Malaspinas AS, Jensen JD, Marques-Bonet T, Alkan C, Prufer K, Meyer M, Burbano HA, Good JM, Schultz R, Aximu-Petri A, Butthof A, Hober B, Hoffner B, Siegemund M, Weihmann A, Nusbaum C, Lander ES, Russ C, Novod N, Affourtit J, Egholm M, Verna C, Rudan P, Brajkovic D, Kucan Z, Gusic I, Doronichev VB, Golovanova LV, Lalueza-Fox C, de la Rasilla M, Fortea J, Rosas A, Schmitz RW, Johnson PLF, Eichler EE, Falush D, Birney E, Mullikin JC, Slatkin M, Nielsen R, Kelso J, Lachmann M, Reich D, Paabo S (2010) A draft sequence of the Neandertal genome. Science 328, 710-722. https://doi.org/10.1126/science.1188021

Guo W, Sun D, Cao Y, Xiao L, Huang X, Ren W, Xu S, Yang G (2021) Extensive Interspecific Gene Flow Shaped Complex Evolutionary History and Underestimated Species Diversity in Rapidly Radiated Dolphins. Journal of Mammalian Evolution 29, 353-367. https://doi.org/10.1007/s10914-021-09581-6

Hibbins MS, Hahn MW (2022) Phylogenomic approaches to detecting and characterizing introgression. Genetics 220, iyab173. https://doi.org/10.1093/genetics/iyab173

Jiao X, Flouri T, Rannala B, Yang Z (2020) The Impact of Cross-Species Gene Flow on Species Tree Estimation. Syst Biol 69, 830-847. https://doi.org/10.1093/sysbio/syaa001

Li H, Durbin R (2011) Inference of human population history from individual whole-genome sequences. Nature 475, 493-496. https://doi.org/10.1038/nature10231

Malinsky M, Matschiner M, Svardal H (2021) Dsuite - Fast D-statistics and related admixture evidence from VCF files. Mol Ecol Resour 21, 584-595. https://doi.org/10.1111/1755-0998.13265

Malinsky M, Svardal H, Tyers AM, Miska EA, Genner MJ, Turner GF, Durbin R (2018) Whole-genome sequences of Malawi cichlids reveal multiple radiations interconnected by gene flow. Nature Ecology & Evolution 2, 1940-1955. https://doi.org/10.1038/s41559-018-0717-x

Mallet J, Besansky N, Hahn MW (2016) How reticulated are species? Bioessays 38, 140-149. https://doi.org/10.1002/bies.201500149 

McGowen MR, Tsagkogeorga G, Alvarez-Carretero S, Dos Reis M, Struebig M, Deaville R, Jepson PD, Jarman S, Polanowski A, Morin PA, Rossiter SJ (2020) Phylogenomic Resolution of the Cetacean Tree of Life Using Target Sequence Capture. Syst Biol 69, 479-501. https://doi.org/10.1093/sysbio/syz068

Moura AE, Shreves K, Pilot M, Andrews KR, Moore DM, Kishida T, Möller L, Natoli A, Gaspari S, McGowen M, Chen I, Gray H, Gore M, Culloch RM, Kiani MS, Willson MS, Bulushi A, Collins T, Baldwin R, Willson A, Minton G, Ponnampalam L, Hoelzel AR (2020) Phylogenomics of the genus Tursiops and closely related Delphininae reveals extensive reticulation among lineages and provides inference about eco-evolutionary drivers. Molecular Phylogenetics and Evolution 146,107047. https://doi.org/10.1016/j.ympev.2020.106756

Müller NF, Ogilvie HA, Zhang C, Fontaine MC, Amaya-Romero JE, Drummond AJ, Stadler T (in press) Joint inference of species histories and gene flow. Syst Biol.

Pease JB, Hahn MW (2015) Detection and Polarization of Introgression in a Five-Taxon Phylogeny. Syst Biol 64, 651-662. https://doi.org/10.1093/sysbio/syv023

Rabier CE, Berry V, Stoltz M, Santos JD, Wang W, Glaszmann JC, Pardi F, Scornavacca C (2021) On the inference of complex phylogenetic networks by Markov Chain Monte-Carlo. PLoS Comput Biol 17, e1008380. https://doi.org/10.1371/journal.pcbi.1008380

Seehausen O (2004) Hybridization and adaptive radiation. Trends Ecol Evol 19, 198-207. https://doi.org/10.1016/j.tree.2004.01.003

Solís-Lemus C, Yang M, Ané C (2016) Inconsistency of Species Tree Methods under Gene Flow. Syst Biol 65, 843-851. https://doi.org/10.1093/sysbio/syw030

Thawornwattana Y, Dalquen D, Yang Z, Tamura K (2018) Coalescent Analysis of Phylogenomic Data Confidently Resolves the Species Relationships in the Anopheles gambiae Species Complex. Molecular Biology and Evolution 35, 2512-2527. https://doi.org/10.1093/molbev/msy158

Wang K, Mathieson I, O’Connell J, Schiffels S (2020) Tracking human population structure through time from whole genome sequences. PLOS Genetics 16, e1008552. https://doi.org/10.1371/journal.pgen.1008552

Westbury MV, Cabrera AA, Rey-Iglesia A, Cahsan BD, Duchêne DA, Hartmann S, Lorenzen ED (2022) A genomic assessment of the marine-speciation paradox within the toothed whale superfamily Delphinoidea. bioRxiv, 2020.10.23.352286, ver. 7 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.10.23.352286

Westbury MV, Petersen B, Lorenzen ED (2019) Genomic analyses reveal an absence of contemporary introgressive admixture between fin whales and blue whales, despite known hybrids. PLoS ONE 14, e0222004. https://doi.org/10.1371/journal.pone.0222004

Yang Z (2015) The BPP program for species tree estimation and species delimitation. Current Zoology 61, 854-865. https://doi.org/10.1093/czoolo/61.5.854

Zhang C, Rabiee M, Sayyari E, Mirarab S (2018) ASTRAL-III: polynomial time species tree reconstruction from partially resolved gene trees. BMC Bioinformatics 19, 153. https://doi.org/10.1186/s12859-018-2129-y