Mapping autophagosome contents identifies interleukin-7 receptor-α as a key cargo modulating CD4+ T cell proliferation.

CD4+ T cells are pivotal cells playing roles in the orchestration of humoral and cytotoxic immune responses. It is known that CD4+ T cell proliferation relies on autophagy, but identification of the autophagosomal cargo involved is missing. Here we create a transgenic mouse model, to enable direct mapping of the proteinaceous content of autophagosomes in primary cells by LC3 proximity labelling. Interleukin-7 receptor-α, a cytokine receptor mostly found in naïve and memory T cells, is reproducibly detected in autophagosomes of activated CD4+ T cells. Consistently, CD4+ T cells lacking autophagy show increased interleukin-7 receptor-α surface expression, while no defect in internalisation is observed. Mechanistically, excessive surface interleukin-7 receptor-α sequestrates the common gamma chain, impairing the interleukin-2 receptor assembly and downstream signalling crucial for T cell proliferation. This study shows that key autophagy substrates can be reliably identified in this mouse model and help mechanistically unravel autophagy's contribution to healthy physiology and disease.

An Integrated Taxonomy for Monogenic Inflammatory Bowel Disease.

BACKGROUND & AIMS: Monogenic forms of inflammatory bowel disease (IBD) illustrate the essential roles of individual genes in pathways and networks safeguarding immune tolerance and gut homeostasis. METHODS: To build a taxonomy model, we assessed 165 disorders. Genes were prioritized based on penetrance of IBD and disease phenotypes were integrated with multi-omics datasets. Monogenic IBD genes were classified by (1) overlapping syndromic features, (2) response to hematopoietic stem cell transplantation, (3) bulk RNA-sequencing of 32 tissues, (4) single-cell RNA-sequencing of >50 cell subsets from the intestine of healthy individuals and patients with IBD (pediatric and adult), and (5) proteomes of 43 immune subsets. The model was validated by addition of newly identified monogenic IBD defects. As a proof-of-concept, we explore the intersection between immunometabolism and antimicrobial activity for a group of disorders (G6PC3/SLC37A4). RESULTS: Our quantitative integrated taxonomy defines the cellular landscape of monogenic IBD gene expression across 102 genes with high and moderate penetrance (81 in the model set and 21 genes in the validation set). We illustrate distinct cellular networks, highlight expression profiles across understudied cell types (e.g., CD8+ T cells, neutrophils, epithelial subsets, and endothelial cells) and define genotype-phenotype associations (perianal disease and defective antimicrobial activity). We illustrate processes and pathways shared across cellular compartments and phenotypic groups and highlight cellular immunometabolism with mammalian target of rapamycin activation as one of the converging pathways. There is an overlap of genes and enriched cell-specific expression between monogenic and polygenic IBD. CONCLUSION: Our taxonomy integrates genetic, clinical and multi-omic data; providing a basis for genomic diagnostics and testable hypotheses for disease functions and treatment responses.

GWAS of stool frequency provides insights into gastrointestinal motility and irritable bowel syndrome.

Gut dysmotility is associated with constipation, diarrhea, and functional gastrointestinal disorders like irritable bowel syndrome (IBS), although its molecular underpinnings are poorly characterized. We studied stool frequency (defined by the number of bowel movements per day, based on questionnaire data) as a proxy for gut motility in a GWAS meta-analysis including 167,875 individuals from UK Biobank and four smaller population-based cohorts. We identify 14 loci associated with stool frequency (p ≤ 5.0 × 10-8). Gene set and pathway analyses detected enrichment for genes involved in neurotransmitter/neuropeptide signaling and preferentially expressed in enteric motor neurons controlling peristalsis. PheWAS identified pleiotropic associations with dysmotility syndromes and the response to their pharmacological treatment. The genetic architecture of stool frequency correlates with that of IBS, and UK Biobank participants from the top 1% of stool frequency polygenic score distribution were associated with 5× higher risk of IBS with diarrhea. These findings pave the way for the identification of actionable pathological mechanisms in IBS and the dysmotility syndromes.

Genome-wide analysis of 53,400 people with irritable bowel syndrome highlights shared genetic pathways with mood and anxiety disorders.

Irritable bowel syndrome (IBS) results from disordered brain-gut interactions. Identifying susceptibility genes could highlight the underlying pathophysiological mechanisms. We designed a digestive health questionnaire for UK Biobank and combined identified cases with IBS with independent cohorts. We conducted a genome-wide association study with 53,400 cases and 433,201 controls and replicated significant associations in a 23andMe panel (205,252 cases and 1,384,055 controls). Our study identified and confirmed six genetic susceptibility loci for IBS. Implicated genes included NCAM1, CADM2, PHF2/FAM120A, DOCK9, CKAP2/TPTE2P3 and BAG6. The first four are associated with mood and anxiety disorders, expressed in the nervous system, or both. Mirroring this, we also found strong genome-wide correlation between the risk of IBS and anxiety, neuroticism and depression (rg > 0.5). Additional analyses suggested this arises due to shared pathogenic pathways rather than, for example, anxiety causing abdominal symptoms. Implicated mechanisms require further exploration to help understand the altered brain-gut interactions underlying IBS.

Comparison of LABORAS with static incapacitance testing for assessing spontaneous pain behaviour in surgically-induced murine osteoarthritis.

OBJECTIVE: Evoked responses following mechanical or thermal stimulation are typically used to assess pain behaviour in murine osteoarthritis (OA). However, there is no consensus on how best to measure spontaneous pain behaviour. METHOD: OA by partial meniscectomy (PMX), or sham surgery was performed in 10-week old C57BL/6 male mice. Collagen-induced arthritis (CIA) was induced in 10 week old DBA1 male mice. Spontaneous pain behaviour, either at the time of active inflammatory disease (CIA), or over the 12 weeks after induction of OA, was assessed by static incapacitance testing (measuring percentage of weight placed through each hindlimb), and Laboratory Animal Behaviour Observation Registration and Analysis System (LABORAS) (translating cage vibrations of singly house animals into specific activities). Data were analysed by repeated measures two way ANOVA with post hoc testing comparing experimental groups with either sham operated or naïve controls. RESULTS: By incapacitance testing, two phases of painful behaviour were evident after PMX: a transient, post-operative phase, which resolved within one week, and a late OA pain phase starting 8 weeks post surgery and reaching statistical significance at week 12 (95% CI: sham 89.51-98.19, PMX 76.18-98.16). LABORAS, was able to detect pain behaviour in mice with CIA, but no statistically significant pain behaviour was observed in OA mice either post operatively (once analgesia had been controlled for) or at any later time points for any activity compared with the sham group. CONCLUSION: Static incapacitance testing is superior to LABORAS for measuring spontaneous pain behaviour in surgically induced murine OA.

Genomic profiling of T-cell activation suggests increased sensitivity of memory T cells to CD28 costimulation.

T-cell activation is a critical driver of immune responses. The CD28 costimulation is an essential regulator of CD4 T-cell responses, however, its relative importance in naive and memory T cells is not fully understood. Using different model systems, we observe that human memory T cells are more sensitive to CD28 costimulation than naive T cells. To deconvolute how the T-cell receptor (TCR) and CD28 orchestrate activation of human T cells, we stimulate cells using varying intensities of TCR and CD28 and profiled gene expression. We show that genes involved in cell cycle progression and division are CD28-driven in memory cells, but under TCR control in naive cells. We further demonstrate that T-helper differentiation and cytokine expression are controlled by CD28. Using chromatin accessibility profiling, we observe that AP1 transcriptional regulation is enriched when both TCR and CD28 are engaged, whereas open chromatin near CD28-sensitive genes is enriched for NF-kB motifs. Lastly, we show that CD28-sensitive genes are enriched in GWAS regions associated with immune diseases, implicating a role for CD28 in disease development. Our study provides important insights into the differential role of costimulation in naive and memory T-cell responses and disease susceptibility.

Graphical Model Selection for Gaussian Conditional Random Fields in the Presence of Latent Variables.

We consider the problem of learning a conditional Gaussian graphical model in the presence of latent variables. Building on recent advances in this field, we suggest a method that decomposes the parameters of a conditional Markov random field into the sum of a sparse and a low-rank matrix. We derive convergence bounds for this estimator and show that it is well-behaved in the high-dimensional regime as well as "sparsistent" (i.e., capable of recovering the graph structure). We then show how proximal gradient algorithms and semi-definite programming techniques can be employed to fit the model to thousands of variables. Through extensive simulations, we illustrate the conditions required for identifiability and show that there is a wide range of situations in which this model performs significantly better than its counterparts, for example, by accommodating more latent variables. Finally, the suggested method is applied to two datasets comprising individual level data on genetic variants and metabolites levels. We show our results replicate better than alternative approaches and show enriched biological signal. Supplementary materials for this article are available online.

Correction: Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

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

Bayesian meta-analysis across genome-wide association studies of diverse phenotypes.

Genome-wide association studies (GWAS) are a powerful tool for understanding the genetic basis of diseases and traits, but most studies have been conducted in isolation, with a focus on either a single or a set of closely related phenotypes. We describe MetABF, a simple Bayesian framework for performing integrative meta-analysis across multiple GWAS using summary statistics. The approach is applicable across a wide range of study designs and can increase the power by 50% compared with standard frequentist tests when only a subset of studies have a true effect. We demonstrate its utility in a meta-analysis of 20 diverse GWAS which were part of the Wellcome Trust Case Control Consortium 2. The novelty of the approach is its ability to explore, and assess the evidence for a range of possible true patterns of association across studies in a computationally efficient framework.

Active immunisation targeting nerve growth factor attenuates chronic pain behaviour in murine osteoarthritis.

OBJECTIVES: Nerve growth factor (NGF) has emerged as a key driver of pain in osteoarthritis (OA) and antibodies to NGF are potent analgesics in human disease. Here, we validate a novel vaccine strategy to generate anti-NGF antibodies for reversal of pain behaviour in a surgical model of OA. METHODS: Virus-like particles were derived from the cucumber mosaic virus (CuMV) and coupled to expressed recombinant NGF to create the vaccine. 10-week-old male mice underwent partial meniscectomy to induce OA or sham-surgery. Spontaneous pain behaviour was measured by Linton incapacitance and OA severity was quantified using OARSI histological scoring. Mice (experimental and a sentinel cohort) were inoculated with CuMVttNGF (Vax) or CuMVttctrl (Mock) either before surgery or once pain was established. Efficacy of anti-NGF from the plasma of sentinel vaccinated mice was measured in vitro using a neurite outgrowth assay in PC12 cells. RESULTS: Anti-NGF titres were readily detectable in the vaccinated but not mock vaccinated mice. Regular boosting with fresh vaccine was required to maintain anti-NGF titres as measured in the sentinel cohort. Both prophylactic and therapeutic vaccination demonstrated a reversal of pain behaviour by incapacitance testing, and a meta-analysis of the two studies showing analgesia at peak anti-NGF titres was highly statistically significant. Serum anti-NGF was able to inhibit neurite outgrowth equivalent to around 150 ug/mL of recombinant monoclonal antibody. CONCLUSIONS: This study demonstrates therapeutic efficacy of a novel NGF vaccine strategy that reversibly alleviates spontaneous pain behaviour in surgically induced murine OA.

Insights into the genetic epidemiology of Crohn's and rare diseases in the Ashkenazi Jewish population.

As part of a broader collaborative network of exome sequencing studies, we developed a jointly called data set of 5,685 Ashkenazi Jewish exomes. We make publicly available a resource of site and allele frequencies, which should serve as a reference for medical genetics in the Ashkenazim (hosted in part at https://ibd.broadinstitute.org, also available in gnomAD at http://gnomad.broadinstitute.org). We estimate that 34% of protein-coding alleles present in the Ashkenazi Jewish population at frequencies greater than 0.2% are significantly more frequent (mean 15-fold) than their maximum frequency observed in other reference populations. Arising via a well-described founder effect approximately 30 generations ago, this catalog of enriched alleles can contribute to differences in genetic risk and overall prevalence of diseases between populations. As validation we document 148 AJ enriched protein-altering alleles that overlap with "pathogenic" ClinVar alleles (table available at https://github.com/macarthur-lab/clinvar/blob/master/output/clinvar.tsv), including those that account for 10-100 fold differences in prevalence between AJ and non-AJ populations of some rare diseases, especially recessive conditions, including Gaucher disease (GBA, p.Asn409Ser, 8-fold enrichment); Canavan disease (ASPA, p.Glu285Ala, 12-fold enrichment); and Tay-Sachs disease (HEXA, c.1421+1G>C, 27-fold enrichment; p.Tyr427IlefsTer5, 12-fold enrichment). We next sought to use this catalog, of well-established relevance to Mendelian disease, to explore Crohn's disease, a common disease with an estimated two to four-fold excess prevalence in AJ. We specifically attempt to evaluate whether strong acting rare alleles, particularly protein-truncating or otherwise large effect-size alleles, enriched by the same founder-effect, contribute excess genetic risk to Crohn's disease in AJ, and find that ten rare genetic risk factors in NOD2 and LRRK2 are enriched in AJ (p < 0.005), including several novel contributing alleles, show evidence of association to CD. Independently, we find that genomewide common variant risk defined by GWAS shows a strong difference between AJ and non-AJ European control population samples (0.97 s.d. higher, p<10-16). Taken together, the results suggest coordinated selection in AJ population for higher CD risk alleles in general. The results and approach illustrate the value of exome sequencing data in case-control studies along with reference data sets like ExAC (sites VCF available via FTP at ftp.broadinstitute.org/pub/ExAC_release/release0.3/) to pinpoint genetic variation that contributes to variable disease predisposition across populations.

Bayesian analysis of genetic association across tree-structured routine healthcare data in the UK Biobank.

Genetic discovery from the multitude of phenotypes extractable from routine healthcare data can transform understanding of the human phenome and accelerate progress toward precision medicine. However, a critical question when analyzing high-dimensional and heterogeneous data is how best to interrogate increasingly specific subphenotypes while retaining statistical power to detect genetic associations. Here we develop and employ a new Bayesian analysis framework that exploits the hierarchical structure of diagnosis classifications to analyze genetic variants against UK Biobank disease phenotypes derived from self-reporting and hospital episode statistics. Our method displays a more than 20% increase in power to detect genetic effects over other approaches and identifies new associations between classical human leukocyte antigen (HLA) alleles and common immune-mediated diseases (IMDs). By applying the approach to genetic risk scores (GRSs), we show the extent of genetic sharing among IMDs and expose differences in disease perception or diagnosis with potential clinical implications.

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.

Exploring the genetic architecture of inflammatory bowel disease by whole-genome sequencing identifies association at ADCY7.

To further resolve the genetic architecture of the inflammatory bowel diseases ulcerative colitis and Crohn's disease, we sequenced the whole genomes of 4,280 patients at low coverage and compared them to 3,652 previously sequenced population controls across 73.5 million variants. We then imputed from these sequences into new and existing genome-wide association study cohorts and tested for association at ∼12 million variants in a total of 16,432 cases and 18,843 controls. We discovered a 0.6% frequency missense variant in ADCY7 that doubles the risk of ulcerative colitis. Despite good statistical power, we did not identify any other new low-frequency risk variants and found that such variants explained little heritability. We detected a burden of very rare, damaging missense variants in known Crohn's disease risk genes, suggesting that more comprehensive sequencing studies will continue to improve understanding of the biology of complex diseases.

Genome-wide association study implicates immune activation of multiple integrin genes in inflammatory bowel disease.

Genetic association studies have identified 215 risk loci for inflammatory bowel disease, thereby uncovering fundamental aspects of its molecular biology. We performed a genome-wide association study of 25,305 individuals and conducted a meta-analysis with published summary statistics, yielding a total sample size of 59,957 subjects. We identified 25 new susceptibility loci, 3 of which contain integrin genes that encode proteins in pathways that have been identified as important therapeutic targets in inflammatory bowel disease. The associated variants are correlated with expression changes in response to immune stimulus at two of these genes (ITGA4 and ITGB8) and at previously implicated loci (ITGAL and ICAM1). In all four cases, the expression-increasing allele also increases disease risk. We also identified likely causal missense variants in a gene implicated in primary immune deficiency, PLCG2, and a negative regulator of inflammation, SLAMF8. Our results demonstrate that new associations at common variants continue to identify genes relevant to therapeutic target identification and prioritization.

Genome-wide association study identifies distinct genetic contributions to prognosis and susceptibility in Crohn's disease.

For most immune-mediated diseases, the main determinant of patient well-being is not the diagnosis itself but instead the course that the disease takes over time (prognosis). Prognosis may vary substantially between patients for reasons that are poorly understood. Familial studies support a genetic contribution to prognosis, but little evidence has been found for a proposed association between prognosis and the burden of susceptibility variants. To better characterize how genetic variation influences disease prognosis, we performed a within-cases genome-wide association study in two cohorts of patients with Crohn's disease. We identified four genome-wide significant loci, none of which showed any association with disease susceptibility. Conversely, the aggregated effect of all 170 disease susceptibility loci was not associated with disease prognosis. Together, these data suggest that the genetic contribution to prognosis in Crohn's disease is largely independent of the contribution to disease susceptibility and point to a biology of prognosis that could provide new therapeutic opportunities.

Genome-wide association study of primary sclerosing cholangitis identifies new risk loci and quantifies the genetic relationship with inflammatory bowel disease.

Primary sclerosing cholangitis (PSC) is a rare progressive disorder leading to bile duct destruction; ∼75% of patients have comorbid inflammatory bowel disease (IBD). We undertook the largest genome-wide association study of PSC (4,796 cases and 19,955 population controls) and identified four new genome-wide significant loci. The most associated SNP at one locus affects splicing and expression of UBASH3A, with the protective allele (C) predicted to cause nonstop-mediated mRNA decay and lower expression of UBASH3A. Further analyses based on common variants suggested that the genome-wide genetic correlation (rG) between PSC and ulcerative colitis (UC) (rG = 0.29) was significantly greater than that between PSC and Crohn's disease (CD) (rG = 0.04) (P = 2.55 × 10-15). UC and CD were genetically more similar to each other (rG = 0.56) than either was to PSC (P < 1.0 × 10-15). Our study represents a substantial advance in understanding of the genetics of PSC.

Resolving TYK2 locus genotype-to-phenotype differences in autoimmunity.

Thousands of genetic variants have been identified, which contribute to the development of complex diseases, but determining how to elucidate their biological consequences for translation into clinical benefit is challenging. Conflicting evidence regarding the functional impact of genetic variants in the tyrosine kinase 2 (TYK2) gene, which is differentially associated with common autoimmune diseases, currently obscures the potential of TYK2 as a therapeutic target. We aimed to resolve this conflict by performing genetic meta-analysis across disorders; subsequent molecular, cellular, in vivo, and structural functional follow-up; and epidemiological studies. Our data revealed a protective homozygous effect that defined a signaling optimum between autoimmunity and immunodeficiency and identified TYK2 as a potential drug target for certain common autoimmune disorders.

Genetic Complexity of Crohn's Disease in Two Large Ashkenazi Jewish Families.

BACKGROUND & AIMS: Crohn's disease (CD) is a highly heritable disease that is particularly common in the Ashkenazi Jewish population. We studied 2 large Ashkenazi Jewish families with a high prevalence of CD in an attempt to identify novel genetic risk variants. METHODS: Ashkenazi Jewish patients with CD and a positive family history were recruited from the University College London Hospital. We used genome-wide, single-nucleotide polymorphism data to assess the burden of common CD-associated risk variants and for linkage analysis. Exome sequencing was performed and rare variants that were predicted to be deleterious and were observed at a high frequency in cases were prioritized. We undertook within-family association analysis after imputation and assessed candidate variants for evidence of association with CD in an independent cohort of Ashkenazi Jewish individuals. We examined the effects of a variant in DUOX2 on hydrogen peroxide production in HEK293 cells. RESULTS: We identified 2 families (1 with >800 members and 1 with >200 members) containing 54 and 26 cases of CD or colitis, respectively. Both families had a significant enrichment of previously described common CD-associated risk variants. No genome-wide significant linkage was observed. Exome sequencing identified candidate variants, including a missense mutation in DUOX2 that impaired its function and a frameshift mutation in CSF2RB that was associated with CD in an independent cohort of Ashkenazi Jewish individuals. CONCLUSIONS: In a study of 2 large Ashkenazi Jewish with multiple cases of CD, we found the genetic basis of the disease to be complex, with a role for common and rare genetic variants. We identified a frameshift mutation in CSF2RB that was replicated in an independent cohort. These findings show the value of family studies and the importance of the innate immune system in the pathogenesis of CD.

Analysis of five chronic inflammatory diseases identifies 27 new associations and highlights disease-specific patterns at shared loci.

We simultaneously investigated the genetic landscape of ankylosing spondylitis, Crohn's disease, psoriasis, primary sclerosing cholangitis and ulcerative colitis to investigate pleiotropy and the relationship between these clinically related diseases. Using high-density genotype data from more than 86,000 individuals of European ancestry, we identified 244 independent multidisease signals, including 27 new genome-wide significant susceptibility loci and 3 unreported shared risk loci. Complex pleiotropy was supported when contrasting multidisease signals with expression data sets from human, rat and mouse together with epigenetic and expressed enhancer profiles. The comorbidities among the five immune diseases were best explained by biological pleiotropy rather than heterogeneity (a subgroup of cases genetically identical to those with another disease, possibly owing to diagnostic misclassification, molecular subtypes or excessive comorbidity). In particular, the strong comorbidity between primary sclerosing cholangitis and inflammatory bowel disease is likely the result of a unique disease, which is genetically distinct from classical inflammatory bowel disease phenotypes.