Identifying transcript-level differential expression in primary human immune cells.

BACKGROUND: Multipotential hematopoietic stem cells differentiate into a wide variety of immune cells with a diversity of functions, including the ability to respond to a variety of stimuli. Importantly, numerous studies have demonstrated the importance of gene transcription in defining cell identity and functions. While these studies have primarily been performed at the level of the gene, it is known that key immune genes such as CD44 and CD45 generate multiple different transcripts that are differentially expressed across different immune cells, and that encode proteins with different sequences and functions. Prior genomic surveys have shown that the mechanisms for generating diversity in expressed transcripts (alternate splicing, alternate transcription start sites, etc.) are very active in immune cells, but have been lacking in terms of identifying genes with multiple transcripts, that are differentially expressed, and likely to affect cell functions. METHODS: We first identified the set of genes that had at least two transcripts expressed in our RNA sequencing dataset generated from purified populations of neutrophils, monocytes and five lymphocyte populations (B, NK, γδ T, CD4 + T and CD8 + T) from twelve healthy donors. Next, we developed a heuristic approach to identify genes where two or more transcripts have distinct expression patterns across lymphoid and/or myeloid populations. We then focused our annotation and interpretation on differentially expressed transcripts that affect the coding sequence. This process was repeated to identify transcripts that were differentially expressed between monocytes and populations of macrophages and LPS-stimulated macrophages derived from these monocytes in vitro. RESULTS: We found that over 55 % of genes had two or more expressed transcripts, with an average ∼3 transcripts per gene, and that 70 % of these had at least two of the transcripts that encoded proteins with different sequences. As expected, we identified a complex pattern of differential expression for multiple transcripts encoding the CD45 transmembrane protein, but we also found similar evidence for ten other genes (CD300A, FYB1, GPI, LITAF, PSMA1, PTMA, RPL32, SEPTIN9, SH3BP2, SH3KBP1) when comparing the expression patterns of transcripts within myeloid and lymphoid cells. We also identified five genes with differentially expressed transcripts associated with the transition from monocytes to macrophages (FNBP1, KLF6, and SEPTIN9) or between macrophages and LPS-stimulated macrophages (CD44, OAZ2, and SEPTIN9). For the most part, we found that the different transcripts of these genes are expected to impact specific biological functions, for example the different transcripts of SEPTIN9 likely regulate the cytoskeleton in immune cells via their interactions with actins filaments and microtubules. CONCLUSIONS: This analytic approach successfully identified multi-transcript genes that are differentially expressed across immune cells and could be applied to other transcriptomic data. DATA AVAILABILITY STATEMENT: Researchers can request access to the individual-level data from the current study by contacting the Montreal Heart Institute ethics committee at the following institutional email address: cer.icm@icm-mhi.org.

Prostaglandins and calprotectin are genetically and functionally linked to the Inflammatory Bowel Diseases.

BACKGROUND: Genome wide association studies (GWAS) have identified and validated more than 200 genomic loci associated with the inflammatory bowel disease (IBD), although for most the causal gene remains unknown. Given the importance of myeloid cells in IBD pathogenesis, the current study aimed to uncover the role of genes within IBD genetic loci that are endogenously expressed in this cell lineage. METHODS: The open reading frames (ORF) of 42 genes from IBD-associated loci were expressed via lentiviral transfer in the THP-1 model of human monocytes and the impact of each of these on the cell's transcriptome was analyzed using a RNA sequencing-based approach. We used a combination of genetic and pharmacologic approaches to validate our findings in the THP-1 line with further validation in human induced pluripotent stem cell (hiPSC)-derived-monocytes. RESULTS: This functional genomics screen provided evidence that genes in four IBD GWAS loci (PTGIR, ZBTB40, SLC39A11 and NFKB1) are involved in controlling S100A8 and S100A9 gene expression, which encode the two subunits of calprotectin (CP). We demonstrated that increasing PTGIR expression and/or stimulating PTGIR signaling resulted in increased CP expression in THP-1. This was further validated in hiPSC-derived monocytes. Conversely, knocking-down PTGIR endogenous expression and/or inhibiting PTGIR signaling led to decreased CP expression. These analyses were extended to the known IBD gene PTGER4, whereby its specific agonist also led to increased CP expression. Furthermore, we demonstrated that the PTGIR and PTGER4 mediated control of CP expression was dependent on signaling via adenylate cyclase and STAT3. Finally, we demonstrated that LPS-mediated increases in CP expression could be potentiated by agonists of PTGIR and PTGER4, and diminished by their antagonists. CONCLUSION: Our results support a causal role for the PTGIR, PTGER4, ZBTB40, SLC39A11 and NFKB1 genes in IBD, with all five genes regulating the expression of CP in myeloid cells, as well as potential roles for the prostacyclin/prostaglandin biogenesis and signaling pathways in IBD susceptibility and pathogenesis.

Large-scale sequencing identifies multiple genes and rare variants associated with Crohn's disease susceptibility.

Genome-wide association studies (GWASs) have identified hundreds of loci associated with Crohn's disease (CD). However, as with all complex diseases, robust identification of the genes dysregulated by noncoding variants typically driving GWAS discoveries has been challenging. Here, to complement GWASs and better define actionable biological targets, we analyzed sequence data from more than 30,000 patients with CD and 80,000 population controls. We directly implicate ten genes in general onset CD for the first time to our knowledge via association to coding variation, four of which lie within established CD GWAS loci. In nine instances, a single coding variant is significantly associated, and in the tenth, ATG4C, we see additionally a significantly increased burden of very rare coding variants in CD cases. In addition to reiterating the central role of innate and adaptive immune cells as well as autophagy in CD pathogenesis, these newly associated genes highlight the emerging role of mesenchymal cells in the development and maintenance of intestinal inflammation.

Epigenetic reader SP140 loss of function drives Crohn's disease due to uncontrolled macrophage topoisomerases.

How mis-regulated chromatin directly impacts human immune disorders is poorly understood. Speckled Protein 140 (SP140) is an immune-restricted PHD and bromodomain-containing epigenetic "reader," and SP140 loss-of-function mutations associate with Crohn's disease (CD), multiple sclerosis (MS), and chronic lymphocytic leukemia (CLL). However, the relevance of these mutations and mechanisms underlying SP140-driven pathogenicity remains unexplored. Using a global proteomic strategy, we identified SP140 as a repressor of topoisomerases (TOPs) that maintains heterochromatin and macrophage fate. In humans and mice, SP140 loss resulted in unleashed TOP activity, de-repression of developmentally silenced genes, and ultimately defective microbe-inducible macrophage transcriptional programs and bacterial killing that drive intestinal pathology. Pharmacological inhibition of TOP1/2 rescued these defects. Furthermore, exacerbated colitis was restored with TOP1/2 inhibitors in Sp140-/- mice, but not wild-type mice, in vivo. Collectively, we identify SP140 as a TOP repressor and reveal repurposing of TOP inhibition to reverse immune diseases driven by SP140 loss.

Human enteric viruses autonomously shape inflammatory bowel disease phenotype through divergent innate immunomodulation.

Altered enteric microorganisms in concert with host genetics shape inflammatory bowel disease (IBD) phenotypes. However, insight is limited to bacteria and fungi. We found that eukaryotic viruses and bacteriophages (collectively, the virome), enriched from non-IBD, noninflamed human colon resections, actively elicited atypical anti-inflammatory innate immune programs. Conversely, ulcerative colitis or Crohn's disease colon resection viromes provoked inflammation, which was successfully dampened by non-IBD viromes. The IBD colon tissue virome was perturbed, including an increase in the enterovirus B species of eukaryotic picornaviruses, not previously detected in fecal virome studies. Mice humanized with non-IBD colon tissue viromes were protected from intestinal inflammation, whereas IBD virome mice exhibited exacerbated inflammation in a nucleic acid sensing-dependent fashion. Furthermore, there were detrimental consequences for IBD patient-derived intestinal epithelial cells bearing loss-of-function mutations within virus sensor MDA5 when exposed to viromes. Our results demonstrate that innate recognition of IBD or non-IBD human viromes autonomously influences intestinal homeostasis and disease phenotypes. Thus, perturbations in the intestinal virome, or an altered ability to sense the virome due to genetic variation, contribute to the induction of IBD. Harnessing the virome may offer therapeutic and biomarker potential.

IBD-associated G protein-coupled receptor 65 variant compromises signalling and impairs key functions involved in inflammation.

BACKGROUND AND AIMS: Inflammatory bowel diseases (IBD) result in chronic inflammation of the gastrointestinal tract. Genetic studies have shown that the GPR65 gene, as well as its missense coding variant, GPR65*Ile231Leu, is associated with IBD. We aimed to define the signalling and biological pathways downstream of GPR65 activation and evaluate the impact of GPR65*231Leu on these. METHODS: We used HEK 293 cells stably expressing GPR65 and deficient for either Gαs, Gαq/11 or Gα12/13, to define GPR65 signalling pathways, IBD patient biopsies and a panel of human tissues, primary immune cells and cell lines to determine biologic context, and genetic modulation of human THP-1-derived macrophages to examine the impact of GPR65 in bacterial phagocytosis and NLRP3 inflammasome activation. RESULTS: We confirmed that GPR65 signals via the Gαs pathway, leading to cAMP accumulation. GPR65 can also signal via the Gα12/13 pathway leading to formation of stress fibers, actin remodeling and RhoA activation; all impaired by the IBD-associated GPR65*231Leu allele. Gene expression profiling revealed greater expression of GPR65 in biopsies from inflamed compared to non-inflamed tissues from IBD patients or control individuals, potentially explained by infiltration of inflammatory immune cells. Decreased GPR65 expression in THP-1-derived macrophages leads to impaired bacterial phagocytosis, increased NLRP3 inflammasome activation and IL-1ß secretion in response to an inflammatory stimulus. CONCLUSIONS: We demonstrate that GPR65 exerts its effects through Gαs- and Gα12/13-mediated pathways, that the IBD-associated GPR65*231Leu allele has compromised interactions with Gα12/13 and that KD of GPR65 leads to impaired bacterial phagocytosis and increased inflammatory signalling via the NLRP3 inflammasome. This work identifies a target for development of small molecule therapies.

Functional screen of inflammatory bowel disease genes reveals key epithelial functions.

BACKGROUND: Genetic studies have been tremendously successful in identifying genomic regions associated with a wide variety of phenotypes, although the success of these studies in identifying causal genes, their variants, and their functional impacts has been more limited. METHODS: We identified 145 genes from IBD-associated genomic loci having endogenous expression within the intestinal epithelial cell compartment. We evaluated the impact of lentiviral transfer of the open reading frame (ORF) of these IBD genes into the HT-29 intestinal epithelial cell line via transcriptomic analyses. By comparing the genes in which expression was modulated by each ORF, as well as the functions enriched within these gene lists, we identified ORFs with shared impacts and their putative disease-relevant biological functions. RESULTS: Analysis of the transcriptomic data for cell lines expressing the ORFs for known causal genes such as HNF4a, IFIH1, and SMAD3 identified functions consistent with what is already known for these genes. These analyses also identified two major clusters of genes: Cluster 1 contained the known IBD causal genes IFIH1, SBNO2, NFKB1, and NOD2, as well as genes from other IBD loci (ZFP36L1, IRF1, GIGYF1, OTUD3, AIRE and PITX1), whereas Cluster 2 contained the known causal gene KSR1 and implicated DUSP16 from another IBD locus. Our analyses highlight how multiple IBD gene candidates can impact on epithelial structure and function, including the protection of the mucosa from intestinal microbiota, and demonstrate that DUSP16 acts a regulator of MAPK activity and contributes to mucosal defense, in part via its regulation of the polymeric immunoglobulin receptor, involved in the protection of the intestinal mucosa from enteric microbiota. CONCLUSIONS: This functional screen, based on expressing IBD genes within an appropriate cellular context, in this instance intestinal epithelial cells, resulted in changes to the cell's transcriptome that are relevant to their endogenous biological function(s). This not only helped in identifying likely causal genes within genetic loci but also provided insight into their biological functions. Furthermore, this work has highlighted the central role of intestinal epithelial cells in IBD pathophysiology, providing a scientific rationale for a drug development strategy that targets epithelial functions in addition to the current therapies targeting immune functions.

A transcriptome-based approach to identify functional modules within and across primary human immune cells.

Genome-wide transcriptomic analyses have provided valuable insight into fundamental biology and disease pathophysiology. Many studies have taken advantage of the correlation in the expression patterns of the transcriptome to infer a potential biologic function of uncharacterized genes, and multiple groups have examined the relationship between co-expression, co-regulation, and gene function on a broader scale. Given the unique characteristics of immune cells circulating in the blood, we were interested in determining whether it was possible to identify functional co-expression modules in human immune cells. Specifically, we sequenced the transcriptome of nine immune cell types from peripheral blood cells of healthy donors and, using a combination of global and targeted analyses of genes within co-expression modules, we were able to determine functions for these modules that were cell lineage-specific or shared among multiple cell lineages. In addition, our analyses identified transcription factors likely important for immune cell lineage commitment and/or maintenance.

Effect of Sex and Underlying Disease on the Genetic Association of QT Interval and Sudden Cardiac Death.

Background Sudden cardiac death (SCD) accounts for ≈300 000 deaths annually in the United States. Men have a higher risk of SCD and are more likely to have underlying coronary artery disease, while women are more likely to have arrhythmic events in the setting of inherited or acquired QT prolongation. Moreover, there is evidence of sex differences in the genetics of QT interval duration. Using sex- and coronary artery disease-stratified analyses, we assess differences in genetic association between longer QT interval and SCD risk. Methods and Results We examined 2282 SCD subjects and 3561 Finnish controls. The SCD subjects were stratified by underlying disease (ischemic versus nonischemic) and by sex. We used logistic regression to test for association between the top QT interval-associated single-nucleotide polymorphism, rs12143842 (in the NOS1AP locus), and SCD risk. We also performed Mendelian randomization to test for causal association of QT interval in the various subgroups. No statistically significant differences were observed between the sexes for associations with rs12143842, despite the odds ratio being higher in females across all subgroup analyses. Consistent with our hypothesis, female non-ischemics had the highest odds ratio point estimate for association between rs12143842 and SCD risk and male ischemics the lowest odds ratio point estimate (P=0.036 for difference). Similar trends were observed for the Mendelian randomization analysis. Conclusions While individual subgroup comparisons did not achieve traditional criteria for statistical significance, this study is consistent with the hypothesis that the causal association of longer QT interval on SCD risk is stronger in women and nonischemic individuals.

A comprehensive evaluation of the genetic architecture of sudden cardiac arrest.

Aims: Sudden cardiac arrest (SCA) accounts for 10% of adult mortality in Western populations. We aim to identify potential loci associated with SCA and to identify risk factors causally associated with SCA. Methods and results: We carried out a large genome-wide association study (GWAS) for SCA (n = 3939 cases, 25 989 non-cases) to examine common variation genome-wide and in candidate arrhythmia genes. We also exploited Mendelian randomization (MR) methods using cross-trait multi-variant genetic risk score associations (GRSA) to assess causal relationships of 18 risk factors with SCA. No variants were associated with SCA at genome-wide significance, nor were common variants in candidate arrhythmia genes associated with SCA at nominal significance. Using cross-trait GRSA, we established genetic correlation between SCA and (i) coronary artery disease (CAD) and traditional CAD risk factors (blood pressure, lipids, and diabetes), (ii) height and BMI, and (iii) electrical instability traits (QT and atrial fibrillation), suggesting aetiologic roles for these traits in SCA risk. Conclusions: Our findings show that a comprehensive approach to the genetic architecture of SCA can shed light on the determinants of a complex life-threatening condition with multiple influencing factors in the general population. The results of this genetic analysis, both positive and negative findings, have implications for evaluating the genetic architecture of patients with a family history of SCA, and for efforts to prevent SCA in high-risk populations and the general community.

Characterization of a Human Induced Pluripotent Stem Cell-Derived Cardiomyocyte Model for the Study of Variant Pathogenicity: Validation of a KCNJ2 Mutation.

BACKGROUND: Long-QT syndrome is a potentially fatal condition for which 30% of patients are without a genetically confirmed diagnosis. Rapid identification of causal mutations is thus a priority to avoid at-risk situations that can lead to fatal cardiac events. Massively parallel sequencing technologies are useful for the identification of sequence variants; however, electrophysiological testing of newly identified variants is crucial to demonstrate causality. Long-QT syndrome could, therefore, benefit from having a standardized platform for functional characterization of candidate variants in the physiological context of human cardiomyocytes. METHODS AND RESULTS: Using a variant in Kir2.1 (Gly52Val) revealed by whole-exome sequencing in a patient presenting with symptoms of long-QT syndrome as a proof of principle, we demonstrated that commercially available human induced pluripotent stem cell-derived cardiomyocytes are a powerful model for screening variants involved in genetic cardiac diseases. Immunohistochemistry experiments and whole-cell current recordings in human embryonic kidney cells expressing the wild-type or the mutant Kir2.1 demonstrated that Kir2.1-52V alters channel cellular trafficking and fails to form a functional channel. Using human induced pluripotent stem cell-derived cardiomyocytes, we not only confirmed these results but also further demonstrated that Kir2.1-52V is associated with a dramatic prolongation of action potential duration with evidence of arrhythmic activity, parameters which could not have been studied using human embryonic kidney cells. CONCLUSIONS: Our study confirms the pathogenicity of Kir2.1-52V in 1 patient with long-QT syndrome and also supports the use of isogenic human induced pluripotent stem cell-derived cardiomyocytes as a physiologically relevant model for the screening of variants of unknown function.

Erratum: Genetic loci associated with heart rate variability and their effects on cardiac disease risk.

This corrects the article DOI: 10.1038/ncomms15805.

Genetic loci associated with heart rate variability and their effects on cardiac disease risk.

Reduced cardiac vagal control reflected in low heart rate variability (HRV) is associated with greater risks for cardiac morbidity and mortality. In two-stage meta-analyses of genome-wide association studies for three HRV traits in up to 53,174 individuals of European ancestry, we detect 17 genome-wide significant SNPs in eight loci. HRV SNPs tag non-synonymous SNPs (in NDUFA11 and KIAA1755), expression quantitative trait loci (eQTLs) (influencing GNG11, RGS6 and NEO1), or are located in genes preferentially expressed in the sinoatrial node (GNG11, RGS6 and HCN4). Genetic risk scores account for 0.9 to 2.6% of the HRV variance. Significant genetic correlation is found for HRV with heart rate (-0.74

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.

TECRL, a new life-threatening inherited arrhythmia gene associated with overlapping clinical features of both LQTS and CPVT.

Genetic causes of many familial arrhythmia syndromes remain elusive. In this study, whole-exome sequencing (WES) was carried out on patients from three different families that presented with life-threatening arrhythmias and high risk of sudden cardiac death (SCD). Two French Canadian probands carried identical homozygous rare variant in TECRL gene (p.Arg196Gln), which encodes the trans-2,3-enoyl-CoA reductase-like protein. Both patients had cardiac arrest, stress-induced atrial and ventricular tachycardia, and QT prolongation on adrenergic stimulation. A third patient from a consanguineous Sudanese family diagnosed with catecholaminergic polymorphic ventricular tachycardia (CPVT) had a homozygous splice site mutation (c.331+1G>A) in TECRL Analysis of intracellular calcium ([Ca2+]i) dynamics in human induced pluripotent stem cell-derived cardiomyocytes (hiPSC-CMs) generated from this individual (TECRLHom-hiPSCs), his heterozygous but clinically asymptomatic father (TECRLHet-hiPSCs), and a healthy individual (CTRL-hiPSCs) from the same Sudanese family, revealed smaller [Ca2+]i transient amplitudes as well as elevated diastolic [Ca2+]i in TECRLHom-hiPSC-CMs compared with CTRL-hiPSC-CMs. The [Ca2+]i transient also rose markedly slower and contained lower sarcoplasmic reticulum (SR) calcium stores, evidenced by the decreased magnitude of caffeine-induced [Ca2+]i transients. In addition, the decay phase of the [Ca2+]i transient was slower in TECRLHom-hiPSC-CMs due to decreased SERCA and NCX activities. Furthermore, TECRLHom-hiPSC-CMs showed prolonged action potentials (APs) compared with CTRL-hiPSC-CMs. TECRL knockdown in control human embryonic stem cell-derived CMs (hESC-CMs) also resulted in significantly longer APs. Moreover, stimulation by noradrenaline (NA) significantly increased the propensity for triggered activity based on delayed afterdepolarizations (DADs) in TECRLHom-hiPSC-CMs and treatment with flecainide, a class Ic antiarrhythmic drug, significantly reduced the triggered activity in these cells. In summary, we report that mutations in TECRL are associated with inherited arrhythmias characterized by clinical features of both LQTS and CPVT Patient-specific hiPSC-CMs recapitulated salient features of the clinical phenotype and provide a platform for drug screening evidenced by initial identification of flecainide as a potential therapeutic. These findings have implications for diagnosis and treatment of inherited cardiac arrhythmias.

A protein-truncating R179X variant in RNF186 confers protection against ulcerative colitis.

Protein-truncating variants protective against human disease provide in vivo validation of therapeutic targets. Here we used targeted sequencing to conduct a search for protein-truncating variants conferring protection against inflammatory bowel disease exploiting knowledge of common variants associated with the same disease. Through replication genotyping and imputation we found that a predicted protein-truncating variant (rs36095412, p.R179X, genotyped in 11,148 ulcerative colitis patients and 295,446 controls, MAF=up to 0.78%) in RNF186, a single-exon ring finger E3 ligase with strong colonic expression, protects against ulcerative colitis (overall P=6.89 × 10(-7), odds ratio=0.30). We further demonstrate that the truncated protein exhibits reduced expression and altered subcellular localization, suggesting the protective mechanism may reside in the loss of an interaction or function via mislocalization and/or loss of an essential transmembrane domain.

Reduction in Na(+) current by angiotensin II is mediated by PKCα in mouse and human-induced pluripotent stem cell-derived cardiomyocytes.

BACKGROUND: Ventricular arrhythmias and sudden cardiac deaths are among the leading causes of mortality in patients with heart failure, and the underlying mechanisms remain incompletely understood. Chronic elevation of angiotensin II (ANGII) is known to be one of the main contributors to heart failure. OBJECTIVE: We tested whether ANGII can alter ventricular conduction and Na(+) current using transgenic mice with cardiomyocyte-restricted overexpression of ANGII type 1 receptor (AT1R). METHODS: We used surface electrocardiograms along with current- and voltage-clamp techniques to characterize the electrophysiological properties of AT1R mice while the underlying regulatory mechanisms were explored using reverse transcription/quantitative polymerase chain reaction, Western blots, and immunofluorescence techniques. RESULTS: Electrophysiological data indicated that chronic AT1R activation in ventricular myocytes caused a 60% reduction in Na(+) current density that slowed the maximal velocity of the action potential upstroke, leading to a prolongation of the QRS complex. These changes occur independently of cardiac hypertrophy, suggesting a direct role for ANGII/AT1R in slowing ventricular conduction. Western blots demonstrated a selective increase in sarcolemmal protein kinase Cα (PKCα) in AT1R mice, indicating PKCα activation. Furthermore, immunofluorescence analysis showed reorganization of PKCα expression to sarcolemma and colocalization with NaV1.5 in AT1R myocytes. The involvement of PKCα in regulating Na(+) current was subsequently demonstrated in human-induced pluripotent stem cell-derived cardiomyocytes where ANGII treatment reduced Na(+) current density. Concomitant treatment with αV5-3, a PKCα translocation inhibitor peptide, blocked the ANGII effect. CONCLUSION: Overall, this study suggests that in mouse and human cardiomyocytes, PKCα is an important mediator of the ANGII-induced reduction in Na(+) current and may contribute to ventricular arrhythmias.

Inherited determinants of Crohn's disease and ulcerative colitis phenotypes: a genetic association study.

BACKGROUND: Crohn's disease and ulcerative colitis are the two major forms of inflammatory bowel disease; treatment strategies have historically been determined by this binary categorisation. Genetic studies have identified 163 susceptibility loci for inflammatory bowel disease, mostly shared between Crohn's disease and ulcerative colitis. We undertook the largest genotype association study, to date, in widely used clinical subphenotypes of inflammatory bowel disease with the goal of further understanding the biological relations between diseases. METHODS: This study included patients from 49 centres in 16 countries in Europe, North America, and Australasia. We applied the Montreal classification system of inflammatory bowel disease subphenotypes to 34,819 patients (19,713 with Crohn's disease, 14,683 with ulcerative colitis) genotyped on the Immunochip array. We tested for genotype-phenotype associations across 156,154 genetic variants. We generated genetic risk scores by combining information from all known inflammatory bowel disease associations to summarise the total load of genetic risk for a particular phenotype. We used these risk scores to test the hypothesis that colonic Crohn's disease, ileal Crohn's disease, and ulcerative colitis are all genetically distinct from each other, and to attempt to identify patients with a mismatch between clinical diagnosis and genetic risk profile. FINDINGS: After quality control, the primary analysis included 29,838 patients (16,902 with Crohn's disease, 12,597 with ulcerative colitis). Three loci (NOD2, MHC, and MST1 3p21) were associated with subphenotypes of inflammatory bowel disease, mainly disease location (essentially fixed over time; median follow-up of 10·5 years). Little or no genetic association with disease behaviour (which changed dramatically over time) remained after conditioning on disease location and age at onset. The genetic risk score representing all known risk alleles for inflammatory bowel disease showed strong association with disease subphenotype (p=1·65 × 10(-78)), even after exclusion of NOD2, MHC, and 3p21 (p=9·23 × 10(-18)). Predictive models based on the genetic risk score strongly distinguished colonic from ileal Crohn's disease. Our genetic risk score could also identify a small number of patients with discrepant genetic risk profiles who were significantly more likely to have a revised diagnosis after follow-up (p=6·8 × 10(-4)). INTERPRETATION: Our data support a continuum of disorders within inflammatory bowel disease, much better explained by three groups (ileal Crohn's disease, colonic Crohn's disease, and ulcerative colitis) than by Crohn's disease and ulcerative colitis as currently defined. Disease location is an intrinsic aspect of a patient's disease, in part genetically determined, and the major driver to changes in disease behaviour over time. FUNDING: International Inflammatory Bowel Disease Genetics Consortium members funding sources (see Acknowledgments for full list).