Rapid evolution of post-zygotic reproductive isolation is widespread in Arctic plant lineages.

BACKGROUND AND AIMS: The Arctic tundra, with its extreme temperatures and short growing season, is evolutionarily young and harbours one of the most species-poor floras on Earth. Arctic species often show little phenotypic and genetic divergence across circumpolar ranges. However, strong intraspecific post-zygotic reproductive isolation (RI) in terms of hybrid sterility has frequently evolved within selfing Arctic species of the genus Draba. Here we assess whether incipient biological species are common in the Arctic flora. METHODS: We conducted an extensive crossing experiment including six species representing four phylogenetically distant families collected across the circumpolar Arctic. We crossed conspecific parental populations representing different spatial scales, raised 740 F1 hybrids to maturity and measured fertility under laboratory conditions. We examined genetic divergence between populations for two of these species (Cardamine bellidifolia and Ranunculus pygmaeus). KEY RESULTS: In five of the six species, we find extensive reduction in pollen fertility and seed set in F1 hybrids; 219 (46 %) of the 477 F1 hybrids generated between parents separated by ≥427 km had <20 % pollen fertility. Isolation with migration (IM) and *BEAST analyses of sequences of eight nuclear genes in C. bellidifolia suggests that reproductively isolated populations of this species diverged during, or even after, the last glaciation. Likewise, Arctic populations of R. pygmaeus were genetically very similar despite exhibiting strongly reduced fertility in crosses, suggesting that RI evolved recently also in this species. CONCLUSION: We show that post-zygotic RI has developed multiple times within taxonomically recognized Arctic species belonging to several distantly related lineages, and that RI may have developed over just a few millennia. Rapid and widespread evolution of incipient biological species in the Arctic flora might be associated with frequent bottlenecks due to glacial cycles, and/or selfing mating systems, which are common in the harsh Arctic environment where pollinators are scarce.

Fifty thousand years of Arctic vegetation and megafaunal diet.

Although it is generally agreed that the Arctic flora is among the youngest and least diverse on Earth, the processes that shaped it are poorly understood. Here we present 50 thousand years (kyr) of Arctic vegetation history, derived from the first large-scale ancient DNA metabarcoding study of circumpolar plant diversity. For this interval we also explore nematode diversity as a proxy for modelling vegetation cover and soil quality, and diets of herbivorous megafaunal mammals, many of which became extinct around 10 kyr bp (before present). For much of the period investigated, Arctic vegetation consisted of dry steppe-tundra dominated by forbs (non-graminoid herbaceous vascular plants). During the Last Glacial Maximum (25-15 kyr bp), diversity declined markedly, although forbs remained dominant. Much changed after 10 kyr bp, with the appearance of moist tundra dominated by woody plants and graminoids. Our analyses indicate that both graminoids and forbs would have featured in megafaunal diets. As such, our findings question the predominance of a Late Quaternary graminoid-dominated Arctic mammoth steppe.

Vicariance, long-distance dispersal, and regional extinction-recolonization dynamics explain the disjunct circumpolar distribution of the arctic-alpine plant Silene acaulis.

PREMISE OF THE STUDY: Many arctic-alpine species have vast geographic ranges, but these may encompass substantial gaps whose origins are poorly understood. Here we address the phylogeographic history of Silene acaulis, a perennial cushion plant with a circumpolar distribution except for a large gap in Siberia. METHODS: We assessed genetic variation in a range-wide sample of 103 populations using plastid DNA (pDNA) sequences and AFLPs (amplified fragment length polymorphisms). We constructed a haplotype network and performed Bayesian phylogenetic analyses based on plastid sequences. We visualized AFLP patterns using principal coordinate analysis, identified genetic groups using the program structure, and estimated genetic diversity and rarity indices by geographic region. KEY RESULTS: The history of the main pDNA lineages was estimated to span several glaciations. AFLP data revealed a distinct division between Beringia/North America and Europe/East Greenland. These two regions shared only one of 17 pDNA haplotypes. Populations on opposite sides of the Siberian range gap (Ural Mountains and Chukotka) were genetically distinct and appear to have resulted from postglacial leading-edge colonizations. We inferred two refugia in North America (Beringia and the southern Rocky Mountains) and two in Europe (central-southern Europe and northern Europe/East Greenland). Patterns in the East Atlantic region suggested transoceanic long-distance dispersal events. CONCLUSIONS: Silene acaulis has a highly dynamic history characterized by vicariance, regional extinction, and recolonization, with persistence in at least four refugia. Long-distance dispersal explains patterns across the Atlantic Ocean, but we found no evidence of dispersal across the Siberian range gap.

Fungal palaeodiversity revealed using high-throughput metabarcoding of ancient DNA from arctic permafrost.

The taxonomic and ecological diversity of ancient fungal communities was assessed by combining next generation sequencing and metabarcoding of DNA preserved in permafrost. Twenty-six sediment samples dated 16 000-32 000 radiocarbon years old from two localities in Siberia were analysed for fungal ITS. We detected 75 fungal OTUs from 21 orders representing three phyla, although rarefaction analyses suggested that the full diversity was not recovered despite generating an average of 6677 ± 3811 (mean ± SD) sequences per sample and that preservation bias likely has considerable effect on the recovered DNA. Most OTUs (75.4%) represented ascomycetes. Due to insufficient sequencing depth, DNA degradation and putative preservation biases in our samples, the recovered taxa probably do not represent the complete historic fungal community, and it is difficult to determine whether the fungal communities varied geographically or experienced a composition shift within the period of 16 000-32 000 bp. However, annotation of OTUs to functional ecological groups provided a wealth of information on the historic communities. About one-third of the OTUs are presumed plant-associates (pathogens, saprotrophs and endophytes) typical of graminoid- and forb-rich habitats. We also detected putative insect pathogens, coprophiles and keratinophiles likely associated with ancient insect and herbivore faunas. The detection of putative insect pathogens, mycoparasites, aquatic fungi and endophytes broadens our previous knowledge of the diversity of fungi present in Beringian palaeoecosystems. A large group of putatively psychrophilic/psychrotolerant fungi was also detected, most likely representing a modern, metabolically active fungal community.

Pervasive Phylogenomic Incongruence Underlies Evolutionary Relationships in Eyebrights (Euphrasia, Orobanchaceae).

Disentangling the phylogenetic relationships of taxonomically complex plant groups is often mired by challenges associated with recent speciation, hybridization, complex mating systems, and polyploidy. Here, we perform a phylogenomic analysis of eyebrights (Euphrasia), a group renowned for taxonomic complexity, with the aim of documenting the extent of phylogenetic discordance at both deep and at shallow phylogenetic scales. We generate whole-genome sequencing data and integrate this with prior genomic data to perform a comprehensive analysis of nuclear genomic, nuclear ribosomal (nrDNA), and complete plastid genomes from 57 individuals representing 36 Euphrasia species. The species tree analysis of 3,454 conserved nuclear scaffolds (46 Mb) reveals that at shallow phylogenetic scales postglacial colonization of North Western Europe occurred in multiple waves from discrete source populations, with most species not being monophyletic, and instead combining genomic variants from across clades. At a deeper phylogenetic scale, the Euphrasia phylogeny is structured by geography and ploidy, and partially by taxonomy. Comparative analyses show Southern Hemisphere tetraploids include a distinct subgenome indicative of independent polyploidy events from Northern Hemisphere taxa. In contrast to the nuclear genome analyses, the plastid genome phylogeny reveals limited geographic structure, while the nrDNA phylogeny is informative of some geographic and taxonomic affinities but more thorough phylogenetic inference is impeded by the retention of ancestral polymorphisms in the polyploids. Overall our results reveal extensive phylogenetic discordance at both deeper and shallower nodes, with broad-scale geographic structure of genomic variation but a lack of definitive taxonomic signal. This suggests that Euphrasia species either have polytopic origins or are maintained by narrow genomic regions in the face of extensive homogenizing gene flow. Moreover, these results suggest genome skimming will not be an effective extended barcode to identify species in groups such as Euphrasia, or many other postglacial species groups.

Elucidating the evolutionary history of the Southeast Asian, holoparasitic, giant-flowered Rafflesiaceae: pliocene vicariance, morphological convergence and character displacement.

The aim of the present study is to elucidate the evolutionary history of the enigmatic holoparasitic Rafflesiaceae. More specifically, floral morphological evolution is interpreted in a molecular phylogenetic context, the biogeographic history of the family is investigated, and the possibility of character displacement to have been operating in this family is assessed. Parsimony and Bayesian methods are used to estimate phylogeny and divergence times among Rafflesiaceae species based on nuclear and mitochondrial DNA sequence data from Barkman et al. (2008) as well as new sequence data from additional samples and an additional genetic marker, the plastid 16S. Ancestral areas are inferred using dispersal-vicariance analysis (DIVA) as well a more recently developed parametric likelihood method (LAGRANGE), now including an update that allows for estimation over the posterior distribution of dated trees. Our extended molecular phylogeny of Rafflesiaceae implies a general lack of morphological synapomorphies as well as a high level of morphological homoplasy. In particular, a high level of floral morphological homoplasy is detected among Rafflesia species suggestive of similar patterns of pollinator-based selection in different geographic areas, and multiple instances of divergent floral size evolution is consistent with a model of character displacement. Initial diversification of Rafflesiaceae during the Late Cretaceous was followed by a long period of no-net diversification, likely due to extinctions caused by a Late Eocene to Miocene dramatic reduction in rainforest cover. A Late Miocene to Early Pliocene rise in sea-level probably caused the vicariant diversification observed between areas of endemism. The most recent species divergences are concordant with Pleistocene changes in climate and sea-levels, but apparently with no successful inter-area migrations, supportive of savannah, rather than rainforest, covered landbridges. An explosive increase in net diversification rate, most pronounced in Rafflesia, may be explained by Mid-Miocene to Pliocene rainforest-favorable conditions as well as natural selection promoting character displacement for Rafflesia flower size.

Afro-alpine flagships revisited: Parallel adaptation, intermountain admixture and shallow genetic structuring in the giant senecios (Dendrosenecio).

Distantly related lineages of the enigmatic giant rosette plants of tropical alpine environments provide classical examples of convergent adaptation. For the giant senecios (Dendrosenecio), the endemic landmarks of the East African sky islands, it has also been suggested that parallel adaptation has been important for within-lineage differentiation. To test this hypothesis and to address potential gene flow and hybridization among the isolated sky islands, we organized field expeditions to all major mountains. We sampled all currently accepted species and all but one subspecies and genotyped 460 plants representing 109 populations. We tested whether genetic structuring corresponds to geography, as predicted by a parallel adaptation hypothesis, or to altitudinal belt and habitat rather than mountains, as predicted by a hypothesis of a single origin of adaptations. Bayesian and Neighbor-Net analyses showed that the main genetic structure is shallow and largely corresponds to geography, supporting a hypothesis of recent, rapid radiation via parallel altitude/habitat adaptation on different mountains. We also found evidence for intermountain admixture, suggesting several long-distance dispersals by wind across vast areas of unsuitable habitat. The combination of parallel adaptation, secondary contact, and hybridization may explain the complex patterns of morphological variation and the contradicting taxonomic treatments of these rare enigmatic giants, supporting the use of wide taxonomic concepts. Notably, the within-population genetic diversity was very low and calls for increased conservation efforts.

Molecular phylogeny and biogeography of the bipolar Euphrasia (Orobanchaceae): recent radiations in an old genus.

Euphrasia includes perennial or annual green parasitic plants, and has a disjunct bipolar distribution except for one transtropical connection across the high mountains of Oceania. The disjunction is coupled with strikingly contrasting patterns of morphological diversity between the southern and northern hemispheres, making it an exciting model to study processes of evolutionary diversification which shaped present floras. We inferred the relationships among 51 species representing 14 of the 15 sections of the genus based on nrDNA ITS and cpDNA trnL intron, trnL-trnF and atpB-rbcL intergenic spacers. Maximum parsimony and Bayesian inference support monophyly of the genus and of several intrageneric groups characterized by morphology, ploidy level, and geographic range. Molecular phylogenetic dating using Bayesian "relaxed" clock methods suggests that the earliest Euphrasia radiations occurred minimum 11-8 Mya with bipolarity being achieved 7-5 Mya. Biogeographic analyses using explicit model-based approach inferred Eurasia as an ancestral area for the genus. The most parsimonious reconstruction found by a dispersal-vicariance analysis requires 17 dispersals to account for the current biogeographic pattern and supports Eurasian origin for Euphrasia. Both long-distance dispersal and across land vicariance can be invoked to explain the diversification in the genus, which experienced rapid radiations driven by new ecological opportunities of the late Pliocene and Pleistocene but also retained a set of local endemic or relict species of an earlier origin.

DNA barcoding a taxonomically complex hemiparasitic genus reveals deep divergence between ploidy levels but lack of species-level resolution.

DNA barcoding is emerging as a useful tool not only for species identification but also for studying evolutionary and ecological processes. Although plant DNA barcodes do not always provide species-level resolution, the generation of large DNA barcode data sets can provide insights into the mechanisms underlying the generation of species diversity. Here, we study evolutionary processes in taxonomically complex British Euphrasia (Orobanchaceae), a group with multiple ploidy levels, frequent self-fertilization, young species divergence and widespread hybridization. We use a phylogenetic approach to investigate the colonization history of British Euphrasia, followed by a DNA barcoding survey and population genetic analyses to reveal the causes of shared sequence variation. Phylogenetic analysis shows Euphrasia have colonized Britain from mainland Europe on multiple occasions. DNA barcoding reveals that no British Euphrasia species has a consistent diagnostic sequence profile, and instead, plastid haplotypes are either widespread across species, or are population specific. The partitioning of nuclear genetic variation suggests differences in ploidy act as a barrier to gene exchange, while the divergence between diploid and tetraploid ITS sequences supports the polyploids being allotetraploid in origin. Overall, these results show that even when lacking species-level resolution, analyses of DNA barcoding data can reveal evolutionary patterns in taxonomically complex genera.

DNA Metabarcoding Reveals Diet Overlap between the Endangered Walia Ibex and Domestic Goats - Implications for Conservation.

Human population expansion and associated degradation of the habitat of many wildlife species cause loss of biodiversity and species extinctions. The small Simen Mountains National Park in Ethiopia is one of the last strongholds for the preservation of a number of afro-alpine mammals, plants and birds, and it is home to the rare endemic Walia ibex, Capra walie. The narrow distribution range of this species as well as potential competition for resources with livestock, especially with domestic goat, Capra hircus, may compromise its future survival. Based on a curated afro-alpine taxonomic reference library constructed for plant taxon identification, we investigated the diet of the Walia ibex and addressed the dietary overlap with domestic goat using DNA metabarcoding of faecal samples. Faeces of both species were collected from different localities in the National Park. We show that both species are browsers, with forbs, shrubs and trees comprising the largest proportion of their diet, supplemented by grasses. There was a considerable overlap in dietary preferences. Several of the preferred diet items of the Walia ibex (Alchemilla sp., Hypericum revolutum, Erica arborea and Rumex sp.) were also among the most preferred diet items of the domestic goat. These results indicate that there is potential for competition between the two species, especially during the dry season, when resources are limited. Our findings, in combination with the expected increase in domestic herbivores, suggest that management plans should consider the potential threat posed by domestic goats to ensure future survival of the endangered Walia ibex.

Highly overlapping winter diet in two sympatric lemming species revealed by DNA metabarcoding.

Sympatric species are expected to minimize competition by partitioning resources, especially when these are limited. Herbivores inhabiting the High Arctic in winter are a prime example of a situation where food availability is anticipated to be low, and thus reduced diet overlap is expected. We present here the first assessment of diet overlap of high arctic lemmings during winter based on DNA metabarcoding of feces. In contrast to previous analyses based on microhistology, we found that the diets of both collared (Dicrostonyx groenlandicus) and brown lemmings (Lemmus trimucronatus) on Bylot Island were dominated by Salix while mosses, which were significantly consumed only by the brown lemming, were a relatively minor food item. The most abundant plant taxon, Cassiope tetragona, which alone composes more than 50% of the available plant biomass, was not detected in feces and can thus be considered to be non-food. Most plant taxa that were identified as food items were consumed in proportion to their availability and none were clearly selected for. The resulting high diet overlap, together with a lack of habitat segregation, indicates a high potential for resource competition between the two lemming species. However, Salix is abundant in the winter habitats of lemmings on Bylot Island and the non-Salix portion of the diets differed between the two species. Also, lemming grazing impact on vegetation during winter in the study area is negligible. Hence, it seems likely that the high potential for resource competition predicted between these two species did not translate into actual competition. This illustrates that even in environments with low primary productivity food resources do not necessarily generate strong competition among herbivores.

Sources of variation in small rodent trophic niche: new insights from DNA metabarcoding and stable isotope analysis.

Intraspecific competition for food is expected to increase the trophic niche width of consumers, defined here as their diet diversity, but this process has been little studied in herbivores. Population densities of small rodents fluctuate greatly, providing a good study model to evaluate effects of competition on trophic niche. We studied resource use in five arctic small rodent populations of four species combining DNA metabarcoding of stomach contents and stable isotope analysis (SIA). Our results suggest that for small rodents, the most pronounced effect of competition on trophic niche is due to increased use of secondary habitats and to habitat-specific diets, rather than an expansion of trophic niche in primary habitat. DNA metabarcoding and SIA provided complementary information about the composition and temporal variation of herbivore diets. Combing these two approaches requires caution, as the underlying processes causing observed patterns may differ between methodologies due to different spatiotemporal scales.

Genetics of cryptic speciation within an Arctic mustard, Draba nivalis.

Crossing experiments indicate that hybrid sterility barriers frequently have developed within diploid, circumpolar plant species of the genus Draba. To gain insight into the rapid evolution of postzygotic reproductive isolation in this system, we augmented the linkage map of one of these species, D. nivalis, and searched for quantitative trait loci (QTLs) associated with reproductive isolation. The map adds 63 new dominant markers to a previously published dataset of 31 co-dominant microsatellites. These markers include 52 amplified fragment length polymorphisms (AFLPs) and 11 sequence-specific amplified polymorphisms (SSAPs) based on retrotransposon sequence. 22 markers displaying transmission ratio distortion were further included in the map. We resolved eight linkage groups with a total map length of 894 cM. Significant genotype-trait associations, or quantitative trait loci (QTL), were detected for reproductive phenotypes including pollen fertility (4 QTLs), seed set (3 QTLs), flowering time (3 QTLs) and number of flowers (4 QTLs). Observed patterns of inheritance were consistent with the influence of both nuclear-nuclear interactions and chromosomal changes on these traits. All seed set QTLs and one pollen fertility QTL displayed underdominant effects suggestive of the involvement of chromosomal rearrangements in hybrid sterility. Interestingly, D. nivalis is predominantly self-fertilizing, which may facilitate the establishment of underdominant loci and contribute to reproductive isolation.