EXCAVATOR: a computer program for efficiently mining gene expression data

Genetic Determinants of EGFR-Driven Lung Cancer Growth and Therapeutic Response In Vivo

In lung adenocarcinoma, oncogenic EGFR mutations co-occur with many tumor suppressor gene alterations; however, the extent to which these contribute to tumor growth and response to therapy in vivo remains largely unknown. By quantifying the effects of inactivating 10 putative tumor suppressor genes in a mouse model of EGFR-driven Trp53-deficient lung adenocarcinoma, we found that Apc, Rb1, or Rbm10 inactivation strongly promoted tumor growth. Unexpectedly, inactivation of Lkb1 or Setd2-the strongest drivers of growth in a KRAS-driven model-reduced EGFR-driven tumor growth. These results are consistent with mutational frequencies in human EGFR- and KRAS-driven lung adenocarcinomas. Furthermore, KEAP1 inactivation reduced the sensitivity of EGFR-driven tumors to the EGFR inhibitor osimertinib, and mutations in genes in the KEAP1 pathway were associated with decreased time on tyrosine kinase inhibitor treatment in patients. Our study highlights how the impact of genetic alterations differs across oncogenic contexts and that the fitness landscape shifts upon treatment. SIGNIFICANCE: By modeling complex genotypes in vivo, this study reveals key tumor suppressors that constrain the growth of EGFR-mutant tumors. Furthermore, we uncovered that KEAP1 inactivation reduces the sensitivity of these tumors to tyrosine kinase inhibitors. Thus, our approach identifies genotypes of biological and therapeutic importance in this disease.This article is highlighted in the In This Issue feature, p. 1601.

Mouse-human co-clinical trials demonstrate superior anti-tumour effects of buparlisib (BKM120) and cetuximab combination in squamous cell carcinoma of head and neck

Background

Recurrent and/or metastatic squamous cell carcinoma of head and neck (R/M SCCHN) is a common cancer with high recurrence and mortality. Current treatments have low response rates (RRs).

Methods

Fifty-three patients with R/M SCCHN received continuous oral buparlisib. In parallel, patient-derived xenografts (PDXs) were established in mice to evaluate resistance mechanisms and efficacy of buparlisib/cetuximab combination. Baseline and on-treatment tumour genomes and transcriptomes were sequenced. Based on the integrated clinical and PDX data, 11 patients with progression under buparlisib monotherapy were treated with a combination of buparlisib and cetuximab.

Results

For buparlisib monotherapy, disease control rate (DCR) was 49%, RR was 3% and median progression-free survival (PFS) and overall survival (OS) were 63 and 143 days, respectively. For combination therapy, DCR was 91%, RR was 18% and median PFS and OS were 111 and 206 days, respectively. Four PDX models were originated from patients enrolled in the current clinical trial. While buparlisib alone did not inhibit tumour growth, combination therapy achieved tumour inhibition in three of seven PDXs. Genes associated with apoptosis and cell-cycle arrest were expressed at higher levels with combination treatment than with buparlisib or cetuximab alone.

Conclusions

The buparlisib/cetuximab combination has significant promise as a treatment strategy for R/M SCCHN.

Clinical Trial Registration

NCT01527877.

Whole Exome and Transcriptome Analyses Integrated with Microenvironmental Immune Signatures of Lung Squamous Cell Carcinoma

The immune microenvironment in lung squamous cell carcinoma (LUSC) is not well understood, with interactions between the host immune system and the tumor, as well as the molecular pathogenesis of LUSC, awaiting better characterization. To date, no molecularly targeted agents have been developed for LUSC treatment. Identification of predictive and prognostic biomarkers for LUSC could help optimize therapy decisions. We sequenced whole exomes and RNA from 101 tumors and matched noncancer control Korean samples. We used the information to predict subtype-specific interactions within the LUSC microenvironment and to connect genomic alterations with immune signatures. Hierarchical clustering based on gene expression and mutational profiling revealed subtypes that were either immune defective or immune competent. We analyzed infiltrating stromal and immune cells to further characterize the tumor microenvironment. Elevated expression of macrophage 2 signature genes in the immune competent subtype confirmed that tumor-associated macrophages (TAM) linked inflammation and mutation-driven cancer. A negative correlation was evident between the immune score and the amount of somatic copy-number variation (SCNV) of immune genes (r = -0.58). The SCNVs showed a potential detrimental effect on immunity in the immune-deficient subtype. Knowledge of the genomic alterations in the tumor microenvironment could be used to guide design of immunotherapy options that are appropriate for patients with certain cancer subtypes. Cancer Immunol Res; 6(7); 848-59. ©2018 AACR.

A Hybrid Approach for Biomarker Discovery from Microarray Gene Expression Data for Cancer Classification

Microarrays allow researchers to monitor the gene expression patterns for tens of thousands of genes across a wide range of cellular responses, phenotype and conditions. Selecting a small subset of discriminate genes from thousands of genes is important for accurate classification of diseases and phenotypes. Many methods have been proposed to find subsets of genes with maximum relevance and minimum redundancy, which can distinguish accurately between samples with different labels. To find the minimum subset of relevant genes is often referred as biomarker discovery. Two main approaches, filter and wrapper techniques, have been applied to biomarker discovery. In this paper, we conducted a comparative study of different biomarker discovery methods, including six filter methods and three wrapper methods. We then proposed a hybrid approach, FR-Wrapper, for biomarker discovery. The aim of this approach is to find an optimum balance between the precision of the biomarker discovery and the computation cost, by taking advantages of both filter method's efficiency and wrapper method's high accuracy. Our hybrid approach applies Fisher's ratio, a simple method easy to understand and implement, to filter out most of the irrelevant genes, then a wrapper method is employed to reduce the redundancy. The performance of FR-Wrapper approach is evaluated over four widely used microarray datasets. Analysis of experimental results reveals that the hybrid approach can achieve the goal of maximum relevance with minimum redundancy.

Express path analysis identifies a tyrosine kinase Src-centric network regulating divergent host responses to Mycobacterium tuberculosis infection

Global gene expression profiling has emerged as a major tool in understanding complex response patterns of biological systems to perturbations. However, a lack of unbiased analytical approaches has restricted the utility of complex microarray data to gain novel system level insights. Here we report a strategy, express path analysis (EPA), that helps to establish various pathways differentially recruited to achieve specific cellular responses under contrasting environmental conditions in an unbiased manner. The analysis superimposes differentially regulated genes between contrasting environments onto the network of functional protein associations followed by a series of iterative enrichments and network analysis. To test the utility of the approach, we infected THP1 macrophage cells with a virulent Mycobacterium tuberculosis strain (H37Rv) or the attenuated non-virulent strain H37Ra as contrasting perturbations and generated the temporal global expression profiles. EPA of the results provided details of response-specific and time-dependent host molecular network perturbations. Further analysis identified tyrosine kinase Src as the major regulatory hub discriminating the responses between wild-type and attenuated Mtb infection. We were then able to verify this novel role of Src experimentally and show that Src executes its role through regulating two vital antimicrobial processes of the host cells (i.e. autophagy and acidification of phagolysosome). These results bear significant potential for developing novel anti-tuberculosis therapy. We propose that EPA could prove extremely useful in understanding complex cellular responses for a variety of perturbations, including pathogenic infections.

HAMSTER: visualizing microarray experiments as a set of minimum spanning trees

Background

Visualization tools allow researchers to obtain a global view of the interrelationships between the probes or experiments of a gene expression (e.g. microarray) data set. Some existing methods include hierarchical clustering and k-means. In recent years, others have proposed applying minimum spanning trees (MST) for microarray clustering. Although MST-based clustering is formally equivalent to the dendrograms produced by hierarchical clustering under certain conditions; visually they can be quite different.

Methods

HAMSTER (Helpful Abstraction using Minimum Spanning Trees for Expression Relations) is an open source system for generating a set of MSTs from the experiments of a microarray data set. While previous works have generated a single MST from a data set for data clustering, we recursively merge experiments and repeat this process to obtain a set of MSTs for data visualization. Depending on the parameters chosen, each tree is analogous to a snapshot of one step of the hierarchical clustering process. We scored and ranked these trees using one of three proposed schemes. HAMSTER is implemented in C++ and makes use of Graphviz for laying out each MST.

Results

We report on the running time of HAMSTER and demonstrate using data sets from the NCBI Gene Expression Omnibus (GEO) that the images created by HAMSTER offer insights that differ from the dendrograms of hierarchical clustering. In addition to the C++ program which is available as open source, we also provided a web-based version (HAMSTER⁺) which allows users to apply our system through a web browser without any computer programming knowledge.

Conclusion

Researchers may find it helpful to include HAMSTER in their microarray analysis workflow as it can offer insights that differ from hierarchical clustering. We believe that HAMSTER would be useful for certain types of gradient data sets (e.g time-series data) and data that indicate relationships between cells/tissues. Both the source and the web server variant of HAMSTER are available from http://hamster.cbrc.jp/.

A hybrid approach for biomarker discovery from microarray gene expression data for cancer classification

Microarrays allow researchers to monitor the gene expression patterns for tens of thousands of genes across a wide range of cellular responses, phenotype and conditions. Selecting a small subset of discriminate genes from thousands of genes is important for accurate classification of diseases and phenotypes. Many methods have been proposed to find subsets of genes with maximum relevance and minimum redundancy, which can distinguish accurately between samples with different labels. To find the minimum subset of relevant genes is often referred as biomarker discovery. Two main approaches, filter and wrapper techniques, have been applied to biomarker discovery. In this paper, we conducted a comparative study of different biomarker discovery methods, including six filter methods and three wrapper methods. We then proposed a hybrid approach, FR-Wrapper, for biomarker discovery. The aim of this approach is to find an optimum balance between the precision of the biomarker discovery and the computation cost, by taking advantages of both filter method's efficiency and wrapper method's high accuracy. Our hybrid approach applies Fisher's ratio, a simple method easy to understand and implement, to filter out most of the irrelevant genes, then a wrapper method is employed to reduce the redundancy. The performance of FR-Wrapper approach is evaluated over four widely used microarray datasets. Analysis of experimental results reveals that the hybrid approach can achieve the goal of maximum relevance with minimum redundancy.

Microarray-MD: a system for exploratory analysis of microarray gene expression data

In this paper, we present Microarray Medical Data explorer (Microarray-MD), a novel software system that is able to assist in the exploratory analysis of gene expression microarray data. It implements a combination scheme of multiple Support Vector Machines, which integrates a variety of gene selection criteria and allows for the discrimination of multiple diseases or subtypes of a disease. The system can be trained and automatically tune its parameters with the provision of pathologically characterized gene expression data to its input. Given a set of new, uncharacterized, patient's data as input, it outputs a decision on the type or the subtype of a disease. A graphical user interface provides easy access to the system operations and direct adjustment of its parameters. It has been tested on various publicly available datasets. The overall accuracy it achieves was estimated to exceed 90%.

Text mining biomedical literature for discovering gene-to-gene relationships: a comparative study of algorithms

Partitioning closely related genes into clusters has become an important element of practically all statistical analyses of microarray data. A number of computer algorithms have been developed for this task. Although these algorithms have demonstrated their usefulness for gene clustering, some basic problems remain. This paper describes our work on extracting functional keywords from MEDLINE for a set of genes that are isolated for further study from microarray experiments based on their differential expression patterns. The sharing of functional keywords among genes is used as a basis for clustering in a new approach called BEA-PARTITION in this paper. Functional keywords associated with genes were extracted from MEDLINE abstracts. We modified the Bond Energy Algorithm (BEA), which is widely accepted in psychology and database design but is virtually unknown in bioinformatics, to cluster genes by functional keyword associations. The results showed that BEA-PARTITION and hierarchical clustering algorithm outperformed k-means clustering and self-organizing map by correctly assigning 25 of 26 genes in a test set of four known gene groups. To evaluate the effectiveness of BEA-PARTITION for clustering genes identified by microarray profiles, 44 yeast genes that are differentially expressed during the cell cycle and have been widely studied in the literature were used as a second test set. Using established measures of cluster quality, the results produced by BEA-PARTITION had higher purity, lower entropy, and higher mutual information than those produced by k-means and self-organizing map. Whereas BEA-PARTITION and the hierarchical clustering produced similar quality of clusters, BEA-PARTITION provides clear cluster boundaries compared to the hierarchical clustering. BEA-PARTITION is simple to implement and provides a powerful approach to clustering genes or to any clustering problem where starting matrices are available from experimental observations.

Protein synthesis profiling in the developing brain: a graph theoretic clustering approach

Mapping regional brain development in terms of protein synthesis (PS) activity yields insight on specific spatio-temporal ontogenetic patterns. The biosynthetic activity of an individual brain nucleus is represented as a time-series object, and clustering of time-series contributes to the problem of inducing indicative patterns of brain developmental events and forming respective PS chronological maps. Clustering analysis of PS chronological maps, in comparison with epigenetic influences of alpha2 adrenoceptors treatment, reveals relationships between distantly located brain structures. Clustering is performed with a novel graph theoretic clustering approach (GTC). The approach is based on the weighted graph arrangement of the input objects and the iterative partitioning of the corresponding minimum spanning tree. The final result is a hierarchical clustering-tree organization of the input objects. Application of GTC on the PS patterns in developing brain revealed five main clusters that correspond to respective brain development indicative profiles. The induced profiles confirm experimental findings, and provide evidence for further experimental studies.