Dynamic landscape of immune cell-specific gene regulation in immune-mediated diseases

Mechanisms of age-related Treg dysfunction in an arthritic environment

Rheumatoid arthritis (RA) is an autoimmune disease characterized by a polyarticular synovitis. In recent years, elderly onset rheumatoid arthritis (EORA) has been increasing. Treg cells in RA have been reported to be dysfunctional, but the relationship between aging and their functional changes is unclear. Here, we found that Treg cells from EORA patients had increased percentages, but decreased activity compared to those from younger onset RA (YORA) patients. In experiments using arthritis model mice, decreased suppressive function and oxygen consumption rate (OCR) were observed in Treg cells only from old arthritic mice. Furthermore, type I interferon (IFN) signaling was upregulated in Treg cells from old GIA mice, and IFN-β decreased the suppressive function of Treg cells. Our findings demonstrate that increased type I IFN signaling in old Treg cells is induced only in the arthritic environment and relates to decreased suppressive function of Treg cells, gets involved in EORA.

SPATS2L is a positive feedback regulator of the type I interferon signaling pathway and plays a vital role in lupus

Through genome-wide association studies (GWAS) and integrated expression quantitative trait locus (eQTL) analyses, numerous susceptibility genes ("eGenes", whose expressions are significantly associated with common variants) associated with systemic lupus erythematosus (SLE) have been identified. Notably, a subset of these eGenes is correlated with disease activity. However, the precise mechanisms through which these genes contribute to the initiation and progression of the disease remain to be fully elucidated. In this investigation, we initially identify SPATS2L as an SLE eGene correlated with disease activity. eSignaling and transcriptomic analyses suggest its involvement in the type I interferon (IFN) pathway. We observe a significant increase in SPATS2L expression following type I IFN stimulation, and the expression levels are dependent on both the concentration and duration of stimulation. Furthermore, through dual-luciferase reporter assays, western blot analysis, and imaging flow cytometry, we confirm that SPATS2L positively modulates the type I IFN pathway, acting as a positive feedback regulator. Notably, siRNA-mediated intervention targeting SPATS2L, an interferon-inducible gene, in peripheral blood mononuclear cells (PBMCs) from patients with SLE reverses the activation of the interferon pathway. In conclusion, our research highlights the pivotal role of SPATS2L as a positive-feedback regulatory molecule within the type I IFN pathway. Our findings suggest that SPATS2L plays a critical role in the onset and progression of SLE and may serve as a promising target for disease activity assessment and intervention strategies.

Deciphering the genetics and mechanisms of predisposition to multiple myeloma

Multiple myeloma (MM) is an incurable malignancy of plasma cells. Epidemiological studies indicate a substantial heritable component, but the underlying mechanisms remain unclear. Here, in a genome-wide association study totaling 10,906 cases and 366,221 controls, we identify 35 MM risk loci, 12 of which are novel. Through functional fine-mapping and Mendelian randomization, we uncover two causal mechanisms for inherited MM risk: longer telomeres; and elevated levels of B-cell maturation antigen (BCMA) and interleukin-5 receptor alpha (IL5RA) in plasma. The largest increase in BCMA and IL5RA levels is mediated by the risk variant rs34562254-A at TNFRSF13B. While individuals with loss-of-function variants in TNFRSF13B develop B-cell immunodeficiency, rs34562254-A exerts a gain-of-function effect, increasing MM risk through amplified B-cell responses. Our results represent an analysis of genetic MM predisposition, highlighting causal mechanisms contributing to MM development.

From genetic associations to genes: methods, applications, and challenges

Genome-wide association studies (GWASs) have identified numerous genetic loci associated with human traits and diseases. However, pinpointing the causal genes remains a challenge, which impedes the translation of GWAS findings into biological insights and medical applications. In this review, we provide an in-depth overview of the methods and technologies used for prioritizing genes from GWAS loci, including gene-based association tests, integrative analysis of GWAS and molecular quantitative trait loci (xQTL) data, linking GWAS variants to target genes through enhancer-gene connection maps, and network-based prioritization. We also outline strategies for generating context-dependent xQTL data and their applications in gene prioritization. We further highlight the potential of gene prioritization in drug repurposing. Lastly, we discuss future challenges and opportunities in this field.

Quantification of escape from X chromosome inactivation with single-cell omics data reveals heterogeneity across cell types and tissues

Several X-linked genes escape from X chromosome inactivation (XCI), while differences in escape across cell types and tissues are still poorly characterized. Here, we developed scLinaX for directly quantifying relative gene expression from the inactivated X chromosome with droplet-based single-cell RNA sequencing (scRNA-seq) data. The scLinaX and differentially expressed gene analyses with large-scale blood scRNA-seq datasets consistently identified the stronger escape in lymphocytes than in myeloid cells. An extension of scLinaX to a 10x multiome dataset (scLinaX-multi) suggested a stronger escape in lymphocytes than in myeloid cells at the chromatin-accessibility level. The scLinaX analysis of human multiple-organ scRNA-seq datasets also identified the relatively strong degree of escape from XCI in lymphoid tissues and lymphocytes. Finally, effect size comparisons of genome-wide association studies between sexes suggested the underlying impact of escape on the genotype-phenotype association. Overall, scLinaX and the quantified escape catalog identified the heterogeneity of escape across cell types and tissues.

Genome-Wide Methylation Analysis Reveals a KCNK3-Prominent Causal Cascade on Hypertension

BACKGROUND

Despite advances in understanding hypertension's genetic structure, how noncoding genetic variants influence it remains unclear. Studying their interaction with DNA methylation is crucial to deciphering this complex disease's genetic mechanisms.

METHODS

We investigated the genetic and epigenetic interplay in hypertension using whole-genome bisulfite sequencing. Methylation profiling in 918 males revealed allele-specific methylation and methylation quantitative trait loci. We engineered rs1275988T/C mutant mice using CRISPR (clustered regularly interspaced short palindromic repeats)/Cas9 (CRISPR-associated protein 9), bred them for homozygosity, and subjected them to a high-salt diet. Telemetry captured their cardiovascular metrics. Protein-DNA interactions were elucidated using DNA pull-downs, mass spectrometry, and Western blots. A wire myograph assessed vascular function, and analysis of the Kcnk3 gene methylation highlighted the mutation's role in hypertension.

RESULTS

We discovered that DNA methylation-associated genetic effects, especially in non-cytosine-phosphate-guanine (non-CpG) island and noncoding distal regulatory regions, significantly contribute to hypertension predisposition. We identified distinct methylation quantitative trait locus patterns in the hypertensive population and observed that the onset of hypertension is influenced by the transmission of genetic effects through the demethylation process. By evidence-driven prioritization and in vivo experiments, we unearthed rs1275988 in a cell type-specific enhancer as a notable hypertension causal variant, intensifying hypertension through the modulation of local DNA methylation and consequential alterations in Kcnk3 gene expression and vascular remodeling. When exposed to a high-salt diet, mice with the rs1275988C/C genotype exhibited exacerbated hypertension and significant vascular remodeling, underscored by increased aortic wall thickness. The C allele of rs1275988 was associated with elevated DNA methylation levels, driving down the expression of the Kcnk3 gene by attenuating Nr2f2 (nuclear receptor subfamily 2 group F member 2) binding at the enhancer locus.

CONCLUSIONS

Our research reveals new insights into the complex interplay between genetic variations and DNA methylation in hypertension. We underscore hypomethylation's potential in hypertension onset and identify rs1275988 as a causal variant in vascular remodeling. This work advances our understanding of hypertension's molecular mechanisms and encourages personalized health care strategies.

CCR4 predicts the efficacy of abatacept in rheumatoid arthritis patients through the estimation of Th17 and Treg cell abundance

OBJECTIVES

Predicting the efficacy of biological disease-modifying antirheumatic drugs is challenging. In this study, we aimed to explore markers that predict the efficacy of abatacept in rheumatoid arthritis (RA) patients.

METHODS

Thirty RA patients receiving abatacept were recruited, and peripheral blood mononuclear cells from the participants were subjected to DNA microarray analysis. The expression of CC chemokine receptor 4 (CCR4), which was selected by the result of DNA microarray, was determined by flow cytometry in 16 newly diagnosed treatment-naïve RA patients. CCR4 expression on each helper T-cell subset was also measured.

RESULTS

CCR4 was upregulated in the abatacept responder. The expression levels of CCR4 were significantly correlated with the improvement of the Clinical Disease Activity Index. CCR4 expression was predominantly observed in CD4+ T cells in peripheral blood mononuclear cells. The percentage of CCR4-expressing CD4+ T cells was significantly higher in RA patients than in healthy individuals. Interestingly, Th17 and Treg cells expressed high levels of CCR4 compared to non-Th17-related helper T cells.

CONCLUSIONS

CCR4 is a Th17- and Treg-related gene, and the high CCR4 expression in peripheral blood samples may predict the efficacy of abatacept in RA.

An atlas of transcribed enhancers across helper T cell diversity for decoding human diseases

Transcribed enhancer maps can reveal nuclear interactions underpinning each cell type and connect specific cell types to diseases. Using a 5' single-cell RNA sequencing approach, we defined transcription start sites of enhancer RNAs and other classes of coding and noncoding RNAs in human CD4+ T cells, revealing cellular heterogeneity and differentiation trajectories. Integration of these datasets with single-cell chromatin profiles showed that active enhancers with bidirectional RNA transcription are highly cell type-specific and that disease heritability is strongly enriched in these enhancers. The resulting cell type-resolved multimodal atlas of bidirectionally transcribed enhancers, which we linked with promoters using fine-scale chromatin contact maps, enabled us to systematically interpret genetic variants associated with a range of immune-mediated diseases.

A unified metric of human immune health

Immunological health has been challenging to characterize but could be defined as the absence of immune pathology. While shared features of some immune diseases and the concept of immunologic resilience based on age-independent adaptation to antigenic stimulation have been developed, general metrics of immune health and its utility for assessing clinically healthy individuals remain ill defined. Here we integrated transcriptomics, serum protein, peripheral immune cell frequency and clinical data from 228 patients with 22 monogenic conditions impacting key immunological pathways together with 42 age- and sex-matched healthy controls. Despite the high penetrance of monogenic lesions, differences between individuals in diverse immune parameters tended to dominate over those attributable to disease conditions or medication use. Unsupervised or supervised machine learning independently identified a score that distinguished healthy participants from patients with monogenic diseases, thus suggesting a quantitative immune health metric (IHM). In ten independent datasets, the IHM discriminated healthy from polygenic autoimmune and inflammatory disease states, marked aging in clinically healthy individuals, tracked disease activities and treatment responses in both immunological and nonimmunological diseases, and predicted age-dependent antibody responses to immunizations with different vaccines. This discriminatory power goes beyond that of the classical inflammatory biomarkers C-reactive protein and interleukin-6. Thus, deviations from health in diverse conditions, including aging, have shared systemic immune consequences, and we provide a web platform for calculating the IHM for other datasets, which could empower precision medicine.

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.

The landscape of innate and adaptive immune cell subsets in patients with adult-onset Still's disease

OBJECTIVE

Adult-onset Still's disease (AOSD) is a systemic autoinflammatory disorder. The understanding of the changes in adaptive immune cells and the crosstalk between innate and adaptive immune systems in AOSD is limited. This study aimed to examine the peripheral immune cell composition and inflammatory protein levels in AOSD patients.

METHODS

Twenty-nine active AOSD patients were enrolled. Flow cytometry was used to analyse the cell populations in peripheral blood. Antibody chips were utilized to detect the protein expression profile in serum.

RESULTS

In active AOSD patients, there was an increase in the percentage of classical and non-classical monocytes among peripheral blood mononuclear cells. The proportion of natural killer (NK) cells decreased, with an increase in CD56dim NK1 cells and a decrease in CD56bright NK2 cells compared with healthy controls (HCs). The percentage of naïve central memory T cells was decreased, while the percentage of effector and effector memory T cells was increased among adaptive lymphocytes. The proportion of naïve B and memory B cells was decreased, while plasma cells were increased in AOSD patients, indicating activation of the adaptive immune system. Additionally, the serum levels of 40 proteins were elevated in AOSD patients, primarily involved in cytokine-cytokine receptor interaction, inflammatory response and regulation of mitogen-activated protein kinase cascade.

CONCLUSION

Our findings showed the activation of the innate and adaptive immune system in AOSD. The protein-protein interaction analysis suggested potential communication between innate and adaptive cell subsets. These findings provide new insights into the pathogenesis of the disease and the development of targeted therapies.

KLF2 enhancer variant rs4808485 increases lupus risk by modulating inflammasome machinery and cellular homoeostasis

OBJECTIVE

A recent genome-wide association study linked KLF2 as a novel Asian-specific locus for systemic lupus erythematosus (SLE) susceptibility. However, the underlying causal functional variant(s), cognate target gene(s) and genetic mechanisms associated with SLE risk are unknown.

METHODS

We used bioinformatics to prioritise likely functional variants and validated the best candidate with diverse experimental techniques, including genome editing. Gene expression was compared between healthy controls (HCs) and patients with SLE with or without lupus nephritis (LN+, LN-).

RESULTS

Through bioinformatics and expression quantitative trait locus analyses, we prioritised rs4808485 in active chromatin, predicted to modulate KLF2 expression. Luciferase reporter assays and chromatin immunoprecipitation-qPCR demonstrated differential allele-specific enhancer activity and binding of active histone marks (H3K27ac, H3K4me3 and H3K4me1), Pol II, CTCF, P300 and the transcription factor PARP1. Chromosome conformation capture-qPCR revealed long-range chromatin interactions between rs4808485 and the KLF2 promoter. These were directly validated by CRISPR-based genetic and epigenetic editing in Jurkat and lymphoblastoid cells. Deleting the rs4808485 enhancer in Jurkat (KO) cells disrupted NLRP3 inflammasome machinery by reducing KLF2 and increasing CASPASE1, IL-1β and GSDMD levels. Knockout cells also exhibited higher proliferation and cell-cycle progression than wild type. RNA-seq validated interplay between KLF2 and inflammasome machinery in HC, LN+ and LN-.

CONCLUSIONS

We demonstrate how rs4808485 modulates the inflammasome and cellular homoeostasis through regulating KLF2 expression. This establishes mechanistic connections between rs4808485 and SLE susceptibility.

Epigenetic targets of Janus kinase inhibitors are linked to genetic risks of rheumatoid arthritis

BACKGROUND

Current strategies that target cytokines (e.g., tumor necrosis factor (TNF)-α), or signaling molecules (e.g., Janus kinase (JAK)) have advanced the management for allergies and autoimmune diseases. Nevertheless, the molecular mechanism that underpins its clinical efficacy have largely remained elusive, especially in the local tissue environment. Here, we aimed to identify the genetic, epigenetic, and immunological targets of JAK inhibitors (JAKis), focusing on their effects on synovial fibroblasts (SFs), the major local effectors associated with destructive joint inflammation in rheumatoid arthritis (RA).

METHODS

SFs were activated by cytokines related to inflammation in RA, and were treated with three types of JAKis or a TNF-α inhibitor (TNFi). Dynamic changes in transcriptome and chromatin accessibility were profiled across samples to identify drug targets. Furthermore, the putative targets were validated using luciferase assays and clustered regularly interspaced short palindromic repeat (CRISPR)-based genome editing.

RESULTS

We found that both JAKis and the TNFi targeted the inflammatory module including IL6. Conversely, specific gene signatures that were preferentially inhibited by either of the drug classes were identified. Strikingly, RA risk enhancers for CD40 and TRAF1 were distinctively regulated by JAKis and the TNFi. We performed luciferase assays and CRISPR-based genome editing, and successfully fine-mapped the single causal variants in these loci, rs6074022-CD40 and rs7021049-TRAF1.

CONCLUSIONS

JAKis and the TNFi had a direct impact on different RA risk enhancers, and we identified nucleotide-resolution targets for both drugs. Distinctive targets of clinically effective drugs could be useful for tailoring the application of these drugs and future design of more efficient treatment strategies.

Combined Genetic Association and Differed Expression Analysis of UBE2L3 Uncovers a Genetic Regulatory Role of (Immuno)proteasome in IgA Nephropathy

Introduction

IgA nephropathy (IgAN) is a leading cause of end-stage renal disease. The exact pathogenesis of IgAN is not well defined, but some genetic studies have led to a novel discovery that the (immuno)proteasome probably plays an important role in IgAN.

Methods

We firstly analyzed the association of variants in the UBE2L3 region with susceptibility to IgAN in 3,495 patients and 9,101 controls, and then analyzed the association between lead variant and clinical phenotypes in 1,803 patients with regular follow-up data. The blood mRNA levels of members of the ubiquitin-proteasome system including UBE2L3 were analyzed in peripheral blood mononuclear cells from 53 patients and 28 healthy controls. The associations between UBE2L3 and the expression levels of genes involved in Gd-IgA1 production were also explored.

Results

The rs131654 showed the most significant association signal in UBE2L3 region (OR: 1.10, 95% CI: 1.04-1.16, p = 2.29 × 10-3), whose genotypes were also associated with the levels of Gd-IgA1 (p = 0.04). The rs131654 was observed to exert cis-eQTL effects on UBE2L3 in various tissues and cell types, particularly in immune cell types in multiple databases. The UBE2L3, LUBAC, and proteasome subunits were highly expressed in patients compared with healthy controls. High expression levels of UBE2L3 were not only associated with higher proteinuria (r = 0.34, p = 0.01) and lower eGFR (r = -0.28, p = 0.04), but also positively correlated with the gene expression of LUBAC and other proteasome subunits. Additionally, mRNA expression levels of UBE2L3 were also positively correlated with IL-6 and RELA, but negatively correlated with the expression levels of the key enzyme in the process of glycosylation including C1GALT1 and C1GALT1C1.

Conclusion

In conclusion, by combined genetic association and differed expression analysis of UBE2L3, our data support a role of genetically conferred dysregulation of the (immuno)proteasome in regulating galactose-deficient IgA1 in the development of IgAN.

Glycoprotein Acetyls Is a Novel Biomarker Predicting Cardiovascular Complications in Rheumatoid Arthritis

The relationship between rheumatoid arthritis (RA) and early onset atherosclerosis is well depicted, each with an important inflammatory component. Glycoprotein acetyls (GlycA), a novel biomarker of inflammation, may play a role in the manifestation of these two inflammatory conditions. The present study examined a potential mediating role of GlycA within the RA-atherosclerosis relationship to determine whether it accounts for the excess risk of cardiovascular disease over that posed by lipid risk factors. The UK Biobank dataset was acquired to establish associations among RA, atherosclerosis, GlycA, and major lipid factors: total cholesterol (TC), high- and low-density lipoprotein (HDL, LDL) cholesterol, and triglycerides (TGs). Genome-wide association study summary statistics were collected from various resources to perform genetic analyses. Causality among variables was tested using Mendelian Randomization (MR) analysis. Genes of interest were identified using colocalization analysis and gene enrichment analysis. MR results appeared to indicate that the genetic relationship between GlycA and RA and also between RA and atherosclerosis was explained by horizontal pleiotropy (p-value = 0.001 and <0.001, respectively), while GlycA may causally predict atherosclerosis (p-value = 0.017). Colocalization analysis revealed several functionally relevant genes shared between GlycA and all the variables assessed. Two loci were apparent in all relationships tested and included the HLA region as well as SLC22A1. GlycA appears to mediate the RA-atherosclerosis relationship through several possible pathways. GlycA, although pleiotropically related to RA, appears to causally predict atherosclerosis. Thus, GlycA is suggested as a significant factor in the etiology of atherosclerosis development in RA.

Long-read sequencing for 29 immune cell subsets reveals disease-linked isoforms

Alternative splicing events are a major causal mechanism for complex traits, but they have been understudied due to the limitation of short-read sequencing. Here, we generate a full-length isoform annotation of human immune cells from an individual by long-read sequencing for 29 cell subsets. This contains a number of unannotated transcripts and isoforms such as a read-through transcript of TOMM40-APOE in the Alzheimer's disease locus. We profile characteristics of isoforms and show that repetitive elements significantly explain the diversity of unannotated isoforms, providing insight into the human genome evolution. In addition, some of the isoforms are expressed in a cell-type specific manner, whose alternative 3'-UTRs usage contributes to their specificity. Further, we identify disease-associated isoforms by isoform switch analysis and by integration of several quantitative trait loci analyses with genome-wide association study data. Our findings will promote the elucidation of the mechanism of complex diseases via alternative splicing.

Fundamentals for predicting transcriptional regulations from DNA sequence patterns

Cell-type-specific regulatory elements, cataloged through extensive experiments and bioinformatics in large-scale consortiums, have enabled enrichment analyses of genetic associations that primarily utilize positional information of the regulatory elements. These analyses have identified cell types and pathways genetically associated with human complex traits. However, our understanding of detailed allelic effects on these elements' activities and on-off states remains incomplete, hampering the interpretation of human genetic study results. This review introduces machine learning methods to learn sequence-dependent transcriptional regulation mechanisms from DNA sequences for predicting such allelic effects (not associations). We provide a concise history of machine-learning-based approaches, the requirements, and the key computational processes, focusing on primers in machine learning. Convolution and self-attention, pivotal in modern deep-learning models, are explained through geometrical interpretations using dot products. This facilitates understanding of the concept and why these have been used for machine learning for DNA sequences. These will inspire further research in this genetics and genomics field.

DISSECT: deep semi-supervised consistency regularization for accurate cell type fraction and gene expression estimation

Cell deconvolution is the estimation of cell type fractions and cell type-specific gene expression from mixed data. An unmet challenge in cell deconvolution is the scarcity of realistic training data and the domain shift often observed in synthetic training data. Here, we show that two novel deep neural networks with simultaneous consistency regularization of the target and training domains significantly improve deconvolution performance. Our algorithm, DISSECT, outperforms competing algorithms in cell fraction and gene expression estimation by up to 14 percentage points. DISSECT can be easily adapted to other biomedical data types, as exemplified by our proteomic deconvolution experiments.

Loop Catalog: a comprehensive HiChIP database of human and mouse samples

HiChIP enables cost-effective and high-resolution profiling of regulatory and structural loops. To leverage the increasing number of publicly available HiChIP datasets from diverse cell lines and primary cells, we developed the Loop Catalog (https://loopcatalog.lji.org), a web-based database featuring HiChIP loop calls for 1319 samples across 133 studies and 44 high-resolution Hi-C loop calls. We demonstrate its utility in interpreting fine-mapped GWAS variants (SNP-to-gene linking), in identifying enriched sequence motifs and motif pairs at loop anchors, and in network-level analysis of loops connecting regulatory elements (community detection). Our comprehensive catalog, spanning over 4M unique 5kb loops, along with the accompanying analysis modalities constitutes an important resource for studies in gene regulation and genome organization.

A genome-wide association study identified PTPN2 as a population-specific susceptibility gene locus for primary biliary cholangitis

BACKGROUND AND AIMS

Previous genome-wide association studies (GWAS) have indicated the involvement of shared (population-nonspecific) and nonshared (population-specific) susceptibility genes in the pathogenesis of primary biliary cholangitis (PBC) among European and East-Asian populations. Although a meta-analysis of these distinct populations has recently identified more than 20 novel PBC susceptibility loci, analyses of population-specific genetic architecture are still needed for a more comprehensive search for genetic factors in PBC.

APPROACH AND RESULTS

Protein tyrosine phosphatase nonreceptor type 2 ( PTPN2) was identified as a novel PBC susceptibility gene locus through GWAS and subsequent genome-wide meta-analysis involving 2181 cases and 2699 controls from the Japanese population (GWAS-lead variant: rs8098858, p = 2.6 × 10 -8 ). In silico and in vitro functional analyses indicated that the risk allele of rs2292758, which is a primary functional variant, decreases PTPN2 expression by disrupting Sp1 binding to the PTPN2 promoter in T follicular helper cells and plasmacytoid dendritic cells. Infiltration of PTPN2-positive T-cells and plasmacytoid dendritic cells was confirmed in the portal area of the PBC liver by immunohistochemistry. Furthermore, transcriptomic analysis of PBC-liver samples indicated the presence of a compromised negative feedback loop in vivo between PTPN2 and IFNG in patients carrying the risk allele of rs2292758.

CONCLUSIONS

PTPN2 , a novel susceptibility gene for PBC in the Japanese population, may be involved in the pathogenesis of PBC through an insufficient negative feedback loop caused by the risk allele of rs2292758 in IFN-γ signaling. This suggests that PTPN2 could be a potential molecular target for PBC treatment.

Genetic associations of cardiovascular risk genes in European patients with coronary artery spasm

BACKGROUND

Coronary artery spasm (CAS) is a frequent finding in patients presenting with angina pectoris. Although the pathogenesis of CAS is incompletely understood, previous studies suggested a genetic contribution. Our study aimed to elucidate genetic variants in a cohort of European patients with angina and unobstructed coronary arteries who underwent acetylcholine (ACh) provocation testing.

METHODS

A candidate association analysis of 208 genes previously associated with cardiovascular conditions was performed using genotyped and imputed variants in patients grouped in epicardial (focal, diffuse) CAS (n = 119) and microvascular CAS (n = 87). Patients with a negative ACh test result (n = 45) served as controls.

RESULTS

We found no association below the genome-wide significance threshold of p < 5 × 10-8, thus not confirming variants in ALDH2, NOS3, and ROCK2 previously reported in CAS patients of Asian ancestry. However, the analysis identified suggestive associations (p < 10-05) for the groups of focal epicardial CAS (CDH13) and diffuse epicardial CAS (HDAC9, EDN1). Downstream analysis of the potential EDN1 risk locus showed that CAS patients have significantly increased plasma endothelin-1 levels (ET-1) compared to controls. An EDN1 haplotype comprising rs9349379 and rs2070698 was significantly associated to ET-1 levels (p = 0.01).

CONCLUSIONS

In summary, we suggest EDN1 as potential genetic risk loci for patients with diffuse epicardial CAS, and European ancestry. Plasma ET-1 levels may serve as a potential cardiac marker.

Novel potential lncRNA biomarker in B cells indicates essential pathogenic pathway activation in patients with SLE

OBJECTIVES

Systemic lupus erythematosus (SLE) is a highly heterogeneous disease, and B cell abnormalities play a central role in the pathogenesis of SLE. Long non-coding RNAs (lncRNAs) have also been implicated in the pathogenesis of SLE. The expression of lncRNAs is finely regulated and cell-type dependent, so we aimed to identify B cell-expressing lncRNAs as biomarkers for SLE, and to explore their ability to reflect the status of SLE critical pathway and disease activity.

METHODS

Weighted gene coexpression network analysis (WGCNA) was used to cluster B cell-expressing genes of patients with SLE into different gene modules and relate them to clinical features. Based on the results of WGCNA, candidate lncRNA levels were further explored in public bulk and single-cell RNA-sequencing data. In another independent cohort, the levels of the candidate were detected by RT-qPCR and the correlation with disease activity was analysed.

RESULTS

WGCNA analysis revealed one gene module significantly correlated with clinical features, which was enriched in type I interferon (IFN) pathway. Among non-coding genes in this module, lncRNA RP11-273G15.2 was differentially expressed in all five subsets of B cells from patients with SLE compared with healthy controls and other autoimmune diseases. RT-qPCR validated that RP11-273G15.2 was highly expressed in SLE B cells and positively correlated with IFN scores (r=0.7329, p<0.0001) and disease activity (r=0.4710, p=0.0005).

CONCLUSION

RP11-273G15.2 could act as a diagnostic and disease activity monitoring biomarker for SLE, which might have the potential to guide clinical management.

A noncoding regulatory variant in IKZF1 increases acute lymphoblastic leukemia risk in Hispanic/Latino children

Hispanic/Latino children have the highest risk of acute lymphoblastic leukemia (ALL) in the US compared to other racial/ethnic groups, yet the basis of this remains incompletely understood. Through genetic fine-mapping analyses, we identified a new independent childhood ALL risk signal near IKZF1 in self-reported Hispanic/Latino individuals, but not in non-Hispanic White individuals, with an effect size of ∼1.44 (95% confidence interval = 1.33-1.55) and a risk allele frequency of ∼18% in Hispanic/Latino populations and <0.5% in European populations. This risk allele was positively associated with Indigenous American ancestry, showed evidence of selection in human history, and was associated with reduced IKZF1 expression. We identified a putative causal variant in a downstream enhancer that is most active in pro-B cells and interacts with the IKZF1 promoter. This variant disrupts IKZF1 autoregulation at this enhancer and results in reduced enhancer activity in B cell progenitors. Our study reveals a genetic basis for the increased ALL risk in Hispanic/Latino children.

Age-associated CD4+ T cells with B cell-promoting functions are regulated by ZEB2 in autoimmunity

Aging is a significant risk factor for autoimmunity, and many autoimmune diseases tend to onset during adulthood. We conducted an extensive analysis of CD4+ T cell subsets from 354 patients with autoimmune disease and healthy controls via flow cytometry and bulk RNA sequencing. As a result, we identified a distinct CXCR3midCD4+ effector memory T cell subset that expands with age, which we designated "age-associated T helper (THA) cells." THA cells exhibited both a cytotoxic phenotype and B cell helper functions, and these features were regulated by the transcription factor ZEB2. Consistent with the highly skewed T cell receptor usage of THA cells, gene expression in THA cells from patients with systemic lupus erythematosus reflected disease activity and was affected by treatment with a calcineurin inhibitor. Moreover, analysis of single-cell RNA sequencing data revealed that THA cells infiltrate damaged organs in patients with autoimmune diseases. Together, our characterization of THA cells may facilitate improved understanding of the relationship between aging and autoimmune diseases.

Genetic regulatory effects in response to a high-cholesterol, high-fat diet in baboons

Steady-state expression quantitative trait loci (eQTLs) explain only a fraction of disease-associated loci identified through genome-wide association studies (GWASs), while eQTLs involved in gene-by-environment (GxE) interactions have rarely been characterized in humans due to experimental challenges. Using a baboon model, we found hundreds of eQTLs that emerge in adipose, liver, and muscle after prolonged exposure to high dietary fat and cholesterol. Diet-responsive eQTLs exhibit genomic localization and genic features that are distinct from steady-state eQTLs. Furthermore, the human orthologs associated with diet-responsive eQTLs are enriched for GWAS genes associated with human metabolic traits, suggesting that context-responsive eQTLs with more complex regulatory effects are likely to explain GWAS hits that do not seem to overlap with standard eQTLs. Our results highlight the complexity of genetic regulatory effects and the potential of eQTLs with disease-relevant GxE interactions in enhancing the understanding of GWAS signals for human complex disease using non-human primate models.

SingleQ: a comprehensive database of single-cell expression quantitative trait loci (sc-eQTLs) cross human tissues

Mapping of expression quantitative trait loci (eQTLs) and other molecular QTLs can help characterize the modes of action of disease-associated genetic variants. However, current eQTL databases present data from bulk RNA-seq approaches, which cannot shed light on the cell type- and environment-specific regulation of disease-associated genetic variants. Here, we introduce our Single-cell eQTL Interactive Database which collects single-cell eQTL (sc-eQTL) datasets and provides online visualization of sc-eQTLs across different cell types in a user-friendly manner. Although sc-eQTL mapping is still in its early stage, our database curates the most comprehensive summary statistics of sc-eQTLs published to date. sc-eQTL studies have revolutionized our understanding of gene regulation in specific cellular contexts, and we anticipate that our database will further accelerate the research of functional genomics. Database URL: http://www.sqraolab.com/scqtl.

Pathophysiology and stratification of treatment-resistant rheumatoid arthritis

Early diagnosis and timely therapeutic intervention are clinical challenges of rheumatoid arthritis (RA), especially for treatment-resistant or difficult-to-treat patients. Little is known about the immunological mechanisms involved in refractory RA. In this review, we summarize previous research findings on the immunological mechanisms of treatment-resistant RA. Genetic prediction of treatment-resistant RA is challenging. Patients with and without anti-cyclic citrullinated peptide autoantibodies are considered part of distinct subgroups, especially regarding long-term clinical prognosis and treatment responses. B cells, T cells and other immune cells and fibroblasts are of pathophysiological importance and are associated with treatment responses. Finally, we propose a new hypothesis that stratifies patients with RA into two subgroups with distinct immunological pathologies based on our recent immunomics analysis of RA. One RA subgroup with a favorable prognosis is characterized by increased interferon signaling. Another subgroup with a worse prognosis is characterized by enhanced acquired immune responses. Increases in dendritic cell precursors and diversified autoreactive anti-modified protein antibodies may have pathophysiological roles, especially in the latter subgroup. These findings that improve treatment response predictions might contribute to future precision medicine for RA.

Three-Dimensional Chromosomal Landscape Revealing miR-146a Dysfunctional Enhancer in Lupus and Establishing a CRISPR-Mediated Approach to Inhibit the Interferon Pathway

OBJECTIVE

The diminished expression of microRNA-146a (miR-146a) in systemic lupus erythematosus (SLE) contributes to the aberrant activation of the interferon pathway. Despite its significance, the underlying mechanism driving this reduced expression remains elusive. Considering the integral role of enhancers in steering gene expression, our study seeks to pinpoint the SLE-affected enhancers responsible for modulating miR-146a expression. Additionally, we aim to elucidate the mechanisms by which these enhancers influence the contribution of miR-146a to the activation of the interferon pathway.

METHODS

Circular chromosome conformation capture sequencing and epigenomic profiles were applied to identify candidate enhancers of miR-146a. CRISPR activation was performed to screen functional enhancers. Differential analysis of chromatin accessibility was used to identify SLE-dysregulated enhancers, and the mechanism underlying enhancer dysfunction was investigated by analyzing transcription factor binding. The therapeutic value of a lupus-related enhancer was further evaluated by targeting it in the peripheral blood mononuclear cells (PBMCs) of patients with SLE through a CRISPR activation approach.

RESULTS

We identified shared and cell-specific enhancers of miR-146a in distinct immune cells. An enhancer 32.5 kb downstream of miR-146a possesses less accessibility in SLE, and its chromatin openness was negatively correlated with SLE disease activity. Moreover, CCAAT/enhancer binding protein α, a down-regulated transcription factor in patients with SLE, binds to the 32.5-kb enhancer and induces the epigenomic change of this locus. Furthermore, CRISPR-based activation of this enhancer in SLE PBMCs could inhibit the activity of interferon pathway.

CONCLUSION

Our work defines a promising target for SLE intervention. We adopted integrative approaches to define cell-specific and functional enhancers of the SLE critical gene and investigated the mechanism underlying its dysregulation mediated by a lupus-related enhancer.

Non-coding DNA variants for risk in lupus

Systemic Lupus Erythematosus (SLE) is a multifactorial autoimmune disease that arises from a dynamic interplay between genetics and environmental triggers. The advent of sophisticated genomics technology has catalyzed a shift in our understanding of disease etiology, spotlighting the pivotal role of non-coding DNA variants in SLE pathogenesis. In this review, we present a comprehensive examination of the non-coding variants associated with SLE, shedding light on their role in influencing disease risk and progression. We discuss the latest methodological advancements that have been instrumental in the identification and functional characterization of these genomic elements, with a special focus on the transformative power of CRISPR-based gene-editing technologies. Additionally, the review probes into the therapeutic opportunities that arise from modulating non-coding regions associated with SLE. Through an exploration of the complex network of non-coding DNA, this review aspires to decode the genetic puzzle of SLE and set the stage for groundbreaking gene-based therapeutic interventions and the advancement of precision medicine strategies tailored to SLE management.

Single-cell transcriptome landscape of circulating CD4+ T cell populations in autoimmune diseases

CD4+ T cells are key mediators of various autoimmune diseases; however, their role in disease progression remains unclear due to cellular heterogeneity. Here, we evaluated CD4+ T cell subpopulations using decomposition-based transcriptome characterization and canonical clustering strategies. This approach identified 12 independent gene programs governing whole CD4+ T cell heterogeneity, which can explain the ambiguity of canonical clustering. In addition, we performed a meta-analysis using public single-cell datasets of over 1.8 million peripheral CD4+ T cells from 953 individuals by projecting cells onto the reference and cataloging cell frequency and qualitative alterations of the populations in 20 diseases. The analyses revealed that the 12 transcriptional programs were useful in characterizing each autoimmune disease and predicting its clinical status. Moreover, genetic variants associated with autoimmune diseases showed disease-specific enrichment within the 12 gene programs. The results collectively provide a landscape of single-cell transcriptomes of CD4+ T cell subpopulations involved in autoimmune disease.

sscNOVA: a semi-supervised convolutional neural network for predicting functional regulatory variants in autoimmune diseases

Genome-wide association studies (GWAS) have identified thousands of variants in the human genome with autoimmune diseases. However, identifying functional regulatory variants associated with autoimmune diseases remains challenging, largely because of insufficient experimental validation data. We adopt the concept of semi-supervised learning by combining labeled and unlabeled data to develop a deep learning-based algorithm framework, sscNOVA, to predict functional regulatory variants in autoimmune diseases and analyze the functional characteristics of these regulatory variants. Compared to traditional supervised learning methods, our approach leverages more variants' data to explore the relationship between functional regulatory variants and autoimmune diseases. Based on the experimentally curated testing dataset and evaluation metrics, we find that sscNOVA outperforms other state-of-the-art methods. Furthermore, we illustrate that sscNOVA can help to improve the prioritization of functional regulatory variants from lead single-nucleotide polymorphisms and the proxy variants in autoimmune GWAS data.

Mapping the functional impact of non-coding regulatory elements in primary T cells through single-cell CRISPR screens

BACKGROUND

Drug targets with genetic evidence are expected to increase clinical success by at least twofold. Yet, translating disease-associated genetic variants into functional knowledge remains a fundamental challenge of drug discovery. A key issue is that the vast majority of complex disease associations cannot be cleanly mapped to a gene. Immune disease-associated variants are enriched within regulatory elements found in T-cell-specific open chromatin regions.

RESULTS

To identify genes and molecular programs modulated by these regulatory elements, we develop a CRISPRi-based single-cell functional screening approach in primary human T cells. Our pipeline enables the interrogation of transcriptomic changes induced by the perturbation of regulatory elements at scale. We first optimize an efficient CRISPRi protocol in primary CD4+ T cells via CROPseq vectors. Subsequently, we perform a screen targeting 45 non-coding regulatory elements and 35 transcription start sites and profile approximately 250,000 T -cell single-cell transcriptomes. We develop a bespoke analytical pipeline for element-to-gene (E2G) mapping and demonstrate that our method can identify both previously annotated and novel E2G links. Lastly, we integrate genetic association data for immune-related traits and demonstrate how our platform can aid in the identification of effector genes for GWAS loci.

CONCLUSIONS

We describe "primary T cell crisprQTL" - a scalable, single-cell functional genomics approach for mapping regulatory elements to genes in primary human T cells. We show how this framework can facilitate the interrogation of immune disease GWAS hits and propose that the combination of experimental and QTL-based techniques is likely to address the variant-to-function problem.

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.

Regulus infers signed regulatory relations from few samples' information using discretization and likelihood constraints

MOTIVATION

Transcriptional regulation is performed by transcription factors (TF) binding to DNA in context-dependent regulatory regions and determines the activation or inhibition of gene expression. Current methods of transcriptional regulatory circuits inference, based on one or all of TF, regions and genes activity measurements require a large number of samples for ranking the candidate TF-gene regulation relations and rarely predict whether they are activations or inhibitions. We hypothesize that transcriptional regulatory circuits can be inferred from fewer samples by (1) fully integrating information on TF binding, gene expression and regulatory regions accessibility, (2) reducing data complexity and (3) using biology-based likelihood constraints to determine the global consistency between a candidate TF-gene relation and patterns of genes expressions and region activations, as well as qualify regulations as activations or inhibitions.

RESULTS

We introduce Regulus, a method which computes TF-gene relations from gene expressions, regulatory region activities and TF binding sites data, together with the genomic locations of all entities. After aggregating gene expressions and region activities into patterns, data are integrated into a RDF (Resource Description Framework) endpoint. A dedicated SPARQL (SPARQL Protocol and RDF Query Language) query retrieves all potential relations between expressed TF and genes involving active regulatory regions. These TF-region-gene relations are then filtered using biological likelihood constraints allowing to qualify them as activation or inhibition. Regulus provides signed relations consistent with public databases and, when applied to biological data, identifies both known and potential new regulators. Regulus is devoted to context-specific transcriptional circuits inference in human settings where samples are scarce and cell populations are closely related, using discretization into patterns and likelihood reasoning to decipher the most robust regulatory relations.

rs10924104 in the expression enhancer motif of CD58 confers susceptibility to human autoimmune diseases

CD58 plays roles in cell adhesion and co-stimulation with antigen presentation from major histocompatibility complex class II on antigen-presenting cells to T-cell antigen receptors on naïve T cells. CD58 reportedly contributes to the development of various human autoimmune diseases. Recently, genome-wide association studies (GWASs) identified CD58 as a susceptibility locus for autoimmune diseases such as systemic lupus erythematosus (SLE), multiple sclerosis (MS), and primary biliary cholangitis (PBC). However, the primary functional variant and molecular mechanisms of susceptibility to autoimmune diseases in the CD58 locus were not clarified. Here, rs10924104, located in the ZNF35-binding motif within the gene expression regulatory motif, was identified as the primary functional variant for SLE, MS, and PBC among genetic variants showing stronger linkage disequilibrium (LD) with GWAS-lead variants in the CD58 locus. Expression-quantitative trait locus (e-QTL) data for each distinct blood cell type and in vitro functional analysis using the CRISPR/Cas9 system corroborated the functional role of rs10924104 in the upregulation of CD58 transcription by the disease-risk allele. Additionally, the strength of disease susceptibility observed in the CD58 locus could be accounted for by the strength of LD between rs10924104 and each GWAS-lead variant. In conclusion, the present study demonstrated for the first time the existence of a shared autoimmune disease-related primary functional variant (i.e., rs10924104) that regulates the expression of CD58. Clarifying the molecular mechanism of disease susceptibility derived from such a shared genetic background is important for understanding human autoimmune diseases and human immunology.

Dissecting the Immune System through Gene Regulation

The immune system plays a dual role in human health, functioning both as a protector against pathogens and, at times, as a contributor to disease. This feature emphasizes the importance to uncover the underlying causes of its malfunctions, necessitating an in-depth analysis in both pathological and physiological conditions to better understand the immune system and immune disorders. Recent advances in scientific technology have enabled extensive investigations into gene regulation, a crucial mechanism governing cellular functionality. Studying gene regulatory mechanisms within the immune system is a promising avenue for enhancing our understanding of immune cells and the immune system as a whole. The gene regulatory mechanisms, revealed through various methodologies, and their implications in the field of immunology are discussed in this chapter.

Single-cell eQTL mapping in yeast reveals a tradeoff between growth and reproduction

Expression quantitative trait loci (eQTLs) provide a key bridge between noncoding DNA sequence variants and organismal traits. The effects of eQTLs can differ among tissues, cell types, and cellular states, but these differences are obscured by gene expression measurements in bulk populations. We developed a one-pot approach to map eQTLs in Saccharomyces cerevisiae by single-cell RNA sequencing (scRNA-seq) and applied it to over 100,000 single cells from three crosses. We used scRNA-seq data to genotype each cell, measure gene expression, and classify the cells by cell-cycle stage. We mapped thousands of local and distant eQTLs and identified interactions between eQTL effects and cell-cycle stages. We took advantage of single-cell expression information to identify hundreds of genes with allele-specific effects on expression noise. We used cell-cycle stage classification to map 20 loci that influence cell-cycle progression. One of these loci influenced the expression of genes involved in the mating response. We showed that the effects of this locus arise from a common variant (W82R) in the gene GPA1, which encodes a signaling protein that negatively regulates the mating pathway. The 82R allele increases mating efficiency at the cost of slower cell-cycle progression and is associated with a higher rate of outcrossing in nature. Our results provide a more granular picture of the effects of genetic variants on gene expression and downstream traits.

ChIP-DIP: A multiplexed method for mapping hundreds of proteins to DNA uncovers diverse regulatory elements controlling gene expression

Gene expression is controlled by the dynamic localization of thousands of distinct regulatory proteins to precise regions of DNA. Understanding this cell-type specific process has been a goal of molecular biology for decades yet remains challenging because most current DNA-protein mapping methods study one protein at a time. To overcome this, we developed ChIP-DIP (ChIP Done In Parallel), a split-pool based method that enables simultaneous, genome-wide mapping of hundreds of diverse regulatory proteins in a single experiment. We demonstrate that ChIP-DIP generates highly accurate maps for all classes of DNA-associated proteins, including histone modifications, chromatin regulators, transcription factors, and RNA Polymerases. Using these data, we explore quantitative combinations of protein localization on genomic DNA to define distinct classes of regulatory elements and their functional activity. Our data demonstrate that ChIP-DIP enables the generation of 'consortium level', context-specific protein localization maps within any molecular biology lab.

COX5A as a potential biomarker for disease activity and organ damage in lupus

Systemic Lupus Erythematosus (SLE) is a complex autoimmune disease with limited therapeutic targets or clinical outcome predictors. This study aimed to gain more insights into the underlying immunological pathways and prognostic biomarkers of SLE. Integrated analyses of RNA-seq data from 64 SLE and 62 healthy controls, examining 27 immune cell types to explore the key pathways and driver genes in SLE pathogenesis. Single-cell RNA sequencing data from the skin and kidney were used to determine the association of COX5A expression with organ damage. The associations of COX5A with SLE phenotypes were further evaluated in two independent cohorts, and receiver operating characteristic (ROC) curves were constructed to assess the value of COX5A as a biomarker for disease activity and organ damage in SLE. We found that oxidative phosphorylation (OXPHOS) is the most significantly altered metabolic pathway in SLE, especially in effector T cells. Notably, we identified an OXPHOS-related enzyme, COX5A, whose expression was significantly higher in effector T cells than in naïve T cells and showed associations with disease activity, organ damage, and steroid treatment of SLE. Furthermore, ROC curves showed that COX5A is a robust biomarker for disease activity, kidney involvement, and new-onset skin lesions, with the area under the curve (AUC) values of 0.880, 0.801, and 0.805, respectively. Our results identified the OXPHOS signature as a prominent feature in SLE T cells, and COX5A as a potential candidate biomarker for disease activity and organ damage in SLE.

eQTL studies: from bulk tissues to single cells

An expression quantitative trait locus (eQTL) is a chromosomal region where genetic variants are associated with the expression levels of specific genes that can be both nearby or distant. The identifications of eQTLs for different tissues, cell types, and contexts have led to a better understanding of the dynamic regulations of gene expressions and implications of functional genes and variants for complex traits and diseases. Although most eQTL studies have been performed on data collected from bulk tissues, recent studies have demonstrated the importance of cell-type-specific and context-dependent gene regulations in biological processes and disease mechanisms. In this review, we discuss statistical methods that have been developed to enable the detection of cell-type-specific and context-dependent eQTLs from bulk tissues, purified cell types, and single cells. We also discuss the limitations of the current methods and future research opportunities.

Mapping the dynamic genetic regulatory architecture of HLA genes at single-cell resolution

The human leukocyte antigen (HLA) locus plays a critical role in complex traits spanning autoimmune and infectious diseases, transplantation and cancer. While coding variation in HLA genes has been extensively documented, regulatory genetic variation modulating HLA expression levels has not been comprehensively investigated. Here we mapped expression quantitative trait loci (eQTLs) for classical HLA genes across 1,073 individuals and 1,131,414 single cells from three tissues. To mitigate technical confounding, we developed scHLApers, a pipeline to accurately quantify single-cell HLA expression using personalized reference genomes. We identified cell-type-specific cis-eQTLs for every classical HLA gene. Modeling eQTLs at single-cell resolution revealed that many eQTL effects are dynamic across cell states even within a cell type. HLA-DQ genes exhibit particularly cell-state-dependent effects within myeloid, B and T cells. For example, a T cell HLA-DQA1 eQTL ( rs3104371 ) is strongest in cytotoxic cells. Dynamic HLA regulation may underlie important interindividual variability in immune responses.

A Non-Coding Variant in SLC15A4 Modulates Enhancer Activity and Lysosomal Deacidification Linked to Lupus Susceptibility

Systemic lupus erythematosus (SLE) is a complex autoimmune disease with a strong genetic basis. Despite the identification of several single nucleotide polymorphisms (SNPs) near the SLC15A4 gene that are significantly associated with SLE across multiple populations, specific causal SNP(s) and molecular mechanisms responsible for disease susceptibility are unknown. To address this gap, we employed bioinformatics, expression quantitative trait loci (eQTLs), and 3D chromatin interaction analysis to nominate a likely functional variant, rs35907548, in an active intronic enhancer of SLC15A4. Through luciferase reporter assays followed by chromatin immunoprecipitation (ChIP)-qPCR, we observed significant allele-specific enhancer effects of rs35907548 in diverse cell lines. The rs35907548 risk allele T is associated with increased regulatory activity and target gene expression, as shown by eQTLs and chromosome conformation capture (3C)-qPCR. The latter revealed long-range chromatin interactions between the rs35907548 enhancer and the promoters of SLC15A4, GLTLD1, and an uncharacterized lncRNA. The enhancer-promoter interactions and expression effects were validated by CRISPR/Cas9 knock-out (KO) of the locus in HL60 promyeloblast cells. KO cells also displayed dramatically dysregulated endolysosomal pH regulation. Together, our data show that the rs35907548 risk allele affects multiple aspects of cellular physiology and may directly contribute to SLE.

Multimodal repertoire analysis unveils B cell biology in immune-mediated diseases

OBJECTIVES

Despite the involvement of B cells in the pathogenesis of immune-mediated diseases (IMDs), biological mechanisms underlying their function are scarcely understood. To overcome this gap, here we constructed and investigated a large-scale repertoire catalogue of five B cell subsets of patients with IMDs.

METHODS

We mapped B cell receptor regions from RNA sequencing data of sorted B cell subsets. Our dataset consisted of 595 donors under IMDs and health. We characterised the repertoire features from various aspects, including their association with immune cell transcriptomes and clinical features and their response to belimumab treatment.

RESULTS

Heavy-chain complementarity-determining region 3 (CDR-H3) length among naïve B cells was shortened among autoimmune diseases. Strong negative correlation between interferon signature strength and CDR-H3 length was observed in naïve B cells and suggested the role for interferon in premature B cell development. VDJ gene usage was skewed especially in plasmablasts and unswitched-memory B cells of patients with systemic lupus erythematosus (SLE). We developed a scoring system to quantify this skewing, and it positively correlated with peripheral helper T cell transcriptomic signatures and negatively correlated with the amount of somatic hyper mutations in plasmablasts, suggesting the association of extrafollicular pathway. Further, this skewing led to high usage of IGHV4-34 gene with 9G4 idiotypes in unswitched-memory B cells, which showed a prominent positive correlation with disease activity in SLE. Gene usage skewing in unswitched-memory B cells was ameliorated after belimumab treatment.

CONCLUSIONS

Our multimodal repertoire analysis enabled us the system-level understanding of B cell abnormality in diseases.

Leveraging molecular quantitative trait loci to comprehend complex diseases/traits from the omics perspective

Comprehending the molecular basis of quantitative genetic variation is a principal goal for complex diseases or traits. Molecular quantitative trait loci (molQTLs) have made it possible to investigate the effects of genetic variants hiding behind large-scale omics data. A deeper understanding of molQTL is urgently required in light of the multi-dimensionalization of omics data to more fully elucidate the pertinent biological mechanisms. Herein, we reviewed molQTLs with the corresponding resource from the omics perspective and further discussed the integrative strategy of GWAS-molQTL to infer their causal effects. Subsequently, we described the opportunities and challenges encountered by molQTL. The case studies showed that molQTL is essential for complex diseases and traits, whether single- or multi-omics QTLs. Overall, we highlighted the functional significance of genetic variants to employ the discovery of molQTL in complex diseases and traits.

Systematic differences in discovery of genetic effects on gene expression and complex traits

Most signals in genome-wide association studies (GWAS) of complex traits implicate noncoding genetic variants with putative gene regulatory effects. However, currently identified regulatory variants, notably expression quantitative trait loci (eQTLs), explain only a small fraction of GWAS signals. Here, we show that GWAS and cis-eQTL hits are systematically different: eQTLs cluster strongly near transcription start sites, whereas GWAS hits do not. Genes near GWAS hits are enriched in key functional annotations, are under strong selective constraint and have complex regulatory landscapes across different tissue/cell types, whereas genes near eQTLs are depleted of most functional annotations, show relaxed constraint, and have simpler regulatory landscapes. We describe a model to understand these observations, including how natural selection on complex traits hinders discovery of functionally relevant eQTLs. Our results imply that GWAS and eQTL studies are systematically biased toward different types of variant, and support the use of complementary functional approaches alongside the next generation of eQTL studies.

A Non-Coding Fc Gamma Receptor Cis-Regulatory Variant within the 1q23 Gene Cluster Is Associated with Plasmodium falciparum Infection in Children Residing in Burkina Faso

Antibodies play a crucial role in activating protective immunity against malaria by interacting with Fc-gamma receptors (FcγRs). Genetic variations in genes encoding FcγRs can affect immune cell responses to the parasite. In this study, our aim was to investigate whether non-coding variants that regulate FcγR expression could influence the prevalence of Plasmodium falciparum infection. Through bioinformatics approaches, we selected expression quantitative trait loci (eQTL) for FCGR2A, FCGR2B, FCGR2C, FCGR3A, and FCGR3B genes encoding FcγRs (FCGR), in whole blood. We prioritized two regulatory variants, rs2099684 and rs1771575, located in open genomic regions. These variants were identified using RegVar, ImmuNexUT, and transcription factor annotations specific to immune cells. In addition to these, we genotyped the coding variants FCGR2A/rs1801274 and FCGR2B/rs1050501 in 234 individuals from a malaria-endemic area in Burkina Faso. We conducted age and family-based analyses to evaluate associations with the prevalence of malarial infection in both children and adults. The analysis revealed that the regulatory rs1771575-CC genotype was predicted to influence FCGR2B/FCGR2C/FCGR3A transcripts in immune cells and was the sole variant associated with a higher prevalence of malarial infection in children. In conclusion, this study identifies the rs1771575 cis-regulatory variant affecting several FcγRs in myeloid and neutrophil cells and associates it with the inter-individual capacity of children living in Burkina Faso to control malarial infection.

Investigation of immune-related diseases using patient-derived induced pluripotent stem cells

The precise pathogenesis of immune-related diseases remains unclear, and new effective therapeutic choices are required for the induction of remission or cure in these diseases. Basic research utilizing immune-related disease patient-derived induced pluripotent stem (iPS) cells is expected to be a promising platform for elucidating the pathogenesis of the diseases and for drug discovery. Since autoinflammatory diseases are usually monogenic, genetic mutations affect the cell function and patient-derived iPS cells tend to exhibit disease-specific phenotypes. In particular, iPS cell-derived monocytic cells and macrophages can be used for functional experiments, such as inflammatory cytokine production, and are often employed in research on patients with autoinflammatory diseases.On the other hand, the utilization of disease-specific iPS cells is less successful for research on autoimmune diseases. One reason for this is that autoimmune diseases are usually polygenic, which makes it challenging to determine which factors cause the phenotypes of patient-derived iPS cells are caused by. Another reason is that protocols for differentiating some lymphocytes associated with autoimmunity, such as CD4+T cells or B cells, from iPS cells have not been well established. Nevertheless, several groups have reported studies utilizing autoimmune disease patient-derived iPS cells, including patients with rheumatoid arthritis, systemic lupus erythematosus (SLE), and systemic sclerosis. Particularly, non-hematopoietic cells, such as fibroblasts and cardiomyocytes, differentiated from autoimmune patient-derived iPS cells have shown promising results for further research into the pathogenesis. Recently, our groups established a method for differentiating dendritic cells that produce interferon-alpha, which can be applied as an SLE pathological model. In summary, patient-derived iPS cells can provide a promising platform for pathological research and new drug discovery in the field of immune-related diseases.

Rheumatoid Arthritis Has Won the Battle but Not the War: How Many Joints Will We Save Tomorrow?

Rheumatoid arthritis refers to joint diseases of unclear etiology whose final stages can lead to unbearable pain and complete immobility of the affected joints. As one of the most widely known diseases of the joints, it serves as a study target for a large number of research groups and pharmaceutical companies. Modern treatment with anti-inflammatory drugs, including janus kinase (JAK) inhibitors, monoclonal antibodies, and botanicals (polyphenols, glycosides, alkaloids, etc.) has achieved some success and hope for improving the course of the disease. However, existing drugs against RA have a number of side effects which push researchers to elaborate on more selective and effective drug candidates. The avant-garde of research, which aims to develop treatment of rheumatoid arthritis using antisense oligonucleotides along with nonsteroidal drugs and corticosteroids against inflammation, increases the chances of success and expands the arsenal of drugs. The primary goal in the treatment of this disease is to find therapies that allow patients with rheumatoid arthritis to move their joints without pain. The main purpose of this review is to show the victories and challenges for the treatment of rheumatoid arthritis and the tortuous but promising path of research that aims to help patients experience the joy of freely moving joints without pain.

RNA circuits and RNA-binding proteins in T cells

RNA is integral to the regulatory circuits that control cell identity and behavior. Cis-regulatory elements in mRNAs interact with RNA-binding proteins (RBPs) that can alter RNA sequence, stability, and translation into protein. Similarly, long noncoding RNAs (lncRNAs) scaffold ribonucleoprotein complexes that mediate transcriptional and post-transcriptional regulation of gene expression. Indeed, cell programming is fundamental to multicellular life and, in this era of cellular therapies, it is of particular interest in T cells. Here, we review key concepts and recent advances in our understanding of the RNA circuits and RBPs that govern mammalian T cell differentiation and immune function.