Comparative genomic analysis and genetic variation within the leafy liverwort genus Scapania

The genus Scapania comprises a group of leafy liverworts distributed throughout many bryophytic assemblages. While many Scapania species grow widely, some are assessed as endangered and appear to be specialists with distinct niche environments. Several are found only in alpine forest communities, inhabiting decaying logs in streams, typical of an environment that is threatened by both logging activity and changes to watercourses. Another species, S. nimbosa, has an unusual Oceanic-Montane distribution across Ireland, Scotland, Norway, China and Nepal. Since gemmae and sexual reproduction are absent the species is hypothesized to be primarily dispersed by fragmentation. In Norway S. nimbosa occupies an area of only 13 x 20 km, at altitudes between 300-980 m, and is frequently found with another more abundant asexual species, S. ornithopodioides. This makes S. nimbosa susceptible to local extinction through climate change or perhaps interspecific competition. Genomics is being increasingly used to infer demography and the evolutionary history of a species. Ascertaining levels of genetic variation can also contribute towards an effective conservation management plan. Besides, very little is known about the genomic organization and sexual determination in leafy liverworts. To generate new knowledge about the genus Scapania we sequenced the genomes of the sexual species S. nemorea (both male and female isolates), S. undulata (a single isolate), and several asexual S. ornithopodiodes and S. nimbosa isolates. Illumina paired-end (2x 300 bp) and Oxford Nanopore long reads were used to create genomic references. Initially organellar genomes were assembled, annotated and genetic variation was discovered. This revealed that variation is indeed present even for S. nimbosa and S. ornithopodioides at Norwegian sites. Next we focussed on creating a high quality nuclear reference genome for S. nemorea using the SPAdes assembler (v3.13). Qualities of each assembly and isolate were assessed with QUAST and BUSCO. While one assembly spans 202.6 Mb (10930 scaffolds; N50 of 66 Kb), other isolates of S. nemorea show larger assembled genome sizes and different Kmer distributions, consistent with the expected alternative sexual chromosome complement. We further analyse genomic synteny and diversity, but emphasize that difficulties in extracting DNA from herbarium specimens really hamper analysis.