Intensity data from the Illumina genotyping platform (namely log R ratio (LRR) and B allele frequency (BAF)) were first normalized by correcting for GC content and Loess normalization, and then copy number variant (CNV)-calling was performed using the cnvHap software55. We analyzed all large CNVs (>100 kb) with a frequency >5% for any significant associations with gene expression.