Multi-species analysis of inflammatory response elements reveals ancient and lineage-specific contributions of transposable elements to NF-κB binding

Transposable elements (TEs) provide a source of transcription factor binding sites that can rewire conserved gene regulatory networks. NF-{kappa}B is an evolutionary conserved transcription factor complex primarily involved in innate immunity and inflammation. The extent to which TEs have contributed to NF-{kappa}B responses during mammalian evolution is not well established. Here we performed a multi-species analysis of TEs bound by the NF-{kappa}B subunit RELA (also known as p65) in response to the proinflammatory cytokine TNF. By comparing RELA ChIP-seq data from TNF-stimulated primary aortic endothelial cells isolated from human, mouse and cow, we found that 55 TE subfamilies were associated with RELA bound regions. These RELA-bound transposons possess active epigenetic features and reside near TNF-responsive genes. A prominent example of lineage-specific contribution of transposons comes from the bovine SINE subfamilies Bov-tA1/2/3 which collectively contributed over 14,000 RELA bound regions in cow. By comparing RELA binding data across species, we also found several examples of RELA motif-bearing TEs that colonized the genome prior to the divergence of the three species and contributed to species-specific RELA binding. For example, we found human RELA bound MER81 instances were enriched for the interferon gamma pathway and demonstrated that one RELA bound MER81 element can control the TNF-induced expression of Interferon Gamma Receptor 2 (IFNGR2). Using ancestral reconstructions, we found that RELA containing MER81 instances rapidly decayed during early primate evolution (> 50 million years ago (MYA)) before stabilizing since the separation of Old World monkeys (< 50 MYA). Taken together, our results suggest ancient and lineage-specific transposon subfamilies contributed to mammalian NF-{kappa}B regulatory networks.