Search

Rates of sharing of genome-wide association studies (GWAS) summary statistics are historically low, limiting potential for scientific discovery. Here we show, using GWAS catalog data, that GWAS papers that share data get on average 81.8% more citations, an effect that is sustained over time., C.W. and G.R. are funded by the Wellcome Trust (WT220788). CW is funded by the Medical Research Council (MRC; MC UU 00002/4) and supported by the NIHR Cambridge BRC (BRC-1215-20014). The views expressed are those of the author(s) and not necessarily those of the NHS, the NIHR or the Department of Health and Social Care. The authors thank all GWAS researchers who proactively share their full summary statistics on GWAS catalog or elsewhere.

MotivationWhile many pipelines have been developed for calling genotypes using RNA-sequencing data, they all have adapted DNA genotype callers that do not model biases specific to RNA-sequencing such as reference panel bias or allele specific expression.ResultsHere, we present BBmix, a Bayesian Beta-Binomial mixture model that first learns the expected distribution of read counts for each genotype, and then deploys those learned parameters to call genotypes probabilistically. We benchmarked our model on a wide variety of datasets and showed that our method generally performed better than competitors, mainly due to an increase of up to 1.4% in the accuracy of heterozygous calls. Moreover, BBmix can be easily incorporated into standard pipelines for calling genotypes. We further show that parameters are generally transferable within datasets, such that a single learning run of less than one hour is sufficient to call genotypes in a large number of samples.AvailabilityWe implemented BBmix as an R package that is available for free under a GPL-2 licence athttps://gitlab.com/evigorito/bbmixand accompanying pipeline athttps://gitlab.com/evigorito/bbmix_pipeline.

Funder: GlaxoSmithKline foundation; doi: http://dx.doi.org/10.13039/501100002066, BACKGROUND: Genome-wide association studies (GWAS) are limited in power to detect associations that exceed the stringent genome-wide significance threshold. This limitation can be alleviated by leveraging relevant auxiliary data, such as functional genomic data. Frameworks utilising the conditional false discovery rate have been developed for this purpose, and have been shown to increase power for GWAS discovery whilst controlling the false discovery rate. However, the methods are currently only applicable for continuous auxiliary data and cannot be used to leverage auxiliary data with a binary representation, such as whether SNPs are synonymous or non-synonymous, or whether they reside in regions of the genome with specific activity states. RESULTS: We describe an extension to the cFDR framework for binary auxiliary data, called "Binary cFDR". We demonstrate FDR control of our method using detailed simulations, and show that Binary cFDR performs better than a comparator method in terms of sensitivity and FDR control. We introduce an all-encompassing user-oriented CRAN R package ( https://annahutch.github.io/fcfdr/ ; https://cran.r-project.org/web/packages/fcfdr/index.html ) and demonstrate its utility in an application to type 1 diabetes, where we identify additional genetic associations. CONCLUSIONS: Our all-encompassing R package, fcfdr, serves as a comprehensive toolkit to unite GWAS and functional genomic data in order to increase statistical power to detect genetic associations.

Background Cluster analysis is an integral part of precision medicine and systems biology, used to define groups of patients or biomolecules. Consensus clustering is an ensemble approach that is widely used in these areas, which combines the output from multiple runs of a non-deterministic clustering algorithm. Here we consider the application of consensus clustering to a broad class of heuristic clustering algorithms that can be derived from Bayesian mixture models (and extensions thereof) by adopting an early stopping criterion when performing sampling-based inference for these models. While the resulting approach is non-Bayesian, it inherits the usual benefits of consensus clustering, particularly in terms of computational scalability and providing assessments of clustering stability/robustness. Results In simulation studies, we show that our approach can successfully uncover the target clustering structure, while also exploring different plausible clusterings of the data. We show that, when a parallel computation environment is available, our approach offers significant reductions in runtime compared to performing sampling-based Bayesian inference for the underlying model, while retaining many of the practical benefits of the Bayesian approach, such as exploring different numbers of clusters. We propose a heuristic to decide upon ensemble size and the early stopping criterion, and then apply consensus clustering to a clustering algorithm derived from a Bayesian integrative clustering method. We use the resulting approach to perform an integrative analysis of three ’omics datasets for budding yeast and find clusters of co-expressed genes with shared regulatory proteins. We validate these clusters using data external to the analysis. Conclustions Our approach can be used as a wrapper for essentially any existing sampling-based Bayesian clustering implementation, and enables meaningful clustering analyses to be performed using such implementations, even when computational Bayesian inference is not feasible, e.g. due to poor exploration of the target density (often as a result of increasing numbers of features) or a limited computational budget that does not along sufficient samples to drawn from a single chain. This enables researchers to straightforwardly extend the applicability of existing software to much larger datasets, including implementations of sophisticated models such as those that jointly model multiple datasets.

Motivation Polygenic scores (PGS) aim to genetically predict complex traits at an individual level. PGS are typically trained on genome-wide association summary statistics and require an independent test dataset to tune parameters. More recent methods allow parameters to be tuned on the training data, removing the need for independent test data, but approaches are computationally intensive. Based on fine-mapping principles, we present RápidoPGS, a flexible and fast method to compute PGS requiring summary-level Genome-wide association studies (GWAS) datasets only, with little computational requirements and no test data required for parameter tuning. Results We show that RápidoPGS performs slightly less well than two out of three other widely used PGS methods (LDpred2, PRScs and SBayesR) for case–control datasets, with median r2 difference: -0.0092, -0.0042 and 0.0064, respectively, but up to 17 000-fold faster with reduced computational requirements. RápidoPGS is implemented in R and can work with user-supplied summary statistics or download them from the GWAS catalog. Availability and implementation Our method is available with a GPL license as an R package from CRAN and GitHub. Supplementary information Supplementary data are available at Bioinformatics online.

Takayasu arteritis is a rare inflammatory disease of large arteries. We performed a genetic study in Takayasu arteritis comprising 6,670 individuals (1,226 affected individuals) from five different populations. We discovered HLA risk factors and four non-HLA susceptibility loci in VPS8, SVEP1, CFL2, and chr13q21 and reinforced IL12B, PTK2B, and chr21q22 as robust susceptibility loci shared across ancestries. Functional analysis proposed plausible underlying disease mechanisms and pinpointed ETS2 as a potential causal gene for chr21q22 association. We also identified >60 candidate loci with suggestive association (p < 5 x 10(-s)) and devised a genetic risk score for Takayasu arteritis. Takayasu arteritis was compared to hundreds of other traits, revealing the closest genetic relatedness to inflammatory bowel disease. Epigenetic patterns within risk loci suggest roles for monocytes and B cells in Takayasu arteritis. This work enhances understanding of the genetic basis and pathophysiology of Takayasu arteritis and provides clues for potential new therapeutic targets., National Institute of Arthritis and Musculoskeletal and Skin Diseases of the National Institutes of HealthUnited States Department of Health & Human ServicesNational Institutes of Health (NIH) - USANIH National Institute of Arthritis & Musculoskeletal & Skin Diseases (NIAMS) [R01 AR070148]; National Institute of Arthritis and Musculoskeletal and Skin DiseasesUnited States Department of Health & Human ServicesNational Institutes of Health (NIH) - USANIH National Institute of Arthritis & Musculoskeletal & Skin Diseases (NIAMS) [U54 AR057319, U01 AR51874 04]; National Center for Research ResourcesUnited States Department of Health & Human ServicesNational Institutes of Health (NIH) - USANIH National Center for Research Resources (NCRR) [U54 RR019497]; Office of Rare Diseases Research of the National Center for Advancing Translational Sciences; Imperial College, National Institute for Health Research, Biomedical Research Centre; Wellcome TrustWellcome TrustEuropean Commission [WT107881]; Medical Research CouncilUK Research & Innovation (UKRI)Medical Research Council UK (MRC)European Commission [MC_UU_00002/4], This work was supported by the National Institute of Arthritis and Musculoskeletal and Skin Diseases of the National Institutes of Health grant R01 AR070148 to A.H.S. The Vasculitis Clinical Research Consortium has received support from the National Institute of Arthritis and Musculoskeletal and Skin Diseases (U54 AR057319 and U01 AR51874 04), the National Center for Research Resources (U54 RR019497), and the Office of Rare Diseases Research of the National Center for Advancing Translational Sciences. J.C.M., A.P.K., and R.M.M. acknowledge support from the Imperial College, National Institute for Health Research, Biomedical Research Centre. C.W. and G.R. acknowledge support from The Wellcome Trust (WT107881) and the Medical Research Council (MC_UU_00002/4). This work was supported by the use of study data downloaded from the dbGaP website, under dbGaP: phs000272.v1.p1, phs000431.v2.p1, phs000583.v1.p1, and phs000444.v1.p1.

Background: Autoimmunity may play a role in idiopathic pulmonary arterial hypertension (IPAH). It is not clear if this is causative, or as a bystander of disease and if it carries any prognostic or treatment significance. Methods: Comprehensive immune phenotyping was performed in IPAH patients, utilising patient samples from the UK National cohort study of Idiopathic and Heritable PAH. Standardised flow cytometric analysis of peripheral blood mononuclear cell types was assessed in IPAH patients and healthy controls. Serum autoantibody biomarkers were assayed using a Protoplex™ autoimmune panel in 473 type I PAH patients and 946 matched controls. Finally, a GST-fusion human proteomic screen was undertaken to identify novel circulating autoantibodies in PAH, and 350 PAH sera samples ELISA assayed for a putative antibody to BMPR2. Findings: Flow cytometric immune profiling demonstrates IPAH is associated with a signature consistent with altered humoral immune response. Increased autoantibody positivity was detected in PAH and patient subgroups associated with worse haemodynamic parameters but better overall survival. We identify a small subset of PAH patients who demonstrate immunoglobulin reactivity to the extracellular domain of the TGFbeta superfamily receptor, BMPR2, a protein which is strongly associated with hereditary causes of the disease. Interpretation: This evidence strongly suggests a proportion of patients with IPAH have an autoimmune component to the disease and that the BMPR2 pathway may be a target. Future trials targeting inflammation need to consider biomarker stratification. Funding Information: NIHR BRC, The Dinosaur Trust, MRC (MC_UU_00002/4) and the Wellcome Trust [WT107881]. Declarations of Interest: Gerry Coghlan reports grants and personal fees for conference attendance from Johnson & Johnson, personal fees for lectures from GlaxoSmithKline, Bayer and MSD, outside the submitted work. Paul Corris reports grants and personal fees for lectures and consultations from Actelion and Bayer, and personal fees for lectures and consultations from MSD, outside the submitted work. Luke Howard reports grants, personal fees for lectures, steering committee and advisory board work, and non-financial support for meeting attendance from Actelion, personal fees for lectures and advisory board work from Bayer and MSD, personal fees for advisory board work from GSK, outside the submitted work. Martin Johnson reports grants and personal fees for meeting attendance, lectures, and advisory board work from Actelion and MSD, outside the submitted work. David Kiely reports grants, personal fees, and other support from Actelion, Bayer and GSK, and personal fees and other support from MSD, outside the submitted work. Allan Lawrie reports grants from British Heart Foundation and Medical Research Council, grants, personal fees, and travel support from Actelion Pharmaceuticals Ltd, grants and personal fees from GlaxoSmithKline, outside the submitted work. James Lordan has received modest honoraria from Actelion and modest support to attend conference from Merck Sharp & Dohme, outside the submitted work. Joanna Pepke Zaba received fees for research grant, support to travel to meetings and honoraria from Actelion, Bayer, GlaxoSmithKline and Merck Sharp & Dohme, outside the submitted work. S.John Wort has received modest honoraria from GlaxoSmithKline, Actelion and Bayer and has received significant research grants from Actelion and Bayer, outside the submitted work. Nicholas Morrell reports personal fees from GSK and Johnson & Johnson/Actelion, outside the submitted work. Mark Toshner reports personal fees from GSK, grants and personal fees from J&J/Actelion, grants from Merck and Bayer, outside the submitted work. All other authors have nothing to disclose Ethics Approval Statement: Ethical approval was obtained, and informed written consent was given. For collection and processing of samples for autoantibody analysis, including controls ethical approval was obtained (An integrated approach to dissect the determinants, risk factors and pathways of ageing in the immune system 15/EE/0327, IRAS # 188383 and Genetics and epidemiology of aging-associated conditions in the Sardinian population, ProgeNIA PROT 2171/CE).

OBJECTIVES: Population-level measures of seropositivity are critical for understanding the epidemiology of an emerging pathogen, yet most antibody tests apply a strict cutoff for seropositivity that is not learnt in a data-driven manner, leading to uncertainty when classifying low-titer responses. To improve upon this, we evaluated cutoff-independent methods for their ability to assign likelihood of SARS-CoV-2 seropositivity to individual samples. METHODS: Using robust ELISAs based on SARS-CoV-2 spike (S) and the receptor-binding domain (RBD), we profiled antibody responses in a group of SARS-CoV-2 PCR+ individuals (n = 138). Using these data, we trained probabilistic learners to assign likelihood of seropositivity to test samples of unknown serostatus (n = 5100), identifying a support vector machines-linear discriminant analysis learner (SVM-LDA) suited for this purpose. RESULTS: In the training data from confirmed ancestral SARS-CoV-2 infections, 99% of participants had detectable anti-S and -RBD IgG in the circulation, with titers differing > 1000-fold between persons. In data of otherwise healthy individuals, 7.2% (n = 367) of samples were of uncertain serostatus, with values in the range of 3-6SD from the mean of pre-pandemic negative controls (n = 595). In contrast, SVM-LDA classified 6.4% (n = 328) of test samples as having a high likelihood (> 99% chance) of past infection, 4.5% (n = 230) to have a 50-99% likelihood, and 4.0% (n = 203) to have a 10-49% likelihood. As different probabilistic approaches were more consistent with each other than conventional SD-based methods, such tools allow for more statistically-sound seropositivity estimates in large cohorts. CONCLUSION: Probabilistic antibody testing frameworks can improve seropositivity estimates in populations with large titer variability.

Funder: GlaxoSmithKline; funder-id: http://dx.doi.org/10.13039/100004330, Genome-wide association studies (GWAS) have identified thousands of genetic variants that are associated with complex traits. However, a stringent significance threshold is required to identify robust genetic associations. Leveraging relevant auxiliary covariates has the potential to boost statistical power to exceed the significance threshold. Particularly, abundant pleiotropy and the non-random distribution of SNPs across various functional categories suggests that leveraging GWAS test statistics from related traits and/or functional genomic data may boost GWAS discovery. While type 1 error rate control has become standard in GWAS, control of the false discovery rate can be a more powerful approach. The conditional false discovery rate (cFDR) extends the standard FDR framework by conditioning on auxiliary data to call significant associations, but current implementations are restricted to auxiliary data satisfying specific parametric distributions, typically GWAS p-values for related traits. We relax these distributional assumptions, enabling an extension of the cFDR framework that supports auxiliary covariates from arbitrary continuous distributions (���Flexible cFDR���). Our method can be applied iteratively, thereby supporting multi-dimensional covariate data. Through simulations we show that Flexible cFDR increases sensitivity whilst controlling FDR after one or several iterations. We further demonstrate its practical potential through application to an asthma GWAS, leveraging various functional genomic data to find additional genetic associations for asthma, which we validate in the larger, independent, UK Biobank data resource.

Funder: “Expanding excellence in England” award from Research England, 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.

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.

In genome-wide association studies (GWAS) it is now common to search for, and find, multiple causal variants located in close proximity. It has also become standard to ask whether different traits share the same causal variants, but one of the popular methods to answer this question, coloc, makes the simplifying assumption that only a single causal variant exists for any given trait in any genomic region. Here, we examine the potential of the recently proposed Sum of Single Effects (SuSiE) regression framework, which can be used for fine-mapping genetic signals, for use with coloc. SuSiE is a novel approach that allows evidence for association at multiple causal variants to be evaluated simultaneously, whilst separating the statistical support for each variant conditional on the causal signal being considered. We show this results in more accurate coloc inference than other proposals to adapt coloc for multiple causal variants based on conditioning. We therefore recommend that coloc be used in combination with SuSiE to optimise accuracy of colocalisation analyses when multiple causal variants exist., Author summary Genetic association studies have found evidence that human disease risk or other traits are under the influence of genetic variants. As results of studies are made publicly available, more research focuses on whether different traits are under influence of the same variants, which may help us understand how variants lead to differences in disease risk. However, one of the popular methods to answer this question, coloc, makes the simplifying assumption that no two members of the set of causal variants for any one trait are close to each other. Here, we examine the potential of the recently proposed Sum of Single Effects (SuSiE) regression framework, for use with coloc. SuSiE is a novel approach that allows evidence for association at multiple causal variants in proximity to be evaluated simultaneously. We show this results in more accurate coloc inference than other proposals to adapt coloc for multiple causal variants based on conditioning. We therefore recommend that coloc be used in combination with SuSiE to optimise accuracy of colocalisation analyses when multiple causal variants exist.

Detecting genetic variants associated with traits (quantitative trait loci, QTL) requires genotyped study individuals. Here we describe BaseQTL, a Bayesian method that exploits allele-specific expression to map molecular QTL from sequencing reads (eQTL for gene expression) even when no genotypes are available. When used with genotypes to map eQTL, BaseQTL has lower error rates and increased power compared with existing QTL mapping methods. Running without genotypes limits how many tests can be performed, but due to the proximity of QTL variants to gene bodies, the 2.8% of variants within a 100 kB window that could be tested contained 26% of eQTL detectable with genotypes. eQTL effect estimates were invariably consistent between analyses performed with and without genotypes. Often, sequencing data may be generated in the absence of genotypes on patients and controls in differential expression studies, and we identified an apparent psoriasis-specific eQTL for GSTP1 in one such dataset, providing new insights into disease-dependent gene regulation. BaseQTL is a Bayesian method to map molecular QTL affecting allele-specific expression even when no genotypes are available. It is well suited to discover eQTLs hidden in a wealth of RNA-seq data to unravel molecular mechanisms underpinning disease.

Gene clustering and sample clustering are commonly used to find patterns in gene expression datasets. However, in heterogeneous samples (e.g. different tissues or disease states), genes may cluster differently. Biclustering algorithms aim to solve this issue by performing sample clustering and gene clustering simultaneously. Existing reviews of biclustering algorithms have yet to include a number of more recent algorithms and have based comparisons on simplistic simulated datasets without specific evaluation of biclusters in real datasets, using less robust metrics.In this study we compared four classes of sparse biclustering algorithms on a range of simulated and real datasets. In particular we use a knockout mouse RNA-seq dataset to evaluate each algorithm’s ability to simultaneously cluster genes and cluster samples across multiple tissues. We found that Bayesian algorithms with strict sparsity constraints had high accuracy on the simulated datasets and didn’t require any post-processing, but were considerably slower than other algorithm classes. We assessed whether non-negative matrix factorisation algorithms can be repurposed for biclustering and found that, although the raw output was poor, after using a sparsity-inducing post-processing procedure we introduce, one such algorithm was one of the most highly ranked on real datasets. We also exhibit the limitations of biclustering algorithms by varying the complexity of simulated datasets. The algorithms generally struggled on simulated datasets with a large number of implanted factors, or with a large number of genes. In real datasets, the algorithms rarely returned clusters containing samples from multiple tissues, which highlights the need for further thought in the design and analysis of multi-tissue studies to avoid differences between tissues dominating the analysis.Code to run the analysis is available at https://github.com/nichollskc/biclust_comp, including wrappers for each algorithm, implementations of evaluation metrics, and code to simulate datasets and perform pre- and post-processing. The full tables of results are available at https://doi.org/10.5281/zenodo.4317556

Background: In Sweden, social restrictions to contain SARS‐CoV‐2 have primarily relied upon voluntary adherence to a set of recommendations. Strict lockdowns have not been enforced, potentially affecting viral dissemination. To understand the levels of past SARS‐CoV‐2 infection in the Stockholm population before the start of mass vaccinations, healthy blood donors and pregnant women (n = 5,100) were sampled at random between 14 March 2020 and 28 February 2021. Methods: In this cross‐sectional prospective study, otherwise‐healthy blood donors (n = 2,600) and pregnant women (n = 2,500) were sampled for consecutive weeks (at four intervals) throughout the study period. Sera from all participants and a cohort of historical (negative) controls (n = 595) were screened for IgG responses against stabilized trimers of the SARS‐CoV‐2 spike (S) glycoprotein and the smaller receptor‐binding domain (RBD). As a complement to standard analytical approaches, a probabilistic (cut‐off independent) Bayesian framework that assigns likelihood of past infection was used to analyse data over time. Setting: Healthy participant samples were randomly selected from their respective pools through Karolinska University Hospital. The study was carried out in accordance with Swedish Ethical Review Authority: registration number 2020–01807. Participants: No participants were symptomatic at sampling, and blood donors were all over the age of 18. No additional metadata were available from the participants. Results: Blood donors and pregnant women showed a similar seroprevalence. After a steep rise at the start of the pandemic, the seroprevalence trajectory increased steadily in approach to the winter second wave of infections, approaching 15% of all individuals surveyed by 13 December 2020. By the end of February 2021, 19% of the population tested seropositive. Notably, 96% of seropositive healthy donors screened (n = 56) developed neutralizing antibody responses at titres comparable to or higher than those observed in clinical trials of SARS‐CoV‐2 spike mRNA vaccination, supporting that mild infection engenders a competent B‐cell response. Conclusions: These data indicate that in the first year since the start of community transmission, seropositivity levels in metropolitan in Stockholm had reached approximately one in five persons, providing important baseline seroprevalence information prior to the start of vaccination.

Many genetic mutations affecting phenotypes are presumed to do so via altering gene expression in particular cells or tissues, but identifying the specific genes involved has been challenging. A transcriptome-wide association study (TWAS) attempts to identify disease associated genes by first learning a predictive model on an eQTL dataset and then imputing gene expression levels into a larger genome-wide association study (GWAS). Finally, associations between predicted gene expressions and GWAS phenotype are identified.Here, we compared tree-based machine learning (ML) method of random forests (RF) with more widely used linear methods of lasso, ridge, and elastic net regression, for prediction of gene expression. We also developed a multi-task learning extension to RF which simultaneously makes use of information from multiple tissues (RF-MTL) and compared it to a multi-dataset version of lasso, the joint lasso, and to a single tissue RF. We found that for prediction of gene expression, RF, in general, outperformed linear approaches on our chosen eQTL dataset and that multi-tissue methods generally outperformed their single-tissue counterparts, with RF-MTL performing the best. Simulations showed that these benefits generally propagated to the next steps of the analysis, although highlighted that joint lasso had a tendency to erroneously identify genes in one tissue if there existed a disease signal for that gene in another.We tested all four methods on type 1 diabetes (T1D) GWAS and expression data for several immune cells and found that 46 genes were identified by at least one method, though only 7 by all methods. Joint lasso discovered the most T1D-associated genes, including 15 unique to that method, but this may reflect its higher false positive rate due to “overborrowing” information across tissues. RF-MTL found more unique associated genes than RF for 3 out 5 tissues. Compared to lasso-based analysis, the RF gene list was more likely to relate to T1D in an analysis of independent data types. We conclude that RF, both single- and multi-task version, is competitive and, for some cell types, superior to linear models conventionally used in the TWAS studies.Author summaryA transcriptome-wide association study (TWAS) is a way of integrating expression data and genome-wide association studies (GWAS), which allows for discovery of genes, rather than mutations, associated to traits of interest. In the TWAS framework, we first train predictive models on an eQTL dataset, then use these models to impute gene expression into a GWAS dataset. Finally, we look for significant associations between predicted gene expression and a GWAS trait. In this work, we compare non-linear method of random forests (RF) to linear models, customarily used in TWAS. Furthermore, we demonstrate that TWAS framework can naturally be extended to, and potentially benefit from, a multi-tissue setting, thereby taking advantage of the correlation between gene expression in different tissue types. We applied the RF, a selection of linear models, and the multi-tissue approaches to an eQTL dataset of monocytes and B cells and a large T1D GWAS. We found that RF outperform lasso in terms of predictive accuracy and the number of differentially expressed genes found, and that multi-dataset version of lasso discovered the most T1D-associated genes. Analysis of the gene lists produced for each method in independent data types (excluding genetic association data) showed all related to T1D, but that the RF methods ranked T1D higher in their lists than the linear methods. We conclude that RF is a useful addition to the TWAS tool box.

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.

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.

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 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. Electronic supplementary material The online version of this article (10.1007/s00125-017-4440-y) contains peer-reviewed but unedited supplementary material, which is available to authorised users.

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., Statistical fine-mapping to pinpoint likely causal variants in a genomic region is complicated by linkage disequilibrium (LD). Here, Asimit et al. compare stepwise and stochastic approaches to fine-mapping and propose a Bayesian multinomial stochastic search method which they apply to six immune-mediated diseases.

Cytokines are essential regulatory components of the immune system, and their aberrant levels have been linked to many disease states. Despite increasing evidence that cytokines operate in concert, many of the physiological interactions between cytokines, and the shared genetic architecture that underlies them, remain unknown. Here, we aimed to identify and characterize genetic variants with pleiotropic effects on cytokines. Using three population-based cohorts (n = 9,263), we performed multivariate genome-wide association studies (GWAS) for a correlation network of 11 circulating cytokines, then combined our results in meta-analysis. We identified a total of eight loci significantly associated with the cytokine network, of which two (PDGFRB and ABO) had not been detected previously. In addition, conditional analyses revealed a further four secondary signals at three known cytokine loci. Integration, through the use of Bayesian colocalization analysis, of publicly available GWAS summary statistics with the cytokine network associations revealed shared causal variants between the eight cytokine loci and other traits; in particular, cytokine network variants at the ABO, SERPINE2, and ZFPM2 loci showed pleiotropic effects on the production of immune-related proteins, on metabolic traits such as lipoprotein and lipid levels, on blood-cell-related traits such as platelet count, and on disease traits such as coronary artery disease and type 2 diabetes.

Methods for analysis of GWAS summary statistics have encouraged data sharing and democratised 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. 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, available as an open source R package, 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.

Background Hi-C and capture Hi-C (CHi-C) are used to map physical contacts between chromatin regions in cell nuclei using high-throughput sequencing. Analysis typically proceeds considering the evidence for contacts between each possible pair of fragments independent from other pairs. This can produce long runs of fragments which appear to all make contact with the same baited fragment of interest. Results We hypothesised that these long runs could result from a smaller subset of direct contacts and propose a new method, based on a Bayesian sparse variable selection approach, which attempts to fine map these direct contacts. Our model is conceptually novel, exploiting the spatial pattern of counts in CHi-C data. Although we use only the CHi-C count data in fitting the model, we show that the fragments prioritised display biological properties that would be expected of true contacts: for bait fragments corresponding to gene promoters, we identify contact fragments with active chromatin and contacts that correspond to edges found in previously defined enhancer-target networks; conversely, for intergenic bait fragments, we identify contact fragments corresponding to promoters for genes expressed in that cell type. We show that long runs of apparently co-contacting fragments can typically be explained using a subset of direct contacts consisting of

Genome-wide association studies have been tremendously successful in identifying genetic variants associated with complex diseases. The majority of association signals are intergenic and evidence is accumulating that a high proportion of signals lie in enhancer regions. We use Capture Hi-C to investigate, for the first time, the interactions between associated variants for four autoimmune diseases and their functional targets in B- and T-cell lines. Here we report numerous looping interactions and provide evidence that only a minority of interactions are common to both B- and T-cell lines, suggesting interactions may be highly cell-type specific; some disease-associated SNPs do not interact with the nearest gene but with more compelling candidate genes (for example, FOXO1, AZI2) often situated several megabases away; and finally, regions associated with different autoimmune diseases interact with each other and the same promoter suggesting common autoimmune gene targets (for example, PTPRC, DEXI and ZFP36L1)., There is evidence that a proportion of the polymorphisms identified by genome-wide association studies lie in enchancer regions. Here the authors use Capture Hi-C to investigate the interaction with targets in autoimmune disease, showing interactions can be long range and cell-type specific.

Identification of alterations in the cellular composition of the human immune system is key to understanding the autoimmune process. Recently, a subset of FOXP3+ cells with low CD25 expression was found to be increased in peripheral blood from systemic lupus erythematosus (SLE) patients, although its functional significance remains controversial. Here we find in comparisons with healthy donors that the frequency of FOXP3+ cells within CD127lowCD25low CD4+ T cells (here defined as CD25lowFOXP3+ T cells) is increased in patients affected by autoimmune disease of varying severity, from combined immunodeficiency with active autoimmunity, SLE to type 1 diabetes. We show that CD25lowFOXP3+ T cells share phenotypic features resembling conventional CD127lowCD25highFOXP3+ Tregs, including demethylation of the Treg-specific epigenetic control region in FOXP3, HELIOS expression, and lack of IL-2 production. As compared to conventional Tregs, more CD25lowFOXP3+HELIOS+ T cells are in cell cycle (33.0% vs 20.7% Ki-67+; P = 1.3 × 10−9) and express the late-stage inhibitory receptor PD-1 (67.2% vs 35.5%; P = 4.0 × 10−18), while having reduced expression of the early-stage inhibitory receptor CTLA-4, as well as other Treg markers, such as FOXP3 and CD15s. The number of CD25lowFOXP3+ T cells is correlated (P = 3.1 × 10−7) with the proportion of CD25highFOXP3+ T cells in cell cycle (Ki-67+). These findings suggest that CD25lowFOXP3+ T cells represent a subset of Tregs that are derived from CD25highFOXP3+ T cells, and are a peripheral marker of recent Treg expansion in response to an autoimmune reaction in tissues., Highlights • FOXP3+ compartment within CD127lowCD25low T cells is expanded in autoimmune patients. • Increased numbers of CD25lowFOXP3+ T cells are a circulating marker of autoimmunity. • CD25lowFOXP3+ HELIOS+ T cells are fully demethylated at the FOXP3 TSDR. • CD25lowFOXP3+ T cells could represent a terminal differentiation stage of regulatory T cells.

BackgroundAutoimmune 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.ResultsWithin four hours, 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.ConclusionsOur systematic experimental framework offers an alternative approach to candidate causal gene identification for variants with cell state-specific functional effects, with achievable sample sizes.

Many common diseases show wide phenotypic variation. We present a statistical method for determining whether phenotypically defined subgroups of disease cases represent different genetic architectures, in which disease-associated variants have different effect sizes in two subgroups. Our method models the genome-wide distributions of genetic association statistics with mixture Gaussians. We apply a global test without requiring explicit identification of disease-associated variants, thus maximizing power in comparison to standard variant-by-variant subgroup analysis. Where evidence for genetic subgrouping is found, we present methods for post hoc identification of the contributing genetic variants. We demonstrate the method on a range of simulated and test data sets, for which expected results are already known. We investigate subgroups of individuals with type 1 diabetes (T1D) defined by autoantibody positivity, establishing evidence for differential genetic architecture with positivity for thyroid-peroxidase-specific antibody, driven generally by variants in known T1D-associated genomic regions.

Summary Long-range interactions between regulatory elements and gene promoters play key roles in transcriptional regulation. The vast majority of interactions are uncharted, constituting a major missing link in understanding genome control. Here, we use promoter capture Hi-C to identify interacting regions of 31,253 promoters in 17 human primary hematopoietic cell types. We show that promoter interactions are highly cell type specific and enriched for links between active promoters and epigenetically marked enhancers. Promoter interactomes reflect lineage relationships of the hematopoietic tree, consistent with dynamic remodeling of nuclear architecture during differentiation. Interacting regions are enriched in genetic variants linked with altered expression of genes they contact, highlighting their functional role. We exploit this rich resource to connect non-coding disease variants to putative target promoters, prioritizing thousands of disease-candidate genes and implicating disease pathways. Our results demonstrate the power of primary cell promoter interactomes to reveal insights into genomic regulatory mechanisms underlying common diseases., Graphical Abstract, Highlights • High-resolution maps of promoter interactions in 17 human primary blood cell types • Interaction patterns are cell type specific and segregate with the hematopoietic tree • Promoter-interacting regions enriched for regulatory chromatin features and eQTLs • Promoter interactions link non-coding GWAS variants with putative target genes, This study deploys a promoter capture Hi-C approach in 17 primary blood cell types to match collaborating regulatory regions and identify genes regulated by noncoding disease-associated variants. Explore this and other papers at the Cell Press IHEC webportal at http://www.cell.com/consortium/IHEC.

Identification of candidate causal variants in regions associated with risk of common diseases is complicated by linkage disequilibrium (LD) and multiple association signals. Nonetheless, accurate maps of these variants are needed, both to fully exploit detailed cell specific chromatin annotation data to highlight disease causal mechanisms and cells, and for design of the functional studies that will ultimately be required to confirm causal mechanisms. We adapted a Bayesian evolutionary stochastic search algorithm to the fine mapping problem, and demonstrated its improved performance over conventional stepwise and regularised regression through simulation studies. We then applied it to fine map the established multiple sclerosis (MS) and type 1 diabetes (T1D) associations in the IL-2RA (CD25) gene region. For T1D, both stepwise and stochastic search approaches identified four T1D association signals, with the major effect tagged by the single nucleotide polymorphism, rs12722496. In contrast, for MS, the stochastic search found two distinct competing models: a single candidate causal variant, tagged by rs2104286 and reported previously using stepwise analysis; and a more complex model with two association signals, one of which was tagged by the major T1D associated rs12722496 and the other by rs56382813. There is low to moderate LD between rs2104286 and both rs12722496 and rs56382813 (r2 ≃ 0:3) and our two SNP model could not be recovered through a forward stepwise search after conditioning on rs2104286. Both signals in the two variant model for MS affect CD25 expression on distinct subpopulations of CD4+ T cells, which are key cells in the autoimmune process. The results support a shared causal variant for T1D and MS. Our study illustrates the benefit of using a purposely designed model search strategy for fine mapping and the advantage of combining disease and protein expression data., Author Summary Genetic association studies have identified many DNA sequence variants that associate with disease risk. By exploiting the known correlation that exists between neighbouring variants in the genome, inference can be extended beyond those individual variants tested to identify sets within which a causal variant is likely to reside. However, this correlation, particularly in the presence of multiple disease causing variants in relative proximity, makes disentangling the specific causal variants difficult. Statistical approaches to this fine mapping problem have traditionally taken a stepwise search approach, beginning with the most associated variant in a region, then iteratively attempting to find additional associated variants. We adapted a stochastic search approach that avoids this stepwise process and is explicitly designed for dealing with highly correlated predictors to the fine mapping problem. We showed in simulated data that it outperforms its stepwise counterpart and other variable selection strategies such as the lasso. We applied our approach to understand the association of two immune-mediated diseases to a region on chromosome 10p15. We identified a model for multiple sclerosis containing two variants, neither of which was found through a stepwise search, and functionally linked both of these to the neighbouring candidate gene, IL2RA, in independent data. Our approach can be used to aid fine mapping of other disease-associated regions, which is critical for design of functional follow-up studies required to understand the mechanisms through which genetic variants influence disease.

The appearance of circulating islet-specific autoantibodies before disease diagnosis is a hallmark of human type 1 diabetes (T1D), and suggests a role for B cells in the pathogenesis of the disease. Alterations in the peripheral B cell compartment have been reported in T1D patients; however, to date, such studies have produced conflicting results and have been limited by sample size. In this study, we have performed a detailed characterization of the B cell compartment in T1D patients (n = 45) and healthy controls (n = 46), and assessed the secretion of the anti-inflammatory cytokine interleukin (IL)-10 in purified B cells from the same donors. Overall, we found no evidence for a profound alteration of the B cell compartment or in the production of IL-10 in peripheral blood of T1D patients. We also investigated age-related changes in peripheral B cell subsets and confirmed the sharp decrease with age of transitional CD19(+) CD27(-) CD24(hi) CD38(hi) B cells, a subset that has recently been ascribed a putative regulatory function. Genetic analysis of the B cell compartment revealed evidence for association of the IL2-IL21 T1D locus with IL-10 production by both memory B cells (P = 6·4 × 10(-4) ) and islet-specific CD4(+) T cells (P = 2·9 × 10(-3) ). In contrast to previous reports, we found no evidence for an alteration of the B cell compartment in healthy individuals homozygous for the non-synonymous PTPN22 Trp(620) T1D risk allele (rs2476601; Arg(620) Trp). The IL2-IL21 association we have identified, if confirmed, suggests a novel role for B cells in T1D pathogenesis through the production of IL-10, and reinforces the importance of IL-10 production by autoreactive CD4(+) T cells.

Diagnosis of the autoimmune disease type 1 diabetes (T1D) is preceded by the appearance of circulating autoantibodies to pancreatic islets. However, almost nothing is known about events leading to this islet autoimmunity. Previous epidemiological and genetic data have associated viral infections and antiviral type I interferon (IFN) immune response genes with T1D. Here, we first used DNA microarray analysis to identify IFN-β–inducible genes in vitro and then used this set of genes to define an IFN-inducible transcriptional signature in peripheral blood mononuclear cells from a group of active systemic lupus erythematosus patients (n = 25). Using this predefined set of 225 IFN signature genes, we investigated the expression of the signature in cohorts of healthy controls (n = 87), patients with T1D (n = 64), and a large longitudinal birth cohort of children genetically predisposed to T1D (n = 109; 454 microarrayed samples). Expression of the IFN signature was increased in genetically predisposed children before the development of autoantibodies (P = 0.0012) but not in patients with established T1D. Upregulation of IFN-inducible genes was transient, temporally associated with a recent history of upper respiratory tract infections (P = 0.0064), and marked by increased expression of SIGLEC-1 (CD169), a lectin-like receptor expressed on CD14+ monocytes. DNA variation in IFN-inducible genes altered T1D risk (P = 0.007), as exemplified by IFIH1, one of the genes in our IFN signature for which increased expression is a known risk factor for disease. These findings identify transient increased expression of type I IFN genes in preclinical diabetes as a risk factor for autoimmunity in children with a genetic predisposition to T1D.

Motivation: Genome-wide association studies (GWAS) have identified many loci implicated in disease susceptibility. Integration of GWAS summary statistics (P-values) and functional genomic datasets should help to elucidate mechanisms. Results: We extended a non-parametric SNP set enrichment method to test for enrichment of GWAS signals in functionally defined loci to a situation where only GWAS P-values are available. The approach is implemented in VSEAMS, a freely available software pipeline. We use VSEAMS to identify enrichment of type 1 diabetes (T1D) GWAS associations near genes that are targets for the transcription factors IKZF3, BATF and ESRRA. IKZF3 lies in a known T1D susceptibility region, while BATF and ESRRA overlap other immune disease susceptibility regions, validating our approach and suggesting novel avenues of research for T1D. Availability and implementation: VSEAMS is available for download (http://github.com/ollyburren/vseams). Contact: chris.wallace@cimr.cam.ac.uk Supplementary information: Supplementary data are available at Bioinformatics online.