James A. Eddy

IL-10 Alters Immunoproteostasis in APP Mice, Increasing Plaque Burden and Worsening Cognitive Behavior

Anti-inflammatory strategies are proposed to have beneficial effects in Alzheimers disease. To explore how anti-inflammatory cytokine signaling affects Aβ pathology, we investigated the effects of adeno-associated virus (AAV2/1)-mediated expression of Interleukin (IL)-10 in the brains of APP transgenic mouse models. IL-10 expression resulted in increased Aβ accumulation and impaired memory in APP mice. A focused transcriptome analysis revealed changes consistent with enhanced IL-10 signaling and increased ApoE expression in IL-10-expressing APP mice. ApoE protein was selectively increased in the plaque-associated insoluble cellular fraction, likely because of direct interaction with aggregated Aβ in the IL-10-expressing APP mice. Ex vivo studies also show that IL-10 and ApoE can individually impair glial Aβ phagocytosis. Our observations that IL-10 has an unexpected negative effect on Aβ proteostasis and cognition in APP mouse models demonstrate the complex interplay between innate immunity and proteostasis in neurodegenerative diseases, an interaction we call immunoproteostasis.

Reconstruction of genome-scale metabolic models for 126 human tissues using mCADRE

BackgroundHuman tissues perform diverse metabolic functions. Mapping out these tissue-specific functions in genome-scale models will advance our understanding of the metabolic basis of various physiological and pathological processes. The global knowledgebase of metabolic functions categorized for the human genome (Human Recon 1) coupled with abundant high-throughput data now makes possible the reconstruction of tissue-specific metabolic models. However, the number of available tissue-specific models remains incomplete compared with the large diversity of human tissues.ResultsWe developed a method called metabolic Context-specificity Assessed by Deterministic Reaction Evaluation (mCADRE). mCADRE is able to infer a tissue-specific network based on gene expression data and metabolic network topology, along with evaluation of functional capabilities during model building. mCADRE produces models with similar or better functionality and achieves dramatic computational speed up over existing methods. Using our method, we reconstructed draft genome-scale metabolic models for 126 human tissue and cell types. Among these, there are models for 26 tumor tissues along with their normal counterparts, and 30 different brain tissues. We performed pathway-level analyses of this large collection of tissue-specific models and identified the eicosanoid metabolic pathway, especially reactions catalyzing the production of leukotrienes from arachidnoic acid, as potential drug targets that selectively affect tumor tissues.ConclusionsThis large collection of 126 genome-scale draft metabolic models provides a useful resource for studying the metabolic basis for a variety of human diseases across many tissues. The functionality of the resulting models and the fast computational speed of the mCADRE algorithm make it a useful tool to build and update tissue-specific metabolic models.

Behavior-specific changes in transcriptional modules lead to distinct and predictable neurogenomic states

Using brain transcriptomic profiles from 853 individual honey bees exhibiting 48 distinct behavioral phenotypes in naturalistic contexts, we report that behavior-specific neurogenomic states can be inferred from the coordinated action of transcription factors (TFs) and their predicted target genes. Unsupervised hierarchical clustering of these transcriptomic profiles showed three clusters that correspond to three ecologically important behavioral categories: aggression, maturation, and foraging. To explore the genetic influences potentially regulating these behavior-specific neurogenomic states, we reconstructed a brain transcriptional regulatory network (TRN) model. This brain TRN quantitatively predicts with high accuracy gene expression changes of more than 2,000 genes involved in behavior, even for behavioral phenotypes on which it was not trained, suggesting that there is a core set of TFs that regulates behavior-specific gene expression in the bee brain, and other TFs more specific to particular categories. TFs playing key roles in the TRN include well-known regulators of neural and behavioral plasticity, e.g., Creb, as well as TFs better known in other biological contexts, e.g., NF-κB (immunity). Our results reveal three insights concerning the relationship between genes and behavior. First, distinct behaviors are subserved by distinct neurogenomic states in the brain. Second, the neurogenomic states underlying different behaviors rely upon both shared and distinct transcriptional modules. Third, despite the complexity of the brain, simple linear relationships between TFs and their putative target genes are a surprisingly prominent feature of the networks underlying behavior.

Metabolic network reconstruction and genome-scale model of butanol-producing strain Clostridium beijerinckii NCIMB 8052

BackgroundSolventogenic clostridia offer a sustainable alternative to petroleum-based production of butanol--an important chemical feedstock and potential fuel additive or replacement. C. beijerinckii is an attractive microorganism for strain design to improve butanol production because it (i) naturally produces the highest recorded butanol concentrations as a byproduct of fermentation; and (ii) can co-ferment pentose and hexose sugars (the primary products from lignocellulosic hydrolysis). Interrogating C. beijerinckii metabolism from a systems viewpoint using constraint-based modeling allows for simulation of the global effect of genetic modifications.ResultsWe present the first genome-scale metabolic model (i CM925) for C. beijerinckii, containing 925 genes, 938 reactions, and 881 metabolites. To build the model we employed a semi-automated procedure that integrated genome annotation information from KEGG, BioCyc, and The SEED, and utilized computational algorithms with manual curation to improve model completeness. Interestingly, we found only a 34% overlap in reactions collected from the three databases--highlighting the importance of evaluating the predictive accuracy of the resulting genome-scale model. To validate i CM925, we conducted fermentation experiments using the NCIMB 8052 strain, and evaluated the ability of the model to simulate measured substrate uptake and product production rates. Experimentally observed fermentation profiles were found to lie within the solution space of the model; however, under an optimal growth objective, additional constraints were needed to reproduce the observed profiles--suggesting the existence of selective pressures other than optimal growth. Notably, a significantly enriched fraction of actively utilized reactions in simulations--constrained to reflect experimental rates--originated from the set of reactions that overlapped between all three databases (P = 3.52 × 10-9, Fishers exact test). Inhibition of the hydrogenase reaction was found to have a strong effect on butanol formation--as experimentally observed.ConclusionsMicrobial production of butanol by C. beijerinckii offers a promising, sustainable, method for generation of this important chemical and potential biofuel. i CM925 is a predictive model that can accurately reproduce physiological behavior and provide insight into the underlying mechanisms of microbial butanol production. As such, the model will be instrumental in efforts to better understand, and metabolically engineer, this microorganism for improved butanol production.

Accomplishments in genome‐scale in silico modeling for industrial and medical biotechnology

Driven by advancements in high‐throughput biological technologies and the growing number of sequenced genomes, the construction of in silico models at the genome scale has provided powerful tools to investigate a vast array of biological systems and applications. Here, we review comprehensively the uses of such models in industrial and medical biotechnology, including biofuel generation, food production, and drug development. While the use of in silico models is still in its early stages for delivering to industry, significant initial successes have been achieved. For the cases presented here, genome‐scale models predict engineering strategies to enhance properties of interest in an organism or to inhibit harmful mechanisms of pathogens. Going forward, genome‐scale in silico models promise to extend their application and analysis scope to become a transformative tool in biotechnology.

In silico models of cancer

Cancer is a complex disease that involves multiple types of biological interactions across diverse physical, temporal, and biological scales. This complexity presents substantial challenges for the characterization of cancer biology, and motivates the study of cancer in the context of molecular, cellular, and physiological systems. Computational models of cancer are being developed to aid both biological discovery and clinical medicine. The development of these in silico models is facilitated by rapidly advancing experimental and analytical tools that generate information‐rich, high‐throughput biological data. Statistical models of cancer at the genomic, transcriptomic, and pathway levels have proven effective in developing diagnostic and prognostic molecular signatures, as well as in identifying perturbed pathways. Statistically inferred network models can prove useful in settings where data overfitting can be avoided, and provide an important means for biological discovery. Mechanistically based signaling and metabolic models that apply a priori knowledge of biochemical processes derived from experiments can also be reconstructed where data are available, and can provide insight and predictive ability regarding the behavior of these systems. At longer length scales, continuum and agent‐based models of the tumor microenvironment and other tissue‐level interactions enable modeling of cancer cell populations and tumor progression. Even though cancer has been among the most‐studied human diseases using systems approaches, significant challenges remain before the enormous potential of in silico cancer biology can be fully realized. Copyright

Identifying Tightly Regulated and Variably Expressed Networks by Differential Rank Conservation (DIRAC)

A powerful way to separate signal from noise in biology is to convert the molecular data from individual genes or proteins into an analysis of comparative biological network behaviors. One of the limitations of previous network analyses is that they do not take into account the combinatorial nature of gene interactions within the network. We report here a new technique, Differential Rank Conservation (DIRAC), which permits one to assess these combinatorial interactions to quantify various biological pathways or networks in a comparative sense, and to determine how they change in different individuals experiencing the same disease process. This approach is based on the relative expression values of participating genes—i.e., the ordering of expression within network profiles. DIRAC provides quantitative measures of how network rankings differ either among networks for a selected phenotype or among phenotypes for a selected network. We examined disease phenotypes including cancer subtypes and neurological disorders and identified networks that are tightly regulated, as defined by high conservation of transcript ordering. Interestingly, we observed a strong trend to looser network regulation in more malignant phenotypes and later stages of disease. At a sample level, DIRAC can detect a change in ranking between phenotypes for any selected network. Variably expressed networks represent statistically robust differences between disease states and serve as signatures for accurate molecular classification, validating the information about expression patterns captured by DIRAC. Importantly, DIRAC can be applied not only to transcriptomic data, but to any ordinal data type.

Cell type-specific genes show striking and distinct patterns of spatial expression in the mouse brain

To characterize gene expression patterns in the regional subdivisions of the mammalian brain, we integrated spatial gene expression patterns from the Allen Brain Atlas for the adult mouse with panels of cell type-specific genes for neurons, astrocytes, and oligodendrocytes from previously published transcriptome profiling experiments. We found that the combined spatial expression patterns of 170 neuron-specific transcripts revealed strikingly clear and symmetrical signatures for most of the brain’s major subdivisions. Moreover, the brain expression spatial signatures correspond to anatomical structures and may even reflect developmental ontogeny. Spatial expression profiles of astrocyte- and oligodendrocyte-specific genes also revealed regional differences; these defined fewer regions and were less distinct but still symmetrical in the coronal plane. Follow-up analysis suggested that region-based clustering of neuron-specific genes was related to (i) a combination of individual genes with restricted expression patterns, (ii) region-specific differences in the relative expression of functional groups of genes, and (iii) regional differences in neuronal density. Products from some of these neuron-specific genes are present in peripheral blood, raising the possibility that they could reflect the activities of disease- or injury-perturbed networks and collectively function as biomarkers for clinical disease diagnostics.

Relative expression analysis for molecular cancer diagnosis and prognosis.

The enormous amount of biomolecule measurement data generated from high-throughput technologies has brought an increased need for computational tools in biological analyses. Such tools can enhance our understanding of human health and genetic diseases, such as cancer, by accurately classifying phenotypes, detecting the presence of disease, discriminating among cancer sub-types, predicting clinical outcomes, and characterizing disease progression. In the case of gene expression microarray data, standard statistical learning methods have been used to identify classifiers that can accurately distinguish disease phenotypes. However, these mathematical prediction rules are often highly complex, and they lack the convenience and simplicity desired for extracting underlying biological meaning or transitioning into the clinic. In this review, we survey a powerful collection of computational methods for analyzing transcriptomic microarray data that address these limitations. Relative Expression Analysis (RXA) is based only on the relative orderings among the expressions of a small number of genes. Specifically, we provide a description of the first and simplest example of RXA, the k-TSP classifier, which is based on k pairs of genes; the case k = 1 is the TSP classifier. Given their simplicity and ease of biological interpretation, as well as their invariance to data normalization and parameter-fitting, these classifiers have been widely applied in aiding molecular diagnostics in a broad range of human cancers. We review several studies which demonstrate accurate classification of disease phenotypes (e.g., cancer vs. normal), cancer subclasses (e.g., AML vs. ALL, GIST vs. LMS), disease outcomes (e.g., metastasis, survival), and diverse human pathologies assayed through blood-borne leukocytes. The studies presented demonstrate that RXA—specifically the TSP and k-TSP classifiers—is a promising new class of computational methods for analyzing high-throughput data, and has the potential to significantly contribute to molecular cancer diagnosis and prognosis.

Integrated Proteomics and Genomics Analysis Reveals a Novel Mesenchymal to Epithelial Reverting Transition in Leiomyosarcoma through Regulation of Slug

Leiomyosarcoma is one of the most common mesenchymal tumors. Proteomics profiling analysis by reverse-phase protein lysate array surprisingly revealed that expression of the epithelial marker E-cadherin (encoded by CDH1) was significantly elevated in a subset of leiomyosarcomas. In contrast, E-cadherin was rarely expressed in the gastrointestinal stromal tumors, another major mesenchymal tumor type. We further sought to 1) validate this finding, 2) determine whether there is a mesenchymal to epithelial reverting transition (MErT) in leiomyosarcoma, and if so 3) elucidate the regulatory mechanism responsible for this MErT. Our data showed that the epithelial cell markers E-cadherin, epithelial membrane antigen, cytokeratin AE1/AE3, and pan-cytokeratin were often detected immunohistochemically in leiomyosarcoma tumor cells on tissue microarray. Interestingly, the E-cadherin protein expression was correlated with better survival in leiomyosarcoma patients. Whole genome microarray was used for transcriptomics analysis, and the epithelial gene expression signature was also associated with better survival. Bioinformatics analysis of transcriptome data showed an inverse correlation between E-cadherin and E-cadherin repressor Slug (SNAI2) expression in leiomyosarcoma, and this inverse correlation was validated on tissue microarray by immunohistochemical staining of E-cadherin and Slug. Knockdown of Slug expression in SK-LMS-1 leiomyosarcoma cells by siRNA significantly increased E-cadherin; decreased the mesenchymal markers vimentin and N-cadherin (encoded by CDH2); and significantly decreased cell proliferation, invasion, and migration. An increase in Slug expression by pCMV6-XL5-Slug transfection decreased E-cadherin and increased vimentin and N-cadherin. Thus, MErT, which is mediated through regulation of Slug, is a clinically significant phenotype in leiomyosarcoma.