Search

Abstract Immunotherapy is now an essential tool for cancer treatment, and the unique features of an individual’s T cell receptor repertoire are known to play a key role in its effectiveness. The repertoire, famously vast due to a cascade of cellular mechanisms, can be quantified using repertoire sequencing. In this study, we sampled the repertoire over several time points following treatment with anti-CTLA-4, in a syngeniec mouse model for colorectal cancer, generating a longitudinal dataset of T cell clones and their abundance. The dynamics of the repertoire can be observed in response to treatment over a period of four weeks, as clonal expansion of specific clones ascends and descends. The data made available here can be used to determine treatment and predict its effect, while also providing a unique look at the behavior of the immune system over time.

Abstract Background In the global effort to discover biomarkers for cancer prognosis, prediction tools have become essential resources. TCR (T cell receptor) repertoires contain important features that differentiate healthy controls from cancer patients or differentiate outcomes for patients being treated with different drugs. Considering, tools that can easily and quickly generate and identify important features out of TCR repertoire data and build accurate classifiers to predict future outcomes are essential. Results This paper introduces GENTLE (GENerator of T cell receptor repertoire features for machine LEarning): an open-source, user-friendly web-application tool that allows TCR repertoire researchers to discover important features; to create classifier models and evaluate them with metrics; and to quickly generate visualizations for data interpretations. We performed a case study with repertoires of TRegs (regulatory T cells) and TConvs (conventional T cells) from healthy controls versus patients with breast cancer. We showed that diversity features were able to distinguish between the groups. Moreover, the classifiers built with these features could correctly classify samples (‘Healthy’ or ‘Breast Cancer’)from the TRegs repertoire when trained with the TConvs repertoire, and from the TConvs repertoire when trained with the TRegs repertoire. Conclusion The paper walks through installing and using GENTLE and presents a case study and results to demonstrate the application’s utility. GENTLE is geared towards any researcher working with TCR repertoire data and aims to discover predictive features from these data and build accurate classifiers. GENTLE is available on https://github.com/dhiego22/gentle and https://share.streamlit.io/dhiego22/gentle/main/gentle.py .

Abstract Upon the development of a therapeutic, a successful response to a global pandemic relies on efficient worldwide distribution, a process constrained by our global shipping network. Most existing strategies seek to maximize the outflow of the therapeutics, hence optimizing for rapid dissemination. Here we find that this intuitive approach is, in fact, counterproductive. The reason is that by focusing strictly on the quantity of disseminated therapeutics, these strategies disregard the way in which this quantity distributes across destinations. Most crucially—they overlook the interplay of the therapeutic spreading patterns with those of the pathogens. This results in a discrepancy between supply and demand, that prohibits efficient mitigation even under optimal conditions of superfluous flow. To solve this, we design a dissemination strategy that naturally follows the predicted spreading patterns of the pathogens, optimizing not just for supply volume, but also for its congruency with the anticipated demand. Specifically, we show that epidemics spread relatively uniformly across all destinations, prompting us to introduce an equality constraint into our dissemination that prioritizes supply homogeneity. This strategy may, at times, slow down the supply rate in certain locations, however, thanks to its egalitarian nature, which mimics the flow of the pathogens, it provides a dramatic leap in overall mitigation efficiency, potentially saving more lives with orders of magnitude less resources.

Predicting an individual's response to therapy is an important goal for precision medicine. Here, the authors use an algorithm that takes into account the interaction type and directionality of signalling pathways in protein–protein interactions and find that their pathway analysis can predict essential genes, which may be a target for therapy.

Abstract Materials and construction methods of nests vary between bird species and at present, very little is known about the relationships between architecture and function in these structures. This study combines computational and experimental techniques to study the structural biology of nests fabricated by the edible nest swiftlet Aerodramus fuciphagus on vertical rock walls using threaded saliva. Utilizing its own saliva as a construction material allows the swiftlets full control over the structural features at a very high resolution in a process similar to additive manufacturing. It was hypothesized that the mechanical properties would vary between the structural regions of the nest (i.e. anchoring to the wall, center of the cup, and rim) mainly by means of architecture to offer structural support and bear the natural loads of birds and eggs. We generated numerical models of swiftlet nests from μCT scans based on collected swiftlet nests, which we loaded with a force of birds and eggs. This was done in order to study and assess the stress distribution that characterizes the specific nest’s architecture, evaluate its strength and weak points if any, as well as to understand the rationale and benefits that underlie this natural structure. We show that macro- and micro-scale structural patterns are identical in all nests, suggesting that their construction is governed by specific design principles. The nests’ response to applied loads of birds and eggs in finite element simulations suggests a mechanical overdesign strategy, which ensures the stresses experienced by its components in any loading scenario are actively minimized to be significantly smaller than the tensile fracture strength of the nests’ material. These findings highlight mechanical overdesign as a biological strategy for resilient, single-material constructions designed to protect eggs and hatchlings.

Abstract Audencel is a dendritic cell (DC)-based cellular cancer immunotherapy against glioblastoma multiforme (GBM). It is characterized by loading of DCs with autologous whole tumor lysate and in vitro maturation via “danger signals”. The recent phase II “GBM-Vax” trial showed no clinical efficacy for Audencel as assessed with progression-free and overall survival in all patients. Here we present immunological research accompanying the trial with a focus on immune system factors related to outcome and Audencel’s effect on the immune system. Methodologically, peripheral blood samples (from apheresis before Audencel or venipuncture during Audencel) were subjected to functional characterization via enzyme-linked immunospot (ELISPOT) assays connected with cytokine bead assays (CBAs) as well as phenotypical characterization via flow cytometry and mRNA quantification. GBM tissue samples (from surgery) were subjected to T cell receptor sequencing and immunohistochemistry. As results we found: Patients with favorable pre-existing anti-tumor characteristics lived longer under Audencel than Audencel patients without them. Pre-vaccination blood CD8+ T cell count and ELISPOT Granzyme B production capacity in vitro upon tumor antigen exposure were significantly correlated with overall survival. Despite Audencel’s general failure to induce a significant clinical response, it nevertheless seemed to have an effect on the immune system. For instance, Audencel led to a significant up-regulation of the Th1-related immunovariables ELISPOT IFNγ, the transcription factor T-bet in the blood and ELISPOT IL-2 in a dose-dependent manner upon vaccination. Post-vaccination levels of ELISPOT IFNγ and CD8+ cells in the blood were indicative of a significantly better survival. In summary, Audencel failed to reach an improvement of survival in the recent phase II clinical trial. No clinical efficacy was registered. Our concomitant immunological work presented here indicates that outcome under Audencel was influenced by the state of the immune system. On the other hand, Audencel also seemed to have stimulated the immune system. Overall, these immunological considerations suggest that DC immunotherapy against glioblastoma should be studied further – with the goal of translating an apparent immunological response into a clinical response. Future research should concentrate on investigating augmentation of immune reactions through combination therapies or on developing meaningful biomarkers.

It is increasingly clear that an extraordinarily diverse range of clinically important conditions—including infections, vaccinations, autoimmune diseases, transplants, transfusion reactions, aging, and cancers—leave telltale signatures in the millions of V(D)J-rearranged antibody and T cell receptor [TR per the Human Genome Organization (HUGO) nomenclature but more commonly known as TCR] genes collectively expressed by a person’s B cells (antibodies) and T cells. We refer to these as the immunome. Because of its diversity and complexity, the immunome provides singular opportunities for advancing personalized medicine by serving as the substrate for a highly multiplexed, near-universal blood test. Here we discuss some of these opportunities, the current state of immunome-based diagnostics, and highlight some of the challenges involved. We conclude with a call to clinicians, researchers, and others to join efforts with the Adaptive Immune Receptor Repertoire Community (AIRR-C) to realize the diagnostic potential of the immunome.

Summary: Biology of the response to anti-CTLA-4 involves the dynamics of specific T cell clones. Reasons for clinical success and failure of this treatment are still largely unknown. Here, we quantified the dynamics of the T cell receptor (TCR) repertoire, throughout 4 weeks involving treatment with anti-CTLA-4, in a syngeneic mouse model for colorectal cancer. These dynamics show an initial increase in clonality in tandem with a decrease in diversity, effects which gradually subside. Furthermore, response to treatment is tightly connected to the shared and public parts of the T cell repertoire. We were able to recognize time-dependent behaviors of specific TCR sequences and cell types and to show the response is dominated by specific motifs. We see that a single, specific time point might be useful to inform a physician of the true response to treatmentThe research further highlights the importance of temporal analyses of the immune response.

The partial success of tumor immunotherapy induced by checkpoint blockade, which is not antigen-specific, suggests that the immune system of some patients contain antigen receptors able to specifically identify tumor cells. Here we focused on T-cell receptor (TCR) repertoires associated with spontaneous breast cancer. We studied the alpha and beta chain CDR3 domains of TCR repertoires of CD4 T cells using deep sequencing of cell populations in mice and applied the results to published TCR sequence data obtained from human patients. We screened peripheral blood T cells obtained monthly from individual mice spontaneously developing breast tumors by 5 months. We then looked at identical TCR sequences in published human studies; we used TCGA data from tumors and healthy tissues of 1,256 breast cancer resections and from 4 focused studies including sequences from tumors, lymph nodes, blood and healthy tissues, and from single cell dataset of 3 breast cancer subjects. We now report that mice spontaneously developing breast cancer manifest shared, Public CDR3 regions in both their alpha and beta and that a significant number of women with early breast cancer manifest identical CDR3 sequences. These findings suggest that the development of breast cancer is associated, across species, with biomarker, exclusive TCR repertoires.

Most breast cancer patients are estrogen receptor positive and thus benefit from treatments that inhibit estrogen production; however, one third of tamoxifen-treated patients develops resistance and relapse. Here the authors show that tamoxifen resistant cells are resistant to chemotherapy because of BARD1 and BRCA1 upregulation.

Here, we outline an overview of the mammalian immune system that updates and extends the classical clonal selection paradigm. Rather than focusing on strict self-not-self discrimination, we propose that the system orchestrates variable inflammatory responses that maintain the body and its symbiosis with the microbiome while eliminating the threat from pathogenic infectious agents and from tumors. The paper makes four points: The immune system classifies healthy and pathologic states of the body—including both self and foreign elements—by deploying individual lymphocytes as cellular computing machines; immune cells transform input signals from the body into an output of specific immune reactions.Rather than independent clonal responses, groups of individually activated immune-system cells co-react in lymphoid organs to make collective decisions through a type of self-organizing swarm intelligence or crowd wisdom.Collective choices by swarms of immune cells, like those of schools of fish, are modified by relatively small numbers of individual regulators responding to shifting conditions—such collective inflammatory responses are dynamically responsive.Self-reactive autoantibody and T-cell receptor (TCR) repertoires shared by healthy individuals function in a biological version of experience-based supervised machine learning. Immune system decisions are primed by formative experience with training sets of self-antigens encountered during lymphocyte development; these initially trained T cell and B cell repertoires form a Wellness Profile that then guides immune responses to test sets of antigens encountered later. This experience-based machine learning strategy is analogous to that deployed by supervised machine-learning algorithms.We propose experiments to test these ideas. This overview of the immune system bears clinical implications for monitoring wellness and for treating autoimmune disease, cancer, and allograft reactions.

The T cell repertoire potentially presents complexity compatible, or greater than, that of the human brain. T cell based immune response is involved with practically every part of human physiology, and high-throughput biology needed to follow the T-cell repertoire has made great leaps with the advent of massive parallel sequencing [1]. Nevertheless, tools to handle and observe the dynamics of this complexity have only recently started to emerge [e.g., 2, 3, 4] in parallel with sequencing technologies. Here, we present a network-based view of the dynamics of the T cell repertoire, during the course of mammary tumors development in a mouse model. The transition from the T cell receptor as a feature, to network-based clustering, followed by network-based temporal analyses, provides novel insights to the workings of the system and provides novel tools to observe cancer progression via the perspective of the immune system. The crux of the approach here is at the network-motivated clustering. The purpose of the clustering step is not merely data reduction and exposing structures, but rather to detect hubs, or attractors, within the T cell receptor repertoire that might shed light on the behavior of the immune system as a dynamic network. The Clone-Attractor is in fact an extension of the clone concept, i.e., instead of looking at particular clones we observe the extended clonal network by assigning clusters to graph nodes and edges to adjacent clusters (editing distance metric). Viewing the system as dynamical brings to the fore the notion of an attractors landscape, hence the possibility to chart this space and map the sample state at a given time to a vector in this large space. Based on this representation we applied two different methods to demonstrate its effectiveness in identifying changes in the repertoire that correlate with changes in the phenotype: (1) network analysis of the TCR repertoire in which two measures were calculated and demonstrated the ability to differentiate control from transgenic samples, and, (2) machine learning classifier capable of both stratifying control and trangenic samples, as well as to stratify pre-cancer and cancer samples.

The T cell receptor (TCR) controls the cellular adaptive immune response to antigens, but our understanding of TCR repertoire diversity and response to challenge is still incomplete. For example, TCR clones shared by different individuals with minimal alteration to germline gene sequences (public clones) are detectable in all vertebrates, but their significance is unknown. Although small in size, the zebrafish TCR repertoire is controlled by processes similar to those operating in mammals. Thus, we studied the zebrafish TCR repertoire and its response to stimulation with self and foreign antigens. We found that cross-reactive public TCRs dominate the T cell response, endowing a limited TCR repertoire with the ability to cope with diverse antigenic challenges. These features of vertebrate public TCRs might provide a mechanism for the rapid generation of protective T cell immunity, allowing a short temporal window for the development of more specific private T cell responses.

Mice overexpressing the longevity protein SIRT6 or deficient for the liver’s most prevalent microRNA miR-122 display a similar set of phenotypes, including improved lipid profile and protection against damage linked to obesity. Here, we show that miR-122 and SIRT6 negatively regulate each other’s expression. SIRT6 downregulates miR-122 by deacetylating H3K56 in the promoter region. MiR-122 binds to three sites on the SIRT6 3′ UTR and reduces its levels. The interplay between SIRT6 and miR-122 is manifested in two physiologically relevant ways in the liver. First, they oppositely regulate a similar set of metabolic genes and fatty acid β-oxidation. Second, in hepatocellular carcinoma patients, loss of a negative correlation between SIRT6 and miR-122 expression is significantly associated with better prognosis. These findings show that SIRT6 and miR-122 negatively regulate each other to control various aspects of liver physiology and SIRT6-miR-122 correlation may serve as a biomarker for hepatocarcinoma prognosis.

The field of computer modeling and simulation of biological systems is rapidly advancing, backed by significant progress in the fields of experimentation techniques, computer hardware, and programming software. The result of a simulation may be delivered in several ways, from numerical results, through graphs of the simulated run, to a visualization of the simulation. The vision of an in-silico experiment mimicking an in-vitro or in-vivo experiment as it is viewed under a microscope is appealing but technically demanding and computationally intensive. Here, we report “Cell Studio,” a generic, hybrid platform to simulate an immune microenvironment with biological and biophysical rules. We use game engines—generic programs for game creation which offer ready-made assets and tools—to create a visualized, interactive 3D simulation. We also utilize a scalable architecture that delegates the computational load to a server. The user may view the simulation, move the “camera” around, stop, fast-forward, and rewind it and inject soluble molecules into the extracellular medium at any point in time. During simulation, graphs are created in real time for a broad view of system-wide processes. The model is parametrized using a user-friendly Graphical User Interface (GUI). We show a simple validation simulation and compare its results with those from a “classical” simulation, validated against a “wet” experiment. We believe that interactive, real-time 3D visualization may aid in generating insights from the model and encourage intuition about the immunological scenario.

Immune tolerance toward “self” is critical in multiple immune disorders. While there are several mechanisms to describe the involvement of immune cells in the process, the role of peripheral tissue cells in that context is not yet clear. The theory of ecoimmunity postulates that interactions between immune and tissue cells represent a predator–prey relationship. A lifelong interaction, shaped mainly during early ontogeny, leads to selection of nonimmune cell phenotypes. Normally, therefore, nonimmune cells that evolve alongside an intact immune system would be phenotypically capable of evading immune responses, and cells whose phenotype falls short of satisfying this steady state would expire under hostile immune responses. This view was supported until recently by experimental evidence showing an inferior endurance of severe combined immunodeficiency (SCID)-derived pancreatic islets when engrafted into syngeneic immune-intact wild-type (WT) mice, relative to islets from WT. Here we extend the experimental exploration of ecoimmunity by searching for the presence of the phenotypic changes suggested by the theory. Immune-related phenotypes of islets, spleen, and bone marrow immune cells were determined, as well as SCID and WT nonlymphocytic cells. Islet submass grafting was performed to depict syngeneic graft functionality. Islet cultures were examined under both resting and inflamed conditions for expression of CD40 and major histocompatibility complex (MHC) class I/II and release of interleukin-1α (IL-1α), IL-1β, IL-6, tumor necrosis factor-α (TNF-α), IL-10, and insulin. Results depict multiple pathways that appear to be related to the sculpting of nonimmune cells by immune cells; 59 SCID islet genes displayed relative expression changes compared with WT islets. SCID cells expressed lower tolerability to inflammation and higher levels of immune-related molecules, including MHC class I. Accordingly, islets exhibited a marked increase in insulin release upon immunocyte depletion, in effect resuming endocrine function that was otherwise suppressed by resident immunocytes. This work provides further support of the ecoimmunity theory and encourages subsequent studies to identify its role in the emergence and treatment of autoimmune pathologies, transplant rejection, and cancer.

Fragile X syndrome (FXS) is the most frequent inherited form of mental retardation. The cause for this X-linked disorder is the silencing of the fragile X mental retardation 1 (fmr1) gene and the absence of the fragile X mental retardation protein (Fmrp). The RNA-binding protein Fmrp represses protein translation, particularly in synapses. In Drosophila, Fmrp interacts with the adenosine deaminase acting on RNA (Adar) enzymes. Adar enzymes convert adenosine to inosine (A-to-I) and modify the sequence of RNA transcripts. Utilizing the fmr1 zebrafish mutant (fmr1-/-), we studied Fmrp-dependent neuronal circuit formation, behavior, and Adar-mediated RNA editing. By combining behavior analyses and live imaging of single axons and synapses, we showed hyperlocomotor activity, as well as increased axonal branching and synaptic density, in fmr1-/- larvae. We identified thousands of clustered RNA editing sites in the zebrafish transcriptome and showed that Fmrp biochemically interacts with the Adar2a protein. The expression levels of the adar genes and Adar2 protein increased in fmr1-/- zebrafish. Microfluidic-based multiplex PCR coupled with deep sequencing showed a mild increase in A-to-I RNA editing levels in evolutionarily conserved neuronal and synaptic Adar-targets in fmr1-/- larvae. These findings suggest that loss of Fmrp results in increased Adar-mediated RNA editing activity on target-specific RNAs, which, in turn, might alter neuronal circuit formation and behavior in FXS.

In the following pages we propose a theory on cognitive systems and the com-mon strategies of perception, which are at the basis of their function. We demon-strate that these strategies are easily seen to be in place in known cognitive sys-tems such as vision and language. Furthermore we show that taking these strat-egies into consideration implies a new outlook on immune function calling for a new appraisal of the immune system as a cognitive system.

The transcriptional networks that regulate gene expression and modifications to this network are at the core of the cancer phenotype. MicroRNAs, a well-studied species of small non-coding RNA molecules, have been shown to have a central role in regulating gene expression as part of this transcriptional network. Further, microRNA deregulation is associated with cancer development and with tumor progression. Glioblastoma Multiform (GBM) is the most common, aggressive and malignant primary tumor of the brain and is associated with one of the worst 5-year survival rates among all human cancers. To study the transcriptional network and its modifications in GBM, we utilized gene expression, microRNA sequencing, whole genome sequencing and clinical data from hundreds of patients from different datasets. Using these data and a novel microRNA-gene association approach we introduce, we have identified unique microRNAs and their associated genes. This unique behavior is composed of the ability of the quantifiable association of the microRNA and the gene expression levels, which we show stratify patients into clinical subgroups of high statistical significance. Importantly, this stratification goes unobserved by other methods and is not affiliated by other subsets or phenotypes within the data. To investigate the robustness of the introduced approach, we demonstrate, in unrelated datasets, robustness of findings. Among the set of identified microRNA-gene associations, we closely study the example of MAF and hsa-miR-330-3p, and show how their co-behavior stratifies patients into prognosis clinical groups and how whole genome sequences tells us more about a specific genomic variation as a possible basis for patient variances. We argue that these identified associations may indicate previously unexplored specific disease control mechanisms and may be used as basis for further study and for possible therapeutic intervention.

High resolution, system-wide characterizations have demonstrated the capacity to identify genomic regions that undergo genomic aberrations. Such research efforts often aim at associating these regions with disease etiology and outcome. Identifying the corresponding biologic processes that are responsible for disease and its outcome remains challenging. Using novel analytic methods that utilize the structure of biologic networks, we are able to identify the specific networks that are highly significantly, nonrandomly altered by regions of copy number amplification observed in a systems-wide analysis. We demonstrate this method in breast cancer, where the state of a subset of the pathways identified through these regions is shown to be highly associated with disease survival and recurrence.

We introduce a novel technique to determine the expression state of a gene from quantitative information measuring its expression. Adopting a productive abstraction from current thinking in molecular biology, we consider two expression states for a gene--Up or Down. We determine this state by using a statistical model that assumes the data behaves as a combination of two biological distributions. Given a cohort of hybridizations, our algorithm predicts, for the single reading, the probability of each gene's being in an Up or a Down state in each hybridization. Using a series of publicly available gene expression data sets, we demonstrate that our algorithm outperforms the prevalent algorithm. We also show that our algorithm can be used in conjunction with expression adjustment techniques to produce a more biologically sound gene-state call. The technique we present here enables a routine update, where the continuously evolving expression level adjustments feed into gene-state calculations. The technique can be applied in almost any multi-sample gene expression experiment, and holds equal promise for protein abundance experiments.

Cancer is recognized to be a family of gene-based diseases whose causes are to be found in disruptions of basic biologic processes. An increasingly deep catalogue of canonical networks details the specific molecular interaction of genes and their products. However, mapping of disease phenotypes to alterations of these networks of interactions is accomplished indirectly and non-systematically. Here we objectively identify pathways associated with malignancy, staging, and outcome in cancer through application of an analytic approach that systematically evaluates differences in the activity and consistency of interactions within canonical biologic processes. Using large collections of publicly accessible genome-wide gene expression, we identify small, common sets of pathways - Trka Receptor, Apoptosis response to DNA Damage, Ceramide, Telomerase, CD40L and Calcineurin - whose differences robustly distinguish diverse tumor types from corresponding normal samples, predict tumor grade, and distinguish phenotypes such as estrogen receptor status and p53 mutation state. Pathways identified through this analysis perform as well or better than phenotypes used in the original studies in predicting cancer outcome. This approach provides a means to use genome-wide characterizations to map key biological processes to important clinical features in disease.

Experiments have generated a plethora of data about the genes, molecules, and cells involved in thymocyte development. Here, we use a computer-driven simulation that uses data about thymocyte development to generate an integrated dynamic representation-a novel technology we have termed reactive animation (RA). RA reveals emergent properties in complex dynamic biological systems. We apply RA to thymocyte development by reproducing and extending the effects of known gene knockouts: CXCR4 and CCR9. RA simulation revealed a previously unidentified role of thymocyte competition for major histocompatability complex presentation. We now report that such competition is required for normal anatomical compartmentalization, can influence the rate of thymocyte velocities within chemokine gradients, and can account for the disproportion between single-positive CD4 and CD8 lineages developing from double-positive precursors.