Phylogenetic trees can be used to extract information about the process of diversification that has generated them. The most common approach to conduct this inference is to rely on a likelihood, defined here as the probability of generating a dated tree T given a diversification model (e.g. a birth-death model), and then use standard maximum likelihood. This idea has been explored extensively in the context of the so-called diversification studies, with many variants for the models and for the questions being asked (diversification rates shifting at certain time points or in the ancestors of particular subclades, trait-dependent diversification rates, etc).
However, all this assumes that the dated tree T is known without error. In practice, trees (that is, both the tree topology and the divergence times) are inferred based on DNA sequences, possibly combined with fossil information for calibrating and informing the divergence times. Molecular dating is a delicate exercise, however, and much more so in fact than reconstructing the tree topology. In particular, a mis-specificied model for the relaxed molecular clock, or a mis-specifiied prior, can have a substantial impact on the estimation of divergence dates - which in turn could severely mislead the inference about the underlying diversification process. This thus raises the following question: would that be possible to conduct inference and testing of diversification models without having to go through the dangerous step of molecular dating?
In his article ""Probabilities of tree topologies with temporal constraints and diversification shifts"" [1], Gilles Didier introduces a recursive method for computing the probability of a tree topology under some diversification model of interest, without knowledge of the exact dates, but only interval constraints on the dates of some of the nodes of the tree. Such interval constraints, which are derived from fossil knowledge, are typically used for molecular dating: they provide the calibrations for the relaxed clock analysis. Thus, what is essentially proposed by Gilles Didier is to use them in combination with the tree topology only, thus bypassing the need to estimates divergence times first, before fitting a diversification model to a phylogenetic tree.
This article, which is primarily a mathematical and algorithmic contribution, is then complemented with several applications: testing for a diversification shift in a given subclade of the phylogeny, just based on the (undated) tree topology, with interval constraints on some of its internal nodes; but also, computing the age distribution of each node and sampling on the joint distribution on node ages, conditional on the interval constraints. The test for the presence of a diversification shift is particularly interesting: an application to simulated data (and without any interval constraint in that case) suggests that the method based on the undated tree performs about as well as the classical method based on a dated tree, and this, even granting the classical approach a perfect knowledge of the dates - given that, in practice, one in fact relies on potentially biased estimates. Finally, an application to a well-known example (rate shifts in cetacean phylogeny) is presented.
This article thus represents a particularly meaningful contribution to the methodology for diversification studies; but also, for molecular dating itself: it is a well known problem in molecular dating that computing and sampling from the conditional distributions on node ages, given fossil constraints, and more generally understanding and visualizing how interval constraints on some nodes of the tree impact the distribution at other nodes, is a particularly difficult exercise. For that reason, the algorithmic routines presented in the present article will be useful in this context as well.

References

[1] Didier, G. (2020) Probabilities of tree topologies with temporal constraints and diversification shifts. bioRxiv, 376756, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/376756

Jean Cury and colleagues (2021) have developed a protocol to simulate bacterial evolution in SLiM. In contrast to existing methods that depend on the coalescent, SLiM simulates evolution forward in time. SLiM has, up to now, mostly been used to simulate the evolution of eukaryotes (Haller and Messer 2019), but has been adapted here to simulate evolution in bacteria. Forward-in-time simulations are usually computationally very costly. To circumvent this issue, bacterial population sizes are scaled down. One would now expect results to become inaccurate, however, Cury et al. show that scaled-down forwards simulations provide very accurate results (similar to those provided by coalescent simulators) that are consistent with theoretical expectations. Simulations were analyzed and compared to existing methods in simple and slightly more complex scenarios where recombination affects evolution. In all scenarios, simulation results from coalescent methods (fastSimBac (De Maio and Wilson 2017), ms (Hudson 2002)) and scaled-down forwards simulations were very similar, which is very good news indeed.

A biologist not aware of the complexities of forwards, backwards simulations and the coalescent, might now naïvely ask why another simulation method is needed if existing methods perform just as well. To address this question the manuscript closes with a very neat example of what exactly is possible with forwards simulations that cannot be achieved using existing methods. The situation modeled is the growth and evolution of a set of 50 bacteria that are randomly distributed on a petri dish. One side of the petri dish is covered in an antibiotic the other is antibiotic-free. Over time, the bacteria grow and acquire antibiotic resistance mutations until the entire artificial petri dish is covered with a bacterial lawn. This simulation demonstrates that it is possible to simulate extremely complex (e.g. real world) scenarios to, for example, assess whether certain phenomena are expected with our current understanding of bacterial evolution, or whether there are additional forces that need to be taken into account. Hence, forwards simulators could significantly help us to understand what current models can and cannot explain in evolutionary biology.  

References  

Cury J, Haller BC, Achaz G, Jay F (2021) Simulation of bacterial populations with SLiM. bioRxiv, 2020.09.28.316869, version 5 peer-reviewed and recommended by Peer community in Evolutionary Biology.  https://doi.org/10.1101/2020.09.28.316869

De Maio N, Wilson DJ (2017) The Bacterial Sequential Markov Coalescent. Genetics, 206, 333–343. https://doi.org/10.1534/genetics.116.198796

Haller BC, Messer PW (2019) SLiM 3: Forward Genetic Simulations Beyond the Wright–Fisher Model. Molecular Biology and Evolution, 36, 632–637. https://doi.org/10.1093/molbev/msy228

Hudson RR (2002) Generating samples under a Wright–Fisher neutral model of genetic variation. Bioinformatics, 18, 337–338. https://doi.org/10.1093/bioinformatics/18.2.337

An individual’s sexual phenotype is determined during development. Understanding which pathways are activated or repressed during the developmental stages leading to a sexually mature individual, for example by studying gene expression and how its level is biased between sexes, allows us to understand the functional aspects of dimorphic phenotypes between the sexes.

Several studies have quantified the differences in transcription between the sexes in mature individuals, showing the extent of this sex-bias and which functions are affected. There is, however, less data available on what occurs during the different phases of development leading to this phenotype, especially in species with specific developmental strategies, such as hemimetabolous insects. While many well-studied insects such as the honey bee, drosophila, and butterflies, exhibit an holometabolous development ("holo" meaning "complete" in reference to their drastic metamorphosis from the juvenile to the adult stage), hemimetabolous insects have juvenile stages that look similar to the adult stage (the hemi prefix meaning "half", referring to the more tissue-specific changes during development), as seen in crickets, cockroaches, and stick insects. Learning more about what happens during development in terms of the identity of genes that are sex-biased (are they the same genes at different developmental stages? What are their function? Do they exhibit specific sequence evolution rates? Is one sex over-represented in the sex-biased genes?) and their quantity over developmental time (gradual or abrupt increase in number, if any?) would allow us to better understand the evolution of sexual dimorphism at the gene expression level and how it relates to dimorphism at the organismic level.

Djordjevic et al (2021) studied the transcriptome during development in an hemimetabolous stick insect, to improve our knowledge of this type of development, where the organismic phenotype is already mostly present in the early life stages. To do this, they quantified whole-genome gene expression levels in whole insects, using RNA-seq at three different developmental stages. One of the interesting results presented by Djordjevic and colleagues is that the increase in the number of genes that were sex-biased in expression is gradual over the three stages of development studied and it is mostly the same genes that stay sex-biased over time, reflecting the gradual change in phenotypes between hatchlings, juveniles and adults. Furthermore, male-biased genes had faster sequence divergence rates than unbiased genes and that female-biased genes.

This new information of sex-bias in gene expression in an hemimetabolous insect allowed the authors to do a comparison of sex-biased genes with what has been found in a well-studied holometabolous insect, Drosophila. The gene expression patterns showed that four times more genes were sex-biased in expression in that species than in stick insects. Furthermore, the increase in the number of sex-biased genes during development was quite abrupt and clearly distinct in the adult stage, a pattern that was not seen in stick insects. As pointed out by the authors, this pattern of a "burst" of sex-biased genes at maturity is more common than the gradual increase seen in stick insects.

With this study, we now know more about the evolution of sex-biased gene expression in an hemimetabolous insect and how it relates to their phenotypic dimorphism. Clearly, the next step will be to sample more hemimetabolous species at different life stages, to see how this pattern is widespread or not in this mode of development in insects.

References

Djordjevic J, Dumas Z, Robinson-Rechavi M, Schwander T, Parker DJ (2021) Dynamics of sex-biased gene expression during development in the stick insect Timema californicum. bioRxiv, 2021.01.23.427895, ver. 6 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.01.23.427895

Exciting times are afoot for those of us interested in the ecology and evolution of invasive populations. Recent years have seen evolutionary process woven firmly into our understanding of invasions (Miller et al. 2020). This integration has inspired a welter of empirical and theoretical work. We have moved from field observations and verbal models to replicate experiments and sophisticated mathematical models. Progress has been rapid, and we have seen science at its best; an intimate discussion between theory and data.
An area currently under very active development is our understanding of pushed invasions. Here a population spreads through space driven, not by dispersal and growth originating at the leading tip of the invasion, but by dispersal and growth originating deeper in the bulk of the population. These pushed invasions may be quite common – they result when per capita growth and dispersal rates are higher in the bulk of the wave than at the leading tip. They result from a range of well-known phenomena, including Allee effects and density-dependent dispersal (Gandhi et al. 2016; Bîrzu et al. 2019). Pushed invasions travel faster than we would expect given growth and dispersal rates on the leading tip, and they lose genetic diversity more slowly than classical pulled invasions (Roques et al. 2012; Haond et al. 2018; Bîrzu et al. 2019).
Well… in theory, anyway. The theory on pushed waves has momentarily streaked ahead of the empirical work, because empirical systems for studying pushed invasions are rare (though see Gandhi et al. 2016; Gandhi, Korolev, and Gore 2019). In this paper, Dahirel and colleagues (2020) make the argument that we may be able to generate pushed invasions in laboratory systems simply by reducing the connectedness of our experimental landscapes. If true, we might have a simple tool for turning many of our established experimental systems into systems for studying pushed dynamics.
It’s a nice idea, and the paper goes to careful lengths to explore the possibility in their lab system (a parasitoid wasp, Trichogramma). They run experiments on replicate wasp populations comparing strongly- v poorly-connected arrays, and estimate the resulting invasion speeds and rate of diversity loss. They also build a simulation model of the system, allowing them to explore in-silico a range of possible processes underlying their results.
As well as developing these parallel systems, Dahirel and colleagues (2020) go to careful lengths to develop statistical analyses that allow inference on key parameters, and they apply these analyses to both the experimental and simulation data. They have been motivated to apply methods that might be used in both laboratory and field settings to help classify invasions.
Ultimately, they found reasonable evidence that their poorly-connected habitat did induce a pushed dynamic. Their poorly connected invasions travelled faster than they should have if they were pulled, they lost diversity more slowly than the highly connected habitat, and replicates with a higher carrying capacity tended to have higher invasion speeds. All in line with expectations of a pushed dynamic. Interestingly, however, their simulation results suggest that they probably got this perfect result for unexpected reasons. The strong hint is that their poorly-connected habitat induced density dependent dispersal in the wasps. Without this effect, their simulations suggest they should have seen diversity decreasing much more rapidly than it did.
There is a nuanced, thoughtful, and carefully argued discussion about all this in the paper, and it is worth reading. There is much of value in this paper. Theirs is not a perfect empirical system in which all the model assumptions are met and in which huge population sizes make stochastic effects negligible. Here is a system one step closer to the messy reality of biology. The struggle to align this system with new theory has been worth the effort. Not only does it give us hope that we might usefully be able to discriminate between classes of invasions using real-world data, but it hints at a rule that Tolstoy might have expressed this way: all pulled invasions are alike, each pushed invasion is pushed in its own way.

References

Bîrzu, G., Matin, S., Hallatschek, O., and Korolev, K. S. (2019). Genetic drift in range expansions is very sensitive to density dependence in dispersal and growth. Ecology Letters, 22(11), 1817-1827. doi: https://doi.org/10.1111/ele.13364
Dahirel, M., Bertin, A., Haond, M., Blin, A., Lombaert, E., Calcagno, V., Fellous, S., Mailleret, L., Malausa, T., and Vercken, E. (2020). Shifts from pulled to pushed range expansions caused by reduction of landscape connectivity. bioRxiv, 2020.05.13.092775, ver. 4 peer-reviewed and recommended by PCI Evolutionary Biology. https://doi.org/10.1101/2020.05.13.092775
Gandhi, S. R., Korolev, K. S., and Gore, J. (2019). Cooperation mitigates diversity loss in a spatially expanding microbial population. Proceedings of the National Academy of Sciences, 116(47), 23582-23587. doi: https://doi.org/10.1073/pnas.1910075116
Gandhi, S. R., Yurtsev, E. A., Korolev, K. S., and Gore, J. (2016). Range expansions transition from pulled to pushed waves as growth becomes more cooperative in an experimental microbial population. Proceedings of the National Academy of Sciences, 113(25), 6922-6927. doi: https://doi.org/10.1073/pnas.1521056113
Haond, M., Morel-Journel, T., Lombaert, E., Vercken, E., Mailleret, L. and Roques, L. (2018). When higher carrying capacities lead to faster propagation (2018), bioRxiv, 307322, ver. 4 peer-reviewed and recommended by Peer Community in Ecology. https://doi.org/10.1101/307322
Miller et al. (2020). Eco‐evolutionary dynamics of range expansion. Ecology, 101(10), e03139. doi: https://doi.org/10.1002/ecy.3139
Roques, L., Garnier, J., Hamel, F., and Klein, E. K. (2012). Allee effect promotes diversity in traveling waves of colonization. Proceedings of the National Academy of Sciences, 109(23), 8828-8833. doi: https://doi.org/10.1073/pnas.1201695109

​​Durability of plant resistance to diploid pathogen Hirohisa Kishino Based on the population genetic and epidemiologic model, Saubin et al. (2021) report that the resistant hosts generated by the breeding based on the gene-for-gene interaction is durable much longer against diploid pathogens than haploid pathogens. The avr allele of pathogen that confers the resistance is genetically recessive. The heterozygotes are not recognized by the resistant hosts and only the avr/avr homozygote is adaptive. As a result, the trajectory of avr allele frequency becomes more stochastic due to genetic drift. Although the paper focuses on the evolution of standing polymorphism, it seems obvious that the adaptive mutations in pathogen have much larger probability of being deleted from the population because the individuals own the avr allele mostly in the form of heterozygote at the initial phase after the mutation. Since only few among many models of plant resistance deployment study the case of diploid pathogen and the contribution of the pathogen life cycle, this work will add an important intellect to the literature (Rimbaud et al. 2021).

From the study of host-parasite interaction in flax rust Melampsora lini, Flor (1942, 1955) showed that the host resistance is formed by the interaction of a host resistance gene and a corresponding pathogen gene. This gene-for-gene hypothesis has been supported by experimental evidence and has served as a basis of the methods of molecular breeding targeting the dominant R genes. However, modern agriculture provides the pathogen populations with the homogeneous environments and laid strong selection pressure on them. As a result, the newly developed resistant plants face the risk of immediate resistance breakdown (Möller and Stukenbrock 2017).

Currently, quantitative resistance is getting attention as characters as a potential target for long-life (mild) resistant breeds (Lannou, 2012). They are polygenic and controlled partly by the same genes that mediate qualitative resistance but mostly by the genes that encode defense-related outputs such as strengthening of the cell wall or defense compound biosynthesis (Corwin and Kliebenstein, 2017). Progress of molecular genetics may overcome the technical difficulty (Bakkeren and Szabo, 2020). Saubin et al. (2021) notes that the pattern of genetic inheritance of the pathogen counterparts that respond to the host traits is crucial regarding with the durability of the resistant hosts. The resistance traits for which avr alleles are predicted to be recessive may be the targets of breeding.

References

Bakkeren, G., and Szabo, L. J. (2020) Progress on molecular genetics and manipulation of rust fungi. Phytopathology, 110, 532-543. https://doi.org/10.1094/PHYTO-07-19-0228-IA

Corwin, J. A., and Kliebenstein, D. J. (2017) Quantitative resistance: more than just perception of a pathogen. The Plant Cell, 29, 655-665. https://doi.org/10.1105/tpc.16.00915

Flor, H. H. (1942) Inheritance of pathogenicity in a cross between physiological races 22 and 24 of Melampsova lini. Phytopathology, 35. Abstract.

Flor, H. H. (1955) Host-parasite interactions in flax rust-its genetics and other implications. Phytopathology, 45, 680-685.

Lannou, C. (2012) Variation and selection of quantitative traits in plant pathogens. Annual review of phytopathology, 50, 319-338. https://doi.org/10.1146/annurev-phyto-081211-173031

Möller, M. and Stukenbrock, E. H. (2017) Evolution and genome architecture in fungal plant pathogens. Nature Reviews Microbiology. 15, 756–771. https://doi.org/10.1038/nrmicro.2017.76

Rimbaud, L., Fabre, F., Papaïx, J., Moury, B., Lannou, C., Barrett, L. G., and Thrall, P. H. (2021) Models of Plant Resistance Deployment. Annual Review of Phytopathology, 59. https://doi.org/10.1146/annurev-phyto-020620-122134

Saubin, M., De Mita, S., Zhu, X., Sudret, B. and Halkett, F. (2021) Impact of ploidy and pathogen life cycle on resistance durability. bioRxiv, 2021.05.28.446112, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.05.28.446112