IMPACT: Genomic Annotation of Cell-State-Specific Regulatory Elements Inferred from the Epigenome of Bound Transcription Factors

Enhancing Disease Risk Gene Discovery by Integrating Transcription Factor-Linked Trans-located Variants into Transcriptome-Wide Association Analyses

Transcriptome-wide association studies (TWAS) have been successful in identifying disease susceptibility genes by integrating cis-variants predicted gene expression with genome-wide association studies (GWAS) data. However, trans-located variants for predicting gene expression remain largely unexplored. Here, we introduce transTF-TWAS, which incorporates transcription factor (TF)-linked trans-located variants to enhance model building. Using data from the Genotype-Tissue Expression project, we predict gene expression and alternative splicing and applied these models to large GWAS datasets for breast, prostate, and lung cancers. We demonstrate that transTF-TWAS outperforms other existing TWAS approaches in both constructing gene prediction models and identifying disease-associated genes, as evidenced by simulations and real data analysis. Our transTF-TWAS approach significantly contributes to the discovery of disease risk genes. Findings from this study have shed new light on several genetically driven key regulators and their associated regulatory networks underlying disease susceptibility.

The regulation and differentiation of regulatory T cells and their dysfunction in autoimmune diseases

The discovery of FOXP3+ regulatory T (Treg) cells as a distinct cell lineage with a central role in regulating immune responses provided a deeper understanding of self-tolerance. The transcription factor FOXP3 serves a key role in Treg cell lineage determination and maintenance, but is not sufficient to enable the full potential of Treg cell suppression, indicating that other factors orchestrate the fine-tuning of Treg cell function. Moreover, FOXP3-independent mechanisms have recently been shown to contribute to Treg cell dysfunction. FOXP3 mutations in humans cause lethal fulminant systemic autoinflammation (IPEX syndrome). However, it remains unclear to what degree Treg cell dysfunction is contributing to the pathophysiology of common autoimmune diseases. In this Review, we discuss the origins of Treg cells in the periphery and the multilayered mechanisms by which Treg cells are induced, as well as the FOXP3-dependent and FOXP3-independent cellular programmes that maintain the suppressive function of Treg cells in humans and mice. Further, we examine evidence for Treg cell dysfunction in the context of common autoimmune diseases such as multiple sclerosis, inflammatory bowel disease, systemic lupus erythematosus and rheumatoid arthritis.

A statistical approach for identifying single nucleotide variants that affect transcription factor binding

Non-coding variants located within regulatory elements may alter gene expression by modifying transcription factor (TF) binding sites, thereby leading to functional consequences. Different TF models are being used to assess the effect of DNA sequence variants, such as single nucleotide variants (SNVs). Often existing methods are slow and do not assess statistical significance of results. We investigated the distribution of absolute maximal differential TF binding scores for general computational models that affect TF binding. We find that a modified Laplace distribution can adequately approximate the empirical distributions. A benchmark on in vitro and in vivo datasets showed that our approach improves upon an existing method in terms of performance and speed. Applications on eQTLs and on a genome-wide association study illustrate the usefulness of our statistics by highlighting cell type-specific regulators and target genes. An implementation of our approach is freely available on GitHub and as bioconda package.

The genetic basis of autoimmunity seen through the lens of T cell functional traits

Autoimmune disease heritability is enriched in T cell-specific regulatory regions of the genome. Modern-day T cell datasets now enable association studies between single nucleotide polymorphisms (SNPs) and a myriad of molecular phenotypes, including chromatin accessibility, gene expression, transcriptional programs, T cell antigen receptor (TCR) amino acid usage, and cell state abundances. Such studies have identified hundreds of quantitative trait loci (QTLs) in T cells that colocalize with genetic risk for autoimmune disease. The key challenge facing immunologists today lies in synthesizing these results toward a unified understanding of the autoimmune T cell: which genes, cell states, and antigens drive tissue destruction?

How can genetics analysis allow for early detection of rheumatoid arthritis

Early detection and diagnosis of rheumatoid arthritis (RA) is crucial for initiating appropriate therapy promptly, thereby leading to increase the likelihood of achieving remission and eventual cure. In this article, we will discuss the contribution of genetic information in predicting and detecting RA.

GWAS for systemic sclerosis identifies six novel susceptibility loci including one in the Fcγ receptor region

Here we report the largest Asian genome-wide association study (GWAS) for systemic sclerosis performed to date, based on data from Japanese subjects and comprising of 1428 cases and 112,599 controls. The lead SNP is in the FCGR/FCRL region, which shows a penetrating association in the Asian population, while a complete linkage disequilibrium SNP, rs10917688, is found in a cis-regulatory element for IRF8. IRF8 is also a significant locus in European GWAS for systemic sclerosis, but rs10917688 only shows an association in the presence of the risk allele of IRF8 in the Japanese population. Further analysis shows that rs10917688 is marked with H3K4me1 in primary B cells. A meta-analysis with a European GWAS detects 30 additional significant loci. Polygenic risk scores constructed with the effect sizes of the meta-analysis suggest the potential portability of genetic associations beyond populations. Prioritizing the top 5% of SNPs of IRF8 binding sites in B cells improves the fitting of the polygenic risk scores, underscoring the roles of B cells and IRF8 in the development of systemic sclerosis. The results also suggest that systemic sclerosis shares a common genetic architecture across populations.

Dynamic regulatory elements in single-cell multimodal data implicate key immune cell states enriched for autoimmune disease heritability

In autoimmune diseases such as rheumatoid arthritis, the immune system attacks the body's own cells. Developing a precise understanding of the cell states where noncoding autoimmune risk variants impart causal mechanisms is critical to developing curative therapies. Here, to identify noncoding regions with accessible chromatin that associate with cell-state-defining gene expression patterns, we leveraged multimodal single-nucleus RNA and assay for transposase-accessible chromatin (ATAC) sequencing data across 28,674 cells from the inflamed synovial tissue of 12 donors. Specifically, we used a multivariate Poisson model to predict peak accessibility from single-nucleus RNA sequencing principal components. For 14 autoimmune diseases, we discovered that cell-state-dependent ('dynamic') chromatin accessibility peaks in immune cell types were enriched for heritability, compared with cell-state-invariant ('cs-invariant') peaks. These dynamic peaks marked regulatory elements associated with T peripheral helper, regulatory T, dendritic and STAT1+CXCL10+ myeloid cell states. We argue that dynamic regulatory elements can help identify precise cell states enriched for disease-critical genetic variation.

Modeling tissue co-regulation estimates tissue-specific contributions to disease

Integrative analyses of genome-wide association studies and gene expression data have implicated many disease-critical tissues. However, co-regulation of genetic effects on gene expression across tissues impedes distinguishing biologically causal tissues from tagging tissues. In the present study, we introduce tissue co-regulation score regression (TCSC), which disentangles causal tissues from tagging tissues by regressing gene-disease association statistics (from transcriptome-wide association studies) on tissue co-regulation scores, reflecting correlations of predicted gene expression across genes and tissues. We applied TCSC to 78 diseases/traits (average n = 302,000) and gene expression prediction models for 48 GTEx tissues. TCSC identified 21 causal tissue-trait pairs at a 5% false discovery rate (FDR), including well-established findings, biologically plausible new findings (for example, aorta artery and glaucoma) and increased specificity of known tissue-trait associations (for example, subcutaneous adipose, but not visceral adipose, and high-density lipoprotein). TCSC also identified 17 causal tissue-trait covariance pairs at 5% FDR. In conclusion, TCSC is a precise method for distinguishing causal tissues from tagging tissues.

Integrative analyses for the identification of idiopathic pulmonary fibrosis-associated genes and shared loci with other diseases

BACKGROUND

Although genome-wide association studies (GWAS) have identified many genomic regions associated with idiopathic pulmonary fibrosis (IPF), the causal genes and functions remain largely unknown. Many single-cell expression data have become available for IPF, and there is increasing evidence suggesting a shared genetic basis between IPF and other diseases.

METHODS

We conducted integrative analyses to improve the power of GWAS. First, we calculated global and local genetic correlations to identify IPF genetically associated traits and local regions. Then, we prioritised candidate genes contributing to local genetic correlation. Second, we performed transcriptome-wide association analysis (TWAS) of 44 tissues to identify candidate genes whose genetically predicted expression level is associated with IPF. To replicate our findings and investigate the regulatory role of the transcription factors (TF) in identified candidate genes, we first conducted the heritability enrichment analysis in TF binding sites. Then, we examined the enrichment of the TF target genes in cell-type-specific differentially expressed genes (DEGs) identified from single-cell expression data of IPF and healthy lung samples.

FINDINGS

We identified 12 candidate genes across 13 genomic regions using local genetic correlation, including the POT1 locus (p value=0.00041), which contained variants with protective effects on lung cancer but increasing IPF risk. We identified another 13 novel genes using TWAS. Two TFs, MAFK and SMAD2, showed significant enrichment in both partitioned heritability and cell-type-specific DEGs.

INTERPRETATION

Our integrative analysis identified new genes for IPF susceptibility and expanded the understanding of the complex genetic architecture and disease mechanism of IPF.

Blood cell traits' GWAS loci colocalization with variation in PU.1 genomic occupancy prioritizes causal noncoding regulatory variants

Genome-wide association studies (GWASs) have uncovered numerous trait-associated loci across the human genome, most of which are located in noncoding regions, making interpretation difficult. Moreover, causal variants are hard to statistically fine-map at many loci because of widespread linkage disequilibrium. To address this challenge, we present a strategy utilizing transcription factor (TF) binding quantitative trait loci (bQTLs) for colocalization analysis to identify trait associations likely mediated by TF occupancy variation and to pinpoint likely causal variants using motif scores. We applied this approach to PU.1 bQTLs in lymphoblastoid cell lines and blood cell trait GWAS data. Colocalization analysis revealed 69 blood cell trait GWAS loci putatively driven by PU.1 occupancy variation. We nominate PU.1 motif-altering variants as the likely shared causal variants at 51 loci. Such integration of TF bQTL data with other GWAS data may reveal transcriptional regulatory mechanisms and causal noncoding variants underlying additional complex traits.

Genetics of skeletal proportions in two different populations

Human height can be divided into sitting height and leg length, reflecting growth of different parts of the skeleton whose relative proportions are captured by the ratio of sitting to total height (as sitting height ratio, SHR). Height is a highly heritable trait, and its genetic basis has been well-studied. However, the genetic determinants of skeletal proportion are much less well-characterized. Expanding substantially on past work, we performed a genome-wide association study (GWAS) of SHR in ∼450,000 individuals with European ancestry and ∼100,000 individuals with East Asian ancestry from the UK and China Kadoorie Biobanks. We identified 565 loci independently associated with SHR, including all genomic regions implicated in prior GWAS in these ancestries. While SHR loci largely overlap height-associated loci (P < 0.001), the fine-mapped SHR signals were often distinct from height. We additionally used fine-mapped signals to identify 36 credible sets with heterogeneous effects across ancestries. Lastly, we used SHR, sitting height, and leg length to identify genetic variation acting on specific body regions rather than on overall human height.

Colocalization of blood cell traits GWAS associations and variation in PU.1 genomic occupancy prioritizes causal noncoding regulatory variants

Genome-wide association studies (GWAS) have uncovered numerous trait-associated loci across the human genome, most of which are located in noncoding regions, making interpretations difficult. Moreover, causal variants are hard to statistically fine-map at many loci because of widespread linkage disequilibrium. To address this challenge, we present a strategy utilizing transcription factor (TF) binding quantitative trait loci (bQTLs) for colocalization analysis to identify trait associations likely mediated by TF occupancy variation and to pinpoint likely causal variants using motif scores. We applied this approach to PU.1 bQTLs in lymphoblastoid cell lines and blood cell traits GWAS data. Colocalization analysis revealed 69 blood cell trait GWAS loci putatively driven by PU.1 occupancy variation. We nominate PU.1 motif-altering variants as the likely shared causal variants at 51 loci. Such integration of TF bQTL data with other GWAS data may reveal transcriptional regulatory mechanisms and causal noncoding variants underlying additional complex traits.

Control of naive and effector CD4 T cell receptor repertoires by rheumatoid-arthritis-risk HLA alleles

OBJECTIVE

Human Leukocyte Antigen (HLA) alleles regulate susceptibility to rheumatoid arthritis (RA) and immune-mediated diseases. This study aims to elucidate the impact of HLA alleles to T cell subsets.

METHODS

We performed genome-wide and HLA allele association analysis for T cell receptor (TCR) beta chain repertoire in 13 purified T cell subsets from the ImmuNexUT database, consisting of 407 donors with ten immune-mediated diseases and healthy controls.

RESULTS

HLA class II alleles were associated with TRBV gene usage and the public clones of CD4 T cells, while HLA class I alleles were associated with CD8 T cells. RA-risk and immune-mediated diseases-risk HLA alleles were associated with TRBV gene usage of naive and effector CD4 T cell subsets and public clones accumulating in Th17. Clonal diversity was independent of HLA alleles and was correlated with transcriptome changes that reflect TCR signaling.

CONCLUSION

This study revealed in vivo evidence that both HLA alleles and environmental factors shape naive and effector TCR repertoires in RA and immune-mediated diseases patients.

Integrating transcription factor occupancy with transcriptome-wide association analysis identifies susceptibility genes in human cancers

Transcriptome-wide association studies (TWAS) have successfully discovered many putative disease susceptibility genes. However, TWAS may suffer from inaccuracy of gene expression predictions due to inclusion of non-regulatory variants. By integrating prior knowledge of susceptible transcription factor occupied elements, we develop sTF-TWAS and demonstrate that it outperforms existing TWAS approaches in both simulation and real data analyses. Under the sTF-TWAS framework, we build genetic models to predict alternative splicing and gene expression in normal breast, prostate and lung tissues from the Genotype-Tissue Expression project and apply these models to data from large genome-wide association studies (GWAS) conducted among European-ancestry populations. At Bonferroni-corrected P < 0.05, we identify 354 putative susceptibility genes for these cancers, including 189 previously unreported in GWAS loci and 45 in loci unreported by GWAS. These findings provide additional insight into the genetic susceptibility of human cancers. Additionally, we show the generalizability of the sTF-TWAS on non-cancer diseases.

Deep learning approaches for noncoding variant prioritization in neurodegenerative diseases

Determining how noncoding genetic variants contribute to neurodegenerative dementias is fundamental to understanding disease pathogenesis, improving patient prognostication, and developing new clinical treatments. Next generation sequencing technologies have produced vast amounts of genomic data on cell type-specific transcription factor binding, gene expression, and three-dimensional chromatin interactions, with the promise of providing key insights into the biological mechanisms underlying disease. However, this data is highly complex, making it challenging for researchers to interpret, assimilate, and dissect. To this end, deep learning has emerged as a powerful tool for genome analysis that can capture the intricate patterns and dependencies within these large datasets. In this review, we organize and discuss the many unique model architectures, development philosophies, and interpretation methods that have emerged in the last few years with a focus on using deep learning to predict the impact of genetic variants on disease pathogenesis. We highlight both broadly-applicable genomic deep learning methods that can be fine-tuned to disease-specific contexts as well as existing neurodegenerative disease research, with an emphasis on Alzheimer's-specific literature. We conclude with an overview of the future of the field at the intersection of neurodegeneration, genomics, and deep learning.

Multi-ancestry genome-wide association analyses identify novel genetic mechanisms in rheumatoid arthritis

Rheumatoid arthritis (RA) is a highly heritable complex disease with unknown etiology. Multi-ancestry genetic research of RA promises to improve power to detect genetic signals, fine-mapping resolution and performances of polygenic risk scores (PRS). Here, we present a large-scale genome-wide association study (GWAS) of RA, which includes 276,020 samples from five ancestral groups. We conducted a multi-ancestry meta-analysis and identified 124 loci (P < 5 × 10-8), of which 34 are novel. Candidate genes at the novel loci suggest essential roles of the immune system (for example, TNIP2 and TNFRSF11A) and joint tissues (for example, WISP1) in RA etiology. Multi-ancestry fine-mapping identified putatively causal variants with biological insights (for example, LEF1). Moreover, PRS based on multi-ancestry GWAS outperformed PRS based on single-ancestry GWAS and had comparable performance between populations of European and East Asian ancestries. Our study provides several insights into the etiology of RA and improves the genetic predictability of RA.

Rebooting Regulatory T Cell and Dendritic Cell Function in Immune-Mediated Inflammatory Diseases: Biomarker and Therapy Discovery under a Multi-Omics Lens

Immune-mediated inflammatory diseases (IMIDs) are a group of autoimmune and chronic inflammatory disorders with constantly increasing prevalence in the modern world. The vast majority of IMIDs develop as a consequence of complex mechanisms dependent on genetic, epigenetic, molecular, cellular, and environmental elements, that lead to defects in immune regulatory guardians of tolerance, such as dendritic (DCs) and regulatory T (Tregs) cells. As a result of this dysfunction, immune tolerance collapses and pathogenesis emerges. Deeper understanding of such disease driving mechanisms remains a major challenge for the prevention of inflammatory disorders. The recent renaissance in high throughput technologies has enabled the increase in the amount of data collected through multiple omics layers, while additionally narrowing the resolution down to the single cell level. In light of the aforementioned, this review focuses on DCs and Tregs and discusses how multi-omics approaches can be harnessed to create robust cell-based IMID biomarkers in hope of leading to more efficient and patient-tailored therapeutic interventions.

Virtual ChIP-seq: predicting transcription factor binding by learning from the transcriptome

Existing methods for computational prediction of transcription factor (TF) binding sites evaluate genomic regions with similarity to known TF sequence preferences. Most TF binding sites, however, do not resemble known TF sequence motifs, and many TFs are not sequence-specific. We developed Virtual ChIP-seq, which predicts binding of individual TFs in new cell types, integrating learned associations with gene expression and binding, TF binding sites from other cell types, and chromatin accessibility data in the new cell type. This approach outperforms methods that predict TF binding solely based on sequence preference, predicting binding for 36 TFs (MCC>0.3).

Dissection of multiple sclerosis genetics identifies B and CD4+ T cells as driver cell subsets

Background

Multiple sclerosis (MS) is an autoimmune condition of the central nervous system with a well-characterized genetic background. Prior analyses of MS genetics have identified broad enrichments across peripheral immune cells, yet the driver immune subsets are unclear.

Results

We utilize chromatin accessibility data across hematopoietic cells to identify cell type-specific enrichments of MS genetic signals. We find that CD4 T and B cells are independently enriched for MS genetics and further refine the driver subsets to T_h17 and memory B cells, respectively. We replicate our findings in data from untreated and treated MS patients and find that immunomodulatory treatments suppress chromatin accessibility at driver cell types. Integration of statistical fine-mapping and chromatin interactions nominate numerous putative causal genes, illustrating complex interplay between shared and cell-specific genes.

Conclusions

Overall, our study finds that open chromatin regions in CD4 T cells and B cells independently drive MS genetic signals. Our study highlights how careful integration of genetics and epigenetics can provide fine-scale insights into causal cell types and nominate new genes and pathways for disease.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02694-y.

Single-cell eQTL models reveal dynamic T cell state dependence of disease loci

Non-coding genetic variants may cause disease by modulating gene expression. However, identifying these expression quantitative trait loci (eQTLs) is complicated by differences in gene regulation across fluid functional cell states within cell types. These states-for example, neurotransmitter-driven programs in astrocytes or perivascular fibroblast differentiation-are obscured in eQTL studies that aggregate cells1,2. Here we modelled eQTLs at single-cell resolution in one complex cell type: memory T cells. Using more than 500,000 unstimulated memory T cells from 259 Peruvian individuals, we show that around one-third of 6,511 cis-eQTLs had effects that were mediated by continuous multimodally defined cell states, such as cytotoxicity and regulatory capacity. In some loci, independent eQTL variants had opposing cell-state relationships. Autoimmune variants were enriched in cell-state-dependent eQTLs, including risk variants for rheumatoid arthritis near ORMDL3 and CTLA4; this indicates that cell-state context is crucial to understanding potential eQTL pathogenicity. Moreover, continuous cell states explained more variation in eQTLs than did conventional discrete categories, such as CD4+ versus CD8+, suggesting that modelling eQTLs and cell states at single-cell resolution can expand insight into gene regulation in functionally heterogeneous cell types.

Immune disease variants modulate gene expression in regulatory CD4+ T cells

Identifying cellular functions dysregulated by disease-associated variants could implicate novel pathways for drug targeting or modulation in cell therapies. However, follow-up studies can be challenging if disease-relevant cell types are difficult to sample. Variants associated with immune diseases point toward the role of CD4+ regulatory T cells (Treg cells). We mapped genetic regulation (quantitative trait loci [QTL]) of gene expression and chromatin activity in Treg cells, and we identified 133 colocalizing loci with immune disease variants. Colocalizations of immune disease genome-wide association study (GWAS) variants with expression QTLs (eQTLs) controlling the expression of CD28 and STAT5A, involved in Treg cell activation and interleukin-2 (IL-2) signaling, support the contribution of Treg cells to the pathobiology of immune diseases. Finally, we identified seven known drug targets suitable for drug repurposing and suggested 63 targets with drug tractability evidence among the GWAS signals that colocalized with Treg cell QTLs. Our study is the first in-depth characterization of immune disease variant effects on Treg cell gene expression modulation and dysregulation of Treg cell function.

Fine-mapping of Parkinson's disease susceptibility loci identifies putative causal variants

Recent genome-wide association studies have identified 78 loci associated with Parkinson's disease susceptibility but the underlying mechanisms remain largely unclear. To identify likely causal variants for disease risk, we fine-mapped these Parkinson's-associated loci using four different fine-mapping methods. We then integrated multi-assay cell type-specific epigenomic profiles to pinpoint the likely mechanism of action of each variant, allowing us to identify Consensus single nucleotide polymorphism (SNPs) that disrupt LRRK2 and FCGR2A regulatory elements in microglia, an MBNL2 enhancer in oligodendrocytes, and a DYRK1A enhancer in neurons. This genome-wide functional fine-mapping investigation of Parkinson's disease substantially advances our understanding of the causal mechanisms underlying this complex disease while avoiding focus on spurious, non-causal mechanisms. Together, these results provide a robust, comprehensive list of the likely causal variants, genes and cell-types underlying Parkinson's disease risk as demonstrated by consistently greater enrichment of our fine-mapped SNPs relative to lead GWAS SNPs across independent functional impact annotations. In addition, our approach prioritized an average of 3/85 variants per locus as putatively causal, making downstream experimental studies both more tractable and more likely to yield disease-relevant, actionable results. Large-scale studies comparing individuals with Parkinson's disease to age-matched controls have identified many regions of the genome associated with the disease. However, there is widespread correlation between different parts of the genome, making it difficult to tell which genetic variants cause Parkinson's and which are simply co-inherited with causal variants. We therefore applied a suite of statistical models to identify the most likely causal genetic variants (i.e. fine-mapping). We then linked these genetic variants with epigenomic and gene expression signatures across a wide variety of tissues and cell types to identify how these variants cause disease. Therefore, this study provides a comprehensive and robust list of cellular and molecular mechanisms that may serve as targets in the development of more effective Parkinson's therapeutics.

Fine mapping with epigenetic information and 3D structure

Since 2005, thousands of genome-wide association studies (GWAS) have been published, identifying hundreds of thousands of genetic variants that increase risk of complex traits such as autoimmune diseases. This wealth of data has the potential to improve patient care, through personalized medicine and the identification of novel drug targets. However, the potential of GWAS for clinical translation has not been fully achieved yet, due to the fact that the functional interpretation of risk variants and the identification of causal variants and genes are challenging. The past decade has seen the development of great advances that are facilitating the overcoming of these limitations, by utilizing a plethora of genomics and epigenomics tools to map and characterize regulatory elements and chromatin interactions, which can be used to fine map GWAS loci, and advance our understanding of the biological mechanisms that cause disease.

echolocatoR: an automated end-to-end statistical and functional genomic fine-mapping pipeline

SUMMARY

echolocatoR integrates a diverse suite of statistical and functional fine-mapping tools to identify, test enrichment in, and visualize high-confidence causal consensus variants in any phenotype. It requires minimal input from users (a summary statistics file), can be run in a single R function, and provides extensive access to relevant datasets (e.g. reference linkage disequilibrium panels, quantitative trait loci, genome-wide annotations, cell-type-specific epigenomics), thereby enabling rapid, robust and scalable end-to-end fine-mapping investigations.

AVAILABILITY AND IMPLEMENTATION

echolocatoR is an open-source R package available through GitHub under the GNU General Public License (Version 3) license: https://github.com/RajLabMSSM/echolocatoR.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Chondroprotective effects of CDK4/6 inhibition via enhanced ubiquitin-dependent degradation of JUN in synovial fibroblasts

OBJECTIVE

Targeting synovial fibroblasts (SF) using a cyclin-dependent kinase (CDK) 4/6 inhibitor (CDKI) could be a potent therapy for rheumatoid arthritis (RA) via inhibition of proliferation and MMP-3 production. This study was designed to elucidate the mechanism of chondroprotective effects on SFs by CDK 4/6 inhibition.

METHODS

CDK4/6 activity was inhibited using CDKI treatment or enhanced by adenoviral gene transduction. Chondroprotective effects were evaluated using a collagen induced arthritis model (CIA). Gene and protein expression were evaluated with quantitative PCR, ELISA, and Western blotting. The binding of nuclear extracts to DNA was assessed with an electrophoresis mobility shift assay. RNA-Seq was performed to identify gene sets affected by CDKI treatment.

RESULTS

CDKI attenuated cartilage destruction and MMP-3 production in CIA. In RASFs, CDKI impaired the binding of AP-1 components to DNA and inhibited the production of MMP-1 and MMP-3, which contain the AP-1 binding sequence in their promoter. CDK4/6 protected JUN from proteasome-dependent degradation by inhibiting ubiquitination. The RNA-Seq analysis identified CDKI-sensitive inflammatory genes, which were associated with the pathway of RA-associated genes, cytokine-cytokine receptor interaction, and IL-17 signalling. Notably, the AP-1 motif was enriched in these genes.

CONCLUSION

The mechanism of chondroprotective effects by CDK4/6 inhibition was achieved by the attenuation of AP-1 transcriptional activity via the impaired stability of JUN. Since the pharmacologic inhibition of CDK4/6 has been established as tolerable in cancer treatment, it could also be beneficial in patients with RA due to its chondroprotective and anti-inflammatory effects.

Multi-ancestral GWAS identifies shared and Asian-specific loci for SLE and links type III interferon signaling and lysosomal function to the disease

OBJECTIVES

Systemic lupus erythematosus (SLE) is a prototype autoimmune disease with differences in prevalence and severity among ancestral groups. This study aims to identify novel genetic components either shared or distinct between Asian and European populations.

METHODS

Both trans-ancestral and ancestry-specific meta-analyses of genome-wide association studies (GWAS) for SLE were performed, involving 30,604 participants of European, Chinese or Thai origin. Using public epigenomic data and expression quantitative trait loci, fine-mapping analyses were conducted to identify putative causal variants and genes for the newly identified loci. Performance of polygenic risk scores (PRS) for the Thai cohort was evaluated comparing different training data.

RESULTS

We identified ten novel SLE susceptibility loci, four of which were found by Asian-specific meta-analyses. A 1bp deletion upstream of IFNLR1 was found associated with SLE, with the risk allele correlated with increased expression of IFNLR1. This gene encodes interferon lambda receptor 1, pointing to the role of type III interferon signaling in SLE. An intronic variant in SLC29A3 was found associated with SLE only in Asians. The putative risk variant may modulate SLC29A3 expression in a monocyte-specific manner. SLC29A3 encodes a lysosomal nucleoside transporter, and subsequent analyses suggest reduced lysosomal function and phagocytosis might be the mechanism underlying this association. In addition, trans-ancestral meta-analysis was proved to be valuable in risk prediction for individuals without ancestry-matched data.

CONCLUSION

Multi-ancestral GWAS identified both shared and Asian-specific loci for SLE, and functional annotation pointed to the involvement of increased type III interferon signaling and reduced lysosomal function in SLE.

Beyond GWAS: from simple associations to functional insights

Each human, when born, has slightly different DNA sequences, which make each of us unique. The variations in DNA sequences are called genetic variants. The primary aim of genome-wide association study (GWAS) is to detect associations between genetic variants and human phenotypes. Since GWAS focuses on germ-line variants, there is no reverse causation. Therefore, GWAS is one of the few tools that can assess the causality of human diseases. In the past 10 years, many large-scale GWAS have been conducted. Although the primary outputs of GWAS are just a series of statistics, its downstream analyses provided many insights beyond simple associations: the causal mechanisms for autoimmune diseases and shared etiology between diseases. Moreover, GWAS downstream analyses generated scores potentially helpful in predicting clinical outcomes of each patient. This review focuses on GWAS for autoimmune diseases and introduces significant achievements of its downstream analyses. We also provide future directions that potentially overcome current limitations. We restrict our discussion to common autoimmune diseases (e.g., rheumatoid arthritis) since rare Mendelian diseases possess distinct genetic etiologies and are not tested by GWAS.

Motif-Raptor: a cell type-specific and transcription factor centric approach for post-GWAS prioritization of causal regulators

MOTIVATION

Genome-wide association studies (GWASs) have identified thousands of common trait-associated genetic variants but interpretation of their function remains challenging. These genetic variants can overlap the binding sites of transcription factors (TFs) and therefore could alter gene expression. However, we currently lack a systematic understanding on how this mechanism contributes to phenotype.

RESULTS

We present Motif-Raptor, a TF-centric computational tool that integrates sequence-based predictive models, chromatin accessibility, gene expression datasets and GWAS summary statistics to systematically investigate how TF function is affected by genetic variants. Given trait-associated non-coding variants, Motif-Raptor can recover relevant cell types and critical TFs to drive hypotheses regarding their mechanism of action. We tested Motif-Raptor on complex traits such as rheumatoid arthritis and red blood cell count and demonstrated its ability to prioritize relevant cell types, potential regulatory TFs and non-coding SNPs which have been previously characterized and validated.

AVAILABILITY AND IMPLEMENTATION

Motif-Raptor is freely available as a Python package at: https://github.com/pinellolab/MotifRaptor.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome editing to define the function of risk loci and variants in rheumatic disease

Discoveries in human genetic studies have revolutionized our understanding of complex rheumatic and autoimmune diseases, including the identification of hundreds of genetic loci and single nucleotide polymorphisms that potentially predispose individuals to disease. However, in most cases, the exact disease-causing variants and their mechanisms of action remain unresolved. Functional follow-up of these findings is most challenging for genomic variants that are in non-coding genomic regions, where the large majority of common disease-associated variants are located, and/or that probably affect disease progression via cell type-specific gene regulation. To deliver on the therapeutic promise of human genetic studies, defining the mechanisms of action of these alleles is essential. Genome editing technology, such as CRISPR-Cas, has created a vast toolbox for targeted genetic and epigenetic modifications that presents unprecedented opportunities to decipher disease-causing loci, genes and variants in autoimmunity. In this Review, we discuss the past 5-10 years of progress in resolving the mechanisms underlying rheumatic disease-associated alleles, with an emphasis on how genomic editing techniques can enable targeted dissection and mechanistic studies of causal autoimmune risk variants.

Deep neural networks identify sequence context features predictive of transcription factor binding

Transcription factors (TFs) bind DNA by recognizing specific sequence motifs, typically of length 6-12bp. A motif can occur many thousands of times in the human genome, but only a subset of those sites are actually bound. Here we present a machine learning framework leveraging existing convolutional neural network architectures and model interpretation techniques to identify and interpret sequence context features most important for predicting whether a particular motif instance will be bound. We apply our framework to predict binding at motifs for 38 TFs in a lymphoblastoid cell line, score the importance of context sequences at base-pair resolution, and characterize context features most predictive of binding. We find that the choice of training data heavily influences classification accuracy and the relative importance of features such as open chromatin. Overall, our framework enables novel insights into features predictive of TF binding and is likely to inform future deep learning applications to interpret non-coding genetic variants.

Improving the trans-ancestry portability of polygenic risk scores by prioritizing variants in predicted cell-type-specific regulatory elements

Poor trans-ancestry portability of polygenic risk scores is a consequence of Eurocentric genetic studies and limited knowledge of shared causal variants. Leveraging regulatory annotations may improve portability by prioritizing functional over tagging variants. We constructed a resource of 707 cell-type-specific IMPACT regulatory annotations by aggregating 5,345 epigenetic datasets to predict binding patterns of 142 transcription factors across 245 cell types. We then partitioned the common SNP heritability of 111 genome-wide association study summary statistics of European (average n ≈ 189,000) and East Asian (average n ≈ 157,000) origin. IMPACT annotations captured consistent SNP heritability between populations, suggesting prioritization of shared functional variants. Variant prioritization using IMPACT resulted in increased trans-ancestry portability of polygenic risk scores from Europeans to East Asians across all 21 phenotypes analyzed (49.9% mean relative increase in R2). Our study identifies a crucial role for functional annotations such as IMPACT to improve the trans-ancestry portability of genetic data.

Functional genomics in autoimmune diseases

Associations between genetic loci and increased susceptibility to autoimmune disease have been well characterized, however, translating this knowledge into mechanistic insight and patient benefit remains a challenge. While improvements in the precision, completeness and accuracy of our genetic understanding of autoimmune diseases will undoubtedly be helpful, meeting this challenge will require two interlinked problems to be addressed: first which of the highly correlated variants at an individual locus is responsible for increased disease risk, and second what are the downstream effects of this variant. Given that the majority of loci are thought to affect non-coding regulatory elements, the second question is often reframed as what are the target gene(s) and pathways affected by causal variants. Currently, these questions are being addressed using a wide variety of novel techniques and datasets. In many cases, these approaches are complementary and it is likely that the most accurate picture will be generated by consolidating information relating to transcription, regulatory activity, chromatin accessibility, chromatin conformation and readouts from functional experiments, such as genome editing and reporter assays. It is clear that it will be necessary to gather this information from disease relevant cell types and conditions and that by doing so our understanding of disease etiology will be improved. This review is focused on the field of autoimmune disease functional genomics with a particular focus on the most exciting and significant research to be published within the last couple of years.

Advances in genetics toward identifying pathogenic cell states of rheumatoid arthritis

Rheumatoid arthritis (RA) risk has a large genetic component (~60%) that is still not fully understood. This has hampered the design of effective treatments that could promise lifelong remission. RA is a polygenic disease with 106 known genome-wide significant associated loci and thousands of small effect causal variants. Our current understanding of RA risk has suggested cell-type-specific contexts for causal variants, implicating CD4 + effector memory T cells, as well as monocytes, B cells and stromal fibroblasts. While these cellular states and categories are still mechanistically broad, future studies may identify causal cell subpopulations. These efforts are propelled by advances in single cell profiling. Identification of causal cell subpopulations may accelerate therapeutic intervention to achieve lifelong remission.

Early progression to active tuberculosis is a highly heritable trait driven by 3q23 in Peruvians

Of the 1.8 billion people worldwide infected with Mycobacterium tuberculosis, 5-15% will develop active tuberculosis (TB). Approximately half will progress to active TB within the first 18 months after infection, presumably because they fail to mount an effective initial immune response. Here, in a genome-wide genetic study of early TB progression, we genotype 4002 active TB cases and their household contacts in Peru. We quantify genetic heritability ([Formula: see text]) of early TB progression to be 21.2% (standard error 0.08). This suggests TB progression has a strong genetic basis, and is comparable to traits with well-established genetic bases. We identify a novel association between early TB progression and variants located in a putative enhancer region on chromosome 3q23 (rs73226617, OR = 1.18; P = 3.93 × 10-8). With in silico and in vitro analyses we identify rs73226617 or rs148722713 as the likely functional variant and ATP1B3 as a potential causal target gene with monocyte specific function.