Genetic load: genomic estimates and applications in non-model animals.

Genetic variation, which is generated by mutation, recombination and gene flow, can reduce the mean fitness of a population, both now and in the future. This 'genetic load' has been estimated in a wide range of animal taxa using various approaches. Advances in genome sequencing and computational techniques now enable us to estimate the genetic load in populations and individuals without direct fitness estimates. Here, we review the classic and contemporary literature of genetic load. We describe approaches to quantify the genetic load in whole-genome sequence data based on evolutionary conservation and annotations. We show that splitting the load into its two components - the realized load (or expressed load) and the masked load (or inbreeding load) - can improve our understanding of the population genetics of deleterious mutations.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes.

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

Gene Expression Shifts in Emperor Penguin Adaptation to the Extreme Antarctic Environment.

Gene expression can accelerate ecological divergence by rapidly tweaking the response of an organism to novel environments, with more divergent environments exerting stronger selection and supposedly, requiring faster adaptive responses. Organisms adapted to extreme environments provide ideal systems to test this hypothesis, particularly when compared to related species with milder ecological niches. The Emperor penguin (Aptenodytes forsteri) is the only endothermic vertebrate breeding in the harsh Antarctic winter, in stark contrast with the less cold-adapted sister species, the King penguin (A. patagonicus). Assembling the first de novo transcriptomes and analysing multi-tissue (brain, kidney, liver, muscle, skin) RNA-Seq data from natural populations of both species, we quantified the shifts in tissue-enhanced genes, co-expression gene networks, and differentially expressed genes characterising Emperor penguin adaptation to the extreme Antarctic. Our analyses revealed the crucial role played by muscle and liver in temperature homeostasis, fasting, and whole-body energy metabolism (glucose/insulin regulation, lipid metabolism, fatty acid beta-oxidation, and blood coagulation). Repatterning at the regulatory level appears as more important in the brain of the Emperor penguin, showing the lowest signature of differential gene expression, but the largest co-expression gene network shift. Nevertheless, over-expressed genes related to mTOR signalling in the brain and the liver support their central role in cold and fasting responses. Besides contributing to understanding the genetics underlying complex traits, like body energy reservoir management, our results provide a first insight into the role of gene expression in adaptation to one of the most extreme environmental conditions endured by an endotherm.

Recurrent allopolyploidizations diversify ecophysiological traits in marsh orchids (Dactylorhiza majalis s.l.).

Whole-genome duplication has shaped the evolution of angiosperms and other organisms, and is important for many crops. Structural reorganization of chromosomes and repatterning of gene expression are frequently observed in allopolyploids, with physiological and ecological consequences. Recurrent origins from different parental populations are widespread among polyploids, resulting in an array of lineages that provide excellent models to uncover mechanisms of adaptation to divergent environments in early phases of polyploid evolution. We integrate here transcriptomic and ecophysiological comparative studies to show that sibling allopolyploid marsh orchid species (Dactylorhiza, Orchidaceae) occur in different habitats (low nutrient fens vs. meadows with mesic soils) and are characterized by a complex suite of intertwined, pronounced ecophysiological differences between them. We uncover distinct features in leaf elemental chemistry, light-harvesting, photoprotection, nutrient transport and stomata activity of the two sibling allopolyploids, which appear to match their specific ecologies, in particular soil chemistry differences at their native sites. We argue that the phenotypic divergence between the sibling allopolyploids has a clear genetic basis, generating ecological barriers that maintain distinct, independent lineages, despite pervasive interspecific gene flow. This suggests that recurrent origins of polyploids bring about a long-term potential to trigger and maintain functional and ecological diversity in marsh orchids and other groups.

Historical dispersal and host-switching formed the evolutionary history of a globally distributed multi-host parasite - The Ligula intestinalis species complex.

Studies on parasite biogeography and host spectrum provide insights into the processes driving parasite diversification. Global geographical distribution and a multi-host spectrum make the tapeworm Ligula intestinalis a promising model for studying both the vicariant and ecological modes of speciation in parasites. To understand the relative importance of host association and biogeography in the evolutionary history of this tapeworm, we analysed mtDNA and reduced-represented genomic SNP data for a total of 139 specimens collected from 18 fish-host genera across a distribution range representing 21 countries. Our results strongly supported the existence of at least 10 evolutionary lineages and estimated the deepest divergence at approximately 4.99-5.05 Mya, which is much younger than the diversification of the fish host genera and orders. Historical biogeography analyses revealed that the ancestor of the parasite diversified following multiple vicariance events and was widespread throughout the Palearctic, Afrotropical, and Nearctic between the late Miocene and early Pliocene. Cyprinoids were inferred as the ancestral hosts for the parasite. Later, from the late Pliocene to Pleistocene, new lineages emerged following a series of biogeographic dispersal and host-switching events. Although only a few of the current Ligula lineages show narrow host-specificity (to a single host genus), almost no host genera, even those that live in sympatry, overlapped between different Ligula lineages. Our analyses uncovered the impact of historical distribution shifts on host switching and the evolution of host specificity without parallel host-parasite co-speciation. Historical biogeography reconstructions also found that the parasite colonized several areas (Afrotropical and Australasian) much earlier than was suggested by only recent faunistic data.

Fitness consequences and ancestry loss in the Apennine brown bear after a simulated genetic rescue intervention.

Reduction in population size, with its predicted effects on population fitness, is the most alarming anthropogenic impact on endangered species. By introducing compatible individuals, genetic rescue (GR) is a promising but debated approach for reducing the genetic load unmasked by inbreeding and for restoring the fitness of declining populations. Although GR can improve genetic diversity and fitness, it can also produce loss of ancestry, hampering local adaptation, or replace with introduced variants the unique genetic pools evolved in endemic groups. We used forward genetic simulations based on empirical genomic data to assess fitness benefits and loss of ancestry risks of GR in the Apennine brown bear (Ursus arctos marsicanus). There are approximately 50 individuals of this isolated subspecies, and they have lower genetic diversity and higher inbreeding than other European brown bears, and GR has been suggested to reduce extinction risks. We compared 10 GR scenarios in which the number and genetic characteristics of migrants varied with a non-GR scenario of simple demographic increase due to nongenetic factors. The introduction of 5 individuals of higher fitness or lower levels of deleterious mutations than the target Apennine brown bear from a larger European brown bear population produced a rapid 10-20% increase in fitness in the subspecies and up to 22.4% loss of ancestry over 30 generations. Without a contemporary demographic increase, fitness started to decline again after a few generations. Doubling the population size without GR gradually increased fitness to a comparable level, but without losing ancestry, thus resulting in the best strategy for the Apennine brown bear conservation. Our results highlight the importance for management of endangered species of realistic forward simulations grounded in empirical whole-genome data.

Consecuencias en la aptitud y pérdida de ascendencia del oso pardo de los Apeninos después de un rescate genético simulado Resumen La reducción del tamaño poblacional, con los previsibles efectos sobre su aptitud, es el impacto antropogénico más alarmante sobre las especies amenazadas. Mediante la introducción de individuos compatibles, el rescate genético (RG) es una estrategia prometedora para reducir la carga genética revelada por la endogamia y restaurar la aptitud de las poblaciones en declive, aunque todavía se debate la eficiencia de esta. Aunque el RG puede mejorar la diversidad genética y la aptitud, también puede producir pérdida de ascendencia, lo que puede dificultar la adaptación local, o sustituir con variantes introducidas por los migrantes los acervos genéticos únicos que han evolucionado en grupos endémicos. En este trabajo realizamos simulaciones genéticas a futuro basadas en datos genómicos empíricos para evaluar los beneficios del RG en términos de aptitud y los riesgos de la pérdida de ascendencia en el oso pardo de los Apeninos (Ursus arctos marsicanus). Quedan aproximadamente 50 individuos de esta subespecie aislada que cuentan con una menor diversidad genética y un mayor nivel de endogamia comparado con otros osos pardos europeos y se ha sugerido que el RG podria reducir el riesgo de extinción de esta población. Comparamos 10 escenarios de RG en los que variaban el número y las características genéticas de los osos migrantes con un escenario sin RG con aumento demográfico causado por factores no genéticos. La introducción de 5 individuos procedentes de una población europea de oso pardo con mayor aptitud o niveles menores de mutaciones deletéreas que el oso pardo de los Apeninos produjo un rápido aumento de la aptitud del 10-20% en la subespecie y hasta un 22.4% de pérdida de ascendencia durante 30 generaciones. En las simulaciones sin un aumento demográfico, la aptitud empezó a disminuir de nuevo después de unas pocas generaciones. La duplicación del tamaño de la población sin RG aumentó gradualmente la aptitud hasta un nivel comparable al de algunos escenarios de RG, pero sin pérdida de ascendencia, por lo que parece ser la mejor estrategia para la conservación del oso pardo de los Apeninos. Nuestros resultados resaltan la importancia que tienen las simulaciones realistas a futuro basadas en datos empíricos del genoma completo para la gestión de especies amenazadas.

A high-quality reference genome for the critically endangered Aeolian wall lizard, Podarcis raffonei.

The Aeolian wall lizard, Podarcis raffonei, is an endangered species endemic to the Aeolian archipelago, Italy, where it is present only in 3 tiny islets and a narrow promontory of a larger island. Because of the extremely limited area of occupancy, severe population fragmentation and observed decline, it has been classified as Critically Endangered by the International Union for the Conservation of Nature (IUCN). Using Pacific Biosciences (PacBio) High Fidelity (HiFi) long-read sequencing, Bionano optical mapping and Arima chromatin conformation capture sequencing (Hi-C), we produced a high-quality, chromosome-scale reference genome for the Aeolian wall lizard, including Z and W sexual chromosomes. The final assembly spans 1.51 Gb across 28 scaffolds with a contig N50 of 61.4 Mb, a scaffold N50 of 93.6 Mb, and a BUSCO completeness score of 97.3%. This genome constitutes a valuable resource for the species to guide potential conservation efforts and more generally for the squamate reptiles that are underrepresented in terms of available high-quality genomic resources.

Impact of past climate warming on genomic diversity and demographic history of collared lemmings across the Eurasian Arctic.

The Arctic climate was warmer than today at the last interglacial and the Holocene thermal optimum. To reveal the impact of past climate-warming events on the demographic history of an Arctic specialist, we examined both mitochondrial and nuclear genomic variation in the collared lemming (Dicrostonyx torquatus, Pallas), a keystone species in tundra communities, across its entire distribution in northern Eurasia. The ancestral phylogenetic position of the West Beringian group and divergence time estimates support the hypothesis of continental range contraction to a single refugial area located in West Beringia during high-magnitude warming of the last interglacial, followed by westward recolonization of northern Eurasia in the last glacial period. The West Beringian group harbors the highest mitogenome diversity and its inferred demography indicates a constantly large effective population size over the Late Pleistocene to Holocene. This suggests that northward forest expansion during recent warming of the Holocene thermal optimum did not affect the gene pool of the collared lemming in West Beringia but reduced genomic diversity and effective population size in all other regions of the Eurasian Arctic. Demographic inference from genomic diversity was corroborated by species distribution modeling showing reduction in species distribution during past climate warming. These conclusions are supported by recent paleoecological evidence suggesting smaller temperature increases and moderate northward forest advances in the extreme northeast of Eurasia during the Late Pleistocene-to-Holocene warming events. This study emphasizes the importance of West Beringia as a potential refugium for cold-adapted Arctic species under ongoing climate warming.

Population Dynamics and Structural Effects at Short and Long Range Support the Hypothesis of the Selective Advantage of the G614 SARS-CoV-2 Spike Variant.

SARS-CoV-2 epidemics quickly propagated worldwide, sorting virus genomic variants in newly established propagules of infections. Stochasticity in transmission within and between countries or an actual selective advantage could explain the global high frequency reached by some genomic variants. Using statistical analyses, demographic reconstructions, and molecular dynamics simulations, we show that the globally invasive G614 spike variant 1) underwent a significant demographic expansion in most countries explained neither by stochastic effects nor by overrepresentation in clinical samples, 2) increases the spike S1/S2 furin-like site conformational plasticity (short-range effect), and 3) modifies the internal motion of the receptor-binding domain affecting its cross-connection with other functional domains (long-range effect). Our results support the hypothesis of a selective advantage at the basis of the spread of the G614 variant, which we suggest may be due to structural modification of the spike protein at the S1/S2 proteolytic site, and provide structural information to guide the design of variant-specific drugs.

Population genomics of a forest fungus reveals high gene flow and climate adaptation signatures.

Genome sequencing of spatially distributed individuals sheds light on how evolution structures genetic variation. Populations of Phellopilus nigrolimitatus, a red-listed wood-inhabiting fungus associated with old-growth coniferous forests, have decreased in size over the last century due to a loss of suitable habitats. We assessed the population genetic structure and investigated local adaptation in P. nigrolimitatus, by establishing a reference genome and genotyping 327 individuals sampled from 24 locations in Northern Europe by RAD sequencing. We revealed a shallow population genetic structure, indicating large historical population sizes and high levels of gene flow. Despite this weak substructuring, two genetic groups were recognized; a western group distributed mostly in Norway and an eastern group covering most of Finland, Poland and Russia. This substructuring may reflect coimmigration with the main host, Norway spruce (Picea abies), into Northern Europe after the last ice age. We found evidence of low levels of genetic diversity in southwestern Finland, which has a long history of intensive forestry and urbanization. Numerous loci were significantly associated with one or more environmental factors, indicating adaptation to specific environments. These loci clustered into two groups with different associations with temperature and precipitation. Overall, our findings indicate that the current population genetic structure of P. nigrolimitatus results from a combination of gene flow, genetic drift and selection. The acquisition of similar knowledge especially over broad geographic scales, linking signatures of adaptive genetic variation to evolutionary processes and environmental variation, for other fungal species will undoubtedly be useful for assessment of the combined effects of habitat fragmentation and climate change on fungi strongly bound to old-growth forests.

Conflicting nuclear and mitogenome phylogenies reveal ancient mitochondrial replacement between two North American species of collared lemmings (Dicrostonyx groenlandicus, D. hudsonius).

Collared lemmings (Dicrostonyx) are cold adapted rodents, keystone animals in the tundra communities and the model taxa in studies of Arctic genetic diversity and Quaternary paleontology. We examined mitochondrial and nuclear genomic variation to reconstruct phylogenetic relationships among the Eurasian D. torquatus and North American D. groenlandicus, D. hudsonius and evaluate biogeographic hypothesis of the two colonization events of North America from Eurasia based on morphological variation in dental traits. The nuclear and mitogenome phylogenies support reciprocal monophyly of each species but reveal conflicting relationships among species. The mitogenome tree likely reflects ancient mitochondrial replacement between currently isolated D. groenlandicus and D. hudsonius. The nuclear genome phylogeny reveals species cladogenesis and supports the hypothesis that D. hudsonius with primitive and distinct molar morphology represents a relic of the first migration event from Eurasia to North America. Species widely distributed in the North American Arctic, D. groenlandicus, with advanced dental morphology originated from a later colonization event across the Bering Land Bridge. This study shows ancient mitochondrial capture between two Arctic species and emphasizes the importance of multilocus approaches for phylogenetic inference.

Selection-driven adaptation to the extreme Antarctic environment in the Emperor penguin.

The eco-evolutionary history of penguins is characterised by shifting from temperate to cold environments. Breeding in Antarctica, the Emperor penguin appears as an extreme outcome of this process, with unique features related to insulation, heat production and energy management. However, whether this species actually diverged from a less cold-adapted ancestor, more ecologically similar to its sister species, the King penguin, is still an open question. As the Antarctic colonisation likely resulted in vast changes in selective pressure experienced by the Emperor penguin, the relative quantification of the genomic signatures of selection, unique to each sister species, could answer this question. Applying phylogeny-based selection tests on 7651 orthologous genes, we identified a more pervasive selection shift in the Emperor penguin than in the King penguin, supporting the hypothesis that its extreme cold adaptation is a derived state. Furthermore, among candidate genes under selection, four (TRPM8, LEPR, CRB1, and SFI1) were identified before in other cold-adapted homeotherms, like the woolly Mammoth, while other 161 genes can be assigned to biological functions relevant to cold adaptation identified in previous studies. Location and structural effects of TRPM8 substitutions in Emperor and King penguin lineages support their functional role with putative divergent effects on thermal adaptation. We conclude that extreme cold adaptation in the Emperor penguin largely involved unique genetic options which, however, affect metabolic and physiological traits common to other cold-adapted homeotherms.

Survival and divergence in a small group: The extraordinary genomic history of the endangered Apennine brown bear stragglers.

About 100 km east of Rome, in the central Apennine Mountains, a critically endangered population of ∼50 brown bears live in complete isolation. Mating outside this population is prevented by several 100 km of bear-free territories. We exploited this natural experiment to better understand the gene and genomic consequences of surviving at extremely small population size. We found that brown bear populations in Europe lost connectivity since Neolithic times, when farming communities expanded and forest burning was used for land clearance. In central Italy, this resulted in a 40-fold population decline. The overall genomic impact of this decline included the complete loss of variation in the mitochondrial genome and along long stretches of the nuclear genome. Several private and deleterious amino acid changes were fixed by random drift; predicted effects include energy deficit, muscle weakness, anomalies in cranial and skeletal development, and reduced aggressiveness. Despite this extreme loss of diversity, Apennine bear genomes show nonrandom peaks of high variation, possibly maintained by balancing selection, at genomic regions significantly enriched for genes associated with immune and olfactory systems. Challenging the paradigm of increased extinction risk in small populations, we suggest that random fixation of deleterious alleles (i) can be an important driver of divergence in isolation, (ii) can be tolerated when balancing selection prevents random loss of variation at important genes, and (iii) is followed by or results directly in favorable behavioral changes.

RADpainter and fineRADstructure: Population Inference from RADseq Data.

Powerful approaches to inferring recent or current population structure based on nearest neighbor haplotype "coancestry" have so far been inaccessible to users without high quality genome-wide haplotype data. With a boom in nonmodel organism genomics, there is a pressing need to bring these methods to communities without access to such data. Here, we present RADpainter, a new program designed to infer the coancestry matrix from restriction-site-associated DNA sequencing (RADseq) data. We combine this program together with a previously published MCMC clustering algorithm into fineRADstructure-a complete, easy to use, and fast population inference package for RADseq data (https://github.com/millanek/fineRADstructure; last accessed February 24, 2018). Finally, with two example data sets, we illustrate its use, benefits, and robustness to missing RAD alleles in double digest RAD sequencing.

Ecological plant epigenetics: Evidence from model and non-model species, and the way forward.

Growing evidence shows that epigenetic mechanisms contribute to complex traits, with implications across many fields of biology. In plant ecology, recent studies have attempted to merge ecological experiments with epigenetic analyses to elucidate the contribution of epigenetics to plant phenotypes, stress responses, adaptation to habitat, and range distributions. While there has been some progress in revealing the role of epigenetics in ecological processes, studies with non-model species have so far been limited to describing broad patterns based on anonymous markers of DNA methylation. In contrast, studies with model species have benefited from powerful genomic resources, which contribute to a more mechanistic understanding but have limited ecological realism. Understanding the significance of epigenetics for plant ecology requires increased transfer of knowledge and methods from model species research to genomes of evolutionarily divergent species, and examination of responses to complex natural environments at a more mechanistic level. This requires transforming genomics tools specifically for studying non-model species, which is challenging given the large and often polyploid genomes of plants. Collaboration among molecular geneticists, ecologists and bioinformaticians promises to enhance our understanding of the mutual links between genome function and ecological processes.

Genomic analyses suggest parallel ecological divergence in Heliosperma pusillum (Caryophyllaceae).

The mosaic distribution of interbreeding taxa with contrasting ecology and morphology offers an opportunity to study microevolutionary dynamics during ecological divergence. We investigate here the evolutionary history of an alpine and a montane ecotype of Heliosperma pusillum (Caryophyllaceae) in the south-eastern Alps. From six pairs of geographically close populations of the two ecotypes (120 individuals) we obtained a high-coverage restriction site associated DNA sequencing (RADseq) dataset that was used for demographic inference to test the hypothesis of parallel evolution of the two ecotypes. The data are consistent with repeated ecological divergence in H. pusillum, uncovering up to five polytopic origins of one ecotype from the other. A complex evolutionary history is evidenced, with local isolation-with-migration in two population pairs and intra-ecotype migration in two others. In all cases, the time of divergence or secondary contact was inferred as postglacial. A metagenomic analysis on exogenous contaminant RAD sequences suggests divergent microbial communities between the ecotypes. The lack of shared genomic regions of high divergence across population pairs illustrates the action of drift and/or local selection in shaping genetic divergence across repeated cases of ecological divergence.

Adaptive sequence evolution is driven by biotic stress in a pair of orchid species (Dactylorhiza) with distinct ecological optima.

The orchid family is the largest in the angiosperms, but little is known about the molecular basis of the significant variation they exhibit. We investigate here the transcriptomic divergence between two European terrestrial orchids, Dactylorhiza incarnata and Dactylorhiza fuchsii, and integrate these results in the context of their distinct ecologies that we also document. Clear signals of lineage-specific adaptive evolution of protein-coding sequences are identified, notably targeting elements of biotic defence, including both physical and chemical adaptations in the context of divergent pools of pathogens and herbivores. In turn, a substantial regulatory divergence between the two species appears linked to adaptation/acclimation to abiotic conditions. Several of the pathways affected by differential expression are also targeted by deviating post-transcriptional regulation via sRNAs. Finally, D. incarnata appears to suffer from insufficient sRNA control over the activity of RNA-dependent DNA polymerase, resulting in increased activity of class I transposable elements and, over time, in larger genome size than that of D. fuchsii. The extensive molecular divergence between the two species suggests significant genomic and transcriptomic shock in their hybrids and offers insights into the difficulty of coexistence at the homoploid level. Altogether, biological response to selection, accumulated during the history of these orchids, appears governed by their microenvironmental context, in which biotic and abiotic pressures act synergistically to shape transcriptome structure, expression and regulation.

Processes Driving the Adaptive Radiation of a Tropical Tree (Diospyros, Ebenaceae) in New Caledonia, a Biodiversity Hotspot.

Due to its special geological history, the New Caledonian Archipelago is a mosaic of soil types, and in combination with climatic conditions this results in a heterogeneous environment across relatively small distances. A group of over 20 endemic species of Diospyros (Ebenaceae) has rapidly and recently radiated on the archipelago after a single long-distance dispersal event. Most of the Diospyros species in the radiating group are morphologically and ecologically well differentiated, but they exhibit low levels of DNA variability. To investigate the processes that shaped the diversification of this group we employed restriction site associated DNA sequencing (RADseq). Over 8400 filtered SNPs generally confirm species delimitations and produce a well-supported phylogenetic tree. Our analyses document local introgression, but only a limited potential for gene flow over longer distances. The phylogenetic relationships point to an early regional clustering among populations and species, indicating that allopatric speciation with respect to macrohabitat (i.e., climatic conditions) may have had a role in the initial differentiation within the group. A later, more rapid radiation involved divergence with respect to microhabitat (i.e., soil preference). Several sister species in the group show a parallel divergence in edaphic preference. Searches for genomic regions that are systematically differentiated in this replicated phenotypic divergence pointed to loci potentially involved in ion binding and cellular transport. These loci appear meaningful in the context of adaptations to soil types that differ in heavy-metal and mineral content. Identical nucleotide changes affected only two of these loci, indicating that introgression may have played a limited role in their evolution. Our results suggest that both allopatric diversification and (parapatric) ecological divergence shaped successive rounds of speciation in the Diospyros radiation on New Caledonia.

Testing Classical Species Properties with Contemporary Data: How "Bad Species" in the Brassy Ringlets (Erebia tyndarus complex, Lepidoptera) Turned Good.

All species concepts are rooted in reproductive, and ultimately genealogical, relations. Genetic data are thus the most important source of information for species delimitation. Current ease of access to genomic data and recent computational advances are blooming a plethora of coalescent-based species delimitation methods. Despite their utility as objective approaches to identify species boundaries, coalescent-based methods (1) rely on simplified demographic models that may fail to capture some attributes of biological species, (2) do not make explicit use of the geographic information contained in the data, and (3) are often computationally intensive. In this article, we present a case of species delimitation in the Erebia tyndarus species complex, a taxon regarded as a classic example of problematic taxonomic resolution. Our approach to species delimitation used genomic data to test predictions rooted in the biological species concept and in the criterion of coexistence in sympatry. We (1) obtained restriction-site associated DNA (RAD) sequencing data from a carefully designed sample, (2) applied two genotype clustering algorithms to identify genetic clusters, and (3) performed within-clusters and between-clusters analyses of isolation by distance as a test for intrinsic reproductive barriers. Comparison of our results with those from a Bayes factor delimitation coalescent-based analysis, showed that coalescent-based approaches may lead to overconfident splitting of allopatric populations, and indicated that incorrect species delimitation is likely to be inferred when an incomplete geographic sample is analyzed. While we acknowledge the theoretical justification and practical usefulness of coalescent-based species delimitation methods, our results stress that, even in the phylogenomic era, the toolkit for species delimitation should not dismiss more traditional, biologically grounded, approaches coupling genomic data with geographic information.

Divergent evolution and niche differentiation within the common peatmoss Sphagnum magellanicum.

PREMISE OF THE STUDY: Populations with phenotypic polymorphism in discrete characters may be good models for investigating genome evolution and speciation. Sphagnum magellanicum Brid. is found throughout the northern hemisphere, and despite considerable variation in morphological characters, it is considered one of the least taxonomically controversial peatmoss species. We have observed two main morphs of the species associated with different microhabitats. Here we investigated the genomic and environmental basis of this intraspecific morphological variation. METHODS: We conducted transplant and common garden experiments to test whether the two morphs are genetically differentiated. We then used RAD-sequencing to quantify the genomic divergence between the morphs and approximate Bayesian computation (ABC) to infer the most likely demographic scenario explaining the genome-wide differentiation of the two morphs. KEY RESULTS: We found that genomic differentiation between the two morphs is unexpectedly high and that several of the differentiated morphological characters have a genetic basis. Using simulation approaches, we found support for a scenario of ancient divergence followed by recent secondary contact. CONCLUSIONS: We show that the two morphs represent the two main genetic clusters previously found worldwide. Our results demonstrate that relatively minor morphological differentiation in a presumed phenotypically plastic peatmoss may be associated with massive divergence across the genome.