Contextual Hub Analysis Tool (CHAT): A Cytoscape app for identifying contextually relevant hubs in biological networks

Molecular insights into the proteomic composition of porcine treated dentin matrix

Background

Human-treated dentin matrix (hTDM) has recently been studied as a natural extracellular matrix-based biomaterial for dentin pulp regeneration. However, porcine-treated dentin matrix (pTDM) is a potential alternative scaffold due to limited availability. However, there is a dearth of information regarding the protein composition and underlying molecular mechanisms of pTDM.Methods: hTDM and pTDM were fabricated using human and porcine teeth, respectively, and their morphological characteristics were examined using scanning electron microscopy. Stem cells derived from human exfoliated deciduous teeth (SHEDs) were isolated and characterized using flow cytometry and multilineage differentiation assays. SHEDs were cultured in three-dimensional environments with hTDM, pTDM, or biphasic hydroxyapatite/tricalcium phosphate. The expression of odontogenesis markers in SHEDs were assessed using real-time polymerase chain reaction and immunochemical staining. Subsequently, SHEDs/TDM and SHEDs/HA/TCP complexes were transplanted subcutaneously into nude mice. The protein composition of pTDM was analyzed using proteomics and compared to previously published data on hTDM.Results: pTDM and hTDM elicited comparable upregulation of odontogenesis-related genes and proteins in SHEDs. Furthermore, both demonstrated the capacity to stimulate root-related tissue regeneration in vivo. Proteomic analysis revealed the presence of 278 protein groups in pTDM, with collagens being the most abundant. Additionally, pTDM and hTDM shared 58 identical proteins, which may contribute to their similar abilities to induce odontogenesis.

Conclusions

Both hTDM and pTDM exhibit comparable capabilities in inducing odontogenesis, potentially owing to their distinctive bioactive molecular networks.

Identification of Proteomic Biomarkers in Proliferative Verrucous Leukoplakia through Liquid Chromatography With Tandem Mass Spectrometry

Proliferative verrucous leukoplakia (PVL) is an oral potentially malignant disorder associated with high risk of malignant transformation. Currently, there is no treatment available, and restrictive follow-up of patients is crucial for a better prognosis. Oral leukoplakia (OL) shares some clinical and microscopic features with PVL but exhibits different clinical manifestations and a lower rate of malignant transformation. This study aimed to investigate the proteomic profile of PVL in tissue and saliva samples to identify potential diagnostic biomarkers with therapeutic implications. Tissue and saliva samples obtained from patients with PVL were compared with those from patients with oral OL and controls. Label-free liquid chromatography with tandem mass spectrometry was employed, followed by qualitative and quantitative analyses, to identify differentially expressed proteins. Potential biomarkers were identified and further validated using immunohistochemistry. Staining intensity scan analyses were performed on tissue samples from patients with PVL, patients with OL, and controls from Brazil, Spain, and Finland. The study revealed differences in the immune system, cell cycle, DNA regulation, apoptosis pathways, and the whole proteome of PVL samples. In addition, liquid chromatography with tandem mass spectrometry analyses showed that calreticulin (CALR), receptor of activated protein C kinase 1 (RACK1), and 14-3-3 Tau-protein (YWHAQ) were highly expressed in PVL samples. Immunohistochemistry validation confirmed increased CARL expression in PVL compared with OL. Conversely, RACK1 and YWHA were highly expressed in oral potentially malignant disorder compared to the control group. Furthermore, significant differences in CALR and RACK1 expression were observed in the OL group when comparing samples with and without oral epithelial dysplasia, unlike the PVL. This research provides insights into the molecular mechanisms underlying these conditions and highlights potential targets for future diagnostic and therapeutic approaches.

IL-35 enhances angiogenic effects of small extracellular vesicles in breast cancer

As an indispensable process for breast cancer metastasis, tumour angiogenesis requires a tight interaction between cancer cells and endothelial cells in tumour microenvironment. Here, we explored the participation of small extracellular vesicles (sEVs) derived from breast cancer cells in modulating angiogenesis and investigated the effect of IL-35 in facilitating this process. Firstly, we characterized breast cancer cells-derived sEVs untreated or treated with IL-35 and visualized the internalization of these sEVs by human umbilical vein endothelial cells (HUVECs). Breast cancer cells-derived sEVs promoted endothelial cell proliferation through facilitating cell cycle progression and enhanced capillary-like structures formation and microvessel formation. Subsequent results proved that IL-35 further reinforced the angiogenic effect induced by breast cancer cells-derived sEVs. Moreover, sEVs from breast cancer cells significantly enhanced tumour growth and microvessel density in breast tumour-bearing mice model. Microarray analysis showed that IL-35 might alter the mRNA profiles of sEVs and activate the Ras/Raf/MEK/ERK signalling pathway. These findings demonstrated that IL-35 indirectly promoted angiogenesis in breast cancer through regulating the content of breast cancer cells-derived sEVs, which could be internalized by HUVECs.

Genome profiles of pathologist-defined cell clusters by multiregional LCM and G&T-seq in one triple-negative breast cancer patient

Pathological examination is the gold standard for cancer diagnosis, and breast tumor cells are often found in clusters. We report a case study on one triple-negative breast cancer (TNBC) patient, analyzing tumor development, metastasis, and prognosis with simultaneous DNA and RNA sequencing of pathologist-defined cell clusters from multiregional frozen sections. The cell clusters are isolated by laser capture microdissection (LCM) from primary tumor tissue, lymphatic vessels, and axillary lymph nodes. Data are reported for a total of 97 cell clusters. A combination of tumor cell-cluster clonality and phylogeny reveals 3 evolutionarily distinct pathways for this patient, each associated with a unique mRNA signature, and each correlated with disparate survival outcomes. Hub gene analysis indicates that extensive downregulation of ribosomal protein mRNA is a potential marker of poor prognosis in breast cancer.

Pathologically diverse cell clusters share genomic and transcriptomic profiles

Transcriptome-defined clones are more complex than genome-defined clones

Three distinct pathways were inferred, each with disparate survival outcomes

Lower expression of ribosomal proteins may be an indicator of poor prognosis

Usefulness of protein-based salivary markers in the diagnosis of oral potentially malignant disorders: A systematic review and meta-analysis

By using a meta-analytical approach, this study aimed to analyse the diagnostic capacity of protein-based biomarkers in saliva for the differential diagnosis of oral potentially malignant disorders (OPMDs) and oral squamous cell carcinoma (OSCC) from healthy individuals as control group (HCG).Articles on protein-based biomarkers in saliva, which provided quantitative expression in individuals with clinical and histopathological diagnosis of OPMD or oral leukoplakia (OL) were considered eligible. Searches were conducted in eight electronic databases. The methodological quality was assessed using the Quality Assessment of Diagnostic Studies tool (QUADAS-2). Functional analysis was also performed. Meta-analyses were performed using the OpenMeta tool (Analyst).Meta-analysis was possible for 4 of the 11 biomarkers studied. Only the carcinoembryonic antigen (CEA) and the soluble fragment of cytokeratin 19 (CYFRA21) were significant for the OSCC/OPMD subgroup, both with a very low heterogeneity. CEA had an OE = 25.854 (CI95%: 13.215-38.492, p< 0.001, I2 = 0) and CYFRA21 had an OE = 9.317 (CI95%: 9.014-9.619, p< 0.001, I2 = 0). For the OPMD/HCG subgroup, only CYFRA21 was significant, with an OE = 3.679 (CI95%: 0.663-6.696, p= 0.017) although with high heterogeneity (I2 = 91.24).The CEA and CYFRA21 markers proved very useful when differentiating OSCC from OPMD. The CYFRA21 was the only protein that was capable of distinguishing between OPMD and healthy controls.

Zika Virus Infection of Human Iris Pigment Epithelial Cells

During recent Zika epidemics, adults infected with Zika virus (ZIKV) have developed organ-specific inflammatory complications. The most serious Zika-associated inflammatory eye disease is uveitis, which is commonly anterior in type, affecting both eyes and responding to corticosteroid eye drops. Mechanisms of Zika-associated anterior uveitis are unknown, but ZIKV has been identified in the aqueous humor of affected individuals. The iris pigment epithelium is a target cell population in viral anterior uveitis, and it acts to maintain immune privilege within the anterior eye. Interactions between ZIKV and human iris pigment epithelial cells were investigated with infectivity assays and RNA-sequencing. Primary cell isolates were prepared from eyes of 20 cadaveric donors, and infected for 24 hours with PRVABC59 strain ZIKV or incubated uninfected as control. Cytoimmunofluorescence, RT-qPCR on total cellular RNA, and focus-forming assays of culture supernatant showed cell isolates were permissive to infection, and supported replication and release of infectious ZIKV. To explore molecular responses of cell isolates to ZIKV infection at the whole transcriptome level, RNA was sequenced on the Illumina NextSeq 500 platform, and results were aligned to the human GRCh38 genome. Multidimensional scaling showed clear separation between transcriptomes of infected and uninfected cell isolates. Differential expression analysis indicated a vigorous molecular response of the cell to ZIKV: 7,935 genes were differentially expressed between ZIKV-infected and uninfected cells (FDR < 0.05), and 99% of 613 genes that changed at least two-fold were up-regulated. Reactome and KEGG pathway and Gene Ontology enrichment analyses indicated strong activation of viral recognition and defense, in addition to biosynthesis processes. A CHAT network included 6275 molecular nodes and 24 contextual hubs in the cell response to ZIKV infection. Receptor-interacting serine/threonine kinase 1 (RIPK1) was the most significantly connected contextual hub. Correlation of gene expression with read counts assigned to the ZIKV genome identified a negative correlation between interferon signaling and viral load across isolates. This work represents the first investigation of mechanisms of Zika-associated anterior uveitis using an in vitro human cell model. The results suggest the iris pigment epithelium mounts a molecular response that limits intraocular pathology in most individuals.

Peptidomics-Driven Strategy Reveals Peptides and Predicted Proteases Associated With Oral Cancer Prognosis

Protease activity has been associated with pathological processes that can lead to cancer development and progression. However, understanding the pathological unbalance in proteolysis is challenging because changes can occur simultaneously at protease, their inhibitor, and substrate levels. Here, we present a pipeline that combines peptidomics, proteomics, and peptidase predictions for studying proteolytic events in the saliva of 79 patients and their association with oral squamous cell carcinoma (OSCC) prognosis. Our findings revealed differences in the saliva peptidome of patients with (pN+) or without (pN0) lymph-node metastasis and delivered a panel of ten endogenous peptides correlated with poor prognostic factors plus five molecules able to classify pN0 and pN+ patients (area under the receiver operating characteristic curve > 0.85). In addition, endopeptidases and exopeptidases putatively implicated in the processing of differential peptides were investigated using cancer tissue gene expression data from public repositories, reinforcing their association with poorer survival rates and prognosis in oral cancer. The dynamics of the OSCC-related proteolysis were further explored via the proteomic profiling of saliva. This revealed that peptidase/endopeptidase inhibitors exhibited reduced levels in the saliva of pN+ patients, as confirmed by selected reaction monitoring-mass spectrometry, while minor changes were detected in the level of saliva proteases. Taken together, our results indicated that proteolytic activity is accentuated in the saliva of patients with OSCC and lymph-node metastasis and, at least in part, is modulated by reduced levels of salivary peptidase inhibitors. Therefore, this integrated pipeline provided better comprehension and discovery of molecular features with implications in the oral cancer metastasis prognosis.

Hypothetical Proteins as Predecessors of Long Non-coding RNAs

Hypothetical Proteins [HP] are the transcripts predicted to be expressed in an organism, but no evidence of it exists in gene banks. On the other hand, long non-coding RNAs [lncRNAs] are the transcripts that might be present in the 5’ UTR or intergenic regions of the genes whose lengths are above 200 bases. With the known unknown [KU] regions in the genomes rapidly existing in gene banks, there is a need to understand the role of open reading frames in the context of annotation. In this commentary, we emphasize that HPs could indeed be the predecessors of lncRNAs.

Identification of key genes and expression profiles in osteoarthritis by co-expressed network analysis

BACKGROUND

The underlying molecular characteristics of osteoarthritis (OA), a common age-related joint disease, remains elusive. Here, we aimed to identify potential early diagnostic biomarkers and elucidate underlying mechanisms of OA using weighted gene co-expression network analysis (WGCNA).

MATERIAL AND METHODS

We obtained the gene expression profile dataset GSE55235, GSE55457, and GSE55584, from the Gene Expression Omnibus. WGCNA was used to investigate the changes in co-expressed genes between normal and OA synovial membrane samples. Modules that were highly correlated to OA were subjected to functional enrichment analysis using the R clusterProfiler package. Differentially expressed genes (DEGs) between the two samples were screened using the "limma" package in R. A Venn diagram was constructed to intersect the genes in significant modules and DEGs. RT -PCR was used to further verify the hub gene expression levels between normal and OA samples.

RESULTS

The preserved significant module was found to be highly associated with OA development and progression (P < 1e-200, correlation = 0.92). Functional enrichment analysis suggested that the antiquewhite4 module was highly correlated to FoxO signaling pathway, and the metabolism of fatty acids and 2-oxocarboxylic acid. A total of 13 hub genes were identified based on significant module network topology and DEG analysis, and RT-PCR confirmed that these genes were significantly increased in OA samples compared with that in normal samples.

CONCLUSIONS

We identified 13 hub genes correlated to the development and progression of OA, which may provide new biomarkers and drug targets for OA.

Hypoxic exosomes facilitate angiogenesis and metastasis in esophageal squamous cell carcinoma through altering the phenotype and transcriptome of endothelial cells

Background

In cancer progression, hypoxia, or low oxygen tension, is a major regulator of tumor aggressiveness and metastasis. However, how cancer cells adapt to the hypoxia and communicate with other mesenchymal cells in microenvironment during tumor development remains to be elucidated. Here, we investigated the involvement of exosomes in modulating angiogenesis and enhancing metastasis in esophageal squamous cell carcinoma (ESCC).

Methods

Differential centrifugation, transmission electron microscopy and nanoparticle tracking analysis were used to isolate and characterize exosomes. Colony formation and transwell assay were performed to assess the proliferation, migration and invasion of human umbilical vein endothelial cells (HUVECs). The tube formation assay and matrigel plug assay were used to evaluate the vascular formation ability of HUVECs in vitro and in vivo respectively. An in vivo nude mice model was established to detect the regulatory role of exosomes in ESCC progression. Microarray analysis was performed to analyze the transcriptome profiles in HUVECs.

Results

Exosomes derived from ESCC cells cultured under hypoxia played a better role in promoting proliferation, migration, invasion and tube formation of HUVECs in vitro and in vivo than exosomes from ESCC cells cultured under normoxia. Moreover, hypoxic exosomes significantly enhanced the tumor growth and lung metastasis compared with normoxic exosomes in nude mice models. Interestingly, endothelial cells were programmed by hypoxic and normoxic exosomes from ESCC cells which altered the transcriptome profile of HUVECs.

Conclusions

Taken together, our data identified an angiogenic role of exosomes from ESCC cells which shed light on the further application of exosomes as valuable therapeutic target for ESCC.

Electronic supplementary material

The online version of this article (10.1186/s13046-019-1384-8) contains supplementary material, which is available to authorized users.

Identification of the iduronate-2-sulfatase proteome in wild-type mouse brain

Iduronate-2-sulfatase (IDS) is a lysosomal enzyme involved in the metabolism of the glycosaminoglycans heparan (HS) and dermatan (DS) sulfate. Mutations on IDS gene produce mucopolysaccharidosis II (MPS II), characterized by the lysosomal accumulation of HS and DS, leading to severe damage of the central nervous system (CNS) and other tissues. In this study, we used a neurochemistry and proteomic approaches to identify the brain distribution of IDS and its interacting proteins on wild-type mouse brain. IDS immunoreactivity showed a robust staining throughout the entire brain, suggesting an intracellular reactivity in nerve cells and astrocytes. By using affinity purification and mass spectrometry we identified 187 putative IDS partners-proteins, mainly hydrolases, cytoskeletal proteins, transporters, transferases, oxidoreductases, nucleic acid binding proteins, membrane traffic proteins, chaperons and enzyme modulators, among others. The interactions with some of these proteins were predicted by using bioinformatics tools and confirmed by co-immunoprecipitation analysis and Blue Native PAGE. In addition, we identified cytosolic IDS-complexes containing proteins from predicted highly connected nodes (hubs), with molecular functions including catalytic activity, redox balance, binding, transport, receptor activity and structural molecule activity. The proteins identified in this study would provide new insights about IDS physiological role into the CNS and its potential role in the brain-specific protein networks.

Immunological Molecular Responses of Human Retinal Pigment Epithelial Cells to Infection With Toxoplasma gondii

Ocular toxoplasmosis is the commonest clinical manifestation of infection with obligate intracellular parasite, Toxoplasma gondii. Active ocular toxoplasmosis is characterized by replication of T. gondii tachyzoites in the retina, with reactive inflammation. The multifunctional retinal pigment epithelium is a key target cell population for T. gondii. Since the global gene expression profile is germane to understanding molecular involvements of retinal pigment epithelial cells in ocular toxoplasmosis, we performed RNA-Sequencing (RNA-Seq) of human cells following infection with T. gondii tachyzoites. Primary cell isolates from eyes of cadaveric donors (n = 3), and the ARPE-19 human retinal pigment epithelial cell line, were infected for 24 h with GT-1 strain T. gondii tachyzoites (multiplicity of infection = 5) or incubated uninfected as control. Total and small RNA were extracted from cells and sequenced on the Illumina NextSeq 500 platform; results were aligned to the human hg19 reference sequence. Multidimensional scaling showed good separation between transcriptomes of infected and uninfected primary cell isolates, which were compared in edgeR software. This differential expression analysis revealed a sizeable response in the total RNA transcriptome-with significantly differentially expressed genes totaling 7,234 (28.9% of assigned transcripts)-but very limited changes in the small RNA transcriptome-totaling 30 (0.35% of assigned transcripts) and including 8 microRNA. Gene ontology and pathway enrichment analyses of differentially expressed total RNA in CAMERA software, identified a strong immunologic transcriptomic signature. We conducted RT-qPCR for 26 immune response-related protein-coding and long non-coding transcripts in epithelial cell isolates from different cadaveric donors (n = 3), extracted by a different isolation protocol but similarly infected with T. gondii, to confirm immunological activity of infected cells. For microRNA, increases in miR-146b and miR-212 were detected by RT-qPCR in 2 and 3 of these independent cell isolates. Biological network analysis in the InnateDB platform, including 735 annotated differentially expressed genes plus 2,046 first-order interactors, identified 10 contextural hubs and 5 subnetworks in the transcriptomic immune response of cells to T. gondii. Our observations provide a solid base for future studies of molecular and cellular interactions between T. gondii and the human retinal pigment epithelium to illuminate mechanisms of ocular toxoplasmosis.

In Silico Protein Interaction Network Analysis of Virulence Proteins Associated with Invasive Aspergillosis for Drug Discovery

BACKGROUND

Protein-Protein interaction (PPI) network analysis of virulence proteins of Aspergillus fumigatus is a prevailing strategy to understand the mechanism behind the virulence of A. fumigatus. The identification of major hub proteins and targeting the hub protein as a new antifungal drug target will help in treating the invasive aspergillosis.

MATERIALS & METHOD

In the present study, the PPI network of 96 virulence (drug target) proteins of A. fumigatus were investigated which resulted in 103 nodes and 430 edges. Topological enrichment analysis of the PPI network was also carried out by using STRING database and Network analyzer a cytoscape plugin app. The key enriched KEGG pathway and protein domains were analyzed by STRING.

CONCLUSION

Manual curation of PPI data identified three proteins (PyrABCN-43, AroM-34, and Glt1- 34) of A. fumigatus possessing the highest interacting partners. Top 10% hub proteins were also identified from the network using cytohubba on the basis of seven algorithms, i.e. betweenness, radiality, closeness, degree, bottleneck, MCC and EPC. Homology model and the active pocket of top three hub proteins were also predicted.

Bioinformatics analysis of microRNA expression between patients with and without latent tuberculosis infections

Tuberculosis (TB) is a globally prevalent infectious disease. The mechanisms of latent TB infection (LTBI) remain to be fully elucidated and may provide novel approaches for diagnosis. As therapeutic targets and molecular diagnostic markers, microRNAs (miRs) have been studied and utilized in various diseases. In the present study, the differentially expressed miRs (DEMs) in LTBI were screened and analyzed to determine the underlying mechanisms and identify potential biomarkers, thereby contributing to the diagnosis of LTBI. The GSE25435 and GSE29190 datasets from Gene Expression Omnibus were selected for analysis. The 2 datasets were analyzed individually using the Bioconductor package to screen the DEMs with specific cut-off criteria [P<0.01 and |log (fold change)|≥1]. Target gene prediction and interaction network construction were performed using Targetscan, the Search Tool for the Retrieval of Interacting Genes and Proteins and Cytoscape individually, and were merged using the latter tool. The hub genes were finally selected based on their degree of connectivity (DC). Gene Ontology and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses were performed using the KEGG and GENCLIP. A total of 144 DEMs were identified from the 2 datasets. By exploring the overlapping miRs in the two datasets, Homo sapiens (hsa)-miR-29a and hsa-miR-15b were identified to be decreased, while hsa-miR-576-5p, hsa-miR-500 and hsa-miR-155 were identified to be upregulated. hsa-miR-500a-3p and hsa-miR-29a-3p, as well as 4 genes, namely cell division cycle (CDC)42, actin α1, skeletal muscle (ACTA1), phosphatase and tensin homolog (PTEN) and fos proto-oncogene (FOS), were selected as the key factors in this regulatory network. A total of 9 signaling pathways, including phosphoinositide-3 kinase (PI3K)/AKT and 11 biological processes, were identified to be associated with LTBI. In conclusion, the present analysis identified hsa-miR-500a-3p and hsa-miR-29a-3p, as well as CDC42, ACTA1, PTEN and FOS, as the most promising biomarkers and therapeutic candidates for LTBI. The PI3K/AKT signaling pathway is the key signaling pathway implicated in LTBI, and an in-depth investigation of the efficiency of PI3K/AKT signaling inhibitors may be used to prevent a chronic state of infection in LTBI.

Systems Genomics of Thigh Adipose Tissue From Asian Indian Type-2 Diabetics Revealed Distinct Protein Interaction Hubs

We performed a systematic analysis of genes implicated in thigh subcutaneous adipose tissue of Asian Indian Type 2 Diabetes Mellitus (AIT2DM) and created a phenome-interactome network. This analysis was performed on 60 subjects specific to limb thigh fat by integrating phenotypic traits and similarity scores associated with AIT2DM. Using a phenotypic attribute, a contextual neighbor was identified across all the traits, viz. body mass index (BMI) statistics, adipocyte size, lipid parameters, homeostatic model assessment- insulin resistance (HOMA-IR), HOMA-ß. In this work, we have attempted to characterize transcription signatures using the phenome-interactome maps where each of the traits under study including the intermediary phenotypes has a distinct set of genes forming the hubs. Furthermore, we have identified various clinical, biochemical, and radiological parameters which show significant correlation with distinct hubs. We observed a number of novel pathways and genes including those that are non-coding RNAs implicated in AIT2DM.We showed that they appear to be associated with pathways, viz. tyrosine kinase JAK2, NOTCH thereby recruiting signaling molecules such as STAT5 and Src family kinases on the cell surface regulated them and our analyses comprising significant hubs suggest that thigh subcutaneous adipose tissue plays a role in pathophysiology of AIT2DM.

Agrin has a pathological role in the progression of oral cancer

Background

The extracellular matrix modulates the hallmarks of cancer. Here we examined the role of agrin—a member of this matrix—in progression of oral squamous cell carcinoma (OSCC).

Methods

We evaluated the immunohistochemical expression of agrin in OSCC and dysplasias. Benign lesions were used as control. In subsequent experiments, we investigated whether the silencing of agrin interferes with tumour expansion both in vitro as well as in vivo. To gain insights into the role of agrin, we identified its protein network (interactome) using mass spectrometry-based proteomics and bioinformatics. Finally, we evaluated the clinical relevance of agrin interactome.

Results

Agrin was elevated in malignant and premalignant lesions. Further, we show that agrin silencing interferes with cancer cell motility, proliferation, invasion, colony and tumour spheroid formation, and it also reduces the phosphorylation of FAK, ERK and cyclin D1 proteins in OSCC cells. In orthotopic model, agrin silencing reduces tumour aggressiveness, like vascular and neural invasion. From a clinical perspective, agrin contextual hubs predict a poor clinical prognosis related with overall survival.

Conclusions

Altogether, our results demonstrate that agrin is a histological marker for the progression of oral cancer and is a strong therapeutic target candidate for both premalignant and OSCC lesions.