Scale is a big topic in ecology [1]. Environmental variation happens at particular scales. The typical scale at which organisms disperse is species-specific, but, as a first approximation, an ensemble of similar species, for instance, trees, could be considered to share a typical dispersal scale. Finally, characteristic spatial scales of species interactions are, in general, different from the typical scales of dispersal and environmental variation. Therefore, conceptually, we can distinguish these three characteristic spatial scales associated with three different processes: species selection for a given environment (E), dispersal (D), and species interactions (I), respectively.  

From the famous species-area relation to the spatial distribution of biomass and species richness, the different macro-ecological patterns we usually study emerge from an interplay between dispersal and local interactions in a physical environment that constrains species establishment and persistence in every location. To make things even more complicated, local environments are often modified by the species that thrive in them, which establishes feedback loops.  It is usually assumed that local interactions are short-range in comparison with species dispersal, and dispersal scales are typically smaller than the scales at which the environment varies (I < D < E, see [2]), but this should not always be the case. 

The authors of this paper [2] relax this typical assumption and develop a theoretical framework to study how diversity and ecosystem functioning are affected by different relations between the typical scales governing interactions, dispersal, and environmental variation. This is a huge challenge. First, diversity and ecosystem functioning across space and time have been empirically characterized through a wide variety of macro-ecological patterns. Second, accommodating local interactions, dispersal and environmental variation and species environmental preferences to model spatiotemporal dynamics of full ecological communities can be done also in a lot of different ways. One can ask if the particular approach suggested by the authors is the best choice in the sense of producing robust results, this is, results that would be predicted by alternative modeling approaches and mathematical analyses [3]. The recommendation here is to read through and judge by yourself.  

The main unusual assumption underlying the model suggested by the authors is non-local species interactions. They introduce interaction kernels to weigh the strength of the ecological interaction with distance, which gives rise to a system of coupled integro-differential equations. This kernel is the key component that allows for control and varies the scale of ecological interactions. Although this is not new in ecology [4], and certainly has a long tradition in physics ---think about the electric or the gravity field, this approach has been widely overlooked in the development of the set of theoretical frameworks we have been using over and over again in community ecology, such as the Lotka-Volterra equations or, more recently, the metacommunity concept [5].

In Physics, classic fields have been revised to account for the fact that information cannot travel faster than light. In an analogous way, a focal individual cannot feel the presence of distant neighbors instantaneously. Therefore, non-local interactions do not exist in ecological communities. As the authors of this paper point out, they emerge in an effective way as a result of non-random movements, for instance, when individuals go regularly back and forth between environments (see [6], for an application to infectious diseases), or even migrate between regions. And, on top of this type of movement, species also tend to disperse and colonize close (or far) environments. Individual mobility and dispersal are then two types of movements, characterized by different spatial-temporal scales in general. Species dispersal, on the one hand, and individual directed movements underlying species interactions, on the other, are themselves diverse across species, but it is clear that they exist and belong to two distinct categories. 

In spite of the long and rich exchange between the authors' team and the reviewers, it was not finally clear (at least, to me and to one of the reviewers) whether the model for the spatio-temporal dynamics of the ecological community (see Eq (1) in [2]) is only presented as a coupled system of integro-differential equations on a continuous landscape for pedagogical reasons, but then modeled on a discrete regular grid for computational convenience. In the latter case, the system represents a regular network of local communities,  becomes a system of coupled ODEs, and can be numerically integrated through the use of standard algorithms. By contrast,  in the former case, the system is meant to truly represent a community that develops on continuous time and space, as in reaction-diffusion systems. In that case, one should keep in mind that numerical instabilities can arise as an artifact when integrating both local and non-local spatio-temporal systems. Spatial patterns could be then transient or simply result from these instabilities. Therefore, when analyzing spatiotemporal integro-differential equations, special attention should be paid to the use of the right numerical algorithms. The authors share all their code at https://zenodo.org/record/5543191, and all this can be checked out. In any case, the whole discussion between the authors and the reviewers has inherent value in itself, because it touches on several limitations and/or strengths of the author's approach,  and I highly recommend checking it out and reading it through.

Beyond these methodological issues, extensive model explorations for the different parameter combinations are presented. Several results are reported, but, in practice, what is then the main conclusion we could highlight here among all of them?  The authors suggest that "it will be difficult to manage landscapes to preserve biodiversity and ecosystem functioning simultaneously, despite their causative relationship", because, first, "increasing dispersal and interaction scales had opposing
effects" on these two patterns, and, second, unexpectedly, "ecosystems attained the highest biomass in scenarios which also led to the lowest levels of biodiversity". If these results come to be fully robust, this is, they pass all checks by other research teams trying to reproduce them using alternative approaches, we will have to accept that we should preserve biodiversity on its own rights and not because it enhances ecosystem functioning or provides particular beneficial services to humans. 

References

[1] Levin, S. A. 1992. The problem of pattern and scale in ecology. Ecology 73:1943–1967. https://doi.org/10.2307/1941447

[2] Yuval R. Zelnik, Matthieu Barbier, David W. Shanafelt, Michel Loreau, Rachel M. Germain. 2023. Linking intrinsic scales of ecological processes to characteristic scales of biodiversity and functioning patterns. bioRxiv, ver. 2 peer-reviewed and recommended by Peer Community in Ecology.  https://doi.org/10.1101/2021.10.11.463913

[3] Baron, J. W. and Galla, T. 2020. Dispersal-induced instability in complex ecosystems. Nature Communications  11, 6032. https://doi.org/10.1038/s41467-020-19824-4

[4] Cushing, J. M. 1977. Integrodifferential equations and delay models in population dynamics 
 Springer-Verlag, Berlin. https://doi.org/10.1007/978-3-642-93073-7

[5] M. A. Leibold, M. Holyoak, N. Mouquet, P. Amarasekare, J. M. Chase, M. F. Hoopes, R. D. Holt, J. B. Shurin, R. Law, D. Tilman, M. Loreau, A. Gonzalez. 2004. The metacommunity concept: a framework for multi-scale community ecology. Ecology Letters, 7(7): 601-613. https://doi.org/10.1111/j.1461-0248.2004.00608.x

[6] M. Pardo-Araujo, D. García-García, D. Alonso, and F. Bartumeus. 2023. Epidemic thresholds and human mobility. Scientific reports 13 (1), 11409. https://doi.org/10.1038/s41598-023-38395-0

In most vertebrate species, the period of life spanning from departure from the growing site until reaching a more advanced life stage (immature or adult) is critical. During this period, juveniles are often highly vulnerable because they have not reached the morphological, physiological and behavioural maturity levels of adults yet and are therefore at high risk of mortality, e.g. through starvation, depredation or competition (e.g. Marchetti & Price 1989, Wunderle 1991, Naef-Daenzer & Grüebler 2016). In line with this, juvenile survival is most often far lower than adult survival (e.g. Wooller et al. 1992). In species with parental care, juveniles have to acquire behavioural independence from their parents and possibly establish their own territory during this period of life. Very often, this is also the period that is least well-known in the life cycle (Cox et al. 2014, Naef-Daenzer & Grüebler 2016) because of reduced accessibility to individuals and/or adoption of low conspicuous behaviours. Therefore, our understanding of how juveniles acquire typical adult behaviours and how this progressively increases their survival prospects is still very limited (Naef-Daenzer & Grüebler 2016), and questions such as the length of this transition period or the cognitive (e.g. learning, memorization) mechanisms involved remain largely unresolved. This is particularly true regarding the acquisition of independent foraging behaviour (Marchetti & Price 1989).

Because direct observations of juvenile behaviours are usually very difficult except in specific situations or at the cost of an enormous effort, the use of remote tracking devices can be particularly appealing in this context (e.g. Ponchon et al. 2013, Kays et al. 2015). Over the past decades, technical advances have allowed the monitoring of not only individuals’ movements at both large and small spatial scales but also their activities and behaviours based on different parameters recording e.g. speed of movement or diving depth (Whitford & Klimley 2019). Device miniaturization has in particular allowed smaller species to be equipped and/or longer periods of time to be monitored (e.g. Naef-Daenzer et al. 2005). This has opened up whole fields of research, and has been particularly used on marine seabirds. In these species, individuals are most often inaccessible when at sea, representing most of the time outside (and even within) the breeding season, and the life cycle of these long-lived species can include an extended immature period (up to many years) during which most of them will remain unseen, until they come back as breeders or pre-breeders (e.g. Wooller et al. 1992, Oro & Martínez-Abraín 2009). Survival has been found to increase gradually with age in these species before reaching high values characteristic of the adult stage. However, the mechanisms underlying this increase are still to be deciphered.

The study by Delord et al. (2023) builds upon the hypothesis that juveniles gradually learn foraging techniques and movement strategies, improving their foraging efficiency, as previous data on flight parameters seemed to show in different long-lived bird species. Yet, these previous studies obtained data over a limited period of time, i.e. a few months at best. Whether these data could capture the whole dynamics of the progressive acquisition of foraging and movement skills can only be assessed by measuring behaviour over a longer time period and comparing it to similar data in adults, to account for seasonal variation in relation to both resource availability and energetic demands, e.g. due to molt.

The present study (Delord et al. 2023) addresses these questions by taking advantage of longer-lasting recordings of the location and activity of juvenile, immature and adult birds obtained simultaneously to investigate changes over time in juvenile behaviour and thereby provide hints about how young progressively acquire foraging skills. This study is performed on Amsterdam albatrosses, a highly endangered long-lived sea bird, with obvious conservation issues (Thiebot et al. 2015). The results show progressive changes in foraging effort over the first two months after departure from the birth colony, but large differences remain between life stages over a much longer time frame. They also reveal strong variations between sexes and over time in the year. Overall, this study, therefore, confirms the need for very long-term data to be collected in order to address the question of progressive behavioural maturation and associated survival consequences in such species with strongly deferred maturity. Ideally, the same individuals should be monitored over different life stages, from the juvenile period up to adulthood, but this would require further technical development to release the issue of powering duration limitation.

As reviewers emphasized in the first review round, one main challenge now remains to ascertain the outcome of the observed behavioural changes in foraging behaviour: we expect them to reflect improvement in foraging skills and thus performance of juveniles over time, but this would need to be tested. Collecting data on foraging efficiency is yet another challenge, that future technical developments may also help overcome. Importantly also, data were available only for individuals that could be caught again because the tracking device had to be retrieved from the bird. Here, a substantial fraction of the loggers (one-fifth) could not be found again (Delord et al. 2023). To what extent the birds for which no data could be obtained are a random sample of the equipped birds would also need to be assessed. The further development of remote tracking techniques allowing data to be downloaded from a long distance should help further exploration of behavioural ontogeny of juveniles while maturing and its survival consequences. Because the maturation process explored here is likely to show very different characteristics (e.g. timing and speed) in smaller / shorter-lived species (see Cox et al. 2014, Naef-Daenzer & Grüebler 2016), the development of miniaturization is also expected to allow further investigation of post-fledging behavioural maturation in a wider range of bird species. Our understanding of this crucial life phase in different types of species should thus continue to progress in the coming years.

References

Cox W. A., Thompson F. R. III, Cox A. S. & Faaborg J. 2014. Post-fledging survival in passerine birds and the value of post-fledging studies to conservation. Journal of Wildlife Management, 78: 183-193. https://doi.org/10.1002/jwmg.670

Delord K., Weimerskirch H. & Barbraud C. 2023. The challenges of independence: ontogeny of at-sea behaviour in a long-lived seabird. bioRxiv, ver. 6 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2021.10.23.465439

Kays R., Crofoot M. C., Jetz W. & Wikelski M. 2015. Terrestrial animal tracking as an eye on life and planet. Science, 348 (6240). https://doi.org/10.1126/science.aaa2478

Marchetti K: & Price T. 1989. Differences in the foraging of juvenile and adult birds: the importance of developmental constraints. Biological Reviews, 64: 51-70. https://doi.org/10.1111/j.1469-185X.1989.tb00638.x

Naef-Daenzer B., Fruh D., Stalder M., Wetli P. & Weise E. 2005. Miniaturization (0.2 g) and evaluation of attachment techniques of telemetry transmitters. The Journal of Experimental Biology, 208: 4063–4068. https://doi.org/10.1242/jeb.01870

Naef-Daenzer B. & Grüebler M. U. 2016. Post-fledging survival of altricial birds: ecological determinants and adaptation. Journal of Field Ornithology, 87: 227-250. https://doi.org/10.1111/jofo.12157

Oro D. & Martínez-Abraín A. 2009. Ecology and behavior of seabirds. Marine Ecology, pp.364-389.

Ponchon A., Grémillet D., Doligez B., Chambert T., Tveera T., Gonzàles-Solìs J & Boulinier T. 2013. Tracking prospecting movements involved in breeding habitat selection: insights, pitfalls and perspectives. Methods in Ecology and Evolution, 4: 143-150. https://doi.org/10.1111/j.2041-210x.2012.00259.x

Thiebot J.-B., Delord K., Barbraud C., Marteau C. & Weimerskirch H. 2015. 167 individuals versus millions of hooks: bycatch mitigation in longline fisheries underlies conservation of Amsterdam albatrosses. Aquatic Conservation 26: 674-688. https://doi.org/10.1002/aqc.2578

Whitford M & Klimley A. P. An overview of behavioral, physiological, and environmental sensors used in animal biotelemetry and biologging studies. Animal Biotelemetry, 7: 26. https://doi.org/10.1186/s40317-019-0189-z

Wooller R.D., Bradley J. S. & Croxall J. P. 1992. Long-term population studies of seabirds. Trends in Ecology and Evolution, 7: 111-114. https://doi.org/10.1016/0169-5347(92)90143-y

Wunderle J. M. 1991. Age-specific foraging proficiency in birds. Current Ornithology, 8: 273-324.

In this study, the authors report a new method for estimating the abundance of the Pyrenean brown bear population. Precisely, the methodology involved aims to apply Pollock's closed robust design (PCRD) capture-recapture models to estimate population abundance and trends over time. Overall, the results encourage the use of PCRD to study populations' demographic rates, while minimizing biases due to inter-individual heterogeneity in detection probabilities.

Estimating the size and trends of animal population over time is essential for informing conservation status and management decision-making (Nichols & Williams 2006). This is particularly the case when the population is small, geographically scattered, and threatened. Although several methods can be used to estimate population abundance, they may be difficult to implement when individuals are rare, elusive, solitary, largely nocturnal, highly mobile, and/or occupy large home ranges in remote and/or rugged habitats. Moreover, in such standard methods,

• the population is assumed to be closed both geographically (no immigration nor emigration) and demographically (no births nor deaths) and
• all individuals are assumed to have identical detection probabilities regardless of their individual attributes (e.g., age, body mass, social status) and habitat features (home-range location and composition) (Otis et al. 1978).

However, these conditions are rarely met in real populations, such as wild mammals (e.g., Bellemain et al. 2005; Solbert et al. 2006), and therefore the risk of underestimating population size can rapidly increase because the assumption of perfect detection of all individuals in the population is violated.

Focusing on the critically endangered Pyrenean brown bear that was close to extinction in the mid-1990s, the study by Vanpe et al. (2022), uses protocol and opportunistic data to describe a statistical modeling exercise to construct mark-recapture histories from 2008 to 2020. Among the data, the authors collected non-invasive samples such as a mixture of hair and scat samples used for genetic identification, as well as photographic trap data of recognized individuals. These data are then analyzed in RMark to provide detection and survival estimates. The final model (i.e. PCRD capture-recapture) is then used to provide Bayesian population estimates. Results show a five-fold increase in population size between 2008 and 2020, from 13 to 66 individuals. Thus, this study represents the first published annual abundance and temporal trend estimates of the Pyrenean brown bear population since 2008.

Then, although the results emphasize that the PCRD estimates were broadly close to the MRS counts and had reasonably narrow associated 95% Credibility Intervals, they also highlight that the sampling effort is different according to individuals. Indeed, as expected, the detection of an individual depends on

• the intraspecific home range size variation that results in individuals that move the most being most likely to be detected and
• the mortality rate which is higher on cubs than on adults and subadults (due to infanticide by males, predation, death of the mother, or abandonment).

Overall, the PCRD capture-recapture modelling approach, involved in this study, provides robust estimates of abundance and demographic rates of the Pyrenean brown bear population (with associated uncertainty) while minimizing and considering bias due to inter-individual heterogeneity in detection probabilities.

The authors conclude that mark-recapture provides useful population estimates and urge wildlife ecologists and managers to use robust approaches, such as the RDPC capture-recapture model, when studying large mammal populations. This information is essential to inform management decisions and assess the conservation status of populations.

References

Bellemain, E.V.A., Swenson, J.E., Tallmon, D., Brunberg, S. and Taberlet, P. (2005). Estimating population size of elusive animals with DNA from hunter-collected feces: four methods for brown bears. Cons. Biol. 19(1), 150-161. https://doi.org/10.1111/j.1523-1739.2005.00549.x

Nichols, J.D. and Williams, B.K. (2006). Monitoring for conservation. Trends Ecol. Evol. 21(12), 668-673. https://doi.org/10.1016/j.tree.2006.08.007

Otis, D.L., Burnham, K.P., White, G.C. and Anderson, D.R. (1978). Statistical inference from capture data on closed animal populations. Wildlife Monographs (62), 3-135.

Solberg, K.H., Bellemain, E., Drageset, O.M., Taberlet, P. and Swenson, J.E. (2006). An evaluation of field and non-invasive genetic methods to estimate brown bear (Ursus arctos) population size. Biol. Conserv. 128(2), 158-168. https://doi.org/10.1016/j.biocon.2005.09.025

Vanpé C, Piédallu B, Quenette P-Y, Sentilles J, Queney G, Palazón S, Jordana IA, Jato R, Elósegui Irurtia MM, de la Torre JS, and Gimenez O (2022) Estimating abundance of a recovering transboundary brown bear population with capture-recapture models. bioRxiv, 2021.12.08.471719, ver. 4 recommended and peer-reviewed by PCI Ecology. https://doi.org/10.1101/2021.12.08.471719

Humans are dramatically modifying many aspects of our planet via increasing concentrations of carbon dioxide in the atmosphere, patterns of land-use change, and unsustainable exploitation of the planet’s resources. These changes impact the abundance of species of wild organisms, with winners and losers. Identifying how different species and groups of species are influenced by anthropogenic activity in different biomes, continents, and habitats, has become a pressing scientific question with many publications reporting analyses of disparate data on species population sizes. Many conclusions are based on the linear analysis of rather short time series of organismal abundances.
 
There has been particular interest in how arthropods are impacted by environmental change, with several recent papers reporting contradictory results. To investigate why these contradictions might arise, Duchenne et al. (2022) conducted an analysis of four published data sets along with a series of experimental analyses of simulated time series to examine the power of widely used statistical analyses to gain inference on temporal trends. Their important paper reveals that accurate inference on dynamics, particularly of species that exhibit large temporal fluctuations in abundance, requires time series that are substantially longer than are typically collected, as well as careful thought as to whether linear models are appropriate. Linear analyses of short time series are susceptible to providing unreliable inference as trends can be strongly influenced by points at either end of the time series. 
 
Duchenne et al.’s paper provides important insight on the conditions when strong inference on temporal trends of arthropod (and other species) abundances can be made, and when they should be treated with caution. They do not doubt that many insect and arachnid species are changing their abundances, and that patterns in these changes may vary spatially. What their results do say is that we should treat grand claims of population recovery or rapid declines apparently to extinction with caution when they are based on short time series, particularly of species that show significant boom and bust dynamics. In many ways, these results are not unexpected, but it is nice to see such careful and thoughtful analyses and interpretation. More data are required for most arthropod species before clear assessments of abundance trends can be made. Given our reliance on many arthropods for food, pollination, and numerous ecosystem services, and the ability of other species to spread devastating human diseases such as dengue and malaria, it is advisable that we slow our modification of their habitats while additional data are collected to allow us to better characterise the trajectory of arthropod populations to understand what the consequences of our actions on the natural world are likely to be.  
 
References

Duchenne F, Porcher E, Mihoub J-B, Loïs G, Fontaine C (2022) Controversy over the decline of arthropods: a matter of temporal baseline? bioRxiv, 2022.02.09.479422, ver. 3 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.02.09.479422

"Good and evil are so close as to be chained together in the soul"
Robert Louis Stevenson, Dr. Jekyll and Mr. Hyde

Maternally inherited symbionts—microorganisms that pass from a female host to her progeny—have two main ways of increasing their own fitness. First, they can increase the fecundity or viability of infected females. This “positive fitness effects” strategy is the one commonly used by mutualistic symbionts, such as Buchnera aphidicola—the bacterial endosymbiont of the pea aphid, Acyrthosiphon pisum [4]. Second, maternally inherited symbionts can manipulate the reproduction of infected females in a way that enhances symbiont transmission at the expense of host fitness. A famous example of this “reproductive parasitism” strategy is the cytoplasmic incompatibility (CI) [3] induced by bacteria of the genus Wolbachia in their arthropod and nematode hosts. CI works as a toxin-antidote system, whereby the sperm of infected males is modified in a lethal way (toxin) that can only be reverted if the egg is also infected (antidote) [1]. As a result, CI imposes a kind of conditional sterility on their hosts: while infected females are compatible with both infected and uninfected males, uninfected females experience high offspring mortality if (and only if) they mate with infected males [7].

These two symbiont strategies (positive fitness effects versus reproductive parasitism) have been traditionally studied separately, both empirically and theoretically. However, it has become clear that the two strategies are not mutually exclusive, and that a reproductive parasite can simultaneously act as a mutualist—an infection type that has been dubbed “Jekyll and Hyde” [6], after the famous novella by Robert Louis Stevenson about kind scientist Dr. Jekyll and his evil alter ego, Mr. Hyde. In important previous work, Zug and Hammerstein [7] analyzed the consequences of positive fitness effects on the dynamics of different kind of infections, including “Jekyll and Hyde” infections characterized by CI and other reproductive parasitism strategies. Building on this and related modeling framework, Karisto et al. [2] re-investigate and expand on the interplay between positive fitness effects and reproductive parasitism in Wolbachia infections by focusing on CI in both diplodiploid and haplodiploid populations, and by paying particular attention to the mathematical assumption structure underlying their results.

Karisto et al. begin by reviewing classic models of Wolbachia infections in diplodiploid populations that assume a “negative fitness effect” (modeled as a fertility penalty on infected females), characteristic of a pure strategy of reproductive parasitism. Together with the positive frequency-dependent effects due to CI (whereby the fitness benefits to symbionts infecting females increase with the proportion of infected males in the population) this results in population dynamics characterized by two stable equilibria (the Wolbachia-free state and an interior equilibrium with a high frequency of Wolbachia-carrying hosts) separated by an unstable interior equilibrium. Wolbachia can then spread once the initial frequency is above a threshold or an invasion barrier, but is prevented from fixing by a proportion of infections failing to be passed on to offspring. Karisto et al. show that, given the assumption of negative fitness effects, the stable interior equilibrium can never feature a Wolbachia prevalence below one-half. Moreover, they convincingly argue that a prevalence greater than but close to one-half is difficult to maintain in the presence of stochastic fluctuations, as in these cases the high-prevalence stable equilibrium would be too close to the unstable equilibrium signposting the invasion barrier.

Karisto et al. then relax the assumption of negative fitness effects and allow for positive fitness effects (modeled as a fertility premium on infected females) in a diplodiploid population. They show that positive fitness effects may result in situations where the original invasion threshold is now absent, the bistable coexistence dynamics are transformed into purely co-existence dynamics, and Wolbachia symbionts can now invade when rare. Karisto et al. conclude that positive fitness effects provide a plausible and potentially testable explanation for the low frequencies of symbiont-carrying hosts that are sometimes observed in nature, which are difficult to reconcile with the assumption of negative fitness effects. 

Finally, Karisto et al. extend their analysis to haplodiploid host populations (where all fertilized eggs develop as females). Here, they investigate two types of cytoplasmic incompatibility: a female-killing effect, similar to the CI effect studied in diplodiploid populations (the “Leptopilina type” of Vavre et al. [5]) and a masculinization effect, where CI leads to the loss of paternal chromosomes and to the development of the offspring as a male (the “Nasonia type” of Vavre et al. [5]). The models are now two-sex, which precludes a complete analytical treatment, in particular regarding the stability of fixed points. Karisto et al. compensate by conducting large numerical analyses that support their claims. Importantly, all main conclusions regarding the interplay between positive fitness effects and reproductive parasitism continue to hold under haplodiploidy. 

All in all, the analysis and results by Karisto et al. suggest that it is not necessary to resort to classical (but depending on the situation, unlikely) mechanisms, such as ongoing invasion or source-sink dynamics, to explain arthropod populations featuring low-prevalent Wolbachia infections. Instead, low-frequency equilibria might be simply due to reproductive parasites conferring beneficial fitness effects, or Wolbachia symbionts playing Dr. Jekyll (positive fitness effects) and Mr. Hyde (cytoplasmatic incompatibility). 

References

[1] Beckmann JF, Bonneau M, Chen H, Hochstrasser M, Poinsot D, Merçot H, Weill M, Sicard M, Charlat S (2019) The Toxin–Antidote Model of Cytoplasmic Incompatibility: Genetics and Evolutionary Implications. Trends in Genetics, 35, 175–185. https://doi.org/10.1016/j.tig.2018.12.004

[2] Karisto P, Duplouy A, Vries C de, Kokko H (2022) Positive fitness effects help explain the broad range of Wolbachia prevalences in natural populations. bioRxiv, 2022.04.11.487824, ver. 5 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/2022.04.11.487824

[3] Laven H (1956) Cytoplasmic Inheritance in Culex. Nature, 177, 141–142. https://doi.org/10.1038/177141a0

[4] Perreau J, Zhang B, Maeda GP, Kirkpatrick M, Moran NA (2021) Strong within-host selection in a maternally inherited obligate symbiont: Buchnera and aphids. Proceedings of the National Academy of Sciences, 118, e2102467118. https://doi.org/10.1073/pnas.2102467118

[5] Vavre F, Fleury F, Varaldi J, Fouillet P, Bouletreau M (2000) Evidence for Female Mortality in Wolbachia-Mediated Cytoplasmic Incompatibility in Haplodiploid Insects: Epidemiologic and Evolutionary Consequences. Evolution, 54, 191–200. https://doi.org/10.1111/j.0014-3820.2000.tb00019.x

[6] Zug R, Hammerstein P (2015) Bad guys turned nice? A critical assessment of Wolbachia mutualisms in arthropod hosts. Biological Reviews, 90, 89–111. https://doi.org/10.1111/brv.12098

[7] Zug R, Hammerstein P (2018) Evolution of reproductive parasites with direct fitness benefits. Heredity, 120, 266–281. https://doi.org/10.1038/s41437-017-0022-5