Genome-wide data obtained from even a small number of individuals can provide unprecedented levels of detail about the evolutionary history of populations and species [1], determinants of genetic diversity [2], species boundaries and the process of speciation itself [3]. Loire and Galtier [4] present a clear example, using the emblematic Galápagos giant tortoise (Chelonoidis nigra), of how multi-species comparative population genomic approaches can provide valuable insights about population structure and species delimitation even when sample sizes are limited but the number of loci is large and distributed across the genome.

Galápagos giant tortoises are endemic to the Galápagos Islands and are currently recognized as an endangered, multi-species complex including both extant and extinct taxa. Taxonomic definitions are based on morphology, geographic isolation and population genetic evidence based on short DNA sequences of the mitochondrial genome (mtDNA) and/or a dozen or so nuclear microsatellite loci [5-8]. The species complex enjoys maximal protection. Population recoveries have been quite successful and spectacular conservation programs based on mitochondrial genes and microsatellites are ongoing. This includes for example individual translocations, breeding program, “hybrid” sterilization or removal, and resurrection of extinct lineages).

In 2013, Loire et al. [9] provided the first population genomic analyses based on genome scale data (~1000 coding loci derived from blood-transcriptomes) from five individuals, encompassing three putative “species”: Chelonnoidis becki, C. porteri and C. vandenburghi. Their results raised doubts about the validity/accuracy of the currently accepted designations of “genetic distinctiveness”. However, the implications for conservation and management have remained unnoticed.

In 2017, Loire and Galtier [4] have re-appraised this issue using an original multi-species comparative population genomic analysis of their previous data set [9]. Based on a comparison of 53 animal species, they show that both the level of genome-wide neutral diversity (πS) and level of population structure estimated using the inbreeding coefficient (F) are much lower than would be expected from a sample covering multiple species. The observed values are more comparable to those typically reported at the “among population” level within a single species such as human (Homo sapiens). The authors go to great length to assess the sensitivity of their method to detect population structure (or lack thereof) and show that their results are robust to potential issues, such as contamination and sequencing errors that can occur with Next Generation Sequencing techniques; and biases related to the small sample size and sub-sampling of individuals. They conclude that published mtDNA and microsatellite-based assessment of population structure and species designations may be biased towards over-splitting.

This manuscript is a very good read as it shows the potential of the now relatively affordable genome-wide data for helping to both resolve and clarify population and species boundaries, illuminate demographic trends, evolutionary trajectories of isolated groups, patterns of connectivity among them, and test for evidence of local adaptation and even reproductive isolation. The comprehensive information provided by genome-wide data can critically inform and assist the development of the best strategies to preserve endangered populations and species. Loire and Galtier [4] make a strong case for applying genomic data to the Galápagos giant tortoises, which is likely to redirect conservation efforts more effectively and at lower cost. The case of the Galápagos giant tortoises is certainly a very emblematic example, which will find an echo in many other endangered species conservation programs.

References

[1] Li H and Durbin R. 2011. Inference of human population history from individual whole-genome sequences. Nature, 475: 493–496. doi: 10.1038/nature10231

[2] Romiguier J, Gayral P, Ballenghien M, Bernard A, Cahais V, Chenuil A, Chiari Y, Dernat R, Duret L, Faivre N, Loire E, Lourenco JM, Nabholz B, Roux C, Tsagkogeorga G, Weber AA-T, Weinert LA, Belkhir K, Bierne N, Glémin S and Galtier N. 2014. Comparative population genomics in animals uncovers the determinants of genetic diversity. Nature, 515: 261–263. doi: 10.1038/nature13685

[3] Roux C, Fraïsse C, Romiguier J, Anciaux Y, Galtier N and Bierne N. 2016. Shedding light on the grey zone of speciation along a continuum of genomic divergence. PLoS Biology, 14: e2000234. doi: 10.1371/journal.pbio.2000234

[4] Loire E and Galtier N. 2017. Lacking conservation genomics in the giant Galápagos tortoise. bioRxiv 101980, ver. 4 of September 26, 2017. doi: 10.1101/101980

[5] Beheregaray LB, Ciofi C, Caccone A, Gibbs JP and Powell JR. 2003. Genetic divergence, phylogeography and conservation units of giant tortoises from Santa Cruz and Pinzón, Galápagos Islands. Conservation Genetics, 4: 31–46. doi: 10.1023/A:1021864214375

[6] Ciofi C, Milinkovitch MC, Gibbs JP, Caccone A and Powell JR. 2002. Microsatellite analysis of genetic divergence among populations of giant Galápagos tortoises. Molecular Ecology, 11: 2265–2283. doi: 10.1046/j.1365-294X.2002.01617.x

[7] Garrick RC, Kajdacsi B, Russello MA, Benavides E, Hyseni C, Gibbs JP, Tapia W and Caccone A. 2015. Naturally rare versus newly rare: demographic inferences on two timescales inform conservation of Galápagos giant tortoises. Ecology and Evolution, 5: 676–694. doi: 10.1002/ece3.1388

[8] Poulakakis N, Edwards DL, Chiari Y, Garrick RC, Russello MA, Benavides E, Watkins-Colwell GJ, Glaberman S, Tapia W, Gibbs JP, Cayot LJ and Caccone A. 2015. Description of a new Galápagos giant tortoise species (Chelonoidis; Testudines: Testudinidae) from Cerro Fatal on Santa Cruz island. PLoS ONE, 10: e0138779. doi: 10.1371/journal.pone.0138779

[9] Loire E, Chiari Y, Bernard A, Cahais V, Romiguier J, Nabholz B, Lourenço JM and Galtier N. 2013. Population genomics of the endangered giant Galápagos tortoise. Genome Biology, 14: R136. doi: 10.1186/gb-2013-14-12-r136

Theoretical work, initiated by Levene (1953) [1] and Dempster (1955) [2], suggests that within a given environment, the way populations are regulated and contribute to the next generation is a key factor for the maintenance of local adaptation polymorphism. In this theoretical context, hard selection describes the situation where the genetic composition of each population affects its contribution to the next generation whereas soft selection describes the case where the contribution of each population is fixed, whatever its genetic composition. Soft selection is able to maintain polymorphism, whereas hard selection invariably leads to the fixation of one of the alleles. Although the specific conditions (e.g. of migration between populations or drift level) in which this prediction holds have been studied in details by theoreticians, experimental tests have mainly failed, usually leading to the conclusion that the allele frequency dynamics was driven by other mechanisms in the experimental systems and conditions used. Gallet, Froissart and Ravigné [3] have set up a bacterial experimental system which allowed them to convincingly demonstrate that soft selection generates the conditions for polymorphism maintenance when hard selection does not, everything else being equal. The key ingredients of their experimental system are (1) the possibility to accurately produce hard and soft selection regimes when daily transferring the populations and (2) the ability to establish artificial well-characterized reproducible trade-offs. To do so, they used two genotypes resisting each one to one antibiotic and combined, across habitats, low antibiotic doses and difference in medium productivity. The experimental approach contains two complementary parts: the first one is looking at changes in the frequencies of two genotypes, initially introduced at around 50% each, over a small number of generations (ca 40) in different environments and selection regimes (soft/hard) and the second one is convincingly showing polymorphism protection by establishing that in soft selection regimes, the lowest fitness genotype is not eliminated even when introduced at low frequency. In this manuscript, a key point is the dialog between theoretical and experimental approaches. The experiments have been thought and designed to be as close as possible to the situations analysed in theoretical work. For example, the experimental polymorphism protection test (experiment 2) closely matches the equivalent analysis classically performed in theoretical approaches. This close fit between theory and experiment is clearly a strength of this study. This said, the experimental system allowing them to realise this close match also has some limitations. For example, changes in allele frequencies could only be monitored over a quite low number of generations because a longer time-scale would have allowed the contribution of de novo mutations and the likely emergence of a generalist genotype resisting to both antibiotics used to generate the local adaptation trade-offs. These limitations, as well as the actual significance of the experimental tests, are discussed in deep details in the manuscript.

References

[1] Levene H. 1953. Genetic equilibrium when more than one niche is available. American Naturalist 87: 331–333. doi: 10.1086/281792

[2] Dempster ER. 1955. Maintenance of genetic heterogeneity. Cold Spring Harbor Symposia on Quantitative Biology. 20: 25–32. doi: 10.1101/SQB.1955.020.01.005

[3] Gallet R, Froissart R, Ravigné V. 2017. Things softly attained are long retained: dissecting the impacts of selection regimes on polymorphism maintenance in experimental spatially heterogeneous environments. bioRxiv 100743; doi: 10.1101/100743

How do phenotypic plasticity and adaptive evolution interact in a novel or changing environment? Does evolution by natural selection generally reinforce initially plastic phenotypic responses, or does it instead oppose them? And to what extent does evolution of a trait involve evolution of its plasticity? These questions have lied at the heart of research on phenotypic evolution in heterogeneous environments ever since it was realized that the environment is likely to affect the expression of many (perhaps most) characters of an individual. Importantly, this broad definition of phenotypic plasticity as change in the average phenotype of a given genotype in response to its environment of development (or expression) does not involve any statement about the adaptiveness of the plastic changes. Theory on the evolution of plasticity has devoted much effort to understanding how reaction norm should evolve under different regimes of environmental change in space and time, and depending on genetic constraints on reaction norm shapes. However on an empirically ground, the questions above have mostly been addressed for individual traits, often chosen a priori for their likeliness to exhibit adaptive plasticity, and we still lack more systematic answers. These can be provided by so-called ‘phenomic’ approaches, where a large number of traits are tracked without prior information on their biological or ecological function. A problem is that the number of phenotypic characters that can be measured in an organism is virtually infinite (and to some extent arbitrary), and that scaling issues makes it difficult to compare different sets of traits. Gene-expression levels offer a partial solution to this dilemma, as they can be considered as a very large number of traits (one per typed gene) that can be measured easily and uniformly (fold change in the number of reads in RNAseq). As for any traits, expression levels of different genes may be genetically correlated, to an extent that depends on their regulation mechanism: cis-regulatory sequences that only affect expression of neighboring genes are likely to cause independent gene expression, while more systematic modifiers of expression (e.g. trans-regulators such as transcription factors) may cause correlated genetic responses of the expression of many genes. Huang and Agrawal [1] have studied plasticity and evolution of gene expression level in young larvae of populations of Drosophila melanogaster that have evolved for about 130 generations under either a constant environment (salt or cadmium), or an environment that is heterogeneous in time or space (combining salt and cadmium). They report a wealth of results, of which we summarize the most striking here. First, among genes that (i) were initially highly plastic and (ii) evolved significant divergence in expression levels between constant environment treatments, the evolved divergence is predominantly in the opposite direction to the initial plastic response. This suggests that either plasticity was initially maladaptive, or the selective pressure changed during the evolutionary process (see below). This somewhat unexpected result strikingly mirrors that from a study published last year in Nature [2], where the same pattern was found for responses of guppies to the presence of predators. However, Huang and Agrawal [1] went beyond this study by deciphering the underlying mechanisms in several interesting ways. First, they showed that change in gene expression often occurred at genes close to SNPs with differentiated frequencies across treatments (but not at genes with differentiated SNPs in their coding sequences), suggesting that cis-regulatory sequences are involved. This is also suggested by the fact that changes in gene expression are mostly caused by the increased expression of only one allele at polymorphic loci, and is a first step towards investigating the genetic underpinnings of (co)variation in gene expression levels. Another interesting set of findings concerns evolution of plasticity in treatments with variable environments. To compare the gene-expression plasticity that evolved in these treatments to an expectation, the authors considered that the expression levels in populations maintained for a long time under constant salt or cadmium had reached an optimum. The differences between these expression levels were thus assumed to predict the level of plasticity that should evolve in a heterogeneous environment (with both cadmium and salt) under perfect environmental predictability. The authors showed that plasticity did evolve more in the expected direction in heterogeneous than in constant environments, resulting in better adapted final expression levels across environments. Taken collectively, these results provide an unprecedented set of patterns that are greatly informative on how plasticity and evolution interact in constant versus changing environments. But of course, interpretations in terms of adaptive versus maladaptive plasticity are more challenging, as the authors themselves admit. Even though environmentally determined gene expression is the basic mechanism underlying the phenotypic plasticity of most traits, it is extremely difficult to relate to more integrated phenotypes for which we can understand the selection pressures, especially in multicellular organisms. The authors have recently investigated evolutionary change of quantitative traits in these selected lines, so it might be possible to establish links between reaction norms for macroscopic traits to those for gene expression levels. Such an approach would also involve tracking gene expression throughout life, rather than only in young larvae as done here, thus putting phenotypic complexity back in the picture also for expression levels. Another difficulty is that a plastic response that was originally adaptive may be replaced by an opposite evolutionary response in the long run, without having to invoke initially maladaptive plasticity. For instance, the authors mention the possibility that a generic stress response is initially triggered by cadmium, but is eventually unnecessary and costly after evolution of genetic mechanisms for cadmium detoxification (a case of so-called genetic accommodation). In any case, this study by Huang and Agrawal [1], together with the one by Ghalambor et al. last year [2], reports novel and unexpected results, which are likely to stimulate researchers interested in plasticity and evolution in heterogeneous environments for the years to come.

References

[1] Huang Y, Agrawal AF. 2016. Experimental Evolution of Gene Expression and Plasticity in Alternative Selective Regimes. PLoS Genetics 12:e1006336. doi: 10.1371/journal.pgen.1006336

[2] Ghalambor CK, Hoke KL, Ruell EW, Fischer EK, Reznick DN, Hughes KA. 2015. Non-adaptive plasticity potentiates rapid adaptive evolution of gene expression in nature. Nature 525: 372-375. doi: 10.1038/nature15256

It is a classical result in the virulence evolution literature that treatments decreasing parasite replication within the host should select for higher replication rates, thus driving increased levels of virulence if the two are correlated. There is some evidence for this in vitro but very little in the field. HIV infections in humans offer a unique opportunity to go beyond the simple predictions that treatments should favour more virulent strains because many details of this host-parasite system are known, especially the link between set-point virus load, transmission rate and virulence.

To tackle this question, Herbeck et al. [1] used a detailed individual-based model. This is original because it allows them to integrate existing knowledge from the epidemiology and evolution of HIV (e.g. recent estimates of the ‘heritability’ of set-point virus load from one infection to the next). This detailed model allows them to formulate predictions regarding the effect of different treatment policies; especially regarding the current policy switch away from treatment initiation based on CD4 counts towards universal treatment.

The results show that, perhaps as expected from the theory, treatments based on the level of remaining host target cells (CD4 T cells) do not affect virulence evolution because they do not strongly affect the virulence level that maximizes HIV’s transmission potential. However, early treatments can lead to moderate increase in virulence within several years if coverage is high enough. These results seem quite robust to variation of all the parameters in realistic ranges.

The great step forward in this model is the ability to obtain quantitative prediction regarding how a virus may evolve in response to public health policies. Here the main conclusion is that given our current knowledge in HIV biology, the risk of virulence evolution is perhaps more limited than expected from a direct application of virulence evolution model. Interestingly, the authors also conclude that recently observed increased in HIV virulence [2-3] cannot be explained by the impact of antiretroviral therapy alone; which raises the question about the main mechanism behind this increase. Finally, the authors make the interesting suggestion that “changing virulence is amenable to being monitored alongside transmitted drug resistance in sentinel surveillance”.

References

[1] Herbeck JT, Mittler JE, Gottlieb GS, Goodreau SM, Murphy JT, Cori A, Pickles M, Fraser C. 2016. Evolution of HIV virulence in response to widespread scale up of antiretroviral therapy: a modeling study. Virus Evolution 2:vew028. doi: 10.1093/ve/vew028

[2] Herbeck JT, Müller V, Maust BS, Ledergerber B, Torti C, et al. 2012. Is the virulence of HIV changing? A meta-analysis of trends in prognostic markers of HIV disease progression and transmission. AIDS 26:193-205. doi: 10.1097/QAD.0b013e32834db418

[3] Pantazis N, Porter K, Costagliola D, De Luca A, Ghosn J, et al. 2014. Temporal trends in prognostic markers of HIV-1 virulence and transmissibility: an observational cohort study. Lancet HIV 1:e119-26. doi: 10.1016/s2352-3018(14)00002-2

Much research effort has been extended toward developing systems for managing soil inhabiting insect pests of crops with entomopathogenic nematodes as biocontrol agents. Although small plot or laboratory experiments may suggest a particular insect pest is vulnerable to management in this way, it is often difficult to scale-up nematode production for application at the field- and farm scale to make such a tactic viable. Part of the problem is that entomopathogenic nematode strains must be propagated by serial passage in vivo, because storage by freezing decreases fitness. At the same time, serial propagation results in loss of virulence (ability to infect) over generations in the laboratory, a phenomenon called attenuation.

To probe the underlying reasons for development of attenuation, as a prerequisite to designing strategies to mitigate it, Shapiro-Ilan and Raymond [1] turned to evolutionary theory of social conflict as a possible explanatory framework. Virulence of entomopathogenic nematodes depends on a combination of virulence factors, like various proteases, secreted by both the nematode and symbiotic bacteria to overcome host defenses. Attenuation is characterized in part by a reduced production of these factors. Invasion of a host involves simultaneous attack by a group of nematodes ("cooperators"), which together neutralize host defenses enough to allow individuals to successfully invade. "Cheaters" in the invading population can avoid the metabolic costs of producing virulence factors while reaping the benefits of infecting the host made vulnerable by the cooperators in the population. The authors hypothesize that an increase in frequency of cheaters may contribute to attenuation of virulence during serial propagation in the laboratory. The evolutionary dynamics of cheater frequency in a population have been explored in many contexts as part of kin selection theory. Cheaters can increase in a population by outcompeting cooperators in a host if overall relatedness within the invading population is low. Conversely, frequency of altruism, or costly cooperation, increases in a population if relatedness is high, which is enhanced by low effective dispersal. However, a population that is too isolated can suffer from inbreeding effects, and competition will occur mainly among relatives, which decreases the fitness benefits of altruism.

Shapiro-Ilan and Raymond [1] tested changes in virulence and reproductive output in a serially propagated entomopathogenic nematode, Heterorhabditis floridensis. They compared lines of high or low relatedness, manipulated via multiplicity of infection (MOI) rates (where a low dose of nematodes gives high relatedness and a high dose gives low relatedness); and under global or local competition, manipulated by pooling populations emerging from all or only two host cadavers per generation, respectively. As predicted, treatments of high relatedness (low MOI) and global competition had the greatest level of reproduction, while all lines of low relatedness (high MOI) evolved decreased reproduction and decreased virulence, which led to extinction. The key finding was that lines in the high relatedness (low MOI) and low (local) competition treatment exhibited the most stable virulence through the 12 generations tested. Thus, to minimize attenuation of virulence while maintaining fitness of recently isolated entomopathogenic nematodes, the authors recommend insect hosts be inoculated with low doses of nematodes from inocula pools from as few cadavers as possible.

The application of evolutionary theory, with a clever experimental design, to an important problem in pest management makes this paper particularly noteworthy.

Reference

[1] Shapiro-Ilan D, Raymond B. 2016. Limiting opportunities for cheating stabilizes virulence in insect parasitic nematodes. Evolutionary Applications 9:462-470. doi: 10.1111/eva.12348