A migration-associated supergene reveals loss of biocomplexity in Atlantic cod

Intraspecific phenotypic diversity is integral to ecological resilience and the provision of ecosystem services1. Chromosome structural variation may underpin intraspecific diversity and complex phenotypes2 by reducing recombination within supergenes containing linked, co-adapted alleles. Connecting ecologically-relevant phenotypes to genomic variation can enable more precise conservation of exploited marine species by protecting important genetic diversity3,4. Here, using genome-wide association analysis of a 12K single nucleotide polymorphism (SNP) array we confirm that an ancient, derived chromosomal rearrangement consisting of two adjacent inversions is strongly associated with migratory phenotype and individual-level genetic structure in Atlantic cod (Gadus morhua) across the Northwest Atlantic. The presence of all identified migration-associated loci within this rearrangement indicates that pervasive variation in migration phenotype is in part controlled by a recombination-resistant supergene, facilitating fine-scale individual phenotypic variation within Northern cod. Furthermore, we reconstruct trends in effective population size over the last century, and find genomic signatures of population collapse, and different patterns of population expansion and decline among individuals based on supergene alleles. We demonstrate declines in effective population size consistent with the onset of industrialized harvest (post 1950) and substantially reduced effective size of individuals homozygous for the derived chromosomal rearrangement relative to heterozygous individuals or those homozygous for the ancestral version of this chromosomal region. These results illustrate how chromosomal structural diversity can mediate fine-scale genetic and phenotypic variation in a highly connected marine species, and suggest a loss of biocomplexity from a migration-associated supergene within Northern cod by overfishing.