Vikodak - A Modular Framework for Inferring Functional Potential of Microbial Communities from 16S Metagenomic Datasets

Background

The overall metabolic/functional potential of any given environmental niche is a function of the sum total of genes/proteins/enzymes that are encoded and expressed by various interacting microbes residing in that niche. Consequently, prior (collated) information pertaining to genes, enzymes encoded by the resident microbes can aid in indirectly (re)constructing/ inferring the metabolic/ functional potential of a given microbial community (given its taxonomic abundance profile). In this study, we present Vikodak—a multi-modular package that is based on the above assumption and automates inferring and/ or comparing the functional characteristics of an environment using taxonomic abundance generated from one or more environmental sample datasets. With the underlying assumptions of co-metabolism and independent contributions of different microbes in a community, a concerted effort has been made to accommodate microbial co-existence patterns in various modules incorporated in Vikodak.

Results

Validation experiments on over 1400 metagenomic samples have confirmed the utility of Vikodak in (a) deciphering enzyme abundance profiles of any KEGG metabolic pathway, (b) functional resolution of distinct metagenomic environments, (c) inferring patterns of functional interaction between resident microbes, and (d) automating statistical comparison of functional features of studied microbiomes. Novel features incorporated in Vikodak also facilitate automatic removal of false positives and spurious functional predictions.

Conclusions

With novel provisions for comprehensive functional analysis, inclusion of microbial co-existence pattern based algorithms, automated inter-environment comparisons; in-depth analysis of individual metabolic pathways and greater flexibilities at the user end, Vikodak is expected to be an important value addition to the family of existing tools for 16S based function prediction.

Availability and Implementation

A web implementation of Vikodak can be publicly accessed at: http://metagenomics.atc.tcs.com/vikodak. This web service is freely available for all categories of users (academic as well as commercial).     

IL-25, IL-33 and TSLP receptor are not critical for development of experimental murine malaria

IL-25, IL-33 and TSLP, which are produced predominantly by epithelial cells, can induce production of Th2-type cytokines such as IL-4, IL-5 and/or IL-13 by various types of cells, suggesting their involvement in induction of Th2-type cytokine-associated immune responses. It is known that Th2-type cytokines contribute to host defense against malaria parasite infection in mice. However, the roles of IL-25, IL-33 and TSLP in malaria parasite infection remain unclear. Thus, to elucidate this, we infected wild-type, IL-25^?/?, IL-33^?/? and TSLP receptor (TSLPR)^?/? mice with Plasmodium berghei (P. berghei) ANKA, a murine malaria strain. The expression levels of IL-25, IL-33 and TSLP mRNA were changed in the brain, liver, lung and spleen of wild-type mice after infection, suggesting that these cytokines are involved in host defense against P. berghei ANKA. However, the incidence of parasitemia and survival in the mutant mice were comparable to in the wild-type mice. These findings indicate that IL-25, IL-33 and TSLP are not critical for host defense against P. berghei ANKA.     

Sampling From the Proteome to the Human Leukocyte Antigen-DR (HLA-DR) Ligandome Proceeds Via High Specificity

Comprehensive analysis of the complex nature of the Human Leukocyte Antigen (HLA) class II ligandome is of utmost importance to understand the basis for CD4⁺ T cell mediated immunity and tolerance. Here, we implemented important improvements in the analysis of the repertoire of HLA-DR-presented peptides, using hybrid mass spectrometry-based peptide fragmentation techniques on a ligandome sample isolated from matured human monocyte-derived dendritic cells (DC). The reported data set constitutes nearly 14 thousand unique high-confident peptides, i.e. the largest single inventory of human DC derived HLA-DR ligands to date. From a technical viewpoint the most prominent finding is that no single peptide fragmentation technique could elucidate the majority of HLA-DR ligands, because of the wide range of physical chemical properties displayed by the HLA-DR ligandome. Our in-depth profiling allowed us to reveal a strikingly poor correlation between the source proteins identified in the HLA class II ligandome and the DC cellular proteome. Important selective sieving from the sampled proteome to the ligandome was evidenced by specificity in the sequences of the core regions both at their N- and C- termini, hence not only reflecting binding motifs but also dominant protease activity associated to the endolysosomal compartments. Moreover, we demonstrate that the HLA-DR ligandome reflects a surface representation of cell-compartments specific for biological events linked to the maturation of monocytes into antigen presenting cells. Our results present new perspectives into the complex nature of the HLA class II system and will aid future immunological studies in characterizing the full breadth of potential CD4⁺ T cell epitopes relevant in health and disease.Human Leukocyte Antigen (HLA)¹ class II molecules on professional antigen presenting cells such as dendritic cells (DC) expose peptide fragments derived from exogenous and endogenous proteins to be screened by CD4⁺ T cells (1, 2). The activation and recruitment of CD4⁺ T cells recognizing disease-related peptide antigens is critical for the development of efficient antipathogen or antitumor immunity. Furthermore, the presentation of self-peptides and their interaction with CD4⁺ T cells is essential to maintain immunological tolerance and homeostasis (3). Knowledge of the nature of HLA class II-presented peptides on DC is of great importance to understand the rules of antigen processing and peptide binding motifs (4), whereas the identity of disease-related antigens may provide new knowledge on immunogenicity and leads for the development of vaccines and immunotherapy (5, 6).Mass spectrometry (MS) has proven effective for the analysis HLA class II-presented peptides (4, 7, 8). MS-based ligandome studies have demonstrated that HLA class II molecules predominantly present peptides derived from exogenous proteins that entered the cells by endocytosis and endogenous proteins that are associated with the endo-lysosomal compartments (4). Yet proteins residing in the cytosol, nucleus or mitochondria can also be presented by HLA class II molecules, primarily through autophagy (9–11). Multiple studies have mapped the HLA class II ligandome of antigen presenting cells in the context of infectious pathogens (12), autoimmune diseases (13–17) or cancer (14, 18, 19), or those that are essential for self-tolerance in the human thymus (3, 20). Notwithstanding these efforts, and certainly not in line with the extensive knowledge on the HLA class I ligandome (21), the nature of the HLA class II-presented peptide repertoire and particular its relationship to the cellular source proteome remains poorly understood.To advance our knowledge on the HLA-DR ligandome on activated DC without having to deal with limitations in cell yield from peripheral human blood (12, 21, 22) or tissue isolates (3), we explored the use of MUTZ-3 cells. This cell line has been used as a model of human monocyte-derived DCs. MUTZ-3 cells can be matured to act as antigen presenting cells and express then high levels of HLA class II molecules, and can be propagated in vitro to large cell densities (23–25). We also evaluated the performance of complementary and hybrid MS fragmentation techniques electron-transfer dissociation (ETD), electron-transfer/higher-energy collision dissociation (EThcD) (26), and higher-energy collision dissociation (HCD) to sequence and identify the HLA class II ligandome. Together this workflow allowed for the identification of an unprecedented large set of about 14 thousand unique peptide sequences presented by DC derived HLA-DR molecules, providing an in-depth view of the complexity of the HLA class II ligandome, revealing underlying features of antigen processing and surface-presentation to CD4⁺ T cells.     

Characterization of the venom of the vermivorous cone snail Conus fulgetrum

Over 200 components with molecular mass ranging mainly from 400 to 4000 Da were characterized from the venom of the vermivorous cone snail Conus fulgetrum that inhabit Egyptian Red Sea. One major component having a molecular mass of 2946 Da was purified by HPLC, and its primary structure was determined by a combination of Edman degradation and MS/MS analysis.     

GRP78 regulates sensitivity of human colorectal cancer cells to DNA targeting agents

This study was carried out to investigate the activation status of unfolded protein response (UPR) in colorectal cancer (CRC) and its contribution to CRC resistance to chemotherapy-induced apoptosis. Chemotherapy-induced apoptosis was assessed by the propidium iodide method. Activation of UPR was evaluated in CRC cell lines using immunoblotting technique and in CRC tissues using immunohistochemistry. Findings of the present study revealed that the UPR is constitutively activated in CRC cell lines and CRC tissues isolated from patients, as evidenced by relatively high levels of the 78-kDa glucose-regulated protein (GRP78) and spliced X-box-binding protein 1 mRNA in tissue samples. In addition, CRC cell lines differentially responded to clinically relevant DNA-targeting agents including cisplatin, and 5-flourouracil. Moreover, the levels of GRP78 were inversely associated with sensitivity of CRC cells to chemotherapy-induced apoptosis. Inhibition of GRP78 by siRNA resulted in increased sensitivity of CRC cells to chemotherapeutic agents. Collectively, current results appear to provide novel insights into the role of UPR in determining sensitivity of CRC cells to chemotherapeutic agents and might have important implications for personalized CRC treatment.