Genomic datasets continue to increase in number due to the ease of production for a wider selection of species including non-model organisms. For many of these species, especially those with large or polyploid genomes, highly contiguous and well-annotated genomes are still rare due to the complexity and cost involved in their assembly. As a result, a common starting point for genomic work in non-model species is the production of a linkage map. Dense linkage maps facilitate the analysis of genomic data in a variety of ways, from broad scale observations regarding genome structure e.g., chromosome number and type or sex-related structural differences, to fine scale patterns e.g., recombination rate variation and co-localization of differentiated regions. Here we present both sex-averaged and sex-specific linkage maps for Coregonus sp. Albock, a member of the European whitefish lineage (C. lavaretus spp. complex), containing 5395 single nucleotide polymorphism (SNP) loci across 40 linkage groups to facilitate future investigation into the genomic basis of whitefish adaptation and speciation. The map was produced using restriction-site associated digestion (RAD) sequencing data from two wild-caught parents and 156 F1 offspring. We discuss the differences between our sex-averaged and sex-specific maps and identify genome-wide synteny between C. sp. Albock and Atlantic Salmon (Salmo salar), which have diverged following the salmonid-specific whole genome duplication. Our analysis confirms that many patterns of synteny observed between Atlantic Salmon and Oncorhynchus and Salvelinus species are also shared by members of the Coregoninae subfamily. We also show that regions known for their species-specific rediploidization history can pose challenges for synteny identification since these regions have diverged independently in each salmonid species following the salmonid-specific whole genome duplication. The European whitefish map provided here will enable future studies to understand the distribution of loci of interest, e.g., FST outliers, along the whitefish genome as well as assisting with the de novo assembly of a whitefish reference genome.

Salmonids are of particular interest to evolutionary biologists due to their incredible diversity of life-history strategies and the speed at which many salmonid species have diversified. In Switzerland alone, over 30 species of Alpine whitefish from the subfamily Coregoninae have evolved since the last glacial maximum, with species exhibiting a diverse range of morphological and behavioural phenotypes. This, combined with the whole genome duplication which occurred in the ancestor of all salmonids, makes the Alpine whitefish radiation a particularly interesting system in which to study the genetic basis of adaptation and speciation and the impacts of ploidy changes and subsequent rediploidization on genome evolution. Although well-curated genome assemblies exist for many species within Salmonidae, genomic resources for the subfamily Coregoninae are lacking. To assemble a whitefish reference genome, we carried out PacBio sequencing from one wild-caught Coregonus sp. Balchen from Lake Thun to ~90× coverage. PacBio reads were assembled independently using three different assemblers, falcon, canu and wtdbg2 and subsequently scaffolded with additional Hi-C data. All three assemblies were highly contiguous, had strong synteny to a previously published Coregonus linkage map, and when mapping additional short-read data to each of the assemblies, coverage was fairly even across most chromosome-scale scaffolds. Here, we present the first de novo genome assembly for the Salmonid subfamily Coregoninae. The final 2.2-Gb wtdbg2 assembly included 40 scaffolds, an N50 of 51.9 Mb and was 93.3% complete for BUSCOs. The assembly consisted of ~52% transposable elements and contained 44,525 genes.

Microbes can have profound effects on their hosts, driving natural selection, promoting speciation and determining species distributions. However, soil-dwelling microbes are rarely investigated as drivers of evolutionary change in plants. We used metabarcoding and experimental manipulation of soil microbiomes to investigate the impact of soil and root microbes in a well-known case of sympatric speciation, the Howea palms of Lord Howe Island (Australia). Whereas H. forsteriana can grow on both calcareous and volcanic soils, H. belmoreana is restricted to, but more successful on, volcanic soil, indicating a trade-off in adaptation to the two soil types. We suggest a novel explanation for this trade-off. Arbuscular mycorrhizal fungi (AMF) are significantly depleted in H. forsteriana on volcanic soil, relative to both H. belmoreana on volcanic soil and H. forsteriana on calcareous soil. This is mirrored by the results of survival experiments, where the sterilization of natural soil reduces Howea fitness in every soil-species combination except H. forsteriana on volcanic soil. Furthermore, AMF-associated genes exhibit evidence of divergent selection between Howea species. These results show a mechanism by which divergent adaptation can have knock-on effects on host-microbe interactions, thereby reducing interspecific competition and promoting the coexistence of plant sister species.

Milkweed butterflies in the genus Danaus are studied in a diverse range of research fields including the neurobiology of migration, biochemistry of plant detoxification, host-parasite interactions, evolution of sex chromosomes, and speciation. We have assembled a nearly chromosomal genome for Danaus chrysippus (known as the African Monarch, African Queen, and Plain Tiger) using long-read sequencing data. This species is of particular interest for the study of genome structural change and its consequences for evolution. Comparison with the genome of the North American Monarch Danaus plexippus reveals generally strong synteny but highlights 3 inversion differences. The 3 chromosomes involved were previously found to carry peaks of intraspecific differentiation in D. chrysippus in Africa, suggesting that these inversions may be polymorphic and associated with local adaptation. The D. chrysippus genome is over 40% larger than that of D. plexippus, and nearly all of the additional ∼100 Megabases of DNA comprises repeats. Future comparative genomic studies within this genus will shed light on the evolution of genome architecture.

Adaptive radiations represent some of the most remarkable explosions of diversification across the tree of life. However, the constraints to rapid diversification and how they are sometimes overcome, particularly the relative roles of genetic architecture and hybridization, remain unclear. Here, we address these questions in the Alpine whitefish radiation, using a whole-genome dataset that includes multiple individuals of each of the 22 species belonging to six ecologically distinct ecomorph classes across several lake-systems. We reveal that repeated ecological and morphological diversification along a common environmental axis is associated with both genome-wide allele frequency shifts and a specific, larger effect, locus, associated with the gene edar. Additionally, we highlight the possible role of introgression between species from different lake-systems in facilitating the evolution and persistence of species with unique trait combinations and ecology. These results highlight the importance of both genome architecture and secondary contact with hybridization in fuelling adaptive radiation.

Supergenes maintain adaptive clusters of alleles in the face of genetic mixing. Although usually attributed to inversions, supergenes can be complex, and reconstructing the precise processes that led to recombination suppression and their timing is challenging. We investigated the origin of the BC supergene, which controls variation in warning coloration in the African monarch butterfly, Danaus chrysippus. By generating chromosome-scale assemblies for all three alleles, we identified multiple structural differences. Most strikingly, we find that a region of more than 1 million bp underwent several segmental duplications at least 7.5 Ma. The resulting duplicated fragments appear to have triggered four inversions in surrounding parts of the chromosome, resulting in stepwise growth of the region of suppressed recombination. Phylogenies for the inversions are incongruent with the species tree and suggest that structural polymorphisms have persisted for at least 4.1 Myr. In addition to the role of duplications in triggering inversions, our results suggest a previously undescribed mechanism of recombination suppression through independent losses of divergent duplicated tracts. Overall, our findings add support for a stepwise model of supergene evolution involving a variety of structural changes. This article is part of the theme issue Genomic architecture of supergenes: causes and evolutionary consequences.

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.