Discovery and characterisation of gene by environment and epistatic genetic effects in a vertebrate model

Phenotypic variation arises from the interplay between genetic and environmental factors. However, disentangling these interactions for complex traits remains challenging in observational cohorts such as human biobanks. Instead, model organisms where genetic (G) and environmental (E) variation can be controlled offer a valuable complement to human studies in the analysis of higher-order genetic effects such as GxE interactions, dominance, and epistasis. Here, we utilized 76 medaka strains of the Medaka Inbred Kiyosu-Karlsruhe (MIKK) panel, to compare heart rate plasticity across temperatures. An F2 segregation analysis identified 16 quantitative trait loci (QTLs), with many exhibiting dominance, GxE, GxG, and GxGxE interactions. We experimentally validated four candidate genes using gene editing, revealing their temperature-sensitive impact on heart function. Finally, we devised simulations to assess how GWAS discovery power is influenced by the choice of statistical models, showing that the apparent additivity in human GWAS is to be expected given study design and sample sizes of current studies. This work demonstrates the value of controlled model organism studies for dissecting the genetics of complex traits and provides guidance on the design of genetic association studies.