We estimated the heritability of gene expression in DGN attributable to SNPs in the vicinity of each gene using a mixed-effects model (see Materials and Methods) and calculated variances using restricted maximum likelihood as implemented in GCTA20. We use only local SNPs since we found that heritability estimates using all genotyped SNPs were too noisy to make meaningful inferences.
