Dissection of a Complex Disease Susceptibility Region Using a Bayesian Stochastic Search Approach to Fine Mapping

Dysfunctional β-cell longevity in diabetes relies on energy conservation and positive epistasis

Long-lived PFKFB3-expressing β-cells are dysfunctional partly because of prevailing glycolysis that compromises metabolic coupling of insulin secretion. Their accumulation in type 2 diabetes (T2D) appears to be related to the loss of apoptotic competency of cell fitness competition that maintains islet function by favoring constant selection of healthy "winner" cells. To investigate how PFKFB3 can disguise the competitive traits of dysfunctional "loser" β-cells, we analyzed the overlap between human β-cells with bona fide "loser signature" across diabetes pathologies using the HPAP scRNA-seq and spatial transcriptomics of PFKFB3-positive β-cells from nPOD T2D pancreata. The overlapping transcriptional profile of "loser" β-cells was represented by down-regulated ribosomal biosynthesis and genes encoding for mitochondrial respiration. PFKFB3-positive "loser" β-cells had the reduced expression of HLA class I and II genes. Gene-gene interaction analysis revealed that PFKFB3 rs1983890 can interact with the anti-apoptotic gene MAIP1 implicating positive epistasis as a mechanism for prolonged survival of "loser" β-cells in T2D. Inhibition of PFKFB3 resulted in the clearance of dysfunctional "loser" β-cells leading to restored glucose tolerance in the mouse model of T2D.

The correlation between CpG methylation and gene expression is driven by sequence variants

Gene promoter and enhancer sequences are bound by transcription factors and are depleted of methylated CpG sites (cytosines preceding guanines in DNA). The absence of methylated CpGs in these sequences typically correlates with increased gene expression, indicating a regulatory role for methylation. We used nanopore sequencing to determine haplotype-specific methylation rates of 15.3 million CpG units in 7,179 whole-blood genomes. We identified 189,178 methylation depleted sequences where three or more proximal CpGs were unmethylated on at least one haplotype. A total of 77,789 methylation depleted sequences (~41%) associated with 80,503 cis-acting sequence variants, which we termed allele-specific methylation quantitative trait loci (ASM-QTLs). RNA sequencing of 896 samples from the same blood draws used to perform nanopore sequencing showed that the ASM-QTL, that is, DNA sequence variability, drives most of the correlation found between gene expression and CpG methylation. ASM-QTLs were enriched 40.2-fold (95% confidence interval 32.2, 49.9) among sequence variants associating with hematological traits, demonstrating that ASM-QTLs are important functional units in the noncoding genome.

Fast and flexible joint fine-mapping of multiple traits via the Sum of Single Effects model

We introduce mvSuSiE, a multi-trait fine-mapping method for identifying putative causal variants from genetic association data (individual-level or summary data). mvSuSiE learns patterns of shared genetic effects from data, and exploits these patterns to improve power to identify causal SNPs. Comparisons on simulated data show that mvSuSiE is competitive in speed, power and precision with existing multi-trait methods, and uniformly improves on single-trait fine-mapping (SuSiE) in each trait separately. We applied mvSuSiE to jointly fine-map 16 blood cell traits using data from the UK Biobank. By jointly analyzing the traits and modeling heterogeneous effect sharing patterns, we discovered a much larger number of causal SNPs (>3,000) compared with single-trait fine-mapping, and with narrower credible sets. mvSuSiE also more comprehensively characterized the ways in which the genetic variants affect one or more blood cell traits; 68% of causal SNPs showed significant effects in more than one blood cell type.

Estimating gene-level false discovery probability improves eQTL statistical fine-mapping precision

Statistical fine-mapping prioritizes putative causal variants from a large number of candidate variants, and is widely used in expression quantitative loci (eQTLs) studies. In eQTL fine-mapping, the existence of causal variants for gene expression is not guaranteed, since the genetic heritability of gene expression explained by nearby (cis-) variants is limited. Here we introduce a refined fine-mapping algorithm, named Knockoff-Finemap combination (KFc). KFc estimates the probability that the causal variant(s) exist in the cis-window of a gene through construction of knockoff genotypes (i.e. a set of synthetic genotypes that resembles the original genotypes), and uses it to adjust the posterior inclusion probabilities (PIPs). Utilizing simulated gene expression data, we show that KFc results in calibrated PIP distribution with improved precision. When applied to gene expression data of 465 genotyped samples from the Japan COVID-19 Task Force (JCTF), KFc resulted in significant enrichment of a functional score as well as reporter assay hits in the top PIP bins. When combined with functional priors derived from an external fine-mapping study (GTEx), KFc resulted in a significantly higher proportion of hematopoietic trait putative causal variants in the top PIP bins. Our work presents improvements in the precision of a major fine-mapping algorithm.

Translating non-coding genetic associations into a better understanding of immune-mediated disease

Genome-wide association studies have identified hundreds of genetic loci that are associated with immune-mediated diseases. Most disease-associated variants are non-coding, and a large proportion of these variants lie within enhancers. As a result, there is a pressing need to understand how common genetic variation might affect enhancer function and thereby contribute to immune-mediated (and other) diseases. In this Review, we first describe statistical and experimental methods to identify causal genetic variants that modulate gene expression, including statistical fine-mapping and massively parallel reporter assays. We then discuss approaches to characterise the mechanisms by which these variants modulate immune function, such as clustered regularly interspaced short palindromic repeats (CRISPR)-based screens. We highlight examples of studies that, by elucidating the effects of disease variants within enhancers, have provided important insights into immune function and uncovered key pathways of disease.

Fine-mapping from summary data with the "Sum of Single Effects" model

In recent work, Wang et al introduced the "Sum of Single Effects" (SuSiE) model, and showed that it provides a simple and efficient approach to fine-mapping genetic variants from individual-level data. Here we present new methods for fitting the SuSiE model to summary data, for example to single-SNP z-scores from an association study and linkage disequilibrium (LD) values estimated from a suitable reference panel. To develop these new methods, we first describe a simple, generic strategy for extending any individual-level data method to deal with summary data. The key idea is to replace the usual regression likelihood with an analogous likelihood based on summary data. We show that existing fine-mapping methods such as FINEMAP and CAVIAR also (implicitly) use this strategy, but in different ways, and so this provides a common framework for understanding different methods for fine-mapping. We investigate other common practical issues in fine-mapping with summary data, including problems caused by inconsistencies between the z-scores and LD estimates, and we develop diagnostics to identify these inconsistencies. We also present a new refinement procedure that improves model fits in some data sets, and hence improves overall reliability of the SuSiE fine-mapping results. Detailed evaluations of fine-mapping methods in a range of simulated data sets show that SuSiE applied to summary data is competitive, in both speed and accuracy, with the best available fine-mapping methods for summary data.

Methods for statistical fine-mapping and their applications to auto-immune diseases

Although genome-wide association studies (GWAS) have identified thousands of loci in the human genome that are associated with different traits, understanding the biological mechanisms underlying the association signals identified in GWAS remains challenging. Statistical fine-mapping is a method aiming to refine GWAS signals by evaluating which variant(s) are truly causal to the phenotype. Here, we review the types of statistical fine-mapping methods that have been widely used to date, with a focus on recently developed functionally informed fine-mapping (FIFM) methods that utilize functional annotations. We then systematically review the applications of statistical fine-mapping in autoimmune disease studies to highlight the value of statistical fine-mapping in biological contexts.

Mutation of Hashimoto’s Thyroiditis and Papillary Thyroid Carcinoma Related Genes and the Screening of Candidate Genes

Clinical studies have shown similarities in the genetic background and biological functional characteristics between Hashimoto’s thyroiditis (HT) and papillary thyroid carcinoma (PTC), and that HT may increase risks of PTC. Here, we set to determine the gene expression specificity of HT and PTC by screening related genes or co-expressed genes and exploring their genetic correlation. Referencing the Oncomine database, HT-related genes were discovered to be expressed in many different types of thyroid cancer, such as TSHR that is highly expressed in thyroid cancer. An in-depth genetic analysis and verification of 35 cancer and paracancerous tissue pairs from patients with thyroid cancer, and 35 tissues and blood cells pairs from patients with Hashimoto’s thyroiditis was conducted. Gene chip technology research showed that TSHR, BACH2, FOXE1, RNASET2, CTLA4, PTPN22, IL2RA and other HT-related genes were all expressed in PTC, in which TSHR was significantly over-expressed in PTC patients sensitive to radioactive iodine therapy, while BACH2 was significantly under-expressed in these patients. The biologically significant candidate Tag SNP highlighted from HT-related genes was screened by the high-throughput detection method. Somatic mutations in patients with HT and PTC were detected by target region capture technique, and 75 mutations were found in patients with HT and PTC. The upstream regulatory factors of the different genes shared by HT and PTC were analyzed based on Ingenuity Pathway Analysis (IPA), and it was found that HIF-1α and PD-L1 could be used as important upstream regulatory signal molecules. These results provide a basis for screening key diagnostic genes of PTC by highlighting the relationship between some HT-related genes and their polymorphisms in the pathogenesis of PTC.

The flashfm approach for fine-mapping multiple quantitative traits

Joint fine-mapping that leverages information between quantitative traits could improve accuracy and resolution over single-trait fine-mapping. Using summary statistics, flashfm (flexible and shared information fine-mapping) fine-maps signals for multiple traits, allowing for missing trait measurements and use of related individuals. In a Bayesian framework, prior model probabilities are formulated to favour model combinations that share causal variants to capitalise on information between traits. Simulation studies demonstrate that both approaches produce broadly equivalent results when traits have no shared causal variants. When traits share at least one causal variant, flashfm reduces the number of potential causal variants by 30% compared with single-trait fine-mapping. In a Ugandan cohort with 33 cardiometabolic traits, flashfm gave a 20% reduction in the total number of potential causal variants from single-trait fine-mapping. Here we show flashfm is computationally efficient and can easily be deployed across publicly available summary statistics for signals in up to six traits.

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

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

Fine-mapping, trans-ancestral and genomic analyses identify causal variants, cells, genes and drug targets for type 1 diabetes

We report the largest and most diverse genetic study of type 1 diabetes (T1D) to date (61,427 participants), yielding 78 genome-wide-significant (P < 5 × 10-8) regions, including 36 that are new. We define credible sets of T1D-associated variants and show that they are enriched in immune-cell accessible chromatin, particularly CD4+ effector T cells. Using chromatin-accessibility profiling of CD4+ T cells from 115 individuals, we map chromatin-accessibility quantitative trait loci and identify five regions where T1D risk variants co-localize with chromatin-accessibility quantitative trait loci. We highlight rs72928038 in BACH2 as a candidate causal T1D variant leading to decreased enhancer accessibility and BACH2 expression in T cells. Finally, we prioritize potential drug targets by integrating genetic evidence, functional genomic maps and immune protein-protein interactions, identifying 12 genes implicated in T1D that have been targeted in clinical trials for autoimmune diseases. These findings provide an expanded genomic landscape for T1D.

Genetic and environmental factors regulate the type 1 diabetes gene CTSH via differential DNA methylation

Cathepsin H (CTSH) is a type 1 diabetes (T1D) risk gene; large-scale genetic and epidemiological studies found that T1D genetic risk correlates with high CTSH expression, rapid decline of beta-cell function, and early onset T1D. Counterintuitively, transcriptional downregulation of CTSH by proinflammatory cytokines has been shown to promote beta-cell apoptosis. Here, we potentially explain these observed contrasting effects, describing a new mechanism where proinflammatory cytokines and T1D genetic risk variants regulate CTSH transcription via differential DNA methylation. We show that, in human islets, CTSH downregulation by the proinflammatory cytokine cocktail interleukin 1β + tumor necrosis factor α + interferon γ was coupled with DNA hypermethylation in an open chromatin region in CTSH intron 1. A luciferase assay in human embryonic kidney 293 cells revealed that methylation of three key cytosine-phosphate-guanine dinucleotide (CpG) residues in intron 1 was responsible for the reduction of promoter activity. We further found that cytokine-induced intron 1 hypermethylation is caused by lowered Tet1/3 activities, suggesting that attenuated active demethylation lowered CTSH transcription. Importantly, individuals who carry the T1D risk variant showed lower methylation variability at the intron 1 CpG residues, presumably making them less sensitive to cytokines, whereas individuals who carry the protective variant showed higher methylation variability, presumably making them more sensitive to cytokines and implying differential responses to environment between the two patient populations. These findings suggest that genetic and environmental influences on a T1D locus are mediated by differential variability and mean of DNA methylation.

Fine-mapping genetic associations

Whilst thousands of genetic variants have been associated with human traits, identifying the subset of those variants that are causal requires a further 'fine-mapping' step. We review the basic fine-mapping approach, which is computationally fast and requires only summary data, but depends on an assumption of a single causal variant per associated region which is recognized as biologically unrealistic. We discuss different ways that the approach has been built upon to accommodate multiple causal variants in a region and to incorporate additional layers of functional annotation data. We further review methods for simultaneous fine-mapping of multiple datasets, either exploiting different linkage disequilibrium (LD) structures across ancestries or borrowing information between distinct but related traits. Finally, we look to the future and the opportunities that will be offered by increasingly accurate maps of causal variants for a multitude of human traits.

Polygenic inheritance, GWAS, polygenic risk scores, and the search for functional variants

The reconciliation between Mendelian inheritance of discrete traits and the genetically based correlation between relatives for quantitative traits was Fisher's infinitesimal model of a large number of genetic variants, each with very small effects, whose causal effects could not be individually identified. The development of genome-wide genetic association studies (GWAS) raised the hope that it would be possible to identify single polymorphic variants with identifiable functional effects on complex traits. It soon became clear that, with larger and larger GWAS on more and more complex traits, most of the significant associations had such small effects, that identifying their individual functional effects was essentially hopeless. Polygenic risk scores that provide an overall estimate of the genetic propensity to a trait at the individual level have been developed using GWAS data. These provide useful identification of groups of individuals with substantially increased risks, which can lead to recommendations of medical treatments or behavioral modifications to reduce risks. However, each such claim will require extensive investigation to justify its practical application. The challenge now is to use limited genetic association studies to find individually identifiable variants of significant functional effect that can help to understand the molecular basis of complex diseases and traits, and so lead to improved disease prevention and treatment. This can best be achieved by 1) the study of rare variants, often chosen by careful candidate assessment, and 2) the careful choice of phenotypes, often extremes of a quantitative variable, or traits with relatively high heritability.

In The Blood: Connecting Variant to Function In Human Hematopoiesis

Genome-wide association studies (GWAS) have identified thousands of genetic variants associated with a range of human diseases and traits. However, understanding the mechanisms by which these genetic variants have an impact on associated diseases and traits, often referred to as the variant-to-function (V2F) problem, remains a significant hurdle. Solving the V2F challenge requires us to identify causative genetic variants, relevant cell types/states, target genes, and mechanisms by which variants can cause diseases or alter phenotypic traits. We discuss emerging functional approaches that are being applied to tackle the V2F problem for blood cell traits, illuminating how human genetic variation can impact on key mechanisms in hematopoiesis, as well as highlighting future prospects for this nascent field.

simGWAS: a fast method for simulation of large scale case-control GWAS summary statistics

Motivation

Methods for analysis of GWAS summary statistics have encouraged data sharing and democratized the analysis of different diseases. Ideal validation for such methods is application to simulated data, where some ‘truth’ is known. As GWAS increase in size, so does the computational complexity of such evaluations; standard practice repeatedly simulates and analyses genotype data for all individuals in an example study.

Results

We have developed a novel method based on an alternative approach, directly simulating GWAS summary data, without individual data as an intermediate step. We mathematically derive the expected statistics for any set of causal variants and their effect sizes, conditional upon control haplotype frequencies (available from public reference datasets). Simulation of GWAS summary output can be conducted independently of sample size by simulating random variates about these expected values. Across a range of scenarios, our method, produces very similar output to that from simulating individual genotypes with a substantial gain in speed even for modest sample sizes. Fast simulation of GWAS summary statistics will enable more complete and rapid evaluation of summary statistic methods as well as opening new potential avenues of research in fine mapping and gene set enrichment analysis.

Availability and implementation

Our method is available under a GPL license as an R package from http://github.com/chr1swallace/simGWAS.

Supplementary information

Supplementary data are available at Bioinformatics online.

Dynamic changes in the regulatory T-cell heterogeneity and function by murine IL-2 mutein

Single-cell RNA-seq analysis reveals that IL-2 mutein treatment expands multiple sub-states of regulatory T cells with superior suppressive function in mice.

The therapeutic expansion of Foxp3⁺ regulatory T cells (Tregs) shows promise for treating autoimmune and inflammatory disorders. Yet, how this treatment affects the heterogeneity and function of Tregs is not clear. Using single-cell RNA-seq analysis, we characterized 31,908 Tregs from the mice treated with a half-life extended mutant form of murine IL-2 (IL-2 mutein, IL-2M) that preferentially expanded Tregs, or mouse IgG Fc as a control. Cell clustering analysis revealed that IL-2M specifically expands multiple sub-states of Tregs with distinct expression profiles. TCR profiling with single-cell analysis uncovered Treg migration across tissues and transcriptional changes between clonally related Tregs after IL-2M treatment. Finally, we identified IL-2M–expanded Tnfrsf9⁺Il1rl1⁺ Tregs with superior suppressive function, highlighting the potential of IL-2M to expand highly suppressive Foxp3⁺ Tregs.

Functional effects of variation in transcription factor binding highlight long-range gene regulation by epromoters

Identifying DNA cis-regulatory modules (CRMs) that control the expression of specific genes is crucial for deciphering the logic of transcriptional control. Natural genetic variation can point to the possible gene regulatory function of specific sequences through their allelic associations with gene expression. However, comprehensive identification of causal regulatory sequences in brute-force association testing without incorporating prior knowledge is challenging due to limited statistical power and effects of linkage disequilibrium. Sequence variants affecting transcription factor (TF) binding at CRMs have a strong potential to influence gene regulatory function, which provides a motivation for prioritizing such variants in association testing. Here, we generate an atlas of CRMs showing predicted allelic variation in TF binding affinity in human lymphoblastoid cell lines and test their association with the expression of their putative target genes inferred from Promoter Capture Hi-C and immediate linear proximity. We reveal >1300 CRM TF-binding variants associated with target gene expression, the majority of them undetected with standard association testing. A large proportion of CRMs showing associations with the expression of genes they contact in 3D localize to the promoter regions of other genes, supporting the notion of 'epromoters': dual-action CRMs with promoter and distal enhancer activity.

Improving the coverage of credible sets in Bayesian genetic fine-mapping

Genome Wide Association Studies (GWAS) have successfully identified thousands of loci associated with human diseases. Bayesian genetic fine-mapping studies aim to identify the specific causal variants within GWAS loci responsible for each association, reporting credible sets of plausible causal variants, which are interpreted as containing the causal variant with some "coverage probability". Here, we use simulations to demonstrate that the coverage probabilities are over-conservative in most fine-mapping situations. We show that this is because fine-mapping data sets are not randomly selected from amongst all causal variants, but from amongst causal variants with larger effect sizes. We present a method to re-estimate the coverage of credible sets using rapid simulations based on the observed, or estimated, SNP correlation structure, we call this the "adjusted coverage estimate". This is extended to find "adjusted credible sets", which are the smallest set of variants such that their adjusted coverage estimate meets the target coverage. We use our method to improve the resolution of a fine-mapping study of type 1 diabetes. We found that in 27 out of 39 associated genomic regions our method could reduce the number of potentially causal variants to consider for follow-up, and found that none of the 95% or 99% credible sets required the inclusion of more variants-a pattern matched in simulations of well powered GWAS. Crucially, our method requires only GWAS summary statistics and remains accurate when SNP correlations are estimated from a large reference panel. Using our method to improve the resolution of fine-mapping studies will enable more efficient expenditure of resources in the follow-up process of annotating the variants in the credible set to determine the implicated genes and pathways in human diseases.

Reworking GWAS Data to Understand the Role of Nongenetic Factors in MS Etiopathogenesis

Genome-wide association studies have identified more than 200 multiple sclerosis (MS)-associated loci across the human genome over the last decade, suggesting complexity in the disease etiology. This complexity poses at least two challenges: the definition of an etiological model including the impact of nongenetic factors, and the clinical translation of genomic data that may be drivers for new druggable targets. We reviewed studies dealing with single genes of interest, to understand how MS-associated single nucleotide polymorphism (SNP) variants affect the expression and the function of those genes. We then surveyed studies on the bioinformatic reworking of genome-wide association studies (GWAS) data, with aggregate analyses of many GWAS loci, each contributing with a small effect to the overall disease predisposition. These investigations uncovered new information, especially when combined with nongenetic factors having possible roles in the disease etiology. In this context, the interactome approach, defined as “modules of genes whose products are known to physically interact with environmental or human factors with plausible relevance for MS pathogenesis”, will be reported in detail. For a future perspective, a polygenic risk score, defined as a cumulative risk derived from aggregating the contributions of many DNA variants associated with a complex trait, may be integrated with data on environmental factors affecting the disease risk or protection.

Genetic Variants Predisposing Most Strongly to Type 1 Diabetes Diagnosed Under Age 7 Years Lie Near Candidate Genes That Function in the Immune System and in Pancreatic β-Cells

OBJECTIVE

Immunohistological analyses of pancreata from patients with type 1 diabetes suggest distinct autoimmune islet β-cell pathology between those diagnosed at <7 years (<7 group) and those diagnosed at age ≥13 years (≥13 group), with both B- and T-lymphocyte islet inflammation common in children in the <7 group, whereas B cells are rare in the ≥13 group. Based on these observations, we sought to identify differences in genetic susceptibility between these prespecified age-at-diagnosis groups to inform on the etiology of the most aggressive form of type 1 diabetes that initiates in the first years of life.

RESEARCH DESIGN AND METHODS

Using multinomial logistic regression models, we tested if known type 1 diabetes loci (17 within the HLA and 55 non-HLA loci) had significantly stronger effect sizes in the <7 group compared with the ≥13 group, using genotype data from 27,071 individuals (18,485 control subjects and 3,121 case subjects diagnosed at <7 years, 3,757 at 7-13 years, and 1,708 at ≥13 years).

RESULTS

Six HLA haplotypes/classical alleles and six non-HLA regions, one of which functions specifically in β-cells (GLIS3) and the other five likely affecting key T-cell (IL2RA, IL10, IKZF3, and THEMIS), thymus (THEMIS), and B-cell development/functions (IKZF3 and IL10) or in both immune and β-cells (CTSH), showed evidence for stronger effects in the <7 group.

CONCLUSIONS

A subset of type 1 diabetes-associated variants are more prevalent in children diagnosed under the age of 7 years and are near candidate genes that act in both pancreatic β- and immune cells.

The Genetic Contribution to Type 1 Diabetes

PURPOSE OF REVIEW

To provide an updated summary of discoveries made to date resulting from genome-wide association study (GWAS) and sequencing studies, and to discuss the latest loci added to the growing repertoire of genetic signals predisposing to type 1 diabetes (T1D).

RECENT FINDINGS

Genetic studies have identified over 60 loci associated with T1D susceptibility. GWAS alone does not specifically inform on underlying mechanisms, but in combination with other sequencing and omics-data, advances are being made in our understanding of T1D genetic etiology and pathogenesis. Current knowledge indicates that genetic variation operating in both pancreatic β cells and in immune cells is central in mediating T1D risk. One of the main challenges is to determine how these recently discovered GWAS-implicated variants affect the expression and function of gene products. Once we understand the mechanism of action for disease-causing variants, we will be well placed to apply targeted genomic approaches to impede the premature activation of the immune system in an effort to ultimately prevent the onset of T1D.

A flexible and parallelizable approach to genome-wide polygenic risk scores

The heritability of most complex traits is driven by variants throughout the genome. Consequently, polygenic risk scores, which combine information on multiple variants genome-wide, have demonstrated improved accuracy in genetic risk prediction. We present a new two-step approach to constructing genome-wide polygenic risk scores from meta-GWAS summary statistics. Local linkage disequilibrium (LD) is adjusted for in Step 1, followed by, uniquely, long-range LD in Step 2. Our algorithm is highly parallelizable since block-wise analyses in Step 1 can be distributed across a high-performance computing cluster, and flexible, since sparsity and heritability are estimated within each block. Inference is obtained through a formal Bayesian variable selection framework, meaning final risk predictions are averaged over competing models. We compared our method to two alternative approaches: LDPred and lassosum using all seven traits in the Welcome Trust Case Control Consortium as well as meta-GWAS summaries for type 1 diabetes (T1D), coronary artery disease, and schizophrenia. Performance was generally similar across methods, although our framework provided more accurate predictions for T1D, for which there are multiple heterogeneous signals in regions of both short- and long-range LD. With sufficient compute resources, our method also allows the fastest runtimes.

Stochastic search and joint fine-mapping increases accuracy and identifies previously unreported associations in immune-mediated diseases

Thousands of genetic variants are associated with human disease risk, but linkage disequilibrium (LD) hinders fine-mapping the causal variants. Both lack of power, and joint tagging of two or more distinct causal variants by a single non-causal SNP, lead to inaccuracies in fine-mapping, with stochastic search more robust than stepwise. We develop a computationally efficient multinomial fine-mapping (MFM) approach that borrows information between diseases in a Bayesian framework. We show that MFM has greater accuracy than single disease analysis when shared causal variants exist, and negligible loss of precision otherwise. MFM analysis of six immune-mediated diseases reveals causal variants undetected in individual disease analysis, including in IL2RA where we confirm functional effects of multiple causal variants using allele-specific expression in sorted CD4+ T cells from genotype-selected individuals. MFM has the potential to increase fine-mapping resolution in related diseases enabling the identification of associated cellular and molecular phenotypes.

Amino acid residues in five separate HLA genes can explain most of the known associations between the MHC and primary biliary cholangitis

Primary Biliary Cholangitis (PBC) is a chronic autoimmune liver disease characterised by progressive destruction of intrahepatic bile ducts. The strongest genetic association is with HLA-DQA1*04:01, but at least three additional independent HLA haplotypes contribute to susceptibility. We used dense single nucleotide polymorphism (SNP) data in 2861 PBC cases and 8514 controls to impute classical HLA alleles and amino acid polymorphisms using state-of-the-art methodologies. We then demonstrated through stepwise regression that association in the HLA region can be largely explained by variation at five separate amino acid positions. Three-dimensional modelling of protein structures and calculation of electrostatic potentials for the implicated HLA alleles/amino acid substitutions demonstrated a correlation between the electrostatic potential of pocket P6 in HLA-DP molecules and the HLA-DPB1 alleles/amino acid substitutions conferring PBC susceptibility/protection, highlighting potential new avenues for future functional investigation.

A long-lived IL-2 mutein that selectively activates and expands regulatory T cells as a therapy for autoimmune disease

Susceptibility to multiple autoimmune diseases is associated with common gene polymorphisms influencing IL-2 signaling and T_reg function, making T_reg-specific expansion by IL-2 a compelling therapeutic approach to treatment. As an in vivo IL-2 half-life enhancer we used a non-targeted, effector-function-silent human IgG1 as a fusion protein. An IL-2 mutein (N88D) with reduced binding to the intermediate affinity IL-2Rβγ receptor was engineered with a stoichiometry of two IL-2N88D molecules per IgG, i.e. IgG-(IL-2N88D)₂. The reduced affinity of IgG-(IL-2N88D)₂ for the IL-2Rβγ receptor resulted in a T_reg-selective molecule in human whole blood pSTAT5 assays. Treatment of cynomolgus monkeys with single low doses of IgG-(IL-2N88D)₂ induced sustained preferential activation of T_regs accompanied by a corresponding 10–14-fold increase in CD4⁺ and CD8⁺ CD25⁺FOXP3⁺ T_regs; conditions that had no effect on CD4⁺ or CD8⁺ memory effector T cells. The expanded cynomolgus T_regs had demethylated FOXP3 and CTLA4 epigenetic signatures characteristic of functionally suppressive cells. Humanized mice had similar selective in vivo responses; IgG-(IL-2N88D)₂ increased T_regs while wild-type IgG-IL-2 increased NK cells in addition to T_regs. The expanded human T_regs had demethylated FOXP3 and CTLA4 signatures and were immunosuppressive. These results describe a next-generation immunotherapy using a long-lived and T_reg-selective IL-2 that activates and expands functional T_regsin vivo. Patients should benefit from restored immune homeostasis in a personalized fashion to the extent that their autoimmune disease condition dictates opening up the possibility for remissions and cures.

•

A human IL-2 molecule mutated to decrease binding to the intermediate affinity IL-2 receptor preferentially activates Tregs.

•

Two IL-2 muteins fused to human IgG1 allow for sustained, preferential expansion of Tregs in cynomolgus and humanized mice.

•

As compared to the wild type IL-2 fusion protein, humanized mice expand fewer NK cells in response to the mutein.

•

The dynamic range of Treg increase based on dose suggests the ability to individualize dosing for particular diseases.

From genome-wide associations to candidate causal variants by statistical fine-mapping

Advancing from statistical associations of complex traits with genetic markers to understanding the functional genetic variants that influence traits is often a complex process. Fine-mapping can select and prioritize genetic variants for further study, yet the multitude of analytical strategies and study designs makes it challenging to choose an optimal approach. We review the strengths and weaknesses of different fine-mapping approaches, emphasizing the main factors that affect performance. Topics include interpreting results from genome-wide association studies (GWAS), the role of linkage disequilibrium, statistical fine-mapping approaches, trans-ethnic studies, genomic annotation and data integration, and other analysis and design issues.

Unravelling the Roles of Susceptibility Loci for Autoimmune Diseases in the Post-GWAS Era

Although genome-wide association studies (GWAS) have identified several hundred loci associated with autoimmune diseases, their mechanistic insights are still poorly understood. The human genome is more complex than single nucleotide polymorphisms (SNPs) that are interrogated by GWAS arrays. Apart from SNPs, it also comprises genetic variations such as insertions-deletions, copy number variations, and somatic mosaicism. Although previous studies suggest that common copy number variations do not play a major role in autoimmune disease risk, it is possible that certain rare genetic variations with large effect sizes are relevant to autoimmunity. In addition, other layers of regulations such as gene-gene interactions, epigenetic-determinants, gene and environmental interactions also contribute to the heritability of autoimmune diseases. This review focuses on discussing why studying these elements may allow us to gain a more comprehensive understanding of the aetiology of complex autoimmune traits.

From lipid locus to drug target through human genomics

In the last decade, over 175 genetic loci have robustly been associated to levels of major circulating blood lipids. Most loci are specific to one or two lipids, whereas some (SUGP1, ZPR1, TRIB1, HERPUD1, and FADS1) are associated to all. While exposing the polygenic architecture of circulating lipids and the underpinnings of dyslipidaemia, these genome-wide association studies (GWAS) have provided further evidence of the critical role that lipids play in coronary heart disease (CHD) risk, as indicated by the 2.7-fold enrichment for macrophage gene expression in atherosclerotic plaques and the association of 25 loci (such as PCSK9, APOB, ABCG5-G8, KCNK5, LPL, HMGCR, NPC1L1, CETP, TRIB1, ABO, PMAIP1-MC4R, and LDLR) with CHD. These GWAS also confirmed known and commonly used therapeutic targets, including HMGCR (statins), PCSK9 (antibodies), and NPC1L1 (ezetimibe). As we head into the post-GWAS era, we offer suggestions for how to move forward beyond genetic risk loci, towards refining the biology behind the associations and identifying causal genes and therapeutic targets. Deep phenotyping through lipidomics and metabolomics will refine and increase the resolution to find causal and druggable targets, and studies aimed at demonstrating gene transcriptional and regulatory effects of lipid associated loci will further aid in identifying these targets. Thus, we argue the need for deeply phenotyped, large genetic association studies to reduce costs and failures and increase the efficiency of the drug discovery pipeline. We conjecture that in the next decade a paradigm shift will tip the balance towards a data-driven approach to therapeutic target development and the application of precision medicine where human genomics takes centre stage.

Approaches and advances in the genetic causes of autoimmune disease and their implications

Genome-wide association studies are transformative in revealing the polygenetic basis of common diseases, with autoimmune diseases leading the charge. Although the field is just over 10 years old, advances in understanding the underlying mechanistic pathways of these conditions, which result from a dense multifactorial blend of genetic, developmental and environmental factors, have already been informative, including insights into therapeutic possibilities. Nevertheless, the challenge of identifying the actual causal genes and pathways and their biological effects on altering disease risk remains for many identified susceptibility regions. It is this fundamental knowledge that will underpin the revolution in patient stratification, the discovery of therapeutic targets and clinical trial design in the next 20 years. Here we outline recent advances in analytical and phenotyping approaches and the emergence of large cohorts with standardized gene-expression data and other phenotypic data that are fueling a bounty of discovery and improved understanding of human physiology.

Fine-mapping of prostate cancer susceptibility loci in a large meta-analysis identifies candidate causal variants

Prostate cancer is a polygenic disease with a large heritable component. A number of common, low-penetrance prostate cancer risk loci have been identified through GWAS. Here we apply the Bayesian multivariate variable selection algorithm JAM to fine-map 84 prostate cancer susceptibility loci, using summary data from a large European ancestry meta-analysis. We observe evidence for multiple independent signals at 12 regions and 99 risk signals overall. Only 15 original GWAS tag SNPs remain among the catalogue of candidate variants identified; the remainder are replaced by more likely candidates. Biological annotation of our credible set of variants indicates significant enrichment within promoter and enhancer elements, and transcription factor-binding sites, including AR, ERG and FOXA1. In 40 regions at least one variant is colocalised with an eQTL in prostate cancer tissue. The refined set of candidate variants substantially increase the proportion of familial relative risk explained by these known susceptibility regions, which highlights the importance of fine-mapping studies and has implications for clinical risk profiling.

The chromosome 6q22.33 region is associated with age at diagnosis of type 1 diabetes and disease risk in those diagnosed under 5 years of age

AIMS/HYPOTHESIS

The genetic risk of type 1 diabetes has been extensively studied. However, the genetic determinants of age at diagnosis (AAD) of type 1 diabetes remain relatively unexplained. Identification of AAD genes and pathways could provide insight into the earliest events in the disease process.

METHODS

Using ImmunoChip data from 15,696 cases, we aimed to identify regions in the genome associated with AAD.

RESULTS

Two regions were convincingly associated with AAD (p < 5 × 10-8): the MHC on 6p21, and 6q22.33. Fine-mapping of 6q22.33 identified two AAD-associated haplotypes in the region nearest to the genes encoding protein tyrosine phosphatase receptor kappa (PTPRK) and thymocyte-expressed molecule involved in selection (THEMIS). We examined the susceptibility to type 1 diabetes at these SNPs by performing a meta-analysis including 19,510 control participants. Although these SNPs were not associated with type 1 diabetes overall (p > 0.001), the SNP most associated with AAD, rs72975913, was associated with susceptibility to type 1 diabetes in those individuals diagnosed at less than 5 years old (p = 2.3 × 10-9).

CONCLUSION/INTERPRETATION

PTPRK and its neighbour THEMIS are required for early development of the thymus, which we can assume influences the initiation of autoimmunity. Non-HLA genes may only be detectable as risk factors for the disease in individuals diagnosed under the age 5 years because, after that period of immune development, their role in disease susceptibility has become redundant.

Chromosome contacts in activated T cells identify autoimmune disease candidate genes

Background

Autoimmune disease-associated variants are preferentially found in regulatory regions in immune cells, particularly CD4⁺ T cells. Linking such regulatory regions to gene promoters in disease-relevant cell contexts facilitates identification of candidate disease genes.

Results

Within 4 h, activation of CD4⁺ T cells invokes changes in histone modifications and enhancer RNA transcription that correspond to altered expression of the interacting genes identified by promoter capture Hi-C. By integrating promoter capture Hi-C data with genetic associations for five autoimmune diseases, we prioritised 245 candidate genes with a median distance from peak signal to prioritised gene of 153 kb. Just under half (108/245) prioritised genes related to activation-sensitive interactions. This included IL2RA, where allele-specific expression analyses were consistent with its interaction-mediated regulation, illustrating the utility of the approach.

Conclusions

Our systematic experimental framework offers an alternative approach to candidate causal gene identification for variants with cell state-specific functional effects, with achievable sample sizes.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-017-1285-0) contains supplementary material, which is available to authorized users.

The Genetic Architecture of Type 1 Diabetes

Type 1 diabetes (T1D) is classically characterised by the clinical need for insulin, the presence of disease-associated serum autoantibodies, and an onset in childhood. The disease, as with other autoimmune diseases, is due to the interaction of genetic and non-genetic effects, which induce a destructive process damaging insulin-secreting cells. In this review, we focus on the nature of this interaction, and how our understanding of that gene-environment interaction has changed our understanding of the nature of the disease. We discuss the early onset of the disease, the development of distinct immunogenotypes, and the declining heritability with increasing age at diagnosis. Whilst Human Leukocyte Antigens (HLA) have a major role in causing T1D, we note that some of these HLA genes have a protective role, especially in children, whilst other non-HLA genes are also important. In adult-onset T1D, the disease is often not insulin-dependent at diagnosis, and has a dissimilar immunogenotype with reduced genetic predisposition. Finally, we discuss the putative nature of the non-genetic factors and how they might interact with genetic susceptibility, including preliminary studies of the epigenome associated with T1D.

The genetics revolution in rheumatology: large scale genomic arrays and genetic mapping

Susceptibility to rheumatic diseases, such as osteoarthritis, rheumatoid arthritis, ankylosing spondylitis, systemic lupus erythematosus, juvenile idiopathic arthritis and psoriatic arthritis, includes a large genetic component. Understanding how an individual's genetic background influences disease onset and outcome can lead to a better understanding of disease biology, improved diagnosis and treatment, and, ultimately, to disease prevention or cure. The past decade has seen great progress in the identification of genetic variants that influence the risk of rheumatic diseases. The challenging task of unravelling the function of these variants is ongoing. In this Review, the major insights from genetic studies, gained from advances in technology, bioinformatics and study design, are discussed in the context of rheumatic disease. In addition, pivotal genetic studies in the main rheumatic diseases are highlighted, with insights into how these studies have changed the way we view these conditions in terms of disease overlap, pathways of disease and potential new therapeutic targets. Finally, the limitations of genetic studies, gaps in our knowledge and ways in which current genetic knowledge can be fully translated into clinical benefit are examined.

Widespread Allelic Heterogeneity in Complex Traits

Recent successes in genome-wide association studies (GWASs) make it possible to address important questions about the genetic architecture of complex traits, such as allele frequency and effect size. One lesser-known aspect of complex traits is the extent of allelic heterogeneity (AH) arising from multiple causal variants at a locus. We developed a computational method to infer the probability of AH and applied it to three GWASs and four expression quantitative trait loci (eQTL) datasets. We identified a total of 4,152 loci with strong evidence of AH. The proportion of all loci with identified AH is 4%-23% in eQTLs, 35% in GWASs of high-density lipoprotein (HDL), and 23% in GWASs of schizophrenia. For eQTLs, we observed a strong correlation between sample size and the proportion of loci with AH (R2 = 0.85, p = 2.2 × 10-16), indicating that statistical power prevents identification of AH in other loci. Understanding the extent of AH may guide the development of new methods for fine mapping and association mapping of complex traits.

Joint Bayesian inference of risk variants and tissue-specific epigenomic enrichments across multiple complex human diseases

Genome wide association studies (GWAS) provide a powerful approach for uncovering disease-associated variants in human, but fine-mapping the causal variants remains a challenge. This is partly remedied by prioritization of disease-associated variants that overlap GWAS-enriched epigenomic annotations. Here, we introduce a new Bayesian model RiVIERA (Risk Variant Inference using Epigenomic Reference Annotations) for inference of driver variants from summary statistics across multiple traits using hundreds of epigenomic annotations. In simulation, RiVIERA promising power in detecting causal variants and causal annotations, the multi-trait joint inference further improved the detection power. We applied RiVIERA to model the existing GWAS summary statistics of 9 autoimmune diseases and Schizophrenia by jointly harnessing the potential causal enrichments among 848 tissue-specific epigenomics annotations from ENCODE/Roadmap consortium covering 127 cell/tissue types and 8 major epigenomic marks. RiVIERA identified meaningful tissue-specific enrichments for enhancer regions defined by H3K4me1 and H3K27ac for Blood T-Cell specifically in the nine autoimmune diseases and Brain-specific enhancer activities exclusively in Schizophrenia. Moreover, the variants from the 95% credible sets exhibited high conservation and enrichments for GTEx whole-blood eQTLs located within transcription-factor-binding-sites and DNA-hypersensitive-sites. Furthermore, joint modeling the nine immune traits by simultaneously inferring and exploiting the underlying epigenomic correlation between traits further improved the functional enrichments compared to single-trait models.

Accuracy of Gene Scores when Pruning Markers by Linkage Disequilibrium

OBJECTIVE

Gene scores are often used to model the combined effects of genetic variants. When variants are in linkage disequilibrium, it is common to prune all variants except the most strongly associated. This avoids duplicating information but discards information when variants have independent effects. However, joint modelling of correlated variants increases the sampling error in the gene score. In recent applications, joint modelling has offered only small improvements in accuracy over pruning. We aimed to quantify the relationship between pruning and joint modelling in relation to sample size.

METHODS

We derived the coefficient of determination R2 for a gene score constructed from pruned markers, and for one constructed from correlated markers with jointly estimated effects.

RESULTS

Pruned scores tend to have slightly lower R2 than jointly modelled scores, but the differences are small at sample sizes up to 100,000. If the proportion of correlated variants is high, joint modelling can obtain modest improvements asymptotically.

CONCLUSIONS

The small gains observed to date from joint modelling can be explained by sample size. As studies become larger, joint modelling will be useful for traits affected by many correlated variants, but the improvements may remain small. Pruning remains a useful heuristic for current studies.

Genetic variation in MHC proteins is associated with T cell receptor expression biases

Within each individual, a highly diverse T cell receptor (TCR) repertoire interacts with peptides presented by major histocompatibility complex (MHC) molecules. Despite extensive research, it remains controversial whether germline-encoded TCR-MHC contacts promote TCR-MHC specificity and if so, whether there exist differences in TCR V-gene compatibilities with different MHC alleles. We applied eQTL mapping to test for associations between genetic variation and TCR V-gene usage in a large human cohort. We report strong trans associations between variation in the MHC locus and TCR V-gene usage. Fine mapping of the association signals reveals specific amino acids in MHC genes that bias V-gene usage, many of which contact or are spatially proximal to the TCR or peptide. Hence, these MHC variants, several of which are linked to autoimmune diseases, can directly affect TCR-MHC interaction. These results provide the first examples of trans-QTLs mediated by protein-protein interactions, and are consistent with intrinsic TCR-MHC specificity.

FINEMAP: efficient variable selection using summary data from genome-wide association studies

MOTIVATION

The goal of fine-mapping in genomic regions associated with complex diseases and traits is to identify causal variants that point to molecular mechanisms behind the associations. Recent fine-mapping methods using summary data from genome-wide association studies rely on exhaustive search through all possible causal configurations, which is computationally expensive.

RESULTS

We introduce FINEMAP, a software package to efficiently explore a set of the most important causal configurations of the region via a shotgun stochastic search algorithm. We show that FINEMAP produces accurate results in a fraction of processing time of existing approaches and is therefore a promising tool for analyzing growing amounts of data produced in genome-wide association studies and emerging sequencing projects.

AVAILABILITY AND IMPLEMENTATION

FINEMAP v1.0 is freely available for Mac OS X and Linux at http://www.christianbenner.com

CONTACT

: christian.benner@helsinki.fi or matti.pirinen@helsinki.fi.

Reducing GWAS Complexity

Genome-wide association studies (GWAS) have revealed numerous genomic 'hits' associated with complex phenotypes. In most cases these hits, along with surrogate genetic variation as measure by numerous single nucleotide polymorphisms (SNPs) that are in linkage disequilibrium, are not in coding genes making assignment of functionality or causality intractable. Here we propose that fine-mapping along with the matching of risk SNPs at chromatin biofeatures lessen this complexity by reducing the number of candidate functional/causal SNPs. For example, we show here that only on average 2 SNPs per prostate cancer risk locus are likely candidates for functionality/causality; we further propose that this manageable number should be taken forward in mechanistic studies. The candidate SNPs can be looked up for each prostate cancer risk region in 2 recent publications in 2015 (1,2) from our groups.