Total RNA was extracted from EBV transformed lymphoblastoid cell line pellets by the TRIzol reagent (Ambion), and mRNA and small RNA sequencing of 465 unique individuals was performed on the Illumina HiSeq2000 platform, with paired-end 75bp mRNA-seq and single-end 36bp small RNA-seq. Five samples were sequenced in replicate in each of the seven sequencing laboratories. The mRNA and small RNA reads were mapped with GEM31 and miraligner32 (link), respectively, with an average of 48.9M mRNA-seq reads and 1.2M miRNA reads per sample after QC. Numerous transcript features were quantified using Gencode v1233 (link) and miRBase v1834 (link) annotations: protein-coding and lincRNA genes (16,084 detected in >50% of samples), transcripts (67,603; with FluxCapacitor7 (link)), exons (146,498), annotated splice junctions (129,805; analyzed in detail in Ferreira et al. submitted), transcribed repetitive elements (47,409), and mature miRNAs (715). Data quality was assessed by sample correlations and read and gene count distributions, and technical variation was removed by PEER normalization35 (link) for the QTL and miRNA-mRNA correlation analyses11 (link). The samples clustered uniformly both before and after normalization. The genotype data was obtained from 1000 Genomes Phase 1 data set for 421 samples (80× average exome and 5× whole genome read depth), and the remaining 41 samples were imputed from Omni 2.5M SNP array data. Furthermore, we did functional reannotation for all the 1000 Genomes variants using Gencode v12. QTL mapping was done with linear regression, using genetic variants with >5% frequency in 1MB window and normalized quantifications transformed to standard normal. Permutations were used to adjust FDR to 5%. Full details are provided in Supplementary Methods.