A landscape of gene expression regulation for synovium in arthritis.

The synovium is an important component of any synovial joint and is the major target tissue of inflammatory arthritis. However, the multi-omics landscape of synovium required for functional inference is absent from large-scale resources. Here we integrate genomics with transcriptomics and chromatin accessibility features of human synovium in up to 245 arthritic patients, to characterize the landscape of genetic regulation on gene expression and the regulatory mechanisms mediating arthritic diseases predisposition. We identify 4765 independent primary and 616 secondary cis-expression quantitative trait loci (cis-eQTLs) in the synovium and find that the eQTLs with multiple independent signals have stronger effects and heritability than single independent eQTLs. Integration of genome-wide association studies (GWASs) and eQTLs identifies 84 arthritis related genes, revealing 38 novel genes which have not been reported by previous studies using eQTL data from the GTEx project or immune cells. We further develop a method called eQTac to identify variants that could affect gene expression by affecting chromatin accessibility and identify 1517 regions with potential regulatory function of chromatin accessibility. Altogether, our study provides a comprehensive synovium multi-omics resource for arthritic diseases and gains new insights into the regulation of gene expression.

Enhancer variants on chromosome 2p14 regulating SPRED2 and ACTR2 act as a signal amplifier to protect against rheumatoid arthritis.

Genome-wide association studies (GWASs) have repeatedly reported multiple non-coding single-nucleotide polymorphisms (SNPs) at 2p14 associated with rheumatoid arthritis (RA), but their functional roles in the pathological mechanisms of RA remain to be explored. In this study, we integrated a series of bioinformatics and functional experiments and identified three intronic RA SNPs (rs1876518, rs268131, and rs2576923) within active enhancers that can regulate the expression of SPRED2 directly. At the same time, SPRED2 and ACTR2 influence each other as a positive feedback signal amplifier to strengthen the protective role in RA by inhibiting the migration and invasion of rheumatoid fibroblast-like synoviocytes (FLSs). In particular, the transcription factor CEBPB preferentially binds to the rs1876518-T allele to increase the expression of SPRED2 in FLSs. Our findings decipher the molecular mechanisms behind the GWAS signals at 2p14 for RA and emphasize SPRED2 as a potential candidate gene for RA, providing a potential target and direction for precise treatment of RA.

Individualized pathway activity algorithm identifies oncogenic pathways in pan-cancer analysis.

BACKGROUND: Accumulative evidences have shown that dysregulation of biological pathways contributed to the initiation and progression of malignant tumours. Several methods for pathway activity measurement have been proposed, but they are restricted to making comparisons between groups or sensitive to experimental batch effects. METHODS: We introduced a novel method for individualized pathway activity measurement (IPAM) that is based on the ranking of gene expression levels in individual sample. Taking advantage of IPAM, we calculated the pathway activity of 318 pathways from KEGG database in the 10528 tumour/normal samples of 33 cancer types from TCGA to identify characteristic dysregulated pathways among different cancer types. FINDINGS: IPAM precisely quantified the level of activity of each pathway in pan-cancer analysis and exhibited better performance in cancer classification and prognosis prediction over five widely used tools. The average ROC-AUC of cancer diagnostic model using tumour-educated platelets (TEPs) reached 92.84%, suggesting the potential of our algorithm in early diagnosis of cancer. We identified several pathways significantly deregulated and associated with patient survival in a large fraction of cancer types, such as tyrosine metabolism, fatty acid degradation, cell cycle, p53 signalling pathway and DNA replication. We also confirmed the dominant role of metabolic pathways in cancer pathway dysregulation and identified the driving factors of specific pathway dysregulation, such as PPARA for branched-chain amino acid metabolism and NR1I2, NR1I3 for fatty acid metabolism. INTERPRETATION: Our study will provide novel clues for understanding the pathological mechanisms of cancer, ultimately paving the way for personalized medicine of cancer. FUNDING: A full list of funding can be found in the Acknowledgements section.

Transcriptome-wide association study identifies multiple genes and pathways associated with thyroid function.

Thyroid dysfunction is a common endocrine disease measured by thyroid-stimulating hormone (TSH) level. Although >70 genetic loci associated with TSH have been reported through genome-wide association studies (GWASs), the variants can only explain a small fraction of the thyroid function heritability. To identify novel candidate genes for thyroid function, we conducted the first large-scale transcriptome-wide association study (TWAS) for thyroid function using GWAS-summary data for TSH levels in up to 119 715 individuals combined with precomputed gene expression weights of six panels from four tissue types. The candidate genes identified by TWAS were further validated by TWAS replication and gene expression profiles. We identified 74 conditionally independent genes significantly associated with thyroid function, such as PDE8B (P = 1.67 × 10-282), PDE10A (P = 7.61 × 10-119), NR3C2 (P = 1.50 × 10-92) and CAPZB (P = 3.13 × 10-79). After TWAS replication using UKBB datasets, 26 genes were replicated for significant associations with thyroid-relevant diseases/traits. Among them, 16 genes were causal for their associations to thyroid-relevant diseases/traits and further validated in differential expression analyses, including two novel genes (MFSD6 and RBM47) that did not implicate in previous GWASs. Enrichment analyses detected several pathways associated with thyroid function, such as the cAMP signaling pathway (P = 7.27 × 10-4), hemostasis (P = 3.74 × 10-4), and platelet activation, signaling and aggregation (P = 9.98 × 10-4). Our study identified multiple candidate genes and pathways associated with thyroid function, providing novel clues for revealing the genetic mechanisms of thyroid function and disease.

Modeling circRNA expression pattern with integrated sequence and epigenetic features demonstrates the potential involvement of H3K79me2 in circRNA expression.

MOTIVATION: CircRNAs are an abundant class of non-coding RNAs with widespread, cell-/tissue-specific patterns. Previous work suggested that epigenetic features might be related to circRNA expression. However, the contribution of epigenetic changes to circRNA expression has not been investigated systematically. Here, we built a machine learning framework named CIRCScan, to predict circRNA expression in various cell lines based on the sequence and epigenetic features. RESULTS: The predicted accuracy of the expression status models was high with area under the curve of receiver operating characteristic (ROC) values of 0.89-0.92 and the false-positive rates of 0.17-0.25. Predicted expressed circRNAs were further validated by RNA-seq data. The performance of expression-level prediction models was also good with normalized root-mean-square errors of 0.28-0.30 and Pearson's correlation coefficient r over 0.4 in all cell lines, along with Spearman's correlation coefficient ρ of 0.33-0.46. Noteworthy, H3K79me2 was highly ranked in modeling both circRNA expression status and levels across different cells. Further analysis in additional nine cell lines demonstrated a significant enrichment of H3K79me2 in circRNA flanking intron regions, supporting the potential involvement of H3K79me2 in circRNA expression regulation. AVAILABILITY AND IMPLEMENTATION: The CIRCScan assembler is freely available online for academic use at https://github.com/johnlcd/CIRCScan. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An Allele-Specific Functional SNP Associated with Two Systemic Autoimmune Diseases Modulates IRF5 Expression by Long-Range Chromatin Loop Formation.

Both systemic lupus erythematosus (SLE) and systemic sclerosis (SSc) are autoimmune diseases sharing similar genetic backgrounds. Genome-wide association studies 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 a series of bioinformatics and functional analyses, we prioritized a potential independent functional single-nucleotide polymorphism (rs13239597) within TNPO3 promoter region, residing in a putative enhancer element and validated that IRF5 is the distal target gene (∼118 kb) of rs13239597, which is a key regulator involved in pathogenic autoantibody dysregulation, increasing risk of both SLE and SSc. We experimentally validated the long-range chromatin interactions between rs13239597 and IRF5 using chromosome conformation capture 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 a mechanistic insight of a noncoding functional variant acting as an allele-specific distal enhancer to directly modulate IRF5 expression, which might obligate in understanding of complex genetic architectures of SLE and SSc pathogenesis.

Multiple Functional Variants at 13q14 Risk Locus for Osteoporosis Regulate RANKL Expression Through Long-Range Super-Enhancer.

RANKL is a key regulator involved in bone metabolism, and a drug target for osteoporosis. The clinical diagnosis and assessment of osteoporosis are mainly based on bone mineral density (BMD). Previous powerful genomewide association studies (GWASs) have identified multiple intergenic single-nucleotide polymorphisms (SNPs) located over 100 kb upstream of RANKL and 65 kb downstream of AKAP11 at 13q14.11 for osteoporosis. Whether these SNPs exert their roles on osteoporosis through RANKL is unknown. In this study, we conducted integrative analyses combining expression quantitative trait locus (eQTL), genomic chromatin interaction (high-throughput chromosome conformation capture [Hi-C]), epigenetic annotation, and a series of functional assays. The eQTL analysis identified six potential functional SNPs (rs9533090, rs9594738, r8001611, rs9533094, rs9533095, and rs9594759) exclusively correlated with RANKL gene expression (p < 0.001) at 13q14.11. Co-localization analyses suggested that eQTL signal for RANKL and BMD-GWAS signal shared the same causal variants. Hi-C analysis and functional annotation further validated that the first five osteoporosis SNPs are located in a super-enhancer region to regulate the expression of RANKL via long-range chromosomal interaction. Particularly, dual-luciferase assay showed that the region harboring rs9533090 in the super-enhancer has the strongest enhancer activity, and rs9533090 is an allele-specific regulatory SNP. Furthermore, deletion of the region harboring rs9533090 using CRISPR/Cas9 genome editing significantly reduced RANKL expression in both mRNA level and protein level. Finally, we found that the rs9533090-C robustly recruits transcription factor NFIC, which efficiently elevates the enhancer activity and increases the RANKL expression. In summary, we provided a feasible method to identify regulatory noncoding SNPs to distally regulate their target gene underlying the pathogenesis of osteoporosis by using bioinformatics data analyses and experimental validation. Our findings would be a potential and promising therapeutic target for precision medicine in osteoporosis. © 2018 American Society for Bone and Mineral Research.