A transcriptome-wide Mendelian randomization study in isolated human immune cells highlights risk genes involved in viral infections and potential drug repurposing opportunities for schizophrenia
===================================================================================================================================================================================================

* David Stacey
* Liam Gaziano
* Preethi Eldi
* Catherine Toben
* Beben Benyamin
* S Hong Lee
* Elina Hyppönen

## ABSTRACT

**Background** Schizophrenia is a neurodevelopmental psychiatric disorder characterised by symptoms of psychosis, thought disorder, and flattened affect. Immune mechanisms are associated with schizophrenia, though the precise nature of this relationship (i.e., causal, correlated, consequential) and the mechanisms involved are not fully understood.

**Methods** To elucidate these mechanisms, we conducted a transcriptome-wide Mendelian randomization study using gene expression exposures from 29 human *cis*-eQTL datasets encompassing 11 unique immune cell types, all publicly available from the eQTL catalogue.

**Results** These analyses highlighted 196 genes, including 67 located within the human leukocyte antigen (HLA) region. Enrichment analyses indicated an over-representation of immune genes, which was driven by the HLA genes. Stringent validation and replication steps retained 61 candidate genes, 27 of which were the sole causal signals at their respective loci, thereby representing strong candidate effector genes at known risk loci. We highlighted *L3HYPDH* as a potential novel schizophrenia risk gene and *DPYD* and *MAPK3* as candidate drug repurposing targets. Futhermore, we performed follow-up analyses focused on one of the candidate effectors, interferon regulatory transcription factor 3 (*IRF3*), which coordinates interferon responses to viral infections. We found evidence of shared genetic aetiology between schizophrenia and autoimmune diseases at the *IRF3* locus, and a significant enrichment of IRF3 chromatin binding at known schizophrenia risk loci.

**Conclusions** Our findings highlight a novel schizophrenia risk gene, potential drug repurposing opportunities, and provide support for *IRF3* as a schizophrenia hub gene, which may play critical roles in mediating schizophrenia-autoimmune comorbidities and the impact of viral infections on schizophrenia risk.

Keywords
*   Schizophrenia
*   inflammation
*   immune system
*   Mendelian randomization
*   transcriptome
*   genetics

## INTRODUCTION

Schizophrenia (SCZ) is a devastating psychiatric illness affecting approximately 1 in 300 people worldwide(1), with onset generally occurring during late adolescence or early adulthood(2). It is characterised by several debilitating symptoms including psychosis (i.e., hallucinations, delusions), disorganised thought and speech, and flattened affect(3). Although SCZ is highly heritable (∼80%)(4), its underlying genetic aetiology is complex and heterogeneous. Moreover, approximately 30% of patients are resistant to currently available antipsychotics(5), which underscores the need to identify novel drug targets. Importantly, recent studies have shown that anti-inflammatory medications can alleviate SCZ symptoms in some patients, including treatment-resistant patients(6, 7).

Immune hypotheses of SCZ have been in existence for decades, supported by several lines of evidence. For example, some people with SCZ have abnormalities in the levels of circulating immune markers(8, 9) and autoantibodies(10). People with SCZ also have an increased risk of comorbid immune-mediated diseases, particularly autoimmune diseases(11); while inflammation of the brain – or encephalitis – can often present with psychiatric symptoms such as psychosis(12). Although the relationship between the immune system and SCZ is highly complex and likely confounded by several lifestyle factors (e.g., smoking) and medication effects (e.g., antipsychotics), findings from Mendelian randomization (MR) studies, which are robust to such confounds, suggest that several circulating immune markers may be causally associated with SCZ risk(13–15).

Of these immune markers, C-reactive protein (CRP) – a key modulator of the complement cascade – has received the most widespread attention, with several independent MR studies consistently showing that a genetic propensity to higher levels of circulating CRP is associated with a lower risk of SCZ(13, 15, 16). The precise mechanisms underlying this association are currently unclear, though they may be mediated by protection against infection(17). Indeed, early life or maternal viral/bacterial infection is a well-established risk factor not only for SCZ(18, 19) but also autoimmune diseases(20), which may explain a proportion of the observed SCZ-autoimmune comorbidities.

Over the past decade, GWASes have identified more than 200 independent risk loci for SCZ, of which the most highly significant (by many orders of magnitude) resides within the human leukocyte antigen (HLA) region(21). Fine-mapping studies have highlighted *C4A*, which encodes complement component 4a, as a candidate effector gene at this locus, potentially mediated by perturbed synaptic pruning during adolescence(22). Nevertheless, pinpointing the precise effector genes at complex disease risk loci is a non-trivial task, and so the extent to which other immune-related genes may mediate some of the known genetic risk is currently unclear. The continued integration of GWAS data with functional genomic data from immune tissues or cells may help to address this knowledge gap.

To this end, we performed a transcriptome-wide MR (TWMR) study of SCZ using publicly available expression quantitative trait locus (eQTL) data from 11 isolated immune cell types of both myeloid and lymphoid origin, including microglia. The TWMR approach utilises *cis*-eQTLs to instrument gene expression exposures at scale, and is a proven method for identifying causal genes for complex human diseases(23, 24). We identified several immune-related genes as causal candidates for SCZ, potentially by modulating immune responses to neurotropic viral infection. We performed follow-up analyses of one of the immune gene candidates, interferon regulatory factor 3 (*IRF3*), which encodes a transcription factor responsible for inducing interferon responses to viral infections. We show that SCZ shares genetic aetiology with autoimmune diseases at the *IRF3* locus, and that IRF3 chromatin binding sites are significantly enriched at SCZ risk loci. We propose *IRF3* to be a SCZ hub gene and a potentially druggable therapeutic target.

## METHODS

### *Cis*-eQTL and GWAS summary datasets

We extracted 29 *cis*-eQTL datasets derived from isolated human immune cells from the eQTL catalogue (**Table 1**) (25, 26). Overall, these datasets represented 11 unique immune cell types. We extracted summary data from the largest and most recently published (at the time of writing) SCZ GWAS(21) from the psychiatric genetics consortium (PGC) ([https://pgc.unc.edu/for-researchers/download-results/](https://pgc.unc.edu/for-researchers/download-results/)), which included up to 130,644 Europeans (53,386 cases). See **Supplementary Information** for a complete list of datasets we used and their sources.

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/07/12/2024.02.18.24303002/T1)

Table 1. Summary of human immune cell *cis*-eQTL datasets downloaded from the eQTL catalogue for transcriptome-wide Mendelian randomization and statistical colocalization analyses.

### Transcriptome-wide Mendelian randomization

We performed TWMR analyses using GSMR(24), as implemented in the Genome-Wide Complex Trait Analysis (GCTA) suite of tools (v1.94.1)(27). Prior to running GSMR, we first extracted all *cis*-eQTLs from the eQTL catalogue(25, 26) datasets with *p*-values lower than the estimated empirical *p*-values (calculated from 1,000 permutations) for the corresponding eGene based on the beta distribution. We obtained all pre-computed empirical *p*-values from the eQTL catalogue(25, 26). We supplied the pre-filtered exposure input files to GSMR and artificially set the GSMR significance threshold to *p*<0.05 (--gwas-thresh 0.05) to ensure all extracted *cis*-eQTLs were available for LD clumping within the GSMR workflow. Because *cis*-eQTL data generally yields fewer independent genetic instruments than well-powered GWASes, we set the minimum number of significant variants required for GSMR analysis to 1 (--gsmr-snp-min 1) and we applied a lenient LD (r2) threshold of 0.1 (--clump-r2 0.1) (GSMR natively adjusts estimates for correlated instruments). We set all other GSMR parameters to default.

To estimate variant correlations, we utilized an individual-level reference genotype panel comprising 4,994 participants from the INTERVAL population study(28). To minimise type I errors, we adjusted significance thresholds by applying a Bonferroni correction per eQTL dataset for the total number of eGenes tested (i.e., 0.05 / *n* eGenes). We plotted the GSMR estimates using the forestplot (v3.1.1) ([https://CRAN.R-project.org/package=forestplot](https://CRAN.R-project.org/package=forestplot)) and corrplot (v0.92)(29) R packages.

### Statistical colocalization

We performed statistical colocalization analyses using two software tools: (i) the Coloc R package (v5.1.0.1)(30); and (ii) a C++ implementation of pair-wise conditional and colocalization (PWCoCo) analysis(31). Coloc and PWCoCo have been described previously, but for brief descriptions, see **Supplementary Information**.

### Gene ontology term enrichment analyses

We conducted gene ontology (GO) term enrichment analyses using the functional annotation chart webtool available at the Database for Annotation, Visualization, and Integrated Discovery (DAVID)(32). To facilitate interpretation, we performed enrichment analyses focused solely on the gene ontology (GO) direct terms for biological process, molecular function, and cellular component. We unselected all other functional annotation sources before beginning analyses. For the background gene set, we extracted all unique eGenes across all eQTL catalogue datasets (**Table 1**) used in TWMR analyses.

### Phenome-wide association studies

We conducted single variant phenome-wide association studies (PheWASs) using the Open Targets Genetics(33) and the integrated epidemiology unit (IEU) OpenGWAS(34) platforms. To search the Open Targets Genetics platform, we used the webtool. For each variant of interest, we extracted all phenotypes from the ‘GWAS lead variants’ module, which contained phenotypes/studies with a lead variant linked with the variant of interest either through linkage disequilibrium (r2>0.5) or fine-mapping (i.e., appears in the 95% credible causal set). To search the IEU OpenGWAS platform, we used the ieugwasr R package (v0.1.5) with the *p*-value threshold set to *p*<5×10-8.

### Distance-based clumping

We ordered the 61 SCZ-associated variants by chromosome and position, and extracted all genes from the xHLA region (chr6:25,726,063-33,400,644; hg38)(35). We first assigned all xHLA region genes to a single mega-clump. For the remaining non-xHLA genes, we used the first gene to initiate the first clump and then sequentially checked the distance between the next gene and the preceding one. If the distance was <500Kb, we added the gene to the current clump; and if the distance was ≥500Kb, we used that gene to initiate the next clump. We continued this process until all genes were assigned to a clump.

### Enrichment analysis of regulatory data

To determine whether IRF3 ChIP-seq peaks were enriched at SCZ risk loci, we utilized the GWAS analysis of regulatory or functional information enrichment with LD correction (GARFIELD) v2 software tool(36), available at [https://www.ebi.ac.uk/birney-srv/GARFIELD/](https://www.ebi.ac.uk/birney-srv/GARFIELD/). For further information, see **Supplementary Information**.

## RESULTS

### Transcriptome-wide Mendelian randomization (TWMR) analyses of isolated human immune cell types highlights 196 candidate schizophrenia-associated genes

**Fig. 1** describes the overall study design and main analyses conducted. To identify immune genes with potential causal roles in SCZ, we performed TWMR analyses using summary data from 29 *cis*-eQTL datasets selected (see **Methods**) from the eQTL catalogue generated by 15 different studies from 11 unique immune cell types (**Table 1**). To facilitate future TWMR efforts focused on exposures from immune cells, we have provided as an immunogenomic resource all the genetic variants used to instrument genes in each of the eQTL catalogue datasets we tested (**Table S1**). Overall, these analyses highlighted a total of 196 significant genes (see **Methods, Fig. 2a, Table S2**), of which 156 (79.59%) were protein-coding and 67 (34.18%) were located within the extended HLA (xHLA) region as defined by Horton et al. (2004) (chr6:25,726,063-33,400,644; hg38)(35). Notably, this gene set did not include *CRP* or any genes encoding major cytokines. The absence of *CRP* was expected since it is specifically expressed in the liver, with only residual expression observed in immune cells(37).

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/12/2024.02.18.24303002/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/07/12/2024.02.18.24303002/F1)

Figure 1. Flow diagram summarising the overall study design and key findings from each analysis step.
GSMR: generalised summary-data-based Mendelian randomization, GTEx: genotype-tissue expression project, SCZ: schizophrenia, TWMR: transcriptome-wide Mendelian randomization; pLi: probability of loss of function intolerance.

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/12/2024.02.18.24303002/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/07/12/2024.02.18.24303002/F2)

Figure 2. Transcriptome-wide Mendelian randomization analyses identifies schizophrenia-associated candidate genes.
(a) Heatmap plotting the bxy (logOR) estimates of the causal effect of isolated immune cell gene expression exposures on schizophrenia risk (y-axis: *cis*-eQTL datasets; x-axis: genes). Gene expression exposures with robust evidence of colocalization (PPH4>0.8) with schizophrenia risk are shown with asterisks. Vertical red dashed lines indicate the location of the human leukocyte antigen region (HLA). (b,c) Gene ontology term enrichment analyses with (b) and without (c) the HLA region genes included.

To characterise this list of 196 genes, we performed gene ontology (GO) term enrichment analyses using the Database for Annotation, Visualization, and Integrated Discovery (DAVID) (see **Methods**)(32). An initial analysis including all 196 genes revealed significant (FDR<0.05) enrichment of 30 GO terms, most of which were immune-related (**Fig. 2b**, **Table S3**). However, none of these terms retained significance after removing the 67 HLA region genes from the test set (**Fig. 2c**, **Table S4**), demonstrating that the immune enrichment observed in our initial analysis was driven by the HLA region genes.

### Statistical colocalization analyses provide validation for 61 schizophrenia-associated genes, 27 of which are the sole causal signals at their respective loci

Since TWMR studies are susceptible to confounding by linkage disequilibrium (LD), we sought to validate our findings by performing statistical colocalization analyses (see **Methods**). Overall, we found robust evidence (PPH4>0.8) of colocalization between SCZ and at least one *cis*-eQTL dataset for 61 of 196 genes (31.12%) (**Fig. 2a**, **Fig. 3**, **Table S5**) – herein, we will refer to these colocalising genes as the validated SCZ-associated gene set.

![Figure 3.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/12/2024.02.18.24303002/F3.medium.gif)

[Figure 3.](http://medrxiv.org/content/early/2024/07/12/2024.02.18.24303002/F3)

Figure 3. The validated SCZ-associated gene set contains candidate SCZ effector genes at 38 non-overlapping genomic regions.
Ideogram of the human genome showing the locations of the 61 validated SCZ-associated gene sets. For brevity, chromosomes without a validated signal are not shown. The blue arrows indicate the extended human leukocyte antigen (xHLA) region, and genes with a probability of loss of function intolerance (pLi) > 0.9 are shown in red font.

Most of these 61 genes replicated in the whole blood *cis*-eQTL dataset from the eQTLgen consortium (**Fig. S1**, **Table** S6, see **Supplementary Information**), and many were mutationally constrained and associated with Mendelian or oligogenic diseases featuring neurodevelopmental, neurological, or immunological features (**Table S7**, see **Supplementary Information**). Moreover, most of the colocalizing genes resided outside of the extended HLA region, and the overall evidence for colocalization within the HLA region was sparse (**Fig. 2a**). Indeed, of the 67 xHLA region genes, only six of them (**Fig. 3**; *H4C13*, *APOM*, *RNF5*, *AGER*, *PBX2*, *GPSM3*) were in the validated SCZ-associated gene set, indicating that most of the associations observed in this region do not reflect causal effects on SCZ (at least not in the cell types tested here).

Horizontal pleiotropy is a major source of potential confounding in Mendelian randomization studies, occuring when an instrumental variable has an effect on the outcome outside of its effect on the exposure in question(52). Due to the extensive co-regulation of neighbouring genes, gene expression exposures can be particularly prone to horizontal pleiotropy(53). To assess this risk across our validated set of 61 SCZ-associated genes, we performed distance-based clumping using a 500kb window (see **Methods**). Including the xHLA region (chr6:25,726,063-33400644) as a single clump, this procedure revealed 38 non-overlapping clumps in total, 27 of which contained just a single gene from the validated SCZ-associated gene set (**Fig. 3**, **Table S8**). These 27 genes therefore represent high confidence effector genes for SCZ that are unlikely to be confounded by horizontal pleiotropy.

Of the remaining 11 bi- or multi-gene clumps, we were able to distinguish between genuine signals and those confounded by horizontal pleiotropy for three of them by manual curation (see **Supplementary Information**). This included the *NAGA*-*WBP2NL*-*SNU13* locus on chromosome 22, which contained two independent SCZ GWAS signals, one of which colocalized solely with *SNU13* mRNA expression in CD8+ T cells (**Fig. S2**).

### *L3HYPDH* is a candidate effector gene at a novel SCZ risk locus

Next, we sought to identify genes from the validated SCZ-associated set that did not reside at a genome-wide significant locus for SCZ, and may therefore represent novel SCZ risk genes. We intersected the genomic coordinates of the validated SCZ-associated gene set with the coordinates of the sentinel variants at the independent SCZ genome-wide significant loci reported by Trubetskoy et al. (2022), and found three genes that did not reside within 1Mb of a SCZ risk locus: *FOXN2*, *L3HYPDH*, *PIK3C2A*. We conducted look-ups in the Open Targets Genetics(33) and the IEU OpenGWAS(34) platforms to check whether any of these three genes or loci had been implicated in SCZ by any previous GWASs (see **Supplementary Results**, **Table S9**, **Table S10**). These look-ups revealed that the *FOXN2* and *PIK3C2A* signals, but not the *L3HYPDH* signal, had been detected in previous SCZ GWASs (see **Supplementary Results**, **Fig. S3**, **Table S11**, **Table S12**). Taken together, these findings suggest that *L3HYPDH* is a candidate effector gene at a novel SCZ risk locus.

### The schizophrenia-associated genes *DPYD* and *MAPK3* may present opportunities for drug repurposing

To identify potentially druggable targets of any of the proteins encoded by the 61 validated SCZ-associated genes, we performed gene-based look ups using the online DrugBank resource(54). This search revealed that the products of five of the genes (*DPYD*, *MAPK3*, *PCCB*, *NT5C2*, *FURIN*) were potentially druggable (**Table S13**).

The most robust evidence was for the dihydropyrimidine dehydrogenase (DPD) (*DPYD*) enzyme, which was directly and specifically targeted by a clinically approved drug called tegafur-uracil (**Table S13**). Tegafur-uracil is a chemotherapy drug combination primarily used in the treatment of bowel cancer. It comprises the prodrug tegafur, which is metabolised into 5-fluorouracil (5-FU) and incorporated into the DNA of dividing cancer cells, thereby killing them. However, since 5-FU is inactivated by DPD, tegafur-uracil also comprises an excess of uracil, which prolongs the half-life of 5-FU by competitively binding to DPD. As an alternative to uracil, tegafur is often combined with more potent inhibitors of DPD such as gimeracil. In our study, TWMR analyses indicated that higher expression levels of *DPYD* in a colocalizing (PPH4>0.8) eQTL macrophage (salmonella-treated) dataset(55) was associated with a higher risk of SCZ (β=0.120, se=0.026, *p*=3.60×10-6), and so compounds acting to reduce the activity of DPD like uracil and gimeracil may represent feasible therapeutic candidates.

Additionally, three approved compounds (arsenic trioxide, sulindac, cholecystokinin) are known to target the mitogen-activated protein kinase 3 (MAPK3), albeit with less specificity than the DPD-targeting compounds (**Table S13**). Most notably, sulindac, which inhibits MAPK3, is a nonsteroidal anti-inflammatory drug (NSAID) used in the treatment of acute and chronic inflammatory conditions. In our study, TWMR analyses indicated that higher *MAPK3* mRNA expression in two distinct colocalizing (PPH4>0.8) monocyte (naïve and LPS-treated) *cis*-eQTL datasets(42, 44) was associated with a higher risk of SCZ (naïve: β=0.229, se=0.046, *p*=5.29×10-7; LPS-treated: β=0.108, se=0.022, *p*=9.72×10-7), suggesting that sulindac may represent an opportunity for immune drug repurposing. Moreover, there may be other opportunities for the repurposing of immune drugs that have not yet been approved. For example, taribavirin and seliciclib, which target cytosolic purine 5’-nucleotidase (*NT5C2*) and MAPK3, respectively, are both currently under investigation for their efficacy in treating viral infections including hepatitis C and herpes simplex virus (**Table S13**).

### Interferon regulatory factor 3 is a schizophrenia hub gene and may contribute towards schizophrenia-autoimmune comorbidities

Next, we sought to further elucidate the immune component of SCZ by performing follow-up analyses focused on one of the genes from the validated SCZ-associated gene set. We selected the *IRF3* gene because it is a strong candidate effector gene at a robust genome-wide significant locus for SCZ on chromosome 19 (**Fig. 4a**)(21) and a well-established immune gene. *IRF3* encodes interferon regulatory factor 3, which is a transcription factor responsible for orchestrating interferon responses to viral infections(56). Moreover, rare loss of function variants affecting *IRF3* are associated with HSV-1-induced encephalitis(57), which is characterised by central inflammation and can present with neuropsychiatric symptoms including psychosis (**Table S7**)(12). In our TWMR analyses, we found an association between lower *IRF3* mRNA expression in natural killer (NK) cells(47) and a higher risk of SCZ (Fig. 4a,b; bxy= -0.32, se=0.07, *p*=5.96×10-6), and we also observed directionally consistent sub-threshold associations in both B cells(43) and neutrophils(39) (**Fig. 4b**).

![Figure 4.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/12/2024.02.18.24303002/F4.medium.gif)

[Figure 4.](http://medrxiv.org/content/early/2024/07/12/2024.02.18.24303002/F4)

Figure 4. 
The rs7251 schizophrenia risk locus is associated with peripheral and central *IRF3* mRNA expression as well as risk of autoimmune diseases. (a) Forest plot depicting the bxy (logOR) estimates of the causal effect of immune cell and whole blood *IRF3* mRNA expression on schizophrenia risk. (b) Stacked regional association plot showing the genetic association signals at the *IRF3* locus for schizophrenia (top) and *IRF3* mRNA expression in natural killer cells (bottom). (c) Forest plot showing the results from a PheWAS of rs7251 and a stacked bar plot depicting PPs of colocalization between schizophrenia and each rs7251-associated phenotype.

To identify additional diseases and traits that might be associated with the *IRF3* locus, we performed a PheWAS of rs7251 using the Open Targets Genetics(33) platform (see **Methods**). This revealed six genome-wide significant associations (*p*<5×10-8) with either rs7251 or a proxy, most of which were autoimmune diseases involving thyroid dysfunction (**Fig. 4c**). Importantly, the same allele (rs7251-C) was associated with a higher risk of both SCZ and autoimmune disease, as well as poorer cognitive performance (**Fig. 4c**). We found robust evidence of colocalization (PPH4>0.8) with SCZ for four of the six phenotypes: autoimmune thyroid disease, cognitive performance, thyroid problems (not cancer), and vitiligo (an autoimmune disease affecting the skin) (**Fig. 4c**). These findings suggest that *IRF3* may contribute towards the observed comorbidity between SCZ and autoimmune diseases, particularly those involving thyroid dysfunction.

Transcription factors are often hub genes responsible for orchestrating complex transcriptional networks that regulate important biological processes. Since *IRF3* is a transcription factor, we sought to assess the extent to which *IRF3* might regulate the expression of other SCZ risk genes. To explore this, we performed enrichment analyses using GARFIELD(36) (GWAS analysis of regulatory or functional enrichment with LD correction) by leveraging (i) full SCZ GWAS summary data from Trubetskoy et al. (2022)(21) and (ii) all available chromatin immunoprecipitation sequence (ChIP-seq) data for IRF3 from the ENCODE project(58) (see **Methods**). We found significant enrichment of SCZ risk loci (defined using genome-wide [*p*<5×10-8] and suggestive [*p*<1×10-5] significance thresholds) at IRF3 peaks from lymphoblastoid, epithelial, and neuroblastoma cells (**Fig. 5a**).

![Figure 5.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/07/12/2024.02.18.24303002/F5.medium.gif)

[Figure 5.](http://medrxiv.org/content/early/2024/07/12/2024.02.18.24303002/F5)

Figure 5. Interferon regulatory transcription factor 3 (IRF3) DNA binding sites are enriched at known schizophrenia risk loci.
(a) Forest plot showing beta estimates (se) from GARFIELD enrichment analyses. (b) Graph showing the 26 schizophrenia risk loci with an overlapping IRF3 DNA-binding site as evidenced by chromatin immunoprecipitation-sequence (ChIP-seq) data. Colocalising (PPH4>0.8) *cis*-eQTL signals were from whole blood (eQTLgen consortium). SCZ: schizophrenia.

Overall, IRF3 peaks overlapped 15 independent genome-wide significant (*p*<5×10-8) and 11 suggestive (*p*<1×10-5) SCZ risk loci in at least one of the three significant cell types (**Fig. 5b**). We found that 18 of these loci (11 genome-wide significant, 7 suggestive) overlapped with IRF3 peaks at one or more gene promoters (+/-1Kb of the transcription start site), several of which were cognate to genes from our set of validated SCZ-associated genes (e.g., *ZNF823* and *CREB3L4*), suggesting *IRF3* may regulate the expression of these and other genes. Accordingly, we found robust evidence of colocalization between the SCZ signals at the 26 regions with whole blood *cis*-eQTL signals from eQTLgen for a total of 49 genes (**Table S14**, **Fig. 5b**). Many of these colocalizing *cis*-eQTL signals were cognate to the IRF3 peak-overlapping promoters, though not all of the, which suggests IRF3 may also be binding to other regulatory elements (e.g., enhancers) at some of these SCZ risk loci. Taken together, these findings are consistent with a model whereby SCZ-associated variants at these 26 risk loci exert regulatory effects on local gene expression either by strengthening or weakening IRF3 binding sites.

## DISCUSSION

In this TWMR study, we sought to identify immune genes and elucidate immune mechanisms involved in SCZ by leveraging 29 *cis*-eQTL datasets generated by 15 different studies from 11 unique human immune cell types, all publicly available from the eQTL catalogue(25, 26). Overall, TWMR analyses highlighted 196 unique candidate genes, 67 of which resided within the extended HLA region. We applied a series of stringent validation and replication steps culminating in the prioritization of 61 validated SCZ-associated genes (**Fig. 3**), 30 of which were unlikely to be confounded by horizontal pleiotropy and therefore represent convincing candidate SCZ effector genes. We highlighted *L3HYPDH* as a potentially novel SCZ risk gene and *DPYD* and *MAPK3* as candidate drug repurposing targets. We also conducted in-depth follow-up analyses of *IRF3*, showing it to be a candidate SCZ hub gene and a potential driver of SCZ-autoimmune comorbidities. Finally, we have made available as an immunogenomic resource all the genetic variants used in this study to instrument gene expression exposures from isolated human immune cells (**Table S1**). We hope this catalogue will facilitate future TWMR screening efforts interested in uncovering immune mechanisms of psychiatric disorders and other immune-mediated phenotypes.

A landmark study led by investigators at GlaxoSmithKline used retrospective data to show that drugs targeting genes and proteins supported by genetic evidence are more likely to progress through the drug development pipeline than those without genetic evidence(59). This has led to several large-scale MR studies investigating gene expression and protein abundance exposures seeking to identify novel drug targets and drug repurposing opportunities(60, 61). In this TWMR study, we found that the DPD (encoded by *DPYD*) inhibitor gimeracil, an adjunctive medicine used to enhance the efficacy of a cancer medicine called tegafur, may represent a drug repurposing opportunity for SCZ. DPD is an enzyme involved in the metabolism of the pyrimidines uracil and thymine, which function as building blocks of nucleic acids, play roles in energy metabolism, and can bind to P2Y receptors(62). Moreover, people with DPD deficiency present with – amongst other symptoms – intellectual disability and impaired social functioning(63), suggesting that gimeracil may help to treat SCZ symptoms relating to cognitive and social impairment.

As a means to further elucidate the immune component of SCZ, we conducted deep-dive analyses focused on *IRF3*, a well characterised immune gene from our validated SCZ-associated gene set. *IRF3* is a downstream effector of the cyclic GMP-AMP synthase (cGAS)-stimulator of interferon genes (STING) pathway, and is responsible for coordinating interferon type I responses following the detection of viral DNA by cGAS(64). Our findings indicated that a genetic propensity to lower *IRF3* expression was associated with a higher risk of SCZ (**Fig. 4b**), which would be consistent with a blunted interferon response to viral infections as a potential SCZ risk factor. A blunted interferon response could lead to the incomplete resolution of viral infections, thereby promoting viral latency and subsequent chronic intermittent or low-grade inflammation(65). Inflammatory mechanisms such as these, particularly those affecting the brain, have long been hypothesised to increase the risk of SCZ or SCZ-like symptoms.

Accordingly, rare loss of function variants affecting *IRF3* lead to HSV-1-induced encephalitis, which is characterised by central inflammation and can present with neuropsychiatric symptoms including psychosis(12, 57). HSV-1 is one of several viruses implicated in the pathophysiology of SCZ, and it is a compelling candidate because it is neurotropic – that is, capable of infecting the brain. Indeed, HSV-1 initially infects epithelial cells at orofacial surfaces(66), after which it can travel afferently along sensory neurons from the skin into the brain and has been shown to specifically infect areas of the limbic system, where it can remain latent(67). Our findings implicated *IRF3* in NK cells, which are amongst the first responders to viral infections(68). We therefore propose a model whereby *IRF3* activity in NK cells might modulate SCZ risk by acting as an early buffering system protecting against the spread of neurotropic viral infections from the periphery into the brain. However, since *IRF3* is also expressed in the brain(69), we cannot rule out a direct brain-centric role for *IRF3* in modulating SCZ risk.

Notably, *IRF3* was not the only gene from the validated SCZ-associated gene set involved in the regulation of viral infections. *ZNF823* and *MAPK3* are both annotated to the HSV-1 infection (hsa05168) KEGG pathway(70), while *MAPK3* (and *IRF3*) was also found to be up-regulated in placentas from pregnancies with SARS-CoV-2 (COVID-19) maternal infection(23). *RNF5* has been shown to inhibit type I interferon responses in HSV keratitis through the cGAS-STING/IRF3 signalling pathway(71), and *FURIN* encodes a protease responsible for cleaving proteins produced by several viruses, including herpesvirus glycoprotein B produced by HSV, as well as the SARS-CoV-2 (COVID-19) spike protein(72). Moreover, a recent study has shown that *L3HYPDH*, potentially a novel SCZ risk gene according to this TWMR study, is induced by type I interferon signalling and inhibits the replication of enterovirus 71 (EV71) by interfering with the translation of EV71 proteins(73). The recurrence of this viral regulation theme suggests that this TWMR study, which focused solely on gene expression exposures from isolated human immune cells, has uncovered genes and mechanisms responsible for mediating the relationship between viral infections and SCZ.

There are several potential limitations to this study. First, our TWMR analysis was conducted with exposures from nearly 30 *cis*-eQTL datasets, posing an increased risk of type I errors. To minimise this risk, we applied a conservative Bonferroni correction to each dataset, (see **Methods**) complemented by stringent validation and replication analyses. Second, most of the *cis*-eQTL datasets from isolated immune cells were derived from relatively small sample sizes, which limits discovery power. Access to larger datasets in the future will facilitate more robust discoveries. Third, certain immune cell types tested were represented by a single *cis*-eQTL study (e.g., NK cells, microglia). Independent *cis*-eQTL datasets for these cell types are necessary for replication and to further explore their role in SCZ.

In conclusion, by conducting MR analyses using gene expression exposures from 29 *cis*-eQTL immune cell datasets followed by stringent validation steps, we highlighted 61 candidate effector genes for SCZ. Overall, response to viral infections was a recurring theme of these genes, and so future studies focusing on them will be well-placed to further elucidate the critical mechanisms mediating the established relationship between viral infections and SCZ. We exemplified this by conducting follow-up analyses of *IRF3*, which encodes a transcription factor responsible for orchestrating interferon responses to viral infections(56), and showed it to be a candidate SCZ hub gene that also modulates risk of autoimmune disease. Based on this, we propose that addressing unmet needs in vaccination (e.g., HSV-1) and the use of antiviral drugs may represent a viable path towards mitigating SCZ risk in the population.

## Supporting information

Supplementary Figures [[supplements/303002_file04.pptx]](pending:yes)

Supplementary Tables [[supplements/303002_file05.xlsx]](pending:yes)

Supplementary Information [[supplements/303002_file06.docx]](pending:yes)

## DATA AVAILABILITY

All data generated by this study are available as Supplementary Tables or will be provided upon request from the corresponding author. All transcriptomic exposure data used in this study are publicly available either from the eQTL catalogue ([https://www.ebi.ac.uk/eqtl/Data_access/](https://www.ebi.ac.uk/eqtl/Data_access/)) or the Genotype-Tissue Expression project ([https://gtexportal.org/home/downloads/adult-gtex/overview](https://gtexportal.org/home/downloads/adult-gtex/overview)). Full GWAS summary data for schizophrenia from Trubetskoy et al. (2022)(21) are publicly available from the Psychiatric Genomics Consortium (PGC) website ([https://pgc.unc.edu/for-researchers/download-results/](https://pgc.unc.edu/for-researchers/download-results/)).

## CONFLICTS OF INTEREST

All authors declare no conflicts of interest.

## ACKNOWLEDGEMENTS

We would like to thank Dr. Kaur Alasoo and colleagues for their excellent eQTL catalogue resource.

## Footnotes

*   Several changes implemented in response to peer-review comments including: (i) significantly reduced word count (particularly in the discussion section), (ii) replaced two results subsections with more pertinent ones (i.e., sections addressing the potential for horizontal pleiotropy, novelty, and drugability), and (iii) added an ideogram as a main figure showing the distribution of the 61 'validated' schizophrenia-associated genes.

*   Received February 18, 2024.
*   Revision received July 11, 2024.
*   Accepted July 12, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial-NoDerivs 4.0 International), CC BY-NC-ND 4.0, as described at [http://creativecommons.org/licenses/by-nc-nd/4.0/](http://creativecommons.org/licenses/by-nc-nd/4.0/)

## REFERENCES

1.  1.van Os J, Kapur S (2009): Schizophrenia. Lancet. 374:635–645.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(09)60995-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19700006&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269222600024&link_type=ISI) 

2.  2.Häfner H, Maurer K, Löffler W, Fätkenheuer B, an der Heiden W, Riecher-Rössler A, et al. (1994): The epidemiology of early schizophrenia. Influence of age and gender on onset and early course. Br J Psychiatry Suppl.29–38.
    
    

3.  3.Patel KR, Cherian J, Gohil K, Atkinson D (2014): Schizophrenia: overview and treatment options. P t. 39:638–645.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25210417&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

4.  4.Hilker R, Helenius D, Fagerlund B, Skytthe A, Christensen K, Werge TM, et al. (2018): Heritability of Schizophrenia and Schizophrenia Spectrum Based on the Nationwide Danish Twin Register. Biol Psychiatry. 83:492–498.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2017.08.017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

5.  5.Potkin SG, Kane JM, Correll CU, Lindenmayer JP, Agid O, Marder SR, et al. (2020): The neurobiology of treatment-resistant schizophrenia: paths to antipsychotic resistance and a roadmap for future research. NPJ Schizophr. 6:1.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41537-020-00114-3&link_type=DOI) 

6.  6.Laan W, Grobbee DE, Selten JP, Heijnen CJ, Kahn RS, Burger H (2010): Adjuvant aspirin therapy reduces symptoms of schizophrenia spectrum disorders: results from a randomized, double-blind, placebo-controlled trial. J Clin Psychiatry. 71:520–527.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.4088/JCP.09m05117yel&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20492850&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000277946400002&link_type=ISI) 

7.  7.Nitta M, Kishimoto T, Müller N, Weiser M, Davidson M, Kane JM, et al. (2013): Adjunctive use of nonsteroidal anti-inflammatory drugs for schizophrenia: a meta-analytic investigation of randomized controlled trials. Schizophr Bull. 39:1230–1241.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schbul/sbt070&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23720576&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000326189700010&link_type=ISI) 

8.  8.Miller BJ, Culpepper N, Rapaport MH (2014): C-reactive protein levels in schizophrenia: a review and meta-analysis. Clin Schizophr Relat Psychoses. 7:223–230.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3371/CSRP.MICU.020813&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23428789&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

9.  9.Wang AK, Miller BJ (2018): Meta-analysis of Cerebrospinal Fluid Cytokine and Tryptophan Catabolite Alterations in Psychiatric Patients: Comparisons Between Schizophrenia, Bipolar Disorder, and Depression. Schizophr Bull. 44:75–83.
    
    

10. 10.Jernbom Falk A, Galletly C, Just D, Toben C, Baune BT, Clark SR, et al. (2021): Autoantibody profiles associated with clinical features in psychotic disorders. Transl Psychiatry. 11:474.
    
    

11. 11.Benros ME, Pedersen MG, Rasmussen H, Eaton WW, Nordentoft M, Mortensen PB (2014): A nationwide study on the risk of autoimmune diseases in individuals with a personal or a family history of schizophrenia and related psychosis. Am J Psychiatry. 171:218–226.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2013.13010086&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24129899&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000330601100016&link_type=ISI) 

12. 12.Ford H, Griffith S, Warren N, Swayne A, Blum S, Butzkueven H, et al. (2022): Psychiatric manifestations of autoimmune encephalitis. Autoimmun Rev. 21:103145.
    
    

13. 13.Hartwig FP, Borges MC, Horta BL, Bowden J, Davey Smith G (2017): Inflammatory Biomarkers and Risk of Schizophrenia: A 2-Sample Mendelian Randomization Study. JAMA Psychiatry. 74:1226–1233.
    
    

14. 14.Perry BI, Upthegrove R, Kappelmann N, Jones PB, Burgess S, Khandaker GM (2021): Associations of immunological proteins/traits with schizophrenia, major depression and bipolar disorder: A bi-directional two-sample mendelian randomization study. Brain Behav Immun. 97:176–185.
    
    

15. 15.Prins BP, Abbasi A, Wong A, Vaez A, Nolte I, Franceschini N, et al. (2016): Investigating the Causal Relationship of C-Reactive Protein with 32 Complex Somatic and Psychiatric Outcomes: A Large-Scale Cross-Consortium Mendelian Randomization Study. PLoS Med. 13:e1001976.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pmed.1001976&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

16. 16.Ligthart S, Vaez A, Võsa U, Stathopoulou MG, de Vries PS, Prins BP, et al. (2018): Genome Analyses of >200,000 Individuals Identify 58 Loci for Chronic Inflammation and Highlight Pathways that Link Inflammation and Complex Disorders. Am J Hum Genet. 103:691–706.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2018.09.009&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30388399&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

17. 17.Ligthart S (2019): Commentary: CRP and schizophrenia: cause, consequence or confounding? Int J Epidemiol. 48:1514–1515.
    
    

18. 18.Brown AS, Derkits EJ (2010): Prenatal infection and schizophrenia: a review of epidemiologic and translational studies. Am J Psychiatry. 167:261–280.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2009.09030361&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=20123911&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000275056300008&link_type=ISI) 

19. 19.Nielsen PR, Benros ME, Mortensen PB (2014): Hospital contacts with infection and risk of schizophrenia: a population-based cohort study with linkage of Danish national registers. Schizophr Bull. 40:1526–1532.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/schbul/sbt200&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24379444&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

20. 20.Smatti MK, Cyprian FS, Nasrallah GK, Al Thani AA, Almishal RO, Yassine HM (2019): Viruses and Autoimmunity: A Review on the Potential Interaction and Molecular Mechanisms. Viruses. 11.
    
    

21. 21.Trubetskoy V, Pardiñas AF, Qi T, Panagiotaropoulou G, Awasthi S, Bigdeli TB, et al. (2022): Mapping genomic loci implicates genes and synaptic biology in schizophrenia. Nature. 604:502–508.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/S41586-022-04434-5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

22. 22.Sekar A, Bialas AR, de Rivera H, Davis A, Hammond TR, Kamitaki N, et al. (2016): Schizophrenia risk from complex variation of complement component 4. Nature. 530:177–183.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature16549&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26814963&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

23. 23.Ursini G, Di Carlo P, Mukherjee S, Chen Q, Han S, Kim J, et al. (2023): Prioritization of potential causative genes for schizophrenia in placenta. Nat Commun. 14:2613.
    
    

24. 24.Zhu Z, Zheng Z, Zhang F, Wu Y, Trzaskowski M, Maier R, et al. (2018): Causal associations between risk factors and common diseases inferred from GWAS summary data. Nat Commun. 9:224.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-017-02317-230&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

25. 25.Kerimov N, Hayhurst JD, Peikova K, Manning JR, Walter P, Kolberg L, et al. (2021): A compendium of uniformly processed human gene expression and splicing quantitative trait loci. Nat Genet. 53:1290–1299.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00924-w&link_type=DOI) 

26. 26.Kerimov N, Tambets R, Hayhurst JD, Rahu I, Kolberg P, Raudvere U, et al. (2023): eQTL Catalogue 2023: New datasets, X chromosome QTLs, and improved detection and visualisation of transcript-level QTLs. PLoS Genet. 19:e1010932.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1010932&link_type=DOI) 

27. 27.Yang J, Lee SH, Goddard ME, Visscher PM (2011): GCTA: a tool for genome-wide complex trait analysis. Am J Hum Genet. 88:76–82.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2010.11.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=21167468&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

28. 28.Di Angelantonio E, Thompson SG, Kaptoge S, Moore C, Walker M, Armitage J, et al. (2017): Efficiency and safety of varying the frequency of whole blood donation (INTERVAL): a randomised trial of 45L000 donors. Lancet. 390:2360–2371.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(17)31928-1&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

29. 29.Wei T, Simko V (2021): R package ’corrplot’: Visualization of a Correlation Matrix (Version 0.92).
    
    

30. 30.Giambartolomei C, Vukcevic D, Schadt EE, Franke L, Hingorani AD, Wallace C, et al. (2014): Bayesian test for colocalisation between pairs of genetic association studies using summary statistics. PLoS Genet. 10:e1004383.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1004383&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24830394&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

31. 31.Zheng J, Haberland V, Baird D, Walker V, Haycock PC, Hurle MR, et al. (2020): Phenome-wide Mendelian randomization mapping the influence of the plasma proteome on complex diseases. Nat Genet. 52:1122–1131.
    
    

32. 32.Sherman BT, Hao M, Qiu J, Jiao X, Baseler MW, Lane HC, et al. (2022): DAVID: a web server for functional enrichment analysis and functional annotation of gene lists (2021 update). Nucleic Acids Res. 50:W216–w221.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/NAR/GKAC194&link_type=DOI) 

33. 33.Ghoussaini M, Mountjoy E, Carmona M, Peat G, Schmidt EM, Hercules A, et al. (2021): Open Targets Genetics: systematic identification of trait-associated genes using large-scale genetics and functional genomics. Nucleic Acids Res. 49:D1311–d1320.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/NAR/GKAA840&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

34. 34.Elsworth B, Lyon M, Alexander T, Liu Y, Matthews P, Hallett J, et al. (2020): The MRC IEU OpenGWAS data infrastructure. bioRxiv.2020.2008.2010.244293.
    
    

35. 35.Horton R, Wilming L, Rand V, Lovering RC, Bruford EA, Khodiyar VK, et al. (2004): Gene map of the extended human MHC. Nat Rev Genet. 5:889–899.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nrg1489&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=15573121&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000225416800011&link_type=ISI) 

36. 36.Iotchkova V, Ritchie GRS, Geihs M, Morganella S, Min JL, Walter K, et al. (2019): GARFIELD classifies disease-relevant genomic features through integration of functional annotations with association signals. Nat Genet. 51:343–353.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0322-6&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30692680&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

37. 37.Pepys MB, Hirschfield GM (2003): C-reactive protein: a critical update. J Clin Invest. 111:1805–1812.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/JCI200318921&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12813013&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000183623600002&link_type=ISI) 

38. 38.Alasoo K, Rodrigues J, Mukhopadhyay S, Knights AJ, Mann AL, Kundu K, et al. (2018): Shared genetic effects on chromatin and gene expression indicate a role for enhancer priming in immune response. Nat Genet. 50:424–431.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0046-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29379200&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

39. 39.Chen L, Ge B, Casale FP, Vasquez L, Kwan T, Garrido-Martín D, et al. (2016): Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells. Cell. 167:1398–1414.e1324.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2016.10.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27863251&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

40. 40.Bossini-Castillo L, Glinos DA, Kunowska N, Golda G, Lamikanra AA, Spitzer M, et al. (2022): Immune disease variants modulate gene expression in regulatory CD4(+) T cells. Cell Genom. 2:None.
    
    

41. 41.Theusch E, Chen YI, Rotter JI, Krauss RM, Medina MW (2020): Genetic variants modulate gene expression statin response in human lymphoblastoid cell lines. BMC Genomics. 21:555.
    
    

42. 42.Momozawa Y, Dmitrieva J, Théâtre E, Deffontaine V, Rahmouni S, Charloteaux B, et al. (2018): IBD risk loci are enriched in multigenic regulatory modules encompassing putative causative genes. Nat Commun. 9:2427.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-018-04365-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29930244&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

43. 43.Fairfax BP, Makino S, Radhakrishnan J, Plant K, Leslie S, Dilthey A, et al. (2012): Genetics of gene expression in primary immune cells identifies cell type-specific master regulators and roles of HLA alleles. Nat Genet. 44:502–510.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2205&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22446964&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

44. 44.Fairfax BP, Humburg P, Makino S, Naranbhai V, Wong D, Lau E, et al. (2014): Innate immune activity conditions the effect of regulatory variants upon monocyte gene expression. Science. 343:1246949.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE2OiIzNDMvNjE3NS8xMjQ2OTQ5IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTIvMjAyNC4wMi4xOC4yNDMwMzAwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

45. 45.Gutierrez-Arcelus M, Lappalainen T, Montgomery SB, Buil A, Ongen H, Yurovsky A, et al. (2013): Passive and active DNA methylation and the interplay with genetic variation in gene regulation. Elife. 2:e00523.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.7554/eLife.00523&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23755361&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

46. 46.Lappalainen T, Sammeth M, Friedländer MR, t Hoen PA, Monlong J, Rivas MA, et al. (2013): Transcriptome and genome sequencing uncovers functional variation in humans. Nature. 501:506–511.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12531&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24037378&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000324826300049&link_type=ISI) 

47. 47.Gilchrist JJ, Makino S, Naranbhai V, Sharma PK, Koturan S, Tong O, et al. (2022): Natural Killer cells demonstrate distinct eQTL and transcriptome-wide disease associations, highlighting their role in autoimmunity. Nat Commun. 13:4073.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-022-31626-4&link_type=DOI) 

48. 48.Kasela S, Kisand K, Tserel L, Kaleviste E, Remm A, Fischer K, et al. (2017): Pathogenic implications for autoimmune mechanisms derived by comparative eQTL analysis of CD4+ versus CD8+ T cells. PLoS Genet. 13:e1006643.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pgen.1006643&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=28248954&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

49. 49.Naranbhai V, Fairfax BP, Makino S, Humburg P, Wong D, Ng E, et al. (2015): Genomic modulators of gene expression in human neutrophils. Nat Commun. 6:7545.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ncomms8545&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26151758&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

50. 50.Buil A, Brown AA, Lappalainen T, Viñuela A, Davies MN, Zheng HF, et al. (2015): Gene-gene and gene-environment interactions detected by transcriptome sequence analysis in twins. Nat Genet. 47:88–91.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3162&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25436857&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

51. 51.Young AMH, Kumasaka N, Calvert F, Hammond TR, Knights A, Panousis N, et al. (2021): A map of transcriptional heterogeneity and regulatory variation in human microglia. Nat Genet. 53:861–868.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00875-2&link_type=DOI) 

52. 52.Verbanck M, Chen CY, Neale B, Do R (2018): Detection of widespread horizontal pleiotropy in causal relationships inferred from Mendelian randomization between complex traits and diseases. Nat Genet. 50:693–698.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0099-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29686387&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

53. 53.Porcu E, Rüeger S, Lepik K, Santoni FA, Reymond A, Kutalik Z (2019): Mendelian randomization integrating GWAS and eQTL data reveals genetic determinants of complex and clinical traits. Nat Commun. 10:3300.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-019-10936-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31341166&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

54. 54.Wishart DS, Knox C, Guo AC, Cheng D, Shrivastava S, Tzur D, et al. (2008): DrugBank: a knowledgebase for drugs, drug actions and drug targets. Nucleic Acids Res. 36:D901–906.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkm958&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18048412&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252545400158&link_type=ISI) 

55. 55.Alasoo K, Rodrigues J, Mukhopadhyay S, Knights AJ, Mann AL, Kundu K, et al. (2018): Shared genetic effects on chromatin and gene expression indicate a role for enhancer priming in immune response. Nature Genetics. 50:424–431.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0046-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29379200&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

56. 56.Hiscott J, Pitha P, Genin P, Nguyen H, Heylbroeck C, Mamane Y, et al. (1999): Triggering the interferon response: the role of IRF-3 transcription factor. J Interferon Cytokine Res. 19:1–13.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/107999099314360&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10048763&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000078437900001&link_type=ISI) 

57. 57.Andersen LL, Mørk N, Reinert LS, Kofod-Olsen E, Narita R, Jørgensen SE, et al. (2015): Functional IRF3 deficiency in a patient with herpes simplex encephalitis. J Exp Med. 212:1371–1379.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MzoiamVtIjtzOjU6InJlc2lkIjtzOjEwOiIyMTIvOS8xMzcxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTIvMjAyNC4wMi4xOC4yNDMwMzAwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

58. 58.Consortium TEP (2012): An integrated encyclopedia of DNA elements in the human genome. Nature. 489:57–74.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature11247&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22955616&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308347000039&link_type=ISI) 

59. 59.Nelson MR, Tipney H, Painter JL, Shen J, Nicoletti P, Shen Y, et al. (2015): The support of human genetic evidence for approved drug indications. Nat Genet. 47:856–860.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3314&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26121088&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

60. 60.Burgess S, Mason AM, Grant AJ, Slob EAW, Gkatzionis A, Zuber V, et al. (2023): Using genetic association data to guide drug discovery and development: Review of methods and applications. Am J Hum Genet. 110:195–214.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2022.12.017&link_type=DOI) 

61. 61.Schmidt AF, Finan C, Gordillo-Marañón M, Asselbergs FW, Freitag DF, Patel RS, et al. (2020): Genetic drug target validation using Mendelian randomisation. Nat Commun. 11:3255.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-16969-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32591531&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

62. 62.Wu HL, Gong Y, Ji P, Xie YF, Jiang YZ, Liu GY (2022): Targeting nucleotide metabolism: a promising approach to enhance cancer immunotherapy. J Hematol Oncol. 15:45.
    
    

63. 63.van Kuilenburg ABP, Meijer J, Meinsma R, Pérez-Dueñas B, Alders M, Bhuiyan ZA, et al. (2019): Dihydropyrimidine Dehydrogenase Deficiency: Homozygosity for an Extremely Rare Variant in DPYD due to Uniparental Isodisomy of Chromosome 1. JIMD Rep. 45:65–69.
    
    

64. 64.Liu S, Cai X, Wu J, Cong Q, Chen X, Li T, et al. (2015): Phosphorylation of innate immune adaptor proteins MAVS, STING, and TRIF induces IRF3 activation. Science. 347:aaa2630.
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE2OiIzNDcvNjIyNy9hYWEyNjMwIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDcvMTIvMjAyNC4wMi4xOC4yNDMwMzAwMi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

65. 65.Channappanavar R, Fehr AR, Zheng J, Wohlford-Lenane C, Abrahante JE, Mack M, et al. (2019): IFN-I response timing relative to virus replication determines MERS coronavirus infection outcomes. J Clin Invest. 129:3625–3639.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1172/jci126363&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31355779&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

66. 66.Arduino PG, Porter SR (2008): Herpes Simplex Virus Type 1 infection: overview on relevant clinico-pathological features. J Oral Pathol Med. 37:107–121.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18197856&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000252323300008&link_type=ISI) 

67. 67.Klein HC (2017): Silencing of Viral Elements: An Available Cure for Schizophrenia? Front Psychiatry. 8:284.
    
    

68. 68.Björkström NK, Strunz B, Ljunggren H-G (2022): Natural killer cells in antiviral immunity. Nature Reviews Immunology. 22:112–123.
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

69. 69.GTExconsortium (2013): The Genotype-Tissue Expression (GTEx) project. Nat Genet. 45:580–585.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2653&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23715323&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 

70. 70.Kanehisa M, Goto S (2000): KEGG: kyoto encyclopedia of genes and genomes. Nucleic Acids Res. 28:27–30.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/28.1.27&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10592173&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F07%2F12%2F2024.02.18.24303002.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000084896300007&link_type=ISI) 

71. 71.Liu Z, Xia L (2022): E3 ligase RNF5 inhibits type I interferon response in herpes simplex virus keratitis through the STING/IRF3 signaling pathway. Front Microbiol. 13:944101.
    
    

72. 72.Braun E, Sauter D (2019): Furin-mediated protein processing in infectious diseases and cancer. Clin Transl Immunology. 8:e1073.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/cti2.1073&link_type=DOI) 

73. 73.Liu J, Liu L, Zeng S, Meng X, Lei N, Yang H, et al. (2024): Inhibition of EV71 replication by an interferon-stimulated gene product L3HYPDH. Virus Res. 342:199336.