Whole genome sequencing reveals how plasticity and genetic differentiation underlie sympatric morphs of Arctic charr.

Now published in Molecular Ecology

Whole genome sequencing reveals how plasticity and genetic differentiation underlie sympatric morphs of Arctic charr.

https://onlinelibrary.wiley.com/doi/10.1111/mec.70085

ABSTRACT

Salmonids have a remarkable ability to form sympatric morphs after postglacial colonisation of freshwater lakes. These morphs often differ in morphology, feeding and spawning behaviour. Here, we explored the genetic basis of morph differentiation in Arctic charr (n = 283) by first establishing a high-quality reference genome and then using this in whole genome sequencing of distinct morphs present in two Norwegian and two Icelandic lakes. The four lakes represent the spectrum of genetic differentiation between morphs from one lake with no genetic differentiation between morphs, implying phenotypic plasticity, to two lakes with locus-specific genetic differentiation, implying incomplete reproductive isolation, and one lake with strong genome-wide divergence consistent with complete reproductive isolation. As many as 12 putative inversions ranging from 0.45 to 3.25 Mbp in size segregated among the four morphs present in one lake, Thingvallavatn, and these contributed significantly to the genetic differentiation among morphs. None of the putative inversions were found in any of the other lakes, but there were cases of partial haplotype sharing in similar morph contrasts in other lakes. Our findings are consistent with a highly polygenic basis of morph differentiation with population-specific selection on alleles linked to the development of similar morph phenotypes. The results support a model where morph differentiation is first established through phenotypic plasticity, leading to niche expansion and separation. This may be followed by gradual development of reproductive isolation, locus-specific differentiation and eventually complete reproductive isolation and genome-wide divergence.

FIGURE 4
Genetic differentiation between small and large benthivorous morphs of Arctic charr in Thingvallavatn. (a) Genome scan based on estimated allele frequencies for individual SNPs. Gene annotations are shown for the top 0.001% of significant SNPs at putative inversions marked by stars. (b–e) Zoom-in profile of four putative inversion regions in (a). Top panel: Linkage disequilibrium (LD), measured as r2, between the top significantly associated SNP and surrounding variants; bottom panel: Pairwise LD (r2) among all genotypes within each region, including the following genomic regions (b) Scaffold 4: 75.25–76.13 Mb, (c) Scaffold 5: 22.30–22.75 Mb, (d) Scaffold 9: 61.30–62.11 Mb and (e) Scaffold 17: 32.45–33.20 Mb. The horizontal red and blue lines indicate Bonferroni corrected significance thresholds with significance levels α = 10−3 and α = 10−8, respectively, after taking into account that 107 SNPs were used. (f) Wilcoxon rank test of theta (θ, nucleotide diversity) distributions across the genome and within putative inversion regions. Significance levels: *p < 0.05, **p < 0.01, ***p < 0.001, ns indicates p > 0.05.

About Arnar Pálsson

Arnar Palsson received his bachelor and Masters degrees from University of Iceland and Ph.D. from the Department of Genetics at North Carolina State University. He worked as post-doctoral fellow at the Department of Ecology and Evolution, University of Chicago.
This entry was posted in Publications. Bookmark the permalink.