The preprint by Pratlong et al. [1] is a well thought quest for genomic regions involved in local adaptation to depth in a species a red coral living the Mediterranean Sea. It first describes a pattern of structuration and then attempts to find candidate genes involved in local adaptation by contrasting deep with shallow populations. Although the pattern of structuration is clear and meaningful, the candidate genomic regions involved in local adaptation remain to be confirmed. Two external reviewers and myself found this preprint particularly interesting regarding the right-mindedness of the authors in front of the difficulties they encounter during their experiments. The discussions on the pros and cons of the approach are very sound and can be easily exported to a large number of studies that hunt for local adaptation. In this sense, the lessons one can learn by reading this well documented manuscript are certainly valuable for a wide range of evolutionary biologists.
More precisely, the authors RAD-sequenced 6 pairs of 'shallow vs deep' samples located in 3 geographical sea areas (Banyuls, Corsica and Marseilles). They were hoping to detect genes involved in the adaptation to depth, if there were any. They start by assessing the patterns of structuration of the 6 samples using PCA and AMOVA [2] and also applied the STRUCTURE [3] assignment software. They show clearly that the samples were mostly differentiated between geographical areas and that only 1 out the 3 areas shows a pattern of isolation by depth (i.e. Marseille). They nevertheless went on and scanned for variants that are highly differentiated in the deep samples when compared to the shallow paired samples in Marseilles, using an Fst outliers approach [4] implemented in the BayeScEnv software [5]. No clear functional signal was in the end detected among the highly differentiated SNPs, leaving a list of candidates begging for complementary data.
The scan for local adaptation using signatures of highly divergent regions is a classical problem of population genetics. It has been applied on many species with various degrees of success. This study is a beautiful example of a well-designed study that did not give full satisfactory answers. Readers will especially appreciate the honesty and the in-depth discussions of the authors while exposing their results and their conclusions step by step.

References

[1] Pratlong, M., Haguenauer, A., Brener, K., Mitta, G., Toulza, E., Garrabou, J., Bensoussan, N., Pontarotti P., & Aurelle, D. (2018). Separate the wheat from the chaff: genomic scan for local adaptation in the red coral Corallium rubrum. bioRxiv, 306456, ver. 3 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/306456
[2] Excoffier, L., Smouse, P. E. & Quattro, J. M. (1992). Analysis of molecular variance inferred from metric distances among DNA haplotypes: application to human mitochondrial DNA restriction data. Genetics, 131(2), 479-491.
[3] Pritchard, J. K., Stephens, M., & Donnelly, P. (2000). Inference of population structure using multilocus genotype data. Genetics, 155(2), 945-959.
[4] Lewontin, R. C., & Krakauer, J. (1973). Distribution of gene frequency as a test of the theory of the selective neutrality of polymorphisms. Genetics, 74(1), 175-195.
[5] de Villemereuil, P., & Gaggiotti, O. E. (2015). A new FST‐based method to uncover local adaptation using environmental variables. Methods in Ecology and Evolution, 6(11), 1248-1258. doi: 10.1111/2041-210X.12418

Artificial selection experiments are key experiments in evolutionary biology. The demonstration that application of selective pressure across multiple generations results in heritable phenotypic changes is a tangible and reproducible proof of the evolution by natural selection.
Artificial selection experiments are used to evaluate the joint effects of selection on multiple traits, their genetic covariances and differences in responses in different environments. Most studies on artificial selection experiments report and base their analyses on phenotypic changes [1]. More recently, changes in allele frequency and other patterns of molecular genetic diversity have been used to identify genomic locations where selection has had an effect. However, so far the changes in gene expression have not been in the focus of artificial selection experiment studies (see [2] for an example though).
In plants, one of the most famous artificial selection experiments is the Illinois Corn Experiment where maize (Zea mays) is selected for oil and protein content [3], but in addition, similar experiments have been conducted also for other traits in maize. In Saclay divergent selection experiment [4] two maize inbred lines (F252 and MBS847) have been selected for early and late flowering for 13 generations, resulting in two week difference in flowering time.
In ”Transcriptomic response to divergent selection for flowering time in maize reveals convergence and key players of the underlying gene regulatory network ” [5] Maud Tenaillon and her coworkers study the gene expression differences among these two independently selected maize populations. Their experiments cover two years in field conditions and they use samples of shoot apical meristem at three different developmental stages: vegetative, transitioning and reproductive. They use RNA-seq transcriptome level differences and qRT-PCR for gene expression pattern investigation. The work is continuation to earlier genetic and phenotypic studies on the same material [4, 6].
The reviewers and I agree that dataset is unique and its major benefit is that it has been obtained from field conditions similar to those that species may face under natural setting during selection. Their tissue sampling is supported by flowering time phenotypic observations and covers the developmental transition stage, making a good effort to identify key transcriptional and phenotypic changes and their timing affected by selection.
Tenaillon et al. [5] identify more than 2000 genes that are differentially expressed among early and late flowering populations. Expectedly, they are enriched for known flowering time genes. As they point out, differential expression of thousands of genes does not mean that they all were independently affected by selection, but rather that the whole transcriptional network has shifted, possibly due to just few upstream or hub-genes. Also, the year-to-year variation had smaller effect in gene expression compared to developmental stage or genetic background, possibly indicating selection for stability across environmental fluctuation for such an important phenotype as flowering time.
Another noteworthy observation is that they find convergent patterns of transcriptional changes among the two selected lines. 115 genes expression patterns are shifted due to selection in both genetic backgrounds. This convergent pattern can be a result of either selection on standing variation or de novo mutations. The data does not allow testing which process is underlying the observed convergence. However, their results show that this is an interesting future question that can be addressed using genotype and gene expression data from the same ancestral and derived material and possibly their hybrids.

References

[1] Hill, W. G., & Caballero, A. (1992). Artificial selection experiments. Annual Review of Ecology and Systematics, 23(1), 287-310. doi: 10.1146/annurev.es.23.110192.001443
[2] Konczal, M., Babik, W., Radwan, J., Sadowska, E. T., & Koteja, P. (2015). Initial molecular-level response to artificial selection for increased aerobic metabolism occurs primarily through changes in gene expression. Molecular biology and evolution, 32(6), 1461-1473. doi: 10.1093/molbev/msv038
[3] Moose, S. P., Dudley, J. W., & Rocheford, T. R. (2004). Maize selection passes the century mark: a unique resource for 21st century genomics. Trends in plant science, 9(7), 358-364. doi: 10.1016/j.tplants.2004.05.005
[4] Durand, E., Tenaillon, M. I., Ridel, C., Coubriche, D., Jamin, P., Jouanne, S., Ressayre, A., Charcosset, A. and Dillmann, C. (2010). Standing variation and new mutations both contribute to a fast response to selection for flowering time in maize inbreds. BMC evolutionary biology, 10(1), 2. doi: 10.1186/1471-2148-10-2
[5] Tenaillon, M. I., Seddiki, K., Mollion, M., Le Guilloux, M., Marchadier, E., Ressayre, A. and Dillmann C. (2019). Transcriptomic response to divergent selection for flowering time in maize reveals convergence and key players of the underlying gene regulatory network. BioRxiv, 461947 ver. 5 peer-reviewed and recommended by PCI Evolutionary Biology. doi: 10.1101/461947
[6] Durand, E., Tenaillon, M. I., Raffoux, X., Thépot, S., Falque, M., Jamin, P., Bourgais A., Ressayre, A. and Dillmann, C. (2015). Dearth of polymorphism associated with a sustained response to selection for flowering time in maize. BMC evolutionary biology, 15(1), 103. doi: 10.1186/s12862-015-0382-5

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

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

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

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

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

References

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

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

Partitioning the effects of vicariance and low dispersal has been a long-standing problem in historical biogeography and phylogeography. While the term “vicariance” refers to divergence in allopatry, caused by some physical (geological, geographical) or climatic barriers (e.g. Rosen 1978), isolation by distance refers to the genetic differentiation of remote populations due to the physical distance separating them, when the latter surpasses the scale of dispersal (Wright 1938, 1940, 1943). 

Vicariance and dispersal have long been considered as separate forces leading to separate scenarii of speciation (e.g. reviewed in Hickerson and Meyer 2008). Nevertheless, these two processes are strongly linked, as, for example, vicariance theory relies on the assumption that ancestral lineages were once linked by dispersal prior to physical or climatic isolation (Rosen 1978). Low dispersal and vicariance are not mutually exclusive, and distinguishing these two processes in heterogeneous landscapes, especially for poor dispersers, remains therefore a severe challenge. For example, low dispersal (and/or small population size) can give rise to geographic patterns consistent with a phylogeographic break and be mistaken for geographic isolation (Irwin 2002, Kuo and Avise 2005).

The study of Rancilliac and colleagues (2023) is at the heart of this issue. It focuses on a nominal lizard species, the red-tailed spiny-footed lizard (Acanthodactylus erythrurus, Squamata: Lacertidae), which has a wide spatial distribution (from the Maghreb to the Iberian Peninsula), is found in a variety of different habitats, and has a wide range of morphological traits that do not always correlate with phylogeny. The main question is the following: have “the morphological and ecological diversification of this group been produced by vicariance and lineage diversification, or by local adaptation in the face of historical gene flow?” To tackle this question, the authors used sequence data from multiple mitochondrial and nuclear markers and a nested analysis workflow integrating phylogeography, multiple correspondence analyses and a relatively novel approach to IBD testing (Hausdorf & Henning, 2020). The latter is based on regression analysis and was shown to be less prone to error than the traditional (partial) Mantel test. 

While this set of methods allowed the partitioning of the effect of isolation by distance and vicariance in shaping contemporary genetic diversity in red-tailed spiny-footed lizards, some of the evolutionary history of this species complex remains blurred by ongoing gene flow and admixture, retention of ancestral polymorphism, or selection. The lack of congruence between mitochondrial and nuclear gene trees once again warns us that proposing evolutionary scenarii based on individual gene trees is a risky business. 

References

Hausdorf B, Hennig C (2020) Species delimitation and geography. Molecular Ecology Resources, 20, 950–960. https://doi.org/10.1111/1755-0998.13184

Hickerson MJ, Meyer CP (2008) Testing comparative phylogeographic models of marine vicariance and dispersal using a hierarchical Bayesian approach. BMC Evolutionary Biology, 8, 322. https://doi.org/10.1186/1471-2148-8-322

Irwin DE (2002) Phylogeographic breaks without geographic barriers to gene flow. Evolution, 56, 2383–2394. https://doi.org/10.1111/j.0014-3820.2002.tb00164.x

Kuo C-H, Avise JC (2005) Phylogeographic breaks in low-dispersal species: the emergence of concordance across gene trees. Genetica, 124, 179–186. https://doi.org/10.1007/s10709-005-2095-y

Rancilhac L, Miralles A, Geniez P, Mendez-Aranda D, Beddek M, Brito JC, Leblois R, Crochet P-A (2023) Phylogeographic breaks and how to find them: An empirical attempt at separating vicariance from isolation by distance in a lizard with restricted dispersal. bioRxiv, 2022.09.30.510256, ver. 4 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.09.30.510256

Rosen DE (1978) Vicariant Patterns and Historical Explanation in Biogeography. Systematic Biology, 27, 159–188. https://doi.org/10.2307/2412970

Wright, S (1938) Size of population and breeding structure in relation to evolution. Science 87:430-431.

Wright S (1940) Breeding Structure of Populations in Relation to Speciation. The American Naturalist, 74, 232–248. https://doi.org/10.1086/280891

Wright S (1943) Isolation by distance. Genetics, 28, 114–138. https://doi.org/10.1093/genetics/28.2.114

Studying the effects of speciation, hybridisation, and evolutionary outcomes following reproduction from divergent populations is a major research area in evolutionary genetics [1]. There are two phenomena that have been the focus of contemporary research. First, a classic concept is the formation of ‘Bateson-Dobzhansky-Muller’ incompatibilities (BDMi) [2–4] that negatively affect hybrid fitness. Here, two diverging populations accumulate mutations over time that are unique to that subpopulation. If they subsequently meet, then these mutations might negatively interact, leading to a loss in fitness or even a complete lack of reproduction. BDMi formation can be complex, involving multiple genes and the fitness changes can depend on the direction of introgression [5]. Second, such secondary contact can instead lead to heterosis, where offspring are fitter than their parental progenitors [6].

Understanding which outcomes are likely to arise require one to know the potential fitness effects of mutations underlying reproductive isolation, to determine whether they are likely to reduce or enhance fitness when hybrids are formed. This is far from an easy task, as it requires one to track mutations at several loci, along with their effects, across a fitness landscape.

The work of De Sanctis et al. [7] neatly fills in this knowledge gap, by creating a general mathematical framework for describing the consequences of a cross from two divergent populations. The derivations are based on Fisher’s Geometric Model, which is widely used to quantify selection acting on a general fitness landscape that is affected by several biological traits [8,9], and has previously been used in theoretical studies of hybridisation [10–12]. By doing so, they are able to decompose how divergence at multiple loci affects offspring fitness through both additive and dominance effects.

A key result arising from their analyses is demonstrating how offspring fitness can be captured by two main functions. The first one is the ‘net effect of evolutionary change’ that, broadly defined, measures how phenotypically divergent two populations are. The second is the ‘total amount of evolutionary change’, which reflects how many mutations contribute to divergence and the effect sizes captured by each of them. The authors illustrate these measurements using simulations covering different scenarios, demonstrating how different parental states can lead to similar fitness outcomes. They also propose experimental methods to measure the underlying mutational effects.

This study neatly demonstrates how complex genetic phenomena underlying hybridisation can be captured using fairly simple mathematical formulae. This powerful approach will thus open the door for future research to investigate hybridisation in more detail, whether it is by expanding on these theoretical models or using the elegant outcomes to quantify fitness effects in experiments.

References

1. Coyne JA, Orr HA. Speciation. Sunderland, Mass: Sinauer Associates; 2004.
2. Bateson W, Seward A. Darwin and modern science. Heredity and variation in modern lights. 1909;85: 101. https://doi.org/10.1017/CBO9780511693953.007
3. Dobzhansky T. Genetics and the Origin of Species. Columbia university press; 1937.
4. Muller HJ. Isolating mechanisms, evolution and temperature. Biol Symp. 1942;6: 71-125.
5. Fraïsse C, Elderfield JAD, Welch JJ. The genetics of speciation: are complex incompatibilities easier to evolve? J Evol Biol. 2014;27: 688-699. https://doi.org/10.1111/jeb.12339
6. Birchler JA, Yao H, Chudalayandi S, Vaiman D, Veitia RA. Heterosis. The Plant Cell. 2010;22: 2105-2112. https://doi.org/10.1105/tpc.110.076133
7. De Sanctis B, Schneemann H, Welch JJ. How does the mode of evolutionary divergence affect reproductive isolation? bioRxiv. 2022. 2022.03.08.483443 version 4. https://doi.org/10.1101/2022.03.08.483443 
8. Fisher RA. The genetical theory of natural selection. Oxford: The Clarendon Press; 1930. https://doi.org/10.5962/bhl.title.27468 
9. Tenaillon O. The Utility of Fisher's Geometric Model in Evolutionary Genetics. Annu Rev Ecol Evol Syst. 2014;45: 179-201. https://doi.org/10.1146/annurev-ecolsys-120213-091846
10. Barton NH. The role of hybridization in evolution. Molecular Ecology. 2001;10: 551-568. https://doi.org/10.1046/j.1365-294x.2001.01216.x 
11. Chevin L-M, Decorzent G, Lenormand T. Niche Dimensionality and The Genetics of Ecological Speciation. Evolution. 2014;68: 1244-1256. https://doi.org/10.1111/evo.12346 
12. Fraïsse C, Gunnarsson PA, Roze D, Bierne N, Welch JJ. The genetics of speciation: Insights from Fisher's geometric model. Evolution. 2016;70: 1450-1464. https://doi.org/10.1111/evo.12968