Phenome-wide association study on miRNA-related sequence variants: the UK Biobank.

BACKGROUND: Genetic variants in the coding region could directly affect the structure and expression levels of genes and proteins. However, the importance of variants in the non-coding region, such as microRNAs (miRNAs), remain to be elucidated. Genetic variants in miRNA-related sequences could affect their biogenesis or functionality and ultimately affect disease risk. Yet, their implications and pleiotropic effects on many clinical conditions remain unknown. METHODS: Here, we utilised genotyping and hospital records data in the UK Biobank (N = 423,419) to investigate associations between 346 genetic variants in miRNA-related sequences and a wide range of clinical diagnoses through phenome-wide association studies. Further, we tested whether changes in blood miRNA expression levels could affect disease risk through colocalisation and Mendelian randomisation analysis. RESULTS: We identified 122 associations for six variants in the seed region of miRNAs, nine variants in the mature region of miRNAs, and 27 variants in the precursor miRNAs. These included associations with hypertension, dyslipidaemia, immune-related disorders, and others. Nineteen miRNAs were associated with multiple diagnoses, with six of them associated with multiple disease categories. The strongest association was reported between rs4285314 in the precursor of miR-3135b and celiac disease risk (odds ratio (OR) per effect allele increase = 0.37, P = 1.8 × 10-162). Colocalisation and Mendelian randomisation analysis highlighted potential causal role of miR-6891-3p in dyslipidaemia. CONCLUSIONS: Our study demonstrates the pleiotropic effect of miRNAs and offers insights to their possible clinical importance.

Insights into suicidal behavior among psoriatic arthritis patients: A systematic review and a genetic linkage disequilibrium analysis.

OBJECTIVES: To systematically assess the magnitude of suicidal behavior among PsA patients and identify associated risk factors. Also identify common genes or coinherited single nucleotide polymorphisms (SNPs) implicated in suicidal behavior and PsA. METHODS: Based on the PRISMA guidelines, we conducted a systematic literature review of the online databases PubMed/Medline, Web of Science, and EMBASE from inception to May 2022. Full-text original articles that describe suicidal behavior in PsA patients were eligible. All registered genome-wide association study (GWAS) data in the GWAS catalog database for PsA and psychiatric traits, such as suicidal behavior, and depression, were downloaded for further analysis. RESULTS: A total of 48 articles were identified, and 6 were relevant to the study question .Among the 122,160 PsA patients, 700 had suicidal behavior (0,57%). The range of age in one study was between 30 and 49 years, and 64% of PsA patients with suicidal behavior were female. Among 13,899 PsA patients 74 had suicidal ideation (0.53%) and 125 suicide attempts occurred (0.9%). In two studies, among 17,383 patients, 13 complete suicides occurred (0.07%). A genetic haplotype on chromosomal region 6p21.1, spanning from 29,597,596 to 32,251,264 Mb, contains predisposing SNPs for PsA and depression. 6p21.1-6p21.3 is the chromosomal region containing the HLA genes of classes I, II and III. CONCLUSION: Suicide behavior in PsA patients was associated with depression and other psychiatric comorbidities. Further evidence supports a genetic origin, at least partly. Awareness of these findings can help clinicians to recognize suicide behavior and prevent suicide attempts.

A Genetic Locus within the FMN1/GREM1 Gene Region Interacts with Body Mass Index in Colorectal Cancer Risk.

Colorectal cancer risk can be impacted by genetic, environmental, and lifestyle factors, including diet and obesity. Gene-environment interactions (G × E) can provide biological insights into the effects of obesity on colorectal cancer risk. Here, we assessed potential genome-wide G × E interactions between body mass index (BMI) and common SNPs for colorectal cancer risk using data from 36,415 colorectal cancer cases and 48,451 controls from three international colorectal cancer consortia (CCFR, CORECT, and GECCO). The G × E tests included the conventional logistic regression using multiplicative terms (one degree of freedom, 1DF test), the two-step EDGE method, and the joint 3DF test, each of which is powerful for detecting G × E interactions under specific conditions. BMI was associated with higher colorectal cancer risk. The two-step approach revealed a statistically significant G×BMI interaction located within the Formin 1/Gremlin 1 (FMN1/GREM1) gene region (rs58349661). This SNP was also identified by the 3DF test, with a suggestive statistical significance in the 1DF test. Among participants with the CC genotype of rs58349661, overweight and obesity categories were associated with higher colorectal cancer risk, whereas null associations were observed across BMI categories in those with the TT genotype. Using data from three large international consortia, this study discovered a locus in the FMN1/GREM1 gene region that interacts with BMI on the association with colorectal cancer risk. Further studies should examine the potential mechanisms through which this locus modifies the etiologic link between obesity and colorectal cancer. SIGNIFICANCE: This gene-environment interaction analysis revealed a genetic locus in FMN1/GREM1 that interacts with body mass index in colorectal cancer risk, suggesting potential implications for precision prevention strategies.

Variational Bayes for high-dimensional proportional hazards models with applications within gene expression.

MOTIVATION: Few Bayesian methods for analyzing high-dimensional sparse survival data provide scalable variable selection, effect estimation and uncertainty quantification. Such methods often either sacrifice uncertainty quantification by computing maximum a posteriori estimates, or quantify the uncertainty at high (unscalable) computational expense. RESULTS: We bridge this gap and develop an interpretable and scalable Bayesian proportional hazards model for prediction and variable selection, referred to as sparse variational Bayes. Our method, based on a mean-field variational approximation, overcomes the high computational cost of Markov chain Monte Carlo, whilst retaining useful features, providing a posterior distribution for the parameters and offering a natural mechanism for variable selection via posterior inclusion probabilities. The performance of our proposed method is assessed via extensive simulations and compared against other state-of-the-art Bayesian variable selection methods, demonstrating comparable or better performance. Finally, we demonstrate how the proposed method can be used for variable selection on two transcriptomic datasets with censored survival outcomes, and how the uncertainty quantification offered by our method can be used to provide an interpretable assessment of patient risk. AVAILABILITY AND IMPLEMENTATION: our method has been implemented as a freely available R package survival.svb (https://github.com/mkomod/survival.svb). SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Higher thyrotropin leads to unfavorable lipid profile and somewhat higher cardiovascular disease risk: evidence from multi-cohort Mendelian randomization and metabolomic profiling.

BACKGROUND: Observational studies suggest interconnections between thyroid status, metabolism, and risk of coronary artery disease (CAD), but causality remains to be proven. The present study aimed to investigate the potential causal relationship between thyroid status and cardiovascular disease and to characterize the metabolomic profile associated with thyroid status. METHODS: Multi-cohort two-sample Mendelian randomization (MR) was performed utilizing genome-wide significant variants as instruments for standardized thyrotropin (TSH) and free thyroxine (fT4) within the reference range. Associations between TSH and fT4 and metabolic profile were investigated in a two-stage manner: associations between TSH and fT4 and the full panel of 161 metabolomic markers were first assessed hypothesis-free, then directional consistency was assessed through Mendelian randomization, another metabolic profile platform, and in individuals with biochemically defined thyroid dysfunction. RESULTS: Circulating TSH was associated with 52/161 metabolomic markers, and fT4 levels were associated with 21/161 metabolomic markers among 9432 euthyroid individuals (median age varied from 23.0 to 75.4 years, 54.5% women). Positive associations between circulating TSH levels and concentrations of very low-density lipoprotein subclasses and components, triglycerides, and triglyceride content of lipoproteins were directionally consistent across the multivariable regression, MR, metabolomic platforms, and for individuals with hypo- and hyperthyroidism. Associations with fT4 levels inversely reflected those observed with TSH. Among 91,810 CAD cases and 656,091 controls of European ancestry, per 1-SD increase of genetically determined TSH concentration risk of CAD increased slightly, but not significantly, with an OR of 1.03 (95% CI 0.99-1.07; p value 0.16), whereas higher genetically determined fT4 levels were not associated with CAD risk (OR 1.00 per SD increase of fT4; 95% CI 0.96-1.04; p value 0.59). CONCLUSIONS: Lower thyroid status leads to an unfavorable lipid profile and a somewhat increased cardiovascular disease risk.

The link between attention deficit hyperactivity disorder (ADHD) symptoms and obesity-related traits: genetic and prenatal explanations.

Attention-deficit/hyperactivity disorder (ADHD) often co-occurs with obesity, however, the potential causality between the traits remains unclear. We examined both genetic and prenatal evidence for causality using Mendelian Randomisation (MR) and polygenic risk scores (PRS). We conducted bi-directional MR on ADHD liability and six obesity-related traits using summary statistics from the largest available meta-analyses of genome-wide association studies. We also examined the shared genetic aetiology between ADHD symptoms (inattention and hyperactivity) and body mass index (BMI) by PRS association analysis using longitudinal data from Northern Finland Birth Cohort 1986 (NFBC1986, n = 2984). Lastly, we examined the impact of the prenatal environment by association analysis of maternal pre-pregnancy BMI and offspring ADHD symptoms, adjusted for PRS of both traits, in NFBC1986 dataset. Through MR analyses, we found evidence for bidirectional causality between ADHD liability and obesity-related traits. PRS association analyses showed evidence for genetic overlap between ADHD symptoms and BMI. We found no evidence for a difference between inattention and hyperactivity symptoms, suggesting that neither symptom subtype is driving the association. We found evidence for association between maternal pre-pregnancy BMI and offspring ADHD symptoms after adjusting for both BMI and ADHD PRS (association p-value = 0.027 for inattention, p = 0.008 for hyperactivity). These results are consistent with the hypothesis that the co-occurrence between ADHD and obesity has both genetic and prenatal environmental origins.

Atopic dermatitis or eczema? Consequences of ambiguity in disease name for biomedical literature mining.

BACKGROUND: Biomedical research increasingly relies on computational approaches to extract relevant information from large corpora of publications. OBJECTIVE: To investigate the consequence of the ambiguity between the use of terms "Eczema" and "Atopic Dermatitis" (AD) from the Information Retrieval perspective, and its impact on meta-analyses, systematic reviews and text mining. METHODS: Articles were retrieved by querying the PubMed using terms 'eczema' (D003876) and "dermatitis, atopic" (D004485). We used machine learning to investigate the differences between the contexts in which each term is used. We used a decision tree approach and trained model to predict if an article would be indexed with eczema or AD tags. We used text-mining tools to extract biological entities associated with eczema and AD, and investigated the discrepancy regarding the retrieval of key findings according to the terminology used. RESULTS: Atopic dermatitis query yielded more articles related to veterinary science, biochemistry, cellular and molecular biology; the eczema query linked to public health, infectious disease and respiratory system. Medical Subject Headings terms associated with "AD" or "Eczema" differed, with an agreement between the top 40 lists of 52%. The presence of terms related to cellular mechanisms, especially allergies and inflammation, characterized AD literature. The metabolites mentioned more frequently than expected in articles with AD tag differed from those indexed with eczema. Fewer enriched genes were retrieved when using eczema compared to AD query. CONCLUSIONS AND CLINICAL RELEVANCE: There is a considerable discrepancy when using text mining to extract bio-entities related to eczema or AD. Our results suggest that any systematic approach (particularly when looking for metabolites or genes related to the condition) should be performed using both terms jointly. We propose to use decision tree learning as a tool to spot and characterize ambiguity, and provide the source code for disambiguation at https://github.com/cfrainay/ResearchCodeBase.

Integrating multi-OMICS data through sparse canonical correlation analysis for the prediction of complex traits: a comparison study.

MOTIVATION: Recent developments in technology have enabled researchers to collect multiple OMICS datasets for the same individuals. The conventional approach for understanding the relationships between the collected datasets and the complex trait of interest would be through the analysis of each OMIC dataset separately from the rest, or to test for associations between the OMICS datasets. In this work we show that integrating multiple OMICS datasets together, instead of analysing them separately, improves our understanding of their in-between relationships as well as the predictive accuracy for the tested trait. Several approaches have been proposed for the integration of heterogeneous and high-dimensional (p≫n) data, such as OMICS. The sparse variant of canonical correlation analysis (CCA) approach is a promising one that seeks to penalize the canonical variables for producing sparse latent variables while achieving maximal correlation between the datasets. Over the last years, a number of approaches for implementing sparse CCA (sCCA) have been proposed, where they differ on their objective functions, iterative algorithm for obtaining the sparse latent variables and make different assumptions about the original datasets. RESULTS: Through a comparative study we have explored the performance of the conventional CCA proposed by Parkhomenko et al., penalized matrix decomposition CCA proposed by Witten and Tibshirani and its extension proposed by Suo et al. The aforementioned methods were modified to allow for different penalty functions. Although sCCA is an unsupervised learning approach for understanding of the in-between relationships, we have twisted the problem as a supervised learning one and investigated how the computed latent variables can be used for predicting complex traits. The approaches were extended to allow for multiple (more than two) datasets where the trait was included as one of the input datasets. Both ways have shown improvement over conventional predictive models that include one or multiple datasets. AVAILABILITY AND IMPLEMENTATION: https://github.com/theorod93/sCCA. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

An Enrichment Analysis for Cardiometabolic Traits Suggests Non-Random Assignment of Genes to microRNAs.

MicroRNAs (miRNAs) regulate the expression of the majority of genes. However, it is not known whether they regulate genes in random or are organized according to their function. To this end, we chose cardiometabolic disorders as an example and investigated whether genes associated with cardiometabolic disorders are regulated by a random set of miRNAs or a limited number of them. Single-nucleotide polymorphisms (SNPs) reaching genome-wide level significance were retrieved from most recent genome-wide association studies on cardiometabolic traits, which were cross-referenced with Ensembl to identify related genes and combined with miRNA target prediction databases (TargetScan, miRTarBase, or miRecords) to identify miRNAs that regulate them. We retrieved 520 SNPs, of which 355 were intragenic, corresponding to 304 genes. While we found a higher proportion of genes reported from all GWAS that were predicted targets for miRNAs in comparison to all protein-coding genes (75.1%), the proportion was even higher for cardiometabolic genes (80.6%). Enrichment analysis was performed within each database. We found that cardiometabolic genes were over-represented in target genes for 29 miRNAs (based on TargetScan) and 3 miRNAs (miR-181a, miR-302d and miR-372) (based on miRecords) after Benjamini-Hochberg correction for multiple testing. Our work provides evidence for non-random assignment of genes to miRNAs and supports the idea that miRNAs regulate sets of genes that are functionally related.

Publisher Correction: Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits.

In the version of this article originally published, the name of author Martin H. de Borst was coded incorrectly in the XML. The error has now been corrected in the HTML version of the paper.

Genetic analysis of over 1 million people identifies 535 new loci associated with blood pressure traits.

High blood pressure is a highly heritable and modifiable risk factor for cardiovascular disease. We report the largest genetic association study of blood pressure traits (systolic, diastolic and pulse pressure) to date in over 1 million people of European ancestry. We identify 535 novel blood pressure loci that not only offer new biological insights into blood pressure regulation but also highlight shared genetic architecture between blood pressure and lifestyle exposures. Our findings identify new biological pathways for blood pressure regulation with potential for improved cardiovascular disease prevention in the future.

Regulatory T Cell Responses in Participants with Type 1 Diabetes after a Single Dose of Interleukin-2: A Non-Randomised, Open Label, Adaptive Dose-Finding Trial.

BACKGROUND: Interleukin-2 (IL-2) has an essential role in the expansion and function of CD4+ regulatory T cells (Tregs). Tregs reduce tissue damage by limiting the immune response following infection and regulate autoreactive CD4+ effector T cells (Teffs) to prevent autoimmune diseases, such as type 1 diabetes (T1D). Genetic susceptibility to T1D causes alterations in the IL-2 pathway, a finding that supports Tregs as a cellular therapeutic target. Aldesleukin (Proleukin; recombinant human IL-2), which is administered at high doses to activate the immune system in cancer immunotherapy, is now being repositioned to treat inflammatory and autoimmune disorders at lower doses by targeting Tregs. METHODS AND FINDINGS: To define the aldesleukin dose response for Tregs and to find doses that increase Tregs physiologically for treatment of T1D, a statistical and systematic approach was taken by analysing the pharmacokinetics and pharmacodynamics of single doses of subcutaneous aldesleukin in the Adaptive Study of IL-2 Dose on Regulatory T Cells in Type 1 Diabetes (DILT1D), a single centre, non-randomised, open label, adaptive dose-finding trial with 40 adult participants with recently diagnosed T1D. The primary endpoint was the maximum percentage increase in Tregs (defined as CD3+CD4+CD25highCD127low) from the baseline frequency in each participant measured over the 7 d following treatment. There was an initial learning phase with five pairs of participants, each pair receiving one of five pre-assigned single doses from 0.04 × 106 to 1.5 × 106 IU/m2, in order to model the dose-response curve. Results from each participant were then incorporated into interim statistical modelling to target the two doses most likely to induce 10% and 20% increases in Treg frequencies. Primary analysis of the evaluable population (n = 39) found that the optimal doses of aldesleukin to induce 10% and 20% increases in Tregs were 0.101 × 106 IU/m2 (standard error [SE] = 0.078, 95% CI = -0.052, 0.254) and 0.497 × 106 IU/m2 (SE = 0.092, 95% CI = 0.316, 0.678), respectively. On analysis of secondary outcomes, using a highly sensitive IL-2 assay, the observed plasma concentrations of the drug at 90 min exceeded the hypothetical Treg-specific therapeutic window determined in vitro (0.015-0.24 IU/ml), even at the lowest doses (0.040 × 106 and 0.045 × 106 IU/m2) administered. A rapid decrease in Treg frequency in the circulation was observed at 90 min and at day 1, which was dose dependent (mean decrease 11.6%, SE = 2.3%, range 10.0%-48.2%, n = 37), rebounding at day 2 and increasing to frequencies above baseline over 7 d. Teffs, natural killer cells, and eosinophils also responded, with their frequencies rapidly and dose-dependently decreased in the blood, then returning to, or exceeding, pretreatment levels. Furthermore, there was a dose-dependent down modulation of one of the two signalling subunits of the IL-2 receptor, the ß chain (CD122) (mean decrease = 58.0%, SE = 2.8%, range 9.8%-85.5%, n = 33), on Tregs and a reduction in their sensitivity to aldesleukin at 90 min and day 1 and 2 post-treatment. Due to blood volume requirements as well as ethical and practical considerations, the study was limited to adults and to analysis of peripheral blood only. CONCLUSIONS: The DILT1D trial results, most notably the early altered trafficking and desensitisation of Tregs induced by a single ultra-low dose of aldesleukin that resolves within 2-3 d, inform the design of the next trial to determine a repeat dosing regimen aimed at establishing a steady-state Treg frequency increase of 20%-50%, with the eventual goal of preventing T1D. TRIAL REGISTRATION: ISRCTN Registry ISRCTN27852285; ClinicalTrials.gov NCT01827735.

Protocol of the adaptive study of IL-2 dose frequency on regulatory T cells in type 1 diabetes (DILfrequency): a mechanistic, non-randomised, repeat dose, open-label, response-adaptive study.

INTRODUCTION: Type 1 diabetes (T1D) is caused by autoimmune destruction of the insulin-producing ß cells in the pancreatic islets, leading to insulinopenia and hyperglycaemia. Genetic analyses indicate that alterations of the interleukin-2 (IL-2) pathway mediating immune activation and tolerance predispose to T1D, specifically the polymorphic expression of the IL-2 receptor-α chain (CD25) on T lymphocytes. Replacement of physiological doses of IL-2 could restore self-tolerance and prevent further autoimmunity by enhancing the function of CD4(+) T regulatory cells (Tregs) to limit the activation of auto reactive T effector cells (Teffs). In this experimental medicine study, we use an adaptive trial design to determine the optimal dosing regimen for IL-2 to improve Treg function while limiting activation of Teffs in participants with T1D. METHODS AND ANALYSIS: The Adaptive study of IL-2 dose frequency on Tregs in type 1 diabetes(DILfrequency) is a mechanistic, non-randomised, repeat dose open-label, response-adaptive study of 36 participants with T1D. The objective is to establish the optimal dose and frequency of ultra-low dose IL-2: to increase Treg frequency within the physiological range, to increase CD25 expression on Tregs, without increasing CD4(+) Teffs. DILfrequency has an initial learning phase where 12 participants are allocated to six different doses and frequencies followed by an interim statistical analysis. After analysis of the learning phase, the Dose and Frequency Committee will select the optimal targets for Treg frequency, Treg CD25 expression and Teff frequency. Three groups of eight participants will be treated consecutively in the confirming phase. Each dose and frequency selected will be based on statistical analysis of all data collected from the previous groups. ETHICS: Ethical approval for DILfrequency was granted on 12 August 2014. RESULTS: The results of this study will be reported, through peer-reviewed journals, conference presentations and an internal organisational report. TRIAL REGISTRATION NUMBERS: NCT02265809, ISRCTN40319192, CRN17571.

Widespread seasonal gene expression reveals annual differences in human immunity and physiology.

Seasonal variations are rarely considered a contributing component to human tissue function or health, although many diseases and physiological process display annual periodicities. Here we find more than 4,000 protein-coding mRNAs in white blood cells and adipose tissue to have seasonal expression profiles, with inverted patterns observed between Europe and Oceania. We also find the cellular composition of blood to vary by season, and these changes, which differ between the United Kingdom and The Gambia, could explain the gene expression periodicity. With regards to tissue function, the immune system has a profound pro-inflammatory transcriptomic profile during European winter, with increased levels of soluble IL-6 receptor and C-reactive protein, risk biomarkers for cardiovascular, psychiatric and autoimmune diseases that have peak incidences in winter. Circannual rhythms thus require further exploration as contributors to various aspects of human physiology and disease.

Effective recruitment of participants to a phase I study using the internet and publicity releases through charities and patient organisations: analysis of the adaptive study of IL-2 dose on regulatory T cells in type 1 diabetes (DILT1D).

BACKGROUND: A barrier to the successful development of new disease treatments is the timely recruitment of participants to experimental medicine studies that are primarily designed to investigate biological mechanisms rather than evaluate clinical efficacy. The aim of this study was to analyse the performance of three recruitment sources and the effect of publicity events during the Adaptive study of IL-2 dose on regulatory T cells in type 1 diabetes (DILT1D). METHODS: The final study outcome, demography, disease duration, residence and the effect of publicity events on the performance of three recruitment sources (clinics, type 1 diabetes (T1D) disease register and the internet) were analysed from a bespoke DILT1D recruitment database. For the internet source, the origin of website hits in relation to publicity events was also evaluated. RESULTS: A total of 735 potentially eligible participants were approached to identify the final 45 DILT1D participants. A total of 477 (64%) were identified via the disease register, but only 59 (12%) responded to contact. A total of 317 individuals registered with the DILT1D study team. Self-referral via the study website generated 170 (54%) registered individuals and was the most popular and successful source, with 88 (28%) sourced from diabetes clinics and 59 (19%) from the disease register. Of those with known T1D duration (N = 272), the internet and clinics sources identified a larger number (57, 21%) of newly diagnosed T1D (<100 days post-diagnosis) compared to the register (1, 0.4%). The internet extended the geographical reach of the study, enabling both national and international participation. Targeted website posts and promotional events from organisations supporting T1D research and treatment during the trial were essential to the success of the internet recruitment strategy. CONCLUSIONS: Analysis of the DILT1D study recruitment outcomes illustrates the utility of an active internet recruitment strategy, supported by patient groups and charities, funding agencies and sponsors, in successfully conducting an early phase study in T1D. This recruitment strategy should now be evaluated in late-stage trials to develop treatments for T1D and other diseases. TRIAL REGISTRATION: NCT01827735 (registered: 4 April 2013); ISRCTN27852285 (registered: 23 March 2013); DRN767 (registered: 21 January 2013).

A method for gene-based pathway analysis using genomewide association study summary statistics reveals nine new type 1 diabetes associations.

Pathway analysis can complement point-wise single nucleotide polymorphism (SNP) analysis in exploring genomewide association study (GWAS) data to identify specific disease-associated genes that can be candidate causal genes. We propose a straightforward methodology that can be used for conducting a gene-based pathway analysis using summary GWAS statistics in combination with widely available reference genotype data. We used this method to perform a gene-based pathway analysis of a type 1 diabetes (T1D) meta-analysis GWAS (of 7,514 cases and 9,045 controls). An important feature of the conducted analysis is the removal of the major histocompatibility complex gene region, the major genetic risk factor for T1D. Thirty-one of the 1,583 (2%) tested pathways were identified to be enriched for association with T1D at a 5% false discovery rate. We analyzed these 31 pathways and their genes to identify SNPs in or near these pathway genes that showed potentially novel association with T1D and attempted to replicate the association of 22 SNPs in additional samples. Replication P-values were skewed (P=9.85×10-11) with 12 of the 22 SNPs showing P<0.05. Support, including replication evidence, was obtained for nine T1D associated variants in genes ITGB7 (rs11170466, P=7.86×10-9), NRP1 (rs722988, 4.88×10-8), BAD (rs694739, 2.37×10-7), CTSB (rs1296023, 2.79×10-7), FYN (rs11964650, P=5.60×10-7), UBE2G1 (rs9906760, 5.08×10-7), MAP3K14 (rs17759555, 9.67×10-7), ITGB1 (rs1557150, 1.93×10-6), and IL7R (rs1445898, 2.76×10-6). The proposed methodology can be applied to other GWAS datasets for which only summary level data are available.

Two novel pathway analysis methods based on a hierarchical model.

MOTIVATION: Over the past few years several pathway analysis methods have been proposed for exploring and enhancing the analysis of genome-wide association data. Hierarchical models have been advocated as a way to integrate SNP and pathway effects in the same model, but their computational complexity has prevented them being applied on a genome-wide scale to date. METHODS: We present two novel methods for identifying associated pathways. In the proposed hierarchical model, the SNP effects are analytically integrated out of the analysis, allowing computationally tractable model fitting to genome-wide data. The first method uses Bayes factors for calculating the effect of the pathways, whereas the second method uses a machine learning algorithm and adaptive lasso for finding a sparse solution of associated pathways. RESULTS: The performance of the proposed methods was explored on both simulated and real data. The results of the simulation study showed that the methods outperformed some well-established association methods: the commonly used Fisher's method for combining P-values and also the recently published BGSA. The methods were applied to two genome-wide association study datasets that aimed to find the genetic structure of platelet function and body mass index, respectively. The results of the analyses replicated the results of previously published pathway analysis of these phenotypes but also identified novel pathways that are potentially involved. AVAILABILITY: An R package is under preparation. In the meantime, the scripts of the methods are available on request from the authors.

Comparison of methods for competitive tests of pathway analysis.

It has been suggested that pathway analysis can complement single-SNP analysis in exploring genomewide association data. Pathway analysis incorporates the available biological knowledge of genes and SNPs and is expected to improve the chances of revealing the underlying genetic architecture of complex traits. Methods for pathway analysis can be classified as competitive (enrichment) or self-contained (association) according to the hypothesis tested. Although association tests are statistically more powerful than enrichment tests they can be difficult to calibrate because biases in analysis accumulate across multiple SNPs or genes. Furthermore, enrichment tests can be more scientifically relevant than association tests, as they detect pathways with relatively more evidence for association than the remaining genes. Here we show how some well known association tests can be simply adapted to test for enrichment, and compare their performance to some established enrichment tests. We propose versions of the Adaptive Rank Truncated Product (ARTP), Tail Strength Measure and Fisher's combination of p-values for testing the enrichment null hypothesis. We compare the behaviour of these proposed methods with the established Hypergeometric Test and Gene-Set Enrichment Analysis (GSEA). The results of the simulation study show that the modified version of the ARTP method has generally the best performance across the situations considered. The methods were also applied for finding enriched pathways for body mass index (BMI) and platelet function phenotypes. The pathway analysis of BMI identified the Vasoactive Intestinal Peptide pathway as significantly associated with BMI. This pathway has been previously reported as associated with BMI and the risk of obesity. The ARTP method was the method that identified the largest number of enriched pathways across all tested pathway databases and phenotypes. The simulation and data application results are in agreement with previous work on association tests and suggests that the ARTP should be preferred for both enrichment and association testing.