Genetic diversity is fundamental to the adaptive potential and survival of species. Although its importance has long been recognized in science, it has a history of neglect within policy, until now. The new Global Biodiversity Framework recently adopted by the Convention on Biological Diversity, states that genetic diversity must be maintained at levels assuring adaptive potential of populations, and includes metrics for systematic monitoring of genetic diversity in so called indicators. Similarly, indicators for genetic diversity are being developed at national levels. Here, we apply new indicators for Swedish national use to one of the northernmost salmonid fishes, the Arctic charr (Salvelinus alpinus). We sequence whole genomes to monitor genetic diversity over four decades in three landlocked populations inhabiting protected alpine lakes in central Sweden. We find levels of genetic diversity, inbreeding and load to differ among lakes but remain stable over time. Effective population sizes are generally small (< 500), suggesting a limited ability to maintain adaptive variability if genetic exchange with nearby populations became eliminated. We identify genomic regions potentially shaped by selection; SNPs exhibiting population divergence exceeding expectations under drift and a putative selective sweep acting within one lake to which the competitive brown trout (Salmo trutta) was introduced during the sampling period. Identified genes appear involved in immunity and salinity tolerance. Present results suggest that genetically vulnerable populations of Arctic charr have maintained neutral and putatively adaptive genetic diversity despite small effective sizes, attesting the importance of continued protection and assurance of gene flow among populations.

Satellite DNA (satDNA) is a fast-evolving portion of eukaryotic genomes. The homogeneous and repetitive nature of such satDNA causes problems during the assembly of genomes, and therefore it is still difficult to study it in detail in nonmodel organisms as well as across broad evolutionary timescales. Here, we combined the use of short- and long-read data to explore the diversity and evolution of satDNA between individuals of the same species and between genera of birds spanning ~40 millions of years of bird evolution using birds-of-paradise (Paradisaeidae) and crow (Corvus) species. These avian species highlighted the presence of a GC-rich Corvoidea satellitome composed of 61 satellite families and provided a set of candidate satDNA monomers for being centromeric on the basis of length, abundance, homogeneity and transcription. Surprisingly, we found that the satDNA of crow species rapidly diverged between closely related species while the satDNA appeared more similar between birds-of-paradise species belonging to different genera.

Background: Many wild species have suffered drastic population size declines over the past centuries, which have led to 'genomic erosion' processes characterized by reduced genetic diversity, increased inbreeding, and accumulation of harmful mutations. Yet, genomic erosion estimates of modern-day populations often lack concordance with dwindling population sizes and conservation status of threatened species. One way to directly quantify the genomic consequences of population declines is to compare genome-wide data from pre-decline museum samples and modern samples. However, doing so requires computational data processing and analysis tools specifically adapted to comparative analyses of degraded, ancient or historical, DNA data with modern DNA data as well as personnel trained to perform such analyses.

Results: Here, we present a highly flexible, scalable, and modular pipeline to compare patterns of genomic erosion using samples from disparate time periods. The GenErode pipeline uses state-of-the-art bioinformatics tools to simultaneously process whole-genome re-sequencing data from ancient/historical and modern samples, and to produce comparable estimates of several genomic erosion indices. No programming knowledge is required to run the pipeline and all bioinformatic steps are well-documented, making the pipeline accessible to users with different backgrounds. GenErode is written in Snakemake and Python3 and uses Conda and Singularity containers to achieve reproducibility on high-performance compute clusters. The source code is freely available on GitHub (https://github.com/NBISweden/GenErode).

Conclusions: GenErode is a user-friendly and reproducible pipeline that enables the standardization of genomic erosion indices from temporally sampled whole genome re-sequencing data.

Supergenes govern multi-trait-balanced polymorphisms in a wide range of systems; however, our understanding of their origins and evolution remains incomplete. The reciprocal placement of stigmas and anthers in pin and thrum floral morphs of distylous species constitutes an iconic example of a balanced polymorphism governed by a supergene, the distyly S-locus. Recent studies have shown that the Primula and Turnera distyly supergenes are both hemizygous in thrums, but it remains unknown whether hemizygosity is pervasive among distyly S-loci. As hemizygosity has major consequences for supergene evolution and loss, clarifying whether this genetic architecture is shared among distylous species is critical. Here, we have characterized the genetic architecture and evolution of the distyly supergene in Linum by generating a chromosome-level genome assembly of Linum tenue, followed by the identification of the S-locus using population genomic data. We show that hemizygosity and thrum-specific expression of S-linked genes, including a pistil-expressed candidate gene for style length, are major features of the Linum S-locus. Structural variation is likely instrumental for recombination suppression, and although the non-recombining dominant haplotype has accumulated transposable elements, S-linked genes are not under relaxed purifying selection. Our findings reveal remarkable convergence in the genetic architecture and evolution of independently derived distyly supergenes, provide a counterexample to classic inversion-based supergenes, and shed new light on the origin and maintenance of an iconic floral polymorphism.

The major histocompatibility complex (MHC) is of central importance to the immune system, and an optimal MHC diversity is believed to maximize pathogen elimination. Birds show substantial variation in MHC diversity, ranging from few genes in most bird orders to very many genes in passerines. Our understanding of the evolutionary trajectories of the MHC in passerines is hampered by lack of data on genomic organization. Therefore, we assembled and annotated the MHC genomic region of the great reed warbler (Acrocephalus arundinaceus), using long-read sequencing and optical mapping. The MHC region is large (>5.5 Mb), characterized by structural changes compared to hitherto investigated bird orders and shows higher repeat content than the genome average. These features were supported by analyses in three additional passerines. MHC genes in passerines are found in two different chromosomal arrangements, either as single copy MHC genes located among non-MHC genes, or as tandemly duplicated tightly linked MHC genes. Some single copy MHC genes are old and putative orthologues among species. In contrast tandemly duplicated MHC genes are monophyletic within species and have evolved by simultaneous gene duplication of several MHC genes. Structural differences in the MHC genomic region among bird orders seem substantial compared to mammals and have possibly been fuelled by clade-specific immune system adaptations. Our study provides methodological guidance in characterizing complex genomic regions, constitutes a resource for MHC research in birds, and calls for a revision of the general belief that avian MHC has a conserved gene order and small size compared to mammals.

The sympatric existence of genetically distinguishable populations of the same species remains a puzzle in ecology. Coexisting salmonid fish populations are known from over 100 freshwater lakes. Most studies of sympatric populations have used limited numbers of genetic markers making it unclear if genetic divergence involves certain parts of the genome. We returned to the first reported case of salmonid sympatry, initially detected through contrasting homozygosity at a single allozyme locus (coding for lactate dehydrogenase A) in brown trout in the small Lakes Bunnersjöarna, Sweden. First, we verified the existence of the two coexisting demes using a 96-SNP fluidigm array. We then applied whole-genome resequencing of pooled DNA to explore genome-wide diversity within and between these demes; nucleotide diversity was higher in deme I than in deme II. Strong genetic divergence is observed with genome-wide F_ST ≈ 0.2. Compared with data from populations of similar small lakes, this divergence is of similar magnitude as that between reproductively isolated populations. Individual whole-genome resequencing of two individuals per deme suggests higher inbreeding in deme II versus deme I, indicating different degree of isolation. We located two gene-copies for LDH-A and found divergence between demes in a regulatory section of one of these genes. However, we did not find a perfect fit between the sequence data and previous allozyme results, and this will require further research. Our data demonstrates genome-wide divergence governed mostly by genetic drift but also by diversifying selection in coexisting populations. This type of hidden biodiversity needs consideration in conservation management.

How the avian sex chromosomes first evolved from autosomes remains elusive as 100 million years (My) of divergence and degeneration obscure their evolutionary history. The Sylvioidea group of songbirds is interesting for understanding avian sex chromosome evolution because a chromosome fusion event 24 Ma formed “neo-sex chromosomes” consisting of an added (new) and an ancestral (old) part. Here, we report the complete female genome (ZW) of one Sylvioidea species, the great reed warbler (Acrocephalus arundinaceus). Our long-read assembly shows that the added region has been translocated to both Z and W, and whereas the added-Z has retained its gene order the added-W part has been heavily rearranged. Phylogenetic analyses show that recombination between the homologous added-Z and -W regions continued after the fusion event, and that recombination suppression across this region took several million years to be completed. Moreover, recombination suppression was initiated across multiple positions over the added-Z, which is not consistent with a simple linear progression starting from the fusion point. As expected following recombination suppression, the added-W show signs of degeneration including repeat accumulation and gene loss. Finally, we present evidence for nonrandom maintenance of slowly evolving and dosage-sensitive genes on both ancestral- and added-W, a process causing correlated evolution among orthologous genes across broad taxonomic groups, regardless of sex linkage.

Highly endangered species like the Sumatran rhinoceros are at risk from inbreeding. Five historical and 16 modern genomes from across the species range show mutational load, but little evidence for local adaptation, suggesting that future inbreeding depression could be mitigated by assisted gene flow among populations. Small populations are often exposed to high inbreeding and mutational load that can increase the risk of extinction. The Sumatran rhinoceros was widespread in Southeast Asia, but is now restricted to small and isolated populations on Sumatra and Borneo, and most likely extinct on the Malay Peninsula. Here, we analyse 5 historical and 16 modern genomes from these populations to investigate the genomic consequences of the recent decline, such as increased inbreeding and mutational load. We find that the Malay Peninsula population experienced increased inbreeding shortly before extirpation, which possibly was accompanied by purging. The populations on Sumatra and Borneo instead show low inbreeding, but high mutational load. The currently small population sizes may thus in the near future lead to inbreeding depression. Moreover, we find little evidence for differences in local adaptation among populations, suggesting that future inbreeding depression could potentially be mitigated by assisted gene flow among populations.

The advent of novel sequencing techniques has unraveled a tremendous diversity on Earth. Genomic data allow us to understand ecology and function of organisms that we would not otherwise know existed. However, major methodological challenges remain, in particular for multicellular organisms with large genomes. Arbuscular mycorrhizal (AM) fungi are important plant symbionts with cryptic and complex multicellular life cycles, thus representing a suitable model system for method development. Here, we report a novel method for large scale, unbiased nuclear sorting, sequencing, and de novo assembling of AM fungal genomes. After comparative analyses of three assembly workflows we discuss how sequence data from single nuclei can best be used for different downstream analyses such as phylogenomics and comparative genomics of single nuclei. Based on analysis of completeness, we conclude that comprehensive de novo genome assemblies can be produced from six to seven nuclei. The method is highly applicable for a broad range of taxa, and will greatly improve our ability to study multicellular eukaryotes with complex life cycles.

The ratio of nonsynonymous over synonymous sequence divergence, dN/dS, is a widely used estimate of the nonsynonymous over synonymous fixation rate ratio ω, which measures the extent to which natural selection modulates protein sequence evolution. Its computation is based on a phylogenetic approach and computes sequence divergence of protein-coding DNA between species, traditionally using a single representative DNA sequence per species. This approach ignores the presence of polymorphisms and relies on the indirect assumption that new mutations fix instantaneously, an assumption which is generally violated and reasonable only for distantly related species. The violation of the underlying assumption leads to a time-dependence of sequence divergence, and biased estimates of ω in particular for closely related species, where the contribution of ancestral and lineage-specific polymorphisms to sequence divergence is substantial. We here use a time-dependent Poisson random field model to derive an analytical expression of dN/dS as a function of divergence time and sample size. We then extend our framework to the estimation of the proportion of adaptive protein evolution α. This mathematical treatment enables us to show that the joint usage of polymorphism and divergence data can assist the inference of selection for closely related species. Moreover, our analytical results provide the basis for a protocol for the estimation of ω and α for closely related species. We illustrate the performance of this protocol by studying a population data set of four corvid species, which involves the estimation of ω and α at different time-scales and for several choices of sample sizes.