Primer sequences used in MPRA cloning.

Interbreeding between anatomically modern humans and archaic hominins has contributed to the genomes of present-day human populations. However, our understanding of the specific gene regulatory consequences of Neanderthal, and particularly, Denisovan introgression is limited. Here, we used a massively parallel reporter assay to investigate the regulatory effects of 25,869 high-confidence introgressed SNPs segregating in present-day individuals of Papuan genetic ancestry in immune cell types. Overall, 8.22% of Denisovan and 8.58% of Neanderthal sequences showed active regulatory activity, and 9.22% of these displayed differential activity between archaic and modern alleles. We found no association between introgressed allele frequency on activity regardless of introgression source, but introgressed Denisovan alleles at higher frequencies were less likely to be differentially active than expected, suggesting introgression is under some degree of selective constraint. Both activity and differentially activity were associated with distance to the nearest transcription start site, while differential activity was additionally associated with differential transcription factor binding. Genes predicted to be regulated by differentially active sequences included IFIH1 and TNFAIP3, key immune genes and known examples of archaic introgression. Overall, this work provides experimental validation of regulatory activity for thousands of archaic variants in populations with the highest levels of Denisovan ancestry worldwide, revealing how human evolutionary history actively shapes present-day genetic diversity and immune function.