Computational approaches in cancer multidrug resistance research: Identification of potential biomarkers, drug targets and drug-target interactions.

Like physics in the 19th century, biology and molecular biology in particular, has been fertilized and enhanced like few other scientific fields, by the incorporation of mathematical methods. In the last decades, a whole new scientific field, bioinformatics, has developed with an output of over 30,000 papers a year (Pubmed search using the keyword "bioinformatics"). Huge databases of mass throughput data have been established, with ArrayExpress alone containing more than 2.7 million assays (October 2019). Computational methods have become indispensable tools in molecular biology, particularly in one of the most challenging areas of cancer research, multidrug resistance (MDR). However, confronted with a plethora of different algorithms, approaches, and methods, the average researcher faces key questions: Which methods do exist? Which methods can be used to tackle the aims of a given study? Or, more generally, how do I use computational biology/bioinformatics to bolster my research? The current review is aimed at providing guidance to existing methods with relevance to MDR research. In particular, we provide an overview on: a) the identification of potential biomarkers using expression data; b) the prediction of treatment response by machine learning methods; c) the employment of network approaches to identify gene/protein regulatory networks and potential key players; d) the identification of drug-target interactions; e) the use of bipartite networks to identify multidrug targets; f) the identification of cellular subpopulations with the MDR phenotype; and, finally, g) the use of molecular modeling methods to guide and enhance drug discovery. This review shall serve as a guide through some of the basic concepts useful in MDR research. It shall give the reader some ideas about the possibilities in MDR research by using computational tools, and, finally, it shall provide a short overview of relevant literature.

Patterns of genomic evolution in advanced melanoma.

Genomic alterations occurring during melanoma progression and the resulting genomic heterogeneity between metastatic deposits remain incompletely understood. Analyzing 86 metastatic melanoma deposits from 53 patients with whole-exome sequencing (WES), we show a low branch to trunk mutation ratio and little intermetastatic heterogeneity, with driver mutations almost completely shared between lesions. Branch mutations consistent with UV damage indicate that metastases may arise from different subclones in the primary tumor. Selective gain of mutated BRAF alleles occurs as an early event, contrasting whole-genome duplication (WGD) occurring as a late truncal event in about 40% of cases. One patient revealed elevated mutational diversity, probably related to previous chemotherapy and DNA repair defects. In another patient having received radiotherapy toward a lymph node metastasis, we detected a radiotherapy-related mutational signature in two subsequent distant relapses, consistent with secondary metastatic seeding. Our findings add to the understanding of genomic evolution in metastatic melanomas.

High number of kinome-mutations in non-small cell lung cancer is associated with reduced immune response and poor relapse-free survival.

Lung cancer is the leading cause of cancer related death, and the past years' improved insight into underlying molecular events has significantly improved outcome for specific subsets of patients. In particular, several new therapies that target protein kinases have been implemented, and many more are becoming available. We have investigated lung cancer specimens for somatic mutations in a targeted panel of 612 human genes, the majority being protein kinases. The somatic mutation profiles were correlated to profiles of immune cell infiltration as well as relapse-free survival. Targeted deep sequencing was performed on 117 tumour/normal pairs using the SureSelect Human Kinome kit (Agilent Technologies), with capture probes targeting 3.2 Mb of the human genome, including exons and untranslated regions of all known kinases, kinase receptors and selected cancer-related genes (612 genes in total). CD8 staining was determined using Ventana Benchmark. Survival analyses were performed using SPSS. The number of mutations per sample ranged from 0 to 50 (within the 612 genes tested), with a median of nine. The prognosis was worse for patients with more than the median number of mutations. A significant correlation was found between mutations in one of selected DNA-repair genes and the total number of mutations in that tumour (p < 0.001). There was a significant inverse correlation between the number of infiltrating stromal CD8+ lymphocytes and the presence of EGFR mutations.

BRAF V600E mutation in early-stage multiple myeloma: good response to broad acting drugs and no relation to prognosis.

In this study, we analyzed the prevalence and clone size of BRAF V600E mutation in 209 patients with multiple myeloma and related the results to clinical phenotype, response and survival. Biopsies were screened for BRAF V600E by allele-specific real-time PCR (AS-PCR). Positive results were confirmed by immunohistochemistry, Sanger sequencing and, in three patients from whom we had stored purified myeloma cells, whole-exome sequencing. Eleven patients (5.3%) were BRAF V600E mutation positive by AS-PCR and at least one other method. The fraction of mutated cells varied from 4 to 100%. BRAF V600E-positive patients had no characteristic clinical phenotype except for significantly higher levels of serum creatinine (125 versus 86 µmol/l) Seven of eleven patients responded with at least very good partial response to alkylators, immunomodulatory agents or proteasome inhibitors. Progression-free and overall survival were similar in patients with and without the mutation. By this integrated approach, we found that patients with BRAF V600E mutation responded very well to broad acting drugs and there was no relation to prognosis in early-stage myeloma. In particular, a large mutated cell fraction did not correlate with aggressive disease.

Increased expression of IRF4 and ETS1 in CD4+ cells from patients with intermittent allergic rhinitis.

BACKGROUND: The transcription factor (TF) IRF4 is involved in the regulation of Th1, Th2, Th9, and Th17 cells, and animal studies have indicated an important role in allergy. However, IRF4 and its target genes have not been examined in human allergy. METHODS: IRF4 and its target genes were examined in allergen-challenged CD4(+) cells from patients with IAR, using combined gene expression microarrays and chromatin immunoprecipitation chips (ChIP-chips), computational target prediction, and RNAi knockdowns. RESULTS: IRF4 increased in allergen-challenged CD4(+) cells from patients with IAR, and functional studies supported its role in Th2 cell activation. IRF4 ChIP-chip showed that IRF4 regulated a large number of genes relevant to Th cell differentiation. However, neither Th1 nor Th2 cytokines were the direct targets of IRF4. To examine whether IRF4 induced Th2 cytokines via one or more downstream TFs, we combined gene expression microarrays, ChIP-chips, and computational target prediction and found a putative intermediary TF, namely ETS1 in allergen-challenged CD4(+) cells from allergic patients. ETS1 increased significantly in allergen-challenged CD4(+) cells from patients compared to controls. Gene expression microarrays before and after ETS1 RNAi knockdown showed that ETS1 induced Th2 cytokines as well as disease-related pathways. CONCLUSIONS: Increased expression of IRF4 in allergen-challenged CD4(+) cells from patients with intermittent allergic rhinitis leads to activation of a complex transcriptional program, including Th2 cytokines.

CLC and IFNAR1 are differentially expressed and a global immunity score is distinct between early- and late-onset colorectal cancer.

Colorectal cancer (CRC) incidence increases with age, and early onset of the disease is an indication of genetic predisposition, estimated to cause up to 30% of all cases. To identify genes associated with early-onset CRC, we investigated gene expression levels within a series of young patients with CRCs who are not known to carry any hereditary syndromes (n=24; mean 43 years at diagnosis), and compared this with a series of CRCs from patients diagnosed at an older age (n=17; mean 79 years). Two individual genes were found to be differentially expressed between the two groups, with statistical significance; CLC was higher and IFNAR1 was less expressed in early-onset CRCs. Furthermore, genes located at chromosome band 19q13 were found to be enriched significantly among the genes with higher expression in the early-onset samples, including CLC. An elevated immune content within the early-onset group was observed from the differentially expressed genes. By application of outlier statistics, H3F3A was identified as a top candidate gene for a subset of the early-onset CRCs. In conclusion, CLC and IFNAR1 were identified to be overall differentially expressed between early- and late-onset CRC, and are important in the development of early-onset CRC.

Connectivity can be used to identify key genes in DNA microarray data: a study based on gene expression in nasal polyps before and after treatment with glucocorticoids.

CONCLUSIONS: The presented analysis of nasal polyposis using connectivity based on the PubGene literature co-citation network demonstrates that this tool can be used to identify key genes in DNA microarray studies of human polygenic diseases. OBJECTIVES: DNA microarray studies of complex diseases may reveal differential expression of hundreds of genes. According to network theory and studies of yeast cells, genes that are connected with several other genes appear to have key regulatory roles. This study aimed to examine if this principle can be translated to DNA microarray studies of human disease, using nasal polyposis as a base for the analysis. MATERIALS AND METHODS: The connectivity of differentially expressed genes from a previously described microarray study of nasal polyposis before and after treatment with glucocorticoids was determined. This was done using the literature co-citation network PubGene. RESULTS: In all, 166 genes were differentially expressed; 39 of these were previously defined as inflammatory and considered important for nasal polyposis. The connectivity of all differentially expressed genes was analysed using the PubGene literature co-citation network. Seventy-four of the 166 genes were connected to other genes. By contrast, the average number of connected genes among 100 sets of 166 randomly chosen genes was 31.5. A small number of the differentially expressed genes were highly connected, while most genes had few or no connections. This indicated a scale-free network. The most connected gene was interleukin-8, an inflammatory gene of known importance for nasal polyposis. Twenty-eight of the 74 connected genes were inflammatory (38%), compared with 11 of the 92 unconnected genes (12%), p < 0.0001. Since most evidence suggests that nasal polyps are inflammatory in their nature, this supports the hypothesis that connected genes have more disease relevance than unconnected genes.

Photochemically induced gene silencing using small interfering RNA molecules in combination with lipid carriers.

Novel strategies for efficient delivery of small interfering RNA (siRNA) molecules with a potential for targeting are required for development of RNA interference (RNAi) therapeutics. Here, we present a strategy that is based on delivery of siRNA molecules through the endocytic pathway, in order to develop a method for site-specific gene silencing. To achieve this, we combined the use of cationic lipids and photochemical internalization (PCI). Using the human S100A4 gene as a model system, we obtained potent gene silencing in four tested human cancer cell lines following PCI induction when using the cationic lipid jetSI-ENDO. Gene silencing was shown at both the RNA and protein levels, with no observed PCI toxicity when using the jetSI reagent and an optimized PCI protocol. This novel induction method opens for in vivo site-specific delivery of siRNA molecules toward a sequence of interest.

Photochemically induced gene silencing using PNA-peptide conjugates.

The potential for exploration of peptide nucleic acid (PNA) as an experimental and therapeutic regulator of gene expression has been hampered by a poor delivery and a lack of site-specific targeting. In the present study, we have developed an efficient strategy for nuclear delivery of PNA by combining cationically charged PNA-peptide conjugates and photochemical internalization (PCI) technology. When using the S100A4 gene as a model system, a consistent downregulation to around 10% remaining protein signal was obtained in three selected cell lines. Furthermore, a dose-dependent and time-dependent inhibition of the S100A4 protein was demonstrated. A main benefit of the strategy proposed is the possibility of site-specific targeting.

BRCA1 mutations in ovarian cancer and borderline tumours in Norway: a nested case-control study.

The aims of the present study were to find the frequency of the most common BRCA1 mutations in women with ovarian tumours identified from a population-based cancer registry and in the general population, to estimate the relative risk of ovarian tumours among the mutation carriers, and to explore the value of using CA125 as a prediagnostic test. The study was designed as a nested case-control study within a cohort mainly consisting of participants in population-based health examinations. The data files of The Cancer Registry of Norway and the Janus serum bank were linked to identify cases with ovarian cancer and borderline tumours. Hereditary BRCA1 mutations were determined using archived serum samples and capillary electrophoresis. Altogether 478 ovarian cancer patients and 190 patients with borderline tumours were identified, and 1421 and 568 matching controls were selected. Odds ratios (OR) of developing ovarian cancer and borderline tumours in the presence of BRCA1 mutations and CA125 level were derived from conditional logistic regression models. Among the 478 ovarian cancer patients, 19 BRCA1 mutations were identified (1675delA, 1135insA, 816delGT and 3347delAG), none among the patients with borderline tumours. Only two of the 1989 controls were BRCA1 mutation carriers (0.10%). The risk of ovarian cancer among the mutation carriers was strongly elevated (OR=29, 95% CI=6.6-120). CA125 was a marker for ovarian cancer, but the sensitivity was low. This study showed that BRCA1 mutation carriers have a very high risk of ovarian cancer. However, since the prevalence of BRCA1 mutations in the Norwegian population was low, the proportion of ovarian cancers due to BRCA1 mutations seemed to be low, about 4%. The sensitivity of using CA125 only as a screening test for ovarian cancer was low.

Interferon-gamma suppresses S100A4 transcription independently of apoptosis or cell cycle arrest.

The S100A4 protein has been associated with increased metastatic capacity of cancer cells, and recent studies have suggested a correlation between the expression level of S100A4 and the prognostic outcome for patients with various types of cancer. The knowledge about the mechanisms underlying the metastasis-promoting effects is still limited, and the aim of the present study was to elucidate signal transduction pathways involved in the regulation of S100A4. After treatment of human carcinoma cells with interferon-gamma (IFN-gamma), we observed downregulation of S100A4 both at mRNA and protein levels. The effect was not dependent on IFN-gamma-induced apoptosis or IFN-gamma-mediated cell cycle arrest. Moreover, IFN-gamma-mediated decrease in mRNA stability could not account for the observed decrease in S100A4 transcript level. Finally, microarray analysis suggests ISGF3G, ETV5, ZNF133 and CEBPG as possible candidate genes involved in IFN-gamma-mediated repression of S100A4.

Associations between gene expressions in breast cancer and patient survival.

We analyzed associations between gene expression in breast cancer and patient survival for 8024 genes from a previously published microarray data set. Analysis of survival, by using the logrank test, was performed automatically for each gene. After correcting for multiple testing, we identified 95 genes whose expression was significantly associated with patient survival. The independent prognostic value of the genes ranking the highest in univariate analysis, together with clinical parameters, was assessed by Cox multivariate regression analysis. The P-values from these logrank tests were also mapped to chromosomal positions and compared with previously reported amplicon regions. We used PubGene web tools to identify groups of genes that had co-occurred in the literature and whose expression patterns were associated with survival. Our analyses demonstrate the comprehensiveness of the microarray technology with respect to measuring gene expression and indicate that the technology may be used to screen for potential clinical markers.

Genetic epidemiology of BRCA1 mutations in Norway.

Familial breast-ovarian cancer has been demonstrated to be frequent but unevenly distributed in Norway. This was assumed to be caused by the reduced population size created by the medieval Bubonic plagues 25 generations ago, and by the following rapid expansion. We have previously reported that four mutations account for 68% of the BRCA1 mutation carriers. Subsequent analysis has resulted in a total of 100 separate families carrying one of these founder mutations. The four mutations occurred on one specific BRCA1 haplotype each. The 1675delA, 816delGT and 3347delAG families originated from the South-West coast of Norway with a few families in the north, while the traceable ancestors of the 1135insA families clustered along the historical inland road from the South-East to mid-Norway. The carriers of each of the four mutations today are descendants of one or a few individuals surviving the plagues. We may identify the majority of BRCA1 mutation carriers in Norway by screening for local founder mutations.

[New knowledge derived from measurement of gene expression with the DNA microarray method]. / Ny kunnskap fra måling av genuttrykk gjennom DNA-mikromatrise.

BACKGROUND: The cDNA microarray method offers the first possibility of obtaining a global understanding of biological processes in living organisms, by simultaneous read-outs of tens of thousands of mRNAs. Initial experiments suggest that genes with similar function have similar expression patterns. MATERIAL AND METHODS: Understanding this level of biological complexity will, however, require completely new approaches to data analysis. Computer science methods, such as data mining and knowledge discovery, can synthesize interpretable if-then rules that model the relation between gene expressions and functions and use the rules to classify unknown genes. The huge body of existing biological and medical knowledge makes it necessary to develop methods for extracting knowledge from such repositories. RESULTS: Models of relations between gene expressions and gene functions in a data set from a publicly available source are synthesized semiautomatically and applied to classify unknown genes. Encouraging results have been achieved. The method is applied in the analysis of data from our microarray system which has recently become operational. INTERPRETATION: The principles are of general importance and will be used to evaluate a wide range of complex data sets like decision support in clinical medicine, for situations in which physicians need to handle a large volume of data for each patient.

Optimization of hammerhead ribozymes for the cleavage of S100A4 (CAPL) mRNA.

Previously, suppression of the S100A4 mRNA by an endogenously expressed ribozyme in osteosarcoma cells was shown to inhibit their metastasis in rats. As a prelude to performing similar studies with exogenous, synthetic ribozymes, we compared a series of hammerhead ribozymes targeted against different sites in the mRNA. The ribozymes differed only in the 7-base flanking sequences complementary to the substrate and were protected against nucleases by chemical modification. Cleavage efficiency varied widely and was not obviously related to the predicted secondary structure of the target RNA. The most active ribozyme of the series was chosen for further optimization. Lengthening its flanking sequences was counterproductive and reduced cleavage even when using excess ribozyme. Using excess substrate (multiple-turnover kinetics), cleavage was fastest with the (6+8) ribozyme having 6 nucleotides (nt) in stem III and 8 nt in stem I. Although these stems strongly influence ribozyme performance, their optimization is still empirical. Faster cleavage was obtained by adding facilitator oligonucleotides to ribozymes with shorter stems of (6+6) and (5+5) nt. Stimulation was particularly strong in the case of the (5+5) ribozyme, which was poorly active by itself. The enhancement caused by different facilitator oligonucleotides paralleled their expected ability to hybridize to RNA as a function of length and chemical modification.

The BRCA1 syndrome and other inherited breast or breast-ovarian cancers in a Norwegian prospective series.

Inherited breast cancer is a heterogenous group of diseases. We examined this heterogeneity in a prospective series of inherited breast and ovarian cancers, previously demonstrated to include 84% of inherited cancers. Ninety-two tumours (65 breast and 27 ovarian) in 82 patients from 70 kindreds were prospectively diagnosed. Fifteen of the breast cancers were in situ, 50 were infiltrating. 40 (49%) of the 82 women carried a BRCA1 mutation, whereas no mutation in BRCA2 was found. Approximately, two-thirds of the BRCA1 mutation carriers had one of the four most frequent Norwegian founder mutations. Ninety-five per cent of the epithelial ovarian cancers occurred in BRCA1 mutation carrying women versus 38% of infiltrating breast cancers and 7% of carcinoma in situ of the breast. The BRCA1 syndrome was phenotypically distinct with invasive, high grade, oestrogen receptor-negative breast cancers and epithelial ovarian cancers. Non-BRCA1/2 inherited breast cancers included carcinoma in situ and lobular carcinoma and were frequently bilateral. Non-BRCA1/2 inherited breast cancer is not associated with epithelial ovarian cancer and in breast cancers has distinct biological characteristics, indicating that the different subgroups of inherited breast cancer may need different healthcare services.

A literature network of human genes for high-throughput analysis of gene expression.

We have carried out automated extraction of explicit and implicit biomedical knowledge from publicly available gene and text databases to create a gene-to-gene co-citation network for 13,712 named human genes by automated analysis of titles and abstracts in over 10 million MEDLINE records. The associations between genes have been annotated by linking genes to terms from the medical subject heading (MeSH) index and terms from the gene ontology (GO) database. The extracted database and accompanying web tools for gene-expression analysis have collectively been named 'PubGene'. We validated the extracted networks by three large-scale experiments showing that co-occurrence reflects biologically meaningful relationships, thus providing an approach to extract and structure known biology. We validated the applicability of the tools by analyzing two publicly available microarray data sets.

[Gene therapy in cancer]. / Genterapi ved kreft.

BACKGROUND: There has been a substantial increase in our understanding of the pathogenesis of cancer over the past decades, and new treatment modalities are being developed on the basis of this knowledge. MATERIAL AND METHODS: The literature is reviewed, and the status of gene therapy protocols approved in Norway is presented. RESULTS: About 70% of the more than 400 clinical gene therapy studies started are targeted at cancer. Several different principles are used, including gene therapy targeted at tumour suppressor genes, oncogenes, or central signaling molecules, as well as "suicide gene" therapy. In addition, various gene therapy protocols aim at strengthening immune responses. Most studies have been early clinical trials primarily designed to study safety, applicability and toxicity. Several of these phase I and II studies have, however, shown partial remission of tumours and, in rare cases, complete remission, although curation has not yet been shown. In some trials, including TP53 gene therapy trials, effects on tumour size have been observed in up to 50% of the patients. Up until now, only two phase III and one phase II/III studies have been initiated, but results from these studies have not yet been published. The two first gene therapy protocols approved in Norway are also targeted at cancer. So far, six patients in Norway have undergone gene therapy. INTERPRETATION: As of today, the results should be seen as promising for some of the principles which are being tried out; their clinical importance must, however, be documented in larger controlled clinical trials.

Detection of mutations by single-strand conformation polymorphism (SSCP) analysis and SSCP-hybrid methods.

Single-strand conformation polymorphism (SSCP) analysis detects mutations based on the fact that single-nucleotide changes in DNA sequences alter the mobility of single-stranded DNA in nondenaturing gels. Four methods for detecting mutations based on SSCP are described here. (1) Traditional SSCP analysis is technically easy and can be used for screening large numbers of samples. SSCP-hybrid methods detect mutations based on either an SSCP effect or an altered component independent of the SSCP effect. (2) Dideoxy fingerprinting (ddF) involves PCR amplification of the target and creation of a set of dideoxy-terminated strands with the mutation. (3) Bi-directional dideoxy fingerprinting (Bi-ddF) involves production of two sets of dideoxy-terminated strands that are generated from two different primers. (4) Restriction endonuclease fingerprinting (REF) involves cleavage of the amplified target with five to six groups of restriction endonucleases.

Detection of mutations by denaturing gradient gel electrophoresis.

This unit describes the procedure for determining the melting profile for a given PCR-amplified sequence by perpendicular denaturing gradient gel electrophoresis (DGGE) and, using that information, for developing a screening assay based on either parallel DGGE, CDGE (Constant Denaturant Gel Electrophoresis), or TTGE (Temporal/Temperature Gradient Electrophoresis). Four support protocols describe techniques for pouring perpendicular and parallel denaturing gradient gels, constant denaturant gels, and temporal temperature gradient gels.