Genomic and phenotypic consequences of two independent secondary contact zones between allopatric lineages of the anadromous ice goby Leucopsarion petersii

Genetic and phenotypic analyses of independent secondary contact zones between certain pairs of divergent populations offer powerful opportunities to assess whether the consequences vary with different environmental backgrounds. Populations of the ice goby Leucopsarion petersii are distributed throughout the Japanese archipelago and comprise genetically and phenotypically divergent groups in the Japan Sea and the Pacific Ocean. In particular, populations in the Japan Sea have a larger body size and numbers of vertebrae than those in the Pacific Ocean. Herein, we performed integrated analyses of genotypes and phenotypes of two independent secondary contact zones and investigated their consequences. Population genetic analyses revealed asymmetric introgression of the mitochondrial genome of either lineage relative to little admixture of nuclear genomes in both secondary contact zones. On phenotype analyses, vertebral numbers were clearly explained by nuclear genomic ancestry in both secondary contact zones, whereas body size was not, suggesting that a little introgression of nuclear genes regulates body size. Actually, we observed biased introgression of a candidate gene, neuropeptide Y (NPY), which potentially controls body size in the ice goby. Moreover, the body size changes in the introgressed populations possibly affect the introgression patterns of mitochondrial genomes across these zones. Collectively, our results demonstrated that genomic and phenotypic consequences of secondary contact varied in marine variable environments.