InnateDB: facilitating systems-level analyses of the mammalian innate immune response

Centenarians as super-controls to assess the biological relevance of genetic risk factors for common age-related diseases: a proof of principle on type 2 diabetes

Genetic association studies of age-related, chronic human diseases often suffer from a lack of power to detect modest effects. Here we propose an alternative approach of including healthy centenarians as a more homogeneous and extreme control group. As a proof of principle we focused on type 2 diabetes (T2D) and assessed allelic/genotypic associations of 31 SNPs associated with T2D, diabetes complications and metabolic diseases and SNPs of genes relevant for telomere stability and age-related diseases. We hypothesized that the frequencies of risk variants are inversely correlated with decreasing health and longevity. We performed association analyses comparing diabetic patients and non-diabetic controls followed by association analyses with extreme phenotypic groups (T2D patients with complications and centenarians). Results drew attention to rs7903146 (TCF7L2 gene) that showed a constant increase in the frequencies of risk genotype (TT) from centenarians to diabetic patients who developed macro-complications and the strongest genotypic association was detected when diabetic patients were compared to centenarians (p_value = 9.066*10⁻⁷). We conclude that robust and biologically relevant associations can be obtained when extreme phenotypes, even with a small sample size, are compared.

In vitro infection of bovine monocytes with Mycoplasma bovis delays apoptosis and suppresses production of gamma interferon and tumor necrosis factor alpha but not interleukin-10

Mycoplasma bovis is one of the major causative pathogens of bovine respiratory complex disease (BRD), which is characterized by enzootic pneumonia, mastitis, pleuritis, and polyarthritis. M. bovis enters and colonizes bovine respiratory epithelial cells through inhalation of aerosol from contaminated air. The nature of the interaction between M. bovis and the bovine innate immune system is not well understood. We hypothesized that M. bovis invades blood monocytes and regulates cellular function to support its persistence and systemic dissemination. We used bovine-specific peptide kinome arrays to identify cellular signaling pathways that could be relevant to M. bovis-monocyte interactions in vitro. We validated these pathways using functional, protein, and gene expression assays. Here, we show that infection of bovine blood monocytes with M. bovis delays spontaneous or tumor necrosis factor alpha (TNF-α)/staurosporine-driven apoptosis, activates the NF-κB p65 subunit, and inhibits caspase-9 activity. We also report that M. bovis-infected bovine monocytes do not produce gamma interferon (IFN-γ) and TNF-α, although the level of production of interleukin-10 (IL-10) is elevated. Our findings suggest that M. bovis takes over the cellular machinery of bovine monocytes to prolong bacterial survival and to possibly facilitate subsequent systemic distribution.

MicroRNA expression profile in human macrophages in response to Leishmania major infection

Background

Leishmania (L.) are intracellular protozoan parasites able to survive and replicate in the hostile phagolysosomal environment of infected macrophages. They cause leishmaniasis, a heterogeneous group of worldwide-distributed affections, representing a paradigm of neglected diseases that are mainly embedded in impoverished populations. To establish successful infection and ensure their own survival, Leishmania have developed sophisticated strategies to subvert the host macrophage responses. Despite a wealth of gained crucial information, these strategies still remain poorly understood. MicroRNAs (miRNAs), an evolutionarily conserved class of endogenous 22-nucleotide non-coding RNAs, are described to participate in the regulation of almost every cellular process investigated so far. They regulate the expression of target genes both at the levels of mRNA stability and translation; changes in their expression have a profound effect on their target transcripts.

Methodology/Principal Findings

We report in this study a comprehensive analysis of miRNA expression profiles in L. major-infected human primary macrophages of three healthy donors assessed at different time-points post-infection (three to 24 h). We show that expression of 64 out of 365 analyzed miRNAs was consistently deregulated upon infection with the same trends in all donors. Among these, several are known to be induced by TLR-dependent responses. GO enrichment analysis of experimentally validated miRNA-targeted genes revealed that several pathways and molecular functions were disturbed upon parasite infection. Finally, following parasite infection, miR-210 abundance was enhanced in HIF-1α-dependent manner, though it did not contribute to inhibiting anti-apoptotic pathways through pro-apoptotic caspase-3 regulation.

Conclusions/Significance

Our data suggest that alteration in miRNA levels likely plays an important role in regulating macrophage functions following L. major infection. These results could contribute to better understanding of the dynamics of gene expression in host cells during leishmaniasis.

Leishmania parasites belong to different species, each one characterized by specific vectors and reservoirs, and causes cutaneous or visceral disease(s) of variable clinical presentation and severity. In its mammalian host, the parasite is an obligate intracellular pathogen infecting the monocyte/macrophage lineage. Leishmania have developed ambiguous relationships with macrophages. Indeed, these cells are the shelter of invading parasites, where they will grow and eventually will reside in a silent state for life. But macrophages are also the cells that participate, through the induction of several pro-inflammatory mediators and antigen presentation, to shape the host immune response and ultimately kill the invader. To subvert these anti-parasite responses, Leishmania manipulate the host machinery for their own differentiation and survival. We aimed to evaluate the impact of L. major (the causative agent of zoonotic cutaneous leishmaniasis) infection on deregulation of non-coding miRNAs, a class of important regulators of gene expression. Our results revealed the implication of several miRNAs on macrophage fate upon parasite infection through regulation of different pathways, including cell death. Our findings provided a new insight for understanding mechanisms governing this miRNA deregulation by parasite infection and will help to provide clues for the development of control strategies for this disease.

Identification of 526 conserved metazoan genetic innovations exposes a new role for cofactor E-like in neuronal microtubule homeostasis

The evolution of metazoans from their choanoflagellate-like unicellular ancestor coincided with the acquisition of novel biological functions to support a multicellular lifestyle, and eventually, the unique cellular and physiological demands of differentiated cell types such as those forming the nervous, muscle and immune systems. In an effort to understand the molecular underpinnings of such metazoan innovations, we carried out a comparative genomics analysis for genes found exclusively in, and widely conserved across, metazoans. Using this approach, we identified a set of 526 core metazoan-specific genes (the ‘metazoanome’), approximately 10% of which are largely uncharacterized, 16% of which are associated with known human disease, and 66% of which are conserved in Trichoplax adhaerens, a basal metazoan lacking neurons and other specialized cell types. Global analyses of previously-characterized core metazoan genes suggest a prevalent property, namely that they act as partially redundant modifiers of ancient eukaryotic pathways. Our data also highlights the importance of exaptation of pre-existing genetic tools during metazoan evolution. Expression studies in C. elegans revealed that many metazoan-specific genes, including tubulin folding cofactor E-like (TBCEL/coel-1), are expressed in neurons. We used C. elegans COEL-1 as a representative to experimentally validate the metazoan-specific character of our dataset. We show that coel-1 disruption results in developmental hypersensitivity to the microtubule drug paclitaxel/taxol, and that overexpression of coel-1 has broad effects during embryonic development and perturbs specialized microtubules in the touch receptor neurons (TRNs). In addition, coel-1 influences the migration, neurite outgrowth and mechanosensory function of the TRNs, and functionally interacts with components of the tubulin acetylation/deacetylation pathway. Together, our findings unveil a conserved molecular toolbox fundamental to metazoan biology that contains a number of neuronally expressed and disease-related genes, and reveal a key role for TBCEL/coel-1 in regulating microtubule function during metazoan development and neuronal differentiation.

The evolution of multicellular animals (metazoans) from their single-celled ancestor required new molecular tools to create and coordinate the various biological functions involved in a communal, or multicellular, lifestyle. This would eventually include the unique cellular and physiological demands of specialized tissues like the nervous system. To identify and understand the genetic bases of such unique metazoan traits, we used a comparative genomics approach to identify 526 metazoan-specific genes which have been evolutionarily conserved throughout the diversification of the animal kingdom. Interestingly, we found that some of those genes are still completely uncharacterized or poorly studied. We used the metazoan model organism C. elegans to examine the expression of some poorly characterized metazoan-specific genes and found that many, including one encoding tubulin folding cofactor E-like (TBCEL; C. elegans COEL-1), are expressed in cells of the nervous system. Using COEL-1 as an example to understand the metazoan-specific character of our dataset, our studies reveal a new role for this protein in regulating the stability of the microtubule cytoskeleton during development, and function of the touch receptor neurons. In summary, our findings help define a conserved molecular toolbox important for metazoan biology, and uncover an important role for COEL-1/TBCEL during development and in the nervous system of the metazoan C. elegans.

Small RNA profiling of influenza A virus-infected cells identifies miR-449b as a regulator of histone deacetylase 1 and interferon beta

The mammalian antiviral response relies on the alteration of cellular gene expression, to induce the production of antiviral effectors and regulate their activities. Recent research has indicated that virus infections can induce the accumulation of cellular microRNA (miRNA) species that influence the stability of host mRNAs and their protein products. To determine the potential for miRNA regulation of cellular responses to influenza A virus infection, small RNA profiling was carried out using next generation sequencing. Comparison of miRNA expression profiles in uninfected human A549 cells to cells infected with influenza A virus strains A/Udorn/72 and A/WSN/33, revealed virus-induced changes in miRNA abundance. Gene expression analysis identified mRNA targets for a cohort of highly inducible miRNAs linked to diverse cellular functions. Experiments demonstrate that the histone deacetylase, HDAC1, can be regulated by influenza-inducible miR-449b, resulting in altered mRNA and protein levels. Expression of miR-449b enhances virus and poly(I:C) activation of the IFNβ promoter, a process known to be negatively regulated by HDAC1. These findings demonstrate miRNA induction by influenza A virus infection and elucidate an example of miRNA control of antiviral gene expression in human cells, defining a role for miR-449b in regulation of HDAC1 and antiviral cytokine signaling.

Neonatal pain-related stress and NFKBIA genotype are associated with altered cortisol levels in preterm boys at school age

Neonatal pain-related stress is associated with elevated salivary cortisol levels to age 18 months in children born very preterm, compared to full-term, suggesting early programming effects. Importantly, interactions between immune/inflammatory and neuroendocrine systems may underlie programming effects. We examined whether cortisol changes persist to school age, and if common genetic variants in the promoter region of the NFKBIA gene involved in regulation of immune and inflammatory responses, modify the association between early experience and later life stress as indexed by hair cortisol levels, which provide an integrated index of endogenous HPA axis activity. Cortisol was assayed in hair samples from 128 children (83 born preterm ≤32 weeks gestation and 45 born full-term) without major sensory, motor or cognitive impairments at age 7 years. We found that hair cortisol levels were lower in preterm compared to term-born children. Downregulation of the HPA axis in preterm children without major impairment, seen years after neonatal stress terminated, suggests persistent alteration of stress system programming. Importantly, the etiology was gender-specific such that in preterm boys but not girls, specifically those with the minor allele for NFKBIA rs2233409, lower hair cortisol was associated with greater neonatal pain (number of skin-breaking procedures from birth to term), independent of medical confounders. Moreover, the minor allele (CT or TT) of NFKBIA rs2233409 was associated with higher secretion of inflammatory cytokines, supporting the hypothesis that neonatal pain-related stress may act as a proinflammatory stimulus that induces long-term immune cell activation. These findings are the first evidence that a long-term association between early pain-related stress and cortisol may be mediated by a genetic variants that regulate the activity of NF-κB, suggesting possible involvement of stress/inflammatory mechanisms in HPA programming in boys born very preterm.

Profiling microRNA expression in bovine alveolar macrophages using RNA-seq

MicroRNAs (miRNAs) are important regulators of gene expression and are known to play a key role in regulating both adaptive and innate immunity. Bovine alveolar macrophages (BAMs) help maintain lung homeostasis and constitute the front line of host defense against several infectious respiratory diseases, such as bovine tuberculosis. Little is known, however, about the role miRNAs play in these cells. In this study, we used a high-throughput sequencing approach, RNA-seq, to determine the expression levels of known and novel miRNAs in unchallenged BAMs isolated from lung lavages of eight different healthy Holstein-Friesian male calves. Approximately 80 million sequence reads were generated from eight BAM miRNA Illumina sequencing libraries, and 80 miRNAs were identified as being expressed in BAMs at a threshold of at least 100 reads per million (RPM). The expression levels of miRNAs varied over a large dynamic range, with a few miRNAs expressed at very high levels (up to 800,000RPM), and the majority lowly expressed. Notably, many of the most highly expressed miRNAs in BAMs have known roles in regulating immunity in other species (e.g. bta-let-7i, bta-miR-21, bta-miR-27, bta-miR-99b, bta-miR-146, bta-miR-147, bta-miR-155 and bta-miR-223). The most highly expressed miRNA in BAMs was miR-21, which has been shown to regulate the expression of antimicrobial peptides in Mycobacterium leprae-infected human monocytes. Furthermore, the predicted target genes of BAM-expressed miRNAs were found to be statistically enriched for roles in innate immunity. In addition to profiling the expression of known miRNAs, the RNA-seq data was also analysed to identify potentially novel bovine miRNAs. One putatively novel bovine miRNA was identified. To the best of our knowledge, this is the first RNA-seq study to profile miRNA expression in BAMs and provides an important reference dataset for investigating the regulatory roles miRNAs play in this important immune cell type.

Comparative phosphoproteomics reveals components of host cell invasion and post-transcriptional regulation during Francisella infection

Francisella tularensis is a facultative intracellular bacterium that causes the deadly disease tularemia. Most evidence suggests that Francisella is not well recognized by the innate immune system that normally leads to cytokine expression and cell death. In previous work, we identified new bacterial factors that were hyper-cytotoxic to macrophages. Four of the identified hyper-cytotoxic strains (lpcC, manB, manC, and kdtA) had an impaired lipopolysaccharide (LPS) synthesis and produced an exposed lipid A lacking the O-antigen. These mutants were not only hyper-cytotoxic but also were phagocytosed at much higher rates compared with the wild type parent strain. To elucidate the cellular signaling underlying this enhanced phagocytosis and cell death, we performed a large-scale comparative phosphoproteomic analysis of cells infected with wild-type and delta-lpcC F. novicida. Our data suggest that not only actin but also intermediate filaments and microtubules are important for F. novicida entry into the host cells. In addition, we observed differential phosphorylation of tristetraprolin, a key component of the mRNA-degrading machinery that controls the expression of a variety of genes including many cytokines. Infection with the delta-lpcC mutant induced the hyper-phosphorylation and inhibition of tristetraprolin, leading to the production of cytokines such as IL-1beta and TNF-alpha that may kill the host cells by triggering apoptosis. Together, our data provide new insights for Francisella invasion and a post-transcriptional mechanism that prevents the expression of host immune response factors that control infection by this pathogen.

A miniaturized chemical proteomic approach for target profiling of clinical kinase inhibitors in tumor biopsies

While targeted therapy based on the idea of attenuating the activity of a preselected, therapeutically relevant protein has become one of the major trends in modern cancer therapy, no truly specific targeted drug has been developed and most clinical agents have displayed a degree of polypharmacology. Therefore, the specificity of anticancer therapeutics has emerged as a highly important but severely underestimated issue. Chemical proteomics is a powerful technique combining postgenomic drug-affinity chromatography with high-end mass spectrometry analysis and bioinformatic data processing to assemble a target profile of a desired therapeutic molecule. Due to high demands on the starting material, however, chemical proteomic studies have been mostly limited to cancer cell lines. Herein, we report a down-scaling of the technique to enable the analysis of very low abundance samples, as those obtained from needle biopsies. By a systematic investigation of several important parameters in pull-downs with the multikinase inhibitor bosutinib, the standard experimental protocol was optimized to 100 μg protein input. At this level, more than 30 well-known targets were detected per single pull-down replicate with high reproducibility. Moreover, as presented by the comprehensive target profile obtained from miniaturized pull-downs with another clinical drug, dasatinib, the optimized protocol seems to be extendable to other drugs of interest. Sixty distinct human and murine targets were finally identified for bosutinib and dasatinib in chemical proteomic experiments utilizing core needle biopsy samples from xenotransplants derived from patient tumor tissue. Altogether, the developed methodology proves robust and generic and holds many promises for the field of personalized health care.

Experimental characterization of the human non-sequence-specific nucleic acid interactome

Background

The interactions between proteins and nucleic acids have a fundamental function in many biological processes, including gene transcription, RNA homeostasis, protein translation and pathogen sensing for innate immunity. While our knowledge of the ensemble of proteins that bind individual mRNAs in mammalian cells has been greatly augmented by recent surveys, no systematic study on the non-sequence-specific engagement of native human proteins with various types of nucleic acids has been reported.

Results

We designed an experimental approach to achieve broad coverage of the non-sequence-specific RNA and DNA binding space, including methylated cytosine, and tested for interaction potential with the human proteome. We used 25 rationally designed nucleic acid probes in an affinity purification mass spectrometry and bioinformatics workflow to identify proteins from whole cell extracts of three different human cell lines. The proteins were profiled for their binding preferences to the different general types of nucleic acids. The study identified 746 high-confidence direct binders, 139 of which were novel and 237 devoid of previous experimental evidence. We could assign specific affinities for sub-types of nucleic acid probes to 219 distinct proteins and individual domains. The evolutionarily conserved protein YB-1, previously associated with cancer and drug resistance, was shown to bind methylated cytosine preferentially, potentially conferring upon YB-1 an epigenetics-related function.

Conclusions

The dataset described here represents a rich resource of experimentally determined nucleic acid-binding proteins, and our methodology has great potential for further exploration of the interface between the protein and nucleic acid realms.

Systems biological approaches to measure and understand vaccine immunity in humans

Recent studies have demonstrated the utility of using systems approaches to identify molecular signatures that can be used to predict vaccine immunity in humans. Such approaches are now being used extensively in vaccinology, and are beginning to yield novel insights about the molecular networks driving vaccine immunity. In this review, we present a broad review of the methodologies involved in these studies, and discuss the promise and challenges involved in this emerging field of "systems vaccinology."

MicroRNA profiling of Sendai virus-infected A549 cells identifies miR-203 as an interferon-inducible regulator of IFIT1/ISG56

The mammalian type I interferon (IFN) response is a primary barrier for virus infection and is essential for complete innate and adaptive immunity. Both IFN production and IFN-mediated antiviral signaling are the result of differential cellular gene expression, a process that is tightly controlled at transcriptional and translational levels. To determine the potential for microRNA (miRNA)-mediated regulation of the antiviral response, small-RNA profiling was used to analyze the miRNA content of human A549 cells at steady state and following infection with the Cantell strain of Sendai virus, a potent inducer of IFN and cellular antiviral responses. While the miRNA content of the cells was largely unaltered by infection, specific changes in miRNA abundance were identified during Sendai virus infection. One miRNA, miR-203, was found to accumulate in infected cells and in response to IFN treatment. Results indicate that miR-203 is an IFN-inducible miRNA that can negatively regulate a number of cellular mRNAs, including an IFN-stimulated gene target, IFIT1/ISG56, by destabilizing its mRNA transcript.

Complementing the Eukaryotic Protein Interactome

Protein interaction networks are important for the understanding of regulatory mechanisms, for the explanation of experimental data and for the prediction of protein functions. Unfortunately, most interaction data is available only for model organisms. As a possible remedy, the transfer of interactions to organisms of interest is common practice, but it is not clear when interactions can be transferred from one organism to another and, thus, the confidence in the derived interactions is low. Here, we propose to use a rich set of features to train Random Forests in order to score transferred interactions. We evaluated the transfer from a range of eukaryotic organisms to S. cerevisiae using orthologs. Directly transferred interactions to S. cerevisiae are on average only 24% consistent with the current S. cerevisiae interaction network. By using commonly applied filter approaches the transfer precision can be improved, but at the cost of a large decrease in the number of transferred interactions. Our Random Forest approach uses various features derived from both the target and the source network as well as the ortholog annotations to assign confidence values to transferred interactions. Thereby, we could increase the average transfer consistency to 85%, while still transferring almost 70% of all correctly transferable interactions. We tested our approach for the transfer of interactions to other species and showed that our approach outperforms competing methods for the transfer of interactions to species where no experimental knowledge is available. Finally, we applied our predictor to score transferred interactions to 83 targets species and we were able to extend the available interactome of B. taurus, M. musculus and G. gallus with over 40,000 interactions each. Our transferred interaction networks are publicly available via our web interface, which allows to inspect and download transferred interaction sets of different sizes, for various species, and at specified expected precision levels.

AVAILABILITY

http://services.bio.ifi.lmu.de/coin-db/.

Optimization criteria and biological process enrichment in homologous multiprotein modules

Biological process enrichment is a widely used metric for evaluating the quality of multiprotein modules. In this study, we examine possible optimization criteria for detecting homologous multiprotein modules and quantify their effects on biological process enrichment. We find that modularity, linear density, and module size are the most important criteria considered, complementary to each other, and that graph theoretical attributes account for 36% of the variance in biological process enrichment. Variations in protein interaction similarity within module pairs have only minor effects on biological process enrichment. As random modules increase in size, both biological process enrichment and modularity tend to improve, although modularity does not show this upward trend in modules with size at most 50 proteins. To adjust for these trends, we recommend a size correction based on random sampling of modules when using biological process enrichment or other attributes to evaluate module boundaries. Characteristics of homologous multiprotein modules optimized for each of the optimization criteria are examined.

Differential response of bovine mammary epithelial cells to Staphylococcus aureus or Escherichia coli agonists of the innate immune system

Mastitis caused by Escherichia coli and Staphylococcus aureus is a major pathology of dairy cows. To better understand the differential response of the mammary gland to these two pathogens, we stimulated bovine mammary epithelial cells (bMEC) with either E. coli crude lipopolysaccharide (LPS) or with S. aureus culture supernatant (SaS) to compare the transcriptomic profiles of the initial bMEC response. By using HEK 293 reporter cells for pattern recognition receptors, the LPS preparation was found to stimulate TLR2 and TLR4 but not TLR5, Nod1 or Nod2, whereas SaS stimulated TLR2. Biochemical analysis revealed that lipoteichoic acid, protein A and α-hemolysin were all present in SaS, and bMEC were found to be responsive to each of these molecules. Transcriptome profiling revealed a core innate immune response partly shared by LPS and SaS. However, LPS induced expression of a significant higher number of genes and the fold changes were of greater magnitude than those induced by SaS. Microarray data analysis suggests that the activation pathways and the early chemokine and cytokine production preceded the defense and stress responses. A major differential response was the activation of the type I IFN pathway by LPS but not by SaS. The higher upregulation of chemokines (Cxcl10, Ccl2, Ccl5 and Ccl20) that target mononuclear leucocytes by LPS than by SaS is likely to be related to the differential activation of the type I IFN pathway, and could induce a different profile of the initial recruitment of leucocytes. The MEC responses to the two stimuli were different, as LPS was associated with NF-κB and Fas signaling pathways, whereas SaS was associated with AP-1 and IL-17A signaling pathways. It is noteworthy that at the protein level secretion of TNF-α and IL-1β was not induced by either stimulus. These results suggest that the response of MEC to diffusible stimuli from E. coli and S. aureus contributes to the onset of the response with differential leucocyte recruitment and distinct inflammatory and innate immune reactions of the mammary gland to infection.

In vivo phosphoproteomics analysis reveals the cardiac targets of β-adrenergic receptor signaling

β-Blockers are widely used to prevent cardiac arrhythmias and to treat hypertension by inhibiting β-adrenergic receptors (βARs) and thus decreasing contractility and heart rate. βARs initiate phosphorylation-dependent signaling cascades, but only a small number of the target proteins are known. We used quantitative in vivo phosphoproteomics to identify 670 site-specific phosphorylation changes in murine hearts in response to acute treatment with specific βAR agonists. The residues adjacent to the regulated phosphorylation sites exhibited a sequence-specific preference (R-X-X-pS/T), and integrative analysis of sequence motifs and interaction networks suggested that the kinases AMPK (adenosine 5'-monophosphate-activated protein kinase), Akt, and mTOR (mammalian target of rapamycin) mediate βAR signaling, in addition to the well-established pathways mediated by PKA (cyclic adenosine monophosphate-dependent protein kinase) and CaMKII (calcium/calmodulin-dependent protein kinase type II). We found specific regulation of phosphorylation sites on six ion channels and transporters that mediate increased ion fluxes at higher heart rates, and we showed that phosphorylation of one of these, Ser(92) of the potassium channel KV7.1, increased current amplitude. Our data set represents a quantitative analysis of phosphorylated proteins regulated in vivo upon stimulation of seven-transmembrane receptors, and our findings reveal previously unknown phosphorylation sites that regulate myocardial contractility, suggesting new potential targets for the treatment of heart disease and hypertension.

The comprehensive antibiotic resistance database

The field of antibiotic drug discovery and the monitoring of new antibiotic resistance elements have yet to fully exploit the power of the genome revolution. Despite the fact that the first genomes sequenced of free living organisms were those of bacteria, there have been few specialized bioinformatic tools developed to mine the growing amount of genomic data associated with pathogens. In particular, there are few tools to study the genetics and genomics of antibiotic resistance and how it impacts bacterial populations, ecology, and the clinic. We have initiated development of such tools in the form of the Comprehensive Antibiotic Research Database (CARD; http://arpcard.mcmaster.ca). The CARD integrates disparate molecular and sequence data, provides a unique organizing principle in the form of the Antibiotic Resistance Ontology (ARO), and can quickly identify putative antibiotic resistance genes in new unannotated genome sequences. This unique platform provides an informatic tool that bridges antibiotic resistance concerns in health care, agriculture, and the environment.

The alveolar microenvironment of patients infected with human immunodeficiency virus does not modify alveolar macrophage interactions with Streptococcus pneumoniae

We tested the hypothesis that HIV infection results in activation of alveolar macrophages and that this might be associated with impaired defense against pneumococcus. We compared alveolar macrophages and lymphocytes in 131 bronchoalveolar lavage samples from HIV-infected and healthy controls using inflammatory gene microarrays, flow cytometry, real-time PCR, and enzyme-linked immunosorbent assay (ELISA) to determine the pattern of macrophage activation associated with HIV infection and the effect of this activation on defense against pneumococcus. We used gamma interferon (IFN-γ) priming to mimic the cellular milieu in HIV-infected lungs. InnateDB and BioLayout 3D were used to analyze the interactions of the upregulated genes. Alveolar macrophages from HIV-infected adults showed increased gene expression and cytokine production in a classical pattern. Bronchoalveolar lavage from HIV-infected subjects showed excess CD8(+) lymphocytes with activated phenotype. Toll-like receptor 4 (TLR4) expression was increased in macrophages from HIV-infected subjects, but function was similar between the groups; lung lavage fluid did not inhibit TLR function in transfected HeLa cells. Alveolar macrophages from HIV-infected subjects showed normal binding and internalization of opsonized pneumococci, with or without IFN-γ priming. Alveolar macrophages from HIV-infected subjects showed classical activation compared to that of healthy controls, but this does not alter macrophage interactions with pneumococci.

Common risk alleles for inflammatory diseases are targets of recent positive selection

Genome-wide association studies (GWASs) have identified hundreds of loci harboring genetic variation influencing inflammatory-disease susceptibility in humans. It has been hypothesized that present day inflammatory diseases may have arisen, in part, due to pleiotropic effects of host resistance to pathogens over the course of human history, with significant selective pressures acting to increase host resistance to pathogens. The extent to which genetic factors underlying inflammatory-disease susceptibility has been influenced by selective processes can now be quantified more comprehensively than previously possible. To understand the evolutionary forces that have shaped inflammatory-disease susceptibility and to elucidate functional pathways affected by selection, we performed a systems-based analysis to integrate (1) published GWASs for inflammatory diseases, (2) a genome-wide scan for signatures of positive selection in a population of European ancestry, (3) functional genomics data comprised of protein-protein interaction networks, and (4) a genome-wide expression quantitative trait locus (eQTL) mapping study in peripheral blood mononuclear cells (PBMCs). We demonstrate that loci for inflammatory-disease susceptibility are enriched for genomic signatures of recent positive natural selection, with selected loci forming a highly interconnected protein-protein interaction network. Further, we identify 21 loci for inflammatory-disease susceptibility that display signatures of recent positive selection, of which 13 also show evidence of cis-regulatory effects on genes within the associated locus. Thus, our integrated analyses highlight a set of susceptibility loci that might subserve a shared molecular function and has experienced selective pressure over the course of human history; today, these loci play a key role in influencing susceptibility to multiple different inflammatory diseases, in part through alterations of gene expression in immune cells.

Functional genetic variation in NFKBIA and susceptibility to childhood asthma, bronchiolitis, and bronchopulmonary dysplasia

Respiratory diseases are the most frequent chronic illnesses in babies and children. Although a vigorous innate immune system is critical for maintaining lung health, a balanced response is essential to minimize damaging inflammation. We investigated the functional and clinical impact of human genetic variants in the promoter of NFKBIA, which encodes IκBα, the major negative regulator of NF-κB. In this study, we quantified the functional impact of NFKBIA promoter polymorphisms (rs3138053, rs2233406, and rs2233409) on promoter-driven protein expression, allele-specific and total NFKBIA mRNA expression, IκBα protein expression, and TLR responsiveness; mapped innate immune regulatory networks active during respiratory syncytial virus infection, asthma, and bronchopulmonary dysplasia; and genotyped and analyzed independent cohorts of children with respiratory syncytial virus infection, asthma, and bronchopulmonary dysplasia. Genetic variants in the promoter of NFKBIA influenced NFKBIA gene expression, IκBα protein expression, and TLR-mediated inflammatory responses. Using a systems biology approach, we demonstrated that NFKBIA/IκBα is a central hub in transcriptional responses of prevalent childhood lung diseases, including respiratory syncytial virus infection, asthma, and bronchopulmonary dysplasia. Finally, by examining independent pediatric lung disease cohorts, we established that this immunologically relevant genetic variation in the promoter of NFKBIA is associated with differential susceptibility to severe bronchiolitis following infection with respiratory syncytial virus, airway hyperresponsiveness, and severe bronchopulmonary dysplasia. These data highlight the importance of negative innate immune regulators, such as NFKBIA, in pediatric lung disease and begin to unravel common aspects in the genetic predisposition to bronchopulmonary dysplasia, bronchiolitis, and childhood asthma.

Next generation sequencing reveals the expression of a unique miRNA profile in response to a gram-positive bacterial infection

MicroRNAs (miRNAs) are short, non-coding RNAs, which post-transcriptionally regulate gene expression and are proposed to play a key role in the regulation of innate and adaptive immunity. Here, we report a next generation sequencing (NGS) approach profiling the expression of miRNAs in primary bovine mammary epithelial cells (BMEs) at 1, 2, 4 and 6 hours post-infection with Streptococcus uberis, a causative agent of bovine mastitis. Analysing over 450 million sequencing reads, we found that 20% of the approximately 1,300 currently known bovine miRNAs are expressed in unchallenged BMEs. We also identified the expression of more than 20 potentially novel bovine miRNAs. There is, however, a significant dynamic range in the expression of known miRNAs. The top 10 highly expressed miRNAs account for >80% of all aligned reads, with the remaining miRNAs showing much lower expression. Twenty-one miRNAs were identified as significantly differentially expressed post-infection with S. uberis. Several of these miRNAs have characterised roles in the immune systems of other species. This miRNA response to the Gram-positive S. uberis is markedly different, however, to lipopolysaccharide (LPS) induced miRNA expression. Of 145 miRNAs identified in the literature as being LPS responsive, only 9 were also differentially expressed in response to S. uberis. Computational analysis has also revealed that the predicted target genes of miRNAs, which are down-regulated in BMEs following S. uberis infection, are statistically enriched for roles in innate immunity. This suggests that miRNAs, which potentially act as central regulators of gene expression responses to a Gram-positive bacterial infection, may significantly regulate the sentinel capacity of mammary epithelial cells to mobilise the innate immune system.

A survey of protein interaction data and multigenic inherited disorders

Background

Multigenic diseases are often associated with protein complexes or interactions involved in the same pathway. We wanted to estimate to what extent this is true given a consolidated protein interaction data set. The study stresses data integration and data representation issues.

Results

We constructed 497 multigenic disease groups from OMIM and tested for overlaps with interaction and pathway data. A total of 159 disease groups had significant overlaps with protein interaction data consolidated by iRefIndex. A further 68 disease overlaps were found only in the KEGG pathway database. No single database contained all significant overlaps thus stressing the importance of data integration. We also found that disease groups overlapped with all three interaction data types: n-ary, spoke-represented complexes and binary data – thus stressing the importance of considering each of these data types separately.

Conclusions

Almost half of our multigenic disease groups could potentially be explained by protein complexes and pathways. However, the fact that no database or data type was able to cover all disease groups suggests that no single database has systematically covered all disease groups for potential related complex and pathway data. This survey provides a basis for further curation efforts to confirm and search for overlaps between diseases and interaction data. The accompanying R script can be used to reproduce the work and track progress in this area as databases change. Disease group overlaps can be further explored using the iRefscape plugin for Cytoscape.

Genome-wide expression profiling of schizophrenia using a large combined cohort

Numerous studies have examined gene expression profiles in post-mortem human brain samples from individuals with schizophrenia compared with healthy controls, to gain insight into the molecular mechanisms of the disease. Although some findings have been replicated across studies, there is a general lack of consensus on which genes or pathways are affected. It has been unclear if these differences are due to the underlying cohorts or methodological considerations. Here, we present the most comprehensive analysis to date of expression patterns in the prefrontal cortex of schizophrenic, compared with unaffected controls. Using data from seven independent studies, we assembled a data set of 153 affected and 153 control individuals. Remarkably, we identified expression differences in the brains of schizophrenics that are validated by up to seven laboratories using independent cohorts. Our combined analysis revealed a signature of 39 probes that are upregulated in schizophrenia and 86 that are downregulated. Some of these genes were previously identified in studies that were not included in our analysis, while others are novel to our analysis. In particular, we observe gene expression changes associated with various aspects of neuronal communication and alterations of processes affected as a consequence of changes in synaptic functioning. A gene network analysis predicted previously unidentified functional relationships among the signature genes. Our results provide evidence for a common underlying expression signature in this heterogeneous disorder.

Identification of a nonsense mutation in CWC15 associated with decreased reproductive efficiency in Jersey cattle

With the recent advent of genomic tools for cattle, several recessive conditions affecting fertility have been identified and selected against, such as deficiency of uridine monophosphate synthase, complex vertebral malformation, and brachyspina. The current report refines the location of a recessive haplotype affecting fertility in Jersey cattle using crossover haplotypes, discovers the causative mutation using whole genome sequencing, and examines the gene's role in embryo loss. In an attempt to identify unknown recessive lethal alleles in the current dairy population, a search using deep Mendelian sampling of 5,288 Jersey cattle was conducted for high-frequency haplotypes that have a deficit of homozygotes at the population level. This search led to the discovery of a putative recessive lethal in Jersey cattle on Bos taurus autosome 15. The haplotype, denoted JH1, was associated with reduced fertility, and further investigation identified one highly-influential Jersey bull as the putative source ancestor. By combining SNP analysis of whole-genome sequences aligned to the JH1 interval and subsequent SNP validation a nonsense mutation in CWC15 was identified as the likely causative mutation underlying the fertility phenotype. No homozygous recessive individuals were found in 749 genotyped animals, whereas all known carriers and carrier haplotypes possessed one copy of the mutant allele. This newly identified lethal has been responsible for a substantial number of spontaneous abortions in Jersey dairy cattle throughout the past half-century. With the mutation identified, selection against the deleterious allele in breeding schemes will aid in reducing the incidence of this defect in the population. These results also show that carrier status can be imputed with high accuracy. Whole-genome resequencing proved to be a powerful strategy to rapidly identify a previously mapped deleterious mutation in a known carrier of a recessive lethal allele.

Synthetic cationic peptide IDR-1018 modulates human macrophage differentiation

Macrophages play a critical role in the innate immune response. To respond in a rapid and efficient manner to challenges in the micro-environment, macrophages are able to differentiate towards classically (M1) or alternatively (M2) activated phenotypes. Synthetic, innate defense regulators (IDR) peptides, designed based on natural host defence peptides, have enhanced immunomodulatory activities and reduced toxicity leading to protection in infection and inflammation models that is dependent on innate immune cells like monocytes/macrophages. Here we tested the effect of IDR-1018 on macrophage differentiation, a process essential to macrophage function and the immune response. Using transcriptional, protein and systems biology analysis, we observed that differentiation in the presence of IDR-1018 induced a unique signature of immune responses including the production of specific pro and anti-inflammatory mediators, expression of wound healing associated genes, and increased phagocytosis of apoptotic cells. Transcription factor IRF4 appeared to play an important role in promoting this IDR-1018-induced phenotype. The data suggests that IDR-1018 drives macrophage differentiation towards an intermediate M1-M2 state, enhancing anti-inflammatory functions while maintaining certain pro-inflammatory activities important to the resolution of infection. Synthetic peptides like IDR-1018, which act by modulating the immune system, could represent a powerful new class of therapeutics capable of treating the rising number of multidrug resistant infections as well as disorders associated with dysregulated immune responses.

Interaction of Salmonella typhimurium with dendritic cells derived from pluripotent embryonic stem cells

Using an in vitro differentiation protocol we isolated cells with the properties of dendritic cells (DCs) from immunologically refractive pluripotent murine embryonic stem cells (ESCs). These ES-derived dendritic cells (ESDCs) expressed cytokines and were able to present antigen to a T cell line. Infection of ESDCs with Salmonella Typhimurium stimulated the expression of immune cell markers and thousands of murine genes, many associated with the immune response. Consequently, this system provides a novel in vitro model, amenable to genetic modification, for monitoring host/pathogen interactions.

Enhanced susceptibility to Citrobacter rodentium infection in microRNA-155-deficient mice

MicroRNAs (miRNAs) are small noncoding molecules that control gene expression posttranscriptionally, with microRNA-155 (miR-155) one of the first to be implicated in immune regulation. Here, we show that miR-155-deficient mice are less able to eradicate a mucosal Citrobacter rodentium infection than wild-type C57BL/6 mice. miR-155-deficient mice exhibited prolonged colonization associated with a higher C. rodentium burden in gastrointestinal tissue and spread into systemic tissues. Germinal center formation and humoral immune responses against C. rodentium were severely impaired in infected miR-155-deficient mice. A similarly susceptible phenotype was observed in μMT mice reconstituted with miR-155-deficient B cells, indicating that miR-155 is required intrinsically for mediating protection against this predominantly luminal bacterial pathogen.

OrtholugeDB: a bacterial and archaeal orthology resource for improved comparative genomic analysis

Prediction of orthologs (homologous genes that diverged because of speciation) is an integral component of many comparative genomics methods. Although orthologs are more likely to have similar function versus paralogs (genes that diverged because of duplication), recent studies have shown that their degree of functional conservation is variable. Also, there are inherent problems with several large-scale ortholog prediction approaches. To address these issues, we previously developed Ortholuge, which uses phylogenetic distance ratios to provide more precise ortholog assessments for a set of predicted orthologs. However, the original version of Ortholuge required manual intervention and was not easily accessible; therefore, we now report the development of OrtholugeDB, available online at http://www.pathogenomics.sfu.ca/ortholugedb. OrtholugeDB provides ortholog predictions for completely sequenced bacterial and archaeal genomes from NCBI based on reciprocal best Basic Local Alignment Search Tool hits, supplemented with further evaluation by the more precise Ortholuge method. The OrtholugeDB web interface facilitates user-friendly and flexible ortholog analysis, from single genes to genomes, plus flexible data download options. We compare Ortholuge with similar methods, showing how it may more consistently identify orthologs with conserved features across a wide range of taxonomic distances. OrtholugeDB facilitates rapid, and more accurate, bacterial and archaeal comparative genomic analysis and large-scale ortholog predictions.

InnateDB: systems biology of innate immunity and beyond--recent updates and continuing curation

InnateDB (http://www.innatedb.com) is an integrated analysis platform that has been specifically designed to facilitate systems-level analyses of mammalian innate immunity networks, pathways and genes. In this article, we provide details of recent updates and improvements to the database. InnateDB now contains >196 000 human, mouse and bovine experimentally validated molecular interactions and 3000 pathway annotations of relevance to all mammalian cellular systems (i.e. not just immune relevant pathways and interactions). In addition, the InnateDB team has, to date, manually curated in excess of 18 000 molecular interactions of relevance to innate immunity, providing unprecedented insight into innate immunity networks, pathways and their component molecules. More recently, InnateDB has also initiated the curation of allergy- and asthma-related interactions. Furthermore, we report a range of improvements to our integrated bioinformatics solutions including web service access to InnateDB interaction data using Proteomics Standards Initiative Common Query Interface, enhanced Gene Ontology analysis for innate immunity, and the availability of new network visualizations tools. Finally, the recent integration of bovine data makes InnateDB the first integrated network analysis platform for this agriculturally important model organism.

Merck Ad5/HIV induces broad innate immune activation that predicts CD8⁺ T-cell responses but is attenuated by preexisting Ad5 immunity

To better understand how innate immune responses to vaccination can lead to lasting protective immunity, we used a systems approach to define immune signatures in humans over 1 wk following MRKAd5/HIV vaccination that predicted subsequent HIV-specific T-cell responses. Within 24 h, striking increases in peripheral blood mononuclear cell gene expression associated with inflammation, IFN response, and myeloid cell trafficking occurred, and lymphocyte-specific transcripts decreased. These alterations were corroborated by marked serum inflammatory cytokine elevations and egress of circulating lymphocytes. Responses of vaccinees with preexisting adenovirus serotype 5 (Ad5) neutralizing antibodies were strongly attenuated, suggesting that enhanced HIV acquisition in Ad5-seropositive subgroups in the Step Study may relate to the lack of appropriate innate activation rather than to increased systemic immune activation. Importantly, patterns of chemoattractant cytokine responses at 24 h and alterations in 209 peripheral blood mononuclear cell transcripts at 72 h were predictive of subsequent induction and magnitude of HIV-specific CD8(+) T-cell responses. This systems approach provides a framework to compare innate responses induced by vectors, as shown here by contrasting the more rapid, robust response to MRKAd5/HIV with that to yellow fever vaccine. When applied iteratively, the findings may permit selection of HIV vaccine candidates eliciting innate immune response profiles more likely to drive HIV protective immunity.

Effects of protein interaction data integration, representation and reliability on the use of network properties for drug target prediction

Background

Previous studies have noted that drug targets appear to be associated with higher-degree or higher-centrality proteins in interaction networks. These studies explicitly or tacitly make choices of different source databases, data integration strategies, representation of proteins and complexes, and data reliability assumptions. Here we examined how the use of different data integration and representation techniques, or different notions of reliability, may affect the efficacy of degree and centrality as features in drug target prediction.

Results

Fifty percent of drug targets have a degree of less than nine, and ninety-five percent have a degree of less than ninety. We found that drug targets are over-represented in higher degree bins – this relationship is only seen for the consolidated interactome and it is not dependent on n-ary interaction data or its representation. Degree acts as a weak predictive feature for drug-target status and using more reliable subsets of the data does not increase this performance. However, performance does increase if only cancer-related drug targets are considered. We also note that a protein’s membership in pathway records can act as a predictive feature that is better than degree and that high-centrality may be an indicator of a drug that is more likely to be withdrawn.

Conclusions

These results show that protein interaction data integration and cleaning is an important consideration when incorporating network properties as predictive features for drug-target status. The provided scripts and data sets offer a starting point for further studies and cross-comparison of methods.

The ConsensusPathDB interaction database: 2013 update

Knowledge of the various interactions between molecules in the cell is crucial for understanding cellular processes in health and disease. Currently available interaction databases, being largely complementary to each other, must be integrated to obtain a comprehensive global map of the different types of interactions. We have previously reported the development of an integrative interaction database called ConsensusPathDB (http://ConsensusPathDB.org) that aims to fulfill this task. In this update article, we report its significant progress in terms of interaction content and web interface tools. ConsensusPathDB has grown mainly due to the integration of 12 further databases; it now contains 215 541 unique interactions and 4601 pathways from overall 30 databases. Binary protein interactions are scored with our confidence assessment tool, IntScore. The ConsensusPathDB web interface allows users to take advantage of these integrated interaction and pathway data in different contexts. Recent developments include pathway analysis of metabolite lists, visualization of functional gene/metabolite sets as overlap graphs, gene set analysis based on protein complexes and induced network modules analysis that connects a list of genes through various interaction types. To facilitate the interactive, visual interpretation of interaction and pathway data, we have re-implemented the graph visualization feature of ConsensusPathDB using the Cytoscape.js library.

Altered Toll-like receptor 9 signaling in Mycobacterium avium subsp. paratuberculosis-infected bovine monocytes reveals potential therapeutic targets

Mycobacterium avium subsp. paratuberculosis is the causative agent of Johne's disease in cattle. The complex, multifaceted interaction of M. avium subsp. paratuberculosis with its host includes dampening the ability of infected cells to respond to stimuli that promote M. avium subsp. paratuberculosis clearance. By disrupting host defenses, M. avium subsp. paratuberculosis creates an intracellular environment that favors the establishment and maintenance of infection. Toll-like receptors (TLRs) are important sensors that initiate innate immune responses to microbial challenge and are also immunotherapeutic targets. For example, TLR9 contributes to host defense against M. avium subsp. paratuberculosis, and its agonists (CpG oligodeoxynucleotides [ODNs]) are under investigation for treatment of Johne's disease and other infections. Here we demonstrate that M. avium subsp. paratuberculosis infection changes the responsiveness of bovine monocytes to TLR9 stimulation. M. avium subsp. paratuberculosis inhibits classical TLR9-mediated responses despite a 10-fold increase in TLR9 expression and maintained uptake of CpG ODNs. Other TLR9-mediated responses, such as oxidative burst, which occur through noncanonical signaling, remain functional. Kinome analysis verifies that classic TLR9 signaling is blocked by M. avium subsp. paratuberculosis infection and that signaling instead proceeds through a Pyk2-mediated mechanism. Pyk2-mediated signaling does not hinder infection, as CpG ODNs fail to promote M. avium subsp. paratuberculosis clearance. Indeed, Pyk2 signaling appears to be an important aspect of M. avium subsp. paratuberculosis infection, as Pyk2 inhibitors significantly reduce the number of intracellular M. avium subsp. paratuberculosis bacteria. The actions of M. avium subsp. paratuberculosis on TLR9 signaling may represent a strategy to generate a host environment which is better suited for infection, revealing potential new targets for therapeutic intervention.

Modeling and dynamical analysis of virus-triggered innate immune signaling pathways

The investigation of the dynamics and regulation of virus-triggered innate immune signaling pathways at a system level will enable comprehensive analysis of the complex interactions that maintain the delicate balance between resistance to infection and viral disease. In this study, we developed a delayed mathematical model to describe the virus-induced interferon (IFN) signaling process by considering several key players in the innate immune response. Using dynamic analysis and numerical simulation, we evaluated the following predictions regarding the antiviral responses: (1) When the replication ratio of virus is less than 1, the infectious virus will be eliminated by the immune system’s defenses regardless of how the time delays are changed. (2) The IFN positive feedback regulation enhances the stability of the innate immune response and causes the immune system to present the bistability phenomenon. (3) The appropriate duration of viral replication and IFN feedback processes stabilizes the innate immune response. The predictions from the model were confirmed by monitoring the virus titer and IFN expression in infected cells. The results suggest that the balance between viral replication and IFN-induced feedback regulation coordinates the dynamical behavior of virus-triggered signaling and antiviral responses. This work will help clarify the mechanisms of the virus-induced innate immune response at a system level and provide instruction for further biological experiments.

Manipulation of innate immunity by a bacterial secreted peptide: lantibiotic nisin Z is selectively immunomodulatory

Innate immunity is triggered by a variety of bacterial molecules, resulting in both protective and potentially harmful pro-inflammatory responses. Further, innate immunity also provides a mechanism for the maintenance of homeostasis between the host immune system and symbiotic or non-pathogenic microorganisms. However, the bacterial factors that mediate these protective effects have been incompletely defined. Here, it was demonstrated that the lantiobiotic nisin Z is able to modulate host immune responses and mediate protective host immunity. Nisin Z induced the secretion of the chemokines MCP-1, IL-8 and Gro-α, and significantly reduced TNF-α induction in response to bacterial LPS in human PBMC. The results correlated with the ability of nisin Z to confer protection against both the Gram-positive organism Staphylococcus aureus, and the Gram-negatives Salmonella enterica sv. Typhimurium and Escherichia coli in murine challenge models. Mechanistic studies revealed that nisin Z modulates host immunity through similar mechanisms as natural host defense peptides, engaging multiple signal transduction pathways and growth factor receptors. The results presented herein demonstrate that, in addition to nisin Z, other bacterial cationic peptides and, in particular, the lantibiotics, could represent a new class of secreted bacterial molecule with immunomodulatory activities.

Atypical activation of the unfolded protein response in cystic fibrosis airway cells contributes to p38 MAPK-mediated innate immune responses

Inflammatory lung disease is the major cause of morbidity and mortality in cystic fibrosis (CF); understanding what produces dysregulated innate immune responses in CF cells will be pivotal in guiding the development of novel anti-inflammatory therapies. To elucidate the molecular mechanisms that mediate exaggerated inflammation in CF following TLR signaling, we profiled global gene expression in immortalized human CF and non-CF airway cells at baseline and after microbial stimulation. Using complementary analysis methods, we observed a signature of increased stress levels in CF cells, specifically characterized by endoplasmic reticulum (ER) stress, the unfolded protein response (UPR), and MAPK signaling. Analysis of ER stress responses revealed an atypical induction of the UPR, characterized by the lack of induction of the PERK-eIF2α pathway in three complementary model systems: immortalized CF airway cells, fresh CF blood cells, and CF lung tissue. This atypical pattern of UPR activation was associated with the hyperinflammatory phenotype in CF cells, as deliberate induction of the PERK-eIF2α pathway with salubrinal attenuated the inflammatory response to both flagellin and Pseudomonas aeruginosa. IL-6 production triggered by ER stress and microbial stimulation were both dependent on p38 MAPK activity, suggesting a molecular link between both signaling events. These data indicate that atypical UPR activation fails to resolve the ER stress in CF and sensitizes the innate immune system to respond more vigorously to microbial challenge. Strategies to restore ER homeostasis and normalize the UPR activation profile may represent a novel therapeutic approach to minimize lung-damaging inflammation in CF.

Genomic expression analyses reveal lysosomal, innate immunity proteins, as disease correlates in murine models of a lysosomal storage disorder

Niemann-Pick Type C (NPC) disease is a rare, genetic, lysosomal disorder with progressive neurodegeneration. Poor understanding of the pathophysiology and a lack of blood-based diagnostic markers are major hurdles in the treatment and management of NPC and several additional, neurological lysosomal disorders. To identify disease severity correlates, we undertook whole genome expression profiling of sentinel organs, brain, liver, and spleen of Balb/c Npc1^−/− mice relative to Npc1^+/− at an asymptomatic stage, as well as early- and late-symptomatic stages. Unexpectedly, we found prominent up regulation of innate immunity genes with age-dependent change in their expression, in all three organs. We shortlisted a set of 12 secretory genes whose expression steadily increased with age in both brain and liver, as potential plasma correlates of neurological and/or liver disease. Ten were innate immune genes with eight ascribed to lysosomes. Several are known to be elevated in diseased organs of murine models of other lysosomal diseases including Gaucher’s disease, Sandhoff disease and MPSIIIB. We validated the top candidate lysozyme, in the plasma of Npc1^−/− as well as Balb/c Npc1^nmf164 mice (bearing a point mutation closer to human disease mutants) and show its reduction in response to an emerging therapeutic. We further established elevation of innate immunity in Npc1^−/− mice through multiple functional assays including inhibition of bacterial infection as well as cellular analysis and immunohistochemistry. These data revealed neutrophil elevation in the Npc1^−/− spleen and liver (where large foci were detected proximal to damaged tissue). Together our results yield a set of lysosomal, secretory innate immunity genes that have potential to be developed as pan or specific plasma markers for neurological diseases associated with lysosomal storage and where diagnosis is a major problem. Further, the accumulation of neutrophils in diseased organs (hitherto not associated with NPC) suggests their role in pathophysiology and disease exacerbation.

Proteomic and bioinformatic profile of primary human oral epithelial cells

Wounding of the oral mucosa occurs frequently in a highly septic environment. Remarkably, these wounds heal quickly and the oral cavity, for the most part, remains healthy. Deciphering the normal human oral epithelial cell (NHOEC) proteome is critical for understanding the mechanism(s) of protection elicited when the mucosal barrier is intact, as well as when it is breached. Combining 2D gel electrophoresis with shotgun proteomics resulted in identification of 1662 NHOEC proteins. Proteome annotations were performed based on protein classes, molecular functions, disease association and membership in canonical and metabolic signaling pathways. Comparing the NHOEC proteome with a database of innate immunity-relevant interactions (InnateDB) identified 64 common proteins associated with innate immunity. Comparison with published salivary proteomes revealed that 738/1662 NHOEC proteins were common, suggesting that significant numbers of salivary proteins are of epithelial origin. Gene ontology analysis showed similarities in the distributions of NHOEC and saliva proteomes with regard to biological processes, and molecular functions. We also assessed the interindividual variability of the NHOEC proteome and observed it to be comparable with other primary cells. The baseline proteome described in this study should serve as a resource for proteome studies of the oral mucosa, especially in relation to disease processes.

Popular computational methods to assess multiprotein complexes derived from label-free affinity purification and mass spectrometry (AP-MS) experiments

Advances in sensitivity, resolution, mass accuracy, and throughput have considerably increased the number of protein identifications made via mass spectrometry. Despite these advances, state-of-the-art experimental methods for the study of protein-protein interactions yield more candidate interactions than may be expected biologically owing to biases and limitations in the experimental methodology. In silico methods, which distinguish between true and false interactions, have been developed and applied successfully to reduce the number of false positive results yielded by physical interaction assays. Such methods may be grouped according to: (1) the type of data used: methods based on experiment-specific measurements (e.g., spectral counts or identification scores) versus methods that extract knowledge encoded in external annotations (e.g., public interaction and functional categorisation databases); (2) the type of algorithm applied: the statistical description and estimation of physical protein properties versus predictive supervised machine learning or text-mining algorithms; (3) the type of protein relation evaluated: direct (binary) interaction of two proteins in a cocomplex versus probability of any functional relationship between two proteins (e.g., co-occurrence in a pathway, sub cellular compartment); and (4) initial motivation: elucidation of experimental data by evaluation versus prediction of novel protein-protein interaction, to be experimentally validated a posteriori. This work reviews several popular computational scoring methods and software platforms for protein-protein interactions evaluation according to their methodology, comparative strengths and weaknesses, data representation, accessibility, and availability. The scoring methods and platforms described include: CompPASS, SAINT, Decontaminator, MINT, IntAct, STRING, and FunCoup. References to related work are provided throughout in order to provide a concise but thorough introduction to a rapidly growing interdisciplinary field of investigation.

The integration of 'omic' disciplines and systems biology in cattle breeding

Enormous progress has been made in the selection of animals, including cattle, for specific traits using traditional quantitative genetics approaches. Nevertheless, considerable variation in phenotypes remains unexplained, and therefore represents potential additional gain for animal production. In addition, the paradigm shift in new disciplines now being applied to animal breeding represents a powerful opportunity to prise open the 'black box' underlying the response to selection and fully understand the genetic architecture controlling the traits of interest. A move away from traditional approaches of animal breeding toward systems approaches using integrative analysis of data from the 'omic' disciplines represents a multitude of exciting opportunities for animal breeding going forward as well as providing alternatives for overcoming some of the limitations of traditional approaches such as the expressed phenotype being an imperfect predictor of the individual's true genetic merit, or the phenotype being only expressed in one gender or late in the lifetime of an animal. This review aims to discuss these opportunities from the perspective of their potential application and contribution to cattle breeding. Harnessing the potential of this paradigm shift also poses some new challenges for animal scientists - and they will also be discussed.

Effective adjunctive therapy by an innate defense regulatory peptide in a preclinical model of severe malaria

Case fatality rates for severe malaria remain high even in the best clinical settings because antimalarial drugs act against the parasite without alleviating life-threatening inflammation. We assessed the potential for host-directed therapy of severe malaria of a new class of anti-inflammatory drugs, the innate defense regulator (IDR) peptides, based on host defense peptides. The Plasmodium berghei ANKA model of experimental cerebral malaria was adapted to use as a preclinical screen by combining late-stage intervention in established infections with advanced bioinformatic analysis of early transcriptional changes in co-regulated gene sets. Coadministration of IDR-1018 with standard first-line antimalarials increased survival of infected mice while down-regulating key inflammatory networks associated with fatality. Thus, IDR peptides provided host-directed adjunctive therapy for severe disease in combination with antimalarial treatment.

Viral immune modulators perturb the human molecular network by common and unique strategies

Viruses must enter host cells to replicate, assemble and propagate. Because of the restricted size of their genomes, viruses have had to evolve efficient ways of exploiting host cell processes to promote their own life cycles and also to escape host immune defence mechanisms. Many viral open reading frames (viORFs) with immune-modulating functions essential for productive viral growth have been identified across a range of viral classes. However, there has been no comprehensive study to identify the host factors with which these viORFs interact for a global perspective of viral perturbation strategies. Here we show that different viral perturbation patterns of the host molecular defence network can be deduced from a mass-spectrometry-based host-factor survey in a defined human cellular system by using 70 innate immune-modulating viORFs from 30 viral species. The 579 host proteins targeted by the viORFs mapped to an unexpectedly large number of signalling pathways and cellular processes, suggesting yet unknown mechanisms of antiviral immunity. We further experimentally verified the targets heterogeneous nuclear ribonucleoprotein U, phosphatidylinositol-3-OH kinase, the WNK (with-no-lysine) kinase family and USP19 (ubiquitin-specific peptidase 19) as vulnerable nodes in the host cellular defence system. Evaluation of the impact of viral immune modulators on the host molecular network revealed perturbation strategies used by individual viruses and by viral classes. Our data are also valuable for the design of broad and specific antiviral therapies.

Induction of ligand-specific PrP (C) signaling in human neuronal cells

Cellular prion protein (PrP (C) ) has attracted considerable attention for its role in transmissible spongiform encephalopathies (TSEs). In spite of being a point of intense research effort critical questions still remain regarding the physiological function of PrP (C) and how these functions may change with the conversion of the protein into the infectious and pathological conformation (PrP (Sc) ). While emerging evidence suggests PrP (C/Sc) are involved in signal transduction there is little consensus on the signaling pathways associated with the normal and diseased states. The purported involvement of PrP (C) in signal transduction, and the association of TSEs with neural pathology, makes kinome analysis of human neurons an interesting and appropriate model to characterize patterns of signal transduction following activation of PrP (C) by two commonly employed experimental ligands; antibody-induced dimerization by 6H4 and the amino acids 106-126 PrP peptide fragment (PrP 106-126). Analysis of the induced kinome responses reveals distinct patterns of signaling activity following each treatment. Specifically, stimulation of human neurons with the 6H4 antibody results in alterations in mitogen activated protein kinase (MAPK) signaling pathways while the 106-126 peptide activates growth factor related signaling pathways including vascular endothelial growth factor (VEGF) signaling and the phosphoinositide-3 kinase (PI3K) pathway. These pathways were validated through independent functional assays. Collectively these results indicate that stimulation of PrP (C) with distinct ligands, even within the same cell type, results in unique patterns of signaling. While this investigation highlights the apparent functional versatility of PrP (C) as a signaling molecule and may offer insight into cellular mechanisms of TSE pathology it also emphasizes the potential dangers associated with attributing activation of specific intracellular events to particular receptors through artificial models of receptor activation.

Proteomic analysis of lysine acetylation sites in rat tissues reveals organ specificity and subcellular patterns

Lysine acetylation is a major posttranslational modification involved in a broad array of physiological functions. Here, we provide an organ-wide map of lysine acetylation sites from 16 rat tissues analyzed by high-resolution tandem mass spectrometry. We quantify 15,474 modification sites on 4,541 proteins and provide the data set as a web-based database. We demonstrate that lysine acetylation displays site-specific sequence motifs that diverge between cellular compartments, with a significant fraction of nuclear sites conforming to the consensus motifs G-AcK and AcK-P. Our data set reveals that the subcellular acetylation distribution is tissue-type dependent and that acetylation targets tissue-specific pathways involved in fundamental physiological processes. We compare lysine acetylation patterns for rat as well as human skeletal muscle biopsies and demonstrate its general involvement in muscle contraction. Furthermore, we illustrate that acetylation of fructose-bisphosphate aldolase and glycerol-3-phosphate dehydrogenase serves as a cellular mechanism to switch off enzymatic activity.

Sets2Networks: network inference from repeated observations of sets

BACKGROUND

The skeleton of complex systems can be represented as networks where vertices represent entities, and edges represent the relations between these entities. Often it is impossible, or expensive, to determine the network structure by experimental validation of the binary interactions between every vertex pair. It is usually more practical to infer the network from surrogate observations. Network inference is the process by which an underlying network of relations between entities is determined from indirect evidence. While many algorithms have been developed to infer networks from quantitative data, less attention has been paid to methods which infer networks from repeated co-occurrence of entities in related sets. This type of data is ubiquitous in the field of systems biology and in other areas of complex systems research. Hence, such methods would be of great utility and value.

RESULTS

Here we present a general method for network inference from repeated observations of sets of related entities. Given experimental observations of such sets, we infer the underlying network connecting these entities by generating an ensemble of networks consistent with the data. The frequency of occurrence of a given link throughout this ensemble is interpreted as the probability that the link is present in the underlying real network conditioned on the data. Exponential random graphs are used to generate and sample the ensemble of consistent networks, and we take an algorithmic approach to numerically execute the inference method. The effectiveness of the method is demonstrated on synthetic data before employing this inference approach to problems in systems biology and systems pharmacology, as well as to construct a co-authorship collaboration network. We predict direct protein-protein interactions from high-throughput mass-spectrometry proteomics, integrate data from Chip-seq and loss-of-function/gain-of-function followed by expression data to infer a network of associations between pluripotency regulators, extract a network that connects 53 cancer drugs to each other and to 34 severe adverse events by mining the FDA's Adverse Events Reporting Systems (AERS), and construct a co-authorship network that connects Mount Sinai School of Medicine investigators. The predicted networks and online software to create networks from entity-set libraries are provided online at http://www.maayanlab.net/S2N.

CONCLUSIONS

The network inference method presented here can be applied to resolve different types of networks in current systems biology and systems pharmacology as well as in other fields of research.

Influenza A virus infection of human respiratory cells induces primary microRNA expression

The cellular response to virus infection is initiated by recognition of the invading pathogen and subsequent changes in gene expression mediated by both transcriptional and translational mechanisms. In addition to well established means of regulating antiviral gene expression, it has been demonstrated that RNA interference (RNAi) can play an important role in antiviral responses. Virus-derived small interfering RNA (siRNA) is a primary antiviral response exploited by plants and invertebrate animals, and host-encoded microRNA (miRNA) species have been clearly implicated in the regulation of innate and adaptive immune responses in mammals and other vertebrates. Examination of miRNA abundance in human lung cell lines revealed endogenous miRNAs, including miR-7, miR-132, miR-146a, miR-187, miR-200c, and miR-1275, to specifically accumulate in response to infection with two influenza A virus strains, A/Udorn/72 and A/WSN/33. Known antiviral response pathways, including Toll-like receptor, RIG-I-like receptor, and direct interferon or cytokine stimulation did not alter the abundance of the tested miRNAs to the extent of influenza A virus infection, which initiates primary miRNA transcription via a secondary response pathway. Gene expression profiling identified 26 cellular mRNAs targeted by these miRNAs, including IRAK1, MAPK3, and other components of innate immune signaling systems.