Next-Generation Grade and Survival Expression Biomarkers of Human Gliomas Based on Algorithmically Reconstructed Molecular Pathways

In gliomas, expression of certain marker genes is strongly associated with survival and tumor type and often exceeds histological assessments. Using a human interactome model, we algorithmically reconstructed 7494 new-type molecular pathways that are centered each on an individual protein. Each single-gene expression and gene-centric pathway activation was tested as a survival and tumor grade biomarker in gliomas and their diagnostic subgroups (IDH mutant or wild type, IDH mutant with 1p/19q co-deletion, MGMT promoter methylated or unmethylated), including the three major molecular subtypes of glioblastoma (proneural, mesenchymal, classical). We used three datasets from The Cancer Genome Atlas and the Chinese Glioma Genome Atlas, which in total include 527 glioblastoma and 1097 low grade glioma profiles. We identified 2724 such gene and 2418 pathway survival biomarkers out of total 17,717 genes and 7494 pathways analyzed. We then assessed tumor grade and molecular subtype biomarkers and with the threshold of AUC > 0.7 identified 1322/982 gene biomarkers and 472/537 pathway biomarkers. This suggests roughly two times greater efficacy of the reconstructed pathway approach compared to gene biomarkers. Thus, we conclude that activation levels of algorithmically reconstructed gene-centric pathways are a potent class of new-generation diagnostic and prognostic biomarkers for gliomas.

Disparity between Inter-Patient Molecular Heterogeneity and Repertoires of Target Drugs Used for Different Types of Cancer in Clinical Oncology

Inter-patient molecular heterogeneity is the major declared driver of an expanding variety of anticancer drugs and personalizing their prescriptions. Here, we compared interpatient molecular heterogeneities of tumors and repertoires of drugs or their molecular targets currently in use in clinical oncology. We estimated molecular heterogeneity using genomic (whole exome sequencing) and transcriptomic (RNA sequencing) data for 4890 tumors taken from The Cancer Genome Atlas database. For thirteen major cancer types, we compared heterogeneities at the levels of mutations and gene expression with the repertoires of targeted therapeutics and their molecular targets accepted by the current guidelines in oncology. Totally, 85 drugs were investigated, collectively covering 82 individual molecular targets. For the first time, we showed that the repertoires of molecular targets of accepted drugs did not correlate with molecular heterogeneities of different cancer types. On the other hand, we found that the clinical recommendations for the available cancer drugs were strongly congruent with the gene expression but not gene mutation patterns. We detected the best match among the drugs usage recommendations and molecular patterns for the kidney, stomach, bladder, ovarian and endometrial cancers. In contrast, brain tumors, prostate and colorectal cancers showed the lowest match. These findings provide a theoretical basis for reconsidering usage of targeted therapeutics and intensifying drug repurposing efforts.

Pathway Based Analysis of Mutation Data Is Efficient for Scoring Target Cancer Drugs

Despite the significant achievements in chemotherapy, cancer remains one of the leading causes of death. Target therapy revolutionized this field, but efficiencies of target drugs show dramatic variation among individual patients. Personalization of target therapies remains, therefore, a challenge in oncology. Here, we proposed molecular pathway-based algorithm for scoring of target drugs using high throughput mutation data to personalize their clinical efficacies. This algorithm was validated on 3,800 exome mutation profiles from The Cancer Genome Atlas (TCGA) project for 128 target drugs. The output values termed Mutational Drug Scores (MDS) showed positive correlation with the published drug efficiencies in clinical trials. We also used MDS approach to simulate all known protein coding genes as the putative drug targets. The model used was built on the basis of 18,273 mutation profiles from COSMIC database for eight cancer types. We found that the MDS algorithm-predicted hits frequently coincide with those already used as targets of the existing cancer drugs, but several novel candidates can be considered promising for further developments. Our results evidence that the MDS is applicable to ranking of anticancer drugs and can be applied for the identification of novel molecular targets.

Pathway Instability Is an Effective New Mutation-Based Type of Cancer Biomarkers

DNA mutations play a crucial role in cancer development and progression. Mutation profiles vary dramatically in different cancer types and between individual tumors. Mutations of several individual genes are known as reliable cancer biomarkers, although the number of such genes is tiny and does not enable differential diagnostics for most of the cancers. We report here a technique enabling dramatically increased efficiency of cancer biomarkers development using DNA mutations data. It includes a quantitative metric termed Pathway instability (PI) based on mutations enrichment of intracellular molecular pathways. This method was tested on 5,956 tumor mutation profiles of 15 cancer types from The Cancer Genome Atlas (TCGA) project. Totally, we screened 2,316,670 mutations in 19,872 genes and 1,748 molecular pathways. Our results demonstrated considerable advantage of pathway-based mutation biomarkers over individual gene mutation profiles, as reflected by more than two orders of magnitude greater numbers by high-quality [ROC area-under-curve (AUC)>0.75] biomarkers. For example, the number of such high-quality mutational biomarkers distinguishing between different cancer types was only six for the individual gene mutations, and already 660 for the pathway-based biomarkers. These results evidence that PI value can be used as a new generation of complex cancer biomarkers significantly outperforming the existing gene mutation biomarkers.

Cytotoxic Potential of the Novel Horseshoe Crab Peptide Polyphemusin III

Biological activity of the new antimicrobial peptide polyphemusin III from the horseshoe crab Limulus polyphemus was examined against bacterial strains and human cancer, transformed, and normal cell cultures. Polyphemusin III has the amino acid sequence RRGCFRVCYRGFCFQRCR and is homologous to other β-hairpin peptides from the horseshoe crab. Antimicrobial activity of the peptide was evaluated and MIC (minimal inhibitory concentration) values were determined. IC₅₀ (half-maximal inhibitory concentration) values measured toward human cells revealed that polyphemusin III showed a potent cytotoxic activity at concentrations of <10 μM. Polyphemusin III caused fast permeabilization of the cytoplasmic membrane of human leukemia cells HL-60, which was measured with trypan blue exclusion assay and lactate dehydrogenase-release assay. Flow cytometry experiments for annexin V-FITC/ propidium iodide double staining revealed that the caspase inhibitor, Z-VAD-FMK, did not abrogate disruption of the plasma membrane by polyphemusin III. Our data suggest that polyphemusin III disrupts the plasma membrane integrity and induces cell death that is apparently not related to apoptosis. In comparison to known polyphemusins and tachyplesins, polyphemusin III demonstrates a similar or lower antimicrobial effect, but significantly higher cytotoxicity against human cancer and transformed cells in vitro.

The contribution of Saccharomyces cerevisiae replicative age to the variations in the levels of Trx2p, Pdr5p, Can1p and Idh isoforms

Asymmetrical division can be a reason for microbial populations heterogeneity. In particular, budding yeast daughter cells are more vulnerable to stresses than the mothers. It was suggested that yeast mother cells could also differ from each other depending on their replicative age. To test this, we measured the levels of Idh1-GFP, Idh2-GFP, Trx2-GFP, Pdr5-GFP and Can1-GFP proteins in cells of the few first, most represented, age cohorts. Pdr5p and Can1p were selected because of the pronounced mother-bud asymmetry for these proteins distributions, Trx2p as indicator of oxidative stress. Isocitrate dehydrogenase subunits Idh1p and Idh2p were assessed because their levels are regulated by mitochondria. We found a small negative correlation between yeast replicative age and Idh1-GFP or Idh2-GFP but not Trx2-GFP levels. Mitochondrial network fragmentation was also confirmed as an early event of replicative aging. No significant difference in the membrane proteins levels Pdr5p and Can1p was found. Moreover, the elder mother cells showed lower coefficient of variation for Pdr5p levels compared to the younger ones and the daughters. Our data suggest that the levels of stress-response proteins Pdr5p and Trx2p in the mother cells are stable during the first few cell cycles regardless of their mother-bud asymmetry.

Early stage of cytomegalovirus infection suppresses host microRNA expression regulation in human fibroblasts

Responses to human cytomegalovirus (HCMV) infection are largely individual and cell type specific. We investigated molecular profiles in 2 primary cell cultures of human fibroblasts, which are highly or marginally sensitive to HCMV infection, respectively. We screened expression of genes and microRNAs (miRs) at the early (3 hours) stage of infection. To assess molecular pathway activation profiles, we applied bioinformatic algorithms OncoFinder and MiRImpact. In both cell types, pathway regulation properties at mRNA and miR levels were markedly different. Surprisingly, in the infected highly sensitive cells, we observed a "freeze" of miR expression profiles compared to uninfected controls. Our results evidence that in the sensitive cells, HCMV blocks intracellular regulation of microRNA expression already at the earliest stage of infection. These data suggest somewhat new functions for HCMV products and demonstrate dependence of miR expression arrest on the host-encoded factors.

Mitochondrial signaling in Saccharomyces cerevisiae pseudohyphae formation induced by butanol

Yeasts growing limited for nitrogen source or treated with fusel alcohols form elongated cells--pseudohyphae. Absence of mitochondrial DNA or anaerobic conditions inhibits this process, but the precise role of mitochondria is not clear. We found that a significant percentage of pseudohyphal cells contained mitochondria with different levels of membrane potential within one cell. An uncoupler carbonyl cyanide p-(trifluoromethoxy) phenylhydrazone (FCCP), but not the ATP-synthase inhibitor oligomycin D, prevented pseudohyphal growth. Interestingly, repression of the MIH1 gene encoding phosphatase activator of the G2/M transition partially restores the ability of yeast to form pseudohyphal cells in the presence of FCCP or in the absence of mitochondrial DNA. At the same time, retrograde signaling (the one triggered by dysfunctional mitochondria) appeared to be a positive regulator of butanol-induced pseudohyphae formation: the deletion of any of the retrograde signaling genes (RTG1, RTG2, or RTG3) partially suppressed pseudohyphal growth. Together, our data suggest that two subpopulations of mitochondria are required for filamentous growth: one with high and another with low transmembrane potential. These mitochondria-activated signaling pathways appear to converge at Mih1p level.