Genetic variants associated with cross-disorder and disorder-specific risk for psychiatric disorders are enriched at epigenetically active sites in peripheral lymphoid cells
=============================================================================================================================================================================

* Mary-Ellen Lynall
* Blagoje Soskic
* James Hayhurst
* Jeremy Schwartzentruber
* Daniel F. Levey
* Gita A. Pathak
* Renato Polimanti
* Joel Gelernter
* Murray B. Stein
* Gosia Trynka
* Menna R. Clatworthy
* Ed Bullmore

## Abstract

Multiple psychiatric disorders have been associated with abnormalities in both the innate and adaptive immune systems. The role of these abnormalities in pathogenesis, and whether they are driven by psychiatric risk variants, remains unclear. We tested for enrichment of GWAS variants associated with multiple psychiatric disorders (cross-disorder or trans-diagnostic risk), or 5 specific disorders (cis-diagnostic risk), in regulatory elements in immune cells. We used three independent epigenetic datasets representing multiple organ systems and immune cell subsets. Trans-diagnostic and cis-diagnostic risk variants (for schizophrenia and depression) were enriched at epigenetically active sites in brain tissues and in lymphoid cells (T, B and NK cells), especially stimulated CD4+ T cells. There was no evidence for enrichment of either trans-risk or cis-risk variants for schizophrenia or depression in myeloid cells. This suggests a model where stimuli (e.g. infection) activate T cells to unmask the effects of psychiatric risk variants, contributing to the pathogenesis of mental health disorders.

## Introduction

Clinical diagnosis for mental health disorders is comprised of multiple, categorically distinct clinical syndromes such as schizophrenia, major depressive disorder (MDD), and bipolar disorder. However, symptoms overlap between some different psychiatric diagnoses, and comparative investigations of psychiatric disorders have revealed both shared and specific genetic 1-3 and environmental risk factors 4,5, and brain transcriptomic profiles 6. These and other data support a general predisposition to psychopathology or ‘p factor’ which captures an individual’s likelihood of developing any psychiatric disorder 7. Thus, some genetic and environmental risks operate trans-diagnostically across multiple psychiatric syndromes, rather than being cis-diagnostically aligned to a specific syndrome, as would be expected if each disorder was a biologically discrete disease entity.

Immune system abnormalities have been observed in case-control studies of many psychiatric disorders, including schizophrenia 8,9, MDD 10, bipolar disorder 11, autism spectrum disorder (ASD) 12, and attention deficit hyperactivity disorder (ADHD) 13. Among the most consistently reported findings, across multiple disorders, are increased C-reactive protein (CRP) 14, increased pro-inflammatory cytokines 14,15, increased white blood cell counts in both myeloid and lymphoid lineages 16-22, and inflammasome activation 23-25. Moreover, environmental exposures that elicit an immune response are risk factors for multiple psychiatric disorders, including *in utero* or parental infections 26,27, childhood and adult infections 28-31, childhood adversity 32, and acute or chronic stress 33. On this basis, it is conceivable that the immune system could be implicated in the pathogenesis of psychiatric disorders; but the direct evidence for a causal role of immune mechanisms is limited. Longitudinal studies have shown that immune dysregulation can be detected prior to onset of psychiatric disorder 34; but this could reflect the coincident effects of risk factors for psychiatric disorder, such as high body mass index (BMI), on both the immune system and the brain. Since germline genetic variants cannot be the consequence of disease, sequence variation associated with a disorder (or disorders) could shed light on the immune processes or cells likely to cause mental health symptoms. There is already some genetic evidence that psychiatric risk is mediated by the immune system. Polygenic risk scores (PRS) for depression, bipolar disorder and schizophrenia are associated with increased lymphocyte counts 35. Immune and psychiatric disorders are genetically correlated 36,37. Pathway analysis of genes trans-diagnostically associated with schizophrenia, bipolar disorder and MDD implicated neuronal, histone and immune pathways 38; although a larger trans-diagnostic analysis did not implicate immune cells or pathways 1.

Most genetic variants associated with psychiatric risk are in non-coding regions of the genome, likely exerting their effects by altering the activity of regulatory elements 39 such as promoters or enhancers; and enhancers can be linearly distant (>10 kilobases) from the genes they regulate 40. Some regulatory elements control gene expression in multiple tissues, but others are specific to particular tissues, or particular cell states. For example, some enhancers are active in stimulated but not resting immune cells 41-43. The locations and activity status of putative enhancers and promoters in a given tissue can be identified through characteristic epigenetic modifications, such as histone modifications.

Epigenetic mechanisms have long been thought to be important in psychiatry, especially in mediating gene-environment interactions 44,45. Epigenetic data from brain tissues have been extensively used to investigate the brain cell types and regions implicated by psychiatric risk variants 46-48, by testing whether risk variants tend to be concentrated, or “enriched”, in regions of the genome that are active in a given tissue. However, the enrichment of psychiatric risk variants in immune cell subsets has not been extensively explored. Studies to date have tended to use functional information from whole blood or immune organs, which obscures and dilutes possible effects in the myeloid and lymphoid immune cell subsets comprising these samples. There is some evidence of enrichment of risk variants for bipolar disorder in genes characteristic of neutrophils, T cells and haematopoietic stem cells; and for schizophrenia at genes in T cells and chromatin marks in T and B cells49. To our knowledge, no studies have demonstrated enrichment of trans-diagnostic risk, or of cis-risk for MDD or ASD, in any immune cell type 1,50-52, or tested if immune cell enrichment is independent of brain tissue enrichment (rather than simply due to coincidental overlap of active genomic regions in brain and immune system cells).

We hypothesized that some genetic risk variants for psychiatric disorders act via their effects on regulatory elements in specific immune cell subsets, thus potentially modulating the response of these cells to infections and other environmental stimuli. We further hypothesized that some of these immunogenetic mechanisms may represent a common pathogenic pathway to multiple psychiatric disorders. To test these hypotheses, we integrated data on common genetic variants associated with trans- and cis-diagnostic risks for psychiatric disorder(s) with data on epigenetically active genomic regions in multiple human cell and tissue types. More formally, we tested the null hypothesis that a given set of risk variants was not co-located with tissue-specific marks of epigenetic activation more frequently than expected by chance in each of multiple tissues (Roadmap/ENCODE 53,54), in 19 sorted immune cell subsets (BLUEPRINT 55), and in *ex vivo* stimulated naïve and memory CD4+ T cells and macrophages (Soskic dataset 42). To contextualise our results, we conducted parallel analyses of three “positive control” disorders: Alzheimer’s disease, a brain disorder for which genetic risk has been associated with myeloid immune cells 56; rheumatoid arthritis, a canonical adaptive autoimmune disorder; and body mass index (BMI), a common comorbidity which may contribute to observed immune abnormalities in psychiatric disorders 57.

## Results

### Trans-diagnostic psychiatric risk is enriched at active chromatin states in T cells

For trans-diagnostic risk of having any one of 8 major psychiatric disorders, we tested for enrichment of genetic risk at active regulatory elements in 88 cells or tissues from the Roadmap consortium, using stratified linkage disequilibrium score regression (s-LDSC). S-LDSC is used to test whether SNP-heritability for a disorder is concentrated or enriched in a genomic annotation58. To generate a single binary annotation of active regulatory elements for each tissue, we combined annotations for active promoters and enhancers based on histone marks (see **Methods**). We found that three main tissue classes were significantly enriched for trans-diagnostic risk variants at regulatory elements, following correction for multiple comparisons: multiple adult and fetal brain regions; T cells; and pancreatic islets (**Figure 1A, S1, Table S4**).

![Figure 1](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/08/2021.08.04.21261606/F1.medium.gif)

[Figure 1](http://medrxiv.org/content/early/2021/08/08/2021.08.04.21261606/F1)

Figure 1 Trans-diagnostic risk enrichment at epigenetically active sites in brain tissue and, independently, in T cells.
(A) Enrichment of trans-diagnostic risk at active regulatory elements (active promoters and enhancers) in 88 tissues from the Roadmap epigenomics consortium. *P*-values are shown for the results of stratified linkage disequilibrium score regression (s-LDSC) analysis, taking the union of active elements in a given cell type as the annotation of interest (see **Methods**). The *P*-values from s-LDSC were used to test the null hypotheses that risk variants were not co-located with epigenetically activated sites more frequently than expected by chance, using the Benjamini-Hochberg false discovery rate to correct for multiple tests across N=88 tissues, *P*FDR < 0.05 (purple). Tissues with nominally significant enrichment (*P* < 0.05, green) are also shown for context. Only tissue groups with at least nominally significant results are shown; for non-significant results in other tissue groups (gastrointestinal; skin; fat; muscle; reproductive; vascular) see **Figure S1**. PFC, prefrontal cortex; HSC, hematopoietic stem cells. (B) Brain-conditioned analysis: right-pointing bars show repeat of the analysis in **Figure 1A**, including the active regulatory annotations for all 10 significantly enriched brain regions as additional terms in each of the s-LDSC models for all other cell types. Left-pointing bars show the original results, as in **Figure 1A**. Following conditional analysis, there was no significant decrease in enrichment significance for any immune tissues (*P* > 0.05). Only immune cell results are shown; see **Figure S2B** for conditional results for other tissues. (C) Enrichment of trans-diagnostic risk in immune cell enhancers, genic enhancers and active promoters in Roadmap immune tissues. Tile size indicates s-LDSC significance: large tiles show results significant at *P*FDR < 0.05, using the Benjamini-Hochberg procedure to correct for the 78 annotations tested; mid-sized tiles show nominally significant results (*P* < 0.05) for context. Tile fill indicates the *P*-value rank within each annotation across cell types.

In the central nervous system (CNS), trans-risk variants were most strongly enriched at regulatory elements in fetal brain tissue samples. There was also significant enrichment (*P*FDR < 0.05) at active regulatory elements in brain structures previously reported as abnormal in neuroimaging studies of psychiatric disorders: dorsolateral prefrontal cortex, angular gyrus, inferior temporal lobe, anterior caudate, cingulate gyrus, hippocampus and substantia nigra. In contrast, trans-risk was not enriched at epigenetically active sites in astroglia, consistent with previous analyses demonstrating that genes mapped to trans-risk genetic variants are mainly expressed by neurons1 rather than glia or non-neuronal cells in the brain.

In the immune system, trans-risk variants were significantly enriched (*P*FDR < 0.05) at epigenetically active genomic sites in multiple T cell subsets, including cytotoxic, helper and regulatory T cells in adult blood and T cells in cord blood. Conversely, there was no enrichment (*P* > 0.05) of trans-risk in myeloid cells (monocytes, neutrophils).

Many regulatory elements are common to multiple tissues, so we reasoned that this pattern of CNS and immune system enrichment for trans-risk variants could be driven by coincidental overlap of brain and T cell active elements. In this case, the genetic risk would be theoretically expected to have pathogenic effect primarily by its modulation of epigenetically active sites in the brain, with no clearly independent pathogenic role mediated by T cells. We therefore repeated the s-LDSC analysis but included the active annotations for all 10 significantly enriched brain regions as extra terms in the s-LDSC models for every other cell type. In this brain-conditioned analysis, both helper and cytotoxic T cells remained strongly enriched for trans-diagnostic genetic risk (**Figure 1B**), while pancreatic islets did not (**Figure S2B**). For enriched immune tissues in the original analysis (at *P*FDR < 0.05), none showed significantly decreased enrichment following brain-conditioned analysis (two-sample *Z*-test, *P* > 0.05). Conversely, including the annotation for male fetal brain (the brain tissue showing strongest enrichment) as an extra term in s-LDSC models significantly reduced trans-risk enrichment in all other brain regions (*Z*-test *P* < 0.05) except the substantia nigra (*P* = 0.09) and hippocampus (*P* = 0.07), reflecting some overlap of active elements between different brain regions at different developmental phases, and validating our statistical approach (**Figure S2A, Table S3**). We showed the same effect for female fetal brain, the second most strongly enriched brain tissue (**Figure S2A**), excluding a potential effect of sex differences in brain development.

The global active annotation used as a binary marker of epigenetic activation combines trans-risk enrichment at three different classes of regulatory elements: active promoters, genic enhancers (enhancers found in gene bodies), and non-genic enhancers. To identify which classes were most enriched for trans-risk, we tested each class separately and found that the enrichment of trans-risk observed in terms of the global active annotation in T cells was not driven by a single class of regulatory element: there was enrichment of trans-risk at both active promoters (*P*FDR < 0.05) and enhancers (*P*FDR < 0.05 for genic enhancers) (**Figure 1C**).

### Cis-diagnostic risk is enriched at active chromatin states in T cells

Using data from the Roadmap Epigenomics Consortium, we next investigated the enrichment of cis-diagnostic risk variants at epigenetically active sites in brain tissues and immune cells for each of 5 mental health or neurodevelopmental disorders (schizophrenia, bipolar disorder, MDD, autism and ADHD) and each of 3 positive control disorders (Alzheimer’s disease, obesity [BMI], and rheumatoid arthritis).

In the CNS, cis-risks for adult-onset mental health disorders (schizophrenia, bipolar disorder, MDD) were enriched in multiple fetal and adult brain tissues, and cis-risks for child mental health or neurodevelopmental disorders (autism, ADHD) were enriched more selectively in fetal brain tissue. Cis-risk for obesity (BMI) was also enriched for active sites across multiple fetal and adult brain tissues; but cis-risk for Alzheimer’s disease was only (nominally) significantly enriched in hippocampus; and cis-risk for rheumatoid arthritis was not enriched in any brain tissue (**Figure 2A, Table S4**).

![Figure 2](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/08/2021.08.04.21261606/F2.medium.gif)

[Figure 2](http://medrxiv.org/content/early/2021/08/08/2021.08.04.21261606/F2)

Figure 2 Cis-diagnostic risk enrichment at epigenetically activated sites in adult and fetal brain tissue and immune cells for 8 specific disorders.
For each of 5 mental health disorders (schizophrenia, bipolar disorder, major depressive disorder [MDD], autism, and attention deficit-hyperactivity disorder [ADHD]), and for each of 3 positive control disorders (obesity, Alzheimer’s disease and rheumatoid arthritis), enrichment of cis-risk variants at active regulatory elements (active promoters and enhancers) was tested in (A) 10 brain tissue samples (3 fetal) and (B) 26 immune cell classes (3 fetal) 54. *P*-values are shown for the results of stratified linkage disequilibrium score regression (s-LDSC) analysis, taking the union of active elements in a given cell type as the annotation of interest. Tile size, from large to small, indicates P-value thresholds from *P*FDR < 0.05 (significant after Benjamini-Hochberg correction for all 88 tissues tested, including those not shown here), through *P* < 0.05 (nominally significant), to *P* ≥ 0.05 (not significant). Tile fill indicates the *P*-value rank within each disorder across all cells/tissues to facilitate comparisons across results from differently-powered genetic association studies.

In the immune system, cis-risks for schizophrenia, bipolar disorder, MDD and autism were enriched at globally activated sites in one or more T cell subsets (mostly at *P* < 0.05; **Figure 2B, Table S4**), with signal driven by both enhancers and promoters (**Figure S3A**). Cis-risk for rheumatoid arthritis was strongly enriched at globally active sites in multiple immune cell subsets; cis-risk for Alzheimer’s disease was significantly enriched in myeloid cells and B cells 56,59; and cis-risk for BMI was only enriched in one T cell class at *P* < 0.05 (**Figure 2B**). The statistical significance of these results depends partly on the sample size of the underlying GWAS and the heritability and polygenicity of the disorder (factors influencing power, and captured by the SNP-based heritability *Z*-score)58; but also on the strength of functional enrichment of the phenotype in that annotation. For the two most enriched immune and brain annotations (naïve cytotoxic and helper T cells; fetal male and female brain), we tested the correlation between heritability *Z*-score and functional enrichment *Z*-score across the 9 disorders included in this study (**Figure S3B**). We found a strong relationship between heritability *Z*-score and brain enrichment (fetal male brain: Spearman’s ρ = 0.87, *P* = 0.005; fetal female brain: ρ = 0.87, *P* = 0.005), but no correlation between heritability *Z*-score and immune enrichment (cytotoxic T cells: Spearman’s ρ = 0, *P =* 1, helper T cells: ρ = 0.03, *P =* 0.9), suggesting that GWAS power was not driving these results.

### Trans- and cis-risk variants are enriched at active enhancers/promoters in lymphoid cells: BLUEPRINT data

To assess the generalizability of these results in an independent dataset, we tested for enrichment of trans- and cis-risk variants at active enhancer/promoter marks (H3K27ac) in sorted immune cell subsets from the BLUEPRINT consortium 55, using the CHEERS algorithm 42 to distinguish effects in epigenetically similar cell types (see **Methods**). We replicated key prior findings in the Roadmap data, i.e., trans-risk and cis-risk for schizophrenia and depression were significantly enriched at epigenetically active sites in lymphoid cells; but not myeloid cells (**Figure 3A**). This was not due to problems with the myeloid data, as we detected the expected enrichment of Alzheimer’s Disease risk variants in macrophages **(Figure 3A**). As well as T cell enrichment, we also find enrichment of trans-risk and cis-risk for schizophrenia and (especially) depression in B cells, as well as enrichment of trans-risk and cis-risk for schizophrenia in NK cells. For ADHD and bipolar disorder (less well-powered GWAS studies with fewer independent significant loci available for analysis, see **Table S1**), no cell types were enriched at *P*FDR < 0.05 (**Figure S4A**). Despite both schizophrenia and depression showing strong lymphoid enrichment, the specific histone peaks overlapped by risk variants for these disorders were not generally shared between them (**Figure 3B, S4B**). This indicates that cis-risks for these two disorders were convergently enriched at a cellular level but distinct at the level of specific regulatory elements. It was also notable that cis-risk variants for obesity overlapped with a set of H3K27ac sites that was largely disjoint with the sets of regulatory elements overlapping with cis-risk variants for psychiatric disorders (**Figure 3B**).

![Figure 3](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/08/2021.08.04.21261606/F3.medium.gif)

[Figure 3](http://medrxiv.org/content/early/2021/08/08/2021.08.04.21261606/F3)

Figure 3 Trans- and cis-diagnostic risk variant enrichment at histone-acetylated marks on adult immune cells in the BLUEPRINT dataset.
(a) Bar plots show enrichment of genetic risk for each disorder at active promoters/enhancers (H3K27ac marks) in unstimulated, sorted immune cells. CHEERS was used to detect enrichment of risk loci at cell-type specific H3K27ac peaks (see **Methods**). The dotted black line marks the nominal significance threshold, *P* < 0.05; the solid black line marks the Bonferroni-corrected significance threshold, *P*FDR < 0.05. Note differing x-axis scales. See **Figure S4A** for ADHD and bipolar disorder results (non-significant after FDR control). (b) Venn diagram shows counts of variant-peak overlaps shared between disorders and unique to each disorder (each peak is only counted once even if overlapping multiple variants). For an upset plot of peak overlaps across all disorders, see **Figure S4B**.

### Trans- and cis-diagnostic risk variants are enriched at histone-acetylated sites in stimulated T cells: Soskic immune stimulation dataset

Given that risk of mental health disorders is affected by both genetic variation and environmental factors, we reasoned that trans- and cis-risk variants could be most significantly enriched at sites that were epigenetically activated in immune cells stimulated by cytokines (mimicking environmental insults) towards different activated cell fates. To investigate this hypothesis, and to assess the robustness of our principal findings in a third independent dataset, we tested whether trans- and cis-risks were enriched at regulatory elements (H3K27ac marks) active during immune cell activation, using a dataset of human naïve and memory CD4+ T (helper) cells and macrophages stimulated *ex vivo* in the presence of 13 cytokine combinations. The chromatin activity was assessed at early and late timepoints after exposure to cytokine stimulations (16h and 5 days for T cells and 6h and 24h for macrophages), as well as in unstimulated cells 42. Both trans-diagnostic risk variants, and cis-risk variants for MDD were most significantly enriched in memory T helper cells at day 5 following T cell stimulation with anti-CD3/anti-CD28 beads that mimic activation occurring with T cell receptor-crosslinking; trans-risk variants and cis-risk variants for schizophrenia were also significantly enriched in memory T helper cells at 16 h and in naïve T helper cells at day 5 only (**Figure 4A**). The loci of histone acetylation that overlapped with cis-risk variants for MDD in late-activated memory T cells were almost completely disjoint with the set of epigenetic marks that overlapped with cis-risk variants for schizophrenia in late-activated memory T cells (**Figure S6A**), again demonstrating convergence of cis-risk enrichment at the immune cell subset level, but divergence at the molecular level of specific regulatory elements. Enrichment was generally much greater for stimulated than unstimulated T cells, with smaller differences in the magnitude of enrichment between different cytokine stimulation conditions (**Figure 4A**). For 9 of the 10 cytokine conditions (all except Th17-cytokine polarizing condition), trans-risk enrichment was significantly greater (*Z*-test *P* < 0.05) in stimulated compared to unstimulated late-activated memory T cells.

![Figure 4](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/08/2021.08.04.21261606/F4.medium.gif)

[Figure 4](http://medrxiv.org/content/early/2021/08/08/2021.08.04.21261606/F4)

Figure 4 Trans- and cis-risk variant enrichment at histone-acetylated marks on experimentally stimulated immune cells in the Soskic immune stimulation dataset.
(A) Bar plots show enrichment of genetic risk for each condition at active promoters/enhancers (H3K27ac marks) in sorted and unstimulated or *ex vivo* stimulated immune cell classes: macrophages, naïve CD4+ (helper) T cells and memory CD4+ T cells, assayed at both early and late timepoints after stimulation with one of several different cytokine cocktails promoting differentiation to different T cell states (as shown in row labels). CHEERS was used to detect enrichment of risk loci at cell-type specific H3K27ac peaks (see **Methods**). The dotted black line marks the nominal significance threshold, *P* < 0.05; the solid black line marks the Bonferroni-corrected significance threshold, *P*Bonferroni < 0.05. Note differing x-axis scales. Results for other disorders are shown in **Figure S5**. (B) Venn diagrams show counts of variant-peak overlaps shared between disorders and unique to each disorder. For an upset plot of peak overlaps across all disorders, see **Figure S6B**. (C) All Soskic immune stimulation dataset peaks overlapped by risk variants for major depressive disorder. Each row corresponds to an H3K27ac peak overlapping a risk variant for MDD; each column corresponds to a different cytokine-induced cell state, ordered and colored as in **Figure 4A** (see legend). The blue fill shade represents how specific each peak is to each cell state (specificity rank of each peak normalized to the mean specificity rank of all peaks). Only 9 of the 108 MDD-associated H3K27ac immune peaks also overlap BMI risk variants. (d) For peaks which were both highly specific to T cells (including both unstimulated and stimulated cells) and overlapped by trans-risk variants, nearest genes were identified and tested for enrichment for curated biological pathways (GO and Reactome) using a hypergeometric test. Only the 10 most significant pathways are shown (all *P*FDR < 0.05). Fill colour indicates gene ratio (number of test genes in the pathway / total number of test genes). See **Figure S7** for results for cis-diagnostic risks.

As in the two prior independent datasets, there was no enrichment for trans-or cis-risk of psychiatric disorders at epigenetically activated sites in myeloid cells, either stimulated or unstimulated, with the exception of enrichment of bipolar disorder risk in IL26-stimulated macrophages (**Figure S5**). Cis-risk variants for obesity were enriched in unstimulated and stimulated T cell states (**Figure 4**), but only 9 of the 108 depression-associated H3K27ac peaks also overlapped with BMI risk variants (**Figure 4C, S6B**), indicating that cis-risks for obesity and depression were enriched at distinct regulatory elements in the same cell subsets.

For disorders showing enrichment in T cells, we performed pathway analysis (overrepresentation analysis) for those genes overlapping or nearest to the T-cell specific histone acetylation peaks overlapped by risk variants. Trans-risk (**Figure 4D**) and cis-risk for schizophrenia (**Figure S7**) showed enrichment of pathways including epigenetic regulation, pre-notch processing and estrogen-dependent gene expression in T cells. In contrast, rheumatoid arthritis showed enrichment of lymphoid cell differentiation, activation, and response to antigenic stimulus (**Figure S7**).

## Discussion

We examined the enrichment of genetic risk variants for psychiatric disorders at epigenetically activated regulatory sites across multiple tissues. Trans-diagnostic risk variants, commonly associated with multiple mental health and neurodevelopmental disorders, were significantly enriched at active regulatory sites in several adult and fetal brain tissue samples. Strikingly, we also found that trans-diagnostic risk variants were significantly enriched at an independent set of regulatory elements in peripheral blood lymphoid cells (but were not enriched in myeloid cells). Our key results – enrichment of trans-risk in T cells and lack of enrichment in myeloid cells – were statistically robust to multiple comparisons and replicated in three independent datasets. Other lymphoid cells (for which fewer datasets were available) are likely also implicated in pathogenesis, as we also found enrichment of trans-risk in B cells and NK cells.

Further investigation of cis-diagnostic risk variants, specifically associated with one of 5 mental health or neurodevelopmental disorders, confirmed significant enrichment of genetic risks for schizophrenia and major depressive disorder at active promoters and enhancers in both brain tissue and peripheral lymphoid cells (but not myeloid cells). Epigenetically activated sites in T cells, especially cytokine-stimulated CD4+ T cells, were most consistently and significantly enriched for sequence variants associated with schizophrenia or MDD; however, the active regulatory elements overlapped by these cis-diagnostic variants were specific to each disorder. This suggests convergence of risk for schizophrenia and depression at a cellular level in the immune system, i.e. activated T cells, and raises questions about how the involvement of different specific regulatory elements in these two disorders might relate to the different phenotypic presentations of schizophrenia and depression. We also found strong enrichment of risk for depression in both naïve and memory B cells. To our knowledge, this is the first demonstration of enrichment of genetic risk for MDD at epigenetically active sites in lymphoid cells (or indeed any immune cell type).

The cis-diagnostic enrichment results for schizophrenia and MDD were statistically robust to multiple comparisons and in clear contrast to the comparable results for 3 positive control disorders. Cis-risks for Alzheimer’s disease were significantly enriched at epigenetically activated sites in myeloid cells (but not lymphoid cells); cis-risks for rheumatoid arthritis were enriched at active sites in myeloid and lymphoid cells (but not brain tissue); and cis-risks for obesity (BMI) were enriched at active sites in brain tissue and (in some analyses) in immune cells, but with effects on regulatory elements distinct from those implicated by psychiatric disorders.

On this basis, we propose that genetic variants associated with increased risk for psychiatric disorders are likely to interact with epigenetic activation of specific and distinct regulatory elements in both the central nervous system and the adaptive immune system. This pathogenic model immediately raises three key questions. What environmental exposures cause epigenetic modification at risk-enriched sites in T cells? How could atypical T cell phenotypes, arising from this gene-by-environment interaction, cause changes in the CNS that are ultimately manifest as mental health or neurodevelopmental disorders? What is the antigen presenting cell which activates atypical CD4+ T cells (our data suggest this may be B cells)?

Infection is the most likely environmental stimulus to induce epigenetic activation in the immune system. And there is abundant epidemiological evidence that fetal and post-natal infections increase the risk for multiple psychiatric disorders 26,27,30,31. However, the immune mechanisms by which early-life infection predisposes to later psychiatric symptoms are not known. We know that fetal or childhood infections can cause long-term changes in adaptive immune cell phenotypes, including T cell memory of antigens and B cell production of antibodies, that are crucial to development of adult immunity. Thus, it is conceivable that the epigenetically activated sites enriched for trans- and cis-risks in T cells and memory B cells in these data were “marked” by exposure to infection or inflammation; and that genetic risk variants modulate the infection-induced activation of regulatory elements, leading to atypical T or B cell phenotypes following infection in people at genetic risk of psychiatric disorder. There is already some epidemiological evidence for gene-by-environment interactions between infection and risk variants for schizophrenia 60-62 and MDD 63. Many aspects of our data are compatible with this concept. For example, our finding that trans- and cis-diagnostic risk variants are enriched at sites epigenetically activated by delayed T cell responses to a wide range of pro-inflammatory cytokine stimuli seems consistent with the epidemiological finding that increased risk of multiple psychiatric disorders is found following a wide range of different infections 12,28,29,64,65.

Atypical T cell phenotypes could conceivably have effects on the brain by at least two broad routes: via stimulus-driven T cell activation and via developmental pathways (**Figure 5**). By this account, the (epi)genetic root cause of mental health-related inflammation is located in the adaptive immune system, rather than the innate immune system. Atypical T cells are supposed to impact on neuronal function via the well-documented adverse effects of soluble inflammatory mediators 66; via contact-dependent mechanisms; or via their effects on other immune or non-immune cells which in turn affect neurons. In contrast with autoimmune diseases, which tend to show greatest enrichment in early T cell activation states 42, the strongest enrichment for psychiatric risk variants was in T cells, especially late-activated memory CD4+ T cells, and memory B cells. This may reflect abnormalities in the resolution (rather than onset) of immune responses, potentially leading to chronic, low-grade peripheral inflammation seen in many psychiatric disorders.

![Figure 5](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2021/08/08/2021.08.04.21261606/F5.medium.gif)

[Figure 5](http://medrxiv.org/content/early/2021/08/08/2021.08.04.21261606/F5)

Figure 5 Schematic of potential pathogenic pathways by which genetic risk variants enriched at epigenetically active sites in T cells could lead to neuronal changes and ultimately psychiatric disorders.
Atypical stimulus-driven activation pathway (light blue boxes): Infection or other stressors can induce epigenetic activation of regulatory elements in T cells that are enriched for trans-or cis-diagnostic risk variants, proximally causing atypical T cell phenotypes, and distally causing increased inflammatory activation of innate immune (myeloid) cells in the periphery and CNS. Atypical activation of T cells resident in the meninges or brain, or trafficking into the meninges and brain from the periphery, could directly have adverse effects on neuronal function. Developmental pathway (light green box): T cells typically control microglial pruning of neuronal synapses as part of normative brain developmental programs in childhood and adolescence. Atypical T cells, potentially induced by infection or stress in genetically-at-risk individuals, could promote atypical microglial pruning of synapses, contributing to the formation of disconnected networks or circuits in the adult brain.

The developmental pathway follows growing evidence that T cells – even in the absence of infection – have an important role in controlling the activity of microglia, the brain’s most numerous resident myeloid cells, in phagocytosing synaptic terminals and neurites as part of normative neurodevelopmental programs of synaptic pruning that are ongoing from early childhood throughout adolescence. Thus atypical T cells in the meninges or brain could lead to atypical synaptic pruning 67,68, mediated by innate immune cells (microglia), with long-term consequences for maturation of large-scale brain networks of inter-connected neurons. There is evidence already from neuroimaging studies that schizophrenia and other psychiatric disorders are associated with disrupted development of connectivity in cortico-subcortical circuits 69. It will be important in future studies to investigate more directly the hypothesis that genetic risks for psychiatric disorders can result in development of dysconnected brain networks via immune as well as neuronal mechanisms.

Perhaps surprisingly, given the prior focus on innate immune abnormalities associated with psychiatric disorders, epigenetically active sites in myeloid cells were not significantly enriched for trans-risk variants or for cis-risk variants for schizophrenia or MDD. What does this mean for the pathogenic role of myeloid cells in these disorders? It is possible that genetic risk variants are indeed enriched at epigenetically active sites in myeloid cells, but only under stimulation conditions that were not represented in the three datasets we analysed. Alternatively, it could be that activation of myeloid cells and associated release of CRP and pro-inflammatory cytokines, although commonly measured as immune biomarkers in previous studies of psychiatric disorders, are downstream consequences of an atypical helper or regulatory T cell response to infection or stress. This would be consistent with the well-known role of T cells in coordinating and modulating innate immune function 70-72, including via inflammasome inhibition 73. Thus atypical T cell phenotypes could promote or sustain low-grade inflammatory states of the innate immune system in the periphery, in the meninges, and in the brain, that have been robustly associated with MDD, schizophrenia and other psychiatric disorders 74-76. A third possibility is that inflammatory responses of myeloid cells have a causal role in pathogenesis of psychiatric disorders but this is entirely driven by environmental factors and independent of inherited genetic risks.

The statistical significance of results for risk variant enrichment at epigenetically active regions reflects in part the sample size of the GWAS’s used and the heritability and polygenicity of the disorders. However, in contrast to our results for brain enrichment, we did not find any correlation between GWAS statistical power and immune enrichment. This suggests that while GWAS power can explain some differences in functional enrichment (as in brain), differences in immune cell enrichment may be more indicative of between-disorder differences in how strongly immune mechanisms contribute to pathogenesis in general, or how frequently immune mechanisms are implicated in individual patients clinically diagnosed with a specific disorder.

We focused here on European ancestry genetic results, as the currently available datasets are from European participants, but the immunogenetics of psychiatric risk should be examined in other ancestries. In addition, the epigenetic datasets used here are predominantly adult: given the role of developmental insults in psychiatric risk, it will be important to investigate genetic enrichment in immune cells sampled at different developmental stages, including adolescence.

Likewise the immune cell states considered in this analysis are the canonical states associated with infection and autoimmunity. It will also be important to explore whether genetic risk variants modulate immune cell phenotypes induced by exposure to non-infectious environmental stimuli e.g. stress, especially given that childhood adversity and other social stressors are known to profoundly increase risk for multiple psychiatric disorders.

In psychiatry, T cell phenotypes have been most investigated in schizophrenia 77, with some evidence of decreased proliferative responses to stimulation 78. Our findings motivate further investigation of T cell and B cell phenotypes across multiple psychiatric conditions, with a focus on how trans-risk affects activated T cells (e.g. stimulated ex-vivo). Functional genomic analysis of T cell subsets from patient cohorts will be particularly important to directly test for disease-associated alterations in DNA accessibility, histone modifications, enhancer-promoter interactions and gene expression. We hypothesize that these alterations will be found at those T cell peaks identified in our analysis as overlapping risk SNPs. However, there may be broader epigenetic consequences of risk variants, especially given that pathway analysis implicated epigenetic regulation processes, which may occur at sites distant from the risk variants.

In conclusion, genetic risk variants for psychiatric disorders were significantly enriched at epigenetically active enhancers/promoters in adaptive immune cells, especially stimulated T cells. This suggests a mechanistic role for T cells in mediating the interaction between environmental exposures and genetic risk variants in the pathogenesis of multiple psychiatric disorders.

### Data sharing statement

All datasets and code used for this analysis are publicly available (see **Table S2**).

## Supporting information

Supplementary material [[supplements/261606_file02.pdf]](pending:yes)

Supplementary Table 4 [[supplements/261606_file03.xlsx]](pending:yes)

## Data Availability

All datasets and code used for this analysis are publicly (see Supplementary Table 2 for full details) Github repository: https://github.com/maryellenlynall/psychimmgen2021 Processed datasets deposited in Zenodo: https://doi.org/10.5281/zenodo.5153661

[https://github.com/maryellenlynall/psychimmgen2021](https://github.com/maryellenlynall/psychimmgen2021) 

[https://doi.org/10.5281/zenodo.5153661](https://doi.org/10.5281/zenodo.5153661) 

## Declaration of interests

E.T.B. serves as a member of the scientific advisory board of Sosei Heptares and as a consultant for GlaxoSmithKline. J.G. and R.P. are paid for their editorial work for Complex Psychiatry journal. M.B.S. in the past 3 years has received consulting income from Actelion, Acadia Pharmaceuticals, Aptinyx, atai Life Sciences, Boehringer Ingelheim, Bionomics, BioXcel Therapeutics, Clexio, EmpowerPharm, Engrail Therapeutics, GW Pharmaceuticals, Janssen, Jazz Pharmaceuticals, and Roche/Genentech. M.B.S. has stock options in Oxeia Biopharmaceuticals and EpiVario. He is paid for his editorial work on *Depression and* Anxiety (Editor-in-Chief), *Biological Psychiatry* (Deputy Editor), and *UpToDate* (Co-Editor-in-Chief for Psychiatry). He has also received research support from NIH, Department of Veterans Affairs, and the Department of Defense. He is on the scientific advisory board for the Brain and Behavior Research Foundation and the Anxiety and Depression Association of America.

## Methods

### Trans- and cis-diagnostic genetic risk variants for psychiatric disorders

The primary GWAS datasets used for the identification of trans- and cis-risk genes are listed in **Table S1**. We used summary statistics from a meta-analysis of trans-diagnostic risk across 8 mental health or neurodevelopmental disorders 1: ASD, bipolar disorder, MDD, obsessive-compulsive disorder, schizophrenia, anorexia nervosa, ADHD, and Tourette syndrome. For analysis of cis-risk, i.e. risk of a specific psychiatric disorder, we separately tested 5 large primary genome-wide association studies (GWAS) of MDD79, bipolar disorder 52, schizophrenia80, autism81, and ADHD 82. For comparative purposes, we analysed GWAS results for BMI 83, Alzheimer’s disease 84, and rheumatoid arthritis 85. For all disorders except MDD, we selected the largest publicly available, predominantly-European GWAS dataset; for MDD, we used a larger recent European GWAS 79.

### Testing for enrichment of genome-wide genetic risk at regulatory elements (ROADMAP data)

Stratified linkage disequilibrium score regression (s-LDSC) can be used to test whether genetic risk is concentrated or enriched in a genomic annotation, e.g. a set of active regulatory elements in a specific cell type 58. S-LDSC hinges on the fact that the disease association statistic for a given genetic variant depends on whether that variant is linked to the disease, but also whether variants in linkage disequilibrium (LD) with that variant are linked to the disease. By testing whether variants in LD with the annotation of interest tend to have higher association scores than variants elsewhere, we can calculate an enrichment score capturing the tendency of SNP-based heritability for that disease to be co-located with that annotation 58. We used this method to test for enrichment of psychiatric risk variants at active regulatory elements in 88 cell or tissue types.

For a given tissue, CHiP-seq data assaying multiple histone marks can be integrated to segment the genome into annotations representing different functional epigenetic states, e.g., enhancers, promoters, repressed regions 54. The IDEAS algorithm 86 leverages shared features across cell types to improve this segmentation. Lacking a strong prior hypothesis about which particular regulatory elements in immune cells would be implicated by psychiatric risk, we generated a simple functional annotation of active states for each tissue in the Roadmap Epigenomics Dataset, which includes samples from all major organ systems including brain, heart, muscle, gut, adipose, skin, reproductive and immune tissues 54. Data for a given tissue or cell type sometimes come from multiple donors – as is the case for most of the brain and immune samples – and sometimes from single donors (see [https://egg2.wustl.edu/roadmap/web\_portal/meta.html](https://egg2.wustl.edu/roadmap/web\_portal/meta.html) for metadata). Immune cell subsets were magnetically sorted from live donor blood samples; brain tissues were homogenized post-mortem samples. For each Roadmap tissue/cell type, we generated a whole genome binary annotation of active regulatory elements (**Figures 1A,1B**) from IDEAS annotations based on 5 epigenetic histone marks (H3K4me3, H3K4me1, H3K36me3, H3K27me3 and H3K9me3). We combined the 6 IDEAS annotations representing active promoters and enhancers to generate a single binary annotation of active regulatory elements for each tissue. More exactly, we merged the IDEAS annotations for active transcription start sites (10\_TssA); regions flanking active TSS (8\_TssAFlnk); weak TSS (14\_TssWk); enhancers (4_Enh); genic enhancers (6_EnhG); and genic enhancers associated with transcription (17_EnhGA), following a previous definition of active states 87. We generated partitioned linkage disequilibrium (LD) scores for each tissue as recommended, using HapMap3 SNPs 58.

We then used s-LDSC to test the enrichment of psychiatric risk variants in each cell type, using a separate model for each cell type, as is standard. Summary statistics were preprocessed using the LDSC recommended script munge_sumstats.py and we performed s-LDSC for each tissue in the Roadmap dataset, using recommended settings, excluding the MHC regions, and including standard non-cell type specific annotations (baseline v1.2) (see [https://storage.googleapis.com/broad-alkesgroup-public/LDSCORE/readme\_baseline_versions](https://storage.googleapis.com/broad-alkesgroup-public/LDSCORE/readme_baseline_versions)).

To further dissect the s-LDSC results for the active annotations, we also performed s-LDSC for the 3 types of genomic element comprising the active annotation: promoters, enhancers, and genic enhancers. We generated partitioned LD scores for the promoters (10\_TssA, 8\_TssAFlnk and 14\_TssWk), enhancers (4_Enh) and genic enhancers (6_EnhG and 17_EnhGA) (**Figure 1C, Figure S3A**) then performed s-LDSC using default settings for each of these annotations in the Roadmap immune tissues. The *P*-values from s-LDSC indicate the significance of the coefficient for the cell type specific annotations. *P*-values were corrected for multiple comparisons across tissues using Benjamini-Hochberg false discovery rate. We coloured heatmaps by *P*-value rank to aid comparison across disorders or annotations which are differently powered.

To account for the possible confounding effect of shared regulatory elements between brain and immune tissues, we also performed brain-conditioned enrichment analyses: for each tissue’s s-LDSC model, we added terms for the active regulatory annotations for possibly confounding brain regions. SNP heritability *Z*-scores (heritability / standard error) and s-LDSC *Z*-scores (enrichment coefficient / standard error) were estimated using LDSC. To compare the results of the original and brain-conditional analyses, we used a one-sided two-sample *Z*-test as follows, where *β*1 is the coefficient for the annotation in the original analysis and *β*2 is the coefficient in the brain-conditional analysis. SE is the standard error of the coefficient for the original (*SE*1) or conditional (*SE*2) analysis. *Z*-scores were converted to *P*-values. ![Formula][1]</img>  

### Testing for enrichment of genetic risk variants in cell-type specific active promoters/enhancers

To compare enrichment of genetic risk at regulatory marks in different immune cell subsets, and immune cells stimulated under different conditions, we used the CHEERS algorithm 42. CHEERS quantifies the overlap of lead (independently significant) genetic risk variants with cell-specific epigenetic peaks. Crucially, CHEERS facilitates the comparison of similar cell types or conditions, which tend to have similar epigenetic profiles, by calculating peak specificity scores, indicating how specific a peak is to that cell type relative to other cell types, then quantifying cell-type enrichment as the specificity-weighted sum of overlaps of disease risk variants with these peaks. While s-LDSC leverages genome wide-information, CHEERS focuses on risk loci which meet genome-wide significance. In brief, CHEERS identifies peaks which overlap lead variants or variants in strong LD (r2<0.8) with lead variants, then calculates the mean cell type specificity score (in that cell type) of those peaks, which captures the degree of enrichment of that cell type for a given disorder. One-sided *P*-values were reported from a discrete uniform distribution (reflecting the ranking of specificity scores within each cell type) and corrected for multiple comparisons across tissues using a Bonferroni correction. To identify lead disease risk loci, all summary statistics were processed consistently: liftover to hg38, harmonization, removal of MHC region, and distance-based clumping (see below for more detail). We applied CHEERS using two human H3K27ac ChIP-seq datasets: (i) BLUEPRINT consortium data from 19 sorted unstimulated immune cells subsets (see **Figure 3A**) 55 and (ii) the Soskic immune stimulation data from sorted and *ex vivo* stimulated immune cells 42. H3K27ac marks active (rather than inactive or poised) enhancer and promoter regions 53,88. In the Soskic immune stimulation experiment, macrophages, naïve CD4+ T cells and memory CD4+ T cells were stimulated using different cytokine cocktails associated with autoimmunity or known to promote different cell fates (see **Figure 4A**). In addition, generic T cell receptor and CD28 co-stimulation signals were provided in all stimulated T cell conditions using beads coated with anti-CD3 and anti-CD28 antibodies. H3K27ac data were processed as described previously 42 to obtain cell-type specificity scores for H3K27ac peaks in each cell type or state. Here, we ran CHEERS using r2 linkage disequilibrium values taken from unrelated European individuals from the 1000 genomes dataset89, calculated using PLINK90.

To compare the enrichment between stimulated and unstimulated cell subsets, we used a one-sided, two-sample Z-test as follows, where *x*1 is the mean specificity rank for the stimulated cell subset and *x*2 is the mean specificity rank for the corresponding unstimulated cell subset. SE is the standard error of the mean and depends on the number of SNPs overlapping peaks. For a given disorder, SE is the same across different annotations, as peaks are called across the dataset as a whole. *Z*-scores were converted to *P*-values. ![Formula][2]</img>  

### Identification of independent risk loci

To identify independently significant loci for each disorder, we reprocessed all summary statistics consistently. Given the lack of well-matched linkage disequilibrium data for the populations underlying these studies, we aimed to conservatively identify independent lead variants without using LD information or conditional analysis within loci. We first lifted over the summary statistics (autosomal chromosomes only) and harmonized variants to the reference strand using the EBI summary statistics snakemake pipeline ([https://github.com/EBISPOT/gwas-sumstats-harmoniser](https://github.com/EBISPOT/gwas-sumstats-harmoniser)). Alleles with a minor allele count <10 were filtered out; where minor allele counts were not available, these were imputed from GnomAD v2.1.1 91 European frequencies lifted over to GRCh38. We then filtered all summary statistics to those variants also present at minor allele frequency > 0.01 in 1000 genomes phase 3 (unrelated European participants) called against GRCh38 89. To find independently significant lead loci, we used the Open Targets genetics finemapping pipeline ([https://github.com/opentargets/genetics-finemapping](https://github.com/opentargets/genetics-finemapping)) to filter summary statistics to variants with *P* < 5 × 10−8 (excluding MHC region chr6:28510120-33480577) and performed distance-based clumping of significant variants with a clumping distance of ±500kb. The number of lead SNPs identified for each disorder is shown in **Table S1**. ASD was excluded from downstream CHEERS analysis as only two significant loci were detected.

### Over-representation analysis

Following CHEERS analysis, to test which biological pathways were implicated in T cells, we identified those T cell-specific peaks overlapped by disease risk variants, selected the genes overlapping or nearest to those peaks, then performed pathway analysis on those genes. More specifically, we selected (for a given disorder) the union of peaks highly specific (CHEERS specificity rank >0.9) to any T cell subset in the Soskic immune stimulation dataset which were also overlapped by risk variants for that disorder. For each peak, we used the ChIPseeker seq2gene function 92 to identify those genes overlapping the peak, or with a promoter overlapping the peak, or (if no promoter overlapped the peak) the nearest gene (a maximum of 100 kilobases away). The selected genes were tested for enrichment of GO biological processes and Reactome pathways using a hypergeometric test via the clusterProfiler enricher function, with Benjamini-Hochberg correction for multiple testing 93.

## Acknowledgements

M.E.L is supported by the Medical Research Council. This work was additionally supported by the NIHR Cambridge Biomedical Research Centre. E.T.B. is supported by an NIHR Senior Investigator award. B.S. is funded by the Open Targets grant OTAR040. G.T. is supported by the Open Targets and the Wellcome Trust (grant WT206194). DFL is supported by a US Veterans Affairs Career Development Award. The authors would also like to acknowledge the crucial input of the Sanger Institute Cellular Genetics Informatics team.

*   Received August 4, 2021.
*   Revision received August 4, 2021.
*   Accepted August 8, 2021.


*   © 2021, Posted by Cold Spring Harbor Laboratory

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

## References

1.  Cross-Disorder Group of the Psychiatric Genomics Consortium. Genomic Relationships, Novel Loci, and Pleiotropic Mechanisms across Eight Psychiatric Disorders. Cell 179, 1469–1482 e1411, doi:10.1016/j.cell.2019.11.020 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2019.11.020&link_type=DOI) 

2.  Peyrot, W. J. & Price, A. L. Identifying loci with different allele frequencies among cases of eight psychiatric disorders using CC-GWAS. Nat Genet 53, 445–454, doi:10.1038/s41588-021-00787-1 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-021-00787-1&link_type=DOI) 

3.  Lee, P. H., Feng, Y. A. & Smoller, J. L. Pleiotropy and Cross-Disorder Genetics Among Psychiatric Disorders. Biol Psychiatry 89, 20–31, doi:10.1016/j.biopsych.2020.09.026 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2020.09.026&link_type=DOI) 

4.  Kessler, R. C., Davis, C. G. & Kendler, K. L. Childhood adversity and adult psychiatric disorder in the US National Comorbidity Survey. Psychol Med 27, 1101–1119, doi:10.1017/s0033291797005588 (1997).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291797005588&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=9300515&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=A1997XV80000012&link_type=ISI) 

5.  Schmitt, A., Malchow, B., Hasan, A. & Falkai, P. The impact of environmental factors in severe psychiatric disorders. Front Neurosci 8, 19, doi:10.3389/fnins.2014.00019 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fnins.2014.00019&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24574956&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

6.  Gandal, M. J. et al. Shared molecular neuropathology across major psychiatric disorders parallels polygenic overlap. Science 359, 693–697, doi:10.1126/science.aad6469 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTkvNjM3Ni82OTMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8wOC8yMDIxLjA4LjA0LjIxMjYxNjA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

7.  Caspi, A. et al. The p Factor: One General Psychopathology Factor in the Structure of Psychiatric Disorders? Clin Psychol Sci 2, 119–137, doi:10.1177/2167702613497473 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/2167702613497473&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25360393&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

8.  Miller, B. J. & Goldsmith, D. L. Towards an Immunophenotype of Schizophrenia: Progress, Potential Mechanisms, and Future Directions. Neuropsychopharmacology 42, 299–317, doi:10.1038/npp.2016.211 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2016.211&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27654215&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

9.  Khandaker, G. M. et al. Inflammation and immunity in schizophrenia: implications for pathophysiology and treatment. Lancet Psychiatry 2, 258–270, doi:10.1016/S2215-0366(14)00122-9 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2215-0366(14)00122-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26359903&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

10. Miller, A. H. & Raison, C. L. The role of inflammation in depression: from evolutionary imperative to modern treatment target. Nat Rev Immunol 16, 22–34, doi:10.1038/nri.2015.5 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nri.2015.5&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26711676&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

11. Rosenblat, J. D. & McIntyre, R. L. Bipolar Disorder and Immune Dysfunction: Epidemiological Findings, Proposed Pathophysiology and Clinical Implications. Brain Sci 7, doi:10.3390/brainsci7110144 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3390/brainsci7110144&link_type=DOI) 

12. Meltzer, A. & Van de Water, J. The Role of the Immune System in Autism Spectrum Disorder. Neuropsychopharmacology 42, 284–298, doi:10.1038/npp.2016.158 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/npp.2016.158&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27534269&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

13. Hoekstra, P. J. Attention-deficit/hyperactivity disorder: is there a connection with the immune system? Eur Child Adolesc Psychiatry 28, 601–602, doi:10.1007/s00787-019-01344-2 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s00787-019-01344-2&link_type=DOI) 

14. Yuan, N., Chen, Y., Xia, Y., Dai, J. & Liu, C. Inflammation-related biomarkers in major psychiatric disorders: a cross-disorder assessment of reproducibility and specificity in 43 meta-analyses. Transl Psychiatry 9, 233, doi:10.1038/s41398-019-0570-y (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41398-019-0570-y&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31534116&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

15. Goldsmith, D. R., Rapaport, M. H. & Miller, B. L. A meta-analysis of blood cytokine network alterations in psychiatric patients: comparisons between schizophrenia, bipolar disorder and depression. Mol Psychiatry 21, 1696–1709, doi:10.1038/mp.2016.3 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2016.3&link_type=DOI) 

16. Lynall, M. E. et al. Peripheral Blood Cell-Stratified Subgroups of Inflamed Depression. Biol Psychiatry, doi:10.1016/j.biopsych.2019.11.017 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2019.11.017&link_type=DOI) 

17. Fernandez-Egea, E. et al. Peripheral immune cell populations associated with cognitive deficits and negative symptoms of treatment-resistant schizophrenia. PLoS One 11, e0155631, doi:10.1371/journal.pone.0155631 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0155631&link_type=DOI) 

18. Grosse, L. et al. Deficiencies of the T and natural killer cell system in major depressive disorder: T regulatory cell defects are associated with inflammatory monocyte activation. Brain Behav Immun 54, 38–44, doi:10.1016/j.bbi.2015.12.003 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2015.12.003&link_type=DOI) 

19. Breunis, M. N. et al. High numbers of circulating activated T cells and raised levels of serum IL-2 receptor in bipolar disorder. Biol Psychiatry 53, 157–165, doi:10.1016/s0006-3223(02)01452-x (2003).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0006-3223(02)01452-X&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=12547472&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000180518900008&link_type=ISI) 

20. Barbosa, I. G. et al. Monocyte and lymphocyte activation in bipolar disorder: a new piece in the puzzle of immune dysfunction in mood disorders. Int J Neuropsychopharmacol 18, doi:10.1093/ijnp/pyu021 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/ijnp/pyu021&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25539506&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

21. Saresella, M. et al. An autistic endophenotype results in complex immune dysfunction in healthy siblings of autistic children. Biol Psychiatry 66, 978–984, doi:10.1016/j.biopsych.2009.06.020 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2009.06.020&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19699471&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000271497400012&link_type=ISI) 

22. Jackson, A. J. & Miller, B. L. Meta-analysis of total and differential white blood cell counts in schizophrenia. Acta Psychiatr Scand 142, 18–26, doi:10.1111/acps.13140 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1111/acps.13140&link_type=DOI) 

23. Alcocer-Gomez, E. et al. NLRP3 inflammasome is activated in mononuclear blood cells from patients with major depressive disorder. Brain Behav Immun 36, 111–117, doi:10.1016/j.bbi.2013.10.017 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2013.10.017&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24513871&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

24. Saresella, M. et al. Multiple inflammasome complexes are activated in autistic spectrum disorders. Brain Behav Immun 57, 125–133, doi:10.1016/j.bbi.2016.03.009 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.bbi.2016.03.009&link_type=DOI) 

25. Kim, H. K., Andreazza, A. C., Elmi, N., Chen, W. & Young, L. L. Nod-like receptor pyrin containing 3 (NLRP3) in the post-mortem frontal cortex from patients with bipolar disorder: A potential mediator between mitochondria and immune-activation. J Psychiatr Res 72, 43–50, doi:10.1016/j.jpsychires.2015.10.015 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jpsychires.2015.10.015&link_type=DOI) 

26. Al-Haddad, B. J. S. et al. Long-term Risk of Neuropsychiatric Disease After Exposure to Infection In Utero. JAMA Psychiatry 76, 594–602, doi:10.1001/jamapsychiatry.2019.0029 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapsychiatry.2019.0029&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30840048&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

27. Lydholm, C. N. et al. Parental Infections Before, During, and After Pregnancy as Risk Factors for Mental Disorders in Childhood and Adolescence: A Nationwide Danish Study. Biol Psychiatry 85, 317–325, doi:10.1016/j.biopsych.2018.09.013 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.biopsych.2018.09.013&link_type=DOI) 

28. Benros, M. E. et al. Autoimmune diseases and severe infections as risk factors for mood disorders: a nationwide study. JAMA Psychiatry 70, 812–820, doi:10.1001/jamapsychiatry.2013.1111 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapsychiatry.2013.1111&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23760347&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000322833600007&link_type=ISI) 

29. Benros, M. E. et al. Autoimmune diseases and severe infections as risk factors for schizophrenia: a 30-year population-based register study. Am J Psychiatry 168, 1303–1310, doi:10.1176/appi.ajp.2011.11030516 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2011.11030516&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22193673&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000297947400013&link_type=ISI) 

30. Breithaupt, L. et al. Association of Exposure to Infections in Childhood With Risk of Eating Disorders in Adolescent Girls. JAMA Psychiatry 76, 800–809, doi:10.1001/jamapsychiatry.2019.0297 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapsychiatry.2019.0297&link_type=DOI) 

31. Kohler-Forsberg, O. et al. A Nationwide Study in Denmark of the Association Between Treated Infections and the Subsequent Risk of Treated Mental Disorders in Children and Adolescents. JAMA Psychiatry 76, 271–279, doi:10.1001/jamapsychiatry.2018.3428 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapsychiatry.2018.3428&link_type=DOI) 

32. Hostinar, C. E., Lachman, M. E., Mroczek, D. K., Seeman, T. E. & Miller, G. L. Additive contributions of childhood adversity and recent stressors to inflammation at midlife: Findings from the MIDUS study. Dev Psychol 51, 1630–1644, doi:10.1037/dev0000049 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1037/dev0000049&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26389605&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

33. Rohleder, N. Stress and inflammation - The need to address the gap in the transition between acute and chronic stress effects. Psychoneuroendocrinology 105, 164–171, doi:10.1016/j.psyneuen.2019.02.021 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.psyneuen.2019.02.021&link_type=DOI) 

34. Khandaker, G. M., Pearson, R. M., Zammit, S., Lewis, G. & Jones, P. L. Association of serum interleukin 6 and C-reactive protein in childhood with depression and psychosis in young adult life: a population-based longitudinal study. JAMA Psychiatry 71, 1121–1128, doi:10.1001/jamapsychiatry.2014.1332 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamapsychiatry.2014.1332&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25133871&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

35. Sewell, M. D. E. et al. Associations Between Major Psychiatric Disorder Polygenic Risk Scores and Blood-Based Markers in UK Biobank. medRxiv, 2020.2011.2006.20227066, doi:10.1101/2020.11.06.20227066 (2020).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoibWVkcnhpdiI7czo1OiJyZXNpZCI7czoyMToiMjAyMC4xMS4wNi4yMDIyNzA2NnYxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMDgvMjAyMS4wOC4wNC4yMTI2MTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

36. Tylee, D. S. et al. Genetic correlations among psychiatric and immune-related phenotypes based on genome-wide association data. Am J Med Genet B Neuropsychiatr Genet 177, 641–657, doi:10.1002/ajmg.b.32652 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ajmg.b.32652&link_type=DOI) 

37. Pouget, J. G. et al. Cross-disorder analysis of schizophrenia and 19 immune-mediated diseases identifies shared genetic risk. Hum Mol Genet 28, 3498–3513, doi:10.1093/hmg/ddz145 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddz145&link_type=DOI) 

38. Network & Pathway Analysis Subgroup of Psychiatric Genomics, C. Psychiatric genome-wide association study analyses implicate neuronal, immune and histone pathways. Nat Neurosci 18, 199–209, doi:10.1038/nn.3922 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn.3922&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25599223&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

39. Gusev, A. et al. Partitioning heritability of regulatory and cell-type-specific variants across 11 common diseases. Am J Hum Genet 95, 535–552, doi:10.1016/j.ajhg.2014.10.004 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2014.10.004&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25439723&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

40. Won, H. et al. Chromosome conformation elucidates regulatory relationships in developing human brain. Nature 538, 523–527, doi:10.1038/nature19847 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature19847&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27760116&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

41. Chow, N. A., Jasenosky, L. D. & Goldfeld, A. L. A distal locus element mediates IFN-gamma priming of lipopolysaccharide-stimulated TNF gene expression. Cell Rep 9, 1718–1728, doi:10.1016/j.celrep.2014.11.011 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.celrep.2014.11.011&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25482561&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

42. Soskic, B. et al. Chromatin activity at GWAS loci identifies T cell states driving complex immune diseases. Nat Genet 51, 1486–1493, doi:10.1038/s41588-019-0493-9 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0493-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31548716&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

43. Alasoo, K. et al. Shared genetic effects on chromatin and gene expression indicate a role for enhancer priming in immune response. Nat Genet 50, 424–431, doi:10.1038/s41588-018-0046-7 (2018).
    
    [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%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

44. O’Donnell, K. J. & Meaney, M. L. Epigenetics, Development, and Psychopathology. Annu Rev Clin Psychol 16, 327–350, doi:10.1146/annurev-clinpsy-050718-095530 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1146/annurev-clinpsy-050718-095530&link_type=DOI) 

45. Klengel, T. & Binder, E. L. Epigenetics of Stress-Related Psychiatric Disorders and Gene x Environment Interactions. Neuron 86, 1343–1357, doi:10.1016/j.neuron.2015.05.036 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.neuron.2015.05.036&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26087162&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

46. Hauberg, M. E. et al. Common schizophrenia risk variants are enriched in open chromatin regions of human glutamatergic neurons. Nat Commun 11, 5581, doi:10.1038/s41467-020-19319-2 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-020-19319-2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=33149216&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

47. Li, M. et al. Integrative functional genomic analysis of human brain development and neuropsychiatric risks. Science 362, doi:10.1126/science.aat7615 (2018).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjE3OiIzNjIvNjQyMC9lYWF0NzYxNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDIxLzA4LzA4LzIwMjEuMDguMDQuMjEyNjE2MDYuYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

48. Girdhar, K. et al. Cell-specific histone modification maps in the human frontal lobe link schizophrenia risk to the neuronal epigenome. Nat Neurosci 21, 1126–1136, doi:10.1038/s41593-018-0187-0 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-018-0187-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30038276&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

49. Finucane, H. K. et al. Heritability enrichment of specifically expressed genes identifies disease-relevant tissues and cell types. Nat Genet 50, 621–629, doi:10.1038/s41588-018-0081-4 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0081-4&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29632380&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

50. Howard, D. M. et al. Genome-wide meta-analysis of depression identifies 102 independent variants and highlights the importance of the prefrontal brain regions. Nat Neurosci 22, 343–352, doi:10.1038/s41593-018-0326-7 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-018-0326-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30718901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

51. Alonso-Gonzalez, A., Calaza, M., Rodriguez-Fontenla, C. & Carracedo, A. Novel Gene-Based Analysis of ASD GWAS: Insight Into the Biological Role of Associated Genes. Front Genet 10, 733, doi:10.3389/fgene.2019.00733 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.3389/fgene.2019.00733&link_type=DOI) 

52. Stahl, E. A. et al. Genome-wide association study identifies 30 loci associated with bipolar disorder. Nat Genet 51, 793–803, doi:10.1038/s41588-019-0397-8 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0397-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=31043756&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

53. Consortium, E. P. An integrated encyclopedia of DNA elements in the human genome. Nature 489, 57–74, doi:10.1038/nature11247 (2012).
    
    [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%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000308347000039&link_type=ISI) 

54. Roadmap Epigenomics, C. et al. Integrative analysis of 111 reference human epigenomes. Nature 518, 317–330, doi:10.1038/nature14248 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature14248&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25693563&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

55. Chen, L. et al. Genetic Drivers of Epigenetic and Transcriptional Variation in Human Immune Cells. Cell 167, 1398–1414 e1324, doi:10.1016/j.cell.2016.10.026 (2016).
    
    [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%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

56. Novikova, G. et al. Integration of Alzheimer’s disease genetics and myeloid genomics identifies disease risk regulatory elements and genes. Nat Commun 12, 1610, doi:10.1038/s41467-021-21823-y (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-21823-y&link_type=DOI) 

57. McLaughlin, A. P. et al. The influence of comorbid depression and overweight status on peripheral inflammation and cortisol levels. Psychol Med, 1-8, doi:10.1017/S0033291721000088 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1017/S0033291721000088&link_type=DOI) 

58. Finucane, H. K. et al. Partitioning heritability by functional annotation using genome-wide association summary statistics. Nat Genet 47, 1228–1235, doi:10.1038/ng.3404 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.3404&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26414678&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

59. Kim, K. et al. Therapeutic B-cell depletion reverses progression of Alzheimer’s disease. Nat Commun 12, 2185, doi:10.1038/s41467-021-22479-4 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41467-021-22479-4&link_type=DOI) 

60. Clarke, M. C., Tanskanen, A., Huttunen, M., Whittaker, J. C. & Cannon, M. Evidence for an interaction between familial liability and prenatal exposure to infection in the causation of schizophrenia. Am J Psychiatry 166, 1025–1030, doi:10.1176/appi.ajp.2009.08010031 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2009.08010031&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19487391&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000269483600013&link_type=ISI) 

61. Demontis, D. et al. Association of GRIN1 and GRIN2A-D with schizophrenia and genetic interaction with maternal herpes simplex virus-2 infection affecting disease risk. Am J Med Genet B Neuropsychiatr Genet 156B, 913–922, doi:10.1002/ajmg.b.31234 (2011).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ajmg.b.31234&link_type=DOI) 

62. Borglum, A. D. et al. Genome-wide study of association and interaction with maternal cytomegalovirus infection suggests new schizophrenia loci. Mol Psychiatry 19, 325–333, doi:10.1038/mp.2013.2 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/mp.2013.2&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23358160&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331951600011&link_type=ISI) 

63. Ye, J. et al. Association between herpes simplex virus 1 exposure and the risk of depression in UK Biobank. Clin Transl Med 10, e108, doi:10.1002/ctm2.108 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/ctm2.108&link_type=DOI) 

64. Tioleco, N. et al. Prenatal maternal infection and risk for autism in offspring: A meta-analysis. Autism Res, doi:10.1002/aur.2499 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/aur.2499&link_type=DOI) 

65. Brown, A. S. & Meyer, U. Maternal Immune Activation and Neuropsychiatric Illness: A Translational Research Perspective. Am J Psychiatry 175, 1073–1083, doi:10.1176/appi.ajp.2018.17121311 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1176/appi.ajp.2018.17121311&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30220221&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

66. Choi, G. B. et al. The maternal interleukin-17a pathway in mice promotes autism-like phenotypes in offspring. Science 351, 933–939, doi:10.1126/science.aad0314 (2016).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzNTEvNjI3Ni85MzMiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8wOC8yMDIxLjA4LjA0LjIxMjYxNjA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

67. Pasciuto, E. et al. Microglia Require CD4 T Cells to Complete the Fetal-to-Adult Transition. Cell 182, 625–640 e624, doi:10.1016/j.cell.2020.06.026 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.cell.2020.06.026&link_type=DOI) 

68. Sekar, A. et al. Schizophrenia risk from complex variation of complement component 4. Nature 530, 177–183, doi:10.1038/nature16549 (2016).
    
    [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%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

69. Morgan, S. E. et al. Cortical patterning of abnormal morphometric similarity in psychosis is associated with brain expression of schizophrenia-related genes. Proc Natl Acad Sci U S A 116, 9604–9609, doi:10.1073/pnas.1820754116 (2019).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMToiMTE2LzE5Lzk2MDQiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8wOC8yMDIxLjA4LjA0LjIxMjYxNjA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

70. Cohen, S. B. et al. CXCR3-dependent CD4(+) T cells are required to activate inflammatory monocytes for defense against intestinal infection. PLoS Pathog 9, e1003706, doi:10.1371/journal.ppat.1003706 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.ppat.1003706&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24130498&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

71. Kim, K. D. et al. Adaptive immune cells temper initial innate responses. Nat Med 13, 1248–1252, doi:10.1038/nm1633 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nm1633&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17891146&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000249980200042&link_type=ISI) 

72. Rauch, P. J. et al. Innate response activator B cells protect against microbial sepsis. Science 335, 597–601, doi:science.1215173 [pii] 10.1126/science.1215173 (2012).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6Mzoic2NpIjtzOjU6InJlc2lkIjtzOjEyOiIzMzUvNjA2OC81OTciO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyMS8wOC8wOC8yMDIxLjA4LjA0LjIxMjYxNjA2LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

73. Guarda, G. et al. T cells dampen innate immune responses through inhibition of NLRP1 and NLRP3 inflammasomes. Nature 460, 269–273, doi:10.1038/nature08100 (2009).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature08100&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=19494813&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000267761000043&link_type=ISI) 

74. Alves de Lima, K. et al. Meningeal gammadelta T cells regulate anxiety-like behavior via IL-17a signaling in neurons. Nat Immunol 21, 1421–1429, doi:10.1038/s41590-020-0776-4 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41590-020-0776-4&link_type=DOI) 

75. Filiano, A. J. et al. Unexpected role of interferon-gamma in regulating neuronal connectivity and social behaviour. Nature 535, 425–429, doi:10.1038/nature18626 (2016).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature18626&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27409813&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

76. Ziv, Y. et al. Immune cells contribute to the maintenance of neurogenesis and spatial learning abilities in adulthood. Nat Neurosci 9, 268–275, doi:10.1038/nn1629 (2006).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nn1629&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16415867&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000234990100023&link_type=ISI) 

77. Debnath, M. Adaptive Immunity in Schizophrenia: Functional Implications of T Cells in the Etiology, Course and Treatment. J Neuroimmune Pharmacol 10, 610–619, doi:10.1007/s11481-015-9626-9 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1007/s11481-015-9626-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26162591&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

78. Craddock, R. M. et al. Altered T-cell function in schizophrenia: a cellular model to investigate molecular disease mechanisms. PLoS One 2, e692, doi:10.1371/journal.pone.0000692 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0000692&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17668069&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

79. Levey, D. F. et al. Bi-ancestral depression GWAS in the Million Veteran Program and meta-analysis in >1.2 million individuals highlight new therapeutic directions. Nat Neurosci, doi:10.1038/s41593-021-00860-2 (2021).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41593-021-00860-2&link_type=DOI) 

80. Ripke, S. et al. Genome-wide association analysis identifies 13 new risk loci for schizophrenia. Nat Genet 45, 1150–1159, doi:10.1038/ng.2742 (2013).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/ng.2742&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23974872&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

81. Grove, J. et al. Identification of common genetic risk variants for autism spectrum disorder. Nat Genet 51, 431–444, doi:10.1038/s41588-019-0344-8 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-019-0344-8&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30804558&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

82. Demontis, D. et al. Discovery of the first genome-wide significant risk loci for attention deficit/hyperactivity disorder. Nat Genet 51, 63–75, doi:10.1038/s41588-018-0269-7 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0269-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30478444&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

83. Pulit, S. L. et al. Meta-analysis of genome-wide association studies for body fat distribution in 694 649 individuals of European ancestry. Hum Mol Genet 28, 166–174, doi:10.1093/hmg/ddy327 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/hmg/ddy327&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30239722&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

84. Jansen, I. E. et al. Genome-wide meta-analysis identifies new loci and functional pathways influencing Alzheimer’s disease risk. Nat Genet 51, 404–413, doi:10.1038/s41588-018-0311-9 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41588-018-0311-9&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=30617256&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

85. Okada, Y. et al. Genetics of rheumatoid arthritis contributes to biology and drug discovery. Nature 506, 376–381, doi:10.1038/nature12873 (2014).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/nature12873&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24390342&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000331477800043&link_type=ISI) 

86. Zhang, Y. & Hardison, R. L. Accurate and reproducible functional maps in 127 human cell types via 2D genome segmentation. Nucleic Acids Res 45, 9823–9836, doi:10.1093/nar/gkx659 (2017).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/nar/gkx659&link_type=DOI) 

87. Backenroth, D. et al. FUN-LDA: A Latent Dirichlet Allocation Model for Predicting Tissue-Specific Functional Effects of Noncoding Variation: Methods and Applications. Am J Hum Genet 102, 920–942, doi:10.1016/j.ajhg.2018.03.026 (2018).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ajhg.2018.03.026&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=29727691&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

88. Creyghton, M. P. et al. Histone H3K27ac separates active from poised enhancers and predicts developmental state. Proc Natl Acad Sci U S A 107, 21931–21936, doi:10.1073/pnas.1016071107 (2010).
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NDoicG5hcyI7czo1OiJyZXNpZCI7czoxMjoiMTA3LzUwLzIxOTMxIjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjEvMDgvMDgvMjAyMS4wOC4wNC4yMTI2MTYwNi5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

89. Lowy-Gallego, E. et al. Variant calling on the GRCh38 assembly with the data from phase three of the 1000 Genomes Project. Wellcome Open Res 4, 50, doi:10.12688/wellcomeopenres.15126.2 (2019).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.12688/wellcomeopenres.15126.2&link_type=DOI) 

90. Purcell, S. et al. PLINK: a tool set for whole-genome association and population-based linkage analyses. Am J Hum Genet 81, 559–575, doi:10.1086/519795 (2007).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1086/519795&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=17701901&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

91. Karczewski, K. J. et al. The mutational constraint spectrum quantified from variation in 141,456 humans. Nature 581, 434–443, doi:10.1038/s41586-020-2308-7 (2020).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1038/s41586-020-2308-7&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=32461654&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

92. Yu, G., Wang, L. G. & He, Q. L. ChIPseeker: an R/Bioconductor package for ChIP peak annotation, comparison and visualization. Bioinformatics 31, 2382–2383, doi:10.1093/bioinformatics/btv145 (2015).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/bioinformatics/btv145&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25765347&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 

93. Yu, G., Wang, L. G., Han, Y. & He, Q. L. clusterProfiler: an R package for comparing biological themes among gene clusters. OMICS 16, 284–287, doi:10.1089/omi.2011.0118 (2012).
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1089/omi.2011.0118&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=22455463&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2021%2F08%2F08%2F2021.08.04.21261606.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000303653300007&link_type=ISI)

 [1]: /embed/graphic-6.gif
 [2]: /embed/graphic-7.gif