Abstract
Age‐at‐maturity and iteroparity are two life history variations of steelhead trout (Oncorhynchus mykiss) that are believed to increase population resilience and stability. While repeat‐spawning individuals are thought to have historically made up a substantial portion of the reproductive population in the Columbia River and the majority of females still attempt outmigration as kelts, return rates of repeat‐spawner are low throughout the basin and below 1% for the furthest migrating stocks. Notably, outmigrating adults exhibit variation in rematuration phenology, displaying either “consecutive” (reproduce immediately the following season) or “skip” (delay spawning for future seasons) spawning patterns. Here, we use low coverage whole genome sequencing of consecutive versus skip spawning female Columbia River steelhead from two populations to test for genomic differences between these two iteroparous phenotypes. We identified genomic regions on several chromosomes which were associated with the phenology of iteroparity, including a region on chromosome 25 containing two genes, estradiol receptor beta (ERβ) and glycoprotein hormone beta‐5 (GPHB5), which, in mammals, are estrogen‐sensitive and expressed in reproductive tissues. Allele frequencies in this ERβ/GPHB5 region differed among female steelhead of different age at maturity, but not males. These genes also shared an island of linkage disequilibrium with the SIX6 gene, 600Kbp away on the same chromosome, a region of known association with age‐at‐maturity. These observations contribute to growing evidence that age‐at‐maturity and the phenology of iteroparity are determined by overlapping physiological processes and genetic pathways.