Genetics of circulating inflammatory proteins identifies drivers of immune-mediated disease risk and therapeutic targets.

Circulating proteins have important functions in inflammation and a broad range of diseases. To identify genetic influences on inflammation-related proteins, we conducted a genome-wide protein quantitative trait locus (pQTL) study of 91 plasma proteins measured using the Olink Target platform in 14,824 participants. We identified 180 pQTLs (59 cis, 121 trans). Integration of pQTL data with eQTL and disease genome-wide association studies provided insight into pathogenesis, implicating lymphotoxin-α in multiple sclerosis. Using Mendelian randomization (MR) to assess causality in disease etiology, we identified both shared and distinct effects of specific proteins across immune-mediated diseases, including directionally discordant effects of CD40 on risk of rheumatoid arthritis versus multiple sclerosis and inflammatory bowel disease. MR implicated CXCL5 in the etiology of ulcerative colitis (UC) and we show elevated gut CXCL5 transcript expression in patients with UC. These results identify targets of existing drugs and provide a powerful resource to facilitate future drug target prioritization.

Interleukin-22 orchestrates a pathological endoplasmic reticulum stress response transcriptional programme in colonic epithelial cells.

OBJECTIVE: The functional role of interleukin-22 (IL22) in chronic inflammation is controversial, and mechanistic insights into how it regulates target tissue are lacking. In this study, we evaluated the functional role of IL22 in chronic colitis and probed mechanisms of IL22-mediated regulation of colonic epithelial cells. DESIGN: To investigate the functional role of IL22 in chronic colitis and how it regulates colonic epithelial cells, we employed a three-dimentional mini-gut epithelial organoid system, in vivo disease models and transcriptomic datasets in human IBD. RESULTS: As well as inducing transcriptional modules implicated in antimicrobial responses, IL22 also coordinated an endoplasmic reticulum (ER) stress response transcriptional programme in colonic epithelial cells. In the colon of patients with active colonic Crohn's disease (CD), there was enrichment of IL22-responsive transcriptional modules and ER stress response modules. Strikingly, in an IL22-dependent model of chronic colitis, targeting IL22 alleviated colonic epithelial ER stress and attenuated colitis. Pharmacological modulation of the ER stress response similarly impacted the severity of colitis. In patients with colonic CD, antibody blockade of IL12p40, which simultaneously blocks IL12 and IL23, the key upstream regulator of IL22 production, alleviated the colonic epithelial ER stress response. CONCLUSIONS: Our data challenge perceptions of IL22 as a predominantly beneficial cytokine in IBD and provide novel insights into the molecular mechanisms of IL22-mediated pathogenicity in chronic colitis. Targeting IL22-regulated pathways and alleviating colonic epithelial ER stress may represent promising therapeutic strategies in patients with colitis. TRIAL REGISTRATION NUMBER: NCT02749630.

Genome3D: exploiting structure to help users understand their sequences.

Genome3D (http://www.genome3d.eu) is a collaborative resource that provides predicted domain annotations and structural models for key sequences. Since introducing Genome3D in a previous NAR paper, we have substantially extended and improved the resource. We have annotated representatives from Pfam families to improve coverage of diverse sequences and added a fast sequence search to the website to allow users to find Genome3D-annotated sequences similar to their own. We have improved and extended the Genome3D data, enlarging the source data set from three model organisms to 10, and adding VIVACE, a resource new to Genome3D. We have analysed and updated Genome3D's SCOP/CATH mapping. Finally, we have improved the superposition tools, which now give users a more powerful interface for investigating similarities and differences between structural models.

A large-scale evaluation of computational protein function prediction.

Automated annotation of protein function is challenging. As the number of sequenced genomes rapidly grows, the overwhelming majority of protein products can only be annotated computationally. If computational predictions are to be relied upon, it is crucial that the accuracy of these methods be high. Here we report the results from the first large-scale community-based critical assessment of protein function annotation (CAFA) experiment. Fifty-four methods representing the state of the art for protein function prediction were evaluated on a target set of 866 proteins from 11 organisms. Two findings stand out: (i) today's best protein function prediction algorithms substantially outperform widely used first-generation methods, with large gains on all types of targets; and (ii) although the top methods perform well enough to guide experiments, there is considerable need for improvement of currently available tools.

DISOPRED3: precise disordered region predictions with annotated protein-binding activity.

MOTIVATION: A sizeable fraction of eukaryotic proteins contain intrinsically disordered regions (IDRs), which act in unfolded states or by undergoing transitions between structured and unstructured conformations. Over time, sequence-based classifiers of IDRs have become fairly accurate and currently a major challenge is linking IDRs to their biological roles from the molecular to the systems level. RESULTS: We describe DISOPRED3, which extends its predecessor with new modules to predict IDRs and protein-binding sites within them. Based on recent CASP evaluation results, DISOPRED3 can be regarded as state of the art in the identification of IDRs, and our self-assessment shows that it significantly improves over DISOPRED2 because its predictions are more specific across the whole board and more sensitive to IDRs longer than 20 amino acids. Predicted IDRs are annotated as protein binding through a novel SVM based classifier, which uses profile data and additional sequence-derived features. Based on benchmarking experiments with full cross-validation, we show that this predictor generates precise assignments of disordered protein binding regions and that it compares well with other publicly available tools.

Evaluation of predictions in the CASP10 model refinement category.

Here we report on the assessment results of the third experiment to evaluate the state of the art in protein model refinement, where participants were invited to improve the accuracy of initial protein models for 27 targets. Using an array of complementary evaluation measures, we find that five groups performed better than the naïve (null) method-a marked improvement over CASP9, although only three were significantly better. The leading groups also demonstrated the ability to consistently improve both backbone and side chain positioning, while other groups reliably enhanced other aspects of protein physicality. The top-ranked group succeeded in improving the backbone conformation in almost 90% of targets, suggesting a strategy that for the first time in CASP refinement is successful in a clear majority of cases. A number of issues remain unsolved: the majority of groups still fail to improve the quality of the starting models; even successful groups are only able to make modest improvements; and no prediction is more similar to the native structure than to the starting model. Successful refinement attempts also often go unrecognized, as suggested by the relatively larger improvements when predictions not submitted as model 1 are also considered.

Protein function prediction by massive integration of evolutionary analyses and multiple data sources.

BACKGROUND: Accurate protein function annotation is a severe bottleneck when utilizing the deluge of high-throughput, next generation sequencing data. Keeping database annotations up-to-date has become a major scientific challenge that requires the development of reliable automatic predictors of protein function. The CAFA experiment provided a unique opportunity to undertake comprehensive 'blind testing' of many diverse approaches for automated function prediction. We report on the methodology we used for this challenge and on the lessons we learnt. METHODS: Our method integrates into a single framework a wide variety of biological information sources, encompassing sequence, gene expression and protein-protein interaction data, as well as annotations in UniProt entries. The methodology transfers functional categories based on the results from complementary homology-based and feature-based analyses. We generated the final molecular function and biological process assignments by combining the initial predictions in a probabilistic manner, which takes into account the Gene Ontology hierarchical structure. RESULTS: We propose a novel scoring function called COmbined Graph-Information Content similarity (COGIC) score for the comparison of predicted functional categories and benchmark data. We demonstrate that our integrative approach provides increased scope and accuracy over both the component methods and the naïve predictors. In line with previous studies, we find that molecular function predictions are more accurate than biological process assignments. CONCLUSIONS: Overall, the results indicate that there is considerable room for improvement in the field. It still remains for the community to invest a great deal of effort to make automated function prediction a useful and routine component in the toolbox of life scientists. As already witnessed in other areas, community-wide blind testing experiments will be pivotal in establishing standards for the evaluation of prediction accuracy, in fostering advancements and new ideas, and ultimately in recording progress.

PSICOV: precise structural contact prediction using sparse inverse covariance estimation on large multiple sequence alignments.

MOTIVATION: The accurate prediction of residue-residue contacts, critical for maintaining the native fold of a protein, remains an open problem in the field of structural bioinformatics. Interest in this long-standing problem has increased recently with algorithmic improvements and the rapid growth in the sizes of sequence families. Progress could have major impacts in both structure and function prediction to name but two benefits. Sequence-based contact predictions are usually made by identifying correlated mutations within multiple sequence alignments (MSAs), most commonly through the information-theoretic approach of calculating mutual information between pairs of sites in proteins. These predictions are often inaccurate because the true covariation signal in the MSA is often masked by biases from many ancillary indirect-coupling or phylogenetic effects. Here we present a novel method, PSICOV, which introduces the use of sparse inverse covariance estimation to the problem of protein contact prediction. Our method builds on work which had previously demonstrated corrections for phylogenetic and entropic correlation noise and allows accurate discrimination of direct from indirectly coupled mutation correlations in the MSA. RESULTS: PSICOV displays a mean precision substantially better than the best performing normalized mutual information approach and Bayesian networks. For 118 out of 150 targets, the L/5 (i.e. top-L/5 predictions for a protein of length L) precision for long-range contacts (sequence separation >23) was ≥ 0.5, which represents an improvement sufficient to be of significant benefit in protein structure prediction or model quality assessment. AVAILABILITY: The PSICOV source code can be downloaded from http://bioinf.cs.ucl.ac.uk/downloads/PSICOV.

Immune checkpoint inhibitor-induced colitis is mediated by polyfunctional lymphocytes and is dependent on an IL23/IFNγ axis.

Immune checkpoint inhibitors (CPIs) are a relatively newly licenced cancer treatment, which make a once previously untreatable disease now amenable to a potential cure. Combination regimens of anti-CTLA4 and anti-PD-1 show enhanced efficacy but are prone to off-target immune-mediated tissue injury, particularly at the barrier surfaces. To probe the impact of immune checkpoints on intestinal homoeostasis, mice are challenged with anti-CTLA4 and anti-PD-1 immunotherapy and manipulation of the intestinal microbiota. The immune profile of the colon of these mice with CPI-colitis is analysed using bulk RNA sequencing, single-cell RNA sequencing and flow cytometry. CPI-colitis in mice is dependent on the composition of the intestinal microbiota and by the induction of lymphocytes expressing interferon-γ (IFNγ), cytotoxicity molecules and other pro-inflammatory cytokines/chemokines. This pre-clinical model of CPI-colitis could be attenuated following blockade of the IL23/IFNγ axis. Therapeutic targeting of IFNγ-producing lymphocytes or regulatory networks, may hold the key to reversing CPI-colitis.

Tumor-activated lymph node fibroblasts suppress T cell function in diffuse large B cell lymphoma.

Recent transcriptomic-based analysis of diffuse large B cell lymphoma (DLBCL) has highlighted the clinical relevance of LN fibroblast and tumor-infiltrating lymphocyte (TIL) signatures within the tumor microenvironment (TME). However, the immunomodulatory role of fibroblasts in lymphoma remains unclear. Here, by studying human and mouse DLBCL-LNs, we identified the presence of an aberrantly remodeled fibroblastic reticular cell (FRC) network expressing elevated fibroblast-activated protein (FAP). RNA-Seq analyses revealed that exposure to DLBCL reprogrammed key immunoregulatory pathways in FRCs, including a switch from homeostatic to inflammatory chemokine expression and elevated antigen-presentation molecules. Functional assays showed that DLBCL-activated FRCs (DLBCL-FRCs) hindered optimal TIL and chimeric antigen receptor (CAR) T cell migration. Moreover, DLBCL-FRCs inhibited CD8+ TIL cytotoxicity in an antigen-specific manner. Notably, the interrogation of patient LNs with imaging mass cytometry identified distinct environments differing in their CD8+ TIL-FRC composition and spatial organization that associated with survival outcomes. We further demonstrated the potential to target inhibitory FRCs to rejuvenate interacting TILs. Cotreating organotypic cultures with FAP-targeted immunostimulatory drugs and a bispecific antibody (glofitamab) augmented antilymphoma TIL cytotoxicity. Our study reveals an immunosuppressive role of FRCs in DLBCL, with implications for immune evasion, disease pathogenesis, and optimizing immunotherapy for patients.

Cytokine responsive networks in human colonic epithelial organoids unveil a molecular classification of inflammatory bowel disease.

Interactions between the epithelium and the immune system are critical in the pathogenesis of inflammatory bowel disease (IBD). In this study, we mapped the transcriptional landscape of human colonic epithelial organoids in response to different cytokines responsible for mediating canonical mucosal immune responses. By profiling the transcriptome of human colonic organoids treated with the canonical cytokines interferon gamma, interleukin-13, -17A, and tumor necrosis factor alpha with next-generation sequencing, we unveil shared and distinct regulation patterns of epithelial function by different cytokines. An integrative analysis of cytokine responses in diseased tissue from patients with IBD (n = 1,009) reveals a molecular classification of mucosal inflammation defined by gradients of cytokine-responsive transcriptional signatures. Our systems biology approach detected signaling bottlenecks in cytokine-responsive networks and highlighted their translational potential as theragnostic targets in intestinal inflammation.

Interleukin-22 regulates neutrophil recruitment in ulcerative colitis and is associated with resistance to ustekinumab therapy.

The function of interleukin-22 (IL-22) in intestinal barrier homeostasis remains controversial. Here, we map the transcriptional landscape regulated by IL-22 in human colonic epithelial organoids and evaluate the biological, functional and clinical significance of the IL-22 mediated pathways in ulcerative colitis (UC). We show that IL-22 regulated pro-inflammatory pathways are involved in microbial recognition, cancer and immune cell chemotaxis; most prominently those involving CXCR2+ neutrophils. IL-22-mediated transcriptional regulation of CXC-family neutrophil-active chemokine expression is highly conserved across species, is dependent on STAT3 signaling, and is functionally and pathologically important in the recruitment of CXCR2+ neutrophils into colonic tissue. In UC patients, the magnitude of enrichment of the IL-22 regulated transcripts in colonic biopsies correlates with colonic neutrophil infiltration and is enriched in non-responders to ustekinumab therapy. Our data provide further insights into the biology of IL-22 in human disease and highlight its function in the regulation of pathogenic immune pathways, including neutrophil chemotaxis. The transcriptional networks regulated by IL-22 are functionally and clinically important in UC, impacting patient trajectories and responsiveness to biological intervention.

TGF-ß1 potentiates Vγ9Vδ2 T cell adoptive immunotherapy of cancer.

Despite its role in cancer surveillance, adoptive immunotherapy using γδ T cells has achieved limited efficacy. To enhance trafficking to bone marrow, circulating Vγ9Vδ2 T cells are expanded in serum-free medium containing TGF-ß1 and IL-2 (γδ[T2] cells) or medium containing IL-2 alone (γδ[2] cells, as the control). Unexpectedly, the yield and viability of γδ[T2] cells are also increased by TGF-ß1, when compared to γδ[2] controls. γδ[T2] cells are less differentiated and yet display increased cytolytic activity, cytokine release, and antitumor activity in several leukemic and solid tumor models. Efficacy is further enhanced by cancer cell sensitization using aminobisphosphonates or Ara-C. A number of contributory effects of TGF-ß are described, including prostaglandin E2 receptor downmodulation, TGF-ß insensitivity, and upregulated integrin activity. Biological relevance is supported by the identification of a favorable γδ[T2] signature in acute myeloid leukemia (AML). Given their enhanced therapeutic activity and compatibility with allogeneic use, γδ[T2] cells warrant evaluation in cancer immunotherapy.

Cancer-associated hypersialylated MUC1 drives the differentiation of human monocytes into macrophages with a pathogenic phenotype.

The tumour microenvironment plays a crucial role in the growth and progression of cancer, and the presence of tumour-associated macrophages (TAMs) is associated with poor prognosis. Recent studies have demonstrated that TAMs display transcriptomic, phenotypic, functional and geographical diversity. Here we show that a sialylated tumour-associated glycoform of the mucin MUC1, MUC1-ST, through the engagement of Siglec-9 can specifically and independently induce the differentiation of monocytes into TAMs with a unique phenotype that to the best of our knowledge has not previously been described. These TAMs can recruit and prolong the lifespan of neutrophils, inhibit the function of T cells, degrade basement membrane allowing for invasion, are inefficient at phagocytosis, and can induce plasma clotting. This macrophage phenotype is enriched in the stroma at the edge of breast cancer nests and their presence is associated with poor prognosis in breast cancer patients.

Evaluation of CASP8 model quality predictions.

The model quality assessment problem consists in the a priori estimation of the overall and per-residue accuracy of protein structure predictions. Over the past years, a number of methods have been developed to address this issue and CASP established a prediction category to evaluate their performance in 2006. In 2008 the experiment was repeated and its results are reported here. Participants were invited to infer the correctness of the protein models submitted by the registered automatic servers. Estimates could apply to both whole models and individual amino acids. Groups involved in the tertiary structure prediction categories were also asked to assign local error estimates to each predicted residue in their own models and their results are also discussed here. The correlation between the predicted and observed correctness measures was the basis of the assessment of the results. We observe that consensus-based methods still perform significantly better than those accepting single models, similarly to what was concluded in the previous edition of the experiment.

Evaluation of template-based models in CASP8 with standard measures.

The strategy for evaluating template-based models submitted to CASP has continuously evolved from CASP1 to CASP5, leading to a standard procedure that has been used in all subsequent editions. The established approach includes methods for calculating the quality of each individual model, for assigning scores based on the distribution of the results for each target and for computing the statistical significance of the differences in scores between prediction methods. These data are made available to the assessor of the template-based modeling category, who uses them as a starting point for further evaluations and analyses. This article describes the detailed workflow of the procedure, provides justifications for a number of choices that are customarily made for CASP data evaluation, and reports the results of the analysis of template-based predictions at CASP8.

Using deep maxout neural networks to improve the accuracy of function prediction from protein interaction networks.

Protein-protein interaction network data provides valuable information that infers direct links between genes and their biological roles. This information brings a fundamental hypothesis for protein function prediction that interacting proteins tend to have similar functions. With the help of recently-developed network embedding feature generation methods and deep maxout neural networks, it is possible to extract functional representations that encode direct links between protein-protein interactions information and protein function. Our novel method, STRING2GO, successfully adopts deep maxout neural networks to learn functional representations simultaneously encoding both protein-protein interactions and functional predictive information. The experimental results show that STRING2GO outperforms other protein-protein interaction network-based prediction methods and one benchmark method adopted in a recent large scale protein function prediction competition.

The evaluation of protein structure prediction results.

Methods for protein structure prediction are flourishing and becoming widely available to both experimentalists and computational biologists. However, how good are they? What is their range of applicability and how can we know which method is better suited for the task at hand? These are the questions that this review tries to address, by describing the worldwide Critical Assessment of techniques for protein Structure Prediction (CASP) initiative and focusing on the specific problems of assessing the quality of a protein 3D model.