8Thousands of unfixed transposable element (TE) insertions segregate in the human 9 population, but little is known about their impact on genome function. Recently, a few 10 studies associated polymorphic TE insertions to mRNA levels of adjacent genes, but the 11 biological significance of these associations, their replicability across cell types, and the 12 mechanisms by which they may regulate genes remain largely unknown. Here we 13 performed a TE-expression QTL analysis of 444 lymphoblastoid cell lines and 294 14 induced pluripotent stem cells using a newly developed set of genotypes for 2,806 15 polymorphic TE insertions. We identified 211 and 324 TE-eQTL acting in cis in each 16 respective cell type. Approximately one fourth were shared across cell types with 17 strongly correlated effects. Furthermore, analysis of chromatin accessibility QTL in a 18 subset of the lymphoblastoid cell lines suggests that unfixed TEs often modulate the 19 activity of enhancers and other distal regulatory DNA elements, which tend to lose 20 accessibility when a TE inserts within them. We also document a case of an unfixed TE 21 likely influencing gene expression at the post-transcriptional level. Our study points to 22 broad and diverse cis-regulatory effects of unfixed TEs in the human population and 23 underscores their plausible contribution to phenotypic variation. 24 25 Keywords: unfixed transposable elements, gene expression, chromatin accessibility, 26 QTL, human, lymphoblastoid cell line, induced pluripotent stem cell 27 48 The regulation of gene expression is central to cellular function and differentiation in 49 development and physiology and changes in gene expression are important drivers of 50 phenotypic variation[20]. Steady state mRNA levels partially govern gene expression in 51 response to inputs integrated by various regulatory sequences acting in cis or 52trans [20,21]. Expression quantitative trait loci (eQTL) studies offer a systematic 53 approach to identify such cis-regulatory elements by correlating the genotypes of 54 genomic variants segregating in individuals with mRNA levels of specific genes, which 55 can be measured on a large scale by RNA-sequencing (RNA-seq) [22,23]. Previous 56 eQTL studies have established that virtually every gene in the human genome has its 57 expression affected by at least one genomic variant (typically a single nucleotide 58 polymorphism, SNP) located in cis (usually defined as within a maximum distance of 1 59Mb) and generally in non-coding regions [21,[24][25][26]. A recent analysis of 44 different 60 tissues by the GTEx consortium indicates that such cis-eQTL fall into two broad 61 categories: those shared across most tissues and those apparently acting in a single or 62 a restricted number of similar tissues [21]. Importantly, most eQTL studies thus far have 63 focused on SNPs, yet other type of genomic variants, such as unfixed TEs insertions, 64 are common in the human population and likely to have more drastic effects on gene 65 expression [27][28][29][30]. 66 67 I...