Both Systemic Lupus Erythematosus (SLE) and Systemic Sclerosis (SSc) are autoimmune diseases sharing similar genetic backgrounds. Genome-wide association studies (GWASs) have constantly disclosed numerous genetic variants conferring to both disease risks at 7q32.1, but the functional mechanisms underlying them are still largely unknown. Through combining fine-mapping and functional epigenomic analyses, we prioritized a potential independent functional SNP (rs13239597) within TNPO3 promoter region, residing in a putative enhancer element. Functional analysis integrating expression quantitative trait locus (eQTL) and high-throughput chromatin interaction (Hi-C) demonstrated that IRF5 is the distal target gene (~118kb) of rs13239597, which is a key regulator of pathogenic autoantibody dysregulation increased risk of both SLE and SSc. We experimentally validated the long-range chromatin interactions between rs13239597 and IRF5 using chromosome conformation capture (3C) assay. We further demonstrated that rs13239597-A acted as an allele-specific enhancer regulating IRF5 expression, independently of TNPO3 by using dual-luciferase reporter assays and CRISPR-Cas9. Particularly, the transcription factor EVI1 could preferentially bind to rs13239597-A allele and increase the enhancer activity to regulate IRF5 expression. Taken together, our results uncovered the mechanistic insight connecting between a noncoding functional variant with a distal immunologically pathogenic gene IRF5, which might obligate in understanding the complex genetic architectures of SLE and SSc pathogenesis.

Circular RNAs (circRNAs) are an abundant class of noncoding RNAs with widespread, cell/tissue specific pattern. Because of their involvement in the pathogenesis of multiple disease, they are receiving increasing attention. Previous work suggested that epigenetic features might be related to circRNA expression. However, current algorithms for circRNAs prediction neglect these features, leading to constant results across different cells. Here we built a machine learning framework named CIRCScan, to predict expression status and expression levels of circRNAs in various cell lines based on sequence and epigenetic features. Both expression status and expression levels can be accurately predicted by different groups of features. For expression status, the top features were similar in different cells. However, the top features for predicting expression levels were different in different cells. Noteworthy, the importance of H3K79me2 ranked high in predicting both circRNAs expression status and levels across different cells, indicating its important role in regulating circRNAs expression. Further validation experiment in K562 confirmed that knock down of H3K79me2 did result in reduction of circRNA production. Our study offers new insights into the regulation of circRNAs by incorporating epigenetic features in prediction models in different cellular contexts.

The genome-wide association studies (GWAS) have identified hundreds of susceptibility loci associated with autoimmune diseases. However, over 90% of risk variants are located in the noncoding regions, leading to great challenges in deciphering the underlying causal functional variants/genes and biological mechanisms. Previous studies focused on developing new scoring method to prioritize functional/disease-relevant variants. However, they principally incorporated annotation data across all cells/tissues while omitted the cell-specific or context-specific regulation. Moreover, limited analyses were performed to dissect the detailed molecular regulatory circuits linking functional GWAS variants to disease etiology. Here we devised a new analysis frame that incorporate hundreds of immune cell-specific multi-omics data to prioritize functional noncoding susceptibility SNPs with gene targets and further dissect their downstream molecular mechanisms and clinical applications for 19 autoimmune diseases. Most prioritized SNPs have genetic associations with transcription factors (TFs) binding, histone modification or chromatin accessibility, indicating their allelic regulatory roles on target genes. Their target genes were significantly enriched in immunologically related pathways and other immunologically related functions. We also detected long-range regulation on 90.7% of target genes including 132 ones exclusively regulated by distal SNPs (eg, CD28 , IL2RA ), which involves several potential key TFs (eg, CTCF), suggesting the important roles of long-range chromatin interaction in autoimmune diseases. Moreover, we identified hundreds of known or predicted druggable genes, and predicted some new potential drug targets for several autoimmune diseases, including two genes ( NFKB1 , SH2B3 ) with known drug indications on other diseases, highlighting their potential drug repurposing opportunities. In summary, our analyses may provide unique resource for future functional follow-up and drug application on autoimmune diseases, which are freely available at http://fngwas.online/.