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

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.

Associations of autozygosity with a broad range of human phenotypes

In many species, the offspring of related parents suffer reduced reproductive success, a phenomenon known as inbreeding depression. In humans, the importance of this effect has remained unclear, partly because reproduction between close relatives is both rare and frequently associated with confounding social factors. Here, using genomic inbreeding coefficients (FROH) for >1.4 million individuals, we show that FROH is significantly associated (p < 0.0005) with apparently deleterious changes in 32 out of 100 traits analysed. These changes are associated with runs of homozygosity (ROH), but not with common variant homozygosity, suggesting that genetic variants associated with inbreeding depression are predominantly rare. The effect on fertility is striking: FROH equivalent to the offspring of first cousins is associated with a 55% decrease [95% CI 44-66%] in the odds of having children. Finally, the effects of FROH are confirmed within full-sibling pairs, where the variation in FROH is independent of all environmental confounding.

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.

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 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<r_g<-0.55) and blood pressure (-0.35<r_g<-0.20). These findings provide clinically relevant biological insight into heritable variation in vagal heart rhythm regulation, with a key role for genetic variants (GNG11, RGS6) that influence G-protein heterotrimer action in GIRK-channel induced pacemaker membrane hyperpolarization.

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 ([Ca²⁺]_i) dynamics in human induced pluripotent stem cell‐derived cardiomyocytes (hiPSC‐CMs) generated from this individual (TECRL_Hom‐hiPSCs), his heterozygous but clinically asymptomatic father (TECRL_Het‐hiPSCs), and a healthy individual (CTRL‐hiPSCs) from the same Sudanese family, revealed smaller [Ca²⁺]_i transient amplitudes as well as elevated diastolic [Ca²⁺]_i in TECRL_Hom‐hiPSC‐CMs compared with CTRL‐hiPSC‐CMs. The [Ca²⁺]_i transient also rose markedly slower and contained lower sarcoplasmic reticulum (SR) calcium stores, evidenced by the decreased magnitude of caffeine‐induced [Ca²⁺]_i transients. In addition, the decay phase of the [Ca²⁺]_i transient was slower in TECRL_Hom‐hiPSC‐CMs due to decreased SERCA and NCX activities. Furthermore, TECRL_Hom‐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 TECRL_Hom‐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.

A Method to Exploit the Structure of Genetic Ancestry Space to Enhance Case-Control Studies

One goal of human genetics is to understand the genetic basis of disease, a challenge for diseases of complex inheritance because risk alleles are few relative to the vast set of benign variants. Risk variants are often sought by association studies in which allele frequencies in case subjects are contrasted with those from population-based samples used as control subjects. In an ideal world we would know population-level allele frequencies, releasing researchers to focus on case subjects. We argue this ideal is possible, at least theoretically, and we outline a path to achieving it in reality. If such a resource were to exist, it would yield ample savings and would facilitate the effective use of data repositories by removing administrative and technical barriers. We call this concept the Universal Control Repository Network (UNICORN), a means to perform association analyses without necessitating direct access to individual-level control data. Our approach to UNICORN uses existing genetic resources and various statistical tools to analyze these data, including hierarchical clustering with spectral analysis of ancestry; and empirical Bayesian analysis along with Gaussian spatial processes to estimate ancestry-specific allele frequencies. We demonstrate our approach using tens of thousands of control subjects from studies of Crohn disease, showing how it controls false positives, provides power similar to that achieved when all control data are directly accessible, and enhances power when control data are limiting or even imperfectly matched ancestrally. These results highlight how UNICORN can enable reliable, powerful, and convenient genetic association analyses without access to the individual-level data.

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

High-density mapping of the MHC identifies a shared role for HLA-DRB1*01:03 in inflammatory bowel diseases and heterozygous advantage in ulcerative colitis

Genome-wide association studies of the related chronic inflammatory bowel diseases (IBD) known as Crohn’s disease and ulcerative colitis have shown strong evidence of association to the major histocompatibility complex (MHC). This region encodes a large number of immunological candidates, including the antigen-presenting classical HLA molecules¹. Studies in IBD have indicated that multiple independent associations exist at HLA and non-HLA genes, but lacked the statistical power to define the architecture of association and causal alleles^2,3. To address this, we performed high-density SNP typing of the MHC in >32,000 patients with IBD, implicating multiple HLA alleles, with a primary role for HLA-DRB1*01:03 in both Crohn’s disease and ulcerative colitis. Significant differences were observed between these diseases, including a predominant role of class II HLA variants and heterozygous advantage observed in ulcerative colitis, suggesting an important role of the adaptive immune response to the colonic environment in the pathogenesis of IBD.

Deep resequencing of GWAS loci identifies rare variants in CARD9, IL23R and RNF186 that are associated with ulcerative colitis

Genome-wide association studies and follow-up meta-analyses in Crohn's disease (CD) and ulcerative colitis (UC) have recently identified 163 disease-associated loci that meet genome-wide significance for these two inflammatory bowel diseases (IBD). These discoveries have already had a tremendous impact on our understanding of the genetic architecture of these diseases and have directed functional studies that have revealed some of the biological functions that are important to IBD (e.g. autophagy). Nonetheless, these loci can only explain a small proportion of disease variance (~14% in CD and 7.5% in UC), suggesting that not only are additional loci to be found but that the known loci may contain high effect rare risk variants that have gone undetected by GWAS. To test this, we have used a targeted sequencing approach in 200 UC cases and 150 healthy controls (HC), all of French Canadian descent, to study 55 genes in regions associated with UC. We performed follow-up genotyping of 42 rare non-synonymous variants in independent case-control cohorts (totaling 14,435 UC cases and 20,204 HC). Our results confirmed significant association to rare non-synonymous coding variants in both IL23R and CARD9, previously identified from sequencing of CD loci, as well as identified a novel association in RNF186. With the exception of CARD9 (OR = 0.39), the rare non-synonymous variants identified were of moderate effect (OR = 1.49 for RNF186 and OR = 0.79 for IL23R). RNF186 encodes a protein with a RING domain having predicted E3 ubiquitin-protein ligase activity and two transmembrane domains. Importantly, the disease-coding variant is located in the ubiquitin ligase domain. Finally, our results suggest that rare variants in genes identified by genome-wide association in UC are unlikely to contribute significantly to the overall variance for the disease. Rather, these are expected to help focus functional studies of the corresponding disease loci.

Identification of heart rate-associated loci and their effects on cardiac conduction and rhythm disorders

Elevated resting heart rate is associated with greater risk of cardiovascular disease and mortality. In a 2-stage meta-analysis of genome-wide association studies in up to 181,171 individuals, we identified 14 new loci associated with heart rate and confirmed associations with all 7 previously established loci. Experimental downregulation of gene expression in Drosophila melanogaster and Danio rerio identified 20 genes at 11 loci that are relevant for heart rate regulation and highlight a role for genes involved in signal transmission, embryonic cardiac development and the pathophysiology of dilated cardiomyopathy, congenital heart failure and/or sudden cardiac death. In addition, genetic susceptibility to increased heart rate is associated with altered cardiac conduction and reduced risk of sick sinus syndrome, and both heart rate-increasing and heart rate-decreasing variants associate with risk of atrial fibrillation. Our findings provide fresh insights into the mechanisms regulating heart rate and identify new therapeutic targets.

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.

Genome-wide expression profiling implicates a MAST3-regulated gene set in colonic mucosal inflammation of ulcerative colitis patients

BACKGROUND

Crohn's disease (CD) and ulcerative colitis (UC) are inflammatory bowel diseases (IBDs) presumably caused by dysregulated immune responses to the gut microbiota. Genetic association studies have implicated dozens of chromosomal regions or loci in IBD susceptibility. The next challenge is to explain the individual role of each of these modest effect loci in the disease state. We have previously identified MAST3 as an IBD susceptibility gene through genetic fine-mapping of the 19p linkage region. Testing MAST3 in a reporter assay provided preliminary evidence that MAST3 modulates the activity of inflammation-related transcription factor nuclear factor kappa B.

METHODS

Here we characterized the function of MAST3 through an examination of the influence of the modulation of MAST3 expression on endogenous genome-wide expression patterns. More specifically, we looked at differential gene expression resulting from overexpression and knockdown of the MAST3 gene in epithelial and macrophage cell lines. From we highlight a group of genes whose expression is modulated by MAST3 and correlate their expression with NF-jB activity. Their expression was found to be enriched in inflamed mucosal tissue of UC patients, confirming the importance of these genes in IBD.

RESULTS

We highlight a group of genes whose expression is modulated by MAST3 and correlate their expression with NF-κB activity. Their expression was found to be enriched in inflamed mucosal tissue of UC patients, confirming the importance of these genes in IBD. These MAST3-regulated genes are central to mucosal immune responses. Among them are proinflammatory cytokines (e.g., CCL20, IL8), regulators of NF-κB (e.g., TNFAIP3, LY96, NFKBIA), genes involved in interferon-induced defense against pathogen invasion (e.g., IFIT1, ISG15), and genes involved in cell adhesion and/or migration (e.g., CD44, TMOD1).

CONCLUSIONS

Taken together, these results confirm MAST3 as a modulator of the inflammatory response through regulation of immune gene expression in the gut of IBD patients.

Expression and functional analysis of intestinal organic cation/L-carnitine transporter (OCTN) in Crohn's disease

BACKGROUND

The IBD5 locus is a genetic risk factor for IBD, particularly Crohn's Disease, coding for the organic cation/carnitine transporters (OCTN1 and 2). Two variants of OCTN are associated with susceptibility to Crohn's Disease. Modified transport of carnitine in vitro has been reported for a polymorphism of OCTN1. The aim was to investigate the function of intestinal OCTNs in IBD in relation to genetic polymorphisms.

METHODS

Intestinal tissue was obtained from endoscopic biopsies and surgical resections from IBD patients (n=33 and 14, resp.) and controls (n=22 and 14, resp.). OCTN protein levels were measured in intestinal biopsies and carnitine transport was quantified in intestinal resections.

RESULTS

OCTN1 protein levels were significantly higher in ileal versus colonic tissue (2.95% ± 0.4 vs 0.66% ± 0.2, resp.; p<0.0002). OCTN1 expression was higher in Crohn's disease patients with mutant homozygous or heterozygous genotypes (0.6% ± 0.1 vs 3% ± 0.8, resp., p<0.02). Carnitine transport was very rapid and Na+ dependent (10s). It was not different comparing Crohn's Disease and control groups (0.45 ± 0.12 vs 0.51 ± 0.12 nM carnitine/mg prot/min, resp.). Carnitine transport tended to be higher in subjects with mutant homozygous and heterozygous OCTN1 and OCTN2 genotypes (0.19 vs 0.59 and 0.25 vs 0.6, respectively).

CONCLUSIONS

The present data reveal that OCTN protein levels appear to be similar in intestinal tissue from Crohn's Disease patients and controls. Overall, ileal carnitine transport appears to as well equal in Crohn's Disease and control groups. However, there was a trend towards higher carnitine transport in subjects with OCTN1 and OCTN2 mutations.

Deep resequencing of GWAS loci identifies independent rare variants associated with inflammatory bowel disease

More than 1,000 susceptibility loci have been identified through genome-wide association studies (GWAS) of common variants; however, the specific genes and full allelic spectrum of causal variants underlying these findings have not yet been defined. Here we used pooled next-generation sequencing to study 56 genes from regions associated with Crohn's disease in 350 cases and 350 controls. Through follow-up genotyping of 70 rare and low-frequency protein-altering variants in nine independent case-control series (16,054 Crohn's disease cases, 12,153 ulcerative colitis cases and 17,575 healthy controls), we identified four additional independent risk factors in NOD2, two additional protective variants in IL23R, a highly significant association with a protective splice variant in CARD9 (P < 1 × 10(-16), odds ratio ≈ 0.29) and additional associations with coding variants in IL18RAP, CUL2, C1orf106, PTPN22 and MUC19. We extend the results of successful GWAS by identifying new, rare and probably functional variants that could aid functional experiments and predictive models.

Identification of a sudden cardiac death susceptibility locus at 2q24.2 through genome-wide association in European ancestry individuals

Sudden cardiac death (SCD) continues to be one of the leading causes of mortality worldwide, with an annual incidence estimated at 250,000–300,000 in the United States and with the vast majority occurring in the setting of coronary disease. We performed a genome-wide association meta-analysis in 1,283 SCD cases and >20,000 control individuals of European ancestry from 5 studies, with follow-up genotyping in up to 3,119 SCD cases and 11,146 controls from 11 European ancestry studies, and identify the BAZ2B locus as associated with SCD (P = 1.8×10⁻¹⁰). The risk allele, while ancestral, has a frequency of ∼1.4%, suggesting strong negative selection and increases risk for SCD by 1.92–fold per allele (95% CI 1.57–2.34). We also tested the role of 49 SNPs previously implicated in modulating electrocardiographic traits (QRS, QT, and RR intervals). Consistent with epidemiological studies showing increased risk of SCD with prolonged QRS/QT intervals, the interval-prolonging alleles are in aggregate associated with increased risk for SCD (P = 0.006).

Family studies have clearly demonstrated a role for genes in modifying risk for sudden cardiac death (SCD), however genetic studies have been limited by available samples. Here we have assembled over 4,400 SCD cases with >30,000 controls, all of European ancestry, and utilize a two-stage study design. In the first stage, we conducted an unbiased genome-wide scan in 1,283 SCD cases and >20,000 controls, and then performed follow-up genotyping in the remainder of the samples. We demonstrate strong association to a region of the genome not previously implicated in SCD, the BAZ2B locus, which contains 3 genes not previously known to play a role in cardiac biology. In addition, we used the genome-wide scan data to test a focused hypothesis that genetic variants that modulate ECG traits associated with SCD (QT, QRS, and RR intervals) also modify risk for SCD, and we demonstrate that QT- and QRS-prolonging alleles are, as a group, associated with increased risk of SCD. Taken together, these findings begin to elucidate the genetic contribution to SCD susceptibility and provide important targets for functional studies to investigate the etiology and pathogenesis of SCD.

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies and candidate gene studies in ulcerative colitis have identified 18 susceptibility loci. We conducted a meta-analysis of six ulcerative colitis genome-wide association study datasets, comprising 6,687 cases and 19,718 controls, and followed up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P < 5 × 10(-8)), increasing the number of ulcerative colitis-associated loci to 47. After annotating associated regions using GRAIL, expression quantitative trait loci data and correlations with non-synonymous SNPs, we identified many candidate genes that provide potentially important insights into disease pathogenesis, including IL1R2, IL8RA-IL8RB, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease risk loci is now 99, including a minimum of 28 shared association signals between Crohn's disease and ulcerative colitis.

Meta-analysis identifies 29 additional ulcerative colitis risk loci, increasing the number of confirmed associations to 47

Genome-wide association studies (GWAS) and candidate gene studies in ulcerative colitis (UC) have identified 18 susceptibility loci. We conducted a meta-analysis of 6 UC GWAS, comprising 6,687 cases and 19,718 controls, and followed-up the top association signals in 9,628 cases and 12,917 controls. We identified 29 additional risk loci (P<5×10^-8), increasing the number of UC associated loci to 47. After annotating associated regions using GRAIL, eQTL data and correlations with non-synonymous SNPs, we identified many candidate genes providing potentially important insights into disease pathogenesis, including IL1R2, IL8RA/B, IL7R, IL12B, DAP, PRDM1, JAK2, IRF5, GNA12 and LSP1. The total number of confirmed inflammatory bowel disease (IBD) risk loci is now 99, including a minimum of 28 shared association signals between Crohn’s disease (CD) and UC.

Genome-wide association identifies multiple ulcerative colitis susceptibility loci

Ulcerative colitis is a chronic, relapsing inflammatory condition of the gastrointestinal tract with a complex genetic and environmental etiology. In an effort to identify genetic variation underlying ulcerative colitis risk, we present two distinct genome-wide association studies of ulcerative colitis and their joint analysis with a previously published scan, comprising, in aggregate, 2,693 individuals with ulcerative colitis and 6,791 control subjects. Fifty-nine SNPs from 14 independent loci attained an association significance of P < 10(-5). Seven of these loci exceeded genome-wide significance (P < 5 x 10(-8)). After testing an independent cohort of 2,009 cases of ulcerative colitis and 1,580 controls, we identified 13 loci that were significantly associated with ulcerative colitis (P < 5 x 10(-8)), including the immunoglobulin receptor gene FCGR2A, 5p15, 2p16 and ORMDL3 (orosomucoid1-like 3). We confirmed association with 14 previously identified ulcerative colitis susceptibility loci, and an analysis of acknowledged Crohn's disease loci showed that roughly half of the known Crohn's disease associations are shared with ulcerative colitis. These data implicate approximately 30 loci in ulcerative colitis, thereby providing insight into disease pathogenesis.

Phenotypic and genotypic characteristics of inflammatory bowel disease in French Canadians: comparison with a large North American repository

OBJECTIVES

Phenotype characteristics of inflammatory bowel disease (IBD) may differ significantly among ethnic subpopulations. The aim of this study was to characterize the IBD phenotype in French Canadians, the most prominent founder population in North America.

METHODS

Using well-characterized phenotype data in the National Institute of Diabetes and Digestive and Kidney Diseases (NIDDK)-IBD Genetics Consortium repository on patients with IBD, we compared phenotypic characteristics of 202 French Canadians with those of 1,287 other Caucasian patients. These included diagnosis, anatomical location, disease behavior, extraintestinal manifestations, surgical history, and family history of IBD.

RESULTS

French-Canadian patients with Crohn's disease (CD) were less likely to have stricturing disease (11 vs. 21%, P=0.005; odds ratio (OR): 0.45, 95% confidence interval (95% CI): 0.24-0.85). Using a stringent definition of ethnicity (three out of four grandparents being French Canadians, as opposed to self-report, n=148), French Canadians had a tendency toward developing fistulizing CD (37 vs. 28%, P=0.07), and there was an increased prevalence of sacroiliitis among those with IBD (4 vs. 2%, P=0.045). Among French Canadians, the numbers of current smokers in CD (40 vs. 25%, P=0.006) and former smokers in ulcerative colitis (UC) (35 vs. 20%, P=0.03) were significantly higher. The prevalence of one of the three main variants of nucleotide-binding oligomerization domain containing 2 (NOD2) single-nucleotide polymorphisms (SNPs) among French-Canadian CD patients was 43.2%. The 3020insC SNP correlated with small bowel disease in French Canadians (25 [corrected] vs. 0%, P=0.006).

CONCLUSIONS

French Canadians show an IBD phenotype profile distinct from other Caucasian IBD populations, with an accentuated association between smoking status and IBD. This unique profile may have implications regarding the need for a different approach to the management of IBD in this population.

Genetic variants in the region harbouring IL2/IL21 associated with ulcerative colitis

OBJECTIVES

Genetic susceptibility is known to play a large part in the predisposition to the inflammatory bowel diseases (IBDs) known as Crohn's disease (CD) and ulcerative colitis (UC). The IL2/IL21 locus on 4q27 is known to be a common risk locus for inflammatory disease (shown in coeliac disease, type 1 diabetes, rheumatoid arthritis, systemic lupus erythematosus and psoriasis), while the roles that interleukin 2 (IL2) and IL21 play in the immune response also make them attractive candidates for IBD. The objective of this study was to test for association between the IL2/IL21 locus and the IBDs.

METHODS

The four single nucleotide polymorphisms (SNPs) in the IL2/IL21 locus most associated with coeliac disease were genotyped in 1590 subjects with IBD and 929 controls from The Netherlands, and then replicated in a North American cohort (2387 cases and 1266 controls) and an Italian cohort (805 cases and 421 controls), yielding a total of 4782 cases (3194 UC, 1588 CD) and 2616 controls. Allelic association testing and a pooled analysis using a Cochran-Mantel-Haenszel test were performed.

RESULTS

All four SNPs were strongly associated with UC in all three cohorts and reached genome-wide significance in the pooled analysis (rs13151961 p = 1.35 x 10(-10), rs13119723 p = 8.60 x 10(-8), rs6840978 p = 3.0 7x 10(-8), rs6822844 p = 2.77 x 10(-9)). A moderate association with CD was also found in the pooled analysis (p value range 0.0016-9.86 x 10(-5)).

CONCLUSIONS

A strong association for the IL2/IL21 locus with UC was found, which also confirms it as a general susceptibility locus for inflammatory disease.

MAST3: a novel IBD risk factor that modulates TLR4 signaling

Inflammatory bowel disease (IBD) is a chronic disorder caused by multiple factors in a genetically susceptible host. Significant advances in the study of genetic susceptibility have highlighted the importance of the innate immune system in this disease. We previously completed a genome-wide linkage study and found a significant locus (IBD6) on chromosome 19p. We were interested in identifying the causal variant in IBD6. We performed a two-stage association mapping study. In stage 1, 1530 single-nucleotide polymorphisms (SNPs) were selected from the HapMap database and genotyped in 761 patients with IBD. Among the SNPs that passed the threshold for replication, 26 were successfully genotyped in 754 additional patients (stage 2). One intronic variant, rs273506, located in the microtubule-associated serine/threonine-protein kinase gene-3 (MAST3), was found to be associated in both stages (pooled P=1.8 x 10(-4)). We identified four MAST3 coding variants, including a non-synonymous SNP rs8108738, correlated to rs273506 and associated with IBD. To test whether MAST3 was expressed in cells of interest, we performed expression assays, which showed abundant expression of MAST3 in antigen-presenting cells and in lymphocytes. The knockdown of MAST3 specifically decreased Toll-like receptor-4-dependent NF-kappaB activity. Our findings are additional proofs of the pivotal role played by modulators of NF-kappaB activity in IBD pathogenesis.

Defining the role of the MHC in autoimmunity: a review and pooled analysis

The major histocompatibility complex (MHC) is one of the most extensively studied regions in the human genome because of the association of variants at this locus with autoimmune, infectious, and inflammatory diseases. However, identification of causal variants within the MHC for the majority of these diseases has remained difficult due to the great variability and extensive linkage disequilibrium (LD) that exists among alleles throughout this locus, coupled with inadequate study design whereby only a limited subset of about 20 from a total of approximately 250 genes have been studied in small cohorts of predominantly European origin. We have performed a review and pooled analysis of the past 30 years of research on the role of the MHC in six genetically complex disease traits – multiple sclerosis (MS), type 1 diabetes (T1D), systemic lupus erythematosus (SLE), ulcerative colitis (UC), Crohn's disease (CD), and rheumatoid arthritis (RA) – in order to consolidate and evaluate the current literature regarding MHC genetics in these common autoimmune and inflammatory diseases. We corroborate established MHC disease associations and identify predisposing variants that previously have not been appreciated. Furthermore, we find a number of interesting commonalities and differences across diseases that implicate both general and disease-specific pathogenetic mechanisms in autoimmunity.

Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis

We present a genome-wide association study of ileal Crohn disease and two independent replication studies that identify several new regions of association to Crohn disease. Specifically, in addition to the previously established CARD15 and IL23R associations, we identified strong and significantly replicated associations (combined P < 10(-10)) with an intergenic region on 10q21.1 and a coding variant in ATG16L1, the latter of which was also recently reported by another group. We also report strong associations with independent replication to variation in the genomic regions encoding PHOX2B, NCF4 and a predicted gene on 16q24.1 (FAM92B). Finally, we demonstrate that ATG16L1 is expressed in intestinal epithelial cell lines and that functional knockdown of this gene abrogates autophagy of Salmonella typhimurium. Together, these findings suggest that autophagy and host cell responses to intracellular microbes are involved in the pathogenesis of Crohn disease.

Genome-wide association study identifies new susceptibility loci for Crohn disease and implicates autophagy in disease pathogenesis

We present a genome-wide association study of ileal Crohn disease and two independent replication studies that identify several new regions of association to Crohn disease. Specifically, in addition to the previously established CARD15 and IL23R associations, we identified strong and significantly replicated associations (combined P < 10(-10)) with an intergenic region on 10q21.1 and a coding variant in ATG16L1, the latter of which was also recently reported by another group. We also report strong associations with independent replication to variation in the genomic regions encoding PHOX2B, NCF4 and a predicted gene on 16q24.1 (FAM92B). Finally, we demonstrate that ATG16L1 is expressed in intestinal epithelial cell lines and that functional knockdown of this gene abrogates autophagy of Salmonella typhimurium. Together, these findings suggest that autophagy and host cell responses to intracellular microbes are involved in the pathogenesis of Crohn disease.

Molecular pathogenesis of inflammatory bowel disease: genotypes, phenotypes and personalized medicine

Crohn's disease (CD) and ulcerative colitis (UC), also known as inflammatory bowel diseases (IBD), are characterized by chronic inflammation of the gastrointestinal tract. IBD is among the few complex diseases for which several genomic regions and specific genes have been identified and confirmed in multiple replication studies. We will review the different loci implicated in disease risk in the context of three proposed mechanisms leading to chronic inflammation of the gut mucosa: 1) deregulation of the innate immune response to enteric microflora or pathogens; 2) increased permeability across the epithelial barrier; and 3) defective regulation of the adaptive immune system. As our knowledge of genetic variation, analytical approaches and technology improves, additional genetic risk factors are expected to be identified. With the identification of novel risk variants, additional pathophysiological mechanisms are likely to emerge. The resulting discoveries will further our molecular understanding of IBD, potentially leading to improved disease classification and rational drug design. Moreover, these approaches and tools can be applied in the context of variable drug response with the goal of providing more personalized clinical management of patients with IBD.

RARgamma acts as a tumor suppressor in mouse keratinocytes

All-trans retinoic acid (RA), the principle biologically active form of vitamin A, is essential for many developmental process as well as homeostasis in the adult. Many lines of evidence also suggest that RA, acting through the RA receptors (RARs), can also suppress growth of tumors of diverse origin. To assess directly the role of the RARs in a model of epidermal tumorigenesis, we investigated the incidence of tumor formation using keratinocytes lacking specific RAR types. Our data suggest that loss of RARgamma, but not RARalpha, predisposed keratinocytes to v-Ha-Ras-induced squamous cell carcinoma. We also found that ablation of RARgamma, but not RARalpha, abolished RA-induced cell cycle arrest and apoptosis in these keratinocytes. Reconstitution of receptor expression into RAR-null cells restored sensitivity to RA, and reversed the tumorigenic potential of receptor-deficient keratinocytes. These data strongly support a tumor suppressor effect for the RARs, in particular endogenous RARgamma, in murine keratinocytes.

Multiple transcription start sites and alternative splicing in the methylenetetrahydrofolate reductase gene result in two enzyme isoforms

Methylenetetrahydrofolate reductase (MTHFR) reduces 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate, the major carbon donor in the remethylation of homocysteine to methionine. Mild MTHFR deficiency, due to a common variant at nucleotide 677, has been reported to alter risk for several disorders including cardiovascular disease, neural tube defects, pregnancy complications, and certain cancers. Little is known about MTHFR regulation, since the complete cDNA and gene sequences have not been determined. In earlier work, we isolated and expressed a 2.2-kb human cDNA comprised of 11 coding exons, and we demonstrated that it encoded an active 70-kDa isoform. However, transcript sizes of approximately 7.5 kb and 9.5 kb and the presence of a second isoform of 77 kDa on Western blots suggested that cDNA sequences were incomplete. In this report, we characterized the complete cDNA and gene structure in human and mouse. Variable 5? and 3? UTR regions were identified, resulting in transcript heterogeneity. The 5? and 3? termini of the MTHFR cDNA were found to overlap with the 5? terminus of a chloride ion channel gene (CLCN-6) and the 3? terminus of an unidentified gene, respectively; this finding has resulted in finer mapping of MTHFR on Chromosome (Chr) 1p36.3. Ribonuclease protection assays identified clusters of transcriptional start sites, suggesting the existence of multiple promoters. MTHFR has several polyadenylation sites creating 3?UTR lengths of 0.2 kb-5.0 kb or 0.6 kb-4.0 kb in human and mouse, respectively. In both species, the previously reported exon 1 was redefined to approximately 3.0 kb in length and shown to be alternatively spliced. An important splice variant contains novel coding sequences; this cDNA was expressed and shown to encode the isozyme of 77 kDa. Our results, which suggest intricate regulation of MTHFR, will facilitate additional regulatory and functional studies of the different isoforms.

Regulation of gli activity by all-trans retinoic acid in mouse keratinocytes

Sonic hedgehog (Shh) signaling is essential for many normal developmental processes. The Shh signal is interpreted by the Gli transcription factors. Elevated Gli-1 expression has been associated with several neoplasms, including basal cell carcinoma. All-trans retinoic acid (RA) has strong effects on epidermal growth and differentiation and has been used for the treatment of various epithelial disorders. In this report, we show that RA can inhibit Gli activity in immortalized murine keratinocytes in a RA receptor-specific manner. This inhibition may occur, at least in part, through sequestration of the transcriptional coactivator cyclic AMP-responsive element-binding protein-binding protein and suggests a novel effect of retinoid excess on Shh signaling.

The thermolabile variant 677C-->T can further reduce activity when expressed in cis with severe mutations for human methylenetetrahydrofolate reductase

Methylenetetrahydrofolate reductase (MTHFR) catalyses the reduction of 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate, a carbon donor for homocysteine remethylation to methionine. Severe MTHFR deficiency is associated with hyperhomocysteinemia and homocystinuria. These patients show a wide variety of neurological and vascular symptoms, with variable age of onset. Residual enzyme activity is usually less than 20% of control values, and correlates reasonably well with age of onset of symptoms. A milder deficiency of MTHFR, with 30%-50% residual enzyme activity and increased enzyme thermolability, has been described as a risk factor for vascular disease and for neural tube defects. In earlier work, we isolated the human cDNA for MTHFR, and reported 14 mutations in severe MTHFR deficiency, as well as a common 677C-->T missense mutation (Ala-->Val) that encodes the thermolabile MTHFR. This variant has also been observed in some patients with severe MTHFR deficiency, in cis with their severe mutations. We report here the in vitro expression of seven severe MTHFR mutations in a bacterial expression system; six of these were expressed in cis with the Val allele to mimic the situation in the patients. We show that three of these constructs have significantly reduced enzyme activity (<10% of control); the presence of the thermolabile variant in these patients in cis is unlikely to affect enzyme function since activity is already low. One mutation causes a dramatic increase in activity when it is expressed in cis with the Ala allele, but is associated with extreme lability when in cis with the Val allele. Three mutations cause moderate decreases in enzyme activity, with a further decrease in activity when they are in cis with the Val allele. We hypothesize that deleterious mutations which alter stability may be compromised to a greater degree when the thermolabile variant is present on the same allele.

Characterization of retinoic acid receptor-deficient keratinocytes

Retinoids are essential for normal epidermal growth and differentiation and show potential for the prevention or treatment of various epithelial neoplasms. The retinoic acid receptors (RARalpha, -beta, and -gamma) are transducers of the retinoid signal. The epidermis expresses RARgamma and RARalpha, both of which are potential mediators of the effects of retinoids in the epidermis. To further investigate the role(s) of these receptors, we derived transformed keratinocyte lines from wild-type, RARalpha, RARgamma, and RARalphagamma null mice and investigated their response to retinoids, including growth inhibition, markers of growth and differentiation, and AP-1 activity. Our results indicate that RARgamma is the principle receptor contributing to all-trans-retinoic acid (RA)-mediated growth arrest in this system. This effect partially correlated with inhibition of AP-1 activity. In the absence of RARs, the synthetic retinoid N-(4-hydroxyphenyl)-retinamide inhibited growth; this was not observed with RA, 9-cis RA, or the synthetic retinoid (E)-4-[2-(5, 5, 8, 8 tetramethyl-5,6,7,8-tetrahydro-2-naphthalenyl)-1-propenyl] benzoic acid. Finally, both RARalpha and RARgamma differently affected the expression of some genes, suggesting both specific and overlapping roles for the RARs in keratinocytes.

Gene structure of human and mouse methylenetetrahydrofolate reductase (MTHFR)

Methylenetetrahydrofolate reductase (MTHFR) catalyzes the conversion of 5,10-methylenetetrahydrofolate to 5-methyltetrahydrofolate, a co-substrate for homocysteine remethylation to methionine. A human cDNA for MTHFR, 2.2 kb in length, has been expressed and shown to result in a catalytically active enzyme of approximately 70 kDa. Fifteen mutations have been identified in the MTHFR gene: 14 rare mutations associated with severe enzymatic deficiency and 1 common variant associated with a milder deficiency. The common polymorphism has been implicated in three multifactorial diseases: occlusive vascular disease, neural tube defects, and colon cancer. The human gene has been mapped to chromosomal region 1p36.3 while the mouse gene has been localized to distal Chromosome (Chr) 4. Here we report the isolation and characterization of the human and mouse genes for MTHFR. A human genomic clone (17 kb) was found to contain the entire cDNA sequence of 2.2 kb; there were 11 exons ranging in size from 102 bp to 432 bp. Intron sizes ranged from 250 bp to 1.5 kb with one exception of 4.2 kb. The mouse genomic clones (19 kb) start 7 kb 5' exon 1 and extend to the end of the coding sequence. The mouse amino acid sequence is approximately 90% identical to the corresponding human sequence. The exon sizes, locations of intronic boundaries, and intron sizes are also quite similar between the two species. The availability of human genomic clones has been useful in designing primers for exon amplification and mutation detection. The mouse genomic clones will be helpful in designing constructs for gene targeting and generation of mouse models for MTHFR deficiency.

Correlation of a common mutation in the methylenetetrahydrofolate reductase gene with plasma homocysteine in patients with premature coronary artery disease

Mild hyperhomocysteinemia, a risk factor for occlusive arterial disease, can be caused by disruptions of homocysteine metabolism. Methylenetetrahydrofolate reductase (MTHFR) catalyzes the synthesis of 5-methyltetrahydrofolate, the methyl donor for homocysteine remethylation to methionine. A common mutation in MTHFR, an alanine-to-valine substitution, may contribute to mild hyperhomocysteinemia in coronary artery disease (CAD). To test this hypothesis, we studied 152 patients with CAD by mutation analysis, MTHFR enzymatic assays, and measurements of plasma homocysteine and several vitamins. The MTHFR mutation was associated with reduced enzymatic activity and increased enzyme thermo-lability in these patients. The difference in the prevalence of the homozygous mutant genotype between the CAD patients (14%) and an unmatched group of healthy subjects (10%) was not significant. However, individuals with the homozygous mutant genotype had higher plasma homocysteine, particularly when plasma folate was below the median value. This genetic-environmental interaction is proposed to be a risk factor for CAD.

Human methionine synthase: cDNA cloning and identification of mutations in patients of the cblG complementation group of folate/cobalamin disorders

Methionine synthase catalyzes the remethylation of homocysteine to methionine in a methylcobalamin-dependent reaction. We used specific regions of homology within the methionine synthase sequences of several lower organisms to clone a human methionine synthase cDNA by a combination of RT-PCR and inverse PCR. The enzyme is 1265 amino acids in length and contains the seven residue structure-based sequence fingerprint identified for cobalamin-containing enzymes. The gene was localized to chromosome 1q43 by the FISH technique. We have identified one missense mutation and a 3 bp deletion in patients of the cblG complementation group of inherited homocysteine/folate disorders by SSCP and sequence analysis, as well as an amino acid substitution present in high frequency in the general population. We discuss the possibility that a mild deficiency of methionine synthase activity could be associated with mild hyperhomocysteinemia, a risk factor for cardiovascular disease and possibly neural tube defects.