Identifying pathways affected by cancer mutations

Mutations in 15 cancers, sourced from the COSMIC Whole Genomes database, and 297 human pathways, arranged into pathway groups based on the processes they orchestrate, and sourced from the KEGG pathway database, have together been used to identify pathways affected by cancer mutations. Genes studied in ≥ 15, and mutated in ≥ 10 samples of a cancer have been considered recurrently mutated, and pathways with recurrently mutated genes have been considered affected in the cancer. Novel doughnut plots have been presented which enable visualization of the extent to which pathways and genes, in each pathway group, are targeted, in each cancer. The ‘organismal systems’ pathway group (including organism-level pathways; e.g., nervous system) is the most targeted, more than even the well-recognized signal transduction, cell-cycle and apoptosis, and DNA repair pathway groups. The important, yet poorly-recognized, role played by the group merits attention. Pathways affected in ≥ 7 cancers yielded insights into processes affected.     

HBV感染者肝细胞样本的关键基因和信号通路的生物信息学挖掘

为探究乙型肝炎病毒(Hepatitis B virus,HBV)感染后肝细胞基因表达和信号通路的改变情况,从Gene Expression Omnibus(GEO)数据库中下载了HBV感染者肝细胞样本制作的基因表达谱数据集GSE83148,进行质量检测、数据标准化后筛选出差异表达基因,进一步做GO和KEGG富集分析以及基因网络相互作用分析,筛选关键基因和信号通路。从HBV感染样本中筛选出fold change≥2,p-value0. 05的上调差异表达基因44个,GO分析获得关键基因BAK1和TP63,差异表达基因网络互作获得5个位于枢纽位置的基因:NDUFS1、NDUFS2、COX7B、ATP5B、OPA1。KEGG分析获得关键信号通路有:乙型肝炎信号通路、病毒癌变信号通路、Fox O信号通路、PI3K-Akt信号通路。本研究筛选出的多数基因与线粒体和氧化呼吸链有关,造成这一现象的具体机制还需进一步探究。     

Identification of key genes and pathways in benign prostatic hyperplasia

Benign prostatic hyperplasia (BPH) is one of the most common causes of lower urinary tract symptoms (LUTS) in elderly man. However, the underlying molecular mechanisms of BPH have not been completely elucidated. We identified the key genes and pathways by using analysis of Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) were identified using edgeR. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed for the DEGs by Database for Annotation, Visualization and Integrated Discovery (DAVID) database and ConsensusPathDB, respectively. Then, protein–protein interaction (PPI) networks were established by the Search Tool for the Retrieval of Interacting Genes (STRING) database and visualized by Cytoscape software. Finally, we identified 660 DEGs ultimately including 268 upregulated genes and 392 downregulated genes. GO analysis revealed that DEGs were mainly enriched in extracellular exosome, identical protein binding, mitochondrial adenosine triphosphate (ATP) synthesis coupled proton transport, extracelluar matrix, focal adhesion, cytosol, Golgi apparatus, cytoplasm, protein binding, and Golgi membrane. Focal adhesion pathway, FoxO signaling pathway, and autophagy pathway were selected. Ubiquitin-conjugating enzyme E2 C (UBE2C), serine/threonine kinase (AKT1), mitogen-activated protein kinase 1 (MAPK1), cyclin B1 (CCNB1), polo-like kinase 1 (PLK1) were filtrated as the hub genes according to the degree of connectivity from the PPI network. The five hub genes including UBE2C, AKT1, MAPK1, CCNB1, PLK1 may play key roles in the pathogenesis of benign prostatic hyperplasia (BPH). Focal adhesion pathway, FoxO signaling pathway, and autophagy pathway may be crucial for the progression of BPH.     

WWOX,large common fragile site genes,and cancer

WWOX is a gene that spans an extremely large chromosomal region. It is derived from within chromosomal band 16q23.2 which is a region with frequent deletions and other alterations in a variety of different cancers. This chromosomal band also contains the FRA16D common fragile site (CFS). CFSs are chromosomal regions found in all individuals which are highly unstable. WWOX has also been demonstrated to function as a tumor suppressor that is involved in the development of many cancers. Two other highly unstable CFSs, FRA3B (3p14.2) and FRA6E (6q26), also span extremely large genes, FHIT and PARK2, respectively, and these two genes are also found to be important tumor suppressors. There are a number of interesting similarities between these three large CFS genes. In spite of the fact that they are derived from some of the most unstable chromosomal regions in the genome, they are found to be highly evolutionarily conserved and the chromosomal region spanning the mouse homologs of both WWOX and FHIT are also CFSs in mice. Many of the other CFSs also span extremely large genes and many of these are very attractive tumor suppressor candidates. WWOX is therefore a member of a very interesting family of very large CFS genes.     

The role of mitogen- and stress-activated protein kinase pathways in melanoma

Recent discoveries have increased our comprehension of the molecular signaling events critical for melanoma development and progression. Many oncogenes driving melanoma have been identified, and most of them exert their oncogenic effects through the activation of the RAF/MEK/ERK mitogen-activated protein kinase (MAPK) pathway. The c-Jun N-terminal kinase (JNK) and p38 MAPK pathways are also important in melanoma, but their precise role is not clear yet. This review summarizes our current knowledge on the role of the three main MAPK pathways, extracellular regulated kinase (ERK), JNK, and p38, and their impact on melanoma biology. Although the results obtained with BRAF inhibitors in melanoma patients are impressive, several mechanisms of acquired resistance have emerged. To overcome this obstacle constitutes the new challenge in melanoma therapy. Given the major role that MAPKs play in melanoma, understanding their functions and the interconnection among them and with other signaling pathways represents a step forward toward this goal.     

Finding off‐targets,biological pathways,and target diseases for chymase inhibitors via structure‐based systems biology approach

Off‐target binding connotes the binding of a small molecule of therapeutic significance to a protein target in addition to the primary target for which it was proposed. Progressively such off‐targeting is emerging to be regular practice to reveal side effects. Chymase is an enzyme of hydrolase class that catalyzes hydrolysis of peptide bonds. A link between heart failure and chymase is ascribed, and a chymase inhibitor is in clinical phase II for treatment of heart failure. However, the underlying mechanisms of the off‐target effects of human chymase inhibitors are still unclear. Here, we develop a robust computational strategy that is applicable to any enzyme system and that allows the prediction of drug effects on biological processes. Putative off‐targets for chymase inhibitors were identified through various structural and functional similarity analyses along with molecular docking studies. Finally, literature survey was performed to incorporate these off‐targets into biological pathways and to establish links between pathways and particular adverse effects. Off‐targets of chymase inhibitors are linked to various biological pathways such as classical and lectin pathways of complement system, intrinsic and extrinsic pathways of coagulation cascade, and fibrinolytic system. Tissue kallikreins, granzyme M, neutrophil elastase, and mesotrypsin are also identified as off‐targets. These off‐targets and their associated pathways are elucidated for the effects of inflammation, cancer, hemorrhage, thrombosis, and central nervous system diseases (Alzheimer's disease). Prospectively, our approach is helpful not only to better understand the mechanisms of chymase inhibitors but also for drug repurposing exercises to find novel uses for these inhibitors. Proteins 2015; 83:1209–1224. © 2014 Wiley Periodicals, Inc.     

Transcriptome analysis of differentially expressed genes and pathways associated with mitoxantrone treatment prostate cancer

The global physiological function of specifically expressed genes of mitoxantrone (MTX)‐resistant prostate cancer (PCa) is unclear. In this study, gene expression pattern from microarray data was investigated for identifying differentially expressed genes (DEGs) in MTX‐resistant PCa xenografts. Human PCa cell lines DU145 and PC3 were cultured in vitro and xenografted into severe combined immunodeficiency (SCID) mice, treated with MTX intragastrically, three times a week until all mice relapsed. Gene expression profiles of the xenografts from castrated mice were performed with Affymetrix human whole genomic oligonucleotide microarray. The Cytoscape software was used to investigate the relationship between proteins and the signalling transduction network. A total of 355 overlapping genes were differentially expressed in MTX‐resistant DU145R and PC3R xenografts. Of these, 16 genes were selected to be validated by quantitative real‐time PCR (qRT‐PCR) in these xenografts, and further tested in a set of formalin‐fixed, paraffin‐embedded and optimal cutting temperature (OCT) clinical tumour samples. Functional and pathway enrichment analyses revealed that these DEGs were closely related to cellular activity, androgen synthesis, DNA damage and repair, also involved in the ERK/MAPK, PI3K/serine‐threonine protein kinase, also known as protein kinase B, PKB (AKT) and apoptosis signalling pathways. This exploratory analysis provides information about potential candidate genes and may bring new insights into the molecular cascade involvement in MTX‐resistant PCa.     

生物信息学分析筛选结直肠癌靶基因及评估预后价值

为寻找与结直肠癌发展和预后相关的潜在关键基因及信号通路.从美国国立信息中心NCBI的GEO数据库获得结直肠癌基因表达数据集GSE106582,通过PCA对样本进行分组,利用GEO2R进行综合分析,筛选结直肠癌与癌旁对照组的差异表达基因;通过DAVID在线工具对差异表达基因进行GO本体分析和KEGG通路富集分析,初步分析...     

Genes acting in longevity-related pathways in the endoparasitoid,Pteromalus puparum

Among insects, lifespans vary over a broad range, from the short-lived mayflies to the 17-year periodical cicadas. Generally, lifespans are determined by a phase in life, the reproductive lifespan, which varies among species. Numerous pathways, such as the insulin/insulin-like growth factor signaling pathway, the target of rapamycin pathway and the mitogen-activated protein kinase/extracellular signal-regulated kinases pathways, influence aging and lifespan. Components of these pathways were identified as lifespan-related genes, including genes mediating growth, metabolism, development, resistance, and other processes. Many age-related genes have been discovered in fruit flies, honeybees, and ants among other insect species. Studies of insect aging and longevity can help understand insect biology and develop new pest management technologies. In this paper, we interrogated the new Pteromalus puparum genome, from which we predicted 133 putative lifespan-related genes based on their homology with known lifespan-related genes of Drosophila melanogaster. These genes function in five signaling pathways and three physiological processes. The conserved domain structures of these genes were predicted and their expression patterns were analyzed. Amino acid sequence alignments and domain structure analysis indicate that most components remain conserved across at least six insect orders. The data in this paper will facilitate future work on parasitoid lifespans, which may have economic value in biocontrol programs.     

Hub genes identification and association of key pathways with hypoxia in cancer cells: A bioinformatics analysis

Three human cancer cell lines (A549, HCT116, and HeLa) were used to investigate the molecular mechanisms and potential prognostic biomarkers associated with hypoxia. We obtained gene expression data from Gene Expression Omnibus (GEO) datasets GSE11704, GSE147384, and GSE38061, which included 5 hypoxic and 8 control samples. Using the GEO2R tool and Venn diagram software, we identified common differentially expressed genes (cDEGs). The cDEGs were then subjected to Gene ontology (GO) and Kyoto Encyclopedia of Gene and Genome (KEGG) pathway analysis by employing DAVID. The hub genes were identified from critical PPI subnetworks through CytoHuba plugin and these genes' prognostic significance and expression were verified using Kaplan-Meier analysis and Gene Expression Profiling Interactive Analysis (GEPIA), respectively. The research showed 676 common DEGs (cDEGs), with 207 upregulated and 469 downregulated genes. The STRING analysis showed 673 nodes and 1446 edges in the PPI network. We identified 4 significant modules and 19 downregulated hub genes. GO analysis revealed all of them were majorly involved in ribosomal large subunit assembly and biogenesis, rRNA processing, ribosome biogenesis, translation, RNA & protein binding frequently at the sites of nucleolus and nucleoplasm while 11 were significantly associated with a better prognosis of hypoxic tumors. Our research sheds light on the molecular mechanisms that underpin hypoxia in human cancer cell lines and identifies potential prognostic biomarkers for hypoxic tumors.     

Engineered single‐ and multi‐cell chemotaxis pathways in E. coli

We have engineered the chemotaxis system of Escherichia coli to respond to molecules that are not attractants for wild‐type cells. The system depends on an artificially introduced enzymatic activity that converts the target molecule into a ligand for an E. coli chemoreceptor, thereby enabling the cells to respond to the new attractant. Two systems were designed, and both showed robust chemotactic responses in semisolid and liquid media. The first incorporates an asparaginase enzyme and the native E. coli aspartate receptor to produce a response to asparagine; the second uses penicillin acylase and an engineered chemoreceptor for phenylacetic acid to produce a response to phenylacetyl glycine. In addition, by taking advantage of a ‘hitchhiker’ effect in which cells producing the ligand can induce chemotaxis of neighboring cells lacking enzymatic activity, we were able to design a more complex system that functions as a simple microbial consortium. The result effectively introduces a logical ‘AND’ into the system so that the population only swims towards the combined gradients of two attractants.     

Development,characterization, and application of a 2‐Compartment system to investigate the impact of pH inhomogeneities in large‐scale CHO‐based processes

Large‐scale bioreactors for the production of monoclonal antibodies reach volumes of up to 25 000 L. With increasing bioreactor size, mixing is however affected negatively, resulting in the formation of gradients throughout the reactor. These gradients can adversely affect process performance at large scale. Since mammalian cells are sensitive to changes in pH, this study investigated the effects of pH gradients on process performance. A 2‐Compartment System was established for this purpose to expose only a fraction of the cell population to pH excursions and thereby mimicking a large‐scale bioreactor. Cells were exposed to repeated pH amplitudes of 0.4 units (pH 7.3), which resulted in decreased viable cell counts, as well as the inhibition of the lactate metabolic shift. These effects were furthermore accompanied by increased absolute lactate levels. Continuous assessment of molecular attributes of the expressed target protein revealed that subunit assembly or N‐glycosylation patterns were only slightly influenced by the pH excursions. The exposure of more cells to the same pH amplitudes further impaired process performance, indicating this is an important factor, which influences the impact of pH inhomogeneity. This knowledge can aid in the design of pH control strategies to minimize the effects of pH inhomogeneity in large‐scale bioreactors.     

Identification of key microRNAs and hub genes in non-small-cell lung cancer using integrative bioinformatics and functional analyses

Non-small-cell lung cancer (NSCLC) is an extremely debilitating respiratory malignancy. However, the pathogenesis of NSCLC has not been fully clarified. The main objective of our study was to identify potential microRNAs (miRNAs) and their regulatory mechanism in NSCLC. Using a systematic review, two NSCLC-associated miRNA data sets (GSE102286 and GSE56036) were obtained from Gene Expression Omnibus, and the differentially expressed miRNAs (DE-miRNAs) were accessed by GEO2R. Survival analysis of candidate DE-miRNAs was conducted using the Kaplan-Meier plotter database. To further illustrate the roles of DE-miRNAs in NSCLC, their potential target genes were predicted by miRNet and were annotated by the Database for Annotation, Visualization and Integrated Discovery (DAVID) program. Moreover, the protein-protein interaction (PPI) and miRNA-hub gene regulatory network were established using the STRING database and Cytoscape software. The function of DE-miRNAs in NSCLC cells was evaluated by transwell assay. Compared with normal tissues, a total of eight DE-miRNAs was commonly changed in two data sets. The survival analysis showed that six miRNAs (miR-21-5p, miR-31-5p, miR-708-5p, miR-30a-5p, miR-451a, and miR-126-3p) were significantly correlated with overall survival. Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis indicated that target genes of upregulated miRNAs were enriched in pathways in cancer, microRNAs in cancer and proteoglycans in cancer, while the target genes of downregulated miRNAs were mainly associated with pathways in cancer, the PI3K-Akt signaling pathway and HTLV-I infection. Based on the miRNA-hub gene network and expression analysis, PTEN, EGFR, STAT3, RHOA, VEGFA, TP53, CTNNB1, and KRAS were identified as potential target genes. Furthermore, all six miRNAs exhibited significant effects on NSCLC cell invasion. These findings indicate that six DE-miRNAs and their target genes may play important roles in the pathogenesis of NSCLC, which will provide novel information for NSCLC treatments.     

Metabolic pathways reconstruction by frequency and amplitude response to forced glycolytic oscillations in yeast

The hypothesis that frequency and amplitude response can be used in a complicated metabolic pathway kinetics model for optimal parameter estimation, as speculated by its successful prior usage for a mechanical oscillator and a heterogeneous chemical system, is tested here. Given the complexity of the glycolysis model of yeast chosen, this question is limited to three kinetics parameters of the 87 in the in vitro model developed in the literature. The direct application of the approach, used with the uninformed selection of operating conditions for the oscillation of external glucose concentration, led to miring the data assimilation process in local minima. Application of linear systems theory, however, identified two natural resonant frequencies that, when excited by external forced oscillations of the same frequency, result in the expression of many harmonics in the Fourier spectra, that is, information-rich experiments. A single such information-rich experiment at one of the resonant frequencies was sufficient to break away from the local minima to find the optimum kinetics parameter estimates. The resonant frequencies themselves represent oscillation modes in glycolysis akin to those previously observed. Furthermore, operation of the bioreactor with large amplitude oscillations of glucose feed (25%) leads to enhanced ethanol average yield by 1.6% at the resonant frequency.     

The role of lipid‐related genes,aging‐related processes,and environment in healthspan

The inherent complexity of aging‐related traits can temper progress in unraveling the genetic origins of healthspan. We focus on two generations in the Framingham Heart Study, the original (FHS) and offspring (FHSO) cohorts, to determine whether aging‐related processes in changing environments can substantially impact the role of lipid‐related genes discovered in candidate gene (the apolipoprotein E (APOE) e2/3/4 polymorphism) and genome‐wide (the APOB rs1042034 (C/T)) studies, in regulation of total cholesterol (TC) and onset of cardiovascular disease (CVD). We demonstrate that the APOE e4 allele and APOB CC genotype can play detrimental, neutral, and protective sex‐specific roles in the etiology of CVD at different ages and in different environments. We document antagonistic roles for the e4 allele in the onset of CVD characterized by detrimental effects at younger ages (RR_{≤ 75 years} = 1.49, P = 7.5 × 10^?4) and protective effects at older ages (RR_76+years = 0.77, P = 0.044) for FHS participants. We found that disregarding the role of aging erroneously nullifies the significant effects of the e4 allele in this sample (RR = 0.92, P = 0.387). The leading biogenetic pathways mediating genetic effects on CVD may be more relevant to lipid metabolism for APOB than APOE. Aging‐related processes can modulate the strength of genetic associations with TC in the same individuals at different chronological ages. We found substantial differences in the effects of the same APOE and APOB alleles on CVD and TC across generations. The results suggest that aging‐related processes in changing environments may play key roles in the genetics of healthspan. Detailed systemic integrative analyses may substantially advance the progress.