Multiple Pleistocene refugia for Arctic White Heather (Cassiope tetragona) supported by population genomics analyses of contemporary and Little-Ice-Age samples

AimArctic plants survived the Pleistocene glaciations in unglaciated refugia, but the number of these refugia is often unclear. We use high-resolution genomic data from present-day and Little-Ice-Age populations of Arctic White Heather (Cassiope tetragona) to re-evaluate the biogeography of this species and determine whether it had multiple independent refugia or a single refugium in Beringia. LocationCircumpolar Arctic and Coastal British Columbia (BC) alpine TaxonCassiope tetragona L., subspecies saximontana and tetragona, outgroup C. mertensiana (Ericaceae) MethodsWe built genotyping-by-sequencing (GBS) libraries using Cassiope tetragona tissue from 36 Arctic locations, including two [~]250-500-year-old populations collected under glacial ice on Ellesmere Island, Canada. We assembled a de novo GBS reference and called variants in dDocent. Population structure, genetic diversity, and demography were inferred from PCA, ADMIXTURE, fastsimcoal2, SplitsTree, and several population genomics statistics. ResultsPopulation structure analyses identified 4-5 clusters that align with geographic locations. Nucleotide diversity was highest in Beringia and decreased eastwards across Canada. Demographic coalescent analysis of the site-frequency-spectrum dated the following splits from Alaska: BC subspecies saximontana (6 mya), Russia (1.5 mya), Europe (>300-600 kya), Greenland (100 kya). Northern Canada populations appear to be from the current interglacial (7-9 kya). Genetic variants from Alaska appeared more frequently in present-day than historic plants on Ellesmere Island. ConclusionsDemographic analyses show BC, Alaska, Russia, Europe, and Greenland all had separate refugia during the last major glaciations. Northern Canadian populations appear to be founded during the current interglacial with genetic contributions from Alaska, Europe, and Greenland. On Ellesmere Island, there is evidence for continued, recent gene flow with foreign variants introduced in the last 250-500 years. These results suggest that a re-analysis of other Arctic species with shallow population structure using higher resolution genomic markers and demographic analyses may help reveal deeper structure and other circumpolar glacial refugia.