Human follicular helper T cell promoter connectomes reveal novel genes and regulatory elements at SLE GWAS loci

Abstract

Systemic lupus erythematosus (SLE) is a complex inflammatory disease mediated by autoreactive antibodies that damages multiple tissues in children and adults. Genome-wide association studies (GWAS) have statistically implicated hundreds of loci in the susceptibility to human disease, including SLE, but the majority have failed to identify the causal variants or the effector genes. As a physicochemical approach to detecting functional variants and connecting them to target genes, we generated comprehensive, high-resolution maps of SLE variant accessibility and gene connectivity in the context of the three-dimensional chromosomal architecture of human tonsillar follicular helper T cells (TFH), a cell type required for the production of anti-nuclear antibodies characteristic of SLE. These spatial epigenomic maps identified a shortlist of over 400 potentially functional variants across 48 GWAS-implicated SLE loci. Twenty percent of these variants were located in open promoters of highly-expressed TFH genes, while 80% reside in non-promoter genomic regions that are connected in 3D to genes that likewise tend to be highly expressed in TFH. Importantly, we find that 90% of SLE-associated variants exhibit spatial proximity to genes that are not nearby in the 1D sequence of the genome, and over 60% of variants ‘skip’ the nearest gene to physically interact only with the promoters of distant genes. Gene ontology confirmed that genes in spatial proximity to SLE variants reside in highly SLE-relevant networks, including accessible variants that loop 200-1000 kb to interact with the promoters of the canonical TFH genes BCL6 and CXCR5. CRISPR-Cas9 genome editing confirmed that these variants reside in novel, distal regulatory elements required for normal BCL6 and CXCR5 expression by T cells. Furthermore, SLE-associated SNP-promoter interactomes implicated a set of novel genes with no known role in TFH or SLE disease biology, including the homeobox-interacting protein kinase HIPK1 and the Ste kinase homolog MINK1. Targeting these kinases in primary human TFH cells inhibited production of IL-21, a requisite cytokine for production of class-switched antibodies by B cells. This 3D-variant-to-gene mapping approach gives mechanistic insight into the SLE-associated regulatory architecture of the human genome.