Submit a preprint

Latest recommendationsrsstwitter

IdTitle * Authors * Abstract * Picture * Thematic fields * RecommenderReviewersSubmission date
08 Feb 2022
article picture

The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia

Experimental evidence for asymmetrical species recognition in East African Ophthalmotilapia cichlids

Recommended by based on reviews by George Turner and 2 anonymous reviewers

I recommend the Van Steenberge et al. study. With over 2000 endemic species, the East African cichlids are a well-established model system in speciation research (Salzburger 2018) and several models have been proposed and tested to explain how these radiations formed (Kocher 2004). Hybridization was shown to be a main driver of the rapid speciation and adaptive radiations of the East African Cichlid fishes (Seehausen 2004). However, it is obvious that unrestrained hybridization also has the potential to reduce taxonomic diversity by erasing species barriers. In the classical model of cichlid evolution, special emphasis was placed on mate preference (Kocher 2004). However, no attention was placed on species recognition, which was implicitly assumed. There is, however, more research needed on what species recognition means, especially in radiating lineages such as cichlids. In a previous study, Nevado et al. 2011 found traces of asymmetrical hybridization between members of the Lake Tanganyika radiation: the genus Ophthalmotilapia. This recommended study by Van Steenberge et al. is based on Nevado et al. (2011), which detected that in one genus of Ophthalmotilapia mitochondrial DNA ‘typical’ for one of the four species (O. nasuta) was also found in three other species (O. ventralis, O. heterodonta, and O. boops). The authors suggested that this could be explained by the fact that females of the three other species accepted O. nasuta males, but that O. nasuta females were more selective and accepted only conspecifc males. This could hence be due to asymmetric mate preferences, or by asymmetric abilities for species recognition. 

This is exactly what the current study by Van Steenberge et al. did. They tested the latter hypothesis by presenting females of two different Ophthalmotilapia species with con- and heterospecific males. This was tested through experiments, making use of wild specimens of two species: O. nasuta and O. ventralis. The authors assumed that if they performed classical “choice-experiments”, they would not notice the recognition effects, given that females would just select preferred, most likely conspecific, males. Instead, specimens were only briefly presented to other fishes since the authors wanted to compare differences in the ability for ‘species recognition’. In this, the authors followed Mendelson and Shaw (2012) who used “a measurable difference in behavioural response towards conspecifics as compared to heterospecifics’’ as a definition for recognition. Instead of the focus on selection/preference, they investigated if females of different species behaved differently, and hence detected the difference between conspecific and heterospecific males. This was tested by a short (15 minutes) exposure to another fish in an isolated part of the aquarium. Recognition was defined as the ‘difference in a particular behaviour between the two conditions’. What was monitored was the swimming behaviour and trajectory (1 image per second) together with known social behaviours of this genus. The selection of these behaviours was further facilitated based on experimental set-ups of reproductive behaviour or the same species previously described by the same research team (Kéver et al. 2018).

The result was that O. nasuta females, for which it was expected that they would not hybridize, showed a different behaviour towards a con- or a heterospecific male. They interacted less with males of the other species. What was unexpected is that there was no difference in behaviour of the females whether they recognized a male or (control) female of their own species. This suggests that they did not detect differences in reproductive behaviour, but rather in the interactions between conspecifics. For females of O. ventralis, for which there are indications for hybridization in the wild, they did not find a difference in behaviour. Females of this species behaved identically with respect to the right and wrong males as well as towards the control females. Interestingly is thus that a complex pattern between species in the wild could be (partially) explained by the behaviour/interaction at first impression of the individuals of these species. 

References

Kéver L, Parmentier E, Derycke S, Verheyen E, Snoeks J, Van Steenberge M, Poncin P (2018) Limited possibilities for prezygotic barriers in the reproductive behaviour of sympatric Ophthalmotilapia species (Teleostei, Cichlidae). Zoology, 126, 71–81. https://doi.org/10.1016/j.zool.2017.12.001

Kocher TD (2004) Adaptive evolution and explosive speciation: the cichlid fish model. Nature Reviews Genetics, 5, 288–298. https://doi.org/10.1038/nrg1316

Mendelson TC, Shaw KL (2012) The (mis)concept of species recognition. Trends in Ecology & Evolution, 27, 421–427. https://doi.org/10.1016/j.tree.2012.04.001

Nevado B, Fazalova V, Backeljau T, Hanssens M, Verheyen E (2011) Repeated Unidirectional Introgression of Nuclear and Mitochondrial DNA Between Four Congeneric Tanganyikan Cichlids. Molecular Biology and Evolution, 28, 2253–2267. https://doi.org/10.1093/molbev/msr043

Salzburger W (2018) Understanding explosive diversification through cichlid fish genomics. Nature Reviews Genetics, 19, 705–717. https://doi.org/10.1038/s41576-018-0043-9

Seehausen O (2004) Hybridization and adaptive radiation. Trends in Ecology & Evolution, 19, 198–207. https://doi.org/10.1016/j.tree.2004.01.003

Steenberge MV, Jublier N, Kéver L, Gresham S, Derycke S, Snoeks J, Parmentier E, Poncin P, Verheyen E (2022) The initial response of females towards congeneric males matches the propensity to hybridise in Ophthalmotilapia. bioRxiv, 2021.08.07.455508, ver. 3 peer-reviewed and recommended by Peer Community in Zoology. https://doi.org/10.1101/2021.08.07.455508

The initial response of females towards congeneric males matches the propensity to hybridise in OphthalmotilapiaMaarten Van Steenberge, Noemie Jublier, Loic Kever, Sophie Gresham, Sofie Derycke, Jos Snoeks, Eric Parmentier, Pascal Poncin, Erik Verheyen<p style="text-align: justify;">Cichlid radiations often harbour closely related species with overlapping niches and distribution ranges. Such species sometimes hybridise in nature, which raises the question how can they coexist. This also holds f...Aquatic, Behavior, Evolution, Fish, Vertebrates, Veterinary entomologyEllen Decaestecker2021-08-09 12:22:49 View
01 Jul 2020
article picture

Sub-lethal insecticide exposure affects host biting efficiency of Kdr-resistant Anopheles gambiae

kdr homozygous resistant An. gambiae displayed enhanced feeding success when exposed to permethrin Insect-Treated Nets

Recommended by based on reviews by Thomas Guillemaud, Niels Verhulst, Etienne Bilgo and 1 anonymous reviewer

Malaria is a vector-borne parasitic disease found in 91 countries with an estimated of 228 million cases occurred worldwide during 2018. The 93% (213 million) of those cases were reported in the African Region (WHO 2019). Six species of Plasmodium parasites can produce the disease but only P. falciparum and P. vivax are the predominant species globally. More than 40 species of Anopheles mosquitoes are important malaria vectors (Asley et al. 2018). Intrinsic (genetic background, parasite susceptibility) and extrinsic (feeding host preference, host diversity and availability, mosquito abundance) factors affect the capacity of mosquitoes to vector the disease (Macdonald 1952). Malaria is prevented by chemoprophylaxis, vaccination, bite-avoidance and vector-control measures. The mainstays of vector control are long-lasting insecticide (pyrethroid) treated nets and indoor residual spraying with insecticides (Asley et al. 2018). The widespread use of pyrethroid insecticides forced the emergence of insecticide resistance in malaria vectors reducing the insecticidal effect. Mosquitoes can modify their behaviour avoiding insecticide contact and so potentially reducing vector control tools efficacy. In this sense, Diop et al. (2020) investigated whether pre-exposure to an Insecticide-Treated Net (ITN) modulates the mosquito ability to take a blood meal in Anopheles gambiae. By means of video recording experiments the authors analyzed how the feeding/bitting behaviour was affected by kdr mutation genotypes (homozygous susceptible – SS-, heterozygotes -RS- and homozygous resistant -RR-) when exposed to two different insecticides (permethrin and deltamethrin). According to the results, the blood-feeding success did not differ between the three genotypes in the absence of insecticide exposure. However, authors observed differences in the feeding duration and blood meal size. In example, RR mosquitoes spent less time taking their blood meal than RS and SS. On the other hand, RS mosquitoes took higher blood volumes than RR females. These differences can affect the mosquito fitness by decreasing/increasing the likelihood to be killed by the host defensive behavior or increase the oogenesis so enhancing fecundity. Regarding the effect of exposition to insecticides authors detected a strong relationship between kdr genotype and Knock Down (KD) phenotype when mosquitoes were exposed to Permethrin. Previously, the authors have evidenced that RR mosquitoes prefer a host protected by a permethrin-treated net rather than an untreated net and that heterozygotes RS mosquitoes have a remarkable ability to find a hole into a bet net (Diop et al. 2015, Porciani et al. 2017). With data here obtained, they demonstrated that kdr homozygous resistant An. gambiae displayed enhanced feeding success when exposed to permethrin ITN. The changes observed in the feeding/biting mosquito behaviour can affect their fitness shaping the evolution of the insecticide resistance in mosquitoes’ natural populations. Moreover, this may also alter parasite transmission dynamics by modifying vector/host interactions and so vector capacity.

References

World Health Organization (2019). World malaria report 2019. Geneva: World Health Organization; 2019. ISBN 978-92-4-156572-1
Ashley EA, Pyae Phyo A, Woodrow CJ (2018). Malaria. Lancet. 391(10130):1608‐1621. doi: 10.1016/S0140-6736(18)30324-6
Macdonald G (1952). The analysis of equilibrium in malaria. Trop Dis Bull 49: 813-828.
Diop MM, Chandre F, Rossignol M, Porciani A, Château M, Moiroux N and Pennetier, C. (2020). Sub-lethal insecticide exposure affects host biting efficiency of Kdr-resistant Anopheles gambiae. bioRxiv 653980, ver. 4 peer-reviewed and recommended by PCI Zoology. doi: 10.1101/653980
Diop MM, Moiroux N, Chandre F, Martin-Herrou H, Milesi P, Boussari O, et al. (2015) Behavioral cost and overdominance in Anopheles gambiae. PLoS ONE. 10(4):e0121755. doi: 10.1371/journal.pone.0121755
Porciani A, Diop M, Moiroux N, Kadoke-Lambi T, Cohuet A, Chandre F, et al. (2017) Influence of pyrethroïd-treated bed net on host seeking behavior of Anopheles gambiae s.s. carrying the kdr allele. PLOS ONE. 12(7):e0164518. doi: 10.1371/journal.pone.0164518

Sub-lethal insecticide exposure affects host biting efficiency of Kdr-resistant Anopheles gambiaeMalal Mamadou Diop, Fabrice Chandre, Marie Rossignol, Angelique Porciani, Mathieu Chateau, Nicolas Moiroux, Cedric Pennetier<p>The massive use of insecticide-treated nets (ITNs) has drastically changed the environment for malaria vector mosquitoes, challenging their host-seeking behaviour and biting success. Here, we investigated the effect of a brief exposure to an IT...Behavior, Ecology, Evolution, Medical entomology, Pesticide resistanceAdrian DiazAnonymous2019-05-29 19:40:25 View
19 Jul 2024
article picture

Museomics of Carabus giant ground beetles shows an Oligocene origin and in situ Alpine diversification

Natural history collections continue to inform ground beetle genetics.

Recommended by based on reviews by Michael Caterino, Julian Dupuis and 1 anonymous reviewer

Some of the biodiversity of our planet now exists only in museums, due to continuing habitat destruction and climate change. With more than 380 million entomological specimens already preserved in museums (Johnson and Owens 2023), there is much work left to document what we already have. Fortunately, new advances in DNA sequencing have given us the opportunity to get enormous amounts of information from dried specimens on pins.
 
One such advance is HyRAD-X, which uses RAD-derived probes originally developed using RNA extracted from a selection of specimens with high RNA-integrity (Schmid et al. 2017). These exome-limited probes can then be used to capture low-integrity DNA extracted from a single leg from museum specimens, followed by Illumina sequencing of the enriched libraries.
 
Ground beetles allow an excellent demonstration of this approach, as their diversity, large size, and charismatic appearance has led to them being well represented in museums. Using a HyRAD-X probe set previously developed for a higher phylogeny within the subfamily Carabinae (Toussaint et al. 2021), the authors have now applied the same probe set to produce a comprehensive phylogeny for Arcifera, a clade of four subgenera and ten species within the genus Carabus (Pauli et al. 2024).
 
Of the 96 specimens that they started out with, 90% were from natural history collections and 40% dropped out immediately due to poor DNA extraction yield. After filtering the resulting sequence reads for minimum coverage and minimum number of samples per locus, they ended up with 35 museum specimens with an average of 793 loci. Phylogenetic analysis of this data supported the current classification of these beetles.
 
Pauli et al’s. (2024) study has effectively shown the power of HyRAD-X methods for applications at the species level. In-house production of probes makes the method accessible, expanding the opportunity to use museum specimens for population genetic research.

References

Johnson KR, Owens, (IFP. 2023) A global approach for natural history museum collections. Science 379,1192-1194(2023). https://doi.org/10.1126/science.adf6434

Pauli MT, Gauthier J, Labédan M, Blanc M, Bilat J, Toussaint EFA (2024) Museomics of Carabus giant ground beetles shows an Oligocene origin and in situ alpine diversification. bioRxiv, ver. 5 peer-reviewed and recommended by Peer Community in Zoology. https://doi.org/10.1101/2024.03.21.586057

Schmid, S., Genevest, R., Gobet, E., Suchan, T., Sperisen, C., Tinner, W. and Alvarez, N. (2017), HyRAD-X, a versatile method combining exome capture and RAD sequencing to extract genomic information from ancient DNA. Methods Ecol Evol, 8: 1374-1388. https://doi.org/10.1111/2041-210X.12785

Toussaint EFA, Gauthier J, Bilat J, Gillett CPDT, Gough HM, Lundkvist H, Blanc M, Muñoz-Ramírez CP, Alvarez N (2021) HyRAD-X Exome Capture Museomics Unravels Giant Ground Beetle Evolution, Genome Biology and Evolution, Volume 13, Issue 7, evab112, https://doi.org/10.1093/gbe/evab112

Museomics of *Carabus* giant ground beetles shows an Oligocene origin and *in situ* Alpine diversificationMarie T. PAULI, Jeremy GAUTHIER, Marjorie LABEDAN, Mickael BLANC, Julia BILAT, Emmanuel F.A. TOUSSAINT<p style="text-align: justify;">The development of museomics represents a major paradigm shift in the use of natural history collection specimens for systematics and evolutionary biology. New approaches in this field allow the sequencing of hundre...Insecta, Phylogeny, SystematicsFelix Sperling2024-03-27 15:30:31 View
26 Aug 2022
article picture

Within and among population differences in cuticular hydrocarbons in the seabird tick Ixodes uriae

Seabird tick diversification and cuticular hydrocarbons

Recommended by based on reviews by 2 anonymous reviewers

Ticks are notorious vectors of diseases in humans and other vertebrates. Much effort has been expended to understand tick diversity and ecology with the aim of managing their populations to alleviate the misery they bring. Further, the fundamental question of whether ticks are usually host generalists or host specialists has been debated at length and is important both for understanding the mechanisms of their diversification as well as for focusing control of ticks [1].

One elegant resolution of this question is to consider most tick species to be global generalists but local specialists [1]. This is well illustrated in a series of studies of the seabird tick, Ixodes uriae, which is comprised of host-specific races that show genetic [2], morphological [3] and host performance [4] differences associated with the seabirds they feed on. Such a pattern has clear ramifications for sympatric speciation; however, the factors that potentially act to drive these differences have remained elusive.

Dupraz et al. [5] have now made intriguing and important steps toward bridging the gap between demonstrating local patterns of tick host association and understanding the physiological mechanisms that may facilitate such divergences. They collected I. uriae ticks from the nests of two seabirds – Atlantic puffins and common guillemots – on the north side of Iceland. Four populations of ticks were sampled, with one island providing both puffin ticks and guillemot ticks, to give two tick populations from each of the two seabird host species. They then washed the ticks in solvent and analyzed the dissolved cuticular hydrocarbons (CHCs) using GC mass spectrometry, revealing 22 different hydrocarbon compounds common to most of these samples. CHCs are known to be important across arthropods for a variety of functions ranging from reducing water loss to facilitating communication and recognition between individuals with species.

Dupraz et al. [5] found three hydrocarbons that distinguished puffin ticks most consistently from guillemot ticks. A cross-validation test for host type also assigned 75% of the tick pools to the seabird host of origin. However, with these limited sample sizes, statistical analysis revealed no significant difference in CHC profiles between the host types, although a tendency was evident. Nonetheless, this study revealed a number of potentially diagnostic CHCs for tick host type, as well as some that may be more diagnostic of locations. This provides a fascinating and actionable foundation for further work using additional sites and host types, as well as an entry point into discerning the mechanisms at play in producing the diversity, complexity and adaptability that make ticks such medical menaces.

References

[1]  McCoy, K.D., Léger, E., Dietrich, M., 2013. Host specialization in ticks and transmission of tick-borne diseases: a review. Front. Cell. Infect. Microbiol. 3. https://doi.org/10.3389/fcimb.2013.00057

[2]  McCoy, K.D., Chapuis, E., Tirard, C., Boulinier, T., Michalakis, Y., Bohec, C.L., Maho, Y.L., Gauthier-Clerc, M., 2005. Recurrent evolution of host-specialized races in a globally distributed parasite. Proc. R. Soc. B Biol. Sci. 272, 2389–2395. https://doi.org/10.1098/rspb.2005.3230

[3]  Dietrich, M., Beati, L., Elguero, E., Boulinier, T., McCoy, K.D., 2013. Body size and shape evolution in host races of the tick Ixodes uriae. Biol. J. Linn. Soc. 108, 323–334. https://doi.org/10.1111/j.1095-8312.2012.02021.x

[4]  Dietrich, M., Lobato, E., Boulinier, T., McCoy, K.D., 2014. An experimental test of host specialization in a ubiquitous polar ectoparasite: a role for adaptation? J. Anim. Ecol. 83, 576–587. https://doi.org/10.1111/1365-2656.12170

[5] Dupraz, M., Leroy, C., Thórarinsson, T. L., d’Ettorre, P. and McCoy, K. D. (2022) Within and among population differences in cuticular hydrocarbons in the seabird tick Ixodes uriae. bioRxiv, 2022.01.21.477272, ver. 5 peer-reviewed and recommended by Peer Community in Zoology. https://doi.org/10.1101/2022.01.21.477272

Within and among population differences in cuticular hydrocarbons in the seabird tick Ixodes uriaeMarlène Dupraz, Chloe Leroy, Thorkell Lindberg Thórarinsson, Patrizia d’Ettorre, Karen D. McCoy<p>The hydrophobic layer of the arthropod cuticle acts to maintain water balance, but can also serve to transmit chemical signals via cuticular hydrocarbons (CHC), essential mediators of arthropod behavior. CHC signatures typically vary qualitativ...Acari, Biology, Ecology, EvolutionFelix Sperling2022-02-08 13:00:52 View
07 Jun 2024
article picture

Relationship between weapon size and six key behavioural and physiological traits in males of the European earwig

The unreliable signal: No correlation between forceps length and male quality in European earwigs

Recommended by ORCID_LOGO based on reviews by Luna Grey and 2 anonymous reviewers

In animals, male weapons such as antlers, horns, spurs, fangs, and tusks typically provide advantages in male contests and increase access to females, thereby enhancing reproductive success. However, such large and extravagant morphological structures are expected to come at a cost, potentially imposing trade-offs with life history traits, physiological functions, or certain behaviors (Emlen, 2001; Emlen, 2008). These costs should be manageable only by males in the best condition. The present study by Blackwell et al. (2024) examines this assumption through a comprehensive study on the European earwig, where males possess forceps-like cerci that vary widely in size within populations.

In the European earwig (Forficula auricularia), male forceps are used in male-male contests as weapons to deter competitors prior to mating (Styrsky & Rhein, 1999) or to interrupt mating pairs by non-copulating males (Forslund, 2000; Walker & Fell, 2001). Despite providing benefits in terms of mating success (Eberhard & Gutierrez, 1991; Tomkins & Brown, 2004), it remains unknown whether long or short forceps are associated with other important life-history traits.

In this laboratory study, Blackwell et al. (2024) investigated two European earwig populations, each divided into two subpopulations: one with the shortest forceps and one with the longest forceps. They examined the potential costs of long forceps on six different traits: one reproductive trait (sperm storage); three non-reproductive behavioral traits such as locomotor performance (involved in search for resources), fleeing reaction face to a risk (long forceps are supposed to be correlated with boldness), and aggregation behavior (European earwigs are facultative group-living organisms); and survival (when deprived of food and subsequently when exposed to an entomopathogenic fungus).

As males in the best condition are supposed to be those that can afford to develop large forceps, Blackwell et al. (2024) predicted that males with long forceps would perform better than those with short forceps across the investigated traits. However, their predictions were not validated, as no correlation between weapon size and male quality was detected in either population. Although the sample size is sometimes limited, the consistency of these results across different populations adds robustness to their conclusions.

By demonstrating that forceps length in the European earwig does not reliably indicate male quality, this paper challenges existing theories and highlights the complexity of evolutionary processes shaping morphological traits. Furthermore, the study raises important questions about the evolutionary mechanisms maintaining weapon size diversity, providing a fresh perspective that could stimulate further research and debate in the field, notably the search for other traits where costs might be incurred.

References

Blackwell, S.E.M., Pasquier, L., Dupont, S., Devers, S., Lécureuil, C. & Meunier, J. (2024). Relationship between weapon size and six key behavioural and physiological traits in males of the European earwig. bioRxiv, ver. 3 peer-reviewed and recommended by Peer Community in Zoology. https://doi.org/10.1101/2024.03.20.585871

Eberhard, W.G., & Gutierrez, E.E. (1991). Male dimorphisms in beetles and earwigs and the question of developmental constraints. Evolution, 45(1), 18–28. https://doi.org/10.2307/2409478

Emlen, D.J. (2001). Costs and the diversification of exaggerated animal structures. Science, 291(5508), 1534–1536. https://doi.org/10.1126/science.1056607

Emlen, D.J. (2008). The evolution of animal weapons. Annual Review of Ecology, Evolution, and Systematics, 39(1), 387–413. https://doi.org/10.1146/annurev.ecolsys.39.110707.173502

Forslund, P. (2000). Male-male competition and large size mating advantage in European earwigs, Forficula auricularia. Animal Behaviour, 59(4), 753–762. https://doi.org/10.1006/anbe.1999.1359

Styrsky, J.D., & Rhein, S.V. (1999). Forceps size does not determine fighting success in European earwigs. Journal of Insect Behavior, 12(4), 475–482. https://doi.org/10.1023/A:1020962606724

Tomkins, J.L., & Brown, G.S. (2004). Population density drives the local evolution of a threshold dimorphism. Nature, 431, 1099–1103. https://doi.org/10.1038/nature02936.1.

Walker, K.A., & Fell, R.D. (2001). Courtship roles of male and female European earwigs, Forficula auricularia L. (Dermaptera: Forficulidae), and sexual use of forceps. Journal of Insect Behavior, 14(1), 1–17. https://doi.org/10.1023/A:1007843227591

Relationship between weapon size and six key behavioural and physiological traits in males of the European earwigSamantha E.M. Blackwell, Laura Pasquier, Simon Dupont, Séverine Devers, Charlotte Lécureuil, *Joël Meunier <p style="text-align: justify;">In many animals, male weapons are large and extravagant morphological structures that typically enhance fighting ability and reproductive success. It is generally assumed that growing and carrying large weapons is c...Behavior, Evolution, Insecta, Invertebrates, Life histories, MorphologyOlivier Roux2024-03-26 08:56:27 View
26 Apr 2023
article picture

Microbiome mediated tolerance to biotic stressors: a case study of the interaction between a toxic cyanobacterium and an oomycete-like infection in Daphnia magna

Multi-stress responses depend on the microbiome in the planktonic crustacean Daphnia

Recommended by and ORCID_LOGO based on reviews by Natacha Kremer and 2 anonymous reviewers

The critical role that gut microbiota play in many aspects of an animal’s life, including pathogen resistance, detoxification, digestion, and nutritional physiology, is becoming more and more apparent (Engel and Moran 2013; Lindsay et al., 2020). Gut microbiota recruitment and maintenance can be largely affected by the surrounding environment (Chandler et al., 2011; Callens et al., 2020). The environment may thus dictate gut microbiota composition and diversity, which in turn can affect organismal responses to stress. Only few studies have, however, taken the gut microbiota into account to estimate life histories in response to multiple stressors in aquatic systems (Macke et al., 2016). 

Houwenhuyse et al., investigate how the microbiome affects life histories in response to ecologically relevant single and multiple biotic stressors (an oomycete-like parasite, and a toxic cyanobacterium) in Daphnia magna (Houwenhuyse et al., 2023). Daphnia is an excellent model, because this aquatic system lends itself extremely well for gut microbiota transplantation and manipulation. This is due to the possibility to sterilize eggs (making them free of bacteria), horizontal transmission of bacteria from the environment, and the relative ease of culturing genetically similar Daphnia clones in large numbers. 

The authors use an elegant experimental design to show that the Daphnia gut microbial community differs when derived from a laboratory versus natural inoculum, the latter being more diverse. The authors subsequently show that key life history traits (survival, fecundity, and body size) depend on the stressors (and combination thereof), the microbiota (structure and diversity), and Daphnia genotype. A key finding is that Daphnia exposed to both biotic stressors show an antagonistic interaction effect on survival (being higher), but only in individuals containing laboratory gut microbiota. The exact mechanism remains to be determined, but the authors propose several interesting hypotheses as to why Daphnia with more diverse gut microbiota do less well. This could be due, for example, to increased inter-microbe competition or an increased chance of contracting opportunistic, parasitic bacteria. For Daphnia with less diverse laboratory gut microbiota, a monopolizing species may be particularly beneficial for stress tolerance. Alongside these interesting findings, the paper also provides extensive information about the gut microbiota composition (available in the supplementary files), which is a very useful resource for other researchers. 

Overall, this study reveals that multiple, interacting factors affect the performance of Daphnia under stressful conditions. Of importance is that laboratory studies may be based on simpler microbiota systems, meaning that stress responses measured in the laboratory may not accurately reflect what is happening in nature. 

REFERENCES

Callens M, De Meester L, Muylaert K, Mukherjee S, Decaestecker E. The bacterioplankton community composition and a host genotype dependent occurrence of taxa shape the Daphnia magna gut bacterial community. FEMS Microbiology Ecology. 2020;96(8):fiaa128. https://doi.org/10.1093/femsec/fiaa128

Chandler JA, Lang JM, Bhatnagar S, Eisen JA, Kopp A. Bacterial communities of diverse Drosophila species: ecological context of a host-microbe model system. PLOS Genetics. 2011;7(9):e1002272. https://doi.org/10.1371/journal.pgen.1002272

Engel P, Moran NA. The gut microbiota of insects - diversity in structure and function. FEMS Microbiology Reviews. 2013;37(5):699-735. https://doi.org/10.1111/1574-6976.12025

Houwenhuyse S, Bulteel L, Vanoverberghe I, Krzynowek A, Goel N et al. Microbiome mediated tolerance to biotic stressors: a case study of the interaction between a toxic cyanobacterium and an oomycete-like infection in Daphnia magna. 2023. OSF, ver. 2 peer-reviewed and recommended by Peer Community in Zoology. https://doi.org/10.31219/osf.io/9n4mg

Lindsay EC, Metcalfe NB, Llewellyn MS. The potential role of the gut microbiota in shaping host energetics and metabolic rate. Journal of Animal Ecology. 2020;89(11):2415-2426. https://doi.org/10.1111/1365-2656.13327

Macke E, Tasiemski A, Massol F, Callens M, Decaestecker E. Life history and eco-evolutionary dynamics in light of the gut microbiota. Oikos. 2017;126(4):508-531. https://doi.org/10.1111/oik.03900

Microbiome mediated tolerance to biotic stressors: a case study of the interaction between a toxic cyanobacterium and an oomycete-like infection in *Daphnia magna*Shira Houwenhuyse*, Lore Bulteel*, Isabel Vanoverberghe, Anna Krzynowek, Naina Goel, Manon Coone, Silke Van den Wyngaert, Arne Sinnesael, Robby Stoks & Ellen Decaestecker<p style="text-align: justify;">Organisms are increasingly facing multiple, potentially interacting stressors in natural populations. The ability of populations coping with combined stressors depends on their tolerance to individual stressors and ...Aquatic, Biology, Crustacea, Ecology, Life histories, SymbiosisBertanne Visser2021-05-17 16:18:18 View
28 Apr 2021
article picture

Inference of the worldwide invasion routes of the pinewood nematode Bursaphelenchus xylophilus using approximate Bayesian computation analysis

Extracting the maximum historical information on pine wood nematode worldwide invasion from genetic data

Recommended by ORCID_LOGO based on reviews by Aude Gilabert and 1 anonymous reviewer

Redistribution of domesticated and non domesticated species by humans profoundly affected earth biogeography and in return human activities. This process accelerated exponentially since human expansion out of Africa, leading to the modern global, highly connected and homogenized, agriculture and trade system (Mack et al. 2000, Jaksic and Castro 2021), that threatens biological diversity and genetic resources. To accompany quarantine and control effort, the reconstruction of invasion routes provides valuable information that help identifying critical nodes and edges in the global networks (Estoup and Guillemaud 2010, Cristescu 2015). Historical records and genetic markers are the two major sources of information of this corpus of knowledge on Anthropocene historical phylogeography. With the advances of molecular genetics tools, the genealogy of these introductions events could be revisited and empowered. Due to their idiosyncrasy and intimate association with the contingency of human trades and activities, understanding the invasion and domestication routes require particular statistical tools (Fraimout et al. 2017).

Because it encompasses all these theoretical, ecological and economical implications, I am pleased to recommend the readers of PCI Zoology this article by Mallez et al. (2021) on pine wood nematode invasion route inference from genetic markers using Approximate Bayesian Computation (ABC) methods.

Economically and ecologically, this pest, is responsible for killing millions of pines worldwide each year. The results show these damages and the global genetic patterns are due to few events of successful introductions. The authors consider that this low probability of introductions success reinforces the idea that quarantine measures are efficient. This is illustrated in Europe where the pine-worm has been quarantined successfully in the Iberian Peninsula since 1999. Another relevant conclusion is that hybridization between invasive populations have not been observed and implied in the invasion process. Finally the present study reinforced the role of Asiatic bridgehead populations in invasion process including in Europe.

Methodologically, for the first time, ABC was applied to this species. A total of 310 individual sequences were added to the Mallez et al. (2015) microsatellite dataset. Fraimoult et al. (2017) showed the interest to apply random forest to improve scenario selection in ABC framework. This method, implemented in the DiYABC software (Collin et al. 2020) for invasion route scenario selection allows to handle more complex scenario alternatives and was used in this study. In this article by Mallez et al. (2021), you will also find a clear illustration of the step-by-step approach to select scenario using ABC techniques (Lombaert et al. 2014). The rationale is to reduce number of scenario to be tested by assuming that most recent invasions cannot be the source of the most ancient invasions and to use posterior results on most ancient routes as prior hypothesis to distinguish following invasions. The other simplification is to perform classical population genetic analysis to characterize genetic units and representative populations prior to invasion routes scenarios selection by ABC.

Yet, even when using the most advanced Bayesian inference methods, it is recognized by the authors that the method can be pushed to its statistical power limits. The method is appropriate when population show strong inter-population genetic structure. But the high number of differentiated populations in native area can be problematic since it is generally associated to incomplete sampling scheme. The hypothesis of ghost populations source allowed to bypass this difficulty, but the authors consider simulation studies are needed to assess the joint effect of genetic diversity and number of genetic markers on the inference results in such situation. Also the need to use a stepwise approach to reduce the number of scenario to test has to be considered with caution. Scenarios that are not selected but have non negligible posterior, cannot be ruled out in the constitution of next step scenarios hypotheses.

Due to its interest to understand this major facet of Anthropocene, reconstruction of invasion routes should be more considered as a guide to damper biological homogenization process.

References

Collin, F.-D., Durif, G., Raynal, L., Lombaert, E., Gautier, M., Vitalis, R., Marin, J.-M. and Estoup, A. (2020) Extending Approximate Bayesian Computation with Supervised Machine Learning to infer demographic history from genetic polymorphisms using DIYABC Random Forest. Authorea. doi: https://doi.org/10.22541/au.159480722.26357192

Cristescu, M.E. (2015) Genetic reconstructions of invasion history. Molecular Ecology, 24, 2212–2225. doi: https://doi.org/10.1111/mec.13117

Estoup, A. and Guillemaud, T., (2010) Reconstructing routes of invasion using genetic data: Why, how and so what? Molecular Ecology, 9, 4113-4130. doi: https://doi.org/10.1111/j.1365-294X.2010.04773.x

Fraimout, A., Debat, V., Fellous, S., Hufbauer, R.A., Foucaud, J., Pudlo, P., Marin, J.M., Price, D.K., Cattel, J., Chen, X., Deprá, M., Duyck, P.F., Guedot, C., Kenis, M., Kimura, M.T., Loeb, G., Loiseau, A., Martinez-Sañudo, I., Pascual, M., Richmond, M.P., Shearer, P., Singh, N., Tamura, K., Xuéreb, A., Zhang, J., Estoup, A. and Nielsen, R. (2017) Deciphering the routes of invasion of Drosophila suzukii by Means of ABC Random Forest. Molecular Biology and Evolution, 34, 980-996. doi: https://doi.org/10.1093/molbev/msx050

Jaksic, F.M. and Castro, S.A. (2021). Biological Invasions in the Anthropocene, in: Jaksic, F.M., Castro, S.A. (Eds.), Biological Invasions in the South American Anthropocene: Global Causes and Local Impacts. Springer International Publishing, Cham, pp. 19-47. doi: https://doi.org/10.1007/978-3-030-56379-0_2

Lombaert, E., Guillemaud, T., Lundgren, J., Koch, R., Facon, B., Grez, A., Loomans, A., Malausa, T., Nedved, O., Rhule, E., Staverlokk, A., Steenberg, T. and Estoup, A. (2014) Complementarity of statistical treatments to reconstruct worldwide routes of invasion: The case of the Asian ladybird Harmonia axyridis. Molecular Ecology, 23, 5979-5997. doi: https://doi.org/10.1111/mec.12989

Mack, R.N., Simberloff, D., Lonsdale, M.W., Evans, H., Clout, M., Bazzaz, F.A. (2000) Biotic Invasions : Causes , Epidemiology , Global Consequences , and Control. Ecological Applications, 10, 689-710. doi: https://doi.org/10.1890/1051-0761(2000)010[0689:BICEGC]2.0.CO;2

Mallez, S., Castagnone, C., Lombaert, E., Castagnone-Sereno, P. and Guillemaud, T. (2021) Inference of the worldwide invasion routes of the pinewood nematode Bursaphelenchus xylophilus using approximate Bayesian computation analysis. bioRxiv, 452326, ver. 6 peer-reviewed and recommended by Peer community in Zoology. doi: https://doi.org/10.1101/452326

Inference of the worldwide invasion routes of the pinewood nematode Bursaphelenchus xylophilus using approximate Bayesian computation analysisSophie Mallez, Chantal Castagnone, Eric Lombaert, Philippe Castagnone-Sereno, Thomas Guillemaud<p>Population genetics have been greatly beneficial to improve knowledge about biological invasions. Model-based genetic inference methods, such as approximate Bayesian computation (ABC), have brought this improvement to a higher level and are now...Biogeography, Biological invasions, Ecology, Evolution, Genetics/Genomics, Herbivores, Invertebrates, Molecular biology, NematologyStéphane Dupas2020-09-15 10:59:41 View
21 Mar 2023
article picture

Population genetics of Glossina fuscipes fuscipes from southern Chad

Population genetics of tsetse, the vector of African Trypanosomiasis, helps informing strategies for control programs

Recommended by ORCID_LOGO based on reviews by 2 anonymous reviewers

Human African Trypanosomiasis (HAT), or sleeping sickness, is caused by trypanosome parasites. In sub-Saharan Africa, two forms are present, Trypanosoma brucei gambiense and T. b. rhodesiense, the former responsible for 95% of reported cases. The parasites are transmitted through a vector, Genus Glossina, also known as tsetse, which means fly in Tswana, a language from southern Africa. Through a blood meal, tsetse picks up the parasite from infected humans or animals (in animals, the parasite causes Animal African Trypanosomiasis or nagana disease). Through medical interventions and vector control programs, the burden of the disease has drastically reduced over the past two decades, so the WHO neglected tropical diseases road map targets the interruption of transmission (zero cases) for 2030 (WHO 2022).

Meaningful vector control programs utilize traps for the removal of animals and for surveillance, along with different methods of spraying insecticides. However, in existing HAT risk areas, it will be essential to understand the ecology of the vector species to implement control programs in a way that areas cleared from the vector will not be reinvaded from other populations. Thus, it will be crucial to understand basic population genetics parameters related to population structure and subdivision, migration frequency and distances, population sizes, and the potential for sex-biased dispersal. The authors utilize genotyping using nine highly polymorphic microsatellite markers of samples from Chad collected in differently affected regions and at different time points (Ravel et al., 2023). Two major HAT zones exist that are targeted by vector control programs, namely Madoul and Maro, while two other areas, Timbéri and Dokoutou, are free of trypanosomes. Samples were taken before vector control programs started.

The sex ratio was female-biased, most strongly in Mandoul and Maro, the zones with the lowest population density. This could be explained by resource limitation, which could be the hosts for a blood meal or the sites for larviposition. Limited resources mean that females must fly further, increasing the chance that more females are caught in traps. 

The effective population densities of Mandoul and Maro were low. However, there was a convergence of population density and trapping density, which might be explained by the higher preservation of flies in the high-density areas of Timbéri and Dokoutou after the first round of sampling, which can only be tested using a second sampling. 

The dispersal distances are the highest recorded so far, especially in Mandoul and Maro, with 20-30 km per generation. However, in Timbéri and Dokoutou, which are 50 km apart, very little exchange occurs (approx. 1-2 individuals every six months). A major contributor to this is the massive destruction of habitat that started in the early 1990s and left patchily distributed and fragmented habitats. The Mandoul zone might be safe from reinvasion after eradication, as for a successful re-establishment, either a pair of a female and male or a pregnant female are required. As the trypanosome prevalence amongst humans was 0.02 and of tsetse 0.06 (Ibrahim et al., 2021) before interventions began, medical interventions and vector control might have further reduced these levels, making a reinvasion and subsequent re-establishment of HAT very unlikely. Maro is close to the border of the Central African Republic, and the area has not been well investigated concerning refugee populations of tsetse, which could contribute to a reinvasion of the Maro zone. The higher level of genetic heterogeneity of the Maro population indicates that invasions from neighboring populations are already ongoing. This immigration could also be the reason for not detecting the bottleneck signature in the Maro population. 

The two HAT areas need different levels of attention while implementing vector eradication programs. While Madoul is relatively safe against reinvasion, Maro needs another type of attention, as frequent and persistent immigration might counteract eradication efforts. Thus, it is recommended that continuous tsetse suppression needs to be implemented in Maro.  

This study shows nicely that an in-depth knowledge of the processes within and between populations is needed to understand how these populations behave. This can be used to extrapolate, make predictions, and inform the organisations implementing vector control programs to include valuable adjustments, as in the case of Maro. Such integrative approaches can help prevent the failure of programs, potentially saving costs and preventing infections of humans and animals who might die if not treated.

References

Ibrahim MAM, Weber JS, Ngomtcho SCH, Signaboubo D, Berger P, Hassane HM, Kelm S (2021) Diversity of trypanosomes in humans and cattle in the HAT foci Mandoul and Maro, Southern Chad- Southern Chad-A matter of concern for zoonotic potential? PLoS Neglected Tropical Diseases, 15, e000 323. https://doi.org/10.1371/journal.pntd.0009323

Ravel S, Mahamat MH, Ségard A, Argiles-Herrero R, Bouyer J, Rayaisse JB, Solano P, Mollo BG, Pèka M, Darnas J, Belem AMG, Yoni W, Noûs C, de Meeûs T (2023) Population genetics of Glossina fuscipes fuscipes from southern Chad. Zenodo, ver. 9 peer-reviewed and recommended by PCI Zoology. https://doi.org/10.5281/zenodo.7763870

WHO (2022) Trypanosomiasis, human African (sleeping sickness). https://www.who.int/news-room/fact-sheets/detail/trypanosomiasis-human-african-(sleeping-sickness), retrieved 17. March 2023

Population genetics of Glossina fuscipes fuscipes from southern ChadSophie Ravel, Mahamat Hissène Mahamat, Adeline Ségard, Rafael Argiles-Herrero, Jérémy Bouyer, Jean-Baptiste Rayaisse, Philippe Solano, Brahim Guihini Mollo, Mallaye Pèka, Justin Darnas, Adrien Marie Gaston Belem, Wilfrid Yoni, Camille Noûs, Thierr...<p>In Subsaharan Africa, tsetse flies (genus Glossina) are vectors of trypanosomes causing Human African Trypanosomiasis (HAT) and Animal African Trypanosomosis (AAT). Some foci of HAT persist in Southern Chad, where a program of tsetse control wa...Biology, Ecology, Evolution, Genetics/Genomics, Insecta, Medical entomology, Parasitology, Pest management, Veterinary entomologyMichael Lattorff Audrey Bras2022-04-22 11:25:24 View
08 Mar 2024
article picture

A comparison of the parasitoid wasp species richness of tropical forest sites in Peru and Uganda – subfamily Rhyssinae (Hymenoptera: Ichneumonidae)

Two sides of tropical richness, parasitoid wasps collected by Malaise traps in tropical rainforests of South America and Africa

Recommended by based on reviews by Mabel Alvarado, Filippo Di Giovanni and 2 anonymous reviewers

Insect species richness and diversity comparisons between samples of the tropics around the world are rare, especially in taxa composed mainly of cryptic species as parasitoid wasps.

The article by Hopkins et al. (2024) compares samples of parasitoid wasps of the subfamily Rhyssinae (Hymenoptera: Ichneumonidae) collected by Malaise traps in tropical rainforests of Perú and Uganda. The samples presented several differences in the time of collecting, covertures, and the sampling number; however, they used the same kind of traps, and the taxonomic process for species delimitation was made for the same team of ichneumonid experts, using equivalent characters.

Publications about this kind of comparative study are difficult to find because cooperative projects on insect richness and diversity from South American and African continents are not frequent. In this sense, this study presented a valuable contrast that shows interesting results about the higher richness and lower abundance of the biota of the American tropics, even with a small sample, in comparison with the biota of the African tropics. The results are supported mainly by the rarefaction curves shown. This pattern of higher species richness and lower specimen abundance, observed in other American tropical taxa such as trees, birds, or butterflies, is observed too in these parasitoid wasps, increasing the body of information that could support the extension of the pattern to the entire biota of the American tropics. The authors recognize the study's limitations, which include strong differences in the size of the forest coverture between places. However, these differences and others are enough described and discussed.

This work is useful because it increases the information about the diversity patterns of the tropics around the world and because study a taxon mainly composed of cryptic species, with a small amount of information in tropical regions.

References

Hopkins T., Tuomisto H., Gómez I.C., Sääksjärvi I. E. 2024. A comparison of the parasitoid wasp species richness of tropical forest sites in Peru and Uganda – subfamily Rhyssinae (Hymenoptera: Ichneumonidae). bioRxiv, ver. 2 peer-reviewed and recommended by Peer Community in Zoology. https://doi.org/10.1101/2023.08.23.554460

A comparison of the parasitoid wasp species richness of tropical forest sites in Peru and Uganda – subfamily Rhyssinae (Hymenoptera: Ichneumonidae)Tapani Hopkins, Hanna Tuomisto, Isrrael C. Gómez, Ilari E. Sääksjärvi<p style="text-align: justify;">The global distribution of parasitoid wasp species richness is poorly known. Past attempts to compare data from different sites have been hampered by small sample sizes and lack of standardisation. During the past d...Biodiversity, Biogeography, InsectaGiovanny Fagua2023-08-24 18:30:26 View
28 Aug 2022
article picture

A simple procedure to detect, test for the presence of stuttering, and cure stuttered data with spreadsheet programs

Improved population genetics parameters through control for microsatellite stuttering

Recommended by ORCID_LOGO based on reviews by Thibaut Malausa, Fabien Halkett and Thierry Rigaud

Molecular markers have drastically changed and improved our understanding of biological processes. In combination with PCR, markers revolutionized the study of all organisms, even tiny insects, and eukaryotic pathogens amongst others. Microsatellite markers were the most prominent and successful ones. Their success started in the early 1990s. They were used for population genetic studies, mapping of genes and genomes, and paternity testing and inference of relatedness. Their popularity is based on some of their characteristics as codominance, the high polymorphism information content, and their ease of isolation (Schlötterer 2004). Still, microsatellites are the marker of choice for a range of non-model organisms as next-generation sequencing technologies produce a huge amount of single nucleotide polymorphisms (SNPs), but often at expense of sample size and higher costs.
 
The high level of polymorphism of microsatellite markers, which consist of one to six base-pair nucleotide motifs replicated up to 10 or 20 times, results from slippage events during DNA replication. Short hairpin loops might shorten the template strand or extend the new strand. However, such slippage events might occur during PCR amplification resulting in additional bands or peaks. Such stutter alleles often appear to differ by one repeat unit and might be hard to interpret but definitively reduce automated scoring of microsatellite results.
 
A standalone software package available to handle stuttering is Microchecker (van Oosterhout et al., 2004, which nowadays faces incompatibilities with updated versions of different operating systems. Thus, de Meeûs and Noûs (2022), in their manuscript, tackled the stuttering issue by developing an OS-independent analysis pipeline based on standard spreadsheet software such as Microsoft Office (Excel) or Apache Open Office (Calc). The authors use simulated populations differing in the mating system (pangamic, selfing (30%), clonal) and a different number of subpopulations and individuals per subpopulation to test for differences among the null model (no stuttering), a test population with 2 out of 20 loci (10%) with stuttering, and the latter with stuttering cured. Further to this, the authors also re-analyse data from previous studies utilising organisms differing in the mating system to understand whether control of stuttering changes major parameter estimates and conclusions of those studies.
 
Stuttering of microsatellite loci might result in increased heterozygote deficits. The authors utilise the FIS (inbreeding coefficient) as a tool to compare the different treatments of the simulated populations. Their method detected stuttering in pangamic and selfing populations, while the detection of stuttering in clonal organisms is more difficult. The cure for stuttering resulted in FIS values similar to those populations lacking stuttering. The re-analysis of four previously published studies indicated that the new method presented here is more accurate than Microchecker (van Oosterhout et al., 2004) in a direct comparison. For the Lyme disease-transmitting tick Ixodes scapularis (De Meeûs et al., 2021), three loci showed stuttering and curing these resulted in data that are in good agreement with pangamic reproduction. In the tsetse fly Glossina palpalis palpalis (Berté et al., 2019), two out of seven loci were detected as stuttering. Curing them resulted in decreased FIS for one locus, while the other showed an increased FIS, an indication of other problems such as the occurrence of null alleles. Overall, in dioecious pangamic populations, the method works well, and the cure of stuttering improves population genetic parameter estimates, although FST and FIS might be slightly overestimated. In monoecious selfers, the detection and cure work well, if other factors such as null alleles do not interfere. In clonal organisms, only loci with extremely high FIS might need a cure to improve parameter estimates.
 
This spreadsheet-based method helps to automate microsatellite analysis at very low costs and thus improves the accuracy of parameter estimates. This might certainly be very useful for a range of non-model organisms, parasites, and their vectors, for which microsatellites are still the marker of choice. 
 
References

Berté D, De Meeus T, Kaba D, Séré M, Djohan V, Courtin F, N'Djetchi KM, Koffi M, Jamonneau V, Ta BTD, Solano P, N’Goran EK, Ravel S (2019) Population genetics of Glossina palpalis palpalis in sleeping sickness foci of Côte d'Ivoire before and after vector control. Infection Genetics and Evolution 75, 103963. https://doi.org/0.1016/j.meegid.2019.103963

de Meeûs T, Chan CT, Ludwig JM, Tsao JI, Patel J, Bhagatwala J, Beati L (2021) Deceptive combined effects of short allele dominance and stuttering: an example with Ixodes scapularis, the main vector of Lyme disease in the U.S.A. Peer Community Journal 1, e40. https://doi.org/10.24072/pcjournal.34

de Meeûs T, Noûs C (2022) A simple procedure to detect, test for the presence of stuttering, and cure stuttered data with spreadsheet programs. Zenodo, v5, peer-reviewed and recommended by PCI Zoology. https://doi.org/10.5281/zenodo.7029324

Schlötterer C (2004) The evolution of molecular markers - just a matter of fashion? Nature Reviews Genetics 5, 63-69. https://doi.org/10.1038/nrg1249

van Oosterhout C, Hutchinson WF, Wills DPM, Shipley P (2004) MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Molecular Ecology Notes 4, 535-538. https://doi.org/10.1111/j.1471-8286.2004.00684.x

A simple procedure to detect, test for the presence of stuttering, and cure stuttered data with spreadsheet programsThierry de Meeûs and Camille Noûs<p>Microsatellite are powerful markers for empirical population genetics, but may be affected by amplification problems like stuttering that produces heterozygote deficits between alleles with one repeat difference. In this paper, we present a sim...Acari, Ecology, Evolution, Genetics/Genomics, Helminthology, Invertebrates, Medical entomology, Molecular biology, Parasitology, Theoretical biology, Veterinary entomologyMichael Lattorff2021-12-06 14:30:47 View