The long-distance transport and introduction of new species by humans is increasingly leading divergent lineages to interact, and sometimes interbreed, even after thousands or millions of years of separation. It is thus of prime importance to understand the consequences of these contemporary admixture events on the evolutionary fitness of interacting organisms, and their ecological implications.

Ciona robusta and Ciona intestinalis are two species of sea squirts that diverged between 1.5 and 2 million years ago and recently came into contact again. This occurred through human-mediated introduction of C. robusta (native to the Northwest Pacific) into the range of C. intestinalis (the English channeled Northeast Atlantic). In this study, Fraïsse et al. (2022) follow up on earlier work by Le Moan et al. (2021), who had identified a long genomic hotspot of introgression of C. robusta ancestry segments in chromosome 5 of C. intestinalis. The hotspot bears suggestive evidence of positive selection and the authors aimed to investigate this further using fully phased whole-genome sequences.

The authors narrow down on the exact boundaries of the introgressed region, and make a compelling case that it has been the likely target of positive selection after introgression, using various complementary approaches based on patterns of population differentiation, haplotype structure and local levels of diversity in the region. Using extensive demographic modeling, they also show that the introgression event was likely recent (approximately 75 years ago), and distinct from other tracts in the C. intestinalis genome that are likely a product of more ancient episodes of interbreeding in the past 30,000 years. Narrowing down on the potential drivers of selection, the authors show that candidate SNPs in the region overlap with the cytochrome family 2 subfamily U gene - involved in the detoxification of exogenous compounds - potentially reflecting adaptation to chemicals encountered in the sea squirt's environment. There also appears to be copy number variation at the candidate SNPs, which provides clues into the adaptation mechanism in the region.

All reviewers agreed that the work carried out by the authors is elegant and the results are robustly supported and well presented. In a round of reviews, various clarifications of the manuscript were suggested by the reviewers, including on the quality of the newly generated sequencing data, and some suggestions for qualifications on the conclusions reached by the authors as well as changes in the figures to increase their clarity. The authors addressed the different concerns of the reviewers, and the new version is much improved. 

This study into human-mediated introgression and its consequences for adaptation is, in my view, both well thought-out and executed. I therefore provide an enthusiastic recommendation of this manuscript.

References

Fraïsse C, Le Moan A, Roux C, Dubois G, Daguin-Thiébaut C, Gagnaire P-A, Viard F and Bierne N (2022) Introgression between highly divergent sea squirt genomes: an adaptive breakthrough? bioRxiv, 2022.03.22.485319, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.03.22.485319

Le Moan A, Roby C, Fraïsse C, Daguin-Thiébaut C, Bierne N, Viard F (2021) An introgression breakthrough left by an anthropogenic contact between two ascidians. Molecular Ecology, 30, 6718–6732. https://doi.org/10.1111/mec.16189

The unique challenges associated with living on an island often result in organisms displaying a specific suite of traits commonly referred to as “island syndrome” (Adler and Levins, 1994; Burns, 2019; Baeckens and Van Damme, 2020). Large phenotypic shifts such as changes in size (e.g. shifts to gigantism or dwarfism, Lomolino, 2005) or coloration (Doutrelant et al., 2016) abound in the literature. However, less obvious phenotypes may also play a key role in adaptation to islands.

One such trait, reduced immune function, has important implications for the future of island populations in the face of anthropogenic-induced changes. Due to lower parasite pressure caused by a less diverse and less virulent parasite population, island hosts may show a decrease in immune defenses (Beadell et al., 2006; Pérez‐Rodríguez et al., 2013). However, this hypothesis has been challenged, as many studies have found ambiguous or conflicting results (Matson, 2006; Illera et al., 2015).

While most previous work has examined various immunological parameters (e.g., antibody concentrations), here, Barthe et al. (2022) take the novel approach of examining molecular signatures of immune genes. Using comparative genomic data from 34 different species of birds the authors examine the ratio of synonymous substitutions (i.e., not changing an amino acid) to non-synonymous substitutions (i.e., changing an amino acid) in innate and acquired immune genes (Pn/Ps ratio). Because population sizes on islands are lower which will affect molecular evolution, they compare these results to data from 97 control genes.  Assuming relaxed selection on islands predicts that the difference between the Pn/Ps ratio of immune genes and of control genes (ΔPn/Ps) is greater in island species compared to mainland ones.

As with previous work the authors found that the results differ depending on the category of immune genes. Both forms of innate defense: beta-defensins and Toll-like receptors did not show higher ΔPn/Ps for island populations. As these genes still have a higher Pn/Ps than control genes, the authors argue these results are in line with these genes being under purifying selection but lacking an “island effect”. Instead, the authors argue that demographic effects (i.e., reductions in Ne) may lead to the decreased immunity documented in other studies. In contrast, there was a reduction in Pn/Ps in MHC II genes, known to be under balancing selection. This reduction was stronger in island species and thus the authors argue that this is the only class of genes where a role for relaxed selection can be invoked. 

Together these results demonstrate that the changes in immunity experienced by island species are complex and that different categories of immune genes can experience different selective pressures. By including control genes in their study, they particularly highlight the importance of accounting for shifts in Ne when examining patterns of island species evolution. Hopefully, this kind of framework will be applied to other taxa to determine if these results are widespread or more specific to birds. 

References

Adler GH, Levins R (1994) The Island Syndrome in Rodent Populations. The Quarterly Review of Biology, 69, 473–490. https://doi.org/10.1086/418744

Baeckens S, Van Damme R (2020) The island syndrome. Current Biology, 30, R338–R339. https://doi.org/10.1016/j.cub.2020.03.029

Barthe M, Doutrelant C, Covas R, Melo M, Illera JC, Tilak M-K, Colombier C, Leroy T, Loiseau C, Nabholz B (2022) Evolution of immune genes in island birds: reduction in population sizes can explain island syndrome. bioRxiv, 2021.11.21.469450, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.11.21.469450

Beadell JS, Ishtiaq F, Covas R, Melo M, Warren BH, Atkinson CT, Bensch S, Graves GR, Jhala YV, Peirce MA, Rahmani AR, Fonseca DM, Fleischer RC (2006) Global phylogeographic limits of Hawaii’s avian malaria. Proceedings of the Royal Society B: Biological Sciences, 273, 2935–2944. https://doi.org/10.1098/rspb.2006.3671

Burns KC (2019) Evolution in Isolation: The Search for an Island Syndrome in Plants. Cambridge University Press, Cambridge. https://doi.org/10.1017/9781108379953

Doutrelant C, Paquet M, Renoult JP, Grégoire A, Crochet P-A, Covas R (2016) Worldwide patterns of bird colouration on islands. Ecology Letters, 19, 537–545. https://doi.org/10.1111/ele.12588

Illera JC, Fernández-Álvarez Á, Hernández-Flores CN, Foronda P (2015) Unforeseen biogeographical patterns in a multiple parasite system in Macaronesia. Journal of Biogeography, 42, 1858–1870. https://doi.org/10.1111/jbi.12548

Lomolino MV (2005) Body size evolution in insular vertebrates: generality of the island rule. Journal of Biogeography, 32, 1683–1699. https://doi.org/10.1111/j.1365-2699.2005.01314.x

Matson KD (2006) Are there differences in immune function between continental and insular birds? Proceedings of the Royal Society B: Biological Sciences, 273, 2267–2274. https://doi.org/10.1098/rspb.2006.3590

Pérez-Rodríguez A, Ramírez Á, Richardson DS, Pérez-Tris J (2013) Evolution of parasite island syndromes without long-term host population isolation: parasite dynamics in Macaronesian blackcaps Sylvia atricapilla. Global Ecology and Biogeography, 22, 1272–1281. https://doi.org/10.1111/geb.12084

Evolutionary genomics, and population genetics in particular, aim to decipher the respective influence of neutral and selective forces shaping genetic polymorphism in a species/population. This is a much-needed requirement before scanning genome data for footprints of species adaptation to their biotic and abiotic environment (Johri et al. 2022). In general, we would like to quantify the proportion of the genome evolving neutrally and under selective (positive, balancing and negative) pressures (Kern and Hahn 2018, Johri et al. 2021). We thus need to understand patterns of linked selection along the genome, that is how the distribution of genetic polymorphisms is shaped by selected sites and the recombination landscape. The present contribution by Pavinato et al. (2022) provides an additional method in the population genomics toolbox to quantify the extent of linked positive and negative selection using temporal data.

The availability of genomics data for model and non-model species has led to improvement of the modeling framework for demography and selection (Johri et al. 2022), but also new inference methods making use of the full genome data based on the Sequential Markovian Coalescent (SMC, Li and Durbin 2011), Approximate Bayesian Computation (ABC, Jay et al. 2019), ABC and machine learning (Pudlo et al. 2016, Raynal et al. 2019) or Deep Learning (Sanchez et al. 2021). These methods are based on one sample in time and the use of the coalescent theory to reconstruct the past (demographic) history. However, it is also possible to obtain for many species temporal data sampled over several time points. For species with short generation time (in experimental evolution or monitored populations), one can sample a population every couple of generations as exemplified with Drosophila melanogaster (Bergland et al. 2010). For species with longer generation times that cannot be easily regularly sampled in time, it becomes possible to sequence available specimens from museums (e.g. Cridland et al. 2018) or ancient DNA samples. Methods using temporal data are based on the classical population genomics assumption that demography (migration, population subdivision, population size changes) leaves a genome-wide signal, while selection leaves a localized signal in the close vicinity of the causal mutation. Several methods do assess the demography of a population (change in effective population size, Ne, in time) using temporal data (e.g. Jorde and Ryman 2007) which can be used to calibrate the detection of loci under strong positive selection (Foll et al. 2014). Recently Buffalo and Coop (2020) used genome-wide covariance between allele frequency changes across time samples (and across replicates) to quantify the effects of linked selection over short timescales. 

In the present contribution, Pavinato et al. (2022) make use of temporal data to draw the joint estimation of demographic and selective parameters using a simulation-based method (ABC-Random Forests). This study by Pavinato et al. (2022) builds a framework allowing to infer the census size of the population in time (N) separately from the effect of genetic drift, which is determined by change in effective population size (Ne) in time, as well estimates of genome-wide parameters of selection. In a nutshell, the authors use a forward simulator and summarize genome data by genomic windows using classic statistics (nucleotide diversity, Tajima’s D, FST, heterozygosity) between time samples and for each sample. They specifically use the distributions (higher moments) of these statistics among all windows. The authors combine as input for the ABC-RF, vectors of summary statistics, model parameters and five latent variables: Ne, the ratio Ne/N, the number of beneficial mutations under strong selection, the average selection coefficient of strongly selected mutations, and the average substitution load. Indeed, the authors are interested in three different types of selection components: 1) the adaptive potential of a population which is estimated as the population mutation rate of beneficial mutations (θb), 2) the number of mutations under strong selection (irrespective of whether they reached fixation or not), and 3) the overall population fitness which is a function of the genetic load. In other words, the novelty of this method is not to focus on the detection of loci under selection, but to infer key parameters/distributions summarizing the genome-wide signal of demography and (positive and negative) selection. As a proof of principle, the authors then apply their method to a dataset of feral populations of honey bees (Apis mellifera) collected in California across many years and recovered from Museum samples (Cridland et al. 2018). The approach yields estimates of Ne which are on the same order of magnitude of previous estimates in hymenopterans, and the authors discuss why the different populations show various values of Ne and N which can be explained by different history of admixture with wild but also domesticated lineages of bees.

This study focuses on quantifying the genome-wide joint footprints of demography, and strong positive and negative selection to determine which proportion of the genome evolves neutrally or not. Further application of this method can be anticipated, for example, to study species with ecological and life-history traits which generate discrepancies between census size and Ne, for example for plants with selfing or seed banking (Sellinger et al. 2020), and for which the genome-wide effect of linked selection is not fully understood.

References

Johri P, Aquadro CF, Beaumont M, Charlesworth B, Excoffier L, Eyre-Walker A, Keightley PD, Lynch M, McVean G, Payseur BA, Pfeifer SP, Stephan W, Jensen JD (2022) Recommendations for improving statistical inference in population genomics. PLOS Biology, 20, e3001669. https://doi.org/10.1371/journal.pbio.3001669

Kern AD, Hahn MW (2018) The Neutral Theory in Light of Natural Selection. Molecular Biology and Evolution, 35, 1366–1371. https://doi.org/10.1093/molbev/msy092

Johri P, Riall K, Becher H, Excoffier L, Charlesworth B, Jensen JD (2021) The Impact of Purifying and Background Selection on the Inference of Population History: Problems and Prospects. Molecular Biology and Evolution, 38, 2986–3003. https://doi.org/10.1093/molbev/msab050

Pavinato VAC, Mita SD, Marin J-M, Navascués M de (2022) Joint inference of adaptive and demographic history from temporal population genomic data. bioRxiv, 2021.03.12.435133, ver. 6 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.03.12.435133

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

Jay F, Boitard S, Austerlitz F (2019) An ABC Method for Whole-Genome Sequence Data: Inferring Paleolithic and Neolithic Human Expansions. Molecular Biology and Evolution, 36, 1565–1579. https://doi.org/10.1093/molbev/msz038

Pudlo P, Marin J-M, Estoup A, Cornuet J-M, Gautier M, Robert CP (2016) Reliable ABC model choice via random forests. Bioinformatics, 32, 859–866. https://doi.org/10.1093/bioinformatics/btv684

Raynal L, Marin J-M, Pudlo P, Ribatet M, Robert CP, Estoup A (2019) ABC random forests for Bayesian parameter inference. Bioinformatics, 35, 1720–1728. https://doi.org/10.1093/bioinformatics/bty867

Sanchez T, Cury J, Charpiat G, Jay F (2021) Deep learning for population size history inference: Design, comparison and combination with approximate Bayesian computation. Molecular Ecology Resources, 21, 2645–2660. https://doi.org/10.1111/1755-0998.13224

Bergland AO, Behrman EL, O’Brien KR, Schmidt PS, Petrov DA (2014) Genomic Evidence of Rapid and Stable Adaptive Oscillations over Seasonal Time Scales in Drosophila. PLOS Genetics, 10, e1004775. https://doi.org/10.1371/journal.pgen.1004775

Cridland JM, Ramirez SR, Dean CA, Sciligo A, Tsutsui ND (2018) Genome Sequencing of Museum Specimens Reveals Rapid Changes in the Genetic Composition of Honey Bees in California. Genome Biology and Evolution, 10, 458–472. https://doi.org/10.1093/gbe/evy007

Jorde PE, Ryman N (2007) Unbiased Estimator for Genetic Drift and Effective Population Size. Genetics, 177, 927–935. https://doi.org/10.1534/genetics.107.075481

Foll M, Shim H, Jensen JD (2015) WFABC: a Wright–Fisher ABC-based approach for inferring effective population sizes and selection coefficients from time-sampled data. Molecular Ecology Resources, 15, 87–98. https://doi.org/10.1111/1755-0998.12280

Buffalo V, Coop G (2020) Estimating the genome-wide contribution of selection to temporal allele frequency change. Proceedings of the National Academy of Sciences, 117, 20672–20680. https://doi.org/10.1073/pnas.1919039117

Sellinger TPP, Awad DA, Moest M, Tellier A (2020) Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data. PLOS Genetics, 16, e1008698. https://doi.org/10.1371/journal.pgen.1008698

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