Prioritization of autoimmune disease-associated genetic variants that perturb regulatory element activity in T cells

Challenges in IBD Research 2024: Preclinical Human IBD Mechanisms

Preclinical human inflammatory bowel disease (IBD) mechanisms is one of 5 focus areas of the Challenges in IBD Research 2024 document, which also includes environmental triggers, novel technologies, precision medicine, and pragmatic clinical research. Herein, we provide a comprehensive overview of current gaps in inflammatory bowel diseases research that relate to preclinical research and deliver actionable approaches to address them with a focus on how these gaps can lead to advancements in IBD interception, remission, and restoration. The document is the result of multidisciplinary input from scientists, clinicians, patients, and funders and represents a valuable resource for patient-centric research prioritization. This preclinical human IBD mechanisms section identifies major research gaps whose investigation will elucidate pathways and mechanisms that can be targeted to address unmet medical needs in IBD. Research gaps were identified in the following areas: genetics, risk alleles, and epigenetics; the microbiome; cell states and interactions; barrier function; IBD complications (specifically fibrosis and stricturing); and extraintestinal manifestations. To address these gaps, we share specific opportunities for investigation for basic and translational scientists and identify priority actions.

Lack of evidence for GWAS signals of exfoliation glaucoma working via monogenic loss-of-function mutation in the nearest gene

PURPOSE

Exfoliation syndrome (XFS) is a systemic disease of elastin-rich tissues involving a deposition of fibrillar exfoliative material (XFM) in the anterior chamber of the eye, which can promote glaucoma. The purpose of this study was to create mice with CRISPR/Cas9-induced variations in candidate genes identified from human genome-wide association studies (GWAS) and screen them for indices of XFS.

METHODS

Variants predicted to be deleterious were sought in the Agpat1, Cacna1a, Loxl1, Pomp, Rbms3, Sema6a, and Tlcd5 genes of C57BL/6J mice using CRISPR/Cas9-based gene editing. Strains were phenotyped by slit-lamp, SD-OCT imaging, and fundus exams at 1-5 mos of age. Smaller cohorts of 12-mos-old mice were also studied.

RESULTS

Deleterious variants were identified in six targets; Pomp was recalcitrant to targeting. Multiple alleles of some targets were isolated, yielding 12 strains. Across all genotypes and ages, 277 mice were assessed by 902 slit-lamp exams, 928 SD-OCT exams, and 358 fundus exams. Homozygosity for Agpat1 or Cacna1a mutations led to early lethality; homozygosity for Loxl1 mutations led to pelvic organ prolapse, preventing aging. Loxl1 homozygotes exhibited a conjunctival phenotype of potential relevance to XFS. Multiple other genotype-specific phenotypes were variously identified. XFM was not observed in any mice.

CONCLUSIONS

This study did not detect XFM in any of the strains. This may have been due to species-specific differences, background dependence, or insufficient aging. Alternatively, it is possible that the current candidates, selected based on proximity to GWAS signals, are not effectors acting via monogenic loss-of-function mechanisms.

Functional dissection of complex and molecular trait variants at single nucleotide resolution

Identifying the causal variants and mechanisms that drive complex traits and diseases remains a core problem in human genetics. The majority of these variants have individually weak effects and lie in non-coding gene-regulatory elements where we lack a complete understanding of how single nucleotide alterations modulate transcriptional processes to affect human phenotypes. To address this, we measured the activity of 221,412 trait-associated variants that had been statistically fine-mapped using a Massively Parallel Reporter Assay (MPRA) in 5 diverse cell-types. We show that MPRA is able to discriminate between likely causal variants and controls, identifying 12,025 regulatory variants with high precision. Although the effects of these variants largely agree with orthogonal measures of function, only 69% can plausibly be explained by the disruption of a known transcription factor (TF) binding motif. We dissect the mechanisms of 136 variants using saturation mutagenesis and assign impacted TFs for 91% of variants without a clear canonical mechanism. Finally, we provide evidence that epistasis is prevalent for variants in close proximity and identify multiple functional variants on the same haplotype at a small, but important, subset of trait-associated loci. Overall, our study provides a systematic functional characterization of likely causal common variants underlying complex and molecular human traits, enabling new insights into the regulatory grammar underlying disease risk.

ITPRIPL1 binds CD3ε to impede T cell activation and enable tumor immune evasion

Cancer immunotherapy has transformed treatment possibilities, but its effectiveness differs significantly among patients, indicating the presence of alternative pathways for immune evasion. Here, we show that ITPRIPL1 functions as an inhibitory ligand of CD3ε, and its expression inhibits T cells in the tumor microenvironment. The binding of ITPRIPL1 extracellular domain to CD3ε on T cells significantly decreased calcium influx and ZAP70 phosphorylation, impeding initial T cell activation. Treatment with a neutralizing antibody against ITPRIPL1 restrained tumor growth and promoted T cell infiltration in mouse models across various solid tumor types. The antibody targeting canine ITPRIPL1 exhibited notable therapeutic efficacy against naturally occurring tumors in pet clinics. These findings highlight the role of ITPRIPL1 (or CD3L1, CD3ε ligand 1) in impeding T cell activation during the critical "signal one" phase. This discovery positions ITPRIPL1 as a promising therapeutic target against multiple tumor types.

MPRAVarDB: an online database and web server for exploring regulatory effects of genetic variants

Massively parallel reporter assay (MPRA) is an important technology to evaluate the impact of genetic variants on gene regulation. Here, we present MPRAVarDB, an online database and web server, for exploring regulatory effects of genetic variants. MPRAVarDB harbors 18 MPRA experiments designed to assess the regulatory effects of genetic variants associated with GWAS loci, eQTLs and various genomic features, resulting in a total of 242,818 variants tested across more than 30 cell lines and 30 human diseases or traits. MPRAVarDB empowers the query of MPRA variants by genomic region, disease and cell line or by any combination of these query terms. Notably, MPRAVarDB offers a suite of pretrained machine learning models tailored to the specific disease and cell line, facilitating the genome-wide prediction of regulatory variants. MPRAVarDB is friendly to use, and users only need a few clicks to receive query and prediction results.

Tissue-specific enhancer-gene maps from multimodal single-cell data identify causal disease alleles

Translating genome-wide association study (GWAS) loci into causal variants and genes requires accurate cell-type-specific enhancer-gene maps from disease-relevant tissues. Building enhancer-gene maps is essential but challenging with current experimental methods in primary human tissues. Here we developed a nonparametric statistical method, SCENT (single-cell enhancer target gene mapping), that models association between enhancer chromatin accessibility and gene expression in single-cell or nucleus multimodal RNA sequencing and ATAC sequencing data. We applied SCENT to 9 multimodal datasets including >120,000 single cells or nuclei and created 23 cell-type-specific enhancer-gene maps. These maps were highly enriched for causal variants in expression quantitative loci and GWAS for 1,143 diseases and traits. We identified likely causal genes for both common and rare diseases and linked somatic mutation hotspots to target genes. We demonstrate that application of SCENT to multimodal data from disease-relevant human tissue enables the scalable construction of accurate cell-type-specific enhancer-gene maps, essential for defining noncoding variant function.

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?

Systematic identification of genotype-dependent enhancer variants in eosinophilic esophagitis

Eosinophilic esophagitis (EoE) is a rare atopic disorder associated with esophageal dysfunction, including difficulty swallowing, food impaction, and inflammation, that develops in a small subset of people with food allergies. Genome-wide association studies (GWASs) have identified 9 independent EoE risk loci reaching genome-wide significance (p < 5 × 10-8) and 27 additional loci of suggestive significance (5 × 10-8 < p < 1 × 10-5). In the current study, we perform linkage disequilibrium (LD) expansion of these loci to nominate a set of 531 variants that are potentially causal. To systematically interrogate the gene regulatory activity of these variants, we designed a massively parallel reporter assay (MPRA) containing the alleles of each variant within their genomic sequence context cloned into a GFP reporter library. Analysis of reporter gene expression in TE-7, HaCaT, and Jurkat cells revealed cell-type-specific gene regulation. We identify 32 allelic enhancer variants, representing 6 genome-wide significant EoE loci and 7 suggestive EoE loci, that regulate reporter gene expression in a genotype-dependent manner in at least one cellular context. By annotating these variants with expression quantitative trait loci (eQTL) and chromatin looping data in related tissues and cell types, we identify putative target genes affected by genetic variation in individuals with EoE. Transcription factor enrichment analyses reveal possible roles for cell-type-specific regulators, including GATA3. Our approach reduces the large set of EoE-associated variants to a set of 32 with allelic regulatory activity, providing functional insights into the effects of genetic variation in this disease.

The role of transcription factors in shaping regulatory T cell identity

Forkhead box protein 3-expressing (FOXP3+) regulatory T cells (Treg cells) suppress conventional T cells and are essential for immunological tolerance. FOXP3, the master transcription factor of Treg cells, controls the expression of multiples genes to guide Treg cell differentiation and function. However, only a small fraction (<10%) of Treg cell-associated genes are directly bound by FOXP3, and FOXP3 alone is insufficient to fully specify the Treg cell programme, indicating a role for other accessory transcription factors operating upstream, downstream and/or concurrently with FOXP3 to direct Treg cell specification and specialized functions. Indeed, the heterogeneity of Treg cells can be at least partially attributed to differential expression of transcription factors that fine-tune their trafficking, survival and functional properties, some of which are niche-specific. In this Review, we discuss the emerging roles of accessory transcription factors in controlling Treg cell identity. We specifically focus on members of the basic helix-loop-helix family (AHR), basic leucine zipper family (BACH2, NFIL3 and BATF), CUT homeobox family (SATB1), zinc-finger domain family (BLIMP1, Ikaros and BCL-11B) and interferon regulatory factor family (IRF4), as well as lineage-defining transcription factors (T-bet, GATA3, RORγt and BCL-6). Understanding the imprinting of Treg cell identity and specialized function will be key to unravelling basic mechanisms of autoimmunity and identifying novel targets for drug development.

Using Synthetic DNA Libraries to Investigate Chromatin and Gene Regulation

Despite the recent explosion in genome-wide studies in chromatin and gene regulation, we are still far from extracting a set of genetic rules that can predict the function of the regulatory genome. One major reason for this deficiency is that gene regulation is a multi-layered process that involves an enormous variable space, which cannot be fully explored using native genomes. This problem can be partially solved by introducing synthetic DNA libraries into cells, a method that can test the regulatory roles of thousands to millions of sequences with limited variables. Here, we review recent applications of this method to study transcription factor (TF) binding, nucleosome positioning, and transcriptional activity. We discuss the design principles, experimental procedures, and major findings from these studies and compare the pros and cons of different approaches.

Non-coding sequence variation reveals fragility within interleukin 2 feedback circuitry and shapes autoimmune disease risk

Genetic variants associated with human autoimmune diseases commonly map to non-coding control regions, particularly enhancers that function selectively in immune cells and fine-tune gene expression within a relatively narrow range of values. How such modest, cell-type-selective changes can meaningfully shape organismal disease risk remains unclear. To explore this issue, we experimentally manipulated species-conserved enhancers within the disease-associated IL2RA locus and studied accompanying changes in the progression of autoimmunity. Perturbing distinct enhancers with restricted activity in conventional T cells (Tconvs) or regulatory T cells (Tregs)-two functionally antagonistic T cell subsets-caused only modest, cell-type-selective decreases in IL2ra expression parameters. However, these same perturbations had striking and opposing effects in vivo , completely preventing or severely accelerating disease in a murine model of type 1 diabetes. Quantitative tissue imaging and computational modelling revealed that each enhancer manipulation impinged on distinct IL-2-dependent feedback circuits. These imbalances altered the intracellular signaling and intercellular communication dynamics of activated Tregs and Tconvs, producing opposing spatial domains that amplified or constrained ongoing autoimmune responses. These findings demonstrate how subtle changes in gene regulation stemming from non-coding variation can propagate across biological scales due to non-linearities in intra- and intercellular feedback circuitry, dramatically shaping disease risk at the organismal level.

Pathophysiological role of BACH transcription factors in digestive system diseases

BTB and CNC homologous (BACH) proteins, including BACH1 and BACH2, are transcription factors that are widely expressed in human tissues. BACH proteins form heterodimers with small musculoaponeurotic fibrosarcoma (MAF) proteins to suppress the transcription of target genes. Furthermore, BACH1 promotes the transcription of target genes. BACH proteins regulate physiological processes, such as the differentiation of B cells and T cells, mitochondrial function, and heme homeostasis as well as pathogenesis related to inflammation, oxidative-stress damage caused by drugs, toxicants, or infections; autoimmunity disorders; and cancer angiogenesis, epithelial-mesenchymal transition, chemotherapy resistance, progression, and metabolism. In this review, we discuss the function of BACH proteins in the digestive system, including the liver, gallbladder, esophagus, stomach, small and large intestines, and pancreas. BACH proteins directly target genes or indirectly regulate downstream molecules to promote or inhibit biological phenomena such as inflammation, tumor angiogenesis, and epithelial-mesenchymal transition. BACH proteins are also regulated by proteins, miRNAs, LncRNAs, labile iron, and positive and negative feedback. Additionally, we summarize a list of regulators targeting these proteins. Our review provides a reference for future studies on targeted drugs in digestive diseases.

Common genetic factors among autoimmune diseases

Autoimmune diseases display a high degree of comorbidity within individuals and families, suggesting shared risk factors. Over the past 15 years, genome-wide association studies have established the polygenic basis of these common conditions and revealed widespread sharing of genetic effects, indicative of a shared immunopathology. Despite ongoing challenges in determining the precise genes and molecular consequences of these risk variants, functional experiments and integration with multimodal genomic data are providing valuable insights into key immune cells and pathways driving these diseases, with potential therapeutic implications. Moreover, genetic studies of ancient populations are shedding light on the contribution of pathogen-driven selection pressures to the increased prevalence of autoimmune disease. This Review summarizes the current understanding of autoimmune disease genetics, including shared effects, mechanisms, and evolutionary origins.

Hacking hematopoiesis - emerging tools for examining variant effects

Hematopoiesis is a continuous process of blood and immune cell production. It is orchestrated by thousands of gene products that respond to extracellular signals by guiding cell fate decisions to meet the needs of the organism. Although much of our knowledge of this process comes from work in model systems, we have learned a great deal from studies on human genetic variation. Considerable insight has emerged from studies on presumed monogenic blood disorders, which continue to provide key insights into the mechanisms critical for hematopoiesis. Furthermore, the emergence of large-scale biobanks and cohorts has uncovered thousands of genomic loci associated with blood cell traits and diseases. Some of these blood cell trait-associated loci act as modifiers of what were once thought to be monogenic blood diseases. However, most of these loci await functional validation. Here, we discuss the validation bottleneck and emerging methods to more effectively connect variant to function. In particular, we highlight recent innovations in genome editing, which have paved the path forward for high-throughput functional assessment of loci. Finally, we discuss existing barriers to progress, including challenges in manipulating the genomes of primary hematopoietic cells.

Widespread perturbation of ETS factor binding sites in cancer

Although >90% of somatic mutations reside in non-coding regions, few have been reported as cancer drivers. To predict driver non-coding variants (NCVs), we present a transcription factor (TF)-aware burden test based on a model of coherent TF function in promoters. We apply this test to NCVs from the Pan-Cancer Analysis of Whole Genomes cohort and predict 2555 driver NCVs in the promoters of 813 genes across 20 cancer types. These genes are enriched in cancer-related gene ontologies, essential genes, and genes associated with cancer prognosis. We find that 765 candidate driver NCVs alter transcriptional activity, 510 lead to differential binding of TF-cofactor regulatory complexes, and that they primarily impact the binding of ETS factors. Finally, we show that different NCVs within a promoter often affect transcriptional activity through shared mechanisms. Our integrated computational and experimental approach shows that cancer NCVs are widespread and that ETS factors are commonly disrupted.

Hereditable variants of classical protein tyrosine phosphatase genes: Will they prove innocent or guilty?

Protein tyrosine phosphatases, together with protein tyrosine kinases, control many molecular signaling steps that control life at cellular and organismal levels. Impairing alterations in the genes encoding the involved proteins is expected to profoundly affect the quality of life-if compatible with life at all. Here, we review the current knowledge on the effects of germline variants that have been reported for genes encoding a subset of the protein tyrosine phosphatase superfamily; that of the thirty seven classical members. The conclusion must be that the newest genome research tools produced an avalanche of data that suggest 'guilt by association' for individual genes to specific disorders. Future research should face the challenge to investigate these accusations thoroughly and convincingly, to reach a mature genotype-phenotype map for this intriguing protein family.

TFSyntax: a database of transcription factors binding syntax in mammalian genomes

In mammals, transcriptional factors (TFs) drive gene expression by binding to regulatory elements in a cooperative manner. Deciphering the rules of such cooperation is crucial to obtain a full understanding of cellular homeostasis and development. Although this is a long-standing topic, there is no comprehensive database for biologists to access the syntax of TF binding sites. Here we present TFSyntax (https://tfsyntax.zhaopage.com), a database focusing on the arrangement of TF binding sites. TFSyntax maps the binding motif of 1299 human TFs and 890 mouse TFs across 382 cells and tissues, representing the most comprehensive TF binding map to date. In addition to location, TFSyntax defines motif positional preference, density and colocalization within accessible elements. Powered by a series of functional modules based on web interface, users can freely search, browse, analyze, and download data of interest. With comprehensive characterization of TF binding syntax across distinct tissues and cell types, TFSyntax represents a valuable resource and platform for studying the mechanism of transcriptional regulation and exploring how regulatory DNA variants cause disease.

Genetic basis of defects in immune tolerance underlying the development of autoimmunity

Genetic variants associated with susceptibility to autoimmune disease have provided important insight into the mechanisms responsible for the loss of immune tolerance and the subsequent development of autoantibodies, tissue damage, and onset of clinical disease. Here, we review how genetic variants shared across multiple autoimmune diseases have contributed to our understanding of global tolerance failure, focusing on variants in the human leukocyte antigen region, PTPN2 and PTPN22, and their role in antigen presentation and T and B cell homeostasis. Variants unique to a specific autoimmune disease such as those in PADI2 and PADI4 that are associated with rheumatoid arthritis are also discussed, addressing their role in disease-specific immunopathology. Current research continues to focus on determining the functional consequences of autoimmune disease-associated variants but has recently expanded to variants in the non-coding regions of the genome using novel approaches to investigate the impact of these variants on mechanisms regulating gene expression. Lastly, studying genetic risk variants in the setting of autoimmunity has clinical implications, helping predict who will develop autoimmune disease and also identifying potential therapeutic targets.