Identification of Common Prognostic Gene Expression Signatures with Biological Meanings from Microarray Gene Expression Datasets

Numerous prognostic gene expression signatures for breast cancer were generated previously with few overlap and limited insight into the biology of the disease. Here we introduce a novel algorithm named SCoR (Survival analysis using Cox proportional hazard regression and Random resampling) to apply random resampling and clustering methods in identifying gene features correlated with time to event data. This is shown to reduce overfitting noises involved in microarray data analysis and discover functional gene sets linked to patient survival. SCoR independently identified a common poor prognostic signature composed of cell proliferation genes from six out of eight breast cancer datasets. Furthermore, a sequential SCoR analysis on highly proliferative breast cancers repeatedly identified T/B cell markers as favorable prognosis factors. In glioblastoma, SCoR identified a common good prognostic signature of chromosome 10 genes from two gene expression datasets (TCGA and REMBRANDT), recapitulating the fact that loss of one copy of chromosome 10 (which harbors the tumor suppressor PTEN) is linked to poor survival in glioblastoma patients. SCoR also identified prognostic genes on sex chromosomes in lung adenocarcinomas, suggesting patient gender might be used to predict outcome in this disease. These results demonstrate the power of SCoR to identify common and biologically meaningful prognostic gene expression signatures.     

A Fungal Sarcolemmal Membrane-Associated Protein (SLMAP) Homolog Plays a Fundamental Role in Development and Localizes to the Nuclear Envelope,Endoplasmic Reticulum,and Mitochondria

Sarcolemmal membrane-associated protein (SLMAP) is a tail-anchored protein involved in fundamental cellular processes, such as myoblast fusion, cell cycle progression, and chromosomal inheritance. Further, SLMAP misexpression is associated with endothelial dysfunctions in diabetes and cancer. SLMAP is part of the conserved striatin-interacting phosphatase and kinase (STRIPAK) complex required for specific signaling pathways in yeasts, filamentous fungi, insects, and mammals. In filamentous fungi, STRIPAK was initially discovered in Sordaria macrospora, a model system for fungal differentiation. Here, we functionally characterize the STRIPAK subunit PRO45, a homolog of human SLMAP. We show that PRO45 is required for sexual propagation and cell-to-cell fusion and that its forkhead-associated (FHA) domain is essential for these processes. Protein-protein interaction studies revealed that PRO45 binds to STRIPAK subunits PRO11 and SmMOB3, which are also required for sexual propagation. Superresolution structured-illumination microscopy (SIM) further established that PRO45 localizes to the nuclear envelope, endoplasmic reticulum, and mitochondria. SIM also showed that localization to the nuclear envelope requires STRIPAK subunits PRO11 and PRO22, whereas for mitochondria it does not. Taken together, our study provides important insights into fundamental roles of the fungal SLMAP homolog PRO45 and suggests STRIPAK-related and STRIPAK-unrelated functions.     

RADIA: RNA and DNA Integrated Analysis for Somatic Mutation Detection

The detection of somatic single nucleotide variants is a crucial component to the characterization of the cancer genome. Mutation calling algorithms thus far have focused on comparing the normal and tumor genomes from the same individual. In recent years, it has become routine for projects like The Cancer Genome Atlas (TCGA) to also sequence the tumor RNA. Here we present RADIA (RNA and DNA Integrated Analysis), a novel computational method combining the patient-matched normal and tumor DNA with the tumor RNA to detect somatic mutations. The inclusion of the RNA increases the power to detect somatic mutations, especially at low DNA allelic frequencies. By integrating an individual’s DNA and RNA, we are able to detect mutations that would otherwise be missed by traditional algorithms that examine only the DNA. We demonstrate high sensitivity (84%) and very high precision (98% and 99%) for RADIA in patient data from endometrial carcinoma and lung adenocarcinoma from TCGA. Mutations with both high DNA and RNA read support have the highest validation rate of over 99%. We also introduce a simulation package that spikes in artificial mutations to patient data, rather than simulating sequencing data from a reference genome. We evaluate sensitivity on the simulation data and demonstrate our ability to rescue back mutations at low DNA allelic frequencies by including the RNA. Finally, we highlight mutations in important cancer genes that were rescued due to the incorporation of the RNA.     

Pgrmc1 (Progesterone Receptor Membrane Component 1) Associates with Epidermal Growth Factor Receptor and Regulates Erlotinib Sensitivity

Tumorigenesis requires the concerted action of multiple pathways, including pathways that stimulate proliferation and metabolism. Epidermal growth factor receptor (EGFR) is a transmembrane receptor-tyrosine kinase that is associated with cancer progression, and the EGFR inhibitors erlotinib/tarceva and tyrphostin/AG-1478 are potent anti-cancer therapeutics. Pgrmc1 (progesterone receptor membrane component 1) is a cytochrome b₅-related protein that is up-regulated in tumors and promotes cancer growth. Pgrmc1 and its homologues have been implicated in cell signaling, and we show here that Pgrmc1 increases susceptibility to AG-1478 and erlotinib, increases plasma membrane EGFR levels, and co-precipitates with EGFR. Pgrmc1 co-localizes with EGFR in cytoplasmic vesicles and co-fractionates with EGFR in high density microsomes. The findings have therapeutic potential because a Pgrmc1 small molecule ligand, which inhibits growth in a variety of cancer cell types, de-stabilized EGFR in multiple tumor cell lines. EGFR is one of the most potent receptor-tyrosine kinases driving tumorigenesis, and our data support a role for Pgrmc1 in promoting several cancer phenotypes at least in part by binding EGFR and stabilizing plasma membrane pools of the receptor.     

The molecular complexity of mu and pi symbiont DNA of Paramecium aurelia

The molecular size of mu and pi symbionts of Parameciumaurelia has been calculated from renaturation kinetic data. Observed values were 0.78 × 10⁹ daltons for mu particle DNA and 0.81 × 10⁹ daltons for pi particle DNA. Estimates of analytical complexity were 4.45 × 10⁹ and 5.05 × 10⁹ daltons respectively. Based on these data, mu and pi symbionts appear to possess multiple genomes and contain a minimum of 5 or 6 copies of each DNA sequence.     

Yu Ping Feng San,an Ancient Chinese Herbal Decoction Containing Astragali Radix,Atractylodis Macrocephalae Rhizoma and Saposhnikoviae Radix,Regulates the Release of Cytokines in Murine Macrophages

Yu Ping Feng San (YPFS), a Chinese herbal decoction, is composed of Astragali Radix (AR; Huangqi), Atractylodis Macrocephalae Rhizoma (AMR; Baizhu) and Saposhnikoviae Radix (SR; Fangfeng) in a weight ratio of 1∶2∶1. Clinically, YPFS has been widely used to regulate immune functions; however, the action mechanism of it is not known. Here, we addressed this issue by providing detail analyses of chemical and biological properties of YPFS. By using rapid resolution liquid chromatography coupled with mass spectrometry, fifteen chemicals deriving from different herbs of YPFS were determined, and which served as a control for the standardization of the herbal extract of YPFS. In general, the amounts of chosen chemical markers were higher in a preparation of YPFS as compared to that of single herb or two-herb compositions. In order to reveal the immune functions of YPFS, the standardized extract was applied onto cultured murine macrophages. The treatment of YPFS stimulated the mRNA and protein expressions of pro-inflammatory cytokines via activation of NF-κB by enhancing IκBα degradation. In contrast, the application of YPFS suppressed the expressions of pro-inflammatory cytokines significantly in the lipopolysaccharide (LPS)-induced chronic inflammation model. In addition, YPFS could up regulate the phagocytic activity in cultured macrophages. These results therefore supported the bi-directional immune-modulatory roles of YPFS in regulating the releases of cytokines from macrophages.     

Evolution of the B3 DNA Binding Superfamily: New Insights into REM Family Gene Diversification

Background

The B3 DNA binding domain includes five families: auxin response factor (ARF), abscisic acid-insensitive3 (ABI3), high level expression of sugar inducible (HSI), related to ABI3/VP1 (RAV) and reproductive meristem (REM). The release of the complete genomes of the angiosperm eudicots Arabidopsis thaliana and Populus trichocarpa, the monocot Orysa sativa, the bryophyte Physcomitrella patens,the green algae Chlamydomonas reinhardtii and Volvox carteri and the red algae Cyanidioschyzon melorae provided an exceptional opportunity to study the evolution of this superfamily.

Methodology

In order to better understand the origin and the diversification of B3 domains in plants, we combined comparative phylogenetic analysis with exon/intron structure and duplication events. In addition, we investigated the conservation and divergence of the B3 domain during the origin and evolution of each family.

Conclusions

Our data indicate that showed that the B3 containing genes have undergone extensive duplication events, and that the REM family B3 domain has a highly diverged DNA binding. Our results also indicate that the founding member of the B3 gene family is likely to be similar to the ABI3/HSI genes found in C. reinhardtii and V. carteri. Among the B3 families, ABI3, HSI, RAV and ARF are most structurally conserved, whereas the REM family has experienced a rapid divergence. These results are discussed in light of their functional and evolutionary roles in plant development.     

The Linker for Activation of B Cells (LAB)/Non-T Cell Activation Linker (NTAL) Regulates Triggering Receptor Expressed on Myeloid Cells (TREM)-2 Signaling and Macrophage Inflammatory Responses Independently of the Linker for Activation of T Cells

Triggering receptor expressed on myeloid cells-2 (TREM-2) is rapidly emerging as a key regulator of the innate immune response via its regulation of macrophage inflammatory responses. Here we demonstrate that proximal TREM-2 signaling parallels other DAP12-based receptor systems in its use of Syk and Src-family kinases. However, we find that the linker for activation of T cells (LAT) is severely reduced as monocytes differentiate into macrophages and that TREM-2 exclusively uses the linker for activation of B cells (LAB encoded by the gene Lat2^−/−) to mediate downstream signaling. LAB is required for TREM-2-mediated activation of Erk1/2 and dampens proximal TREM-2 signals through a novel LAT-independent mechanism resulting in macrophages with proinflammatory properties. Thus, Lat2^−/− macrophages have increased TREM-2-induced proximal phosphorylation, and lipopolysaccharide stimulation of these cells leads to increased interleukin-10 (IL-10) and decreased IL-12p40 production relative to wild type cells. Together these data identify LAB as a critical, LAT-independent regulator of TREM-2 signaling and macrophage development capable of controlling subsequent inflammatory responses.     

Origin and Evolution of Sulfadoxine Resistant Plasmodium falciparum

The Thailand-Cambodia border is the epicenter for drug-resistant falciparum malaria. Previous studies have shown that chloroquine (CQ) and pyrimethamine resistance originated in this region and eventually spread to other Asian countries and Africa. However, there is a dearth in understanding the origin and evolution of dhps alleles associated with sulfadoxine resistance. The present study was designed to reveal the origin(s) of sulfadoxine resistance in Cambodia and its evolutionary relationship to African and South American dhps alleles. We sequenced 234 Cambodian Plasmodium falciparum isolates for the dhps codons S436A/F, A437G, K540E, A581G and A613S/T implicated in sulfadoxine resistance. We also genotyped 10 microsatellite loci around dhps to determine the genetic backgrounds of various alleles and compared them with the backgrounds of alleles prevalent in Africa and South America. In addition to previously known highly-resistant triple mutant dhps alleles SGEGA and AGEAA (codons 436, 437, 540, 581, 613 are sequentially indicated), a large proportion of the isolates (19.3%) contained a 540N mutation in association with 437G/581G yielding a previously unreported triple mutant allele, SGNGA. Microsatellite data strongly suggest the strength of selection was greater on triple mutant dhps alleles followed by the double and single mutants. We provide evidence for at least three independent origins for the double mutants, one each for the SGKGA, AGKAA and SGEAA alleles. Our data suggest that the triple mutant allele SGEGA and the novel allele SGNGA have common origin on the SGKGA background, whereas the AGEAA triple mutant was derived from AGKAA on multiple, albeit limited, genetic backgrounds. The SGEAA did not share haplotypes with any of the triple mutants. Comparative analysis of the microsatellite haplotypes flanking dhps alleles from Cambodia, Kenya, Cameroon and Venezuela revealed an independent origin of sulfadoxine resistant alleles in each of these regions.     

HEK293T Cells Are Heterozygous for CCR5 Delta 32 Mutation

C-C chemokine receptor 5 (CCR5) is a receptor for chemokines and a co-receptor for HIV-1 entry into the target CD4+ cells. CCR5 delta 32 deletion is a loss-of-function mutation, resistant to HIV-1 infection. We tried to induce the CCR5 delta 32 mutation harnessing the genome editing technique, CRISPR-Cas9 (Clustered Regularly Interspaced Short Palindromic Repeats, CRISPR and CRISPR associated protein 9, Cas9) in the commonly used cell line human embryonic kidney HEK 293T cells. Surprisingly, we found that HEK293T cells are heterozygous for CCR5 delta 32 mutation, in contrast to the wild type CCR5 cells, human acute T cell leukemia cell line Jurkat and human breast adenocarcinoma cell line MDA-MB-231 cells. This finding indicates that at least one human cell line is heterozygous for the CCR5 delta 32 mutation. We also found that in PCR amplification, wild type CCR5 DNA and mutant delta 32 DNA can form mismatched heteroduplex and move slowly in gel electrophoresis.     

PASNet: pathway-associated sparse deep neural network for prognosis prediction from high-throughput data

Background

Predicting prognosis in patients from large-scale genomic data is a fundamentally challenging problem in genomic medicine. However, the prognosis still remains poor in many diseases. The poor prognosis may be caused by high complexity of biological systems, where multiple biological components and their hierarchical relationships are involved. Moreover, it is challenging to develop robust computational solutions with high-dimension, low-sample size data.

Results

In this study, we propose a Pathway-Associated Sparse Deep Neural Network (PASNet) that not only predicts patients’ prognoses but also describes complex biological processes regarding biological pathways for prognosis. PASNet models a multilayered, hierarchical biological system of genes and pathways to predict clinical outcomes by leveraging deep learning. The sparse solution of PASNet provides the capability of model interpretability that most conventional fully-connected neural networks lack. We applied PASNet for long-term survival prediction in Glioblastoma multiforme (GBM), which is a primary brain cancer that shows poor prognostic performance. The predictive performance of PASNet was evaluated with multiple cross-validation experiments. PASNet showed a higher Area Under the Curve (AUC) and F1-score than previous long-term survival prediction classifiers, and the significance of PASNet’s performance was assessed by Wilcoxon signed-rank test. Furthermore, the biological pathways, found in PASNet, were referred to as significant pathways in GBM in previous biology and medicine research.

Conclusions

PASNet can describe the different biological systems of clinical outcomes for prognostic prediction as well as predicting prognosis more accurately than the current state-of-the-art methods. PASNet is the first pathway-based deep neural network that represents hierarchical representations of genes and pathways and their nonlinear effects, to the best of our knowledge. Additionally, PASNet would be promising due to its flexible model representation and interpretability, embodying the strengths of deep learning. The open-source code of PASNet is available at https://github.com/DataX-JieHao/PASNet.

     

Preferential Utilization of d-Lactate by Shewanella oneidensis

Shewanella oneidensis couples oxidation of lactate to respiration of many substrates. Here we report that llpR (l-lactate-positive regulator, SO_3460) encodes a positive regulator of l-lactate utilization distinct from previously studied regulators. We also demonstrate d-lactate inhibition of l-lactate utilization in S. oneidensis, resulting in preferential utilization of the d isomer.     

Quantitation of the common components of deoxyribonucleic acids by mass spectrometry: application to the analysis of DNAs of unusual composition

A mass spectral method for the quantitation of the percentages of deoxyadenosine, deoxyguanosine, deoxycytidine, and thymidine in intact DNAs has been devised. Standard curves for each nucleoside have been constructed which are based upon the observation that a direct correlation exists between the heights (% deflection) of diagnostic peaks from these nucleosides in a mass spectrum and the published percent composition of specific DNAs. Analyses of DNA from Clostridiumperfringens, Micrococcusluteus, Escherichiacoli, Bacillussubtilis, Pseudomonasfluorescens, Drosophilamelanogaster, salmon sperm, and bacteriophage lambda were used to determine standard curves. The validity of the method was demonstrated by comparison of the results from the mass spectral procedure with results from the chemical analyses of the DNAs from calf thymus and wheat germ. Analysis of ØX-174 DNA yielded values consistent with the published values obtained via sequence analysis and indicated that the method is applicable to both single and double-stranded DNAs. Results from T2 DNA, which contains no cytidine, exhibited artificially high values for adenosine, guanosine and thymidine with concomitant alteration in the A/T and G/C molar ratios. Such skewed results are useful in predicting the presence of modified nucleosides. The extreme sensitivity of the method has been exploited in the analysis of subnanogram quantities of restriction endonuclease fragments from DNA.     

Life without tRNAArg–adenosine deaminase TadA: evolutionary consequences of decoding the four CGN codons as arginine in Mycoplasmas and other Mollicutes

In most bacteria, two tRNAs decode the four arginine CGN codons. One tRNA harboring a wobble inosine (tRNA^Arg_ICG) reads the CGU, CGC and CGA codons, whereas a second tRNA harboring a wobble cytidine (tRNA^Arg_CCG) reads the remaining CGG codon. The reduced genomes of Mycoplasmas and other Mollicutes lack the gene encoding tRNA^Arg_CCG. This raises the question of how these organisms decode CGG codons. Examination of 36 Mollicute genomes for genes encoding tRNA^Arg and the TadA enzyme, responsible for wobble inosine formation, suggested an evolutionary scenario where tadA gene mutations first occurred. This allowed the temporary accumulation of non-deaminated tRNA^Arg_ACG, capable of reading all CGN codons. This hypothesis was verified in Mycoplasma capricolum, which contains a small fraction of tRNA^Arg_ACG with a non-deaminated wobble adenosine. Subsets of Mollicutes continued to evolve by losing both the mutated tRNA^Arg_CCG and tadA, and then acquired a new tRNA^Arg_UCG. This permitted further tRNA^Arg_ACG mutations with tRNA^Arg_GCG or its disappearance, leaving a single tRNA^Arg_UCG to decode the four CGN codons. The key point of our model is that the A-to-I deamination activity had to be controlled before the loss of the tadA gene, allowing the stepwise evolution of Mollicutes toward an alternative decoding strategy.