Large scale functional screen identifies genetic variants with splicing effects in modern and archaic humans

Humans co-existed and interbred with other hominins which later became extinct. These archaic hominins are known to us only through fossil records and for two cases, genome sequences. Here we engineer Neanderthal and Denisovan sequences into thousands of artificial genes to reconstruct the pre-mRNA processing patterns of these extinct populations. Of the 5,224 alleles tested in this massively parallel splicing reporter assay (MaPSy), we report 969 exonic splicing mutations (ESMs) that correspond to differences in exon recognition between extant and extinct hominins. Using MaPSy splicing variants, predicted splicing variants, and splicing quantitative trait loci, we show that splice-disrupting variants experienced greater purifying selection in anatomically modern humans than in Neanderthals. Adaptively introgressed variants were enriched for moderate effect splicing variants, consistent with positive selection for alternative spliced alleles following introgression. As particularly compelling examples, we characterized a novel tissue-specific alternative splicing variant at the adaptively introgressed innate immunity gene TLR1, as well as a novel Neanderthal introgressed alternative splicing variant in the gene HSPG2 that encodes perlecan. We further identified potentially pathogenic splicing variants found only in Neanderthals and Denisovans in genes related to sperm maturation and immunity. Finally, we found splicing variants that may contribute to variation among modern humans in total bilirubin, balding, hemoglobin levels, and lung capacity. Our findings provide novel insights into natural selection acting on splicing in human evolution and demonstrate how functional assays can be used to identify candidate causal variants underlying differences in gene regulation and phenotype. Overall design: Illumina 2x150 paired-end sequencing data for massively parallel splicing assay (MaPSy) experiment with two separately transfected minigene pools (labeled N2A, N2B) each with three input replicates and three output library replicates (labeled a,b,c).