STRING v11: protein-protein association networks with increased coverage, supporting functional discovery in genome-wide experimental datasets

Quantitative proteomic analysis reveals high interference on protein expression of H9c2 cells activated with glucose and cardiotonic steroids

In recent decades, the incidence of death and morbidity due to diabetes has increased worldwide, causing a high social and economic impact. Diabetes is a major cause of blindness, kidney failure, heart attack, stroke and lower limb amputation. However, the molecular mechanisms that make the heart and kidneys the main targets of diabetes are not completely understood. To better understand the complex biochemical mechanism of diabetic cardiomyopathy, we investigated the effects of hyperglycemia with concomitant digoxin and ouabain stimulation in H9c2 cells. Total extracted proteins were analyzed by label-free LC-MS/MS, quantified by Scaffold software and validated by parallel reaction monitoring (PRM) methodology. Here, we show that the eukaryotic initiation factors (Eifs) and elongation factors (Eefs) Eif3f, Eef2 and Eif4a1 are overexpressed following cardiotonic steroid (CTS) stimulation. Similarly, the expression of four 14-3-3 proteins that play a key role in cardiac ventricular compaction was altered after CTS stimulation. In total, the expression of nine protein groups was altered in response to the stimulation of H9c2 cells. Here, the biological consequences of these changes are discussed in depth. SIGNIFICANCE: Hyperglycemia is the main physiological condition that provokes tissue and vascular injuries in heart of diabetic patients. However, the changings at large scale in the expression of proteins of cardiomyocytes generated by this condition was not yet studied. Here we report for the first time the altered biosynthesis of nine groups of proteins of H9c2 cells activated by high glucose concentrations and by cardiotonic steroids (CTS). Furthermore, the increased biosynthesis of Eifs, Eefs and 14-3-3 protein groups by CTS, which play a crucial role in cardiomyopathies are original data reported in this work. These findings not only enhance our knowledge concerning to the effects of hyperglycemia and CTS on H9c2 cells but also indicate potential molecular targets to interfere in diabetes cardiomyopathy progression.

Bioinformatics analyses of publicly available NEPCa datasets

Gene expression profiles are valuable resources for the identification of key players that driver disease progression. However, neuroendocrine prostate cancer (NEPCa) specimens are rare, limiting research on this aggressive disease. In this study, we generated a 12-gene signature of NEPCa and used this signature to differentiate NEPCa from prostate adenocarcinoma (AdPCa) samples in publicly available datasets. From these samples, we identified genes that were differentially expressed in NEPCa and AdPCa. Gene ontology and network analyses revealed key players in the pathogenesis of NEPCa, including E2Fs, members of MHC class II, and factors involved in neuron differentiation, neurogenesis, and stem cell signaling. In conclusion, we identified a 12-gene signature of NEPCa and found pathways that are important for the pathologic development of NEPCa.

Phylogenetic correlations can suffice to infer protein partners from sequences

Determining which proteins interact together is crucial to a systems-level understanding of the cell. Recently, algorithms based on Direct Coupling Analysis (DCA) pairwise maximum-entropy models have allowed to identify interaction partners among paralogous proteins from sequence data. This success of DCA at predicting protein-protein interactions could be mainly based on its known ability to identify pairs of residues that are in contact in the three-dimensional structure of protein complexes and that coevolve to remain physicochemically complementary. However, interacting proteins possess similar evolutionary histories. What is the role of purely phylogenetic correlations in the performance of DCA-based methods to infer interaction partners? To address this question, we employ controlled synthetic data that only involve phylogeny and no interactions or contacts. We find that DCA accurately identifies the pairs of synthetic sequences that share evolutionary history. While phylogenetic correlations confound the identification of contacting residues by DCA, they are thus useful to predict interacting partners among paralogs. We find that DCA performs as well as phylogenetic methods to this end, and slightly better than them with large and accurate training sets. Employing DCA or phylogenetic methods within an Iterative Pairing Algorithm (IPA) allows to predict pairs of evolutionary partners without a training set. We further demonstrate the ability of these various methods to correctly predict pairings among real paralogous proteins with genome proximity but no known direct physical interaction, illustrating the importance of phylogenetic correlations in natural data. However, for physically interacting and strongly coevolving proteins, DCA and mutual information outperform phylogenetic methods. We finally discuss how to distinguish physically interacting proteins from proteins that only share a common evolutionary history.

Many biologically important protein-protein interactions are conserved over evolutionary time scales. This leads to two different signals that can be used to computationally predict interactions between protein families and to identify specific interaction partners. First, the shared evolutionary history leads to highly similar phylogenetic relationships between interacting proteins of the two families. Second, the need to keep the interaction surfaces of partner proteins biophysically compatible causes a correlated amino-acid usage of interface residues. Employing simulated data, we show that the shared history alone can be used to detect partner proteins. Similar accuracies are achieved by algorithms comparing phylogenetic relationships and by methods based on Direct Coupling Analysis (DCA), which are primarily known for their ability to detect the second type of signal. Using natural sequence data, we show that in cases with shared evolutionary history but without known physical interactions, both methods work with similar accuracy, while for some physically interacting systems, DCA and mutual information outperform phylogenetic methods. We propose methods allowing both to predict interactions between protein families and to find interacting partners among paralogs.

High gene expression levels of VEGFA and CXCL8 in the peritumoral brain zone are associated with the recurrence of glioblastoma: A bioinformatics analysis

The present study aimed to identify differentially regulated genes between the peritumoral brain zone (PBZ) and tumor core (TC) of glioblastoma (GBM), to elucidate the underlying molecular mechanisms and provide a target for the treatment of tumors. The GSE13276 and GSE116520 datasets were downloaded from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) for the PBZ and TC were obtained using the GEO2R tool. The bioinformatics and evolutionary genomics online tool Venn was used to identify common DEGs between the two datasets. The Database for Annotation, Visualization, and Integrated Discovery online tool was used to analyze enriched pathways of the Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases. The Search Tool for the Retrieval of Interacting Genes/Proteins online tool was used to construct a protein-protein interaction (PPI) network of DEGs. Hub genes were identified using Cytohubba, a plug-in for Cytoscape. The Gene Expression Profiling Interactive Analysis (GEPIA) database was utilized to perform survival analysis. In total, 75 DEGs, including 12 upregulated and 63 downregulated genes, were identified. In the GO term analysis, these DEGs were mainly enriched in ‘regulation of angiogenesis’ and ‘central nervous system development’. Furthermore, in the KEGG pathway analysis, the DEGs were mainly enriched in ‘bladder cancer’ and ‘endocytosis’. When filtering the results of the PPI network analysis using Cytohubba, a total of 10 hub genes, including proteolipid protein 1, myelin associated oligodendrocyte basic protein, contactin 2, myelin oligodendrocyte glycoprotein, myelin basic protein, myelin associated glycoprotein, SRY-box transcription factor 10, C-X-C motif chemokine ligand 8 (CXCL8), vascular endothelial growth factor A (VEGFA) and plasmolipin, were identified. These hub genes were further subjected to GO term and KEGG pathway analysis, and were revealed to be enriched in ‘central nervous system development’, ‘bladder cancer’ and ‘rheumatoid arthritis’. These hub genes were used to perform survival analysis using the GEPIA database, and it was determined that VEGFA and CXCL8 were significantly associated with a reduction in the overall survival of patients with GBM. In conclusion, the results suggest that the recurrence of GBM is associated with high gene expression levels VEGFA and CXCL8, and the development of the central nervous system.

Biological Network Approaches and Applications in Rare Disease Studies

Network biology has the capability to integrate, represent, interpret, and model complex biological systems by collectively accommodating biological omics data, biological interactions and associations, graph theory, statistical measures, and visualizations. Biological networks have recently been shown to be very useful for studies that decipher biological mechanisms and disease etiologies and for studies that predict therapeutic responses, at both the molecular and system levels. In this review, we briefly summarize the general framework of biological network studies, including data resources, network construction methods, statistical measures, network topological properties, and visualization tools. We also introduce several recent biological network applications and methods for the studies of rare diseases.

The transcriptome difference between colorectal tumor and normal tissues revealed by single-cell sequencing

The previous cancer studies were difficult to reproduce since the tumor tissues were analyzed directly. But the tumor tissues were actually a mixture of different cancer cells. The transcriptome of single-cell was much robust than the transcriptome of a mixed tissue. The single-cell transcriptome had much smaller variance. In this study, we analyzed the single-cell transcriptome of 272 colorectal cancer (CRC) epithelial cells and 160 normal epithelial cells and identified 342 discriminative transcripts using advanced machine learning methods. The most discriminative transcripts were LGALS4, PHGR1, C15orf48, HEPACAM2, PERP, FABP1, FCGBP, MT1G, TSPAN1 and CKB. We further clustered the 342 transcripts into two categories. The upregulated transcripts in CRC epithelial cells were significantly enriched in Ribosome, Protein processing in endoplasmic reticulum, Antigen processing and presentation and p53 signaling pathway. The downregulated transcripts in CRC epithelial cells were significantly enriched in Mineral absorption, Aldosterone-regulated sodium reabsorption and Oxidative phosphorylation pathways. The biological analysis of the discriminative transcripts revealed the possible mechanism of colorectal cancer.

Genome-wide investigation of superoxide dismutase (SOD) gene family and their regulatory miRNAs reveal the involvement in abiotic stress and hormone response in tea plant (Camellia sinensis)

Superoxide dismutases (SODs), as a family of metalloenzymes related to the removal of reactive oxygen species (ROS), have not previously been investigated at genome-wide level in tea plant. In this study, 10 CsSOD genes were identified in tea plant genome, including 7 Cu/Zn-SODs (CSDs), 2 Fe-SODs (FSDs) and one Mn-SOD (MSD), and phylogenetically classified in three subgroups, respectively. Physico-chemical characteristic, conserved motifs and potential protein interaction analyses about CsSOD proteins were carried out. Exon-intron structures and codon usage bias about CsSOD genes were also examined. Exon-intron structures analysis revealed that different CsSOD genes contained various number of introns. On the basis of the prediction of regulatory miRNAs of CsSODs, a modification 5’ RNA ligase-mediated (RLM)-RACE was performed and validated that csn-miR398a-3p-1 directly cleaves CsCSD4. By prediction of cis-acting elements, the expression patterns of 10 CsSOD genes and their regulatory miRNAs were detected under cold, drought, exogenous methyl jasmonate (MeJA) and gibberellin (GA₃) treatments. The results showed that most of CsSODs except for CsFSD2 were induced under cold stress and CsCSDs may play primary roles under drought stress; exogenous GA₃ and MeJA could also stimulated/inhibited distinct CsSODs at different stages. In addition, we found that csn-miR398a-3p-1 negatively regulated the expression of CsCSD4 may be a crucial regulatory mechanism under cold stress. This study provides a certain basis for the studies about stress resistance in tea plants, even provide insight into comprehending the classification, evolution, diverse functions and influencing factors of expression patterns for CsSOD genes.

Analysis of predicted amino acid biosynthesis in Rv3344c in Mycobacterium tuberculosis H37 Rv using bioinformatics tools

According to World Health Organization (WHO) report, Mycobacterium tuberculosis H₃₇ Rv (M. tuberculosis) affects one-third population of the world. Emergence of effective treatment/research against this disease is need of the hour. Therefore, we present some important aspects of Rv3344c, which is a PE_PGRS protein. Evidence shows that PE_PGRS proteins show fibronectin binding activity. This protein has affinity for calcium and also shows motifs of GTP-binding protein. It also shows the presence of sites for ribose-5-phosphate binding and motifs of aspartate-beta-semialdehyde dehydrogenase, both of which are involved in amino acid biosynthesis. Thus, this protein might be targeted to block the amino acid biosynthesis in M. tuberculosis. This article takes into consideration some important aspects of Rv3344c protein as its function is still unknown. This study includes retrieval of protein sequence database, multiple sequence alignment, protein-protein interaction, epitope prediction, localization, function prediction, phosphorylation site prediction, model building and its validation, ligand-binding prediction along with mutational analysis. Hence, this study might be an important step in the development of new drugs and treatment of tuberculosis.

Asparagine Deprivation Causes a Reversible Inhibition of Human Cytomegalovirus Acute Virus Replication

As obligate intracellular pathogens, viruses rely on the host cell machinery to replicate efficiently, with the host metabolism extensively manipulated for this purpose. High-throughput small interfering RNA (siRNA) screens provide a systematic approach for the identification of novel host-virus interactions. Here, we report a large-scale screen for host factors important for human cytomegalovirus (HCMV), consisting of 6,881 siRNAs. We identified 47 proviral factors and 68 antiviral factors involved in a wide range of cellular processes, including the mediator complex, proteasome function, and mRNA splicing. Focused characterization of one of the hits, asparagine synthetase (ASNS), demonstrated a strict requirement for asparagine for HCMV replication which leads to an early block in virus replication before the onset of DNA amplification. This effect is specific to HCMV, as knockdown of ASNS had little effect on herpes simplex virus 1 or influenza A virus replication, suggesting that the restriction is not simply due to a failure in protein production. Remarkably, virus replication could be completely rescued 7 days postinfection with the addition of exogenous asparagine, indicating that while virus replication is restricted at an early stage, it maintains the capacity for full replication days after initial infection. This study represents the most comprehensive siRNA screen for the identification of host factors involved in HCMV replication and identifies the nonessential amino acid asparagine as a critical factor in regulating HCMV virus replication. These results have implications for control of viral latency and the clinical treatment of HCMV in patients.IMPORTANCE HCMV accounts for more than 60% of complications associated with solid organ transplant patients. Prophylactic or preventative treatment with antivirals, such as ganciclovir, reduces the occurrence of early onset HCMV disease. However, late onset disease remains a significant problem, and prolonged treatment, especially in patients with suppressed immune systems, greatly increases the risk of antiviral resistance. Very few antivirals have been developed for use against HCMV since the licensing of ganciclovir, and of these, the same viral genes are often targeted, reducing the usefulness of these drugs against resistant strains. An alternative approach is to target host genes essential for virus replication. Here we demonstrate that HCMV replication is highly dependent on levels of the amino acid asparagine and that knockdown of a critical enzyme involved in asparagine synthesis results in severe attenuation of virus replication. These results suggest that reducing asparagine levels through dietary restriction or chemotherapeutic treatment could limit HCMV replication in patients.

Boronate affinity-based photoactivatable magnetic nanoparticles for the oriented and irreversible conjugation of Fc-fused lectins and antibodies

The utilization of immuno-magnetic nanoparticles (MNPs) for the selective capture, enrichment, and separation of specific glycoproteins from complicated biological samples is appealing for the discovery of disease biomarkers. Herein, MNPs were designed and anchored with abundant boronic acid (BA) and photoreactive alkyl diazirine (Diaz) functional groups to obtain permanently tethered Fc-fused Siglec-2 and antiserum amyloid A (SAA) mAb with the assistance of reversible boronate affinity and UV light activation in an orientation-controlled manner. The Siglec-2-Fc-functionalized MNPs showed excellent stability in fetal bovine serum (FBS) and excellent efficiency in the extraction of cell membrane glycoproteins. The anti-SAA mAb-functionalized MNPs maintained active Ab orientation and preserved antigen recognition capability in biological samples. Thus, the BA-Diaz-based strategy holds promise for the immobilization of glycoproteins, such as antibodies, with the original protein binding activity maintained, which can provide better enrichment for the sensitive detection of target proteins.

R2R3 MYB Transcription Factors - Functions outside the DNA-Binding Domain

Several transcription factor (TF) families, including the MYB family, regulate a wide array of biological processes. TFs contain DNA-binding domains (DBDs) and regulatory regions; although information on protein structure is scarce for plant MYB TFs, various in silico methods suggest that the non-MYB regions contain extensive intrinsically disordered regions (IDRs). Although IDRs do not fold into stable globular structures, they comprise functional regions including interaction motifs, and recent research has shown that IDRs perform crucial biological roles. We map here domain organization, disorder predictions, and functional regions across the entire Arabidopsis thaliana R2R3 MYB TF family, and highlight where an increased research focus will be necessary to shape a new understanding of structure-function relationships in plant TFs.

Redefining malignant pleural mesothelioma types as a continuum uncovers immune-vascular interactions

BACKGROUND

Malignant Pleural Mesothelioma (MPM) is an aggressive disease related to asbestos exposure, with no effective therapeutic options.

METHODS

We undertook unsupervised analyses of RNA-sequencing data of 284 MPMs, with no assumption of discreteness. Using immunohistochemistry, we performed an orthogonal validation on a subset of 103 samples and a biological replication in an independent series of 77 samples.

FINDINGS

A continuum of molecular profiles explained the prognosis of the disease better than any discrete model. The immune and vascular pathways were the major sources of molecular variation, with strong differences in the expression of immune checkpoints and pro-angiogenic genes; the extrema of this continuum had specific molecular profiles: a "hot" bad-prognosis profile, with high lymphocyte infiltration and high expression of immune checkpoints and pro-angiogenic genes; a "cold" bad-prognosis profile, with low lymphocyte infiltration and high expression of pro-angiogenic genes; and a "VEGFR2+/VISTA+" better-prognosis profile, with high expression of immune checkpoint VISTA and pro-angiogenic gene VEGFR2. We validated the gene expression levels at the protein level for a subset of five selected genes belonging to the immune and vascular pathways (CD8A, PDL1, VEGFR3, VEGFR2, and VISTA), in the validation series, and replicated the molecular profiles as well as their prognostic value in the replication series.

INTERPRETATION

The prognosis of MPM is best explained by a continuous model, which extremes show specific expression patterns of genes involved in angiogenesis and immune response.

FGB and FGG derived from plasma exosomes as potential biomarkers to distinguish benign from malignant pulmonary nodules

Previous proteomic analysis (label-free) of plasma exosomes revealed that the expression of FGG and FGB was significantly higher in the malignant pulmonary nodules group, compared to the benign pulmonary nodules group. The present study was performed to evaluate the role of plasma exosomal proteins FGB and FGG in the diagnosis of benign and malignant pulmonary nodules. We examined the expression levels of FGB and FGG in plasma exosomes from 63 patients before surgery. Postoperative pathological diagnosis confirmed that 43 cases were malignant and 20 cases were benign. The ROC curve was used to describe the sensitivity, specificity, area under the curve (AUC) of the biomarker and the corresponding 95% confidence interval. We confirmed that the expression levels of FGB and FGG were higher in the plasma exosomes of malignant group than in the benign group. The sensitivity and AUC of FGB combined with FGG detection to determine the nature of pulmonary nodules are superior to single FGB or FGG detection. FGB and FGG might represent novel and sensitive biomarker to distinguish benign from malignant pulmonary nodules.

Overexpression of GmCAMTA12 Enhanced Drought Tolerance in Arabidopsis and Soybean

Fifteen transcription factors in the CAMTA (calmodulin binding transcription activator) family of soybean were reported to differentially regulate in multiple stresses; however, their functional analyses had not yet been attempted. To characterize their role in stresses, we first comprehensively analyzed the GmCAMTA family in silico and thereafter determined their expression pattern under drought. The bioinformatics analysis revealed multiple stress-related cis-regulatory elements including ABRE, SARE, G-box and W-box, 10 unique miRNA (microRNA) targets in GmCAMTA transcripts and 48 proteins in GmCAMTAs' interaction network. We then cloned the 2769 bp CDS (coding sequence) of GmCAMTA12 in an expression vector and overexpressed in soybean and Arabidopsis through Agrobacterium-mediated transformation. The T3 (Transgenic generation 3) stably transformed homozygous lines of Arabidopsis exhibited enhanced tolerance to drought in soil as well as on MS (Murashige and Skoog) media containing mannitol. In their drought assay, the average survival rate of transgenic Arabidopsis lines OE5 and OE12 (Overexpression Line 5 and Line 12) was 83.66% and 87.87%, respectively, which was ~30% higher than that of wild type. In addition, the germination and root length assays as well as physiological indexes such as proline and malondialdehyde contents, catalase activity and leakage of electrolytes affirmed the better performance of OE lines. Similarly, GmCAMTA12 overexpression in soybean promoted drought-efficient hairy roots in OE chimeric plants as compare to that of VC (Vector control). In parallel, the improved growth performance of OE in Hoagland-PEG (polyethylene glycol) and on MS-mannitol was revealed by their phenotypic, physiological and molecular measures. Furthermore, with the overexpression of GmCAMTA12, the downstream genes including AtAnnexin5, AtCaMHSP, At2G433110 and AtWRKY14 were upregulated in Arabidopsis. Likewise, in soybean hairy roots, GmELO, GmNAB and GmPLA1-IId were significantly upregulated as a result of GmCAMTA12 overexpression and majority of these upregulated genes in both plants possess CAMTA binding CGCG/CGTG motif in their promoters. Taken together, we report that GmCAMTA12 plays substantial role in tolerance of soybean against drought stress and could prove to be a novel candidate for engineering soybean and other plants against drought stress. Some research gaps were also identified for future studies to extend our comprehension of Ca-CaM-CAMTA-mediated stress regulatory mechanisms.

Integration of Bioinformatics Resources Reveals the Therapeutic Benefits of Gemcitabine and Cell Cycle Intervention in SMAD4-Deleted Pancreatic Ductal Adenocarcinoma

Pancreatic ductal adenocarcinoma (PDAC) is the most common and aggressive type of pancreatic cancer. The five-year survival rate of PDAC is very low (less than 8%), which is associated with the late diagnosis, high metastatic potential, and resistance to therapeutic agents. The identification of better prognostic or therapeutic biomarker may have clinical benefits for PDAC treatment. SMAD4, a central mediator of transforming growth factor beta (TGFβ) signaling pathway, is considered a tumor suppressor gene. SMAD4 inactivation is frequently found in PDAC. However, its role in prognosis and therapeutics of PDAC is still unclear. In this study, we applied bioinformatics approaches, and integrated publicly available resources, to investigate the role of SMAD4 gene deletion in PDAC. We found that SMAD4 deletion was associated with poorer disease-free, but not overall, survival in PDAC patients. Cancer hallmark enrichment and pathway analysis suggested that the upregulation of cell cycle-related genes in SMAD4-deleted PDAC. Chemotherapy response profiling of PDAC cell lines and patient-derived organoids revealed that SMAD4-deleted PDAC was sensitive to gemcitabine, the first-line treatment for PDAC, and specific cell cycle-targeting drugs. Taken together, our study provides an insight into the prognostic and therapeutic roles of SMAD4 gene deletion in PDAC, and SMAD4 gene copy numbers may be used as a therapeutic biomarker for PDAC treatment.

Enhancing cardiac glycolysis causes an increase in PDK4 content in response to short-term high-fat diet

The healthy heart has a dynamic capacity to respond and adapt to changes in nutrient availability. Metabolic inflexibility, such as occurs with diabetes, increases cardiac reliance on fatty acids to meet energetic demands, and this results in deleterious effects, including mitochondrial dysfunction, that contribute to pathophysiology. Enhancing glucose usage may mitigate metabolic inflexibility and be advantageous under such conditions. Here, we sought to identify how mitochondrial function and cardiac metabolism are affected in a transgenic mouse model of enhanced cardiac glycolysis (Glyco^Hi) basally and following a short-term (7-day) high-fat diet (HFD). Glyco^Hi mice constitutively express an active form of phosphofructokinase-2, resulting in elevated levels of the PFK-1 allosteric activator fructose 2,6-bisphosphate. We report that basally Glyco^Hi mitochondria exhibit augmented pyruvate-supported respiration relative to fatty acids. Nevertheless, both WT and Glyco^Hi mitochondria had a similar shift toward increased rates of fatty acid-supported respiration following HFD. Metabolic profiling by GC-MS revealed distinct features based on both genotype and diet, with a unique increase in branched-chain amino acids in the Glyco^Hi HFD group. Targeted quantitative proteomics analysis also supported both genotype- and diet-dependent changes in protein expression and uncovered an enhanced expression of pyruvate dehydrogenase kinase 4 (PDK4) in the Glyco^Hi HFD group. These results support a newly identified mechanism whereby the levels of fructose 2,6-bisphosphate promote mitochondrial PDK4 levels and identify a secondary adaptive response that prevents excessive mitochondrial pyruvate oxidation when glycolysis is sustained after a high-fat dietary challenge.

The Construction and Comprehensive Analysis of ceRNA Networks and Tumor-Infiltrating Immune Cells in Bone Metastatic Melanoma

Background/Aims: As a malignant and melanocytic tumor, cutaneous melanoma is the devastating skin tumor with high rates of recurrence and metastasis. Bone is the common metastatic location, and bone metastasis may result in pathologic fracture, neurologic damage, and severe bone pain. Although metastatic melanoma was reported to get benefits from immunotherapy, molecular mechanisms and immune microenviroment underlying the melanoma bone metastasis and prognostic factors are still unknown.

Methods: Gene expression profiling of 112 samples, including 104 primary melanomas and 8 bone metastatic melanomas from The Cancer Genome Atlas database, was assayed to construct a ceRNA network associated with bone metastases. Besides, we detected the fraction of 22 immune cell types in melanoma via the algorithm of “cell type identification by estimating relative subsets of RNA transcripts (CIBERSORT).” Based on the significant ceRNAs or immune cells, we constructed nomograms to predict the prognosis of patients with melanoma. Ultimately, correlation analysis was implemented to discover the relationship between the significant ceRNA and immune cells to reveal the potential signaling pathways.

Results: We constructed a ceRNA network based on the interaction among 8 pairs of long noncoding RNA–microRNA and 15 pairs of microRNA–mRNA. CIBERSORT and ceRNA integration analysis discovered that AL118506.1 has both significant prognostic value (P = 0.002) and high correlation with T follicular helper cells (P = 0.033). Meanwhile, T cells CD8 and macrophages M2 were negatively correlated (P < 0.001). Moreover, we constructed two satisfactory nomograms (area under curve of 3-year survival: 0.899; 5-year survival: 0.885; and concordance index: 0.780) with significant ceRNAs or immune cells, to predict the prognosis of patients.

Conclusions: In this study, we suggest that bone metastasis in melanoma might be related to AL118506.1 and its role in regulating thrombospondin 2 and T follicular helper cells. Two nomograms were constructed to predict the prognosis of patients with melanoma and demonstrated their value in improving the personalized management.

Proteomics-Based Monitoring of Pathway Activity Reveals that Blocking Diacylglycerol Biosynthesis Rescues from Alpha-Synuclein Toxicity

Proteinaceous inclusions containing alpha-synuclein (α-Syn) have been implicated in neuronal toxicity in Parkinson's disease, but the pathways that modulate toxicity remain enigmatic. Here, we used a targeted proteomic assay to simultaneously measure 269 pathway activation markers and proteins deregulated by α-Syn expression across a panel of 33 Saccharomyces cerevisiae strains that genetically modulate α-Syn toxicity. Applying multidimensional linear regression analysis to these data predicted Pah1, a phosphatase that catalyzes conversion of phosphatidic acid to diacylglycerol at the endoplasmic reticulum membrane, as an effector of rescue. Follow-up studies demonstrated that inhibition of Pah1 activity ameliorates the toxic effects of α-Syn, indicate that the diacylglycerol branch of lipid metabolism could enhance α-Syn neuronal cytotoxicity, and suggest a link between α-Syn toxicity and the biology of lipid droplets.

Identification of the lncRNA-miRNA-mRNA network associated with gastric cancer via integrated bioinformatics analysis

The aim of the present study was to investigate the long non-coding RNA (lncRNA)-microRNA (miRNA)-mRNA regulatory network in gastric cancer (GC) using bioinformatics analysis. Two mRNA gene expression profiles, GSE79973 and GSE54129, and two miRNA expression profiles, GSE93415 and GSE78091, were downloaded from the Gene Expression Omnibus database. The differentially expressed mRNAs (DEMs) and the differentially expressed miRNAs (DEMis) were merged separately. Gene ontology and pathway enrichment analysis were conducted using the Database for Annotation, Visualization and Integrated Discovery. A protein-protein interaction (PPI) network was then constructed and the 10 top hub genes in the network were analyzed using the Search Tool for the Retrieval of Interacting Genes. The lncRNA-miRNA-mRNA networks were visualized using Cytoscape software. As a result, 158 shared DEMs (40 upregulated and 118 downregulated) were identified from two mRNA datasets. A total of 30 upregulated miRNAs and 1 downregulated miRNA functioned as DEMis. The PPI network consisted of 129 nodes and 572 interactions. The 10 top hub genes were selected by degree using Cytohubba, including Jun proto-oncogene, mitogen-activated protein kinase (MAPK)3, transforming growth factor-β1, Fos proto-oncogene, AP-1 transcription factor subunit, interleukin (IL)-8, MAPK1, RELA proto-oncogene nuclear factor-κB subunit, interferon regulatory factor 7, ubiquitin like modifier and vascular endothelial growth factor A. In the lncRNA-miRNA-mRNA network, a total of 1,215 regulatory associations were constructed using Cytoscape. In conclusion, the present study provides a novel perspective of the molecular mechanisms underlying GC by identifying the lncRNA-miRNA-mRNA regulatory network via bioinformatics analysis.

Induction of Lysosome-associated Protein Transmembrane 4 Beta via Sulfatase 2 Enhances Autophagic Flux in Liver Cancer Cells

Autophagy has been shown to be a key cellular event controlling tumor growth in different neoplasms including hepatocellular carcinoma (HCC). Although this biological role of autophagy has been clearly established, the mechanism underlying its regulation remains elusive. Here, we demonstrate a role of sulfatase 2 (SULF2), a 6‐O‐endosulfatase modulating various growth factors and cytokine‐related signaling pathways controlling tumor cell proliferation and survival, in the regulation of autophagy in HCC cells. SULF2 increased autophagosome formation, shown by increased LC3B‐II protein and green fluorescent protein–LC3 puncta. Increased fusion between autophagosomes and lysosomes/lysosomal enzymes, higher expression of lysosomal membrane protein, and an increase in autolysosomes were also shown by western blot, immunofluorescence, and electron microscopy of SULF2‐expressing cells, indicating enhanced autophagic flux. In contrast, RNA‐interference silencing of SULF2 in Huh7 cells induced lysosomal membrane permeabilization with diffuse cytosolic staining of cathepsin D and punctate staining of galectin‐3. Analysis of the mechanism showed that inhibition of lysosome‐associated protein transmembrane 4 beta (LAPTM4B), a gene induced by SULF2, resulted in decreased autophagosome formation, decreased fusion between autophagosomes and lysosomes, and increased lysosomal membrane permeabilization. Interestingly, down‐regulation of LAPTM4B also phenocopies the knockdown of SULF2, significantly reducing cell viability and colony formation. Conclusion: Our results demonstrate a role for SULF2 in the regulation of autophagic flux that is mediated through LAPTM4B induction in HCC cells, and provide a foundation for future translational efforts targeting autophagy in liver malignancies.

High Expression of FGF5 Is an Independent Prognostic Factor for Poor Overall Survival and Relapse-Free Survival in Lung Adenocarcinoma

Lung cancer is not only a serious disease but also a public problem threatening human health all over the world. Nonsmall cell lung cancer-which accounts for the majority of lung cancer-is mainly composed of lung adenocarcinoma (LUAD) and squamous cell carcinoma (LUSC). FGF5 is a gene located in q21.21. In the past years, research on FGF5 is mainly focused on hair length and hereditary spherocytosis. In our study, bioinformatics analysis of FGF5 was performed through multiple databases. Expression of FGF5 was compared between tumor and normal tissues, association between gene expression and clinical outcomes was investigated in LUAD and LUSC separately, and potential signaling pathways were predicted. FGF5 expression was upregulated in lung cancer tissues compared with normal tissues. What is more, the high FGF5 expression group had significantly lower proportions of lymph node negative (N0) patients (77/144, 53.5%, vs. 253/358, 70.7%, p = 0.000), and is associated with worse overall survival (OS) (p < 0.0001) and relapse-free survival (RFS) (p = 0.024) in LUAD patients, which could not be seen in LUSC. The following analysis confirmed that high FGF5 expression could be an independent prognostic factor for poor OS (HR: 0.431, 95% CI: 0.312-0.597, p = 0.001) and RFS (HR: 0.678, 95% CI: 0.471-0.977, p = 0.037) in LUAD, but not in LUSC. Coexpression genes related to FGF5 were explored and potential pathways were investigated for further research. FGF5 is a tumor-associated gene that upregulated in lung cancer tissues, and could be an independent prognostic factor that have potential value for further research.

Human-specific LAIR2 contributes to the high invasiveness of human extravillous trophoblast cells

The placenta is a temporary vital organ for intra-uterine development and growth. The anatomical structure of the placenta has evolved substantially, resulting in broad inter-species diversity. In particular, human placental extravillous trophoblast cells (EVTs) have evolved aggressive features, although the mechanism underlying this aggressiveness remains elusive. In the present study, we compared the human and mouse homologous gene databases and obtained 2272 human-specific genes, 807 of which are expressed in the placenta according to the UniGene database. Using the human trophoblast cell line HTR8/SVneo, we further verified the expression and function of one of these genes, the leukocyte-associated immunoglobulin-like receptor 2 (LAIR2). This gene shows increased expression during pregnancy and its reduced expression is associated with pregnancy complications. Although LAIR2 was expressed in the human placenta villus and decidua in the first trimester of pregnancy, it was not expressed in mouse tissues. Knockdown of LAIR2 markedly improved cell viability and inhibited the invasive ability of HTR8/SVneo cells. These data suggest that species-specific genes are pivotal to the evolution of a more aggressive human placenta to match the physiological demands of human development. Further investigation is required to obtain evidence on the function of LAIR2 and other specific genes in the placenta, providing insight on the mechanism, properties, and possible applications of this in humans.

Role of miR-10b-5p in the prognosis of breast cancer

Breast cancer is the leading cause of cancer-related death in women worldwide. Aberrant expression levels of miR-10b-5p in breast cancer has been reported while the molecular mechanism of miR-10b-5p in tumorigenesis remains elusive. Therefore, this study was aimed to investigate the role of miR-10b-5p in breast cancer and the network of its target genes using bioinformatics analysis. In this study, the expression profiles and prognostic value of miR-10b-5p in breast cancer were analyzed from public databases. Association between miR-10b-5p and clinicopathological parameters were analyzed by non-parametric test. Moreover, the optimal target genes of miR-10b-5p were obtained and their expression patterns were examined using starBase and HPA database. Additionally, the role of these target genes in cancer development were explored via Cancer Hallmarks Analytics Tool (CHAT). The protein–protein interaction (PPI) networks were constructed to further investigate the interactive relationships among these genes. Furthermore, GO, KEGG pathway and Reactome pathway analyses were carried out to decipher functions of these target genes. Results demonstrated that miR-10b-5p was down-regulated in breast cancer and low expression of miR-10b-5p was significantly correlated to worse outcome. Five genes, BIRC5, E2F2, KIF2C, FOXM1, and MCM5, were considered as potential key target genes of miR-10b-5p. As expected, higher expression levels of these genes were observed in breast cancer tissues than in normal tissues. Moreover, analysis from CHAT revealed that these genes were mainly involved in sustaining proliferative signaling in cancer development. In addition, PPI networks analysis revealed strong interactions between target genes. GO, KEGG, and Reactome pathway analysis suggested that these target genes of miR-10b-5p in breast cancer were significantly involved in cell cycle. Predicted target genes were further validated by qRT-PCR analysis in human breast cancer cell line MDA-MB-231 transfected with miR-10b mimic or antisense inhibitors. Taken together, our data suggest that miR-10b-5p functions to impede breast carcinoma progression via regulation of its key target genes and hopefully serves as a potential diagnostic and prognostic marker for breast cancer.

Serum-Derived Extracellular Vesicles from African Swine Fever Virus-Infected Pigs Selectively Recruit Viral and Porcine Proteins

: African swine fever is a devastating hemorrhagic infectious disease, which affects domestic and wild swines (Susscrofa) of all breeds and ages, with a high lethality of up to 90-100% in naïve animals. The causative agent, African swine fever virus (ASFV), is a large and complex double-stranded DNA arbovirus which is currently spreading worldwide, with serious socioeconomic consequences. There is no treatment or effective vaccine commercially available, and most of the current research is focused on attenuated viral models, with limited success so far. Thus, new strategies are under investigation. Extracellular vesicles (EVs) have proven to be a promising new vaccination platform for veterinary diseases in situations in which conventional approaches have not been completely successful. Here, serum extracellular vesicles from infected pigs using two different ASFV viruses (OURT 88/3 and Benin ΔMGF), corresponding to a naturally attenuated virus and a deletion mutant, respectively, were characterized in order to determine possible differences in the content of swine and viral proteins in EV-enriched fractions. Firstly, EVs were characterized by their CD5, CD63, CD81 and CD163 surface expression. Secondly, ASFV proteins were detected on the surface of EVs from ASFV-infected pig serum. Finally, proteomic analysis revealed few specific proteins from ASFV in the EVs, but 942 swine proteins were detected in all EV preparations (negative controls, and OURT 88/3 and Benin ΔMGF-infected preparations). However, in samples from OURT 88/3-infected animals, only a small number of proteins were differentially identified compared to control uninfected animals. Fifty-six swine proteins (Group Benin) and seven proteins (Group OURT 88/3) were differentially detected on EVs when compared to the EV control group. Most of these were related to coagulation cascades. The results presented here could contribute to a better understanding of ASFV pathogenesis and immune/protective responses in the host.

Altered Fecal Small RNA Profiles in Colorectal Cancer Reflect Gut Microbiome Composition in Stool Samples

The characteristics of microbial small RNA transcription are largely unknown, while it is of primary importance for a better identification of molecules with functional activities in the gut niche under both healthy and disease conditions. By performing combined analyses of metagenomic and small RNA sequencing (sRNA-Seq) data, we characterized both the human and microbial small RNA contents of stool samples from healthy individuals and from patients with colorectal carcinoma or adenoma. With the integrative analyses of metagenomic and sRNA-Seq data, we identified a human and microbial small RNA signature which can be used to improve diagnosis of the disease. Our analysis of human and gut microbiome small RNA expression is relevant to generation of the first hypotheses about the potential molecular interactions occurring in the gut of CRC patients, and it can be the basis for further mechanistic studies and clinical tests.

Dysbiotic configurations of the human gut microbiota have been linked to colorectal cancer (CRC). Human small noncoding RNAs are also implicated in CRC, and recent findings suggest that their release in the gut lumen contributes to shape the gut microbiota. Bacterial small RNAs (bsRNAs) may also play a role in carcinogenesis, but their role has been less extensively explored. Here, we performed small RNA and shotgun sequencing on 80 stool specimens from patients with CRC or with adenomas and from healthy subjects collected in a cross-sectional study to evaluate their combined use as a predictive tool for disease detection. We observed considerable overlap and a correlation between metagenomic and bsRNA quantitative taxonomic profiles obtained from the two approaches. We identified a combined predictive signature composed of 32 features from human and microbial small RNAs and DNA-based microbiome able to accurately classify CRC samples separately from healthy and adenoma samples (area under the curve [AUC] = 0.87). In the present study, we report evidence that host-microbiome dysbiosis in CRC can also be observed by examination of altered small RNA stool profiles. Integrated analyses of the microbiome and small RNAs in the human stool may provide insights for designing more-accurate tools for diagnostic purposes.

IMPORTANCE The characteristics of microbial small RNA transcription are largely unknown, while it is of primary importance for a better identification of molecules with functional activities in the gut niche under both healthy and disease conditions. By performing combined analyses of metagenomic and small RNA sequencing (sRNA-Seq) data, we characterized both the human and microbial small RNA contents of stool samples from healthy individuals and from patients with colorectal carcinoma or adenoma. With the integrative analyses of metagenomic and sRNA-Seq data, we identified a human and microbial small RNA signature which can be used to improve diagnosis of the disease. Our analysis of human and gut microbiome small RNA expression is relevant to generation of the first hypotheses about the potential molecular interactions occurring in the gut of CRC patients, and it can be the basis for further mechanistic studies and clinical tests.

Network Pharmacology Reveals the Molecular Mechanism of Cuyuxunxi Prescription in Promoting Wound Healing in Patients with Anal Fistula

Background

The healing process of the surgical wound of anal fistulotomy is much slower because of the presence of stool within the wound. Cuyuxunxi (CYXX) prescription is a Chinese herbal fumigant that is being used to wash surgical wound after anal fistulotomy. This study aimed at investigating the molecular mechanism of CYXX prescription using a network pharmacology-based strategy.

Materials and Methods

The active compounds in each herbal medicine were retrieved from the traditional Chinese medicine systems pharmacology (TCMSP) database and in Traditional Chinese Medicine Integrated Database (TCMID) analysis platform based on the criteria of oral bioavailability ≥40% and drug-likeness ≥0.2. The disease-related target genes were extracted from the Comparative Toxicogenomics Database. Protein-protein interaction network was built for the overlapped genes as well as functional enrichment analysis. Finally, an ingredient-target genes-pathway network was built by integrating all information.

Results

A total of 375 chemical ingredients of the 5 main herbal medicines in CYXX prescription were retrieved from TCMSP database and TCMID. Among the 375 chemical ingredients, 59 were active compounds. Besides, 325 target genes for 16 active compounds in 3 herbal medicines were obtained. Functional enrichment analysis revealed that these overlapped genes were significantly related with immune response, biosynthesis of antibiotics, and complement and coagulation cascades. A comprehensive network which contains 133 nodes (8 disease nodes, 3 drug nodes, 8 ingredients, 103 target gene nodes, 7 GO nodes, and 4 pathway nodes) was built.

Conclusion

The network built in this study might aid in understanding the action mechanism of CYXX prescription at molecular level to pathway level.

Association between promoter DNA methylation and gene expression in the pathogenesis of ischemic stroke

To assess DNA methylation sites as well as gene expression related to ischemic stroke (IS) and comprehensively reveal their correlation and possible pathological mechanisms, we implemented (1) genome-wide DNA methylation profiling from the GEO repository related to IS with and without symptoms; (2) identification of differentially methylation positions (DMPs) and genes (DMGs), functional enrichment analysis along with DMG regulatory network construction; (3) validation tests of 2 differential methylation positions of interest as well as analogous gene expression in other datasets and in IS patients and controls; and (4) correlation analysis of DNA methylation and mRNA expression data. In total, 870 DMPs were physically located within 693 DMGs. After disease ontology (DO), Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway, gene ontology (GO), protein-protein interaction (PPI) network construction as well as module analysis, HLA-DRB1 and HLA-DQB1 were identified. Their expression was validated in 4 other datasets but was significant in only 1, and the expression was lower in the IS group (P < 0.05). After validation in IS patients and controls, we found that these two genes showed more hypermethylation and lower expression levels in the IS group (P < 0.001). The methylation of genes was negatively associated with their expression (P < 0.05). The current study recognized a connection among DNA methylation and gene expression and emphasized the prominence of HLA-DRB1 and HLA-DQB1 in IS pathogenesis.

Global Transcriptomic Analysis of the Candida albicans Response to Treatment with a Novel Inhibitor of Filamentation

These results from whole-genome transcriptional profiling provide further insights into the biological activity and mode of action of a small-molecule inhibitor of C. albicans filamentation. This information will assist in the development of novel antivirulence strategies against C. albicans infections.

The opportunistic pathogenic fungus Candida albicans can cause devastating infections in immunocompromised patients. Its ability to undergo a morphogenetic transition from yeast to filamentous forms allows it to penetrate tissues and damage tissues, and the expression of genes associated with a number of pathogenetic mechanisms is also coordinately regulated with the yeast-to-hypha conversion. Therefore, it is widely considered that filamentation represents one of the main virulence factors of C. albicans. We have previously identified N-[3-(allyloxy)-phenyl]-4-methoxybenzamide (compound 9029936) as the lead compound in a series of small-molecule inhibitors of C. albicans filamentation and characterized its activity both in vitro and in vivo. This compound appears to be a promising candidate for the development of alternative antivirulence strategies for the treatment of C. albicans infections. In this study, we performed RNA sequencing analysis of samples obtained from C. albicans cells grown under filament-inducing conditions in the presence or absence of this compound. Overall, treatment with compound 9029936 resulted in 618 upregulated and 702 downregulated genes. Not surprisingly, some of the most downregulated genes included well-characterized genes associated with filamentation and virulence such as SAP5, ECE1 (candidalysin), and ALS3, as well as genes that impact metal chelation and utilization. Gene ontology analysis revealed an overrepresentation of cell adhesion, iron transport, filamentation, biofilm formation, and pathogenesis processes among the genes downregulated during treatment with this leading compound. Interestingly, the top upregulated genes suggested an enhancement of vesicular transport pathways, particularly those involving SNARE interactions.

IMPORTANCE These results from whole-genome transcriptional profiling provide further insights into the biological activity and mode of action of a small-molecule inhibitor of C. albicans filamentation. This information will assist in the development of novel antivirulence strategies against C. albicans infections.

Prognostic significance of TOP2A in non-small cell lung cancer revealed by bioinformatic analysis

Background

Lung cancer has been a common malignant tumor with a leading cause of morbidity and mortality, current molecular targets are woefully lacking comparing to the highly progressive cancer. The study is designed to identify new prognostic predictors and potential gene targets based on bioinformatic analysis of Gene Expression Omnibus (GEO) database.

Methods

Four cDNA expression profiles GSE19188, GSE101929, GSE18842 and GSE33532 were chosen from GEO database to analyze the differently expressed genes (DEGs) between non-small cell lung cancer (NSCLC) and normal lung tissues. After the DEGs functions were analyzed, the protein-protein interaction network (PPI) of DEGs were constructed, and the core gene in the network which has high connectivity degree with other genes was identified. We analyzed the association of the gene with the development of NSCLC as well as its prognosis. Lastly we explored the conceivable signaling mechanism of the gene regulation during the development of NSCLC.

Results

A total of 92 up regulated and 214 down regulated DEGs were shared in four cDNA expression profiles. Based on their PPI network, TOP2A was connected with most of other genes and was selected for further analysis. Kaplan-Meier overall survival analysis (OS) revealed that TOP2A was associated with worse NSCLC patients survival. And both GEPIA analysis and immunohistochemistry experiment (IHC) confirmed that TOP2A was aberrant gain of expression in cancer comparing to normal tissues. The clinical significance of TOP2A and probable signaling pathways it involved in were further explored, and a positive correlation between TOP2A and TPX2 expression was found in lung cancer tissues.

Conclusion

Using bioinformatic analysis, we revealed that TOP2A could be adopted as a prognostic indicator of NSCLC and it potentially regulate cancer development through co-work with TPX2. However, more detailed experiments are needed to clarify its drug target role in clinical medical use.

Phenotypical microRNA screen reveals a noncanonical role of CDK2 in regulating neutrophil migration

Neutrophil migration is essential for inflammatory responses to kill pathogens; however, excessive neutrophilic inflammation also leads to tissue injury and adverse effects. To discover novel therapeutic targets that modulate neutrophil migration, we performed a neutrophil-specific microRNA (miRNA) overexpression screen in zebrafish and identified 8 miRNAs as potent suppressors of neutrophil migration. Among those, miR-199 decreases neutrophil chemotaxis in zebrafish and human neutrophil-like cells. Intriguingly, in terminally differentiated neutrophils, miR-199 alters the cell cycle-related pathways and directly suppresses cyclin-dependent kinase 2 (Cdk2), whose known activity is restricted to cell cycle progression and cell differentiation. Inhibiting Cdk2, but not DNA replication, disrupts cell polarity and chemotaxis of zebrafish neutrophils without inducing cell death. Human neutrophil-like cells deficient in CDK2 fail to polarize and display altered signaling downstream of the formyl peptide receptor. Chemotaxis of primary human neutrophils is also reduced upon CDK2 inhibition. Furthermore, miR-199 overexpression or CDK2 inhibition significantly improves the outcome of lethal systemic inflammation challenges in zebrafish. Our results therefore reveal previously unknown functions of miR-199 and CDK2 in regulating neutrophil migration and provide directions in alleviating systemic inflammation.

Clinical, pathophysiologic, and genomic analysis of the outcomes of primary head and neck malignancy after pulmonary metastasectomy

The median overall survival (OS) of some head and neck malignancies, such as head and neck squamous cell carcinoma (HNSCC), with metastatic lesions was only 12 months. Whether aggressive pulmonary metastasectomy (PM) improves survival is controversial. Patients with primary head and neck malignancy undergoing PM were enrolled. Clinical outcomes were compared among different histological types. Whole-exome sequencing was used for matched pulmonary metastatic samples. The genes where genetic variants have been identified were sent for analysis by DAVID, IPA, and STRING. Forty-nine patients with primary head and neck malignancies were enrolled. Two-year postmetastasectomy survival (PMS) rates of adenoid cystic carcinoma, thyroid carcinoma, nasopharyngeal carcinoma, and HNSCC were 100%, 88.2%, 71.4%, and 59.2%, respectively (P = 0.024). In HNSCC, the time to distant metastasis was an independent predictive factor of the efficacy of PM. Several pathways, such as branched-chain amino acid (BCAA) consumption, were significantly associated with the progression of HNSCC [P < 0.001, fold enrichment (FE) = 5.45]. Moreover, metabolism-associated signaling pathways also seemed to be involved in cancer metastasis. Histological types and time to distant metastasis were important factors influencing the clinical outcomes of PM. For HNSCC, metabolic-associated signaling pathways were significantly associated with tumor progression and distant metastasis. Future validations are warranted.

Identification of Transcriptional Markers and microRNA-mRNA Regulatory Networks in Colon Cancer by Integrative Analysis of mRNA and microRNA Expression Profiles in Colon Tumor Stroma

The aberrant expression of microRNAs (miRNAs) and genes in tumor microenvironment (TME) has been associated with the pathogenesis of colon cancer. An integrative exploration of transcriptional markers (gene signatures) and miRNA–mRNA regulatory networks in colon tumor stroma (CTS) remains lacking. Using two datasets of mRNA and miRNA expression profiling in CTS, we identified differentially expressed miRNAs (DEmiRs) and differentially expressed genes (DEGs) between CTS and normal stroma. Furthermore, we identified the transcriptional markers which were both gene targets of DEmiRs and hub genes in the protein–protein interaction (PPI) network of DEGs. Moreover, we investigated the associations between the transcriptional markers and tumor immunity in colon cancer. We identified 17 upregulated and seven downregulated DEmiRs in CTS relative to normal stroma based on a miRNA expression profiling dataset. Pathway analysis revealed that the downregulated DEmiRs were significantly involved in 25 KEGG pathways (such as TGF-β, Wnt, cell adhesion molecules, and cytokine–cytokine receptor interaction), and the upregulated DEmiRs were involved in 10 pathways (such as extracellular matrix (ECM)-receptor interaction and proteoglycans in cancer). Moreover, we identified 460 DEGs in CTS versus normal stroma by a meta-analysis of two gene expression profiling datasets. Among them, eight upregulated DEGs were both hub genes in the PPI network of DEGs and target genes of the downregulated DEmiRs. We found that three of the eight DEGs were negative prognostic factors consistently in two colon cancer cohorts, including COL5A2, EDNRA, and OLR1. The identification of transcriptional markers and miRNA–mRNA regulatory networks in CTS may provide insights into the mechanism of tumor immune microenvironment regulation in colon cancer.

The Spectrinome: The Interactome of a Scaffold Protein Creating Nuclear and Cytoplasmic Connectivity and Function

We provide a review of Spectrin isoform function in the cytoplasm, the nucleus, the cell surface, and in intracellular signaling. We then discuss the importance of Spectrin’s E2/E3 chimeric ubiquitin conjugating and ligating activity in maintaining cellular homeostasis. Finally we present spectrin isoform subunit specific human diseases. We have created the Spectrinome, from the Human Proteome, Human Reactome and Human Atlas data and demonstrated how it can be a useful tool in visualizing and understanding spectrins myriad of cellular functions.

Impact statement

Spectrin was for the first 12 years after its discovery thought to be found only in erythrocytes. In 1981, Goodman and colleagues1found that spectrin-like molecules were ubiquitously found in non-erythroid cells leading to a great multitude of publications over the next thirty eight years. The discovery of multiple spectrin isoforms found associated with every cellular compartment, and representing 2-3% of cellular protein, has brought us to today’s understanding that spectrin is a scaffolding protein, with its own E2/E3 chimeric ubiquitin conjugating ligating activity that is involved in virtually every cellular function. We cover the history, localized functions of spectrin isoforms, human diseases caused by mutations, and provide the spectrinome: a useful tool for understanding the myriad of functions for one of the most important proteins in all eukaryotic cells.

Non-coding RNA regulatory networks

It is well established that the vast majority of human RNA transcripts do not encode for proteins and that non-coding RNAs regulate cell physiology and shape cellular functions. A subset of them is involved in gene regulation at different levels, from epigenetic gene silencing to post-transcriptional regulation of mRNA stability. Notably, the aberrant expression of many non-coding RNAs has been associated with aggressive pathologies. Rapid advances in network biology indicates that the robustness of cellular processes is the result of specific properties of biological networks such as scale-free degree distribution and hierarchical modularity, suggesting that regulatory network analyses could provide new insights on gene regulation and dysfunction mechanisms. In this study we present an overview of public repositories where non-coding RNA-regulatory interactions are collected and annotated, we discuss unresolved questions for data integration and we recall existing resources to build and analyse networks.

Identification of Aloperine as an anti-apoptotic Bcl2 protein inhibitor in glioma cells

OBJECTIVE

Aloperine (ALO), an alkaloid isolated from the leaves of Sophora alopecuroides, has been suggested to exhibit anti-inflammatory and anti-tumor properties and is traditionally used to treat various human diseases, including cancer. However, limited information is available about the mechanisms that determine the anti-tumor activities of ALO.

METHODS

Herein, through comprehensive bioinformatics methods and in vitro functional analyses, we evaluated the detailed anti-tumor mechanisms of ALO.

RESULTS

Using the databases Bioinformatics analysis tool for molecular mechanism of traditional Chinese medicine and PubChem Project, we identified the potential targets of ALO. A protein-protein interaction network was constructed to determine the relationship among these probable targets. Functional enrichment analysis revealed that ALO is potentially involved in the induction of apoptosis. In addition, molecular docking demonstrated that ALO expectedly docks into the active pocket of the Bcl2 protein, suggesting Bcl2 as a direct target of ALO. Moreover, western blot and qPCR analysis showed that ALO downregulated Bcl2 expression in human glioma cell lines, SK-N-AS and U118. Using flow cytometry methods, we further confirmed that ALO significantly promotes apoptosis in SK-N-AS and U118 cell lines, similar to the effect induced by ABT-737, a well-known Bcl2 inhibitor. In addition, Bcl-2 overexpression could rescue ALO-induced Bcl-2 inhibition and suppress pro-apoptotic effects in glioma cells.

CONCLUSION

Taken together, these findings suggest that the natural agent ALO effectively enhances apoptosis by acting as a potential Bcl2 inhibitor in human glioma cells.

Replication timing networks reveal a link between transcription regulatory circuits and replication timing control

DNA replication occurs in a defined temporal order known as the replication timing (RT) program and is regulated during development, coordinated with 3D genome organization and transcriptional activity. However, transcription and RT are not sufficiently coordinated to predict each other, suggesting an indirect relationship. Here, we exploit genome-wide RT profiles from 15 human cell types and intermediate differentiation stages derived from human embryonic stem cells to construct different types of RT regulatory networks. First, we constructed networks based on the coordinated RT changes during cell fate commitment to create highly complex RT networks composed of thousands of interactions that form specific functional subnetwork communities. We also constructed directional regulatory networks based on the order of RT changes within cell lineages, and identified master regulators of differentiation pathways. Finally, we explored relationships between RT networks and transcriptional regulatory networks (TRNs) by combining them into more complex circuitries of composite and bipartite networks. Results identified novel trans interactions linking transcription factors that are core to the regulatory circuitry of each cell type to RT changes occurring in those cell types. These core transcription factors were found to bind cooperatively to sites in the affected replication domains, providing provocative evidence that they constitute biologically significant directional interactions. Our findings suggest a regulatory link between the establishment of cell-type-specific TRNs and RT control during lineage specification.

Structures of the substrate-binding protein YfeA in apo and zinc-reconstituted holo forms

In the structural biology of bacterial substrate-binding proteins (SBPs), a growing number of comparisons between substrate-bound and substrate-free forms of metal atom-binding (cluster A-I) SBPs have revealed minimal structural differences between forms. These observations contrast with SBPs that bind substrates such as amino acids or nucleic acids and may undergo >60° rigid-body rotations. Substrate transfer in these SBPs is described by a Venus flytrap model, although this model may not apply to all SBPs. In this report, structures are presented of substrate-free (apo) and reconstituted substrate-bound (holo) YfeA, a polyspecific cluster A-I SBP from Yersinia pestis. It is demonstrated that an apo cluster A-I SBP can be purified by fractionation when co-expressed with its cognate transporter, adding an alternative strategy to the mutagenesis or biochemical treatment used to generate other apo cluster A-I SBPs. The apo YfeA structure contains 111 disordered protein atoms in a mobile helix located in the flexible carboxy-terminal lobe. Metal binding triggers a 15-fold reduction in the solvent-accessible surface area of the metal-binding site and reordering of the 111 protein atoms in the mobile helix. The flexible lobe undergoes a 13.6° rigid-body rotation that is driven by a spring-hammer metal-binding mechanism. This asymmetric rigid-body rotation may be unique to metal atom-binding SBPs (i.e. clusters A-I, A-II and D-IV).

Autograft microskin combined with adipose-derived stem cell enhances wound healing in a full-thickness skin defect mouse model

OBJECTIVE

Autograft microskin transplantation has been widely used as a skin graft therapy in full-thickness skin defect. However, skin grafting failure can lead to a pathological delay wound healing due to a poor vascularization bed. Considering the active role of adipose-derived stem cell (ADSC) in promoting angiogenesis, we intend to investigate the efficacy of autograft microskin combined with ADSC transplantation for facilitating wound healing in a full-thickness skin defect mouse model.

MATERIAL AND METHODS

An in vivo full-thickness skin defect mouse model was used to evaluate the contribution of transplantation microskin and ADSC in wound healing. The angiogenesis was detected by immunohistochemistry staining. In vitro paracrine signaling pathway was evaluated by protein array and Gene Ontology, Kyoto Encyclopedia of Genes and Genomes pathway, and protein-protein interaction network analysis.

RESULTS

Co-transplantation of microskin and ADSC potentiated the wound healing with better epithelization, smaller scar thickness, and higher angiogenesis (CD31) in the subcutaneous layer. We found both EGF and VEGF cytokines were secreted by microskin in vitro. Additionally, secretome proteomic analysis in a co-culture system of microskin and ADSC revealed that ADSC could secrete a wide range of important molecules to form a reacting network with microskin, including VEGF, IL-6, EGF, uPAR, MCP-3, G-CSF, and Tie-2, which most likely supported the angiogenesis effect as observed.

CONCLUSION

Overall, we concluded that the use of ADSC partially modulates microskin function and enhances wound healing by promoting angiogenesis in a full-thickness skin defect mouse model.

Polyamine Synthesis Effects Capsule Expression by Reduction of Precursors in Streptococcus pneumoniae

Streptococcus pneumoniae (pneumococcus, Spn) colonizes the human nasopharynx asymptomatically but can cause infections such as otitis media, and invasive pneumococcal disease such as community-acquired pneumonia, meningitis, and sepsis. Although the success of Spn as a pathogen can be attributed to its ability to synthesize and regulate capsular polysaccharide (CPS) for survival in the host, the mechanisms of CPS regulation are not well-described. Recent studies from our lab demonstrate that deletion of a putative polyamine biosynthesis gene (ΔcadA) in Spn TIGR4 results in the loss of the capsule. In this study, we characterized the transcriptome and metabolome of ΔcadA and identified specific mechanisms that could explain the regulatory role of polyamines in pneumococcal CPS biosynthesis. Our data indicate that impaired polyamine synthesis impacts galactose to glucose interconversion via the Leloir pathway which limits the availability of UDP-galactose, a precursor of serotype 4 CPS, and UDP-N-acetylglucosamine (UDP-GlcNAc), a nucleotide sugar precursor that is at the intersection of CPS and peptidoglycan repeat unit biosynthesis. Reduced carbon flux through glycolysis, coupled with altered fate of glycolytic intermediates further supports impaired synthesis of UDP-GlcNAc. A significant increase in the expression of transketolases indicates a potential shift in carbon flow toward the pentose phosphate pathway (PPP). Higher PPP activity could constitute oxidative stress responses in ΔcadA which warrants further investigation. The results from this study clearly demonstrate the potential of polyamine synthesis, targeted for cancer therapy in human medicine, for the development of novel prophylactic and therapeutic strategies for treating bacterial infections.

Differential impact of classical and non-canonical NF-κB pathway-related gene expression on the survival of breast cancer patients

Inflammation is a well-known driver of carcinogenesis and cancer progression, often attributed to the tumor microenvironment. However, tumor cells themselves are capable of secreting a variety of inflammatory molecules, leading to the activation of specific signaling pathways that promote tumor progression. The NF-κB signaling pathway is one of the most important connections between inflammation and tumorigenesis. NF-κB is a superfamily of transcription factors that plays an important role in several types of hematological and solid tumors, including breast cancer. However, the role of the NF-κB pathway in the survival of breast cancer patients is poorly studied. In this study, we analyzed and related the expression of both canonical and alternative NF-κB pathways and selected target genes with the relapse-free and overall survival of breast cancer patients. We used the public database Kaplan-Meier plotter (KMplot) which includes gene expression data and survival information of 3951 breast cancer patients. We found that the expression of IKKα was associated with poor relapse-free survival in patients with ER-positive tumors. Moreover, the expression of IL-8 and MMP-1 was associated with poor relapse-free and overall survival. In contrast, expression of IKKβ, p50, and p65 from the canonical pathway, and NIK and RELB from the alternative pathway correlated with better relapse-free survival also when the patients were classified by their hormonal and nodal status. Our study suggests that the expression of genes of the canonical and alternative NF-κB pathways is ultimately critical for tumor persistence. Understanding the communication between both pathways would help to find better therapeutic and prophylactic targets to prevent breast cancer progression and relapse.

Molecular characterization and computational structure prediction of activin receptor type IIB in aseel and broiler chicken

The present study was formulated to characterize and comprehend the molecular structural characteristics of ACTRIIB receptor in Aseel and control broiler (CB) populations. The full length coding sequence (1539 bp) of the receptor was amplified, cloned, sequenced and analyzed using bioinformatic tools. The physico chemical properties of protein and structural features like secondary structure, solvent accessibility and disorder regions were computed. The 3D structure was predicted by I-TASSER and evaluated by Ramachandran Plot and tools under Structural Analysis and Verification Server. The nucleotides differences between CB and Aseel were c. [156G > A; 210 T > C; 493C > T; c.520G > C; 665A > C; 686G > A; 937C > G; 1011A > C; 1130A > G; 1208 T > A; 1326 T > C; 1433 T > C]. The amino acid substitutions between CB and Aseel were p. [(Pro165Ser; Glu174Gln; Gln222Pro; Ser229Asn; His313Asp; Gln377Arg; Val403Asp; and Ile478Thr)]. While, the silent changes includes p. [(Lys53=; Glu71=; Leu337=; Asp442=)]. The molecular weight of mature protein was predicted to be 55.51 kDa and 57.80 kDa in Aseel and CB, respectively. The higher rank 3D model had a C-score of -1.60 in Aseel and - 1.41 in CB, while the estimated TM-score (0.54 ± 0.14) and RMSD (5.8 ± 1.2 Å) were found to be similar in Aseel and CB. Among the 512 residues, >90% were in favored region, 4.7% in allowed region and <1.5% in disallowed region in both Aseel and CB. The pattern of contact map was comparable in Aseel and CB. The Hydrogen bond plots of the Aseel and CB shared similar secondary structure pattern. The ACTRIIB protein was predicted to interact with ACVR1B, ACVR1C, INHBA, SMAD 1,2,5,7 & 9 and BMPR1A&B. Clustal and phylogenetic analysis implied that both the lines were closely related and formed a sub cluster with in avian cluster. The current research provides insights about structural and functional aspects of the receptor and also aids in understanding the evolutionary history of ACTRIIB.

Long non-coding RNA small nucleolar RNA host gene 7 is upregulated and promotes cell proliferation in thyroid cancer

Thyroid cancer (THCA) is one of the most common types of endocrine cancer worldwide. However, the mechanisms underlying THCA progression have not been fully elucidated. Recent studies have demonstrated that long non-coding RNAs (lncRNAs) are dysregulated in human diseases, and are involved in regulating various biological processes. Furthermore, several reports have indicated that lncRNAs serve important roles in THCA. In the present study, a dataset from The Cancer Genome Atlas was used to analyze the expression levels and the clinical information of small nucleolar RNA host gene 7 (SNHG7) in THCA. Starbase was used to construct the competing endogenous RNA network. The Molecule Annotation System was used to analyze the data from Gene Ontology and Kyoto Encyclopedia of Genes and Genomes databases. Furthermore, cell proliferation and cell cycle assays were used to detect the functions of SNHG7 in THCA. The present study revealed for the first time, to the best of our knowledge, that SNHG7 is markedly upregulated in THCA samples following analysis of The Cancer Genome Atlas datasets. SNHG7 expression was higher in advanced stage compared with early stage THCA samples. In addition, high expression levels of SNHG7 were associated with shorter survival times in THCA patients compared with low expression levels. Bioinformatics analysis revealed that SNHG7 was associated with the processes of ‘protein translation’, ‘viral life cycle’, ‘RNA processing’, ‘mRNA splicing’, ‘histone ubiquitination’, ‘endoplasmic reticulum-to-Golgi vesicle-mediated transport’, ‘sister chromatid cohesion’, ‘DNA damage checkpoint regulation’, ‘translation’ and ‘the spliceosome’. Additionally, knockdown of SNHG7 significantly inhibited thyroid cancer cell proliferation and cell cycle progression in vitro. Taken together, the results obtained in the present study suggested that SNHG7 may serve as a novel therapeutic and prognostic target for THCA.

Profiling the lncRNA-miRNA-mRNA ceRNA network to reveal potential crosstalk between inflammatory bowel disease and colorectal cancer

Background

Because of the increasing dysplasia rate in the lifelong course of inflammatory bowel disease (IBD) patients, it is imperative to characterize the crosstalk between IBD and colorectal cancer (CRC). However, there have been no reports revealing the occurrence of the ceRNA network in IBD-related CRC.

Methods

In this study, we conducted gene expression profile studies of databases and performed an integrated analysis to detect the potential of lncRNA-miRNA-mRNA ceRNA in regulating disease transformation. R packages were used to screen differentially expressed mRNA, lncRNA and miRNA among CRC, IBD and normal tissue. The lncRNA-miRNA-mRNA network was constructed based on predicted miRNA-targeted lncRNAs and miRNA-targeted mRNAs. Functional analyses were then conducted to identify genes involved in the ceRNA network, and key lncRNAs were evaluated based on several clinical outcomes.

Results

A total of three lncRNAs, 15 miRNAs, and 138 mRNAs were identified as potential mediators in the pathophysiological processes of IBD-related CRC. Gene Ontology annotation enrichment analysis confirmed that the dysplasia process was strongly associated with immune response, response to lipopolysaccharide, and inflammatory response. Survival analysis showed that LINC01106 (HR = 1.7; p < 0.05) were strongly associated with overall survival of colorectal cancer patients. The current study identified a series of IBD-related mRNAs, miRNA, and lncRNAs, and highlighted the important role of ceRNAs in the pathogenesis of IBD-related CRC. Among them, the LINC01106-miRNA-mRNA axis was identified as vital targets for further research.