Nonlinear analysis of biological sequences

This is the final report of a three-year, Laboratory Directed Research and Development (LDRD) project at the Los Alamos National Laboratory (LANL). The main objectives of this project involved deriving new capabilities for analyzing biological sequences. The authors focused on tabulating the statistical properties exhibited by Human coding DNA sequences and on techniques of inferring the phylogenetic relationships among protein sequences related by descent

The influence of T cell development on pathogen specificity and autoreactivity

T cells orchestrate adaptive immune responses upon activation. T cell activation requires sufficiently strong binding of T cell receptors on their surface to short peptides derived from foreign proteins bound to protein products of the major histocompatibility (MHC) gene products, which are displayed on the surface of antigen presenting cells. T cells can also interact with peptide-MHC complexes, where the peptide is derived from host (self) proteins. A diverse repertoire of relatively self-tolerant T cell receptors is selected in the thymus. We study a model, computationally and analytically, to describe how thymic selection shapes the repertoire of T cell receptors, such that T cell receptor recognition of pathogenic peptides is both specific and degenerate. We also discuss the escape probability of autoimmune T cells from the thymus.Comment: 12 pages, 7 figure

Gene signature of the post-Chernobyl papillary thyroid cancer

Purpose: Following the nuclear accidents in Chernobyl and later in Fukushima, the nuclear community has been faced with important issues concerning how to search for and diagnose biological consequences of low-dose internal radiation contamination. Although after the Chernobyl accident an increase in childhood papillary thyroid cancer (PTC) was observed, it is still not clear whether the molecular biology of PTCs associated with low-dose radiation exposure differs from that of sporadic PTC. Methods: We investigated tissue samples from 65 children/young adults with PTC using DNA microarray (Affymetrix, Human Genome U133 2.0 Plus) with the aim of identifying molecular differences between radiation-induced (exposed to Chernobyl radiation, ECR) and sporadic PTC. All participants were resident in the same region so that confounding factors related to genetics or environment were minimized. Results: There were small but significant differences in the gene expression profiles between ECR and non-ECR PTC (global test, p < 0.01), with 300 differently expressed probe sets (p < 0.001) corresponding to 239 genes. Multifactorial analysis of variance showed that besides radiation exposure history, the BRAF mutation exhibited independent effects on the PTC expression profile; the histological subset and patient age at diagnosis had negligible effects. Ten genes (PPME1, HDAC11, SOCS7, CIC, THRA, ERBB2, PPP1R9A, HDGF, RAD51AP1, and CDK1) from the 19 investigated with quantitative RT-PCR were confirmed as being associated with radiation exposure in an independent, validation set of samples. Conclusion: Significant, but subtle, differences in gene expression in the post-Chernobyl PTC are associated with previous low-dose radiation exposure

Absence of a specific radiation signature in post-Chernobyl thyroid cancers

Thyroid cancers have been the main medical consequence of the Chernobyl accident. On the basis of their pathological features and of the fact that a large proportion of them demonstrate RET-PTC translocations, these cancers are considered as similar to classical sporadic papillary carcinomas, although molecular alterations differ between both tumours. We analysed gene expression in post-Chernobyl cancers, sporadic papillary carcinomas and compared to autonomous adenomas used as controls. Unsupervised clustering of these data did not distinguish between the cancers, but separates both cancers from adenomas. No gene signature separating sporadic from post-Chernobyl PTC (chPTC) could be found using supervised and unsupervised classification methods although such a signature is demonstrated for cancers and adenomas. Furthermore, we demonstrate that pooled RNA from sporadic and chPTC are as strongly correlated as two independent sporadic PTC pools, one from Europe, one from the US involving patients not exposed to Chernobyl radiations. This result relies on cDNA and Affymetrix microarrays. Thus, platform-specific artifacts are controlled for. Our findings suggest the absence of a radiation fingerprint in the chPTC and support the concept that post-Chernobyl cancer data, for which the cancer-causing event and its date are known, are a unique source of information to study naturally occurring papillary carcinomas

Genome-wide gene expression profiling suggests distinct radiation susceptibilities in sporadic and post-Chernobyl papillary thyroid cancers

Papillary thyroid cancers (PTCs) incidence dramatically increased in the vicinity of Chernobyl. The cancer-initiating role of radiation elsewhere is debated. Therefore, we searched for a signature distinguishing radio-induced from sporadic cancers. Using microarrays, we compared the expression profiles of PTCs from the Chernobyl Tissue Bank (CTB, n=12) and from French patients with no history of exposure to ionising radiations (n=14). We also compared the transcriptional responses of human lymphocytes to the presumed aetiological agents initiating these tumours, γ-radiation and H2O2. On a global scale, the transcriptomes of CTB and French tumours are indistinguishable, and the transcriptional responses to γ-radiation and H2O2 are similar. On a finer scale, a 118 genes signature discriminated the γ-radiation and H2O2 responses. This signature could be used to classify the tumours as CTB or French with an error of 15–27%. Similar results were obtained with an independent signature of 13 genes involved in homologous recombination. Although sporadic and radio-induced PTCs represent the same disease, they are distinguishable with molecular signatures reflecting specific responses to γ-radiation and H2O2. These signatures in PTCs could reflect the susceptibility profiles of the patients, suggesting the feasibility of a radiation susceptibility test

Balancing Robustness against the Dangers of Multiple Attractors in a Hopfield-Type Model of Biological Attractors

Background: Many chronic human diseases are of unclear origin, and persist long beyond any known insult or instigating factor. These diseases may represent a structurally normal biologic network that has become trapped within the basin of an abnormal attractor. Methodology/Principal Findings: We used the Hopfield net as the archetypical example of a dynamic biological network. By progressively removing the links of fully connected Hopfield nets, we found that a designated attractor of the nets could still be supported until only slightly more than 1 link per node remained. As the number of links approached this minimum value, the rate of convergence to this attractor from an arbitrary starting state increased dramatically. Furthermore, with more than about twice the minimum of links, the net became increasingly able to support a second attractor. Conclusions/Significance: We speculate that homeostatic biological networks may have evolved to assume a degree of connectivity that balances robustness and agility against the dangers of becoming trapped in an abnormal attractor

Meta-analysis of muscle transcriptome data using the MADMuscle database reveals biologically relevant gene patterns

<p>Abstract</p> <p>Background</p> <p>DNA microarray technology has had a great impact on muscle research and microarray gene expression data has been widely used to identify gene signatures characteristic of the studied conditions. With the rapid accumulation of muscle microarray data, it is of great interest to understand how to compare and combine data across multiple studies. Meta-analysis of transcriptome data is a valuable method to achieve it. It enables to highlight conserved gene signatures between multiple independent studies. However, using it is made difficult by the diversity of the available data: different microarray platforms, different gene nomenclature, different species studied, etc.</p> <p>Description</p> <p>We have developed a system tool dedicated to muscle transcriptome data. This system comprises a collection of microarray data as well as a query tool. This latter allows the user to extract similar clusters of co-expressed genes from the database, using an input gene list. Common and relevant gene signatures can thus be searched more easily. The dedicated database consists in a large compendium of public data (more than 500 data sets) related to muscle (skeletal and heart). These studies included seven different animal species from invertebrates (<it>Drosophila melanogaster, Caenorhabditis elegans</it>) and vertebrates (<it>Homo sapiens, Mus musculus, Rattus norvegicus, Canis familiaris, Gallus gallus</it>). After a renormalization step, clusters of co-expressed genes were identified in each dataset. The lists of co-expressed genes were annotated using a unified re-annotation procedure. These gene lists were compared to find significant overlaps between studies.</p> <p>Conclusions</p> <p>Applied to this large compendium of data sets, meta-analyses demonstrated that conserved patterns between species could be identified. Focusing on a specific pathology (Duchenne Muscular Dystrophy) we validated results across independent studies and revealed robust biomarkers and new pathways of interest. The meta-analyses performed with MADMuscle show the usefulness of this approach. Our method can be applied to all public transcriptome data.</p

Representation in the (Artificial) Immune System

Much of contemporary research in Artificial Immune Systems (AIS) has partitioned into either algorithmic machine learning and optimisation, or, modelling biologically plausible dynamical systems, with little overlap between. We propose that this dichotomy is somewhat to blame for the lack of significant advancement of the field in either direction and demonstrate how a simplistic interpretation of Perelson’s shape-space formalism may have largely contributed to this dichotomy. In this paper, we motivate and derive an alternative representational abstraction. To do so we consider the validity of shape-space from both the biological and machine learning perspectives. We then take steps towards formally integrating these perspectives into a coherent computational model of notions such as life-long learning, degeneracy, constructive representations and contextual recognition—rhetoric that has long inspired work in AIS, while remaining largely devoid of operational definition

Prediction of proteasome cleavage motifs by neural networks

We present a predictive method that can simulate an essential step in the antigen presentation in higher vertebrates, namely the step involving the proteasomal degradation of polypeptides into fragments which have the potential to bind to MHC Class I molecules. Proteasomal cleavage prediction algorithms published so far were trained on data from in vitro digestion experiments with constitutive proteasomes. As a result, they did not take into account the characteristics of the structurally modified proteasomes--often called immunoproteasomes--found in cells stimulated by gamma-interferon under physiological conditions. Our algorithm has been trained not only on in vitro data, but also on MHC Class I ligand data, which reflect a combination of immunoproteasome and constitutive proteasome specificity. This feature, together with the use of neural networks, a non-linear classification technique, make the prediction of MHC Class I ligand boundaries more accurate: 65% of the cleavage sites and 85% of the non-cleavage sites are correctly determined. Moreover, we show that the neural networks trained on the constitutive proteasome data learns a specificity that differs from that of the networks trained on MHC Class I ligands, i.e. the specificity of the immunoproteasome is different than the constitutive proteasome. The tools developed in this study in combination with a predictor of MHC and TAP binding capacity should give a more complete prediction of the generation and presentation of peptides on MHC Class I molecules. Here we demonstrate that such an approach produces an accurate prediction of the CTL the epitopes in HIV Nef. The method is available at www.cbs.dtu.dk/services/NetChop/.Journal Articleinfo:eu-repo/semantics/publishe