Each individual multicellular organism, each of our bodies, is a small universe. Every living surface -skin, cuticle, bark, mucosa- is the home place to milliards of bacteria, fungi and viruses. They constitute our microbiota. Some of them are essential for certain organisms. Other could not live without their hosts. For many species, the relationship between host and microbiota is so close that their histories are inseparable. The recognition of this biological inextricability has led to the notion of holobiont as the organism ensemble of host and microbiota. When individuals of a particular animal or plant species expand their geographical range, it is the holobiont that expands. And these processes of migration, expansion and colonization are often accompanied by evolutionary and ecological innovations in the interspecies relationships, at the macroscopic level (e.g. novel predator-prey or host-parasite interactions) and at the microscopic level (e.g. changes in the microbiota composition). From the human point of view, these novel interactions can be economically disastrous if they involve and threaten important crop or cattle species. And this is especially worrying in the present context of genetic standardization and intensification for mass-production on the one hand, and of climate change on the other.

With this perspective, the international team led by Amandine Cornille presents a study aiming at understanding the evolutionary history of the rosy apple aphid Dysaphis plantaginea Passerini, a major pest of the cultivated apple tree Malus domestica Borkh (1). The apple tree was probably domesticated in Central Asia, and later disseminated by humans over the world in different waves, and it was probably introduced in Europe by the Greeks. It is however unclear when and where D. plantaginea started parasitizing the cultivated apple tree. The ancestral D. plantaginea could have already infected the wild ancestor of current cultivated apple trees, but the aphid is not common in Central Asia. Alternatively, it may have gained access only later to the plant, possibly via a host jump, from Pyrus to Malus that may have occurred in Asia Minor or in the Caucasus. In the present preprint, Olvera-Vázquez and coworkers have analysed over 650 D. plantaginea colonies from 52 orchards in 13 countries, in Western, Central and Eastern Europe as well as in Morocco and the USA. The authors have analysed the genetic diversity in the sampled aphids, and have characterized as well the composition of the associated endosymbiont bacteria. The analyses detect substantial recent admixture, but allow to identify aphid subpopulations slightly but significantly differentiated and isolated by distance, especially those in Morocco and the USA, as well as to determine the presence of significant gene flow. This process of colonization associated to gene flow is most likely indirectly driven by human interactions. Very interestingly, the data show that this genetic diversity in the aphids is not reflected by a corresponding diversity in the associated microbiota, largely dominated by a few Buchnera aphidicola variants. In order to determine polarity in the evolutionary history of the aphid-tree association, the authors have applied approximate Bayesian computing and machine learning approaches. Albeit promising, the results are not sufficiently robust to assess directionality nor to confidently assess the origin of the crop pest. Despite the large effort here communicated, the authors point to the lack of sufficient data (in terms of aphid isolates), especially originating from Central Asia. Such increased sampling will need to be implemented in the future in order to elucidate not only the origin and the demographic history of the interaction between the cultivated apple tree and the rosy apple aphid. This knowledge is needed to understand how this crop pest struggles with the different seasonal and geographical selection pressures while maintaining high genetic diversity, conspicuous gene flow, differentiated populations and low endosymbiontic diversity.

References

1. Olvera-Vazquez SG, Remoué C, Venon A, Rousselet A, Grandcolas O, Azrine M, Momont L, Galan M, Benoit L, David GM, Alhmedi A, Beliën T, Alins G, Franck P, Haddioui A, Jacobsen SK, Andreev R, Simon S, Sigsgaard L, Guibert E, Tournant L, Gazel F, Mody K, Khachtib Y, Roman A, Ursu TM, Zakharov IA, Belcram H, Harry M, Roth M, Simon JC, Oram S, Ricard JM, Agnello A, Beers EH, Engelman J, Balti I, Salhi-Hannachi A, Zhang H, Tu H, Mottet C, Barrès B, Degrave A, Razmjou J, Giraud T, Falque M, Dapena E, Miñarro M, Jardillier L, Deschamps P, Jousselin E, Cornille A (2021) Large-scale geographic survey provides insights into the colonization history of a major aphid pest on its cultivated apple host in Europe, North America and North Africa. bioRxiv, 2020.12.11.421644, ver. 3 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2020.12.11.421644

The total amount of DNA utilized to store hereditary information varies immensely among eukaryotic organisms. Single copy genome sizes – disregarding differences due to ploidy - differ by more than three orders of magnitude ranging from a few million nucleotides (Mb) to hundreds of billions (Gb). With the ever-increasing availability of fully sequenced genomes we now know that most of the difference is due either to whole genome duplication or to variation in the abundance of repetitive elements. Regarding repetitive elements, the evolutionary forces underlying the large variation 'allowing' more or less elements in a genome remain largely elusive. A tentative correlation between an organism's complexity (however this may be adequately measured) and genome size, the so called C-value paradox [1], has long been dismissed. Studies testing for selection on secondary phenotypic effects associated with genome size (cell size, metabolic rates, nutrient availability) have yielded mixed results. Nonadaptive theories capitalizing on a role of deleterious insertion-deletion mutations and genetic drift as the main drivers have likewise received mixed support [2-3]. Overall, most evidence was derived from analyses across broad taxonomical scales [4-6].

Lefébure and colleagues [7] take a different approach. They confine their considerations to a homogeneous, restricted taxonomical group, isopod crustaceans of the superfamily Aselloidea. This taxonomic focus allows the authors to circumvent many of the confounding factors such as phylogenetic inertia, life history divergence and mutation rate variation that tend to trouble analyses across broad taxonomic timescales. Another important feature of the chosen system is the evolutionary independent transition of habitat use that has occurred at least 11 times. One group of species inhabits subterranean ecosystems (groundwater), another group thrives on surface water. Populations of the former live in low-energy habitats and are expected to be outnumbered by their surface dwelling relatives. Interestingly – and a precondition for the study - the groundwater species have significantly larger genomes (up to 137%). With this unique set-up, the authors are able to investigate the link between genome size and evolutionary forces related to a proxy of long-term population size by removing many of the confounding factors a priori.

Upfront, we learn that the dN/dS ratio is higher in the groundwater species. This may either suggest prevalent positive selection or lower efficacy of purifying selection (relaxed constraint) in the group of species in which population sizes are expected to be low. Using a series of population genetic analyses the authors provide compelling evidence for the latter. Analyses are carefully conducted and include models for estimating the intensity and frequency of purifying and positive selection, the DoS (direction of selection) and α statistic. Next the authors also exclude the possibility that increased dN/dS of the subterranean groundwater species may be due to nonfunctionalization, which may result from the subterranean lifestyle.

Overall, these analyses suggest relaxed constraint in smaller populations as the most plausible alternative to explain increased dN/dS ratios. In addition to the efficacy of selection, the authors estimate the timing of the ecological transition under the rationale that the amount of time a species may have been exposed to the subterranean habitat may reflect long term population sizes. To calibrate the 'colonization clock' they apply a neat trick based on the degree of degeneration of the opsin gene (as vision tends to get lost in these habitats). When finally testing which parameters may explain differences in genome size all factors – ecological status, selection efficiency as measured by dN/dS and colonization time - turned out to be significant predictors. Direct estimates of the short term effective population size Ne from polymorphism data, however, did not correlate with genome size. Ruling out the effect of other co-variates such as body size and growth rate the authors conclude that genome size was overall best predicted by long-term population size change upon habitat shift. In that the authors provide convincing evidence that the increase in genome size is linked to a decrease in long-term reduction of selection efficiency of subterranean species. Assuming a bias for insertion mutations over deletion mutations (which is usually the case in eukaryotes) this result is in agreement with the theory of mutational hazard [4-6]. This theory proposed by Michael Lynch postulates that the accumulation of non-functional DNA has a weak deleterious effect that can only be efficiently opposed by natural selection in species with high Ne.

In conclusion, Lefébure and colleagues provide novel and welcome evidence supporting a 'neutralist' hypothesis of genome size evolution without the need to invoke an adaptive component. Methodologically, the study cautions against the common use of polymorphism-based estimates of Ne which are often obfuscated by transitory demographic change. Instead, alternative measures of selection efficacy linked to long-term population size may serve as better predictors of genome size. We hope that this study will stimulate additional work testing the link between Ne and genome size variation in other taxonomical groups [8-9]. Using genome sequences instead of the transcriptome approach applied here may concomitantly further our understanding of the molecular mechanisms underlying genome size change.

References

[1] Thomas, CA Jr. 1971. The genetic organization of chromosomes. Annual Review of Genetics 5: 237–256. doi: 10.1146/annurev.ge.05.120171.001321

[2] Ågren JA, Greiner S, Johnson MTJ, Wright SI. 2015. No evidence that sex and transposable elements drive genome size variation in evening primroses. Evolution 69: 1053–1062. doi: 10.1111/evo.12627

[3] Bast J, Schaefer I, Schwander T, Maraun M, Scheu S, Kraaijeveld K. 2016. No accumulation of transposable elements in asexual arthropods. Molecular Biology and Evolution 33: 697–706. doi: 10.1093/molbev/msv261

[4] Lynch M. 2007. The Origins of Genome Architecture. Sinauer Associates.

[5] Lynch M, Bobay LM, Catania F, Gout JF, Rho M. 2011. The repatterning of eukaryotic genomes by random genetic drift. Annual Review of Genomics and Human Genetics 12: 347–366. doi: 10.1146/annurev-genom-082410-101412

[6] Lynch M, Conery JS. 2003. The origins of genome complexity. Science 302: 1401–1404. doi: 10.1126/science.1089370

[7] Lefébure T, Morvan C, Malard F, François C, Konecny-Dupré L, Guéguen L, Weiss-Gayet M, Seguin-Orlando A, Ermini L, Der Sarkissian C, Charrier NP, Eme D, Mermillod-Blondin F, Duret L, Vieira C, Orlando L, and Douady CJ. 2017. Less effective selection leads to larger genomes. Genome Research 27: 1016-1028. doi: 10.1101/gr.212589.116

[8] Lower SS, Johnston JS, Stanger-Hall KF, Hjelmen CE, Hanrahan SJ, Korunes K, Hall D. 2017. Genome size in North American fireflies: Substantial variation likely driven by neutral processes. Genome Biolology and Evolution 9: 1499–1512. doi: 10.1093/gbe/evx097

[9] Sessegolo C, Burlet N, Haudry A. 2016. Strong phylogenetic inertia on genome size and transposable element content among 26 species of flies. Biology Letters 12: 20160407. doi: 10.1098/rsbl.2016.0407

The human genome not only encodes for biological functions and for what makes us human, it also encodes the population history of our ancestors. Changes in past population sizes, for example, affect the distribution of times to the most recent common ancestor (tMRCA) of genomic segments, which in turn can be inferred by sophisticated modelling along the genome.
A key framework for such modelling of local tMRCA tracts along genomes is the Sequentially Markovian Coalescent (SMC) (McVean and Cardin 2005, Marjoram and Wall 2006) . The problem that the SMC solves is that the mosaic of local tMRCAs along the genome is unknown, both in their actual ages and in their positions along the genome. The SMC allows to effectively sum across all possibilities and handle the uncertainty probabilistically. Several important tools for inferring the demographic history of a population have been developed built on top of the SMC, including PSMC (Li and Durbin 2011), diCal (Sheehan et al 2013), MSMC (Schiffels and Durbin 2014), SMC++ (Terhorst et al 2017), eSMC (Sellinger et al. 2020) and others.
In this paper, Sellinger, Abu Awad and Tellier (2020) review these SMC-based methods and provide a coherent simulation design to comparatively assess their strengths and weaknesses in a variety of demographic scenarios (Sellinger, Abu Awad and Tellier 2020). In addition, they used these simulations to test how breaking various key assumptions in SMC methods affects estimates, such as constant recombination rates, or absence of false positive SNP calls.
As a result of this assessment, the authors not only provide practical guidance for researchers who want to use these methods, but also insights into how these methods work. For example, the paper carefully separates sources of error in these methods by observing what they call “Best-case convergence” of each method if the data behaves perfectly and separating that from how the method applies with actual data. This approach provides a deeper insight into the methods than what we could learn from application to genomic data alone.
In the age of genomics, computational tools and their development are key for researchers in this field. All the more important is it to provide the community with overviews, reviews and independent assessments of such tools. This is particularly important as sometimes the development of new methods lacks primary visibility due to relevant testing material being pushed to Supplementary Sections in papers due to space constraints. As SMC-based methods have become so widely used tools in genomics, I think the detailed assessment by Sellinger et al. (2020) is timely and relevant.
In conclusion, I recommend this paper because it bridges from a mere review of the different methods to an in-depth assessment of performance, thereby addressing both beginners in the field who just seek an initial overview, as well as experienced researchers who are interested in theoretical boundaries and assumptions of the different methods.

References

[1] Li, H., and Durbin, R. (2011). Inference of human population history from individual whole-genome sequences. Nature, 475(7357), 493-496. doi: https://doi.org/10.1038/nature10231
[2] Marjoram, P., and Wall, J. D. (2006). Fast"" coalescent"" simulation. BMC genetics, 7(1), 16. doi: https://doi.org/10.1186/1471-2156-7-16
[3] McVean, G. A., and Cardin, N. J. (2005). Approximating the coalescent with recombination. Philosophical Transactions of the Royal Society B: Biological Sciences, 360(1459), 1387-1393. doi: https://doi.org/10.1098/rstb.2005.1673
[4] Schiffels, S., and Durbin, R. (2014). Inferring human population size and separation history from multiple genome sequences. Nature genetics, 46(8), 919-925. doi: https://doi.org/10.1038/ng.3015
[5] Sellinger, T. P. P., Awad, D. A., Moest, M., and Tellier, A. (2020). Inference of past demography, dormancy and self-fertilization rates from whole genome sequence data. PLoS Genetics, 16(4), e1008698. doi: https://doi.org/10.1371/journal.pgen.1008698
[6] Sellinger, T. P. P., Awad, D. A. and Tellier, A. (2020) Limits and Convergence properties of the Sequentially Markovian Coalescent. bioRxiv, 2020.07.23.217091, ver. 3 peer-reviewed and recommended by PCI Evolutionary Biology. doi: https://doi.org/10.1101/2020.07.23.217091
[7] Sheehan, S., Harris, K., and Song, Y. S. (2013). Estimating variable effective population sizes from multiple genomes: a sequentially Markov conditional sampling distribution approach. Genetics, 194(3), 647-662. doi: https://doi.org/10.1534/genetics.112.149096
[8] Terhorst, J., Kamm, J. A., and Song, Y. S. (2017). Robust and scalable inference of population history from hundreds of unphased whole genomes. Nature genetics, 49(2), 303-309. doi: https://doi.org/10.1038/ng.3748

Dating nodes in a phylogeny is an important problem in evolution and is typically performed by using molecular clocks and fossil age estimates [1]. The manuscript by Chauve et al. [2] reports a novel method, which uses lateral gene transfers to help ordering nodes in a species tree. The idea is that a lateral gene transfer can only occur between two species living at the same time, which indirectly informs on node relative ages in a phylogeny: the donor species cannot be more recent than the recipient species. Horizontal gene transfers are increasingly recognized as frequent, even in eukaryotes, and especially in micro-organisms that have little fossil records [3-7]. Yet, such an important source of information has been very rarely used so far for inferring relative node ages in phylogenies. In this context, the method by Chauve et al. [2] represents an innovative and original approach to a difficult problem. An obvious limitation of the approach is that it relies on inferences of horizontal transfers, which detection is in itself a difficult problem. Incomplete taxon sampling, or the extinction of the true donor lineage may render patterns difficult to interpret in a temporary fashion. Yet, for clades with no fossils this may be the only piece of information we have at hand, and the growing amount of sequence data is likely to minimize issues derived from incomplete sampling.

The developed method, MaxTiC (for Maximal Time Consistency) [2], represents a very nice application of theoretical developments on the well-known « Feedback Arc Set » computer science problem to the evolutionary question of ordering nodes in a phylogeny. MaxTiC uses as input a species tree and a set of time constraints based on lateral gene transfers inferred using other softwares, and minimizes conflicts between node ordering and these time constraints. The application of MaxTiC on simulated datasets indicated that node ordering was fairly accurate [2]. MaxTiC is implemented in a freely available software, which represents original and relevant contribution to the field of evolutionary biology.

References

[1] Donoghue P and Smith M, editors. 2003. Telling the evolutionary time. CRC press.

[2] Chauve C, Rafiey A, Davin AA, Scornavacca C, Veber P, Boussau B, Szöllősi GJ, Daubin V and Tannier E. 2017. MaxTiC: Fast ranking of a phylogenetic tree by Maximum Time Consistency with lateral gene transfers. bioRxiv 127548, ver. 6 of 6th November 2017. doi: 10.1101/127548

[3] Ropars J, Rodríguez de la Vega RC, Lopez-Villavicencio M, Gouzy J, Sallet E, Debuchy R, Dupont J, Branca A and Giraud T. 2015. Adaptive horizontal gene transfers between multiple cheese-associated fungi. Current Biology 19, 2562–2569. doi: 10.1016/j.cub.2015.08.025

[4] Novo M, Bigey F, Beyne E, Galeote V, Gavory F, Mallet S, Cambon B, Legras JL, Wincker P, Casaregola S and Dequin S. 2009. Eukaryote-to-eukaryote gene transfer events revealed by the genome sequence of the wine yeast Saccharomyces cerevisiae EC1118. Proceeding of the National Academy of Science USA, 106, 16333–16338. doi: 10.1073/pnas.0904673106

[5] Naranjo-Ortíz MA, Brock M, Brunke S, Hube B, Marcet-Houben M, Gabaldón T. 2016. Widespread inter- and intra-domain horizontal gene transfer of d-amino acid metabolism enzymes in Eukaryotes. Frontiers in Microbiology 7, 2001. doi: 10.3389/fmicb.2016.02001

[6] Alexander WG, Wisecaver JH, Rokas A, Hittinger CT. 2016. Horizontally acquired genes in early-diverging pathogenic fungi enable the use of host nucleosides and nucleotides. Proceeding of the National Academy of Science USA. 113, 4116–4121. doi: 10.1073/pnas.1517242113

[7] Marcet-Houben M, Gabaldón T. 2010. Acquisition of prokaryotic genes by fungal genomes. Trends in Genetics. 26, 5–8. doi: 10.1016/j.tig.2009.11.007

The genome of eukaryotic species is a complex structure that experience many different interactions within itself and with the surrounding environment. The genetic architecture of a phenotype (that is, the set of genetic elements affecting a trait of the organism) plays a fundamental role in understanding the adaptation process of a species to, for example, different climate environments, or to its interaction with other species. Thus, it is fundamental to study the different aspects of the genetic architecture of the species and its relationship with its surronding environment. Aspects such as modularity (the number of genetic units and the degree to which each unit is affecting a trait of the organism), pleiotropy (the number of different effects that a genetic unit can have on an organism) or linkage (the degree of association between the different genetic units) are essential to understand the genetic architecture and to interpret the effects of selection on the genome. Indeed, the knowledge of the different aspects of the genetic architecture could clarify whether genes are affected by multiple aspects of the environment or, on the contrary, are affected by only specific aspects [1,2].

The work performed by Lotterhos et al. [3] sought to understand the genetic architecture of the adaptation to different environments in lodgepole pine (Pinus contorta), considering as candidate SNPs those previously detected as a result of its extreme association patterns to different environmental variables or to extreme population differentiation. This consideration is very important because the study is only relevant if the studied markers are under the effect of selection. Otherwise, the genetic architecture of the adaptation to different environments would be masked by other (neutral) kind of associations that would be difficult to interpret [4,5]. In order to understand the relationship between genetic architecture and adaptation, it is relevant to detect the association networks of the candidate SNPs with climate variables (a way to measure modularity) and if these SNPs (and loci) are affected by single or multiple environments (a way to measure pleiotropy).

The authors used co-association networks, an innovative approach in this field, to analyse the interaction between the environmental information and the genetic polymorphism of each individual. This methodology is more appropriate than other multivariate methods - such as analysis based on principal components - because it is possible to cluster SNPs based on associations with similar environmental variables. In this sense, the co-association networks allowed to both study the genetic and physical linkage between different co-associations modules but also to compare two different models of evolution: a Modular environmental response architecture (specific genes are affected by specific aspects of the environment) or a Universal pleiotropic environmental response architecture (all genes are affected by all aspects of the environment). The representation of different correlations between allelic frequency and environmental factors (named galaxy biplots) are especially informative to understand the effect of the different clusters on specific aspects of the environment (for example, the co-association network ‘Aridity’ shows strong associations with hot/wet versus cold/dry environments).

The analysis performed by Lotterhos et al. [3], although it has some unavoidable limitations (e.g., only extreme candidate SNPs are selected, limiting the results to the stronger effects; the genetic and physical map is incomplete in this species), includes relevant results and also implements new methodologies in the field. To highlight some of them: the preponderance of a Modular environmental response architecture (evolution in separated modules), the detection of physical linkage among SNPs that are co-associated with different aspects of the environment (which was unexpected a priori), the implementation of co-association networks and galaxy biplots to see the effect of modularity and pleiotropy on different aspects of environment. Finally, this work contains remarkable introductory Figures and Tables explaining unambiguously the main concepts [6] included in this study. This work can be treated as a starting point for many other future studies in the field.

References

[1] Hancock AM, Brachi B, Faure N, Horton MW, Jarymowycz LB, Sperone FG, Toomajian C, Roux F & Bergelson J. 2011. Adaptation to climate across the Arabidopsis thaliana genome. Science 334: 83–86. doi: 10.1126/science.1209244
[2] Wagner GP & Zhang J. The pleiotropic structure of the genotype­phenotype map: the evolvability of complex organisms. Nature Review Genetics 12: 204–213. doi: 10.1038/nrg2949
[3] Lotterhos KE, Yeaman S, Degner J, Aitken S, Hodgins K. 2018. Modularity of genes involved in local adaptation to climate despite physical linkage. bioRxiv 202481, ver. 4 peer-reviewed by Peer Community In Evolutionary Biology. doi: 10.1101/202481
[4] Lotterhos KE & Whitlock MC. 2014. Evaluation of demographic history and neutral parameterization on the performance of FST outlier tests. Molecular Ecology 23: 2178–2192. doi: 10.1111/mec.12725
[5] Lotterhos KE & Whitlock MC. 2015. The relative power of genome scans to detect local adaptation depends on sampling design and statistical method. Molecular Ecology 24: 1031–1046. doi: 10.1111/mec.13100
[6] Paaby AB & Rockman MV. 2013. The many faces of pleiotropy. Trends in Genetics 29: 66-73. doi: 10.1016/j.tig.2012.10.010