Genetic diversity in temperate and boreal forests tree species has been strongly affected by late Pleistocene climate oscillations [2,3,5], but also by anthropogenic forces. Particularly in Europe, where a long history of human intervention has re-distributed species and populations, it can be difficult to know if a given forest arose through natural regeneration and gene flow or through some combination of natural and human-mediated processes. This uncertainty can confound inferences of the causes and consequences of standing genetic variation, which may impact our interpretation of demographic events that shaped species before humans became dominant on the landscape. In their paper entitled "Genomic data provides new insights on the demographic history and the extent of recent material transfers in Norway spruce", Chen et al. [1] used a genome-wide dataset of 400k SNPs to infer the demographic history of Picea abies (Norway spruce), the most widespread and abundant spruce species in Europe, and to understand its evolutionary relationship with two other spruces (Picea obovata [Siberian spruce] and P. omorika [Serbian spruce]). Three major Norway spruce clusters were identified, corresponding to central Europe, Russia and the Baltics, and Scandinavia, which agrees with previous studies. The density of the SNP data in the present paper enabled inference of previously uncharacterized admixture between these groups, which corresponds to the timing of postglacial recolonization following the last glacial maximum (LGM). This suggests that multiple migration routes gave rise to the extant distribution of the species, and may explain why Chen et al.'s estimates of divergence times among these major Norway spruce groups were older (15mya) than those of previous studies (5-6mya) – those previous studies may have unknowingly included admixed material [4]. Treemix analysis also revealed extensive admixture between Norway and Siberian spruce over the last ~100k years, while the geographically-restricted Serbian spruce was both isolated from introgression and had a dramatically smaller effective population size (Ne) than either of the other two species. This small Ne resulted from a bottleneck associated with the onset of the iron age ~3000 years ago, which suggests that anthropogenic depletion of forest resources has severely impacted this species. Finally, ancestry of Norway spruce samples collected in Sweden and Denmark suggest their recent transfer from more southern areas of the species range. This northward movement of genotypes likely occurred because the trees performed well relative to local provenances, which is a common observation when trees from the south are planted in more northern locations (although at the potential cost of frost damage due to inappropriate phenology). While not the reason for the transfer, the incorporation of southern seed sources into the Swedish breeding and reforestation program may lead to more resilient forests under climate change. Taken together, the data and analysis presented in this paper allowed inference of the intra- and interspecific demographic histories of a tree species group at a very high resolution, and suggest caveats regarding sampling and interpretation of data from areas with a long history of occupancy by humans.

References

[1] Chen, J., Milesi, P., Jansson, G., Berlin, M., Karlsson, B., Aleksić, J. M., Vendramin, G. G., Lascoux, M. (2018). Genomic data provides new insights on the demographic history and the extent of recent material transfers in Norway spruce. BioRxiv, 402016. ver. 3 peer-reviewed and recommended by PCI Evol Biol. doi: 10.1101/402016
[2] Holliday, J. A., Yuen, M., Ritland, K., & Aitken, S. N. (2010). Postglacial history of a widespread conifer produces inverse clines in selective neutrality tests. Molecular Ecology, 19(18), 3857–3864. doi: 10.1111/j.1365-294X.2010.04767.x
[3] Ingvarsson, P. K. (2008). Multilocus patterns of nucleotide polymorphism and the demographic history of Populus tremula. Genetics, 180, 329-340. doi: 10.1534/genetics.108.090431
[4] Lockwood, J. D., Aleksić, J. M., Zou, J., Wang, J., Liu, J., & Renner, S. S. (2013). A new phylogeny for the genus Picea from plastid, mitochondrial, and nuclear sequences. Molecular Phylogenetics and Evolution, 69(3), 717–727. doi: 10.1016/j.ympev.2013.07.004
[5] Pyhäjärvi, T., Garcia-Gil, M. R., Knürr, T., Mikkonen, M., Wachowiak, W., & Savolainen, O. (2007). Demographic history has influenced nucleotide diversity in European Pinus sylvestris populations. Genetics, 177(3), 1713–1724. doi: 10.1534/genetics.107.077099 "

The genomic revolution has given us access to inexpensive genetic data for any species. Simultaneously we have lost the ability to easily identify chimerism in samples or some unusual deviations from standard Mendelian genetics. Methods have been developed to identify sex chromosomes, characterise the ploidy, or understand the exact form of parthenogenesis from genomic data. However, we rarely consider that the tissues we extract DNA from could be a mixture of cells with different genotypes or karyotypes. This can nonetheless happen for a variety of (fascinating) reasons such as somatic chromosome elimination, transmissible cancer, or parental genome elimination. Without a dedicated analysis, it is very easy to miss it.

In this preprint, Jaron et al. (2022) used an ingenious analysis of whole individual NGS data to test the hypothesis of paternal genome elimination in the globular springtail Allacma fusca. The authors suspected that a high fraction of the whole body of males is made of sperm in this species and if this species undergoes paternal genome elimination, we would expect that sperm would only contain maternally inherited chromosomes. Given the reference genome was highly fragmented, they developed a two-tissue model to analyse Kmer spectra and obtained confirmation that around one-third of the tissue was sperm in males. This allowed them to test whether coverage patterns were consistent with the species exhibiting paternal genome elimination. They combined their estimation of the fraction of haploid tissue with allele coverages in autosomes and the X chromosome to obtain support for a bias toward one parental allele, suggesting that all sperm carries the same parental haplotype. It could be the maternal or the paternal alleles, but paternal genome elimination is most compatible with the known biology of Arthropods. SNP calling was used to confirm conclusions based on the analysis of the raw pileups.

I found this study to be a good example of how a clever analysis of Kmer spectra and allele coverages can provide information about unusual modes of reproduction in a species, even though it does not have a well-assembled genome yet. As advocated by the authors, routine inspection of Kmer spectra and allelic read-count distributions should be included in the best practice of NGS data analysis. They provide the method to identify paternal genome elimination but also the way to develop similar methods to detect another kind of genetic chimerism in the avalanche of sequence data produced nowadays.

References

Jaron KS, Hodson CN, Ellers J, Baird SJ, Ross L (2022) Genomic evidence of paternal genome elimination in the globular springtail Allacma fusca. bioRxiv, 2021.11.12.468426, ver. 5 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2021.11.12.468426

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

Modelling the evolution of complete genome sequences in populations requires accounting for the recombination process, as a single tree can no longer describe the underlying genealogy. The sequentially Markov coalescent (SMC, McVean and Cardin 2005; Marjoram and Wall 2006) approximates the standard coalescent with recombination process and permits estimating population genetic parameters (e.g., population sizes, recombination rates) using population genomic datasets. As such datasets become available for an increasing number of species, more fine-tuned models are needed to encompass the diversity of life cycles of organisms beyond the model species on which most methods have been benchmarked.

The work by Korfmann et al. (Korfmann et al. 2024) represents a significant step forward as it accounts for multiple mergers in SMC models. Multiple merger models account for simultaneous coalescence events so that more than two lineages find a common ancestor in a given generation. This feature is not allowed in standard coalescent models and may result from selection or skewed offspring distributions, conditions likely met by a broad range of species, particularly microbial.

Yet, this work goes beyond extending the SMC, as it introduces several methodological innovations. The "classical" SMC-based inference approaches rely on hidden Markov models to compute the likelihood of the data while efficiently integrating over the possible ancestral recombination graphs (ARG). Following other recent works (e.g. Gattepaille et al. 2016), Korfmann et al. propose to separate the ARG inference from model parameter estimation under maximum likelihood (ML). They introduce a procedure where the ARG is first reconstructed from the data and then taken as input in the model fitting step. While this approach does not permit accounting for the uncertainty in the ARG reconstruction (which is typically large), it potentially allows for the extraction of more information from the ARG, such as the occurrence of multiple merging events. Going away from maximum likelihood inference, the authors trained a graph neural network (GNN) on simulated ARGs, introducing a new, flexible way to estimate population genomic parameters.

The authors used simulations under a beta-coalescent model with diverse demographic scenarios and showed that the ML and GNN approaches introduced can reliably recover the simulated parameter values. They further show that when the true ARG is given as input, the GNN outperforms the ML approach, demonstrating its promising power as ARG reconstruction methods improve. In particular, they showed that trained GNNs can disentangle the effects of selective sweeps and skewed offspring distributions while inferring past population size changes.

This work paves the way for new, exciting applications, though many questions must be answered. How frequent are multiple mergers? As the authors showed that these events "erase" the record of past demographic events, how many genomes are needed to conduct reliable inference, and can the methods computationally cope with the resulting (potentially large) amounts of required data? This is particularly intriguing as micro-organisms, prone to strong selection and skewed offspring distributions, also tend to carry smaller genomes.

References

Gattepaille L, Günther T, Jakobsson M. 2016. Inferring Past Effective Population Size from Distributions of Coalescent Times. Genetics 204:1191-1206.
https://doi.org/10.1534/genetics.115.185058
 
Korfmann K, Sellinger T, Freund F, Fumagalli M, Tellier A. 2024. Simultaneous Inference of Past Demography and Selection from the Ancestral Recombination Graph under the Beta Coalescent. bioRxiv, 2022.09.28.508873. ver. 5 peer-reviewed and recommended by Peer Community in Evolutionary Biology. https://doi.org/10.1101/2022.09.28.508873
 
Marjoram P, Wall JD. 2006. Fast "coalescent" simulation. BMC Genet. 7:16.
https://doi.org/10.1186/1471-2156-7-16
 
McVean GAT, Cardin NJ. 2005. Approximating the coalescent with recombination. Philos. Trans. R. Soc. Lond. B. Biol. Sci. 360:1387-1393.
https://doi.org/10.1098/rstb.2005.1673

Two back-to-back papers published earlier this year in Nature Genetics provide compelling evidence for the control of a male reproductive polymorphism in a wading bird by a "supergene", a cluster of tightly linked genes [1-2]. The bird in question, the ruff (Philomachus pugnax), has a rather unusual reproductive system that consists of three distinct types of males ("reproductive morphs"): aggressive "independents" who represent the majority of males; a smaller fraction of non-territorial "satellites" who are submissive towards "independents"; and "faeders" who mimic females and are rare. Previous work has shown that the male morphs differ in major aspects of mating and aggression behavior, plumage coloration and body size, and that – intriguingly – this complex multi-trait polymorphism is apparently controlled by a single autosomal Mendelian locus with three alleles [3]. To uncover the genetic control of this polymorphism two independent teams, led by Terry Burke [1] and Leif Andersson [2], have set out to analyze the genomes of male ruffs. Using a combination of genomics and genetics, both groups managed to pin down the supergene locus and map it to a non-recombining, 4.5 Mb large inversion which arose 3.8 million years ago. While "independents" are homozygous for the ancestral uninverted sequence, "satellites" and "faeders" carry evolutionarily divergent, dominant alternative haplotypes of the inversion. Thus, as in several other notable cases, for example the supergene control of disassortative mating, aggressiveness and plumage color in white-throated sparrows [4], of mimicry in Heliconius and Papilio butterflies [5-6], or of social structure in ants [7], an inversion – behaving as a single "locus" – underpins the mechanistic basis of the supergene. More generally, and beyond inversions, a growing number of studies now shows that selection can favor the evolution of suppressed recombination, thereby leading to the emergence of clusters of tightly linked loci which can then control – presumably due to polygenic gene action – a suite of complex phenotypes [8-10]. A largely unresolved question in this field concerns the identity of the causative alleles and loci within a given supergene. Recent progress on this question has been made for example in Papilio polytes butterflies where a mimicry supergene has been found to involve – surprisingly – only a single but large gene: multiple mimicry alleles in the doublesex gene are maintained in strong linkage disequilibrium via an inversion. It will clearly be of great interest to see future examples of such a fine-scale genetic dissection of supergenes. In conclusion, we were impressed by the data and analyses of Küpper et al. [1] and Lamichhaney et al. [2]: both papers beautifully illustrate how genomics and evolutionary ecology can be combined to make new, exciting discoveries. Both papers will appeal to readers with an interest in supergenes, inversions, the interplay of selection and recombination, or the genetic control of complex phenotypes.

References

[1] Küpper C, Stocks M, Risse JE, dos Remedios N, Farrell LL, McRae SB, Morgan TC, Karlionova N, Pinchuk P, Verkuil YI, et al. 2016. A supergene determines highly divergent male reproductive morphs in the ruff. Nature Genetics 48:79-83. doi: 10.1038/ng.3443

[2] Lamichhaney S, Fan G, Widemo F, Gunnarsson U, Thalmann DS, Hoeppner MP, Kerje S, Gustafson U, Shi C, Zhang H, et al. 2016. Structural genomic changes underlie alternative reproductive strategies in the ruff (Philomachus pugnax). Nature Genetics 48:84-88. doi: 10.1038/ng.3430

[3] Lank DB, Smith CM, Hanotte O, Burke T, Cooke F. 1995. Genetic polymorphism for alternative mating behaviour in lekking male ruff Philomachus pugnax. Nature 378:59-62. doi: 10.1038/378059a0

[4] Tuttle Elaina M, Bergland Alan O, Korody Marisa L, Brewer Michael S, Newhouse Daniel J, Minx P, Stager M, Betuel A, Cheviron Zachary A, Warren Wesley C, et al. 2016. Divergence and Functional Degradation of a Sex Chromosome-like Supergene. Current Biology 26:344-350. doi: 10.1016/j.cub.2015.11.069

[5] Joron M, Frezal L, Jones RT, Chamberlain NL, Lee SF, Haag CR, Whibley A, Becuwe M, Baxter SW, Ferguson L, et al. 2011. Chromosomal rearrangements maintain a polymorphic supergene controlling butterfly mimicry. Nature 477:203-206. doi: 10.1038/nature10341

[6] Kunte K, Zhang W, Tenger-Trolander A, Palmer DH, Martin A, Reed RD, Mullen SP, Kronforst MR. 2014. doublesex is a mimicry supergene. Nature 507:229-232. doi: 10.1038/nature13112

[7] Wang J, Wurm Y, Nipitwattanaphon M, Riba-Grognuz O, Huang Y-C, Shoemaker D, Keller L. 2013. A Y-like social chromosome causes alternative colony organization in fire ants. Nature 493:664-668. doi: 10.1038/nature11832

[8] Thompson MJ, Jiggins CD. 2014. Supergenes and their role in evolution. Heredity 113:1-8. doi: 10.1038/hdy.2014.20

[9] Schwander T, Libbrecht R, Keller L. 2014. Supergenes and Complex Phenotypes. Current Biology 24:R288-R294. doi: 10.1016/j.cub.2014.01.056

[10] Charlesworth D. 2015. The status of supergenes in the 21st century: recombination suppression in Batesian mimicry and sex chromosomes and other complex adaptations. Evolutionary Applications 9:74-90. doi: 10.1111/eva.12291