An assessment of gene-by-gene interactions as a tool to unfold missing heritability in dyslexia

Even if substantial heritability has been reported and candidate genes have been identified extensively, all known marker associations explain only a small proportion of the phenotypic variance of developmental dyslexia (DD) and related quantitative phenotypes. Gene-by-gene interaction (also known as “epistasis”—G × G) triggers a non-additive effect of genes at different loci and should be taken into account in explaining part of the missing heritability of this complex trait. We assessed potential G × G interactions among five DD candidate genes, i.e., DYX1C1, DCDC2, KIAA0319, ROBO1, and GRIN2B, upon DD-related neuropsychological phenotypes in 493 nuclear families with DD, by implementing two complementary regression-based approaches: (1) a general linear model equation whereby the trait is predicted by the main effect of the number of rare alleles of the two genes and by the effect of the interaction between them, and (2) a family-based association test to detect G × G interactions between two unlinked markers by splitting up the association effect into a between- and a within-family genetic orthogonal components. After applying 500,000 permutations and correcting for multiple testing, both methods show that G × G effects between markers within the DYX1C1, KIAA0319/TTRAP, and GRIN2B genes lower the memory letters composite z-score of on average 0.55 standard deviation. We provided initial evidence that the effects of familial transmission of synergistic interactions between genetic risk variants can be exploited in the study of the etiology of DD, explain part of its missing heritability, and assist in designing customized charts of individualized neurocognitive impairments in complex disorders, such as DD.