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

Genetic architecture of the inflammatory bowel diseases across East Asian and European ancestries

Inflammatory bowel diseases (IBDs) are chronic disorders of the gastrointestinal tract with the following two subtypes: Crohn's disease (CD) and ulcerative colitis (UC). To date, most IBD genetic associations were derived from individuals of European (EUR) ancestries. Here we report the largest IBD study of individuals of East Asian (EAS) ancestries, including 14,393 cases and 15,456 controls. We found 80 IBD loci in EAS alone and 320 when meta-analyzed with ~370,000 EUR individuals (~30,000 cases), among which 81 are new. EAS-enriched coding variants implicate many new IBD genes, including ADAP1 and GIT2. Although IBD genetic effects are generally consistent across ancestries, genetics underlying CD appears more ancestry dependent than UC, driven by allele frequency (NOD2) and effect (TNFSF15). We extended the IBD polygenic risk score (PRS) by incorporating both ancestries, greatly improving its accuracy and highlighting the importance of diversity for the equitable deployment of PRS.

Genetic coding variant in complement factor B (CFB) is associated with increased risk for perianal Crohn's disease and leads to impaired CFB cleavage and phagocytosis

OBJECTIVE

Perianal Crohn's disease (pCD) occurs in up to 40% of patients with CD and is associated with poor quality of life, limited treatment responses and poorly understood aetiology. We performed a genetic association study comparing CD subjects with and without perianal disease and subsequently performed functional follow-up studies for a pCD associated SNP in Complement Factor B (CFB).

DESIGN

Immunochip-based meta-analysis on 4056 pCD and 11 088 patients with CD from three independent cohorts was performed. Serological and clinical variables were analysed by regression analyses. Risk allele of rs4151651 was introduced into human CFB plasmid by site-directed mutagenesis. Binding of recombinant G252 or S252 CFB to C3b and its cleavage was determined in cell-free assays. Macrophage phagocytosis in presence of recombinant CFB or serum from CFB risk, or protective CD or healthy subjects was assessed by flow cytometry.

RESULTS

Perianal complications were associated with colonic involvement, OmpC and ASCA serology, and serology quartile sum score. We identified a genetic association for pCD (rs4151651), a non-synonymous SNP (G252S) in CFB, in all three cohorts. Recombinant S252 CFB had reduced binding to C3b, its cleavage was impaired, and complement-driven phagocytosis and cytokine secretion were reduced compared with G252 CFB. Serine 252 generates a de novo glycosylation site in CFB. Serum from homozygous risk patients displayed significantly decreased macrophage phagocytosis compared with non-risk serum.

CONCLUSION

pCD-associated rs4151651 in CFB is a loss-of-function mutation that impairs its cleavage, activation of alternative complement pathway, and pathogen phagocytosis thus implicating the alternative complement pathway and CFB in pCD aetiology.

An Extremes of Phenotype Approach Confirms Significant Genetic Heterogeneity in Patients with Ulcerative Colitis

BACKGROUND AND AIMS

Ulcerative colitis [UC] is a major form of inflammatory bowel disease globally. Phenotypic heterogeneity is defined by several variables including age of onset and disease extent. The genetics of disease severity remains poorly understood. To further investigate this, we performed a genome wide association [GWA] study using an extremes of phenotype strategy.

METHODS

We conducted GWA analyses in 311 patients with medically refractory UC [MRUC], 287 with non-medically refractory UC [non-MRUC] and 583 controls. Odds ratios [ORs] were calculated for known risk variants comparing MRUC and non-MRUC, and controls.

RESULTS

MRUC-control analysis had the greatest yield of genome-wide significant single nucleotide polymorphisms [SNPs] [2018], including lead SNP = rs111838972 [OR = 1.82, p = 6.28 × 10-9] near MMEL1 and a locus in the human leukocyte antigen [HLA] region [lead SNP = rs144717024, OR = 12.23, p = 1.7 × 10-19]. ORs for the lead SNPs were significantly higher in MRUC compared to non-MRUC [p < 9.0 × 10-6]. No SNPs reached significance in the non-MRUC-control analysis (top SNP, rs7680780 [OR 2.70, p = 5.56 × 10-8). We replicate findings for rs4151651 in the Complement Factor B [CFB] gene and demonstrate significant changes in CFB gene expression in active UC. Detailed HLA analyses support the strong associations with MHC II genes, particularly HLA-DQA1, HLA-DQB1 and HLA-DRB1 in MRUC.

CONCLUSIONS

Our MRUC subgroup replicates multiple known UC risk variants in contrast to non-MRUC and demonstrates significant differences in effect sizes compared to those published. Non-MRUC cases demonstrate lower ORs similar to those published. Additional risk and prognostic loci may be identified by targeted recruitment of individuals with severe disease.

NOMe-HiC: joint profiling of genetic variant, DNA methylation, chromatin accessibility, and 3D genome in the same DNA molecule

Cis-regulatory elements are coordinated to regulate the expression of their targeted genes. However, the joint measurement of cis-regulatory elements' activities and their interactions in spatial proximity is limited by the current sequencing approaches. We describe a method, NOMe-HiC, which simultaneously captures single-nucleotide polymorphisms, DNA methylation, chromatin accessibility (GpC methyltransferase footprints), and chromosome conformation changes from the same DNA molecule, together with the transcriptome, in a single assay. NOMe-HiC shows high concordance with state-of-the-art mono-omic assays across different molecular measurements and reveals coordinated chromatin accessibility at distal genomic segments in spatial proximity and novel types of long-range allele-specific chromatin accessibility.

Development of antidrug antibodies against adalimumab maps to variation within the HLA-DR peptide-binding groove

Targeted biologic therapies can elicit an undesirable host immune response characterized by the development of antidrug antibodies (ADA), an important cause of treatment failure. The most widely used biologic across immune-mediated diseases is adalimumab, a tumor necrosis factor inhibitor. This study aimed to identify genetic variants that contribute to the development of ADA against adalimumab, thereby influencing treatment failure. In patients with psoriasis on their first course of adalimumab, in whom serum ADA had been evaluated 6-36 months after starting treatment, we observed a genome-wide association with ADA against adalimumab within the major histocompatibility complex (MHC). The association signal mapped to the presence of tryptophan at position 9 and lysine at position 71 of the HLA-DR peptide-binding groove, with both residues conferring protection against ADA. Underscoring their clinical relevance, these residues were also protective against treatment failure. Our findings highlight antigenic peptide presentation via MHC class II as a critical mechanism in the development of ADA against biologic therapies and downstream treatment response.

Proteomics: Its Promise and Pitfalls in Shaping Precision Medicine in Solid Organ Transplantation

Solid organ transplantation is an established treatment of choice for end-stage organ failure. However, all transplant patients are at risk of developing complications, including allograft rejection and death. Histological analysis of graft biopsy is still the gold standard for evaluation of allograft injury, but it is an invasive procedure and prone to sampling errors. The past decade has seen an increased number of efforts to develop minimally invasive procedures for monitoring allograft injury. Despite the recent progress, limitations such as the complexity of proteomics-based technology, the lack of standardization, and the heterogeneity of populations that have been included in different studies have hindered proteomic tools from reaching clinical transplantation. This review focuses on the role of proteomics-based platforms in biomarker discovery and validation in solid organ transplantation. We also emphasize the value of biomarkers that provide potential mechanistic insights into the pathophysiology of allograft injury, dysfunction, or rejection. Additionally, we forecast that the growth of publicly available data sets, combined with computational methods that effectively integrate them, will facilitate a generation of more informed hypotheses for potential subsequent evaluation in preclinical and clinical studies. Finally, we illustrate the value of combining data sets through the integration of 2 independent data sets that pinpointed hub proteins in antibody-mediated rejection.

Type 2 diabetes candidate genes, including PAX5, cause impaired insulin secretion in human pancreatic islets

Type 2 diabetes (T2D) is caused by insufficient insulin secretion from pancreatic β cells. To identify candidate genes contributing to T2D pathophysiology, we studied human pancreatic islets from approximately 300 individuals. We found 395 differentially expressed genes (DEGs) in islets from individuals with T2D, including, to our knowledge, novel (OPRD1, PAX5, TET1) and previously identified (CHL1, GLRA1, IAPP) candidates. A third of the identified expression changes in islets may predispose to diabetes, as expression of these genes associated with HbA1c in individuals not previously diagnosed with T2D. Most DEGs were expressed in human β cells, based on single-cell RNA-Seq data. Additionally, DEGs displayed alterations in open chromatin and associated with T2D SNPs. Mouse KO strains demonstrated that the identified T2D-associated candidate genes regulate glucose homeostasis and body composition in vivo. Functional validation showed that mimicking T2D-associated changes for OPRD1, PAX5, and SLC2A2 impaired insulin secretion. Impairments in Pax5-overexpressing β cells were due to severe mitochondrial dysfunction. Finally, we discovered PAX5 as a potential transcriptional regulator of many T2D-associated DEGs in human islets. Overall, we have identified molecular alterations in human pancreatic islets that contribute to β cell dysfunction in T2D pathophysiology.

Therapeutic interventions target the NLRP3 inflammasome in ulcerative colitis: Comprehensive study

One of the significant health issues in the world is the prevalence of ulcerative colitis (UC). UC is a chronic disorder that mainly affects the colon, beginning with the rectum, and can progress from asymptomatic mild inflammation to extensive inflammation of the entire colon. Understanding the underlying molecular mechanisms of UC pathogenesis emphasizes the need for innovative therapeutic approaches based on identifying molecular targets. Interestingly, in response to cellular injury, the NLR family pyrin domain containing 3 (NLRP3) inflammasome is a crucial part of the inflammation and immunological reaction by promoting caspase-1 activation and the release of interleukin-1β. This review discusses the mechanisms of NLRP3 inflammasome activation by various signals and its regulation and impact on UC.

Transcriptomic analyses of treatment-naïve pediatric ulcerative colitis patients and exploration of underlying disease pathogenesis

BACKGROUND

Ulcerative colitis (UC) is a form of chronic inflammatory bowel disease of nonspecific origin. This study used an RNA-Sequencing (RNA-Seq) approach to evaluate the transcriptomic landscape of a well-stratified treatment-naïve pediatric UC patient population by comparing them with healthy control children. The data were analyzed to evaluate the mechanisms driving UC-related intestinal inflammation and fibrosis.

METHODS

Intestinal mucosal samples from five pediatric UC patients and five healthy controls were analyzed by RNA-Seq, and results were verified by qPCR. A CRISPR/Cas9 approach was used to knock out the expression of HLA-DRB5, and molecular biology techniques were used for additional mechanistic studies.

RESULTS

In these analyses, 2290 genes were found to be differentially expressed between the UC and control samples, of which 1258 and 1032 were upregulated and downregulated, respectively. Gene Ontology analysis showed that these genes were enriched in extracellular matrix (ECM)-related processes and that 7 of 8 differentially expressed genes of interest (PIK3CD, IL1β, IL1α, TIMP1, MMP1, MMP12, COL6A3, and HLADRB5) were upregulated and involved in ECM-receptor interaction and inflammatory bowel disease-related pathways. Increased HLA-DRB5 expression driven by intestinal bacteria was found to promote IL-1α secretion, leading to intestinal inflammation and fibrosis, suggesting a possible target for the treatment of UC.

CONCLUSION

These data suggest that intestinal inflammation is present in pediatric UC patients for extended periods before the onset of symptoms, and intestinal fibrosis begins even during the early stages of UC. Intestinal bacteria were also found to trigger intestinal inflammation and fibrosis, with HLA-DRB5 playing a central role in this process.

Long-term Temporal Stability of Peripheral Blood DNA Methylation Profiles in Patients With Inflammatory Bowel Disease

BACKGROUND & AIMS

There is great current interest in the potential application of DNA methylation alterations in peripheral blood leukocytes (PBLs) as biomarkers of susceptibility, progression, and treatment response in inflammatory bowel disease (IBD). However, the intra-individual stability of PBL methylation in IBD has not been characterized. Here, we studied the long-term stability of all probes located on the Illumina HumanMethylation EPIC BeadChip array.

METHODS

We followed a cohort of 46 adult patients with IBD (36 Crohn's disease [CD], 10 ulcerative colitis [UC]; median age, 44 years; interquartile range [IQR] 27-56 years; 50% female) that received standard care follow-up at the Amsterdam University Medical Centers. Paired PBL samples were collected at 2 time points with a median of 7 years (range, 2-9 years) in between. Differential methylation and intra-class correlation (ICC) analyses were used to identify time-associated differences and temporally stable CpGs, respectively.

RESULTS

Around 60% of all EPIC array loci presented poor intra-individual stability (ICC <0.50); 78.114 (≈9%) showed good (ICC, 0.75-0.89), and 41.274 (≈5%) showed excellent (ICC ≥0.90) stability, between both measured time points. Focusing on previously identified consistently differentially methylated positions indicated that 22 CD-, 11 UC-, and 24 IBD-associated loci demonstrated high stability (ICC ≥0.75) over time; of these, we observed a marked stability of CpG loci associated to the HLA genes.

CONCLUSIONS

Our data provide insight into the long-term stability of the PBL DNA methylome within an IBD context, facilitating the selection of biologically relevant and robust IBD-associated epigenetic biomarkers with increased potential for independent validation. These data also have potential implications in understanding disease pathogenesis.

The missing link between genetic association and regulatory function

The genetic basis of most traits is highly polygenic and dominated by non-coding alleles. It is widely assumed that such alleles exert small regulatory effects on the expression of cis-linked genes. However, despite the availability of gene expression and epigenomic datasets, few variant-to-gene links have emerged. It is unclear whether these sparse results are due to limitations in available data and methods, or to deficiencies in the underlying assumed model. To better distinguish between these possibilities, we identified 220 gene-trait pairs in which protein-coding variants influence a complex trait or its Mendelian cognate. Despite the presence of expression quantitative trait loci near most GWAS associations, by applying a gene-based approach we found limited evidence that the baseline expression of trait-related genes explains GWAS associations, whether using colocalization methods (8% of genes implicated), transcription-wide association (2% of genes implicated), or a combination of regulatory annotations and distance (4% of genes implicated). These results contradict the hypothesis that most complex trait-associated variants coincide with homeostatic expression QTLs, suggesting that better models are needed. The field must confront this deficit and pursue this 'missing regulation.'