Abstract
Ornamental medaka strains derived from wild Japanese medaka (Oryzias latipes species complex) are bred worldwide. Over 200 years of selective breeding have produced over 700 strains with a wide variety of phenotypes, including diverse body coloration, scales, eyeball morphology, and fin and body shapes. In this study, we first identified and described 34 phenotypes in ornamental medaka strains. To understand the genomic basis of this phenotypic diversity and the domestication process, we performed whole-genome sequencing on 181 individuals of 86 ornamental Japanese medaka strains. Population genomic analyses revealed that modern ornamental medaka strains are genetically closer to the wild Southern Japan population of the Kansai-Setouchi regions, suggesting the origin of ornamental strains. In addition, the gene loci poc1a, tyr, nme2a, and gabrr2b have undergone selection during domestication. We performed genome-wide association studies analysis for 29 phenotypes observed in ornamental medaka strains and identified strong candidate genes for some phenotypes, including kcnq5a for hirenaga and swallow, bmp5 for deme, adcy5 for orochi, and kitlga for aurora, respectively. We found that loss of exon 8 of adcy5 caused melanism, a dark body color phenotype, in medaka, providing a molecular insight into this phenomenon in vertebrates and human familial dyskinesia. In addition, we uncovered the predominant candidate peaks of genome-wide association studies, including a total of 3,328 genes associated with 26 phenotypes. Our findings highlight the potential of population genomics to explore genotype-phenotype correlations and the genomic basis of body coloration and morphogenesis in medaka.