Gene coexpression networks reveal a broad role for lncRNAs in inflammatory bowel disease.

The role of long noncoding RNAs (lncRNAs) in disease is incompletely understood, but their regulation of inflammation is increasingly appreciated. We addressed the extent of lncRNA involvement in inflammatory bowel disease (IBD) using biopsy-derived RNA-sequencing data from a large cohort of deeply phenotyped patients with IBD. Weighted gene correlation network analysis revealed gene modules of lncRNAs coexpressed with protein-coding genes enriched for biological pathways, correlated with epithelial and immune cell signatures, or correlated with distal colon expression. Correlation of modules with clinical features uncovered a module correlated with disease severity, with an enriched interferon response signature containing the hub lncRNA IRF1-AS1. Connecting genes to IBD-associated single nucleotide polymorphisms (SNPs) revealed an enrichment of SNP-adjacent lncRNAs in biologically relevant modules. Ulcerative colitis-specific SNPs were enriched in distal colon-related modules, suggesting that disease-specific mechanisms may result from altered lncRNA expression. The function of the IBD-associated SNP-adjacent lncRNA IRF1-AS1 was explored in human myeloid cells, and our results suggested IRF1-AS1 promoted optimal production of TNF-α, IL-6, and IL-23. A CRISPR/Cas9-mediated activation screen in THP-1 cells revealed several lncRNAs that modulated LPS-induced TNF-α responses. Overall, this study uncovered the expression patterns of lncRNAs in IBD that identify functional, disease-relevant lncRNAs.

Biopsy and blood-based molecular biomarker of inflammation in IBD.

OBJECTIVE: IBD therapies and treatments are evolving to deeper levels of remission. Molecular measures of disease may augment current endpoints including the potential for less invasive assessments. DESIGN: Transcriptome analysis on 712 endoscopically defined inflamed (Inf) and 1778 non-inflamed (Non-Inf) intestinal biopsies (n=498 Crohn's disease, n=421 UC and 243 controls) in the Mount Sinai Crohn's and Colitis Registry were used to identify genes differentially expressed between Inf and Non-Inf biopsies and to generate a molecular inflammation score (bMIS) via gene set variance analysis. A circulating MIS (cirMIS) score, reflecting intestinal molecular inflammation, was generated using blood transcriptome data. bMIS/cirMIS was validated as indicators of intestinal inflammation in four independent IBD cohorts. RESULTS: bMIS/cirMIS was strongly associated with clinical, endoscopic and histological disease activity indices. Patients with the same histologic score of inflammation had variable bMIS scores, indicating that bMIS describes a deeper range of inflammation. In available clinical trial data sets, both scores were responsive to IBD treatment. Despite similar baseline endoscopic and histologic activity, UC patients with lower baseline bMIS levels were more likely treatment responders compared with those with higher levels. Finally, among patients with UC in endoscopic and histologic remission, those with lower bMIS levels were less likely to have a disease flare over time. CONCLUSION: Transcriptionally based scores provide an alternative objective and deeper quantification of intestinal inflammation, which could augment current clinical assessments used for disease monitoring and have potential for predicting therapeutic response and patients at higher risk of disease flares.

Integrative Analysis of the Inflammatory Bowel Disease Serum Metabolome Improves Our Understanding of Genetic Etiology and Points to Novel Putative Therapeutic Targets.

BACKGROUND & AIMS: Polygenic and environmental factors are underlying causes of inflammatory bowel disease (IBD). We hypothesized that integration of the genetic loci controlling a metabolite's abundance, with known IBD genetic susceptibility loci, may help resolve metabolic drivers of IBD. METHODS: We measured the levels of 1300 metabolites in the serum of 484 patients with ulcerative colitis (UC) and 464 patients with Crohn's disease (CD) and 365 controls. Differential metabolite abundance was determined for disease status, subtype, clinical and endoscopic disease activity, as well as IBD phenotype including disease behavior, location, and extent. To inform on the genetic basis underlying metabolic diversity, we integrated metabolite and genomic data. Genetic colocalization and Mendelian randomization analyses were performed using known IBD risk loci to explore whether any metabolite was causally associated with IBD. RESULTS: We found 173 genetically controlled metabolites (metabolite quantitative trait loci, 9 novel) within 63 non-overlapping loci (7 novel). Furthermore, several metabolites significantly associated with IBD disease status and activity as defined using clinical and endoscopic indexes. This constitutes a resource for biomarker discovery and IBD biology insights. Using this resource, we show that a novel metabolite quantitative trait locus for serum butyrate levels containing ACADS was not supported as causal for IBD; replicate the association of serum omega-6 containing lipids with the fatty acid desaturase 1/2 locus and identify these metabolites as causal for CD through Mendelian randomization; and validate a novel association of serum plasmalogen and TMEM229B, which was predicted as causal for CD. CONCLUSIONS: An exploratory analysis combining genetics and unbiased serum metabolome surveys can reveal novel biomarkers of disease activity and potential mediators of pathology in IBD.

Molecular Characterization of Limited Ulcerative Colitis Reveals Novel Biology and Predictors of Disease Extension.

BACKGROUND AND AIMS: Disease extent varies in ulcerative colitis (UC) from proctitis to left-sided colitis to pancolitis and is a major prognostic factor. When the extent of UC is limited there is often a sharp demarcation between macroscopically involved and uninvolved areas and what defines this or subsequent extension is unknown. We characterized the demarcation site molecularly and determined genes associated with subsequent disease extension. METHODS: We performed RNA sequence analysis of biopsy specimens from UC patients with endoscopically and histologically confirmed limited disease, of which a subset later extended. Biopsy specimens were obtained from the endoscopically inflamed upper (proximal) limit of disease, immediately adjacent to the uninvolved colon, as well as at more proximal, endoscopically uninflamed colonic segments. RESULTS: Differentially expressed genes were identified in the endoscopically inflamed biopsy specimens taken at each patient's most proximal diseased site relative to healthy controls. Expression of these genes in the more proximal biopsy specimens transitioned back to control levels abruptly or gradually, the latter pattern supporting the concept that disease exists beyond the endoscopic disease demarcation site. The gradually transitioning genes were associated with inflammation, angiogenesis, glucuronidation, and homeodomain pathways. A subset of these genes in inflamed biopsy specimens was found to predict disease extension better than clinical features and were responsive to biologic therapies. Network analysis revealed critical roles for interferon signaling in UC inflammation and poly(ADP-ribose) polymerase 14 (PARP14) was a predicted key driver gene of extension. Higher PARP14 protein levels were found in inflamed biopsy specimens of patients with limited UC that subsequently extended. CONCLUSION: Molecular predictors of disease extension reveal novel strategies for disease prognostication and potential therapeutic targeting.

Deep Analysis of the Peripheral Immune System in IBD Reveals New Insight in Disease Subtyping and Response to Monotherapy or Combination Therapy.

BACKGROUND: Inflammatory bowel disease (IBD) is a complex disease with variable presentation, progression, and response to therapies. Current disease classification is based on subjective clinical phenotypes. The peripheral blood immunophenome can reflect local inflammation, and thus we measured 39 circulating immune cell types in a large cohort of IBD and control subjects and performed immunotype:phenotype associations. METHODS: We performed fluorescence-activated cell sorting or CyTOF analysis on blood from 728 Crohn's disease, 464 ulcerative colitis, and 334 non-IBD patients, with available demographics, endoscopic and clinical examinations and medication use. RESULTS: We observed few immune cell types commonly affected in IBD (lowered natural killer cells, B cells, and CD45RA- CD8 T cells). Generally, the immunophenome was distinct between ulcerative colitis and Crohn's disease. Within disease subtype, there were further distinctions, with specific immune cell types associating with disease duration, behavior, and location. Thiopurine monotherapy altered abundance of many cell types, often in the same direction as disease association, while anti-tumor necrosis factor (anti-TNF) monotherapy demonstrated an opposing pattern. Concomitant use of an anti-TNF and thiopurine was not synergistic, but rather was additive. For example, thiopurine monotherapy use alone or in combination with anti-TNF was associated with a dramatic reduction in major subclasses of B cells. CONCLUSIONS: We present a peripheral map of immune cell changes in IBD related to disease entity and therapies as a resource for hypothesis generation. We propose the changes in B cell subsets could affect antibody formation and potentially explain the mechanism behind the superiority of combination therapy through the impact of thiopurines on pharmacokinetics of anti-TNFs.

Polygenic risk score for alcohol drinking behavior improves prediction of inflammatory bowel disease risk.

Epidemiological studies have long recognized risky behaviors as potentially modifiable factors for the onset and flares of inflammatory bowel disease (IBD); yet, the underlying mechanisms are largely unknown. Recently, the genetic susceptibilities to cigarette smoking, alcohol and cannabis use [i.e. substance use (SU)] have been characterized by well-powered genome-wide association studies (GWASs). We aimed to assess the impact of genetic determinants of SU on IBD risk. Using Mount Sinai Crohn's and Colitis Registry (MSCCR) cohort of 1058 IBD cases and 188 healthy controls, we computed the polygenic risk score (PRS) for SU and correlated them with the observed IBD diagnoses, while adjusting for genetic ancestry, PRS for IBD and SU behavior at enrollment. The results were validated in a pediatric cohort with no SU exposure. PRS of alcohol consumption (DrnkWk), smoking cessation and age of smoking initiation, were associated with IBD risk in MSCCR even after adjustment for PRSIBD and actual smoking status. One interquartile range decrease in PRSDrnkWk was significantly associated to higher IBD risk (i.e. inverse association) (with odds ratio = 1.65 and 95% confidence interval: 1.32, 2.06). The association was replicated in a pediatric Crohn's disease cohort. Colocalization analysis identified a locus on chromosome 16 with polymorphisms in IL27, SULT1A2 and SH2B1, which reached genome-wide statistical significance in GWAS (P < 7.7e-9) for both alcohol consumption and IBD risk. This study demonstrated that the genetic predisposition to SU was associated with IBD risk, independent of PRSIBD and in the absence of SU behaviors. Our study may help further stratify individuals at risk of IBD.

Gut microbiota density influences host physiology and is shaped by host and microbial factors.

To identify factors that regulate gut microbiota density and the impact of varied microbiota density on health, we assayed this fundamental ecosystem property in fecal samples across mammals, human disease, and therapeutic interventions. Physiologic features of the host (carrying capacity) and the fitness of the gut microbiota shape microbiota density. Therapeutic manipulation of microbiota density in mice altered host metabolic and immune homeostasis. In humans, gut microbiota density was reduced in Crohn's disease, ulcerative colitis, and ileal pouch-anal anastomosis. The gut microbiota in recurrent Clostridium difficile infection had lower density and reduced fitness that were restored by fecal microbiota transplantation. Understanding the interplay between microbiota and disease in terms of microbiota density, host carrying capacity, and microbiota fitness provide new insights into microbiome structure and microbiome targeted therapeutics. Editorial note: This article has been through an editorial process in which the authors decide how to respond to the issues raised during peer review. The Reviewing Editor's assessment is that all the issues have been addressed (see decision letter).

Relationship between hysterectomy and admixture in African American women.

OBJECTIVE: Most studies suggest that hysterectomies are more common in African American women than in other ethnic groups. To assess this ethnic surgical disparity in a novel way, our main goal was to determine whether admixture (the proportion of sub-Saharan African or European origin in individuals) is associated with hysterectomy frequency in African American women in the Women's Health Initiative. STUDY DESIGN: In this retrospective study, we used ancestry informative single nucleotide polymorphisms to estimate admixture proportions in >10,000 African American women from the Women's Health Initiative. Logistic regression models were used to assess the association between admixture and self-reported history of hysterectomy with and without controls for relevant covariates. Multinomial logistic regression models were used to assess the association between admixture and self-reported age of hysterectomy. We also considered other potential risk factors (adiposity, hypertension, and education) for hysterectomy accounting for admixture. RESULTS: African admixture was a strong risk factor after the adjustment for multiple covariates (odds ratio, 1.85; P < .0001). The admixture risk for hysterectomy was highest for those procedures that were performed in the 35-39 age range (odds ratio, 3.08; P < .0001) and least evident in oldest ages (≥45 years old). Our analyses also suggest that adiposity, hypertension, and education were associated independently with hysterectomy in this population group. CONCLUSION: These results suggest that higher African admixture is associated with higher frequencies of hysterectomy and that genetic studies that specifically target African American women and diseases that are associated with hysterectomy may be especially useful in understanding the pathogenesis and underlying cause of this disparity in health outcome.

Risk for myasthenia gravis maps to a (151) Pro→Ala change in TNIP1 and to human leukocyte antigen-B*08.

OBJECTIVE: The objective of this study is to comprehensively define the genetic basis of early onset myasthenia gravis (EOMG). METHODS: We have carried out a 2-stage genome-wide association study on a total of 649 North European EOMG patients. Cases were matched 1:4 with controls of European ancestry. We performed imputation and conditional analyses across the major histocompatibility complex, as well as in the top regions of association outside the human leukocyte antigen (HLA) region. RESULTS: We observed the strongest association in the HLA class I region at rs7750641 (p = 1.2 × 10(-92) ; odds ratio [OR], 6.25). By imputation and conditional analyses, HLA-B*08 proves to be the major associated allele (p = 2.87 × 10(-113) ; OR, 6.41). In addition to the expected association with PTPN22 (rs2476601; OR, 1.71; p = 8.2 × 10(-10) ), an imputed coding variant (rs2233290) at position 151 (Pro→Ala) in the TNFAIP3-interacting protein 1, TNIP1, confers even stronger risk than PTPN22 (OR, 1.91; p = 3.2 × 10(-10) ). INTERPRETATION: The association at TNIP1 in EOMG implies disease mechanisms involving ubiquitin-dependent dysregulation of NF-κB signaling. The localization of the major HLA signal to the HLA-B*08 allele suggests that CD8(+) T cells may play a key role in disease initiation or pathogenesis.

Relationship between gallbladder surgery and ethnic admixture in African American and Hispanic American women.

OBJECTIVES: The objective of this study was to investigate whether differences in admixture in African American and Hispanic American adult women are associated with differences in gallbladder surgery. METHODS: Gallbladder surgery history on entry to the Women's Health Initiative's (WHI) study was used as a dichotomous outcome measure for this study. The proportion of European, sub-Saharan African, and Amerindian (AMI) admixture was estimated for 10,841 African American and 4,620 Hispanic American women in WHI using 92 ancestry informative markers. Logistic regression analyses assessed the relationship between admixture and gallbladder surgery in WHI women (enrollment at ages >50, mean age 61 years) with or without adjusting for multiple covariates, including measures of adiposity, parity, alcohol use, and education. RESULTS: There was a significant positive association between AMI admixture and the frequency of gallbladder surgery in Hispanic Americans. The odds ratio (OR) and 95% confidence intervals (CIs) for AMI admixture group was OR=2.97, CI=2.01-4.38, P<10(-4). Although there were strong positive associations with parity and adiposity, and negative associations with alcohol consumption and education, accounting for these covariates did not remove the admixture association (OR=2.46, CI=1.62-3.73). In contrast, the effect of African admixture was nearly indistinguishable from that of the European admixture, both of which were protective in the Hispanic American group, and African admixture had a marginal association with decreased gallbladder surgery in the African American group. Measures of adiposity were associated with increased risk for gallbladder surgery and remained significant after accounting for admixture and each of the other covariates. Education level and alcohol use were associated with decreased risk for gallbladder disease. However, after accounting for the other covariates these variably remained significant. CONCLUSIONS: AMI admixture is strongly associated with gallbladder surgery in women, even after adjustment for selected risk factors for cholelithiasis. Additional studies to ascertain the specific genetic risk factors underlying these associations are warranted.

Genome-wide association study identifies novel loci predisposing to cutaneous melanoma.

We performed a multistage genome-wide association study of melanoma. In a discovery cohort of 1804 melanoma cases and 1026 controls, we identified loci at chromosomes 15q13.1 (HERC2/OCA2 region) and 16q24.3 (MC1R) regions that reached genome-wide significance within this study and also found strong evidence for genetic effects on susceptibility to melanoma from markers on chromosome 9p21.3 in the p16/ARF region and on chromosome 1q21.3 (ARNT/LASS2/ANXA9 region). The most significant single-nucleotide polymorphisms (SNPs) in the 15q13.1 locus (rs1129038 and rs12913832) lie within a genomic region that has profound effects on eye and skin color; notably, 50% of variability in eye color is associated with variation in the SNP rs12913832. Because eye and skin colors vary across European populations, we further evaluated the associations of the significant SNPs after carefully adjusting for European substructure. We also evaluated the top 10 most significant SNPs by using data from three other genome-wide scans. Additional in silico data provided replication of the findings from the most significant region on chromosome 1q21.3 rs7412746 (P = 6 × 10(-10)). Together, these data identified several candidate genes for additional studies to identify causal variants predisposing to increased risk for developing melanoma.

A large-scale replication study identifies TNIP1, PRDM1, JAZF1, UHRF1BP1 and IL10 as risk loci for systemic lupus erythematosus.

Genome-wide association studies have recently identified at least 15 susceptibility loci for systemic lupus erythematosus (SLE). To confirm additional risk loci, we selected SNPs from 2,466 regions that showed nominal evidence of association to SLE (P < 0.05) in a genome-wide study and genotyped them in an independent sample of 1,963 cases and 4,329 controls. This replication effort identified five new SLE susceptibility loci (P < 5 x 10(-8)): TNIP1 (odds ratio (OR) = 1.27), PRDM1 (OR = 1.20), JAZF1 (OR = 1.20), UHRF1BP1 (OR = 1.17) and IL10 (OR = 1.19). We identified 21 additional candidate loci with P< or = 1 x 10(-5). A candidate screen of alleles previously associated with other autoimmune diseases suggested five loci (P < 1 x 10(-3)) that may contribute to SLE: IFIH1, CFB, CLEC16A, IL12B and SH2B3. These results expand the number of confirmed and candidate SLE susceptibility loci and implicate several key immunologic pathways in SLE pathogenesis.