Next-Generation Immunosequencing Reveals Pathological T-Cell Architecture in Autoimmune Hepatitis.

BACKGROUND AND AIMS: Autoimmune hepatitis (AIH) is a chronic liver disease that regularly relapses when immunosuppression is tapered. It is thought to be driven by T-cells, whereas the etiologic impact of an apparently deregulated B lineage system, as evidenced by hypergammaglobulinemia and autoantibodies, remains elusive. We set out to investigate T and B cell repertoires supporting chronic inflammation in AIH. APPROACH AND RESULTS: T and B cell receptor (TCR/BCR) and human leukocyte antigen (HLA) next-generation immunosequencing were used to record immune signatures from a cohort of 60 patients with AIH and disease controls. Blood and liver B lineage immune metrics were not indicative of a dominant directional antigen selection apart from a slight skewing of IGHV-J genes. More importantly, we found strong AIH-specific TRBV-J skewing not attributable to the HLA-DRB1 specificities of the cohort. This TCR repertoire bias was generated as a result of peripheral T cell (de)selection and persisted in disease remission. Using a clustering algorithm according to antigenic specificity, we identified liver TCR clusters that were shared between patients with AIH but were absent or deselected in patients with other liver pathologies. CONCLUSIONS: Patients with AIH show profound and persisting T-cell architectural changes that may explain high relapse rates after tapering immunosuppression. Liver T-cell clusters shared between patients may mediate liver damage and warrant further study.

Deep sequencing of bone marrow microenvironments of patients with del(5q) myelodysplastic syndrome reveals imprints of antigenic selection as well as generation of novel T-cell clusters as a response pattern to lenalidomide.

In myelodysplastic syndromes with a partial deletion of the long arm of chromosome 5, del(5q), lenalidomide is believed to reverse anergic T-cell immunity in the bone marrow resulting in suppression of the del(5q) clone. In this study we used next-generation sequencing of immunoglobulin heavy chain (IGH) and T-cell receptor beta (TRB) rearrangements in bone marrow-residing and peripheral blood-circulating lymphocytes of patients with del(5q) myelodysplastic syndromes to assess the immune architecture and track adaptive immune responses during treatment with lenalidomide. The baseline bone marrow B-cell space in patients was comparable to that of age-matched healthy controls in terms of gene usage and IGH clonality, but showed a higher percentage of hypermutated IGH sequences, indicating an expanded number of antigen-experienced B lineage cells. Bone marrow B lineage clonality decreased significantly and hypermutated IGH clones normalized upon lenalidomide treatment, well in line with the proliferative effect on healthy antigen-inexperienced B-cell precursors previously described for this drug. The T-cell space in bone marrow of patients with del(5q) myelodysplastic syndromes showed higher TRB clonality compared to that of healthy controls. Upon lenalidomide treatment, myelodysplastic syndrome-specific T-cell clusters with low to medium spontaneous generation probabilities emerged; these clusters were shared across patients, indicating a common antigen-driven T-cell response pattern. Hence, we observed B lineage diversification and generation of new, antigen-dependent T-cell clusters, compatible with a model of adaptive immunity induced against the del(5q) clone by lenalidomide. Overall, this supports the concept that lenalidomide not only alters the functional T-cell state, but also the composition of the T- and B-cell repertoires in del(5q) myelodysplastic syndromes.

ConBind: motif-aware cross-species alignment for the identification of functional transcription factor binding sites.

Eukaryotic gene expression is regulated by transcription factors (TFs) binding to promoter as well as distal enhancers. TFs recognize short, but specific binding sites (TFBSs) that are located within the promoter and enhancer regions. Functionally relevant TFBSs are often highly conserved during evolution leaving a strong phylogenetic signal. While multiple sequence alignment (MSA) is a potent tool to detect the phylogenetic signal, the current MSA implementations are optimized to align the maximum number of identical nucleotides. This approach might result in the omission of conserved motifs that contain interchangeable nucleotides such as the ETS motif (IUPAC code: GGAW). Here, we introduce ConBind, a novel method to enhance alignment of short motifs, even if their mutual sequence similarity is only partial. ConBind improves the identification of conserved TFBSs by improving the alignment accuracy of TFBS families within orthologous DNA sequences. Functional validation of the Gfi1b + 13 enhancer reveals that ConBind identifies additional functionally important ETS binding sites that were missed by all other tested alignment tools. In addition to the analysis of known regulatory regions, our web tool is useful for the analysis of TFBSs on so far unknown DNA regions identified through ChIP-sequencing.

BioASF: a framework for automatically generating executable pathway models specified in BioPAX.

MOTIVATION: Biological pathways play a key role in most cellular functions. To better understand these functions, diverse computational and cell biology researchers use biological pathway data for various analysis and modeling purposes. For specifying these biological pathways, a community of researchers has defined BioPAX and provided various tools for creating, validating and visualizing BioPAX models. However, a generic software framework for simulating BioPAX models is missing. Here, we attempt to fill this gap by introducing a generic simulation framework for BioPAX. The framework explicitly separates the execution model from the model structure as provided by BioPAX, with the advantage that the modelling process becomes more reproducible and intrinsically more modular; this ensures natural biological constraints are satisfied upon execution. The framework is based on the principles of discrete event systems and multi-agent systems, and is capable of automatically generating a hierarchical multi-agent system for a given BioPAX model. RESULTS: To demonstrate the applicability of the framework, we simulated two types of biological network models: a gene regulatory network modeling the haematopoietic stem cell regulators and a signal transduction network modeling the Wnt/ß-catenin signaling pathway. We observed that the results of the simulations performed using our framework were entirely consistent with the simulation results reported by the researchers who developed the original models in a proprietary language. AVAILABILITY AND IMPLEMENTATION: The framework, implemented in Java, is open source and its source code, documentation and tutorial are available at http://www.ibi.vu.nl/programs/BioASF CONTACT: j.heringa@vu.nl.

Hard-wired heterogeneity in blood stem cells revealed using a dynamic regulatory network model.

MOTIVATION: Combinatorial interactions of transcription factors with cis-regulatory elements control the dynamic progression through successive cellular states and thus underpin all metazoan development. The construction of network models of cis-regulatory elements, therefore, has the potential to generate fundamental insights into cellular fate and differentiation. Haematopoiesis has long served as a model system to study mammalian differentiation, yet modelling based on experimentally informed cis-regulatory interactions has so far been restricted to pairs of interacting factors. Here, we have generated a Boolean network model based on detailed cis-regulatory functional data connecting 11 haematopoietic stem/progenitor cell (HSPC) regulator genes. RESULTS: Despite its apparent simplicity, the model exhibits surprisingly complex behaviour that we charted using strongly connected components and shortest-path analysis in its Boolean state space. This analysis of our model predicts that HSPCs display heterogeneous expression patterns and possess many intermediate states that can act as 'stepping stones' for the HSPC to achieve a final differentiated state. Importantly, an external perturbation or 'trigger' is required to exit the stem cell state, with distinct triggers characterizing maturation into the various different lineages. By focusing on intermediate states occurring during erythrocyte differentiation, from our model we predicted a novel negative regulation of Fli1 by Gata1, which we confirmed experimentally thus validating our model. In conclusion, we demonstrate that an advanced mammalian regulatory network model based on experimentally validated cis-regulatory interactions has allowed us to make novel, experimentally testable hypotheses about transcriptional mechanisms that control differentiation of mammalian stem cells. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

The role of proteosome-mediated proteolysis in modulating potentially harmful transcription factor activity in Saccharomyces cerevisiae.

MOTIVATION: The appropriate modulation of the stress response to variable environmental conditions is necessary to maintain sustained viability in Saccharomyces cerevisiae. Particularly, controlling the abundance of proteins that may have detrimental effects on cell growth is crucial for rapid recovery from stress-induced quiescence. RESULTS: Prompted by qualitative modeling of the nutrient starvation response in yeast, we investigated in vivo the effect of proteolysis after nutrient starvation showing that, for the Gis1 transcription factor at least, proteasome-mediated control is crucial for a rapid return to growth. Additional bioinformatics analyses show that potentially toxic transcriptional regulators have a significantly lower protein half-life, a higher fraction of unstructured regions and more potential PEST motifs than the non-detrimental ones. Furthermore, inhibiting proteasome activity tends to increase the expression of genes induced during the Environmental Stress Response more than those in the rest of the genome. Our combined results suggest that proteasome-mediated proteolysis of potentially toxic transcription factors tightly modulates the stress response in yeast. CONTACT: jasmin.fisher@microsoft.com

Executing multicellular differentiation: quantitative predictive modelling of C.elegans vulval development.

MOTIVATION: Understanding the processes involved in multi-cellular pattern formation is a central problem of developmental biology, hopefully leading to many new insights, e.g. in the treatment of various diseases. Defining suitable computational techniques for development modelling, able to perform in silico simulation experiments, is an open and challenging problem. RESULTS: Previously, we proposed a coarse-grained, quantitative approach based on the basic Petri net formalism, to mimic the behaviour of the biological processes during multicellular differentiation. Here, we apply our modelling approach to the well-studied process of Caenorhabditis elegans vulval development. We show that our model correctly reproduces a large set of in vivo experiments with statistical accuracy. It also generates gene expression time series in accordance with recent biological evidence. Finally, we modelled the role of microRNA mir-61 during vulval development and predict its contribution in stabilizing cell pattern formation.

Response of metastatic mouse invasive lobular carcinoma to mTOR inhibition is partly mediated by the adaptive immune system.

Effective treatment of invasive lobular carcinoma (ILC) of the breast is hampered by late detection, invasive growth, distant metastasis, and poor response to chemotherapy. Phosphoinositide 3-kinase (PI3K) signaling, one of the major druggable oncogenic signaling networks, is frequently activated in ILC. We investigated treatment response and resistance to AZD8055, an inhibitor of mammalian target of rapamycin (mTOR), in the K14-cre;Cdh1Flox/Flox;Trp53Flox/Flox (KEP) mouse model of metastatic ILC. Inhibition of mTOR signaling blocked the growth of primary KEP tumors as well as the progression of metastatic disease. However, primary tumors and distant metastases eventually acquired resistance after long-term AZD8055 treatment, despite continued effective suppression of mTOR signaling in cancer cells. Interestingly, therapeutic responses were associated with increased expression of genes related to antigen presentation. Consistent with this observation, increased numbers of tumor-infiltrating major histocompatibility complex class II-positive (MHCII+) immune cells were observed in treatment-responsive KEP tumors. Acquisition of treatment resistance was associated with loss of MHCII+ cells and reduced expression of genes related to the adaptive immune system. The therapeutic efficacy of mTOR inhibition was reduced in Rag1-/- mice lacking mature T and B lymphocytes, compared to immunocompetent mice. Furthermore, therapy responsiveness could be partially rescued by transplanting AZD8055-resistant KEP tumors into treatment-naïve immunocompetent hosts. Collectively, these data indicate that the PI3K signaling pathway is an attractive therapeutic target in invasive lobular carcinoma, and that part of the therapeutic effect of mTOR inhibition is mediated by the adaptive immune system.

High-Throughput Immunogenetics Reveals a Lack of Physiological T Cell Clusters in Patients With Autoimmune Cytopenias.

Autoimmune cytopenias (AIC) such as immune thrombocytopenia or autoimmune hemolytic anemia are claimed to be essentially driven by a dysregulated immune system. Using next-generation immunosequencing we profiled 59 T and B cell repertoires (TRB and IGH) of 25 newly diagnosed patients with primary or secondary (lymphoma-associated) AIC to test the hypothesis if these patients present a disease-specific immunological signature that could reveal pathophysiological clues and eventually be exploited as blood-based biomarker. Global TRB and IGH repertoire metrics as well as VJ gene usage distribution showed uniform characteristics for all lymphoma patients (high clonality and preferential usage of specific TRBV- and TRBJ genes), but no AIC-specific signature. Since T cell immune reactions toward antigens are unique and polyclonal, we clustered TCRß clones in-silico based on target recognition using the GLIPH (grouping of lymphocyte interactions by paratope hotspots) algorithm. This analysis revealed a considerable lack of physiological T cell clusters in patients with primary AIC. Interestingly, this signature did not discriminate between the different subentities of AIC and was also found in an independent cohort of 23 patients with active autoimmune hepatitis. Taken together, our data suggests that the identified T cell cluster signature could represent a blood biomarker of autoimmune conditions in general and should be functionally validated in future studies.

PTEN Loss in E-Cadherin-Deficient Mouse Mammary Epithelial Cells Rescues Apoptosis and Results in Development of Classical Invasive Lobular Carcinoma.

Invasive lobular carcinoma (ILC) is an aggressive breast cancer subtype with poor response to chemotherapy. Besides loss of E-cadherin, a hallmark of ILC, genetic inactivation of PTEN is frequently observed in patients. Through concomitant Cre-mediated inactivation of E-cadherin and PTEN in mammary epithelium, we generated a mouse model of classical ILC (CLC), the main histological ILC subtype. While loss of E-cadherin induced cell dissemination and apoptosis, additional PTEN inactivation promoted cell survival and rapid formation of invasive mammary tumors that recapitulate the histological and molecular features, estrogen receptor (ER) status, growth kinetics, metastatic behavior, and tumor microenvironment of human CLC. Combined inactivation of E-cadherin and PTEN is sufficient to cause CLC development. These CLCs showed significant tumor regression upon BEZ235-mediated inhibition of PI3K signaling. In summary, this mouse model provides important insights into CLC development and suggests inhibition of phosphatidylinositol 3-kinase (PI3K) signaling as a potential therapeutic strategy for targeting CLC.

An experimentally validated network of nine haematopoietic transcription factors reveals mechanisms of cell state stability.

Transcription factor (TF) networks determine cell-type identity by establishing and maintaining lineage-specific expression profiles, yet reconstruction of mammalian regulatory network models has been hampered by a lack of comprehensive functional validation of regulatory interactions. Here, we report comprehensive ChIP-Seq, transgenic and reporter gene experimental data that have allowed us to construct an experimentally validated regulatory network model for haematopoietic stem/progenitor cells (HSPCs). Model simulation coupled with subsequent experimental validation using single cell expression profiling revealed potential mechanisms for cell state stabilisation, and also how a leukaemogenic TF fusion protein perturbs key HSPC regulators. The approach presented here should help to improve our understanding of both normal physiological and disease processes.