Transcriptomic effects of rs4845604, an IBD and allergy-associated RORC variant, in stimulated ex vivo CD4+ T cells

RORγt is an isoform of RORC, preferentially expressed in Th17 cells, that functions as a critical regulator of type 3 immunity. As murine Th17-driven inflammatory disease models were greatly diminished in RORC knock-out mice, this receptor was prioritised as an attractive therapeutic target for the treatment of several autoimmune diseases. Human genetic studies indicate a significant contributory role for RORC in several human disease conditions. Furthermore, genome-wide association studies (GWAS) report a significant association between inflammatory bowel disease (IBD) and the RORC regulatory variant rs4845604. To investigate if the rs4845604 variant may affect CD4⁺ T cell differentiation events, naïve CD4⁺ T cells were isolated from eighteen healthy subjects homozygous for the rs4845604 minor (A) or major (G) allele). Isolated cells from each subject were differentiated into distinct T cell lineages by culturing in either T cell maintenance medium or Th17 driving medium conditions for six days in the presence of an RORC inverse agonist (to prevent constitutive receptor activity) or an inactive diastereomer (control). Our proof of concept study indicated that genotype had no significant effect on the mean number of naïve CD4 T cells isolated, nor the frequency of Th1-like and Th17-like cells following six days of culture in any of the four culture conditions. Analysis of the derived RNA-seq count data identified genotype-driven transcriptional effects in each of the four culture conditions. Subsequent pathway enrichment analysis of these profiles reported perturbation of metabolic signalling networks, with the potential to affect the cellular detoxification response. This investigation reveals that rs4845604 genotype is associated with transcriptional effects in CD4⁺ T cells that may perturb immune and metabolic pathways. Most significantly, the rs4845604 GG, IBD risk associated, genotype may be associated with a differential detoxification response. This observation justifies further investigation in a larger cohort of both healthy and IBD-affected individuals.

Genomic atlas of the human plasma proteome

Although plasma proteins have important roles in biological processes and are the direct targets of many drugs, the genetic factors that control inter-individual variation in plasma protein levels are not well understood. Here we characterize the genetic architecture of the human plasma proteome in healthy blood donors from the INTERVAL study. We identify 1,927 genetic associations with 1,478 proteins, a fourfold increase on existing knowledge, including trans associations for 1,104 proteins. To understand the consequences of perturbations in plasma protein levels, we apply an integrated approach that links genetic variation with biological pathway, disease, and drug databases. We show that protein quantitative trait loci overlap with gene expression quantitative trait loci, as well as with disease-associated loci, and find evidence that protein biomarkers have causal roles in disease using Mendelian randomization analysis. By linking genetic factors to diseases via specific proteins, our analyses highlight potential therapeutic targets, opportunities for matching existing drugs with new disease indications, and potential safety concerns for drugs under development.

Exome-wide association study reveals novel psoriasis susceptibility locus at TNFSF15 and rare protective alleles in genes contributing to type I IFN signalling

Psoriasis is a common inflammatory skin disorder for which multiple genetic susceptibility loci have been identified, but few resolved to specific functional variants. In this study, we sought to identify common and rare psoriasis-associated gene-centric variation. Using exome arrays we genotyped four independent cohorts, totalling 11 861 psoriasis cases and 28 610 controls, aggregating the dataset through statistical meta-analysis. Single variant analysis detected a previously unreported risk locus at TNFSF15 (rs6478108; P = 1.50 × 10-8, OR = 1.10), and association of common protein-altering variants at 11 loci previously implicated in psoriasis susceptibility. We validate previous reports of protective low-frequency protein-altering variants within IFIH1 (encoding an innate antiviral receptor) and TYK2 (encoding a Janus kinase), in each case establishing a further series of protective rare variants (minor allele frequency < 0.01) via gene-wide aggregation testing (IFIH1: pburden = 2.53 × 10-7, OR = 0.707; TYK2: pburden = 6.17 × 10-4, OR = 0.744). Both genes play significant roles in type I interferon (IFN) production and signalling. Several of the protective rare and low-frequency variants in IFIH1 and TYK2 disrupt conserved protein domains, highlighting potential mechanisms through which their effect may be exerted.

Fine-mapping inflammatory bowel disease loci to single-variant resolution

Inflammatory bowel diseases are chronic gastrointestinal inflammatory disorders that affect millions of people worldwide. Genome-wide association studies have identified 200 inflammatory bowel disease-associated loci, but few have been conclusively resolved to specific functional variants. Here we report fine-mapping of 94 inflammatory bowel disease loci using high-density genotyping in 67,852 individuals. We pinpoint 18 associations to a single causal variant with greater than 95% certainty, and an additional 27 associations to a single variant with greater than 50% certainty. These 45 variants are significantly enriched for protein-coding changes (n = 13), direct disruption of transcription-factor binding sites (n = 3), and tissue-specific epigenetic marks (n = 10), with the last category showing enrichment in specific immune cells among associations stronger in Crohn's disease and in gut mucosa among associations stronger in ulcerative colitis. The results of this study suggest that high-resolution fine-mapping in large samples can convert many discoveries from genome-wide association studies into statistically convincing causal variants, providing a powerful substrate for experimental elucidation of disease mechanisms.

Genome-Wide Association Studies Suggest Limited Immune Gene Enrichment in Schizophrenia Compared to 5 Autoimmune Diseases

There has been intense debate over the immunological basis of schizophrenia, and the potential utility of adjunct immunotherapies. The major histocompatibility complex is consistently the most powerful region of association in genome-wide association studies (GWASs) of schizophrenia and has been interpreted as strong genetic evidence supporting the immune hypothesis. However, global pathway analyses provide inconsistent evidence of immune involvement in schizophrenia, and it remains unclear whether genetic data support an immune etiology per se. Here we empirically test the hypothesis that variation in immune genes contributes to schizophrenia. We show that there is no enrichment of immune loci outside of the MHC region in the largest genetic study of schizophrenia conducted to date, in contrast to 5 diseases of known immune origin. Among 108 regions of the genome previously associated with schizophrenia, we identify 6 immune candidates (DPP4, HSPD1, EGR1, CLU, ESAM, NFATC3) encoding proteins with alternative, nonimmune roles in the brain. While our findings do not refute evidence that has accumulated in support of the immune hypothesis, they suggest that genetically mediated alterations in immune function may not play a major role in schizophrenia susceptibility. Instead, there may be a role for pleiotropic effects of a small number of immune genes that also regulate brain development and plasticity. Whether immune alterations drive schizophrenia progression is an important question to be addressed by future research, especially in light of the growing interest in applying immunotherapies in schizophrenia.

Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci

We simultaneously investigated the genetic landscape of ankylosing spondylitis, Crohn's disease, psoriasis, primary sclerosing cholangitis and ulcerative colitis to investigate pleiotropy and the relationship between these clinically related diseases. Using high-density genotype data from more than 86,000 individuals of European ancestry, we identified 244 independent multidisease signals, including 27 new genome-wide significant susceptibility loci and 3 unreported shared risk loci. Complex pleiotropy was supported when contrasting multidisease signals with expression data sets from human, rat and mouse together with epigenetic and expressed enhancer profiles. The comorbidities among the five immune diseases were best explained by biological pleiotropy rather than heterogeneity (a subgroup of cases genetically identical to those with another disease, possibly owing to diagnostic misclassification, molecular subtypes or excessive comorbidity). In particular, the strong comorbidity between primary sclerosing cholangitis and inflammatory bowel disease is likely the result of a unique disease, which is genetically distinct from classical inflammatory bowel disease phenotypes.

Strategies for fine-mapping complex traits

Genome-wide association studies (GWAS) have identified thousands of robust and replicable genetic associations for complex disease. However, the identification of the causal variants that underlie these associations has been more difficult. This problem of fine-mapping association signals predates GWAS, but the last few years have seen a surge of studies aimed at pinpointing causal variants using both statistical evidence from large association data sets and functional annotations of genetic variants. Combining these two approaches can often determine not only the causal variant but also the target gene. Recent contributions include analyses of custom genotyping arrays, such as the Immunochip, statistical methods to identify credible sets of causal variants and the addition of functional genomic annotations for coding and non-coding variation to help prioritize variants and discern functional consequence and hence the biological basis of disease risk.

Do Genetic Susceptibility Variants Associate with Disease Severity in Early Active Rheumatoid Arthritis?

OBJECTIVE

Genetic variants affect both the development and severity of rheumatoid arthritis (RA). Recent studies have expanded the number of RA susceptibility variants. We tested the hypothesis that these associated with disease severity in a clinical trial cohort of patients with early, active RA.

METHODS

We evaluated 524 patients with RA enrolled in the Combination Anti-Rheumatic Drugs in Early RA (CARDERA) trials. We tested validated susceptibility variants - 69 single-nucleotide polymorphisms (SNP), 15 HLA-DRB1 alleles, and amino acid polymorphisms in 6 HLA molecule positions - for their associations with progression in Larsen scoring, 28-joint Disease Activity Scores, and Health Assessment Questionnaire (HAQ) scores over 2 years using linear mixed-effects and latent growth curve models.

RESULTS

HLA variants were associated with joint destruction. The *04:01 SNP (rs660895, p = 0.0003), *04:01 allele (p = 0.0002), and HLA-DRβ1 amino acids histidine at position 13 (p = 0.0005) and valine at position 11 (p = 0.0012) significantly associated with radiological progression. This association was only significant in anticitrullinated protein antibody (ACPA)-positive patients, suggesting that while their effects were not mediated by ACPA, they only predicted joint damage in ACPA-positive RA. Non-HLA variants did not associate with radiograph damage (assessed individually and cumulatively as a weighted genetic risk score). Two SNP - rs11889341 (STAT4, p = 0.0001) and rs653178 (SH2B3-PTPN11, p = 0.0004) - associated with HAQ scores over 6-24 months.

CONCLUSION

HLA susceptibility variants play an important role in determining radiological progression in early, active ACPA-positive RA. Genome-wide and HLA-wide analyses across large populations are required to better characterize the genetic architecture of radiological progression in RA.

Enhanced meta-analysis and replication studies identify five new psoriasis susceptibility loci

Psoriasis is a chronic autoimmune disease with complex genetic architecture. Previous genome-wide association studies (GWAS) and a recent meta-analysis using Immunochip data have uncovered 36 susceptibility loci. Here, we extend our previous meta-analysis of European ancestry by refined genotype calling and imputation and by the addition of 5,033 cases and 5,707 controls. The combined analysis, consisting of over 15,000 cases and 27,000 controls, identifies five new psoriasis susceptibility loci at genome-wide significance (P<5 × 10(-8)). The newly identified signals include two that reside in intergenic regions (1q31.1 and 5p13.1) and three residing near PLCL2 (3p24.3), NFKBIZ (3q12.3) and CAMK2G (10q22.2). We further demonstrate that NFKBIZ is a TRAF3IP2-dependent target of IL-17 signalling in human skin keratinocytes, thereby functionally linking two strong candidate genes. These results further integrate the genetics and immunology of psoriasis, suggesting new avenues for functional analysis and improved therapies.

Investigating the genetic variation underlying episodicity in major depressive disorder: suggestive evidence for a bipolar contribution

BACKGROUND

Highly recurrent major depressive disorder (MDD) has reportedly increased risk of shifting to bipolar disorder; high recurrence frequency has, therefore, featured as evidence of 'soft bipolarity'. We aimed to investigate the genetic underpinnings of total depressive episode count in recurrent MDD.

METHODS

Our primary sample included 1966 MDD cases with negative family history of bipolar disorder from the RADIANT studies. Total episode count was adjusted for gender, age, MDD duration, study and center before being tested for association with genotype in two separate genome-wide analyses (GWAS), in the full set and in a subset of 1364 cases with positive family history of MDD (FH+). We also calculated polygenic scores from the Psychiatric Genomics Consortium MDD and bipolar disorder studies.

RESULTS

Episodicity (especially intermediate episode counts) was an independent index of MDD familial aggregation, replicating previous reports. The GWAS produced no genome-wide significant findings. The strongest signals were detected in the full set at MAGI1 (p=5.1×10(-7)), previously associated with bipolar disorder, and in the FH+ subset at STIM1 (p=3.9×10(-6) after imputation), a calcium channel signaling gene. However, these findings failed to replicate in an independent Munich cohort. In the full set polygenic profile analyses, MDD polygenes predicted episodicity better than bipolar polygenes; however, in the FH+ subset, both polygenic scores performed similarly.

LIMITATIONS

Episode count was self-reported and, therefore, subject to recall bias.

CONCLUSIONS

Our findings lend preliminary support to the hypothesis that highly recurrent MDD with FH+ is part of a 'soft bipolar spectrum' but await replication in larger cohorts.

Germline mutations affecting the proofreading domains of POLE and POLD1 predispose to colorectal adenomas and carcinomas

Many individuals with multiple or large colorectal adenomas or early-onset colorectal cancer (CRC) have no detectable germline mutations in the known cancer predisposition genes. Using whole-genome sequencing, supplemented by linkage and association analysis, we identified specific heterozygous POLE or POLD1 germline variants in several multiple-adenoma and/or CRC cases but in no controls. The variants associated with susceptibility, POLE p.Leu424Val and POLD1 p.Ser478Asn, have high penetrance, and POLD1 mutation was also associated with endometrial cancer predisposition. The mutations map to equivalent sites in the proofreading (exonuclease) domain of DNA polymerases ɛ and δ and are predicted to cause a defect in the correction of mispaired bases inserted during DNA replication. In agreement with this prediction, the tumors from mutation carriers were microsatellite stable but tended to acquire base substitution mutations, as confirmed by yeast functional assays. Further analysis of published data showed that the recently described group of hypermutant, microsatellite-stable CRCs is likely to be caused by somatic POLE mutations affecting the exonuclease domain.

Host-microbe interactions have shaped the genetic architecture of inflammatory bowel disease

Crohn’s disease (CD) and ulcerative colitis (UC), the two common forms of inflammatory bowel disease (IBD), affect over 2.5 million people of European ancestry with rising prevalence in other populations¹. Genome-wide association studies (GWAS) and subsequent meta-analyses of CD and UC^2,3 as separate phenotypes implicated previously unsuspected mechanisms, such as autophagy⁴, in pathogenesis and showed that some IBD loci are shared with other inflammatory diseases⁵. Here we expand knowledge of relevant pathways by undertaking a meta-analysis of CD and UC genome-wide association scans, with validation of significant findings in more than 75,000 cases and controls. We identify 71 new associations, for a total of 163 IBD loci that meet genome-wide significance thresholds. Most loci contribute to both phenotypes, and both directional and balancing selection effects are evident. Many IBD loci are also implicated in other immune-mediated disorders, most notably with ankylosing spondylitis and psoriasis. We also observe striking overlap between susceptibility loci for IBD and mycobacterial infection. Gene co-expression network analysis emphasizes this relationship, with pathways shared between host responses to mycobacteria and those predisposing to IBD.

Identification of 15 new psoriasis susceptibility loci highlights the role of innate immunity

To gain further insight into the genetic architecture of psoriasis, we conducted a meta-analysis of 3 genome-wide association studies (GWAS) and 2 independent data sets genotyped on the Immunochip, including 10,588 cases and 22,806 controls. We identified 15 new susceptibility loci, increasing to 36 the number associated with psoriasis in European individuals. We also identified, using conditional analyses, five independent signals within previously known loci. The newly identified loci shared with other autoimmune diseases include candidate genes with roles in regulating T-cell function (such as RUNX3, TAGAP and STAT3). Notably, they included candidate genes whose products are involved in innate host defense, including interferon-mediated antiviral responses (DDX58), macrophage activation (ZC3H12C) and nuclear factor (NF)-κB signaling (CARD14 and CARM1). These results portend a better understanding of shared and distinctive genetic determinants of immune-mediated inflammatory disorders and emphasize the importance of the skin in innate and acquired host defense.

Conditional analysis identifies three novel major histocompatibility complex loci associated with psoriasis

Psoriasis is a common, chronic, inflammatory skin disorder. A number of genetic loci have been shown to confer risk for psoriasis. Collectively, these offer an integrated model for the inherited basis for susceptibility to psoriasis that combines altered skin barrier function together with the dysregulation of innate immune pathogen sensing and adap-tive immunity. The major histocompatibility complex (MHC) harbours the psoriasis susceptibility region which exhibits the largest effect size, driven in part by variation contained on the HLA-Cw*0602 allele. However, the resolution of the number and genomic location of potential independent risk loci are hampered by extensive linkage disequilibrium across the region. We leveraged the power of large psoriasis case and control data sets and the statistical approach of conditional analysis to identify potential further association signals distributed across the MHC. In addition to the major loci at HLA-C (P = 2.20 × 10(-236)), we observed and replicated four additional independent signals for disease association, three of which are novel. We detected evidence for association at SNPs rs2507971 (P = 6.73 × 10(-14)), rs9260313 (P = 7.93 × 10(-09)), rs66609536 (P = 3.54 × 10(-07)) and rs380924 (P = 6.24 × 10(-06)), located within the class I region of the MHC, with each observation replicated in an independent sample (P ≤ 0.01). The previously identified locus is close to MICA, the other three lie near MICB, HLA-A and HCG9 (a non-coding RNA gene). The identification of disease associations with both MICA and MICB is particularly intriguing, since each encodes an MHC class I-related protein with potent immunological function.

Conditional analysis identifies three novel major histocompatibility complex loci associated with psoriasis

Psoriasis is a common, chronic, inflammatory skin disorder. A number of genetic loci have been shown to confer risk for psoriasis. Collectively, these offer an integrated model for the inherited basis for susceptibility to psoriasis that combines altered skin barrier function together with the dysregulation of innate immune pathogen sensing and adap-tive immunity. The major histocompatibility complex (MHC) harbours the psoriasis susceptibility region which exhibits the largest effect size, driven in part by variation contained on the HLA-Cw*0602 allele. However, the resolution of the number and genomic location of potential independent risk loci are hampered by extensive linkage disequilibrium across the region. We leveraged the power of large psoriasis case and control data sets and the statistical approach of conditional analysis to identify potential further association signals distributed across the MHC. In addition to the major loci at HLA-C (P = 2.20 × 10(-236)), we observed and replicated four additional independent signals for disease association, three of which are novel. We detected evidence for association at SNPs rs2507971 (P = 6.73 × 10(-14)), rs9260313 (P = 7.93 × 10(-09)), rs66609536 (P = 3.54 × 10(-07)) and rs380924 (P = 6.24 × 10(-06)), located within the class I region of the MHC, with each observation replicated in an independent sample (P ≤ 0.01). The previously identified locus is close to MICA, the other three lie near MICB, HLA-A and HCG9 (a non-coding RNA gene). The identification of disease associations with both MICA and MICB is particularly intriguing, since each encodes an MHC class I-related protein with potent immunological function.

Refinement of the associations between risk of colorectal cancer and polymorphisms on chromosomes 1q41 and 12q13.13

In genome-wide association studies (GWASs) of colorectal cancer, we have identified two genomic regions in which pairs of tagging-single nucleotide polymorphisms (tagSNPs) are associated with disease; these comprise chromosomes 1q41 (rs6691170, rs6687758) and 12q13.13 (rs7163702, rs11169552). We investigated these regions further, aiming to determine whether they contain more than one independent association signal and/or to identify the SNPs most strongly associated with disease. Genotyping of additional sample sets at the original tagSNPs showed that, for both regions, the two tagSNPs were unlikely to identify a single haplotype on which the functional variation lay. Conversely, one of the pair of SNPs did not fully capture the association signal in each region. We therefore undertook more detailed analyses, using imputation, logistic regression, genealogical analysis using the GENECLUSTER program and haplotype analysis. In the 1q41 region, the SNP rs11118883 emerged as a strong candidate based on all these analyses, sufficient to account for the signals at both rs6691170 and rs6687758. rs11118883 lies within a region with strong evidence of transcriptional regulatory activity and has been associated with expression of PDGFRB mRNA. For 12q13.13, a complex situation was found: SNP rs7972465 showed stronger association than either rs11169552 or rs7136702, and GENECLUSTER found no good evidence for a two-SNP model. However, logistic regression and haplotype analyses supported a two-SNP model, in which a signal at the SNP rs706793 was added to that at rs11169552. Post-GWAS fine-mapping studies are challenging, but the use of multiple tools can assist in identifying candidate functional variants in at least some cases.

Colorectal cancer risk is not associated with increased levels of homozygosity in a population from the United Kingdom

Genome-wide association studies have identified several common single nucleotide polymorphisms (SNP) associated with an increased risk of colorectal cancer (CRC), although they have failed to identify any recessively acting alleles that contribute to disease risk. However, two recent studies have suggested that inbreeding and runs of homozygosity (ROH) increase the risk of developing cancer, perhaps by exposing recessive alleles as a result of autozygosity. To examine these results in a relatively large case-control series, we analyzed samples from a cohort in the United Kingdom comprising 921 colorectal tumor cases and 929 controls. Individuals were genotyped using a 550,000 tagging SNP panel. Additionally, we identified from these SNPs a set of approximately 30,000 SNPs in low pairwise linkage disequilibrium. To determine whether homozygosity was associated with CRC, we performed multiple tests to assess homozygosity at individual SNPs and ROHs in cases and controls. No association was found between CRC and (i) homozygosity at any individual SNP, (ii) overall homozygosity or level of inbreeding, (iii) total length or number of ROHs per individual, or (iv) a ROH at any particular genomic location. In conclusion, our results from a large case-control series do not replicate those of previous studies and suggest that homozygosity/autozygosity is not a major risk factor for CRC in an outbred population.

Disease severity and genetic pathways in attenuated familial adenomatous polyposis vary greatly but depend on the site of the germline mutation

BACKGROUND

Attenuated familial adenomatous polyposis (AFAP) is associated with germline mutations in the 5', 3', and exon 9 of the adenomatous polyposis coli (APC) gene. These mutations probably encode a limited amount of functional APC protein.

METHODS AND RESULTS

We found that colonic polyp number varied greatly among AFAP patients but members of the same family tended to have more similar disease severity. 5' Mutants generally had more polyps than other patients. We analysed somatic APC mutations/loss of heterozygosity (LOH) in 235 tumours from 35 patients (16 families) with a variety of AFAP associated germline mutations. In common with two previous studies of individual kindreds, we found biallelic changes ("third hits") in some polyps. We found that the "third hit" probably initiated tumorigenesis. Somatic mutation spectra were similar in 5' and 3' mutant patients, often resembling classical FAP. In exon 9 mutants, in contrast, "third hits" were more common. Most "third hits" left three 20 amino acid repeats (20AARs) on the germline mutant APC allele, with LOH (or proximal somatic mutation) of the wild-type allele; but some polyps had loss of the germline mutant with mutation leaving one 20AAR on the wild-type allele.

CONCLUSIONS

We propose that mutations, such as nt4661insA, that leave three 20AARs are preferentially selected in cis with some AFAP mutations because the residual protein function is near optimal for tumorigenesis. Not all AFAP polyps appear to need "three hits" however. AFAP is phenotypically and genetically heterogeneous. In addition to effects of different germline mutations, modifier genes may be acting on the AFAP phenotype, perhaps influencing the quantity of functional protein produced by the germline mutant allele.