Mapping effector genes at lupus GWAS loci using promoter Capture-C in follicular helper T cells

Dynamic chromatin architecture identifies new autoimmune-associated enhancers for IL2 and novel genes regulating CD4+ T cell activation

Genome-wide association studies (GWAS) have identified hundreds of genetic signals associated with autoimmune disease. The majority of these signals are located in non-coding regions and likely impact cis-regulatory elements (cRE). Because cRE function is dynamic across cell types and states, profiling the epigenetic status of cRE across physiological processes is necessary to characterize the molecular mechanisms by which autoimmune variants contribute to disease risk. We localized risk variants from 15 autoimmune GWAS to cRE active during TCR-CD28 co-stimulation of naïve human CD4+ T cells. To characterize how dynamic changes in gene expression correlate with cRE activity, we measured transcript levels, chromatin accessibility, and promoter-cRE contacts across three phases of naive CD4+ T cell activation using RNA-seq, ATAC-seq, and HiC. We identified ~1200 protein-coding genes physically connected to accessible disease-associated variants at 423 GWAS signals, at least one-third of which are dynamically regulated by activation. From these maps, we functionally validated a novel stretch of evolutionarily conserved intergenic enhancers whose activity is required for activation-induced IL2 gene expression in human and mouse, and is influenced by autoimmune-associated genetic variation. The set of genes implicated by this approach are enriched for genes controlling CD4+ T cell function and genes involved in human inborn errors of immunity, and we pharmacologically validated eight implicated genes as novel regulators of T cell activation. These studies directly show how autoimmune variants and the genes they regulate influence processes involved in CD4+ T cell proliferation and activation.

Implicating type 2 diabetes effector genes in relevant metabolic cellular models using promoter-focused Capture-C

AIMS/HYPOTHESIS

Genome-wide association studies (GWAS) have identified hundreds of type 2 diabetes loci, with the vast majority of signals located in non-coding regions; as a consequence, it remains largely unclear which 'effector' genes these variants influence. Determining these effector genes has been hampered by the relatively challenging cellular settings in which they are hypothesised to confer their effects.

METHODS

To implicate such effector genes, we elected to generate and integrate high-resolution promoter-focused Capture-C, assay for transposase-accessible chromatin with sequencing (ATAC-seq) and RNA-seq datasets to characterise chromatin and expression profiles in multiple cell lines relevant to type 2 diabetes for subsequent functional follow-up analyses: EndoC-BH1 (pancreatic beta cell), HepG2 (hepatocyte) and Simpson-Golabi-Behmel syndrome (SGBS; adipocyte).

RESULTS

The subsequent variant-to-gene analysis implicated 810 candidate effector genes at 370 type 2 diabetes risk loci. Using partitioned linkage disequilibrium score regression, we observed enrichment for type 2 diabetes and fasting glucose GWAS loci in promoter-connected putative cis-regulatory elements in EndoC-BH1 cells as well as fasting insulin GWAS loci in SGBS cells. Moreover, as a proof of principle, when we knocked down expression of the SMCO4 gene in EndoC-BH1 cells, we observed a statistically significant increase in insulin secretion.

CONCLUSIONS/INTERPRETATION

These results provide a resource for comparing tissue-specific data in tractable cellular models as opposed to relatively challenging primary cell settings.

DATA AVAILABILITY

Raw and processed next-generation sequencing data for EndoC-BH1, HepG2, SGBS_undiff and SGBS_diff cells are deposited in GEO under the Superseries accession GSE262484. Promoter-focused Capture-C data are deposited under accession GSE262496. Hi-C data are deposited under accession GSE262481. Bulk ATAC-seq data are deposited under accession GSE262479. Bulk RNA-seq data are deposited under accession GSE262480.

Variant-to-function mapping of late-onset Alzheimer's disease GWAS signals in human microglial cell models implicates RTFDC1 at the CASS4 locus

Late-onset Alzheimer's disease (LOAD) research has principally focused on neurons over the years due to their known role in the production of amyloid beta plaques and neurofibrillary tangles. In contrast, recent genomic studies of LOAD have implicated microglia as culprits of the prolonged inflammation exacerbating the neurodegeneration observed in patient brains. Indeed, recent LOAD genome-wide association studies (GWAS) have reported multiple loci near genes related to microglial function, including TREM2, ABI3, and CR1. However, GWAS alone cannot pinpoint underlying causal variants or effector genes at such loci, as most signals reside in non-coding regions of the genome and could presumably confer their influence frequently via long-range regulatory interactions. We elected to carry out a combination of ATAC-seq and high-resolution promoter-focused Capture-C in two human microglial cell models (iPSC-derived microglia and HMC3) in order to physically map interactions between LOAD GWAS-implicated candidate causal variants and their corresponding putative effector genes. Notably, we observed consistent evidence that rs6024870 at the GWAS CASS4 locus contacted the promoter of nearby gene, RTFDC1. We subsequently observed a directionallly consistent decrease in RTFDC1 expression with the the protective minor A allele of rs6024870 via both luciferase assays in HMC3 cells and expression studies in primary human microglia. Through CRISPR-Cas9-mediated deletion of the putative regulatory region harboring rs6024870 in HMC3 cells, we observed increased pro-inflammatory cytokine secretion and decreased DNA double strand break repair related, at least in part, to RTFDC1 expression levels. Our variant-to-function approach therefore reveals that the rs6024870-harboring regulatory element at the LOAD 'CASS4' GWAS locus influences both microglial inflammatory capacity and DNA damage resolution, along with cumulative evidence implicating RTFDC1 as a novel candidate effector gene.

3D genomic features across >50 diverse cell types reveal insights into the genomic architecture of childhood obesity

The prevalence of childhood obesity is increasing worldwide, along with the associated common comorbidities of type 2 diabetes and cardiovascular disease in later life. Motivated by evidence for a strong genetic component, our prior genome-wide association study (GWAS) efforts for childhood obesity revealed 19 independent signals for the trait; however, the mechanism of action of these loci remains to be elucidated. To molecularly characterize these childhood obesity loci we sought to determine the underlying causal variants and the corresponding effector genes within diverse cellular contexts. Integrating childhood obesity GWAS summary statistics with our existing 3D genomic datasets for 57 human cell types, consisting of high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq, and RNA-seq, we applied stratified LD score regression and calculated the proportion of genome-wide SNP heritability attributable to cell type-specific features, revealing pancreatic alpha cell enrichment as the most statistically significant. Subsequent chromatin contact-based fine-mapping was carried out for genome-wide significant childhood obesity loci and their linkage disequilibrium proxies to implicate effector genes, yielded the most abundant number of candidate variants and target genes at the BDNF, ADCY3, TMEM18 and FTO loci in skeletal muscle myotubes and the pancreatic beta-cell line, EndoC-BH1. One novel implicated effector gene, ALKAL2 - an inflammation-responsive gene in nerve nociceptors - was observed at the key TMEM18 locus across multiple immune cell types. Interestingly, this observation was also supported through colocalization analysis using expression quantitative trait loci (eQTL) derived from the Genotype-Tissue Expression (GTEx) dataset, supporting an inflammatory and neurologic component to the pathogenesis of childhood obesity. Our comprehensive appraisal of 3D genomic datasets generated in a myriad of different cell types provides genomic insights into pediatric obesity pathogenesis.

3D chromatin-based variant-to-gene maps across 57 human cell types reveal the cellular and genetic architecture of autoimmune disease susceptibility

A portion of the genetic basis for many common autoimmune disorders has been uncovered by genome-wide association studies (GWAS), but GWAS do not reveal causal variants, effector genes, or the cell types impacted by disease-associated variation. We have generated 3D genomic datasets consisting of promoter-focused Capture-C, Hi-C, ATAC-seq, and RNA-seq and integrated these data with GWAS of 16 autoimmune traits to physically map disease-associated variants to the effector genes they likely regulate in 57 human cell types. These 3D maps of gene cis-regulatory architecture are highly powered to identify the cell types most likely impacted by disease-associated genetic variation compared to 1D genomic features, and tend to implicate different effector genes than eQTL approaches in the same cell types. Most of the variants implicated by these cis-regulatory architectures are highly trait-specific, but nearly half of the target genes connected to these variants are shared across multiple autoimmune disorders in multiple cell types, suggesting a high level of genetic diversity and complexity among autoimmune diseases that nonetheless converge at the level of target gene and cell type. Substantial effector gene sharing led to the common enrichment of similar biological networks across disease and cell types. However, trait-specific pathways representing potential areas for disease-specific intervention were identified. To test this, we pharmacologically validated squalene synthase, a cholesterol biosynthetic enzyme encoded by the FDFT1 gene implicated by our approach in MS and SLE, as a novel immunomodulatory drug target controlling inflammatory cytokine production by human T cells. These data represent a comprehensive resource for basic discovery of gene cis-regulatory mechanisms, and the analyses reported reveal mechanisms by which autoimmune-associated variants act to regulate gene expression, function, and pathology across multiple, distinct tissues and cell types.

Shared and unique 3D genomic features of substance use disorders across multiple cell types

Recent genome-wide association studies (GWAS) have revealed shared genetic components among alcohol, opioid, tobacco and cannabis use disorders. However, the extent of the underlying shared causal variants and effector genes, along with their cellular context, remain unclear. We leveraged our existing 3D genomic datasets comprising high-resolution promoter-focused Capture-C/Hi-C, ATAC-seq and RNA-seq across >50 diverse human cell types to focus on genomic regions that coincide with GWAS loci. Using stratified LD regression, we determined the proportion of genomewide SNP heritability attributable to the features assayed across our cell types by integrating recent GWAS summary statistics for the relevant traits: alcohol use disorder (AUD), tobacco use disorder (TUD), opioid use disorder (OUD) and cannabis use disorder (CanUD). Statistically significant enrichments (P<0.05) were observed in 14 specific cell types, with heritability reaching 9.2-fold for iPSC-derived cortical neurons and neural progenitors, confirming that they are crucial cell types for further functional exploration. Additionally, several pancreatic cell types, notably pancreatic beta cells, showed enrichment for TUD, with heritability enrichments up to 4.8-fold, suggesting genomic overlap with metabolic processes. Further investigation revealed significant positive genetic correlations between T2D with both TUD and CanUD (FDR<0.05) and a significant negative genetic correlation with AUD. Interestingly, after partitioning the heritability for each cell type's cis-regulatory elements, the correlation between T2D and TUD for pancreatic beta cells was greater (r=0.2) than the global genetic correlation value. Our study provides new genomic insights into substance use disorders and implicates cell types where functional follow-up studies could reveal causal variant-gene mechanisms underpinning these disorders.

Perturbation of the insomnia WDR90 genome-wide association studies locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q

Although genome-wide association studies (GWAS) have identified loci for sleep-related traits, they do not directly uncover the underlying causal variants and corresponding effector genes. The majority of such variants reside in non-coding regions and are therefore presumed to impact cis-regulatory elements. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated phosphatidyl inositol glycan (PIG)-Q as a functionally relevant gene at the insomnia "WDR90" GWAS locus. However, importantly that effort did not characterize the corresponding underlying causal variant. Specifically, our previous 3D genomic datasets nominated a shortlist of three neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium within an intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. We sought to investigate the influence of these SNPs collectively and then individually on PIG-Q modulation to pinpoint the causal "regulatory" variant. Starting with gross level perturbation, deletion of the entire region in NPCs via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from individual luciferase reporter assays for each SNP in iPSCs revealed that the region with the rs3752495 risk allele (RA) induced a ~2.5-fold increase in luciferase expression. Importantly, rs3752495 also exhibited an allele-specific effect, with the RA increasing the luciferase expression by ~2-fold versus the non-RA. In conclusion, our variant-to-function approach and in vitro validation implicate rs3752495 as a causal insomnia variant embedded within WDR90 while modulating the expression of the distally located PIG-Q.

Systematic decoding of cis gene regulation defines context-dependent control of the multi-gene costimulatory receptor locus in human T cells

Cis-regulatory elements (CREs) interact with trans regulators to orchestrate gene expression, but how transcriptional regulation is coordinated in multi-gene loci has not been experimentally defined. We sought to characterize the CREs controlling dynamic expression of the adjacent costimulatory genes CD28, CTLA4 and ICOS, encoding regulators of T cell-mediated immunity. Tiling CRISPR interference (CRISPRi) screens in primary human T cells, both conventional and regulatory subsets, uncovered gene-, cell subset- and stimulation-specific CREs. Integration with CRISPR knockout screens and assay for transposase-accessible chromatin with sequencing (ATAC-seq) profiling identified trans regulators influencing chromatin states at specific CRISPRi-responsive elements to control costimulatory gene expression. We then discovered a critical CCCTC-binding factor (CTCF) boundary that reinforces CRE interaction with CTLA4 while also preventing promiscuous activation of CD28. By systematically mapping CREs and associated trans regulators directly in primary human T cell subsets, this work overcomes longstanding experimental limitations to decode context-dependent gene regulatory programs in a complex, multi-gene locus critical to immune homeostasis.

Foxp3 depends on Ikaros for control of regulatory T cell gene expression and function

Ikaros is a transcriptional factor required for conventional T cell development, differentiation, and anergy. While the related factors Helios and Eos have defined roles in regulatory T cells (Treg), a role for Ikaros has not been established. To determine the function of Ikaros in the Treg lineage, we generated mice with Treg-specific deletion of the Ikaros gene (Ikzf1). We find that Ikaros cooperates with Foxp3 to establish a major portion of the Treg epigenome and transcriptome. Ikaros-deficient Treg exhibit Th1-like gene expression with abnormal production of IL-2, IFNg, TNFa, and factors involved in Wnt and Notch signaling. While Ikzf1-Treg-cko mice do not develop spontaneous autoimmunity, Ikaros-deficient Treg are unable to control conventional T cell-mediated immune pathology in response to TCR and inflammatory stimuli in models of IBD and organ transplantation. These studies establish Ikaros as a core factor required in Treg for tolerance and the control of inflammatory immune responses.

GWAS-informed data integration and non-coding CRISPRi screen illuminate genetic etiology of bone mineral density

Over 1,100 independent signals have been identified with genome-wide association studies (GWAS) for bone mineral density (BMD), a key risk factor for mortality-increasing fragility fractures; however, the effector gene(s) for most remain unknown. Informed by a variant-to-gene mapping strategy implicating 89 non-coding elements predicted to regulate osteoblast gene expression at BMD GWAS loci, we executed a single-cell CRISPRi screen in human fetal osteoblast 1.19 cells (hFOBs). The BMD relevance of hFOBs was supported by heritability enrichment from cross-cell type stratified LD-score regression involving 98 cell types grouped into 15 tissues. 24 genes showed perturbation in the screen, with four (ARID5B, CC2D1B, EIF4G2, and NCOA3) exhibiting consistent effects upon siRNA knockdown on three measures of osteoblast maturation and mineralization. Lastly, additional heritability enrichments, genetic correlations, and multi-trait fine-mapping revealed that many BMD GWAS signals are pleiotropic and likely mediate their effects via non-bone tissues that warrant attention in future screens.

Common variable immunodeficiency, cross currents, and prevailing winds

Common variable immunodeficiency (CVID) is a heterogenous disease category created to distinguish late-onset antibody deficiencies from early-onset diseases like agammaglobulinemia or more expansively dysfunctional combined immunodeficiencies. Opinions vary on which affected patients should receive a CVID diagnosis which confuses clinicians and erects reproducibility barriers for researchers. Most experts agree that CVID's most indeliable feature is defective germinal center (GC) production of isotype-switched, affinity-maturated antibodies. Here, we review the biological factors contributing to CVID-associated GC dysfunction including genetic, epigenetic, tolerogenic, microbiome, and regulatory abnormalities. We also discuss the consequences of these biological phenomena to the development of non-infectious disease complications. Finally, we opine on topics and lines of investigation we think hold promise for expanding our mechanistic understanding of this protean condition and for improving the lives of affected patients.

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.

Proteomic analyses reveal cystatin c is a promising biomarker for evaluation of systemic lupus erythematosus

BACKGROUND

Systemic lupus erythematosus (SLE) is an autoimmune disease with multiple organ involvement, especially the kidneys. However, the underlying mechanism remains unclear, and accurate biomarkers are still lacking. This study aimed to identify biomarkers to assess organ damage and disease activity in patients with SLE using quantitative proteomics.

METHODS

Proteomic analysis was performed using mass spectrometry in 15 patients with SLE and 15 age-matched healthy controls. Proteomic profiles were compared in four main subtypes: SLE with proteinuria (SLE-PN), SLE without proteinuria (SLE-non-PN), SLE with anti-dsDNA positivity (SLE-DP), and SLE with anti-dsDNA negativity (SLE-non-DP). Gene ontology biological process analysis revealed differentially expressed protein networks. Cystatin C (CysC) levels were measured in 200 patients with SLE using an immunoturbidimetric assay. Clinical and laboratory data were collected to assess their correlation with serum CysC levels.

RESULTS

Proteomic analysis showed that upregulated proteins in both the SLE-PN and SLE-DP groups were mainly mapped to neutrophil activation networks. Moreover, CysC from neutrophil activation networks was upregulated in both the SLE-PN and SLE-DP groups. The associations of serum CysC level with proteinuria, anti-dsDNA positivity, lower complement C3 levels, and SLE disease activity index score in patients with SLE were further validated in a large independent cohort.

CONCLUSIONS

Neutrophil activation is more prominent in SLE with proteinuria and anti-dsDNA positivity, and CysC is a promising marker for monitoring organ damage and disease activity in SLE.

Integrative Genomic Analyses Identify LncRNA Regulatory Networks across Pediatric Leukemias and Solid Tumors

Long noncoding RNAs (lncRNA) play an important role in gene regulation and contribute to tumorigenesis. While pan-cancer studies of lncRNA expression have been performed for adult malignancies, the lncRNA landscape across pediatric cancers remains largely uncharted. Here, we curated RNA sequencing data for 1,044 pediatric leukemia and extracranial solid tumors and integrated paired tumor whole genome sequencing and epigenetic data in relevant cell line models to explore lncRNA expression, regulation, and association with cancer. A total of 2,657 lncRNAs were robustly expressed across six pediatric cancers, including 1,142 exhibiting histotype-elevated expression. DNA copy number alterations contributed to lncRNA dysregulation at a proportion comparable to protein coding genes. Application of a multidimensional framework to identify and prioritize lncRNAs impacting gene networks revealed that lncRNAs dysregulated in pediatric cancer are associated with proliferation, metabolism, and DNA damage hallmarks. Analysis of upstream regulation via cell type-specific transcription factors further implicated distinct histotype-elevated and developmental lncRNAs. Integration of these analyses prioritized lncRNAs for experimental validation, and silencing of TBX2-AS1, the top-prioritized neuroblastoma-specific lncRNA, resulted in significant growth inhibition of neuroblastoma cells, confirming the computational predictions. Taken together, these data provide a comprehensive characterization of lncRNA regulation and function in pediatric cancers and pave the way for future mechanistic studies.

SIGNIFICANCE

Comprehensive characterization of lncRNAs in pediatric cancer leads to the identification of highly expressed lncRNAs across childhood cancers, annotation of lncRNAs showing histotype-specific elevated expression, and prediction of lncRNA gene regulatory networks.

Epromoters are new players in the regulatory landscape with potential pleiotropic roles

Precise spatiotemporal control of gene expression during normal development and cell differentiation is achieved by the combined action of proximal (promoters) and distal (enhancers) cis-regulatory elements. Recent studies have reported that a subset of promoters, termed Epromoters, works also as enhancers to regulate distal genes. This new paradigm opened novel questions regarding the complexity of our genome and raises the possibility that genetic variation within Epromoters has pleiotropic effects on various physiological and pathological traits by differentially impacting multiple proximal and distal genes. Here, we discuss the different observations pointing to an important role of Epromoters in the regulatory landscape and summarize the evidence supporting a pleiotropic impact of these elements in disease. We further hypothesize that Epromoter might represent a major contributor to phenotypic variation and disease.

Perturbation of the insomnia WDR90 GWAS locus pinpoints rs3752495 as a causal variant influencing distal expression of neighboring gene, PIG-Q

Although genome wide association studies (GWAS) have been crucial for the identification of loci associated with sleep traits and disorders, the method itself does not directly uncover the underlying causal variants and corresponding effector genes. The overwhelming majority of such variants reside in non-coding regions and are therefore presumed to impact the activity of cis-regulatory elements, such as enhancers. Our previously reported 'variant-to-gene mapping' effort in human induced pluripotent stem cell (iPSC)-derived neural progenitor cells (NPCs), combined with validation in both Drosophila and zebrafish, implicated PIG-Q as a functionally relevant gene at the insomnia 'WDR90' locus. However, importantly that effort did not characterize the corresponding underlying causal variant at this GWAS signal. Specifically, our genome-wide ATAC-seq and high-resolution promoter-focused Capture C datasets generated in this cell setting brought our attention to a shortlist of three tightly neighboring single nucleotide polymorphisms (SNPs) in strong linkage disequilibrium in a candidate intronic enhancer region of WDR90 that contacted the open PIG-Q promoter. The objective of this study was to investigate the influence of the proxy SNPs collectively and then individually on PIG-Q modulation and to pinpoint the causal "regulatory" variant among the three SNPs. Starting at a gross level perturbation, deletion of the entire region harboring all three SNPs in human iPSC-derived neural progenitor cells via CRISPR-Cas9 editing and subsequent RNA sequencing revealed expression changes in specific PIG-Q transcripts. Results from more refined individual luciferase reporter assays for each of the three SNPs in iPSCs revealed that the intronic region with the rs3752495 risk allele induced a ~2.5-fold increase in luciferase expression (n=10). Importantly, rs3752495 also exhibited an allele specific effect, with the risk allele increasing the luciferase expression by ~2-fold compared to the non-risk allele. In conclusion, our variant-to-function approach and subsequent in vitro validation implicates rs3752495 as a causal insomnia risk variant embedded at the WDR90-PIG-Q locus.

3D genome organization and epigenetic regulation in autoimmune diseases

Three-dimensional (3D) genomics is an emerging field of research that investigates the relationship between gene regulatory function and the spatial structure of chromatin. Chromatin folding can be studied using chromosome conformation capture (3C) technology and 3C-based derivative sequencing technologies, including chromosome conformation capture-on-chip (4C), chromosome conformation capture carbon copy (5C), and high-throughput chromosome conformation capture (Hi-C), which allow scientists to capture 3D conformations from a single site to the entire genome. A comprehensive analysis of the relationships between various regulatory components and gene function also requires the integration of multi-omics data such as genomics, transcriptomics, and epigenomics. 3D genome folding is involved in immune cell differentiation, activation, and dysfunction and participates in a wide range of diseases, including autoimmune diseases. We describe hierarchical 3D chromatin organization in this review and conclude with characteristics of C-techniques and multi-omics applications of the 3D genome. In addition, we describe the relationship between 3D genome structure and the differentiation and maturation of immune cells and address how changes in chromosome folding contribute to autoimmune diseases.

Variant-to-gene, gene-to-trait: finding new sleep genes through GWAS

Distilling insomnia genome-wide association study (GWAS) variants, Palermo and colleagues identified several genes that participate in sleep regulation in two different model organisms. This workflow sets off an innovative strategy to extract biological relevance from large human genomic databases.

3D genome alterations in T cells associated with disease activity of systemic lupus erythematosus

OBJECTIVES

Three-dimensional (3D) genome alterations can dysregulate gene expression by rewiring physical interactions within chromosomes in a tissue-specific or cell-specific manner and lead to diseases. We aimed to elucidate the 3D genome structure and its role in gene expression networks dysregulated in systemic lupus erythematosus (SLE).

METHODS

We performed Hi-C experiments using CD4+ T cells from 7 patients with SLE and 5 age-matched and sex-matched healthy controls (HCs) combined with RNA sequencing analysis. Further integrative analyses, including transcription factor motif enrichment, SPI1 knockdown and histone modifications (H3K27ac, H3K4me1, H3K4me3), were performed for altered loop-associated gene loci in SLE.

RESULTS

We deciphered the 3D chromosome organisation in T cells of patients with SLE and found it was clearly distinct from that of HCs and closely associated with the disease activity of SLE. Importantly, we identified loops within chromosomes associated with the disease activity of SLE and differentially expressed genes and found some key histone modifications close to these loops. Moreover, we demonstrated the contribution of the transcription factor SPI1, whose motif is located in the altered loop in SLE, to the overexpression of interferon pathway gene. In addition, we identified the potential influences of genetic variations in 3D genome alterations in SLE.

CONCLUSIONS

Our results highlight the 3D genome structure alterations associated with SLE development and provide a foundation for future interrogation of the relationships between chromosome structure and gene expression control in SLE.

Variant-to-gene mapping followed by cross-species genetic screening identifies GPI-anchor biosynthesis as a regulator of sleep

Genome-wide association studies (GWAS) in humans have identified loci robustly associated with several heritable diseases or traits, yet little is known about the functional roles of the underlying causal variants in regulating sleep duration or quality. We applied an ATAC-seq/promoter focused Capture C strategy in human iPSC-derived neural progenitors to carry out a "variant-to-gene" mapping campaign that identified 88 candidate sleep effector genes connected to relevant GWAS signals. To functionally validate the role of the implicated effector genes in sleep regulation, we performed a neuron-specific RNA interference screen in the fruit fly, Drosophila melanogaster, followed by validation in zebrafish. This approach identified a number of genes that regulate sleep including a critical role for glycosylphosphatidylinositol (GPI)-anchor biosynthesis. These results provide the first physical variant-to-gene mapping of human sleep genes followed by a model organism-based prioritization, revealing a conserved role for GPI-anchor biosynthesis in sleep regulation.

T follicular helper cells in autoimmune diseases

T follicular helper (Tfh) cells with the phenotype of mainly expressing surface molecules C-X-C motif chemokine receptor type 5 (CXCR5), inducible co-stimulator (ICOS), secreting cytokine interleukin-21 (IL-21) and requiring the transcription factor B cell lymphoma 6 (BCL-6) have been recently defined as a new subset of CD4⁺ T cells. They exist in germinal centers (GCs) of lymphoid organs and in peripheral blood. With the ability to promote B cell development, GC formation and antibody production, Tfh cells play critical roles in the pathogenesis of many autoimmune diseases, such as systemic lupus erythematosus (SLE), rheumatoid arthritis (RA), primary Sjögren's syndrome (pSS), etc. The aberrant proliferation and function of Tfh cells will cause the pathological process like autoantibody production and tissue injury. In this paper, we review the recent advances in Tfh cell biology and their roles in autoimmune diseases, with a mention of their use as therapeutic targets, which will shed more light on the pathogenesis and treatment of certain autoimmune diseases.

3D genome organization links non-coding disease-associated variants to genes

Genome sequencing has revealed over 300 million genetic variations in human populations. Over 90% of variants are single nucleotide polymorphisms (SNPs), the remainder include short deletions or insertions, and small numbers of structural variants. Hundreds of thousands of these variants have been associated with specific phenotypic traits and diseases through genome wide association studies which link significant differences in variant frequencies with specific phenotypes among large groups of individuals. Only 5% of disease-associated SNPs are located in gene coding sequences, with the potential to disrupt gene expression or alter of the function of encoded proteins. The remaining 95% of disease-associated SNPs are located in non-coding DNA sequences which make up 98% of the genome. The role of non-coding, disease-associated SNPs, many of which are located at considerable distances from any gene, was at first a mystery until the discovery that gene promoters regularly interact with distal regulatory elements to control gene expression. Disease-associated SNPs are enriched at the millions of gene regulatory elements that are dispersed throughout the non-coding sequences of the genome, suggesting they function as gene regulation variants. Assigning specific regulatory elements to the genes they control is not straightforward since they can be millions of base pairs apart. In this review we describe how understanding 3D genome organization can identify specific interactions between gene promoters and distal regulatory elements and how 3D genomics can link disease-associated SNPs to their target genes. Understanding which gene or genes contribute to a specific disease is the first step in designing rational therapeutic interventions.

Biological insights into systemic lupus erythematosus through an immune cell-specific transcriptome-wide association study

OBJECTIVE

Genome-wide association studies (GWAS) have identified >100 risk loci for systemic lupus erythematosus (SLE), but the disease genes at most loci remain unclear, hampering translation of these genetic discoveries. We aimed to prioritise genes underlying the 110 SLE loci that were identified in the latest East Asian GWAS meta-analysis.

METHODS

We built gene expression predictive models in blood B cells, CD4+ and CD8+ T cells, monocytes, natural killer cells and peripheral blood cells of 105 Japanese individuals. We performed a transcriptome-wide association study (TWAS) using data from the latest genome-wide association meta-analysis of 208 370 East Asians and searched for candidate genes using TWAS and three data-driven computational approaches.

RESULTS

TWAS identified 171 genes for SLE (p<1.0×10-5); 114 (66.7%) showed significance only in a single cell type; 127 (74.3%) were in SLE GWAS loci. TWAS identified a strong association between CD83 and SLE (p<7.7×10-8). Meta-analysis of genetic associations in the existing 208 370 East Asian and additional 1498 cases and 3330 controls found a novel single-variant association at rs72836542 (OR=1.11, p=4.5×10-9) around CD83. For the 110 SLE loci, we identified 276 gene candidates, including 104 genes at recently-identified SLE novel loci. We demonstrated in vitro that putative causal variant rs61759532 exhibited an allele-specific regulatory effect on ACAP1, and that presence of the SLE risk allele decreased ACAP1 expression.

CONCLUSIONS

Cell-level TWAS in six types of immune cells complemented SLE gene discovery and guided the identification of novel genetic associations. The gene findings shed biological insights into SLE genetic associations.

Role of CNSs Conserved Distal Cis-Regulatory Elements in CD4 + T Cell Development and Differentiation

Naïve CD4⁺ T cells differentiate into diverse subsets of effector cells and perform various homeostatic and immune functions. The differentiation and maintenance of these different subsets are controlled through the upregulation and silencing of master genes. Mechanistic studies of the regulation of these master genes identified conserved and distal intronic regulatory elements, which are accessible subsets of conserved non-coding sequences (CNSs), acting as cis-regulatory elements in a lineage-specific manner that controls the function of CD4⁺ T cells. Abnormal CNS activity is associated with incorrect expression of master genes and development of autoimmune diseases or immune suppression. Here, we describe the function of several conserved, distal cis-regulatory elements at the Foxp3, Rorc, Il-4, Il-10 and Il-17 gene locus were shown to play important roles in CD4⁺ T cells differentiation. Together, this review briefly outlines currently known CNSs, with a focus on their regulations and functions in complexes modulating the differentiation and maintenance of various CD4⁺ T cells subsets, in health and disease contexts, as well as during the conversion of T regulatory cells to T helper 17 cells. This article will provide a comprehensive view of CNSs conserved distal cis-regulatory elements at a few loci that control aspects of CD4⁺ T cells function.

Implicating effector genes at COVID-19 GWAS loci using promoter-focused Capture-C in disease-relevant immune cell types

BACKGROUND

SARS-CoV-2 infection results in a broad spectrum of COVID-19 disease, from mild or no symptoms to hospitalization and death. COVID-19 disease severity has been associated with some pre-existing conditions and the magnitude of the adaptive immune response to SARS-CoV-2, and a recent genome-wide association study (GWAS) of the risk of critical illness revealed a significant genetic component. To gain insight into how human genetic variation attenuates or exacerbates disease following SARS-CoV-2 infection, we implicated putatively functional COVID risk variants in the cis-regulatory landscapes of human immune cell types with established roles in disease severity and used high-resolution chromatin conformation capture to map these disease-associated elements to their effector genes.

RESULTS

This functional genomic approach implicates 16 genes involved in viral replication, the interferon response, and inflammation. Several of these genes (PAXBP1, IFNAR2, OAS1, OAS3, TNFAIP8L1, GART) were differentially expressed in immune cells from patients with severe versus moderate COVID-19 disease, and we demonstrate a previously unappreciated role for GART in T cell-dependent antibody-producing B cell differentiation in a human tonsillar organoid model.

CONCLUSIONS

This study offers immunogenetic insight into the basis of COVID-19 disease severity and implicates new targets for therapeutics that limit SARS-CoV-2 infection and its resultant life-threatening inflammation.

Deciphering the Role of 3D Genome Organization in Breast Cancer Susceptibility

Cancer risk by environmental exposure is modulated by an individual's genetics and age at exposure. This age-specific period of susceptibility is referred to as the "Window of Susceptibility" (WOS). Rats have a similar WOS for developing breast cancer. A previous study in rat identified an age-specific long-range regulatory interaction for the cancer gene, Pappa, that is associated with breast cancer susceptibility. However, the global role of three-dimensional genome organization and downstream gene expression programs in the WOS is not known. Therefore, we generated Hi-C and RNA-seq data in rat mammary epithelial cells within and outside the WOS. To systematically identify higher-order changes in 3D genome organization, we developed NE-MVNMF that combines network enhancement followed by multitask non-negative matrix factorization. We examined three-dimensional genome organization dynamics at the level of individual loops as well as higher-order domains. Differential chromatin interactions tend to be associated with differentially up-regulated genes with the WOS and recapitulate several human SNP-gene interactions associated with breast cancer susceptibility. Our approach identified genomic blocks of regions with greater overall differences in contact count between the two time points when the cluster assignments change and identified genes and pathways implicated in early carcinogenesis and cancer treatment. Our results suggest that WOS-specific changes in 3D genome organization are linked to transcriptional changes that may influence susceptibility to breast cancer.

Super interactive promoters provide insight into cell type-specific regulatory networks in blood lineage cell types

Existing studies of chromatin conformation have primarily focused on potential enhancers interacting with gene promoters. By contrast, the interactivity of promoters per se, while equally critical to understanding transcriptional control, has been largely unexplored, particularly in a cell type-specific manner for blood lineage cell types. In this study, we leverage promoter capture Hi-C data across a compendium of blood lineage cell types to identify and characterize cell type-specific super-interactive promoters (SIPs). Notably, promoter-interacting regions (PIRs) of SIPs are more likely to overlap with cell type-specific ATAC-seq peaks and GWAS variants for relevant blood cell traits than PIRs of non-SIPs. Moreover, PIRs of cell-type-specific SIPs show enriched heritability of relevant blood cell trait (s), and are more enriched with GWAS variants associated with blood cell traits compared to PIRs of non-SIPs. Further, SIP genes tend to express at a higher level in the corresponding cell type. Importantly, SIP subnetworks incorporating cell-type-specific SIPs and ATAC-seq peaks help interpret GWAS variants. Examples include GWAS variants associated with platelet count near the megakaryocyte SIP gene EPHB3 and variants associated lymphocyte count near the native CD4 T-Cell SIP gene ETS1. Interestingly, around 25.7% ~ 39.6% blood cell traits GWAS variants residing in SIP PIR regions disrupt transcription factor binding motifs. Importantly, our analysis shows the potential of using promoter-centric analyses of chromatin spatial organization data to identify biologically important genes and their regulatory regions.

Recent advances in understanding the genetic basis of systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is a polygenic chronic autoimmune disease leading to multiple organ damage. A large heritability of up to 66% is estimated in SLE, with roughly 180 reported susceptibility loci that have been identified mostly by genome-wide association studies (GWASs) and account for approximately 30% of genetic heritability. A vast majority of risk variants reside in non-coding regions, which makes it quite challenging to interpret their functional implications in the SLE-affected immune system, suggesting the importance of understanding cell type-specific epigenetic regulation around SLE GWAS variants. The latest genetic studies have been highly fruitful as several dozens of SLE loci were newly discovered in the last few years and many loci have come to be understood in systemic approaches integrating GWAS signals with other biological resources. In this review, we summarize SLE-associated genetic variants in both the major histocompatibility complex (MHC) and non-MHC loci, examining polygenetic risk scores for SLE and their associations with clinical features. Finally, variant-driven pathogenetic functions underlying genetic associations are described, coupled with discussion about challenges and future directions in genetic studies on SLE.

Integrated single-cell transcriptomics and epigenomics reveals strong germinal center-associated etiology of autoimmune risk loci

The germinal center (GC) response is critical for both effective adaptive immunity and establishing peripheral tolerance by limiting autoreactive B cells. Dysfunction in these processes can lead to defective immune responses to infection or contribute to autoimmune disease. To understand the gene regulatory principles underlying the GC response, we generated a single-cell transcriptomic and epigenomic atlas of the human tonsil, a widely studied and representative lymphoid tissue. We characterize diverse immune cell subsets and build a trajectory of dynamic gene expression and transcription factor activity during B cell activation, GC formation, and plasma cell differentiation. We subsequently leverage cell type–specific transcriptomic and epigenomic maps to interpret potential regulatory impact of genetic variants implicated in autoimmunity, revealing that many exhibit their greatest regulatory potential in GC-associated cellular populations. These included gene loci linked with known roles in GC biology (IL21, IL21R, IL4R, and BCL6) and transcription factors regulating B cell differentiation (POU2AF1 and HHEX). Together, these analyses provide a powerful new cell type–resolved resource for the interpretation of cellular and genetic causes underpinning autoimmune disease.

Clinical relevance of T follicular helper cells in systemic lupus erythematosus

INTRODUCTION

T helper cells regulate a variety of immune responses and are involved in the pathogenesis of infection, allergy and autoimmune diseases. T follicular helper (Tfh) cells, which induce B cell maturation, play an important role in the production of the extremely diverse autoantibodies found in systemic lupus erythematosus (SLE).

AREA COVERED

We provide an overview of the plasticity and diversity of Tfh cells in humans and their involvement in the pathology and pathogenesis of SLE. Our review outlines the potential of Tfh cells as a therapeutic target for SLE.

EXPERT OPINION

Tfh cells are involved in the pathogenesis of SLE based on their plasticity and diversity. Tfh cell differentiation and function are variably regulated by cytokines (IL-12, interferons, IL-2, etc), co-stimulatory molecules (ICOS, CD40L, OX40, etc), and intracellular signals (JAK-STAT, etc). Elucidation of the mechanisms underlying Tfh cell differentiation and function may lead to the development of new therapies for SLE.

Detecting chromosomal interactions in Capture Hi-C data with CHiCAGO and companion tools

Capture Hi-C is widely used to obtain high-resolution profiles of chromosomal interactions involving, at least on one end, regions of interest such as gene promoters. Signal detection in Capture Hi-C data is challenging and cannot be adequately accomplished with tools developed for other chromosome conformation capture methods, including standard Hi-C. Capture Hi-C Analysis of Genomic Organization (CHiCAGO) is a computational pipeline developed specifically for Capture Hi-C analysis. It implements a statistical model accounting for biological and technical background components, as well as bespoke normalization and multiple testing procedures for this data type. Here we provide a step-by-step guide to the CHiCAGO workflow that is aimed at users with basic experience of the command line and R. We also describe more advanced strategies for tuning the key parameters for custom experiments and provide guidance on data preprocessing and downstream analysis using companion tools. In a typical experiment, CHiCAGO takes ~2-3 h to run, although pre- and postprocessing steps may take much longer.

Identification of 22 susceptibility loci associated with testicular germ cell tumors

Testicular germ cell tumors (TGCT) are the most common tumor in young white men and have a high heritability. In this study, the international Testicular Cancer Consortium assemble 10,156 and 179,683 men with and without TGCT, respectively, for a genome-wide association study. This meta-analysis identifies 22 TGCT susceptibility loci, bringing the total to 78, which account for 44% of disease heritability. Men with a polygenic risk score (PRS) in the 95th percentile have a 6.8-fold increased risk of TGCT compared to men with median scores. Among men with independent TGCT risk factors such as cryptorchidism, the PRS may guide screening decisions with the goal of reducing treatment-related complications causing long-term morbidity in survivors. These findings emphasize the interconnected nature of two known pathways that promote TGCT susceptibility: male germ cell development within its somatic niche and regulation of chromosomal division and structure, and implicate an additional biological pathway, mRNA translation.

Strong association of common variants in the miRNA-binding site of NOD2 gene with clinicopathological characteristics and disease activity of systemic lupus erythematosus

INTRODUCTION/OBJECTIVES

Systemic lupus erythematosus (SLE) is a multifactorial systemic autoimmune disease, in which genetic susceptibility plays a pivotal role. The nucleotide oligomerization domain 2 (NOD2) gene is one of the main regulators of chronic inflammatory conditions and could be involved in SLE pathogenesis. Single nucleotide polymorphisms (SNPs) in miRNA binding sites which are located in 3'UTR of the NOD2 gene could be associated with SLE risk by dysregulation of NOD2 expression. In the present study, we assessed the possible association between SNPs rs3135500 and rs3135499 in the NOD2 gene with SLE risk in the Iranian population.

METHODS

A case-control study using 110 SLE patients and 120 control subjects was undertaken to estimate rs3135500 (G > A) and rs3135499 (A > C) genotypes via real-time PCR high-resolution melting method (HRM).

RESULTS

No significant association was observed between allele and genotype frequencies of rs3135500 and rs3135499 polymorphisms and SLE risk in this population (P > 0.05). However, there was an obvious association between rs3135500 (A allele) with laboratory factors that are associated with disease activity (P < 0.05) and some clinical manifestations that are associated with disease severity such as neurological symptoms, skin manifestations, renal involvements, and higher serum concentration of creatinine (P < 0.05). Besides, rs3135499 (C allele) was correlated with renal involvement and also the concentration of creatinine (P < 0.05). Moreover, in the patients group, the risk alleles in these polymorphisms were associated with lower age of onset (P < 0.05).

CONCLUSIONS

Our results suggest a substantial association between NOD2 polymorphisms with clinicopathological characteristics and SLE disease activity. Key Points • Single nucleotide polymorphisms (SNPs) in miRNA binding sites which are located in 3'UTR of the NOD2 gene could be associated with SLE risk by dysregulation of NOD2 expression. • Our results suggested that two miRSNPs (rs3135500 and rs3135499) in the NOD2 gene were meaningfully correlated with clinicopathological characteristics and disease activity of SLE.

T Cells in Systemic Lupus Erythematosus

T-cell dysregulation has been implicated in the loss of tolerance and overactivation of B cells in systemic lupus erythematosus (SLE). Recent studies have identified T-cell subsets and genetic, epigenetic, and environmental factors that contribute to pathogenic T-cell differentiation, as well as disease pathogenesis and clinical phenotypes in SLE. Many therapeutics targeting T-cell pathways are under development, and although many have not progressed in clinical trials, the recent approval of the calcineurin inhibitor voclosporin is encouraging. Further study of T-cell subsets and biomarkers of T-cell action may pave the way for specific targeting of pathogenic T-cell populations in SLE.

3D promoter architecture re-organization during iPSC-derived neuronal cell differentiation implicates target genes for neurodevelopmental disorders

Neurodevelopmental disorders are thought to arise from interrupted development of the brain at an early age. Genome-wide association studies (GWAS) have identified hundreds of loci associated with susceptibility to neurodevelopmental disorders; however, which noncoding variants regulate which genes at these loci is often unclear. To implicate neuronal GWAS effector genes, we performed an integrated analysis of transcriptomics, epigenomics and chromatin conformation changes during the development from Induced pluripotent stem cell-derived neuronal progenitor cells (NPCs) into neurons using a combination of high-resolution promoter-focused Capture-C, ATAC-seq and RNA-seq. We observed that gene expression changes during the NPC-to-neuron transition were highly dependent on both promoter accessibility changes and long-range interactions which connect distal cis-regulatory elements (enhancer or silencers) to developmental-stage-specific genes. These genome-scale promoter-cis-regulatory-element atlases implicated 454 neurodevelopmental disorder-associated, putative causal variants mapping to 600 distal targets. These putative effector genes were significantly enriched for pathways involved in the regulation of neuronal development and chromatin organization, with 27 % expressed in a stage-specific manner. The intersection of open chromatin and chromatin conformation revealed development-stage-specific gene regulatory architectures during neuronal differentiation, providing a rich resource to aid characterization of the genetic and developmental basis of neurodevelopmental disorders.

Variant to Gene Mapping to Discover New Targets for Immune Tolerance

The breakdown of immunological tolerance leads to autoimmune disease, and the mechanisms that maintain self-tolerance, especially in humans, are not fully understood. Genome-wide association studies (GWAS) have identified hundreds of human genetic loci statistically linked to autoimmune disease risk, and epigenetic modifications of DNA and chromatin at these loci have been associated with autoimmune disease risk. Because the vast majority of these signals are located far from genes, identifying causal variants, and their functional consequences on the correct effector genes, has been challenging. These limitations have hampered the translation of GWAS findings into novel drug targets and clinical interventions, but recent advances in understanding the spatial organization of the genome in the nucleus have offered mechanistic insights into gene regulation and answers to questions left open by GWAS. Here we discuss the potential for 'variant-to-gene mapping' approaches that integrate GWAS with 3D functional genomic data to identify human genes involved in the maintenance of tolerance.

Multi-cell type gene coexpression network analysis reveals coordinated interferon response and cross-cell type correlations in systemic lupus erythematosus

Systemic lupus erythematosus (SLE) is an incurable autoimmune disease disproportionately affecting women. A major obstacle in finding targeted therapies for SLE is its remarkable heterogeneity in clinical manifestations as well as in the involvement of distinct cell types. To identify cell-specific targets as well as cross-correlation relationships among expression programs of different cell types, we here analyze six major circulating immune cell types from SLE patient blood. Our results show that presence of an interferon response signature stratifies patients into two distinct groups (IFNneg vs. IFNpos). Comparing these two groups using differential gene expression and differential gene coexpression analysis, we prioritize a relatively small list of genes from classical monocytes including two known immune modulators: TNFSF13B/BAFF (target of belimumab, an approved therapeutic for SLE) and IL1RN (the basis of anakinra, a therapeutic for rheumatoid arthritis). We then develop a multi-cell type extension of the weighted gene coexpression network analysis (WGCNA) framework, termed mWGCNA. Applying mWGCNA to RNA-seq data from six sorted immune cell populations (15 SLE, 10 healthy donors), we identify a coexpression module with interferon-stimulated genes (ISGs) among all cell types and a cross-cell type correlation linking expression of specific T helper cell markers to B cell response as well as to TNFSF13B expression from myeloid cells, all of which in turn correlates with disease severity of IFNpos patients. Our results demonstrate the power of a hypothesis-free and data-driven approach to discover drug targets and to reveal novel cross-correlation across cell types in SLE with implications for other autoimmune diseases.

Nonmechanical parfocal and autofocus features based on wave propagation distribution in lensfree holographic microscopy

Performing long-term cell observations is a non-trivial task for conventional optical microscopy, since it is usually not compatible with environments of an incubator and its temperature and humidity requirements. Lensless holographic microscopy, being entirely based on semiconductor chips without lenses and without any moving parts, has proven to be a very interesting alternative to conventional microscopy. Here, we report on the integration of a computational parfocal feature, which operates based on wave propagation distribution analysis, to perform a fast autofocusing process. This unique non-mechanical focusing approach was implemented to keep the imaged object staying in-focus during continuous long-term and real-time recordings. A light-emitting diode (LED) combined with pinhole setup was used to realize a point light source, leading to a resolution down to 2.76 μm. Our approach delivers not only in-focus sharp images of dynamic cells, but also three-dimensional (3D) information on their (x, y, z)-positions. System reliability tests were conducted inside a sealed incubator to monitor cultures of three different biological living cells (i.e., MIN6, neuroblastoma (SH-SY5Y), and Prorocentrum minimum). Altogether, this autofocusing framework enables new opportunities for highly integrated microscopic imaging and dynamic tracking of moving objects in harsh environments with large sample areas.

Roles of follicular helper and regulatory T cells in allergic diseases and allergen immunotherapy

Allergic diseases are characterized by overactive type 2 immune responses to allergens and immunoglobulin E (IgE)-mediated hypersensitivity. Emerging evidence suggests that follicular helper T (T_FH ) cells, rather than type 2 T-helper (T_H 2) cells, play a crucial role in controlling IgE production. However, follicular regulatory T (T_FR ) cells, a specialized subset of regulatory T (T_REG ) cells resident in B-cell follicles, restricts T_FH cell-mediated help in extrafollicular antibody production, germinal center (GC) formation, immunoglobulin affinity maturation, and long-lived, high-affinity plasma and memory B-cell differentiation. In mouse models of allergic asthma and food allergy, CXCR5⁺ T_FH cells, not CXCR5^- conventional T_H 2 cells, are needed to support IgE production, otherwise exacerbated by CXCR5⁺ T_FR cell deletion. Upregulation of T_FH cell activities, including a skewing toward type 2 T_FH (T_FH 2) and IL-13 producing T_FH (T_FH 13) phenotypes, and defects in T_FR cells have been identified in patients with allergic diseases. Allergen immunotherapy (AIT) reinstates the balance between T_FH and T_FR cells in patients with allergic diseases, resulting in clinical benefits. Collectively, further understanding of T_FH and T_FR cells and their role in the immunopathogenesis of allergic diseases creates opportunities to develop novel therapeutic approaches.

Biological constraints on GWAS SNPs at suggestive significance thresholds reveal additional BMI loci

To uncover novel significant association signals (p<5×10-8), genome-wide association studies (GWAS) requires increasingly larger sample sizes to overcome statistical correction for multiple testing. As an alternative, we aimed to identify associations among suggestive signals (5 × 10-8≤p<5×10-4) in increasingly powered GWAS efforts using chromatin accessibility and direct contact with gene promoters as biological constraints. We conducted retrospective analyses of three GIANT BMI GWAS efforts using ATAC-seq and promoter-focused Capture C data from human adipocytes and embryonic stem cell (ESC)-derived hypothalamic-like neurons. This approach, with its extremely low false-positive rate, identified 15 loci at p<5×10-5 in the 2010 GWAS, of which 13 achieved genome-wide significance by 2018, including at NAV1, MTIF3, and ADCY3. Eighty percent of constrained 2015 loci achieved genome-wide significance in 2018. We observed similar results in waist-to-hip ratio analyses. In conclusion, biological constraints on sub-significant GWAS signals can reveal potentially true-positive loci for further investigation in existing data sets without increasing sample size.

Next steps in the identification of gene targets for type 1 diabetes

The purpose of this review is to provide a view of the future of genomics and other omics approaches in defining the genetic contribution to all stages of risk of type 1 diabetes and the functional impact and clinical implementations of the associated variants. From the recognition nearly 50 years ago that genetics (in the form of HLA) distinguishes risk of type 1 diabetes from type 2 diabetes, advances in technology and sample acquisition through collaboration have identified over 60 loci harbouring SNPs associated with type 1 diabetes risk. Coupled with HLA region genes, these variants account for the majority of the genetic risk (~50% of the total risk); however, relatively few variants are located in coding regions of genes exerting a predicted protein change. The vast majority of genetic risk in type 1 diabetes appears to be attributed to regions of the genome involved in gene regulation, but the target effectors of those genetic variants are not readily identifiable. Although past genetic studies clearly implicated immune-relevant cell types involved in risk, the target organ (the beta cell) was left untouched. Through emergent technologies, using combinations of genetics, gene expression, epigenetics, chromosome conformation and gene editing, novel landscapes of how SNPs regulate genes have emerged. Furthermore, both the immune system and the beta cell and their biological pathways have been implicated in a context-specific manner. The use of variants from immune and beta cell studies distinguish type 1 diabetes from type 2 diabetes and, when they are combined in a genetic risk score, open new avenues for prediction and treatment. Graphical abstract.