Sex-Dimorphic Analyses Identify Novel and Sex-Specific Genetic Associations in Inflammatory Bowel Disease.

BACKGROUND: Sex is an integral variable often overlooked in complex disease genetics. Differences between sexes have been reported in natural history, disease complications, and age of onset in inflammatory bowel disease (IBD). While association studies have identified >230 IBD loci, there have been a limited number of studies investigating sex differences underlying these genetic associations. METHODS: We report the first investigation of sex-dimorphic associations via meta-analysis of a sex-stratified association study (34 579 IBD cases, 39 125 controls). In addition, we performed chromosome (chr) X-specific analyses, considering models of X inactivation (XCI) and XCI escape. Demographic and clinical characteristics were also compared between sexes. RESULTS: We identified significant differences between sexes for disease location and perianal complication in Crohn's disease and disease extent in ulcerative colitis. We observed genome-wide-significant sex-dimorphic associations (P < 5 × 10-8) at loci not previously reported in large-scale IBD genetic studies, including at chr9q22, CARMIL1, and UBASH3A. We identified variants in known IBD loci, including in chr2p15 and within the major histocompatibility complex on chr6, exhibiting sex-specific patterns of association (P < 5 × 10-7 in one sex only). We identified 3 chrX associations with IBD, including a novel Crohn's disease susceptibility locus at Xp22. CONCLUSIONS: These analyses identified novel IBD loci, in addition to characterizing sex-specific patterns of associations underlying sex-dimorphic associations. By elucidating the role of sex in IBD genetics, our study will help enhance our understanding of the differences between the sexes in IBD biology and underscores a need to move beyond conventional sex-combined analyses to appreciate the genetic architecture of IBD more comprehensively.

Sex-dimorphic meta-analyses of sex-stratified case-control (n = 73 704) regression identified 3 novel inflammatory bowel disease loci reaching genome-wide significance and highlighted chromosome 2 and major histocompatibility complex variants exhibiting sex-specific association. In addition, a novel chromosome X Crohn's disease susceptibility locus was identified.

Serum Lipidomic Screen Identifies Key Metabolites, Pathways, and Disease Classifiers in Crohn's Disease.

BACKGROUND: There is an unmet medical need for biomarkers that capture host and environmental contributions in inflammatory bowel diseases (IBDs). This study aimed at testing the potential of circulating lipids as disease classifiers given their major roles in inflammation. METHODS: We applied a previously validated comprehensive high-resolution liquid chromatography-mass spectrometry-based untargeted lipidomic workflow covering 25 lipid subclasses to serum samples from 100 Crohn's disease (CD) patients and 100 matched control subjects. Findings were replicated and expanded in another 200 CD patients and 200 control subjects. Key metabolites were tested for associations with disease behavior and location, and classification models were built and validated. Their association with disease activity was tested using an independent cohort of 42 CD patients. RESULTS: We identified >70 metabolites with strong association (P < 1 × 10-4, q < 5 × 10-4) to CD. Highly performing classification models (area under the curve > 0.84-0.97) could be built with as few as 5 to 9 different metabolites, representing 6 major correlated lipid clusters. These classifiers included a phosphatidylethanolamine ether (O-16:0/20:4), a sphingomyelin (d18:1/21:0) and a cholesterol ester (14:1), a very long-chain dicarboxylic acid [28:1(OH)] and sitosterol sulfate. These classifiers and correlated lipids indicate a dysregulated metabolism in host cells, notably in peroxisomes, as well as dysbiosis, oxidative stress, compromised inflammation resolution, or intestinal membrane integrity. A subset of these were associated with disease behavior or location. CONCLUSIONS: Untargeted lipidomic analyses uncovered perturbations in the circulating human CD lipidome, likely resulting from multiple pathogenic mechanisms. Models using as few as 5 biomarkers had strong disease classifier characteristics, supporting their potential use in diagnosis or prognosis.

This study reports a comprehensive untargeted lipidomic analysis of 600 serum samples from patients with Crohn's disease and matched control subjects, identified and replicated ~70 metabolites associated with Crohn's disease, and developed highly performing classification models (area under the curve > 0.84-0.97) with as few as 5 metabolites.

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.

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.

Serum Analyte Profiles Associated With Crohn's Disease and Disease Location.

BACKGROUND: Crohn's disease (CD) can affect any segment of the digestive tract but is most often localized in the ileal, ileocolonic, and colorectal regions of the intestines. It is believed that the chronic inflammation in CD is a result of an imbalance between the epithelial barrier, the immune system, and the intestinal microbiota. The aim of the study was to identify circulating markers associated with CD and/or disease location in CD patients. METHODS: We tested 49 cytokines, chemokines, and growth factors in serum samples from 300 patients with CD and 300 controls. After quality control, analyte levels were tested for association with CD and disease location. RESULTS: We identified 13 analytes that were higher in CD patients relative to healthy controls and that remained significant after conservative Bonferroni correction (P < 0.0015). In particular, CXCL9, CXCL1, and interleukin IL-6 had the greatest effect and were highly significant (P < 5 × 10-7). We also identified 9 analytes that were associated with disease location, with VEGF, IL-12p70, and IL-6 being elevated in patients with colorectal disease (P < 3 × 10-4). CONCLUSIONS: Multiple serum analytes are elevated in CD. These implicate the involvement of multiple cell types from the immune, epithelial, and endothelial systems, suggesting that circulating analytes reflect the inflammatory processes that are ongoing within the gut. Moreover, the identification of distinct profiles according to disease location supports the existence of a biological difference between ileal and colonic CD, consistent with previous genetic and clinical observations.

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.

Transethnic analysis of the human leukocyte antigen region for ulcerative colitis reveals not only shared but also ethnicity-specific disease associations.

Inflammatory bowel disease (IBD) is a chronic inflammatory disease of the gut. Genetic association studies have identified the highly variable human leukocyte antigen (HLA) region as the strongest susceptibility locus for IBD and specifically DRB1*01:03 as a determining factor for ulcerative colitis (UC). However, for most of the association signal such as delineation could not be made because of tight structures of linkage disequilibrium within the HLA. The aim of this study was therefore to further characterize the HLA signal using a transethnic approach. We performed a comprehensive fine mapping of single HLA alleles in UC in a cohort of 9272 individuals with African American, East Asian, Puerto Rican, Indian and Iranian descent and 40 691 previously analyzed Caucasians, additionally analyzing whole HLA haplotypes. We computationally characterized the binding of associated HLA alleles to human self-peptides and analyzed the physicochemical properties of the HLA proteins and predicted self-peptidomes. Highlighting alleles of the HLA-DRB1*15 group and their correlated HLA-DQ-DR haplotypes, we not only identified consistent associations (regarding effects directions/magnitudes) across different ethnicities but also identified population-specific signals (regarding differences in allele frequencies). We observed that DRB1*01:03 is mostly present in individuals of Western European descent and hardly present in non-Caucasian individuals. We found peptides predicted to bind to risk HLA alleles to be rich in positively charged amino acids. We conclude that the HLA plays an important role for UC susceptibility across different ethnicities. This research further implicates specific features of peptides that are predicted to bind risk and protective HLA proteins.

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.

Construction and benchmarking of a multi-ethnic reference panel for the imputation of HLA class I and II alleles.

Genotype imputation of the human leukocyte antigen (HLA) region is a cost-effective means to infer classical HLA alleles from inexpensive and dense SNP array data. In the research setting, imputation helps avoid costs for wet lab-based HLA typing and thus renders association analyses of the HLA in large cohorts feasible. Yet, most HLA imputation reference panels target Caucasian ethnicities and multi-ethnic panels are scarce. We compiled a high-quality multi-ethnic reference panel based on genotypes measured with Illumina's Immunochip genotyping array and HLA types established using a high-resolution next generation sequencing approach. Our reference panel includes more than 1,300 samples from Germany, Malta, China, India, Iran, Japan and Korea and samples of African American ancestry for all classical HLA class I and II alleles including HLA-DRB3/4/5. Applying extensive cross-validation, we benchmarked the imputation using the HLA imputation tool HIBAG, our multi-ethnic reference and an independent, previously published data set compiled of subpopulations of the 1000 Genomes project. We achieved average imputation accuracies higher than 0.924 for the commonly studied HLA-A, -B, -C, -DQB1 and -DRB1 genes across all ethnicities. We investigated allele-specific imputation challenges in regard to geographic origin of the samples using sensitivity and specificity measurements as well as allele frequencies and identified HLA alleles that are challenging to impute for each of the populations separately. In conclusion, our new multi-ethnic reference data set allows for high resolution HLA imputation of genotypes at all classical HLA class I and II genes including the HLA-DRB3/4/5 loci based on diverse ancestry populations.

Comprehensive and Reproducible Untargeted Lipidomic Workflow Using LC-QTOF Validated for Human Plasma Analysis.

The goal of this work was to develop a label-free, comprehensive, and reproducible high-resolution liquid chromatography-mass spectrometry (LC-MS)-based untargeted lipidomic workflow using a single instrument, which could be applied to biomarker discovery in both basic and clinical studies. For this, we have (i) optimized lipid extraction and elution to enhance coverage of polar and nonpolar lipids as well as resolution of their isomers, (ii) ensured MS signal reproducibility and linearity, and (iii) developed a bioinformatic pipeline to correct remaining biases. Workflow validation is reported for 48 replicates of a single human plasma sample: 1124 reproducible LC-MS signals were extracted (median signal intensity RSD = 10%), 50% of which are redundant due to adducts, dimers, in-source fragmentation, contaminations, or positive and negative ion duplicates. From the resulting 578 unique compounds, 428 lipids were identified by MS/MS, including acyl chain composition, of which 394 had RSD < 30% inside their linear intensity range, thereby enabling robust semiquantitation. MS signal intensity spanned 4 orders of magnitude, covering 16 lipid subclasses. Finally, the power of our workflow is illustrated by a proof-of-concept study in which 100 samples from healthy human subjects were analyzed and the data set was investigated using three different statistical testing strategies in order to compare their capacity in identifying the impact of sex and age on circulating lipids.

C1orf106 is a colitis risk gene that regulates stability of epithelial adherens junctions.

Polymorphisms in C1orf106 are associated with increased risk of inflammatory bowel disease (IBD). However, the function of C1orf106 and the consequences of disease-associated polymorphisms are unknown. Here we demonstrate that C1orf106 regulates adherens junction stability by regulating the degradation of cytohesin-1, a guanine nucleotide exchange factor that controls activation of ARF6. By limiting cytohesin-1-dependent ARF6 activation, C1orf106 stabilizes adherens junctions. Consistent with this model, C1orf106-/- mice exhibit defects in the intestinal epithelial cell barrier, a phenotype observed in IBD patients that confers increased susceptibility to intestinal pathogens. Furthermore, the IBD risk variant increases C1orf106 ubiquitination and turnover with consequent functional impairments. These findings delineate a mechanism by which a genetic polymorphism fine-tunes intestinal epithelial barrier integrity and elucidate a fundamental mechanism of cellular junctional control.

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.

Genetic Factors Interact With Tobacco Smoke to Modify Risk for Inflammatory Bowel Disease in Humans and Mice.

BACKGROUND & AIMS: The role of tobacco smoke in the etiology of inflammatory bowel disease (IBD) is unclear. We investigated interactions between genes and smoking (gene-smoking interactions) that affect risk for Crohn's disease (CD) and ulcerative colitis (UC) in a case-only study of patients and in mouse models of IBD. METHODS: We used 55 Immunochip-wide datasets that included 19,735 IBD cases (10,856 CD cases and 8879 UC cases) of known smoking status. We performed 3 meta-analyses each for CD, UC, and IBD (CD and UC combined), comparing data for never vs ever smokers, never vs current smokers, and never vs former smokers. We studied the effects of exposure to cigarette smoke in Il10-/- and Nod2-/- mice, as well as in Balb/c mice without disruption of these genes (wild-type mice). Mice were exposed to the smoke of 5 cigarettes per day, 5 days a week, for 8 weeks, in a ventilated smoking chamber, or ambient air (controls). Intestines were collected and analyzed histologically and by reverse transcription polymerase chain reaction. RESULTS: We identified 64 single nucleotide polymorphisms (SNPs) for which the association between the SNP and IBD were modified by smoking behavior (meta-analysis Wald test P < 5.0 × 10-5; heterogeneity Cochrane Q test P > .05). Twenty of these variants were located within the HLA region at 6p21. Analysis of classical HLA alleles (imputed from SNP genotypes) revealed an interaction with smoking. We replicated the interaction of a variant in NOD2 with current smoking in relation to the risk for CD (frameshift variant fs1007insC; rs5743293). We identified 2 variants in the same genomic region (rs2270368 and rs17221417) that interact with smoking in relation to CD risk. Approximately 45% of the SNPs that interact with smoking were in close vicinity (≤1 Mb) to SNPs previously associated with IBD; many were located near or within genes that regulate mucosal barrier function and immune tolerance. Smoking modified the disease risk of some variants in opposite directions for CD vs UC. Exposure of Interleukin 10 (il10)-deficient mice to cigarette smoke accelerated development of colitis and increased expression of interferon gamma in the small intestine compared to wild-type mice exposed to smoke. NOD2-deficient mice exposed to cigarette smoke developed ileitis, characterized by increased expression of interferon gamma, compared to wild-type mice exposed to smoke. CONCLUSIONS: In an analysis of 55 Immunochip-wide datasets, we identified 64 SNPs whose association with risk for IBD is modified by tobacco smoking. Gene-smoking interactions were confirmed in mice with disruption of Il10 and Nod2-variants of these genes have been associated with risk for IBD. Our findings from mice and humans revealed that the effects of smoking on risk for IBD depend on genetic variants.

Genetic Predictors of Benign Course of Ulcerative Colitis-A North American Inflammatory Bowel Disease Genetics Consortium Study.

BACKGROUND: A subset of patients with ulcerative colitis (UC) have a benign course and an overall favorable prognosis. Early identification of these low-risk patients may allow for a less aggressive therapeutic approach and possible reduction of therapy-associated risks. The aim of this project was to identify the genetic predictors of benign UC phenotype. METHODS: UC patients were selected from the National Institute of Diabetes and Digestive and Kidney Diseases Inflammatory Bowel Disease Genetics Consortium. Benign phenotype was defined as no need for immunomodulatory or biological therapy, hospitalizations, or colectomy. The association between benign UC phenotype and known loci linked to the risk of inflammatory bowel disease (IBD) was evaluated. The results for 156 index single-nucleotide polymorphisms (SNPs) from the known IBD loci were extracted for the main analysis. The association of the benign phenotype to a genetic burden score was also evaluated. RESULTS: None of the index SNPs from the IBD loci reached the predefined threshold of 1 × 10. In the exploratory analysis of the remaining Immunochip SNPs and imputed major histocompatibility complex data, 5 distinct suggestive association signals are identified (rs1697950, rs2523639, rs17836409, rs11742854, and rs75001121). CONCLUSIONS: No SNPs from IBD susceptibility loci were found to be associated (at our predefined threshold of 1 × 10) with a benign UC disease course. The rs11742570 variant on chromosome 5 was the one with the greatest association to benign disease although the association did not reach the predefined significant threshold. Given the modest power of our study, the findings suggested on the exploratory analysis merit extension to larger discovery cohorts.

Correction: The Influence of Age and Sex on Genetic Associations with Adult Body Size and Shape: A Large-Scale Genome-Wide Interaction Study.

[This corrects the article DOI: 10.1371/journal.pgen.1005378.].

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).