Abstract
Background Developing functional understanding into the causal molecular drivers of immunological disease is a critical challenge in genomic medicine. Here we systematically apply Mendelian randomization (MR), genetic colocalization, immune cell-type enrichment and phenome-wide association methods to investigate the effect of genetically predicted gene expression on 12 autoimmune and 4 cancer outcomes.
Results Using whole blood derived estimates for regulatory variants from the eQTLGen consortium (n=31,684) we constructed genetic risk scores (r2<0.1) for 10,104 genes. Applying the inverse-variance weighted Mendelian randomization method transcriptome-wide whilst accounting for linkage disequilibrium structure identified 773 unique genes with evidence of a genetically predicted effect on at least one disease outcome (P<4.81 × 10−5). We next undertook genetic colocalization to investigate whether these effects may be confined to specific cell-types using gene expression data derived from 18 types of immune cells. This highlighted many cell-type dependent effects, such as PRKCQ expression and asthma risk (posterior probability of association (PPA)=0.998), which was T-cell specific, as well as TPM3 expression and prostate cancer risk (PPA=0.821), which was restricted to monocytes. Phenome-wide analyses on 320 complex traits allowed us to explore the shared genetic architecture and prioritize key drivers of disease risk, such as CASP10 which provided evidence of an effect on 7 cancer-related outcomes. Similarly, these evaluations of pervasive pleiotropy may be valuable for evaluations of therapeutic targets to help identify potential adverse effects.
Conclusions Our atlas of results can be used to characterize known and novel loci in autoimmune disease and cancer susceptibility, both in terms of developing insight into cell-type dependent effects as well as dissecting shared genetic architecture and disease pathways. As exemplar, we have highlighted several key findings in this study, although similar evaluations can be conducted interactively at http://mrcieu.mrsoftware.org/immuno_MR/.
Competing Interest Statement
The authors have declared no competing interest.
Funding Statement
CP is supported by a Wellcome Trust PhD studentship in Molecular, Genetic and Lifecourse Epidemiology [108902/B/15/Z]. TGR is a UKRI Innovation Research Fellow (MR/S003886/1). This work was supported by the MRC Integrative Epidemiology Unit which receives funding from the UK Medical Research Council and the University of Bristol (MC_UU_00011/1).
Author Declarations
I confirm all relevant ethical guidelines have been followed, and any necessary IRB and/or ethics committee approvals have been obtained.
Yes
The details of the IRB/oversight body that provided approval or exemption for the research described are given below:
All ethical approval for datasets used in this study has been described in the corresponding studies they were obtained from.
All necessary patient/participant consent has been obtained and the appropriate institutional forms have been archived.
Yes
I understand that all clinical trials and any other prospective interventional studies must be registered with an ICMJE-approved registry, such as ClinicalTrials.gov. I confirm that any such study reported in the manuscript has been registered and the trial registration ID is provided (note: if posting a prospective study registered retrospectively, please provide a statement in the trial ID field explaining why the study was not registered in advance).
Yes
I have followed all appropriate research reporting guidelines and uploaded the relevant EQUATOR Network research reporting checklist(s) and other pertinent material as supplementary files, if applicable.
Yes
Data Availability
All the results from this study are available in our web application (http://mrcieu.mrsoftware.org/immuno_MR/), and code for analyses are available from the 'TwoSampleMR', 'MendelianRandomization' and 'coloc' R packages. Summary statistics for GWAS traits were obtained from publicly available resources as described in Supplementary Table 1 and also from the OpenGWAS platform (https://gwas.mrcieu.ac.uk/). Data on gene expression data from whole blood was downloaded from https://www.eqtlgen.org/ and cell-type specific gene expression was obtained from https://dice-database.org/ and http://dcc.blueprint-epigenome.eu/). The reference panel used in this study was created using data from the UK Biobank (app 15825).