Implications of Human Antimicrobial Peptide Defensin Beta-1 in Clinical Oral Squamous Cell Carcinoma Patients via an Integrated Bioinformatics Approach

BACKGROUND

The human antimicrobial peptide defensin beta 1 (DEFB1) has been found to play antimicrobial and anti-inflammatory roles in oral diseases; however, its tumor-regulating role in oral squamous cell carcinoma (OSCC) has not yet been researched by using an integrative bioinformatics approach.

OBJECTIVE

To investigate the regulating mechanisms of the DEFB1 gene in OSCC in terms of its expression patterns, prognostic values, biological functions, and implication for tumor immunity.

METHODS

The DEFB1 gene expression pattern and regulatory involvement in OSCC were investigated using publically accessible data from TCGA database. R software tools and public web servers were utilized to conduct statistical analysis of data from cancer and noncancerous samples.

RESULTS

DEFB1 was found to be significantly downregulated in OSCC tumor samples compared with healthy control oral samples. The DEFB1 gene was found associated with the prognostic outcomes of OSCC, and its upregulation represented better survival outcome. Gene set enrichment analysis (GSEA) results showed that DEFB1-significantly correlated genes were mainly enriched in four signaling pathways mediating the antitumor role of DEFB1 in OSCC, including extracellular matrix-related pathway, RTK/PI3K/AKT/mTOR pathway, keratinization, and cytokine-related pathway. The gene-gene interaction network showed that DEFB1 was closely correlated with several genes, for example, CCR6 (C-C motif chemokine receptor 6), CXCL1 (C-X-C motif chemokine ligand 1), MAP4K2 (mitogen-activated protein kinase kinase kinase kinase 2), PTGER3 (prostaglandin E receptor 3), and MMP7 (matrix metallopeptidase 7). Moreover, DEFB1 was found to be involved in the tumor immunity of OSCC by regulating the function of tumor macrophage cells, mast cells, T cells, and NK cells.

CONCLUSIONS

Given the dysregulation, prognostic value, and tumor progression-related biological pathway alteration, indicating the tumor immune-modulatory role of DEFB1 in OSCC, the DEFB1 gene should be regarded as a potential therapeutic target for treating oral cancer.

Identification of Therapeutic Targets and Prognostic Biomarkers among Genes from the Mediator Complex Family in the Hepatocellular Carcinoma Tumour-Immune Microenvironment

BACKGROUND

Hepatocellular carcinoma (HCC) is predominant among all types of primary liver cancers characterised by high morbidity and mortality. Genes in the mediator complex (MED) family are engaged in the tumour-immune microenvironment and function as regulatory hubs mediating carcinogenesis and progression across diverse cancer types. Whereas research studies have been conducted to examine the mechanisms in several cancers, studies that systematically focused on the therapeutic and prognostic values of MED in patients with HCC are limited.

METHODS

The online databases ONCOMINE, GEPIA, UALCAN, GeneMANIA, cBioPortal, OmicStudio, STING, Metascape, and TIMER were used in this study.

RESULTS

The transcriptional levels of all members of the MED family in HCC presented an aberrant high expression pattern. Significant correlations were found between the MED1, MED6, MED8, MED10, MED12, MED15, MED17, MED19, MED20, MED21, MED22, MED23, MED24, MED25, MED26, and MED27 expression levels and the pathological stage in the patients with HCC. The patients with high expression levels of MED6, MED8, MED10, MED17, MED19, MED20, MED21, MED22, MED24, and MED25 were significantly associated with poor prognosis. Functional enrichment analysis revealed that the members of the MED family were mainly enriched in the nucleobase-containing compound catabolic process, regulation of chromosome organisation, and transcriptional regulation by TP53. Significant correlations were found between the MED6, MED8, MED10, MED17, MED19, MED20, MED21, MED22, MED24, and MED25 expression levels and all types of immune cells (B cells, CD8+ T cells, CD4+ T cells, macrophages, neutrophils, and dendritic cells). B cells and MED8 were independent predictors of overall survival. We found significant correlations between the somatic copy number alterations of the MED6, MED8, MED10, MED20, MED21, MED22, MED24, and MED25 molecules and the abundance of immune infiltrates.

CONCLUSIONS

Our study delineated a thorough landscape to investigate the therapeutic and prognostic potentials of the MED family for HCC cases, which yielded promising results for the development of immunotherapeutic drugs and construction of a prognostic stratification model.

RNA-Binding Protein Expression Alters Upon Differentiation of Human B Cells and T Cells

B cells and T cells are key players in the defence against infections and malignancies. To exert their function, B cells and T cells differentiate into effector and memory cells. Tight regulation of these differentiation processes is key to prevent their malfunction, which can result in life-threatening disease. Lymphocyte differentiation relies on the appropriate timing and dosage of regulatory molecules, and post-transcriptional gene regulation (PTR) is a key player herein. PTR includes the regulation through RNA-binding proteins (RBPs), which control the fate of RNA and its translation into proteins. To date, a comprehensive overview of the RBP expression throughout lymphocyte differentiation is lacking. Using transcriptome and proteome analyses, we here catalogued the RBP expression for human B cells and T cells. We observed that even though the overall RBP expression is conserved, the relative RBP expression is distinct between B cells and T cells. Differentiation into effector and memory cells alters the RBP expression, resulting into preferential expression of different classes of RBPs. For instance, whereas naive T cells express high levels of translation-regulating RBPs, effector T cells preferentially express RBPs that modulate mRNA stability. Lastly, we found that cytotoxic CD8+ and CD4+ T cells express a common RBP repertoire. Combined, our study reveals a cell type-specific and differentiation-dependent RBP expression landscape in human lymphocytes, which will help unravel the role of RBPs in lymphocyte function.

Cytokines and chemokines expression pattern in herpes simplex virus type-1 encephalitis

To explore the differently expressed cytokines and chemokines to understand the pathways that lead to herpes simplex encephalitis (HSE). Mice in the experimental group were inoculated intracranially with HSV-1. A high-throughput cytokine chip assay was employed to assess the expression of cytokines/chemokines in the mice brain. GO, KEGG, and PPIs analyses were used to investigate the biological process (BP), pathways and interaction network of the differently expressed proteins (DEPs) in HSE. 13 DEPs and various proteins-related signal pathways were identified in HSE, including three new factors (IL-1α, MIP-1γ, and sTNF RI). The proteins were mainly implicated in leukocyte activation and chemotaxis. Additionally, the DEPs constituted a pivotal protein interaction network where IL-6 might be a mediator. 13 DEPs and a series of related signal pathways were associated with the pathophysiological mechanisms responsible for HSE. IL-6 might be a key mediator in the inflammatory responses to the disease.

Novel DNA methylation markers of PRRSV-specific antibodies and their intergenerational transmission from pregnant sows to piglets

The main strategy for preventing porcine reproductive and respiratory syndrome (PRRS) is vaccination. However, current commercial porcine reproductive and respiratory syndrome virus (PRRSV) vaccines have limited effectiveness and may even cause infections in pigs. The identification of stable molecular markers associated with immune responses to PRRSV vaccination in pigs provides a new approach for PRRS prevention. DNA methylation, the most stable epigenetic molecular marker related to PRRSV vaccination, has not been investigated. In the current research, we used whole genome bisulfite sequencing (WGBS) to investigate DNA methylation in pregnant sows that received PRRSV vaccination and their piglets with high and low PRRSV-specific antibody levels. By performing methylation data analysis and basing on our previous transcriptomic studies, we identified several differentially methylated genes (DMGs) that are involved in the pathways of inflammatory and immune responses. Among the DMGs, ISG15, MX1, SERPINE1, GNG11 and IFIT3 were common hub genes in the two generations. MX1 and GNG11 were located in quantitative trait loci related with PRRSV antibody titer and PRRSV susceptibility, respectively. These results suggest that PRRSV vaccination in sows induces DNA methylation changes in genes and DNA methylation changes occur through intergenerational transmission. The novel DNA methylation markers and target genes observed in our study provide new insights into the molecular mechanisms of immune responses to PRRSV vaccination across two pig generations.

Preliminary therapeutic and mechanistic evaluation of S-allylmercapto-N-acetylcysteine in the treatment of pulmonary emphysema

The objective of this work was to study the effects and mechanisms of S-allylmercapto-N-acetylcysteine (ASSNAC) in the treatment of pulmonary emphysema based on network pharmacology analysis and other techniques. Firstly, the potential targets associated with ASSNAC and COPD were integrated using public databases. Then, a protein-protein interaction network was constructed using String database and Cytoscape software. The Gene Ontology analysis and Kyoto Encyclopedia of Genes and Genomes pathway analysis were performed on DAVID platform. The molecular docking of ASSNAC with some key disease targets was implemented on the SwissDock platform. To verify the results of the network pharmacology, a pulmonary emphysema mice model was established and treated with ASSNAC. Besides, the expressions of the predicted targets were detected by immunohistochemistry, Western blot analysis or enzyme-linked immunosorbent assay. Results showed that 33 overlapping targets are achieved, including CXCL8, ICAM1, MAP2K1, PTGS2, ACE and so on. The critical pathways of ASSNAC against COPD involved arachidonic acid metabolism, chemokine pathway, MAPK pathway, renin-angiotensin system, and others. Pharmacodynamic experiments demonstrated that ASSNAC decreased the pulmonary emphysema and inflammation in the pulmonary emphysema mice. Therefore, these results confirm the perspective of network pharmacology in the target verification, and indicate the treatment potential of ASSNAC against COPD.

Prognostic value and immune infiltration of a novel stromal/immune score-related P2RY12 in lung adenocarcinoma microenvironment

OBJECTIVE

Increasing evidence highlights the clinical implications of P2RY12 that belongs to the family of G-protein coupled receptors in carcinogenesis. Here, this study was designed to explore the associations between P2RY12 and tumor immune microenvironment of lung adenocarcinoma (LUAD).

METHODS

Based on 352 LUAD samples from The Cancer Genome Atlas (TCGA), stromal and immune scores of each sample were estimated by ESTIMATE algorithm. Differential expression analysis was presented between stromal/immune high- and low-score groups. Protein-protein interaction (PPI) was then constructed by STRING database. Univariate and multivariate Cox regression analysis was utilized to screen prognosis-related factors. Co-expressed genes of P2RY12 were analyzed, followed by functional enrichment analysis. Furthermore, the correlation between P2RY12 and immune cell infiltrations was estimated using the TIMER database. P2RY12 expression was validated between 37 pairs of LUAD and normal tissues using RT-qPCR and immunohistochemistry. After overexpressing P2RY12, the proliferation and migration of A549 cells was detected by CCK-8 and scratch test.

RESULTS

145 up- and 102 down-regulated stromal- and immune-related mRNAs were identified for LUAD. Based on 145 up-regulated mRNAs, a PPI network was conducted, consisting of 95 nodes and 210 relationship pairs. Ten hub genes were then identified. Among them, CCR8, CNR2, CXCR5, GPR18, GPR31 and P2RY12 were in association with overall survival of patients with LUAD. After adjusting other variables, P2RY12 expression was an independent prognostic factor for LUAD. 288 co-expressed genes of P2RY12 were determined and these co-expressed genes were primarily involved in immune-related biological processes or pathways. Moreover, P2RY12 was significantly correlated with M2 macrophage and dendritic cell infiltration. After validation, P2RY12 expression was significantly decreased in LUAD than normal tissues. Its overexpression distinctly suppressed proliferation and migration of A549 cells.

CONCLUSION

Our findings suggest that P2RY12 is an immune infiltration-related prognostic marker for LUAD.

Proteomics of mucosal exosomes of Cynoglossus semilaevis altered when infected by Vibrio harveyi

The cargo of exosomes contains proteins with various functions, which might be promising biomarkers for disease diagnosis and prognosis. To explore the impact of the Vibrio harveyi pathogen on Cynoglossus semilaevis from a different perspective and develop promising biomarkers for infection, the exosomes from epidermal mucus of healthy controls（EC）and sick fish（ES）were extracted and identified, coupled with proteomic screening through iTRAQ followed with LC-MS/MS. 1531 credible proteins were obtained relating to structural, metabolic and immunological functions. 359 different expressed proteins (DEPs) (FC > 2 or FC < 0.5) were found, with 161 up-regulated and 198 down-regulated in ES. Based on the database of C. semilaevis on Uniprot, 71 proteins were characterized as concrete names, including 19 up-regulated proteins and 52 down-regulated proteins, and were selected as subjects for further studies. Ferritin, Toll-like receptor 5S protein and Calcium-transporting ATPase were upregulated, while Histone H2B and Eukaryotic translation initiation factor 5A were downregulated, consistent with the expression levels of related mRNAs in skin tissue verified by qRT-PCR. The integrated analysis between miRomics and proteomics also provided possible regulatory relationships mediated by mucous exosomes during infection. The signature proteins in mucosal exosomes could make sense in the explanation of the infection defending mechanism and the development of biomarkers which can differentiate diseased and healthy C. semilaevis individuals.

Identification of prognostic genes and construction of a novel gene signature in the skin melanoma based on the tumor microenvironment

Skin melanoma remains a highly prevalent and yet deadly form of cancer, with the exact degree of melanoma-associated mortality being strongly dependent upon the local tumor microenvironment. The exact composition of stromal and immune cells within this microenvironmental region has the potential to profoundly impact melanoma progression and prognosis. As such, the present study was designed with the goal of clarifying the predictive relevance of stromal and immune cell-related genes in melanoma patients through comprehensive bioinformatics analyses. We therefore analyzed melanoma sample gene expression within The Cancer Genome Atlas database and employed the ESTIMATE algorithm as a means of calculating both stromal and immune scores that were in turn used for identifying differentially expressed genes (DEGs). Subsequently, univariate analyses were used to detect DEGs associated with melanoma patient survival, and through additional functional enrichment analyses, we determined that these survival-related DEGs are largely related to inflammatory and immune responses. A prognostic signature comprised of 10 genes (IL15, CCL8, CLIC2, SAMD9L, TLR2, HLA.DQB1, IGHV1-18, RARRES3, GBP4, APOBEC3G) was generated. This 10-gene signature effectively separated melanoma patients into low- and high-risk groups based upon their survival. These low- and high-risk groups also exhibited distinct immune statuses and differing degrees of immune cell infiltration. In conclusion, our results offer novel insights into a number of microenvironment-associated genes that impact survival outcomes in melanoma patients, potentially highlighting these genes as viable therapeutic targets.

Identification of 4-genes model in papillary renal cell tumor microenvironment based on comprehensive analysis

BACKGROUND

The tumor microenvironment acts a pivotal part in the occurrence and development of tumor. However, there are few studies on the microenvironment of papillary renal cell carcinoma (PRCC). Our study aims to explore prognostic genes related to tumor microenvironment in PRCC.

METHODS

PRCC expression profiles and clinical data were extracted from The Cancer Gene Atlas (TCGA) and Gene Expression Omnibus (GEO) database. Immune/stromal scores were performed utilizing the ESTIMATE algorithm. Three hundred fifty-seven samples were split into two groups on the basis of median immune/stromal score, and comparison of gene expression was conducted. Intersect genes were obtained by Venn diagrams. Hub genes were selected through protein-protein interaction (PPI) network construction, and relevant functional analysis was conducted by DAVID. We used Kaplan-Meier analysis to identify the correlations between genes and overall survival (OS) and progression-free survival (PFS). Univariate and multivariate cox regression analysis were employed to construct survival model. Cibersort was used to predict the immune cell composition of high and low risk group. Combined nomograms were built to predict PRCC prognosis. Immune properties of PRCC were validated by The Cancer Immunome Atlas (TCIA).

RESULTS

We found immune/stromal score was correlated with T pathological stages and PRCC subtypes. Nine hundred eighty-nine differentially expressed genes (DEGs) and 1169 DEGs were identified respectively on the basis of immune and stromal score. Venn diagrams indicated that 763 co-upregulated genes and 4 co-downregulated genes were identified. Kaplan-Meier analysis revealed that 120 genes were involved in tumor prognosis. Then PPI network analysis identified 22 hub genes, and four of which were significantly related to OS in patients with PRCC confirmed by cox regression analysis. Finally, we constructed a prognostic nomogram which combined with influence factors.

CONCLUSIONS

Four tumor microenvironment-related genes (CD79A, CXCL13, IL6 and CCL19) were identified as biomarkers for PRCC prognosis.

Identification of key genes and pathways in discoid lupus skin via bioinformatics analysis

ABSTRACT

Discoid lupus erythematosus (DLE) is the most common skin manifestation of lupus; however, the molecular mechanisms underlying DLE remain unknown. Therefore, we aimed to identify key differentially expressed genes (DEGs) in discoid lupus skin and investigate their potential pathways.To identify candidate genes involved in the occurrence and development of the disease, we downloaded the microarray datasets GSE52471 and GSE72535 from the Gene Expression Database (GEO). DEGs between discoid lupus skin and normal controls were selected using the GEO2R tool and Venn diagram software (http://bioinformatics.psb.ugent.be/webtools/Venn/). The Database for Annotation, Visualization, and Integrated Discovery (DAVID), Enrichr, and Cytoscape ClueGo were used to analyze the Kyoto Encyclopedia of Gene and Genome pathways and gene ontology. Protein-protein interactions (PPIs) of these DEGs were further assessed using the Search Tool for the Retrieval Interacting Genes version 10.0.Seventy three DEGs were co-expressed in both datasets. DEGs were predominantly upregulated in receptor signaling pathways of the immune response. In the PPI network, 69 upregulated genes were selected. Furthermore, 4 genes (CXCL10, ISG15, IFIH1, and IRF7) were found to be significantly upregulated in the RIG-I-like receptor signaling pathway, from analysis of Enrichr and Cytoscape ClueGo.The results of this study may provide new insights into the potential molecular mechanisms of DLE. However, further experimentation is required to confirm these findings.

The identification of neutrophils-mediated mechanisms and potential therapeutic targets for the management of sepsis-induced acute immunosuppression using bioinformatics

Neutrophils have crucial roles in defensing against infection and adaptive immune responses. This study aimed to investigate the genetic mechanism in neutrophils in response to sepsis-induced immunosuppression.The GSE64457 dataset was downloaded from the Gene Expression Omnibus database and the neutrophil samples (D3-4 and D6-8 post sepsis shock) were assigned into two groups. The differentially expressed genes (DEGs) were identified. The Short Time-series Expression Miner (STEM) clustering analysis was conducted to select the consistently changed DEGs post sepsis shock. The overlapping genes between the DEGs and the deposited genes associated with immune, sepsis, and immunosuppression in the AmiGO2 and Comparative Toxicogenomics Database were screened out and used for the construction of the protein-protein interaction (PPI) network. The expression of several hub genes in sepsis patients was validated using the PCR analysis. The drugs targeting the hub genes and the therapy strategies for sepsis or immunosuppression were reviewed and used to construct the drug-gene-therapy-cell network to illustrate the potential therapeutic roles of the hub genes.A total of 357 overlapping DEGs between the two groups were identified and were used for the STEM clustering analysis, which generated four significant profiles with 195 upregulated (including annexin A1, ANXA1; matrix metallopeptidase 9, MMP9; and interleukin 15, IL-15) and 151 downregulated DEGs (including, AKT1, IFN-related genes, and HLA antigen genes). Then, a total of 34 of the 151 downregulated DEGs and 39 of the 195 upregulated DEGs were shared between the databases and above DEGs, respectively. The PPI network analysis identified a downregulated module including IFN-related genes. The deregulation of DEGs including AKT1 (down), IFN-inducible protein 6 (IFI6, down), IL-15 (up), and ANXA1 (up) was verified in the neutrophils from patients with sepsis-induced immunosuppression as compared with controls. Literature review focusing on the therapy showed that the upregulation of IL-15, IFN, and HLA antigens are the management targets. Besides, the AKT1 gene was targeted by gemcitabine.These findings provided additional clues for understanding the mechanisms of sepsis-induced immunosuppression. The drugs targeting AKT1 might provide now clues for the management strategy of immunosuppression with the intention to prevent neutrophil infiltration.

Immunological Hallmarks of Inflammatory Status in Vaso-Occlusive Crisis of Sickle Cell Anemia Patients

Sickle Cell Anemia (SCA) is the most common genetic disorder around the world. The mutation in the β-globin gene is responsible for a higher hemolysis rate, with further involvement of immunological molecules, especially cytokines, chemokines, growth factors, and anaphylatoxins. These molecules are responsible for inducing and attracting immune cells into circulation, thus contributing to increases in leukocytes and other pro-inflammatory mediators, and can culminate in a vaso-occlusive crisis (VOC). This study aimed to characterize the levels of these molecules in SCA patients in different clinical conditions in order to identify potential hallmarks of inflammation in these patients. An analytical prospective study was conducted using the serum of SCA patients in steady-state (StSt; n = 27) and VOC (n = 22), along with 53 healthy donors (HD). Samples from the VOC group were obtained on admission and on discharge, in the convalescent phase (CV). Levels of chemokines (CXCL8, CXCL10, CL2, CLL3, CCL4, CL5, and CCL11), cytokines (IL-1β, IL-1ra, IL-2, IL-4, IL-5, IL-6, IL-7, IL-10, IL-12p70, IL-13, IL-17A, TNF-α, and IFN-γ) and growth factors (VEGF, FGFb, PDGF-BB, GM-CSF, and G-CSF) were measured using a Luminex assay, and anaphylatoxins (C3a, C4a, and C5a) were measured using Cytometric Bead Array. SCA patients in StSt showed a pro-inflammatory profile, and were indicated as being higher producers of CCL2, IL-1β, IL-12p70, IFN-γ, IL-17A, and GM-CSF, while VOC is highlighted by molecules IL-4 and IL-5, but also IL-2, IL-7, PDGF-BB, and G-CSF. PDGF-BB and IL-1ra seemed to be two important hallmarks for the acute-to-chronic stage, due to their significant decrease after crisis inflammation and statistical difference in VOC and CV groups. These molecules show higher levels and a strong correlation with other molecules in VOC. Furthermore, they remain at higher levels even after crisis recovery, which suggest their importance in the role of inflammation during crisis and participation in immune cell adhesion and activation. These results support a relevant role of cytokines, neutrophil and monocytes, since these may act as markers of VOC inflammation in SCA patients.

Identification of Potential Risk Genes and the Immune Landscape of Idiopathic Pulmonary Arterial Hypertension via Microarray Gene Expression Dataset Reanalysis

Gene dysfunction and immune cell infiltration play an essential role in the pathogenesis of idiopathic pulmonary arterial hypertension (IPAH). We aimed to investigate the immune landscape and novel differentially expressed genes (DEGs) of IPAH. In addition, potential druggable molecular targets for IPAH were also explored. In this study, the GSE117261 dataset was reanalyzed to explore the immune landscape and hub DEGs of IPAH. Lasso Cox regression analysis and receiver operating characteristic curve analysis were performed to detect the predictive value of IPAH. Additionally, the underlying drug targets for IPAH treatment were determined by drug-gene analysis. IPAH was significantly associated with the transforming growth factor-β (TGF-β) signaling pathway and Wnt signaling pathway as well as energetic metabolism dysfunction. We identified 31 upregulated and 39 downregulated DEGs in IPAH patients. Six hub genes, namely, SAA1, CCL5, CXCR1, CXCR2, CCR1, and ADORA3, were related to IPAH pathogenesis regardless of sex differences. Prediction model analysis showed that the area under the curve values of the hub DEGs except CXCR2 were all above 0.9 for distinguishing IPAH patients. In addition, the relative proportions of 5 subtypes of immune cells, namely, CD₈ ⁺ T cells, CD₄ ⁺ memory resting T cells, γ delta T cells, M1 macrophages, and resting mast cells, were significantly upregulated in the IPAH samples, while 6 subtypes of immune cells, namely, CD₄ ⁺ naive T cells, resting NK cells, monocytes, M0 macrophages, activated mast cells, and neutrophils, were downregulated. Additionally, a total of 17 intersecting drugs targeting 5 genes, CCL5, CXCR1, CXCR2, CCR1, and ADORA3, were generated as potential druggable molecular targets for IPAH. Our study revealed the underlying correlations between genes and immune cells in IPAH and demonstrated for the first time that SAA1, CCL5, CXCR1, CCR1, and ADORA3 may be novel genetic targets for IPAH.

A live attenuated Edwardsiella tarda vaccine induces immunological expression pattern in Japanese flounder (Paralichthys olivaceus) in the early phase of immunization

A previous study showed that an attenuated Edwardsiella tarda strain, TXhfq, as a live vaccine could elicit protective immune effects in fish against E. tarda infection. In the current study, in order to clarify the molecular mechanism of fish immune response at the early stage after TXhfq vaccination, RNA-Seq technology was used to compare the transcriptomes of skin, intestine, and spleen between bath-vaccinated and unvaccinated Japanese flounder (Paralichthys olivaceus). An average of 46.6 million clean reads per library was obtained, ~88.04% of which were successfully mapped to the reference genome, and approximately 24,600 genes were detected in each sample. A total of 565, 878, and 1258 differential expression genes (DEGs) were found in skin, intestine, and spleen, respectively, including 1263 up-regulated genes and 1438 down-regulated genes. The DEGs exhibited different characteristics in each tissue. One hundred and sixteen DEGs belonging to six immune related categories were scrutinized, i.e., inflammatory factors, cytokines, complement and coagulation system, mucins, phagocytosis, and antigen processing and presentation. A protein-protein interaction network was constructed to get the interaction network between immune genes during the early stage of immunization. The top six hub genes highly regulated by TXhfq formed complicated interaction relationship with each other, which were involved in immune processes, notably inflammation and phagocytosis. Our results provide valuable information for the understanding of the immune mechanism underlying the protection of live attenuated vaccines in fish.

SECAT: Quantifying Protein Complex Dynamics across Cell States by Network-Centric Analysis of SEC-SWATH-MS Profiles

Protein-protein interactions (PPIs) play critical functional and regulatory roles in cellular processes. They are essential for macromolecular complex formation, which in turn constitutes the basis for protein interaction networks that determine the functional state of a cell. We and others have previously shown that chromatographic fractionation of native protein complexes in combination with bottom-up mass spectrometric analysis of consecutive fractions supports the multiplexed characterization and detection of state-specific changes of protein complexes. In this study, we extend co-fractionation and mass spectrometric data analysis to perform quantitative, network-based studies of proteome organization, via the size-exclusion chromatography algorithmic toolkit (SECAT). This framework explicitly accounts for the dynamic nature and rewiring of protein complexes across multiple cell states and samples, thus, elucidating molecular mechanisms that are differentially implemented across different experimental settings. Systematic analysis of multiple datasets shows that SECAT represents a highly scalable and effective methodology to assess condition/state-specific protein-network state. A record of this paper's transparent peer review process is included in the Supplemental Information.

Transcriptome analysis based on RNA-seq of common innate immune responses of flounder cells to IHNV, VHSV, and HIRRV

Viral hemorrhagic septicemia virus (VHSV) and hirame rhabdovirus (HIRRV) belong to the genus Novirhabdovirus and are the causative agents of a serious disease in cultured flounder. However, infectious hematopoietic necrosis virus (IHNV), a prototype of the genus Novirhabdovirus, does not cause disease in flounder. To determine whether IHNV growth is restricted in flounder cells, we compared the growth of IHNV with that of VHSV and HIRRV in hirame natural embryo (HINAE) cells infected with novirhabdoviruses at 1 multiplicity of infection. Unexpectedly, we found that IHNV grew as well as VHSV and HIRRV. For successful growth in host cells, viruses modulate innate immune responses exerted by virus-infected cells. Our results suggest that IHNV, like VHSV and HIRRV, has evolved the ability to overcome the innate immune response of flounder cells. To determine the innate immune response genes of virus-infected HINAE cells which are commonly modulated by the three novirhabdoviruses, we infected HINAE cells with novirhabdoviruses at multiplicity of infection (MOI) 1 and performed an RNA sequencing-based transcriptome analysis at 24 h post-infection. We discovered ~12,500 unigenes altered by novirhabdovirus infection and found that many of these were involved in multiple cellular pathways. After novirhabdovirus infection, 170 genes involved in the innate immune response were differentially expressed compared to uninfected cells. Among them, 9 genes changed expression by more than 2-fold and were commonly modulated by all three novirhabdoviruses. Interferon regulatory factor 8 (IRF8), C-X-C motif chemokine receptor 1 (CXCR1), Toll/interleukin-1 receptor domain-containing adapter protein (TIRAP), cholesterol 25-hydroxylase (CH25H), C-X-C motif chemokine ligand 11, duplicate 5 (CXCL11.5), and Toll-like receptor 2 (TLR2) were up-regulated, whereas C-C motif chemokine receptor 6a (CCR6a), interleukin-12a (IL12a), and Toll-like receptor 1 (TLR1) were down-regulated. These genes have been reported to be involved in antiviral responses and, thus, their modulation may be critical for the growth of novirhabdovirus in flounder cells. This is the first report to identify innate immune response genes in flounder that are commonly modulated by IHNV, VHSV, and HIRRV. These data will provide new insights into how novirhabdoviruses survive the innate immune response of flounder cells.

Network pharmacology-based study of the protective mechanism of conciliatory anti-allergic decoction on asthma

BACKGROUND

This study aimed to explore the underlying anti-asthma pharmacological mechanisms of conciliatory anti-allergic decoction (CAD) with a network pharmacology approach.

METHODS

Traditional Chinese medicine related databases were utilized to screen the active ingredients of CAD. Targets of CAD for asthma treatment were also identified based on related databases. The protein-protein interaction network, biological function and KEGG pathway enrichment analysis, and molecular docking of the targets were performed. Furthermore, an asthma mouse model experiment involving HE staining, AB-PAS staining, and ELISA was also performed to assess the anti-asthma effect of CAD.

RESULTS

There were 77 active ingredients in CAD, including quercetin, kaempferol, stigmasterol, luteolin, cryptotanshinone, beta-sitosterol, acacetin, naringenin, baicalin, and 48 related targets for asthma treatment, mainly including TNF, IL4, IL5, IL10, IL13 and IFN-γ, were identified with ideal molecular docking binding scores by network pharmacology analysis. KEGG pathway analysis revealed that these targets were directly involved in the asthma pathway, Th1 and Th2 cell differentiation, and signaling pathways correlated with asthma (NF-κB, IL17, T cell receptor, TNF, JAK-STAT signaling pathways, etc.). Animal experiments also confirmed that CAD could attenuate inflammatory cell invasion, goblet cell hyperplasia and mucus secretion. The levels of the major targets TNF-α, IL4, IL5, and IL13 can also be regulated by CAD in an asthma mouse model.

CONCLUSION

The anti-asthma mechanism of CAD possibly stemmed from the active ingredients targeting asthma-related targets, which are involved in the asthma pathway and signaling pathways to exhibit therapeutic effects.

Natural killer cell immunotypes related to COVID-19 disease severity

Understanding innate immune responses in COVID-19 is important to decipher mechanisms of host responses and interpret disease pathogenesis. Natural killer (NK) cells are innate effector lymphocytes that respond to acute viral infections but might also contribute to immunopathology. Using 28-color flow cytometry, we here reveal strong NK cell activation across distinct subsets in peripheral blood of COVID-19 patients. This pattern was mirrored in scRNA-seq signatures of NK cells in bronchoalveolar lavage from COVID-19 patients. Unsupervised high-dimensional analysis of peripheral blood NK cells furthermore identified distinct NK cell immunotypes that were linked to disease severity. Hallmarks of these immunotypes were high expression of perforin, NKG2C, and Ksp37, reflecting increased presence of adaptive NK cells in circulation of patients with severe disease. Finally, arming of CD56bright NK cells was observed across COVID-19 disease states, driven by a defined protein-protein interaction network of inflammatory soluble factors. This study provides a detailed map of the NK cell activation landscape in COVID-19 disease.

Natural killer cell immunotypes related to COVID-19 disease severity

Understanding innate immune responses in COVID-19 is important to decipher mechanisms of host responses and interpret disease pathogenesis. Natural killer (NK) cells are innate effector lymphocytes that respond to acute viral infections but might also contribute to immunopathology. Using 28-color flow cytometry, we here reveal strong NK cell activation across distinct subsets in peripheral blood of COVID-19 patients. This pattern was mirrored in scRNA-seq signatures of NK cells in bronchoalveolar lavage from COVID-19 patients. Unsupervised high-dimensional analysis of peripheral blood NK cells furthermore identified distinct NK cell immunotypes that were linked to disease severity. Hallmarks of these immunotypes were high expression of perforin, NKG2C, and Ksp37, reflecting increased presence of adaptive NK cells in circulation of patients with severe disease. Finally, arming of CD56^bright NK cells was observed across COVID-19 disease states, driven by a defined protein-protein interaction network of inflammatory soluble factors. This study provides a detailed map of the NK cell activation landscape in COVID-19 disease.

Excessive Neutrophils and Neutrophil Extracellular Traps in COVID-19

Background: Cases of excessive neutrophil counts in the blood in severe coronavirus disease (COVID-19) patients have drawn significant attention. Neutrophil infiltration was also noted on the pathological findings from autopsies. It is urgent to clarify the pathogenesis of neutrophils leading to severe pneumonia in COVID-19. Methods: A retrospective analysis was performed on 55 COVID-19 patients classified as mild (n = 22), moderate (n = 25), and severe (n = 8) according to the Guidelines released by the National Health Commission of China. Trends relating leukocyte counts and lungs examined by chest CT scan were quantified by Bayesian inference. Transcriptional signatures of host immune cells of four COVID19 patients were analyzed by RNA sequencing of lung specimens and BALF. Results: Neutrophilia occurred in 6 of 8 severe patients at 7-19 days after symptom onset, coinciding with lesion progression. Increasing neutrophil counts paralleled lesion CT values (slope: 0.8 and 0.3-1.2), reflecting neutrophilia-induced lung injury in severe patients. Transcriptome analysis revealed that neutrophil activation was correlated with 17 neutrophil extracellular trap (NET)-associated genes in COVID-19 patients, which was related to innate immunity and interacted with T/NK/B cells, as supported by a protein-protein interaction network analysis. Conclusion: Excessive neutrophils and associated NETs could explain the pathogenesis of lung injury in COVID-19 pneumonia.

A robust signature associated with patient prognosis and tumor immune microenvironment based on immune-related genes in lung squamous cell carcinoma

BACKGROUND

Lung squamous cell carcinoma (LUSC) is one common type of lung cancer. Immune-related genes (IRGs) are closely associated with cancer prognosis. This study aims to screen the key genes associated with LUSC and establish an immune-related prognostic model.

METHODS

Based on the Cancer Genome Atlas (TCGA) database, we screened the differentially expressed genes (DEGs) between LUSC and normal samples. Intersecting the DEGs with the immune-related genes (IRGs), we obtained the differentially expressed IRGs (DEIRGs). Univariate as well as multivariate Cox regression analyses were performed to identify the survival-associated IRGs and establish an immune-related prognostic model. The relationship between the prognostic model and tumor-infiltrating immune cells was analyzed by TIMER and CIBERSORT.

RESULTS

A total of 229 DEIRGs were screened, and 14 IRGs associated with survival were identified using univariate Cox analysis. Among the 14 IRGs, six genes were selected out using Lasso and multivariate Cox analyses, and they were used to build the prognostic model. Further analysis indicated that overall survival (OS) of high-risk groups was lower than that of low-risk groups. High risk score was independently related to worse OS. Moreover, the risk score was positively correlated with several immune infiltration cells. Finally, the efficacy of the prognostic model was validated by another independent cohort GSE73403.

CONCLUSION

The DEIRGs described in the study may have the potential to be the prognostic molecular markers for LUSC. In addition, the risk score model could predict the OS and provides more information for the immunotherapy of patients with LUSC.

Personalized therapy design for systemic lupus erythematosus based on the analysis of protein-protein interaction networks

We analyzed protein expression data for Lupus patients, which have been obtained from publicly available databases. A combination of systems biology and statistical thermodynamics approaches was used to extract topological properties of the associated protein-protein interaction networks for each of the 291 patients whose samples were used to provide the molecular data. We have concluded that among the many proteins that appear to play critical roles in this pathology, most of them are either ribosomal proteins, ubiquitination pathway proteins or heat shock proteins. We propose some of the proteins identified in this study to be considered for drug targeting.

Prediction of triptolide targets in rheumatoid arthritis using network pharmacology and molecular docking

Network pharmacology is a novel approach that uses bioinformatics to predict and identify multiple drug targets and interactions in disease. Here, we used network pharmacology to investigate the mechanism by which triptolide acts in rheumatoid arthritis (RA). We first searched public databases for genes and proteins known to be associated with RA, as well as those predicted to be targets of triptolide, and then used Ingenuity Pathway Analysis (IPA) to identify enriched gene pathways and networks. Networks and pathways that overlapped between RA-associated proteins and triptolide target proteins were then used to predict candidate protein targets of triptolide in RA. The following proteins were found to occur in both RA-associated networks and triptolide target networks: CD274, RELA, MCL1, MAPK8, CXCL8, STAT1, STAT3, c-JUN, JNK, c-Fos, NF-κB, and TNF-α. Docking studies suggested that triptolide can fit in the binding pocket of the six top candidate triptolide target proteins (CD274, RELA, MCL1, MAPK8, CXCL8 and STAT1). The overlapping pathways were activation of Th1 and Th2 cells, macrophages, fibroblasts and endothelial cells in RA, while the overlapping networks were involved in cellular movement, hematological system development and function, immune cell trafficking, cell-to-cell signaling and interaction, inflammatory response, cellular function and maintenance, and cell death and survival. These results show that network pharmacology can be used to generate hypotheses about how triptolide exerts therapeutic effects in RA. Network pharmacology may be a useful method for characterizing multi-target drugs in complex diseases.

Understanding allergic multimorbidity within the non-eosinophilic interactome

BACKGROUND

The mechanisms explaining multimorbidity between asthma, dermatitis and rhinitis (allergic multimorbidity) are not well known. We investigated these mechanisms and their specificity in distinct cell types by means of an interactome-based analysis of expression data.

METHODS

Genes associated to the diseases were identified using data mining approaches, and their multimorbidity mechanisms in distinct cell types were characterized by means of an in silico analysis of the topology of the human interactome.

RESULTS

We characterized specific pathomechanisms for multimorbidities between asthma, dermatitis and rhinitis for distinct emergent non-eosinophilic cell types. We observed differential roles for cytokine signaling, TLR-mediated signaling and metabolic pathways for multimorbidities across distinct cell types. Furthermore, we also identified individual genes potentially associated to multimorbidity mechanisms.

CONCLUSIONS

Our results support the existence of differentiated multimorbidity mechanisms between asthma, dermatitis and rhinitis at cell type level, as well as mechanisms common to distinct cell types. These results will help understanding the biology underlying allergic multimorbidity, assisting in the design of new clinical studies.

Assessment of the genes and molecular mechanisms of B cells activation through systems biology approaches

Two classes of T helper lymphocytes, Th1 and Th2, have different roles in B cell activation based on specific cytokines. To understand the difference of molecular mechanisms of B cell activation, the microarray dataset of B cells co-cultured with type 1 and 2 T helper, Be1 and Be2, were investigated. After quality assessment, using the GEO2R tool, the GSE84948 dataset was re-analyzed. Genes with adjusted p-value ⩽ 0.05 were assumed as differentially expressed (DE). The protein-protein interaction (PPI) networks were constructed using CluePedia plugin of Cytoscape, and analyzed by NetworkAnalyzer tool and MCODE plugin. Using ClueGO plugin of Cytoscape software, gene ontology (GO) analysis was performed. The comparison of Be1 and Be2 cells with naive B cells revealed 8742 and 8748 DE genes, respectively. The topology analysis of PPI networks predicted central genes. Among these, Jak3, Actrt3, and Pik3cb genes were determined as crucial genes in Be1 network. Prkx, Smarca4, and Jak2 genes were defined as critical genes in Be2 PPI network. GO analysis with PPI networks genes resulted in the promotion of immune system activation. In conclusion, we explored holistic methods for molecular assay of the difference between B cell activation mechanisms with Th1 and Th2.

Identification of hepatitis B virus and liver cancer bridge molecules based on functional module network

BACKGROUND

The potential role of chronic inflammation in the development of cancer has been widely recognized. However, there has been little research fully and thoroughly exploring the molecular link between hepatitis B virus (HBV) and hepatocellular carcinoma (HCC).

AIM

To elucidate the molecular links between HBV and HCC through analyzing the molecular processes of HBV-HCC using a multidimensional approach.

METHODS

First, maladjusted genes shared between HBV and HCC were identified by disease-related differentially expressed genes. Second, the protein-protein interaction network based on dysfunctional genes identified a series of dysfunctional modules and significant crosstalk between modules based on the hypergeometric test. In addition, key regulators were detected by pivot analysis. Finally, targeted drugs that have regulatory effects on diseases were predicted by modular methods and drug target information.

RESULTS

The study found that 67 genes continued to increase in the HBV-HCC process. Moreover, 366 overlapping genes in the module network participated in multiple functional blocks. It could be presumed that these genes and their interactions play an important role in the relationship between inflammation and cancer. Correspondingly, significant crosstalk constructed a module level bridge for HBV-HCC molecular processes. On the other hand, a series of non-coding RNAs and transcription factors that have potential pivot regulatory effects on HBV and HCC were identified. Among them, some of the regulators also had persistent disorders in the process of HBV-HCC including microRNA-192, microRNA-215, and microRNA-874, and early growth response 2, FOS, and Kruppel-like factor 4. Therefore, the study concluded that these pivots are the key bridge molecules outside the module. Last but not least, a variety of drugs that may have some potential pharmacological or toxic side effects on HBV-induced HCC were predicted, but their mechanisms still need to be further explored.

CONCLUSION

The results suggest that the persistent inflammatory environment of HBV can be utilized as an important risk factor to induce the occurrence of HCC, which is supported by molecular evidence.

GPCR-specific autoantibody signatures are associated with physiological and pathological immune homeostasis

Autoantibodies have been associated with autoimmune diseases. However, studies have identified autoantibodies in healthy donors (HD) who do not develop autoimmune disorders. Here we provide evidence of a network of immunoglobulin G (IgG) autoantibodies targeting G protein-coupled receptors (GPCR) in HD compared to patients with systemic sclerosis, Alzheimer's disease, and ovarian cancer. Sex, age and pathological conditions affect autoantibody correlation and hierarchical clustering signatures, yet many of the correlations are shared across all groups, indicating alterations to homeostasis. Furthermore, we identify relationships between autoantibodies targeting structurally and functionally related molecules, such as vascular, neuronal or chemokine receptors. Finally, autoantibodies targeting the endothelin receptor type A (EDNRA) exhibit chemotactic activity, as demonstrated by neutrophil migration toward HD-IgG in an EDNRA-dependent manner and in the direction of IgG from EDNRA-immunized mice. Our data characterizing the in vivo signatures of anti-GPCR autoantibodies thus suggest that they are a physiological part of the immune system.

Analysis of the transcriptome data in Litopenaeus vannamei reveals the immune basis and predicts the hub regulation-genes in response to high-pH stress

Soil salinization erodes the farmlands and poses a serious threat to human life, reuse of the saline-alkali lands as cultivated resources becomes increasingly prominent. Pacific white shrimp (Litopenaeus vannamei) is an important farmed aquatic species for the development and utilization of the saline-alkali areas. However, little is known about the adaptation mechanism of this species in terms of high-pH stress. In the present study, a transcriptome analysis on the gill tissues of L. vannamei in response to high-pH stress (pH 9.3 ± 0.1) was conducted. After analyzing, the cyclic nucleotide gated channel-Ca2+ (CNGC-Ca2+) and patched 1 (Ptc1) were detected as the majority annotated components in the cAMP signaling pathway (KO04024), indicating that the CNGC-Ca2+ and Ptc1 might be the candidate components for transducing and maintaining the high-pH stress signals, respectively. The immunoglobulin superfamily (IgSF), heat shock protein (HSP), glutathione s-transferase (GST), prophenoloxidase/phenoloxidase (proPO/PO), superoxide dismutase (SOD), anti-lipopolysaccharide factor (ALF) and lipoprotein were discovered as the major transcribed immune factors in response to high-pH stress. To further detect hub regulation-genes, protein-protein interaction (PPI) networks were constructed; the genes/proteins "Polymerase (RNA) II (DNA directed) polypeptide A" (POLR2A), "Histone acetyltransferase p300" (EP300) and "Heat shock 70kDa protein 8" (HSPA8) were suggested as the top three hub regulation-genes in response to acute high-pH stress; the genes/proteins "Heat shock 70kDa protein 4" (HSPA4), "FBJ murine osteosarcoma viral oncogene homolog" (FOS) and "Nucleoporin 54kDa" (NUP54) were proposed as the top three hub regulation-genes involved in adapting endurance high-pH stress; the protein-interactions of "EP300-HSPA8" and "HSPA4-NUP54" were detected as the most important biological interactions in response to the high-pH stress; and the HSP70 family genes might play essential roles in the adaptation of the high-pH stress environment in L. vannamei. These findings provide the first insight into the molecular and immune basis of L. vannamei in terms of high-pH environments, and the construction of a PPI network might improve our understanding in revealing the hub regulation-genes in response to abiotic stress in shrimp species and might be beneficial for further studies.

Defining Host Responses during Systemic Bacterial Infection through Construction of a Murine Organ Proteome Atlas

Group A Streptococcus (GAS) remains one of the top 10 deadliest human pathogens worldwide despite its sensitivity to penicillin. Although the most common GAS infection is pharyngitis (strep throat), it also causes life-threatening systemic infections. A series of complex networks between host and pathogen drive invasive infections, which have not been comprehensively mapped. Attempting to map these interactions, we examined organ-level protein dynamics using a mouse model of systemic GAS infection. We quantified over 11,000 proteins, defining organ-specific markers for all analyzed tissues. From this analysis, an atlas of dynamically regulated proteins and pathways was constructed. Through statistical methods, we narrowed organ-specific markers of infection to 34 from the defined atlas. We show these markers are trackable in blood of infected mice, and a subset has been observed in plasma samples from GAS-infected clinical patients. This proteomics-based strategy provides insight into host defense responses, establishes potentially useful targets for therapeutic intervention, and presents biomarkers for determining affected organs during bacterial infection.

Genetic Variants in Immune-Related Pathways and Breast Cancer Risk in African American Women in the AMBER Consortium

Background: Constitutional immunity shaped by exposure to endemic infectious diseases and parasitic worms in Sub-Saharan Africa may play a role in the etiology of breast cancer among African American (AA) women.Methods: A total of 149,514 gene variants in 433 genes across 45 immune pathways were analyzed in the AMBER consortium among 3,663 breast cancer cases and 4,687 controls. Gene-based pathway analyses were conducted using the adaptive rank truncated product statistic for overall breast cancer risk, and risk by estrogen receptor (ER) status. Unconditional logistic regression analysis was used to estimate ORs and 95% confidence intervals (CIs) for single variants.Results: The top pathways were Interleukin binding (P = 0.01), Biocarta TNFR2 (P = 0.005), and positive regulation of cytokine production (P = 0.024) for overall, ER+, and ER- cancers, respectively. The most significant gene was IL2RB (P = 0.001) for overall cancer, with rs228952 being the top variant identified (OR = 0.85; 95% CI, 0.79-0.92). Only BCL3 contained a significant variant for ER+ breast cancer. Variants in IL2RB, TLR6, IL8, PRKDC, and MAP3K1 were associated with ER- disease. The only genes showing heterogeneity between ER- and ER+ cancers were TRAF1, MAP3K1, and MAPK3 (P ≤ 0.02). We also noted genes associated with autoimmune and atopic disorders.Conclusions: Findings from this study suggest that genetic variants in immune pathways are relevant to breast cancer susceptibility among AA women, both for ER+ and ER- breast cancers.Impact: Results from this study extend our understanding of how inherited genetic variation in immune pathways is relevant to breast cancer susceptibility. Cancer Epidemiol Biomarkers Prev; 27(3); 321-30. ©2018 AACR.

Network ethnopharmacological evaluation of the immunomodulatory activity of Withania somnifera

ETHNOPHARMACOLOGICAL RELEVANCE

Withania somnifera (L.) Dunal (Ashwagandha, WS) is one of the extensively explored Ayurvedic botanicals. Several properties including immunomodulation, anti-cancer and neuro-protection of the botanical have been reported. Even though, in indigenous medicine, WS is well known for its immunomodulatory activity, the molecular mechanism of immunomodulation has not been elucidated.

AIM OF THE STUDY

This study aimed the evaluation of the immunomodulatory effect of WS using network ethnopharmacology technique to elucidate the in silico molecular mechanism.

MATERIALS AND METHODS

Databases- DPED, UNPD, PubChem, Binding DB, ChEMBL, KEGG and STRING were used to gather information to develop the networks. The networks were constructed using Cytoscape 3.2.1. Data analysis was performed with the help of Excel pivot table and Cytoscape network analyzer tool.

RESULTS

Investigation for WS immune modulation mechanism identified five bioactives that are capable of regulating 15 immune system pathways through 16 target proteins by bioactive-target and protein-protein interactions. The study also unveils the potential of withanolide-phytosterol combination to achieve effective immunomodulation and seven novel bioactive-immune target combinations.

CONCLUSION

The study elucidated an in silico molecular mechanism of immunomodulation of WS. It unveils the potential of withanolide-phytosterol combination to achieve a better immunomodulation. Experimental validation of the network findings would aid in understanding the rationale behind WS immunomodulation as well as aid in bioactive formulation based drug discovery.

Global gene expression profile of peripheral blood mononuclear cells challenged with Theileria annulata in crossbred and indigenous cattle

Bovine tropical theileriosis is an important haemoprotozoan disease associated with high rates of morbidity and mortality particularly in exotic and crossbred cattle. It is one of the major constraints of the livestock development programmes in India and Southeast Asia. Indigenous cattle (Bos indicus) are reported to be comparatively less affected than exotic and crossbred cattle. However, genetic basis of resistance to tropical theileriosis in indigenous cattle is not well documented. Recent studies incited an idea that differentially expressed genes in exotic and indigenous cattle play significant role in breed specific resistance to tropical theileriosis. The present study was designed to determine the global gene expression profile in peripheral blood mononuclear cells derived from indigenous (Tharparkar) and cross-bred cattle following in vitro infection of T. annulata (Parbhani strain). Two separate microarray experiments were carried out each for cross-bred and Tharparkar cattle. The cross-bred cattle showed 1082 differentially expressed genes (DEGs). Out of total DEGs, 597 genes were down-regulated and 485 were up-regulated. Their fold change varied from 2283.93 to -4816.02. Tharparkar cattle showed 875 differentially expressed genes including 451 down-regulated and 424 up-regulated. The fold change varied from 94.93 to -19.20. A subset of genes was validated by qRT-PCR and results were correlated well with microarray data indicating that microarray results provided an accurate report of transcript level. Functional annotation study of DEGs confirmed their involvement in various pathways including response to oxidative stress, immune system regulation, cell proliferation, cytoskeletal changes, kinases activity and apoptosis. Gene network analysis of these DEGs plays an important role to understand the interaction among genes. It is therefore, hypothesized that the different susceptibility to tropical theileriosis exhibited by indigenous and crossbred cattle is due to breed-specific differences in the dealing of infected cells with other immune cells, which ultimately influence the immune response responded against T. annulata infection.

The Relationship between B-cell Epitope and Mimotope Sequences

B-cell epitope is a group of residues which is on the surface of an antigen. It invokes humoral responses. Locating B-cell epitope is important for effective vaccine design, and the development of diagnostic reagents. Mimotope-based B-cell epitope prediction method is a kind of conformational B-cell epitope prediction, and the core idea of the method is mapping the mimotope sequences which are obtained from a random phage display library. However, current mimotope-based B-cell epitope prediction methods cannot maintain a high degree of satisfaction in the circumstances of employing only mimotope sequences. In this study, we did a multi-perspective analysis on parameters for conformational B-cell epitopes and characteristics between epitope and mimotope on a benchmark datasets which contains 67 mimotope sets, corresponding to 40 unique complex structures. In these 67 cases, there are 25 antigen-antibody complexes and 42 protein-protein interactions. We analyzed the two parts separately. The results showed the mimotope sequences do have some epitope features, but there are also some epitope properties that mimotope sequences do not contain. In addition, the numbers of epitope segments with different lengths were obviously different between the antigen-antibody complexes and the protein-protein interactions. This study reflects how similar do mimotope sequence and genuine epitopes have; and evaluates existing mimotope-based B-cell epitope prediction methods from a novel viewpoint.

Targeted in vivo inhibition of specific protein-protein interactions using recombinant antibodies

With the growing availability of genomic sequence information, there is an increasing need for gene function analysis. Antibody-mediated "silencing" represents an intriguing alternative for the precise inhibition of a particular function of biomolecules. Here, we describe a method for selecting recombinant antibodies with a specific purpose in mind, which is to inhibit intrinsic protein-protein interactions in the cytosol of plant cells. Experimental procedures were designed for conveniently evaluating desired properties of recombinant antibodies in consecutive steps. Our selection method was successfully used to develop a recombinant antibody inhibiting the interaction of ARABIDOPSIS HISTIDINE PHOSPHOTRANSFER PROTEIN 3 with such of its upstream interaction partners as the receiver domain of CYTOKININ INDEPENDENT HISTIDINE KINASE 1. The specific down-regulation of the cytokinin signaling pathway in vivo demonstrates the validity of our approach. This selection method can serve as a prototype for developing unique recombinant antibodies able to interfere with virtually any biomolecule in the living cell.

In silico and in vitro studies on the protein-protein interactions between Brugia malayi immunomodulatory protein calreticulin and human C1q

Filarial parasites modulate effective immune response of their host by releasing a variety of immunomodulatory molecules, which help in the long persistence of the parasite within the host. The present study was aimed to characterize an immunomodulatory protein of Brugia malayi and its interaction with the host immune component at the structural and functional level. Our findings showed that Brugia malayi Calreticulin (BmCRT) is responsible for the prevention of classical complement pathway activation via its interaction with the first component C1q of the human host. This was confirmed by inhibition of C1q dependent lysis of immunoglobulin-sensitized Red Blood Cells (S-RBCs). This is possibly the first report which predicts CRT-C1q interaction on the structural content of proteins to explain how BmCRT inhibits this pathway. The molecular docking of BmCRT-C1q complex indicated that C1qB chain (IgG/M and CRP binding sites on C1q) played a major role in the interaction with conserved and non-conserved regions of N and P domain of BmCRT. Out of 37 amino acids of BmCRT involved in the interaction, nine amino acids (Pro(126), Glu(132), His(147), Arg(151), His(153), Met(154), Lys(156), Ala(196) and Lys(212)) are absent in human CRT. Both ELISA and in silico analysis showed the significant role of Ca(+2) in BmCRT-HuC1q complex formation and deactivation of C1r2-C1s2. Molecular dynamics studies of BmCRT-HuC1q complex showed a deviation from ∼ 0.4 nm to ∼ 1.0 nm. CD analyses indicated that BmCRT is composed of 49.6% α helix, 9.6% β sheet and 43.6% random coil. These findings provided valuable information on the architecture and chemistry of BmCRT-C1q interaction and supported the hypothesis that BmCRT binds with huC1q at their targets (IgG/M, CRP) binding sites. This interaction enables the parasite to interfere with the initial stage of host complement activation, which might be helpful in parasites establishment. These results might be utilized for help in blocking the C1q/CRT interaction and preventing parasite infection.

Proteomic comparison of 3D and 2D glioma models reveals increased HLA-E expression in 3D models is associated with resistance to NK cell-mediated cytotoxicity

Three-dimensional cell culture techniques can better reflect the in vivo characteristics of tumor cells compared with traditional monolayer cultures. Compared with their 2D counterparts, 3D-cultured tumor cells showed enhanced resistance to the cytotoxic T cell-mediated immune response. However, it remains unclear whether 3D-cultured tumor cells have an enhanced resistance to NK cell cytotoxicity. In this study, a total of 363 differentially expressed proteins were identified between the 2D- and 3D-cultured U251 cells by comparative proteomics, and an immune-associated protein-protein interaction (PPI) network based on these differential proteins was constructed by bioinformatics. Within the network, HLA-E, as a molecule for inhibiting NK cell activation, was significantly up-regulated in the 3D-cultured tumor cells. Then, we found that the 3D-cultured U251 cells exhibited potent resistance to NK cell cytotoxicity in vitro and were prone to tumor formation in vivo. The resistance of the 3D-cultured tumor cells to NK cell lysis was mediated by the HLA-E/NKG2A interaction because the administration of antibodies that block either HLA-E or NKG2A completely eliminated this resistance and significantly decreased tumor formation. Taken together, our findings indicate that HLA-E up-regulation in 3D-cultured cells may result in enhanced tumor resistance to NK cell-mediated immune response.

The structural pathway of interleukin 1 (IL-1) initiated signaling reveals mechanisms of oncogenic mutations and SNPs in inflammation and cancer

Interleukin-1 (IL-1) is a large cytokine family closely related to innate immunity and inflammation. IL-1 proteins are key players in signaling pathways such as apoptosis, TLR, MAPK, NLR and NF-κB. The IL-1 pathway is also associated with cancer, and chronic inflammation increases the risk of tumor development via oncogenic mutations. Here we illustrate that the structures of interfaces between proteins in this pathway bearing the mutations may reveal how. Proteins are frequently regulated via their interactions, which can turn them ON or OFF. We show that oncogenic mutations are significantly at or adjoining interface regions, and can abolish (or enhance) the protein-protein interaction, making the protein constitutively active (or inactive, if it is a repressor). We combine known structures of protein-protein complexes and those that we have predicted for the IL-1 pathway, and integrate them with literature information. In the reconstructed pathway there are 104 interactions between proteins whose three dimensional structures are experimentally identified; only 15 have experimentally-determined structures of the interacting complexes. By predicting the protein-protein complexes throughout the pathway via the PRISM algorithm, the structural coverage increases from 15% to 71%. In silico mutagenesis and comparison of the predicted binding energies reveal the mechanisms of how oncogenic and single nucleotide polymorphism (SNP) mutations can abrogate the interactions or increase the binding affinity of the mutant to the native partner. Computational mapping of mutations on the interface of the predicted complexes may constitute a powerful strategy to explain the mechanisms of activation/inhibition. It can also help explain how an oncogenic mutation or SNP works.

Structural pathways are important because they provide insight into signaling mechanisms; help understand the mechanism of disease-related mutations; and help in drug discovery. While extremely useful, common pathway diagrams lacking structural data are unable to provide mechanistic insight to explain oncogenic mutations or SNPs. Here we focus on the construction of the IL-1 structural pathway and map oncogenic mutations and SNPs to complexes in this pathway. Our results indicate that computational modeling of protein-protein interactions on a large scale can provide accurate, structural atom-level detail of signaling pathways in the human cell and help delineate the mechanism through which a mutation leads to disease. We show that the mutations either thwart the interactions, activating the proteins even in their absence or stabilize them, leading to the same uncontrolled outcome. Computational mapping of mutations on the interface of the predicted complexes may constitute an effective strategy to explain the mechanisms of mutations- constitutive activation or deactivation.

Linking transcriptional changes over time in stimulated dendritic cells to identify gene networks activated during the innate immune response

The innate immune response is primarily mediated by the Toll-like receptors functioning through the MyD88-dependent and TRIF-dependent pathways. Despite being widely studied, it is not yet completely understood and systems-level analyses have been lacking. In this study, we identified a high-probability network of genes activated during the innate immune response using a novel approach to analyze time-course gene expression profiles of activated immune cells in combination with a large gene regulatory and protein-protein interaction network. We classified the immune response into three consecutive time-dependent stages and identified the most probable paths between genes showing a significant change in expression at each stage. The resultant network contained several novel and known regulators of the innate immune response, many of which did not show any observable change in expression at the sampled time points. The response network shows the dominance of genes from specific functional classes during different stages of the immune response. It also suggests a role for the protein phosphatase 2a catalytic subunit α in the regulation of the immunoproteasome during the late phase of the response. In order to clarify the differences between the MyD88-dependent and TRIF-dependent pathways in the innate immune response, time-course gene expression profiles from MyD88-knockout and TRIF-knockout dendritic cells were analyzed. Their response networks suggest the dominance of the MyD88-dependent pathway in the innate immune response, and an association of the circadian regulators and immunoproteasomal degradation with the TRIF-dependent pathway. The response network presented here provides the most probable associations between genes expressed in the early and the late phases of the innate immune response, while taking into account the intermediate regulators. We propose that the method described here can also be used in the identification of time-dependent gene sub-networks in other biological systems.

Epitope mapping and in silico characterization of interactions between Der p 7 allergen and MoAb WH9

Der p 7 is an important house dust mite allergen. However, antigenic determinants of Der p 7 are largely unknown. The purpose of this study is to analyze the determinants of Der p 7 and determine the structural basis of interactions between Der p 7 and WH9, an IgE-binding inhibition mouse monoclonal antibody (MoAb). IgE and WH9-reactive determinant(s) was identified by immunoblot using allergen mutants. A 3-D binary complex structure of Der p 7 and WH9 was simulated with homology modeling and docking methods. Our results obtained showed that among the five Der p 7 mutants (S156A, I157A, L158A, D159A, P160A), serum no. 1045 with IgE-binding against Der p 7 exhibited a reduced IgE immunoblot reactivity against Der p 7 L158A and D159A mutants. WH9 showed reduced immunoblot reactivity against S156A, L158A, D159A and P160A and the observation was confirmed by immunoblot inhibition. The WH9-binding determinant on Der p 7 containing S156, L158, D159 and P160 assumes a loop-like structure. The structural model of the Der p 7-WH9 complex suggests residues S156, I157, L158, D159 and P160 of Der p 7 contribute to WH9 binding via potential hydrogen bonds, electrostatic and hydrophobic interactions. In conclusion, MoAb WH9 interacts with critical residues L158 and D159 of Der p 7 and inhibits IgE-binding to Der p 7. Results obtained advance our understanding on molecular and structural bases of the antigenicity of Der p 7, its interactions with MoAb WH9 and facilitate the design of safer immunotherapy of human atopic disorders.

Autoimmunoreactive IgGs from patients with postural orthostatic tachycardia syndrome

PURPOSE

Autoantibodies are implicated in the pathogenesis of cardiovascular diseases and cardiac arrhythmias. In this pilot study, we tested the hypothesis that autoantibodies are present in patients with postural orthostatic tachycardia syndrome (POTS).

EXPERIMENTAL DESIGN

Seven control subjects (6 F:1 M, average age 36.1 years) and ten patients with the diagnosis of POTS (7 F: 3 M, average age 35.1 years) provided informed consent and 30 mL of venous blood. Human heart membrane proteins were resolved by 2DE and immunoblotted against purified IgGs from controls and patients.

RESULTS

Eighteen protein spots immunoreactive specifically against patient IgGs were detected and they were excised from gels, trypsin-digested, and analyzed by nanoLC-electrospray MS/MS. Forty unique proteins were identified and these include proteins that are associated with cardiac hypertrophy (mimecan, myozenin), cardiac remodeling (periostin), cardiomyopathy (desmin, desmoplakin), cell survival (laminin), structural integrity (filamin), chaperone proteins (crystalline, HSP70), mitochondrial enzymes, and channel proteins. Ingenuity Pathway Analysis showed multiple pathways were involved including those that regulate energy metabolism, redox, fibrosis, cardiac hypertrophy, and degeneration.

CONCLUSIONS AND CLINICAL RELEVANCE

Autoantibodies are present in patients with POTS. These autoantibodies cross-react with a wide range of cardiac proteins and may induce alterations in cardiac function. Autoimmune pathogenetic mechanisms should be further explored in these patients.

Effector-triggered immunity signaling: from gene-for-gene pathways to protein-protein interaction networks

In its simplicity and testability, Flor's gene-for-gene hypothesis has been a powerful driver in plant immunity research for decades. Once the molecular underpinnings of gene-for-gene resistance had come into sharper focus, there was a reassessment of Flor's hypothesis and a name change to effector-triggered immunity. As implied by the name change and exemplified by pioneering studies, plant immunity is increasingly described in terms of protein rather than genetic interactions. This progress leads to a reinterpretation of old concepts of pathogen recognition and resistance signaling and, of course, opens up new questions. Here, we provide a brief historical overview of resistance gene function and how a new focus on protein interactions can lead to a deeper understanding of the logic of plant innate immunity signaling.