Pathway analysis of glioblastoma tissue after preoperative treatment with the EGFR tyrosine kinase inhibitor gefitinib--a phase II trial.

Amplification of the epidermal growth factor receptor (EGFR) gene is one of the most common oncogenic alterations in glioblastoma (45%) making it a prime target for therapy. However, small molecule inhibitors of the EGFR tyrosine kinase showed disappointing efficacy in clinical trials for glioblastoma. Here we aimed at investigating the molecular effects of the tyrosine kinase inhibitor gefitinib on the EGFR signaling pathway in human glioblastoma. Twenty-two patients selected for reoperation of recurrent glioblastoma were treated within a phase II trial for 5 days with 500 mg gefitinib before surgery followed by postoperative gefitinib until recurrence. Resected glioblastoma tissues exhibited high concentrations of gefitinib (median, 4.1 μg/g), 20 times higher than respective plasma. EGFR-pathway activity was evaluated with phosphorylation-specific assays. The EGFR was efficiently dephosphorylated in treated patients as compared to a control cohort of 12 patients. However, no significant effect on 12 pathway constituents was detected. In contrast, in vitro treatment of a glioblastoma cell line, BS-153, with endogenous EGFRwt amplification and EGFRvIII expression resulted not only in dephosphorylation of the EGFR, but also of key regulators in the pathway such as AKT. Treating established xenografts of the same cell line as an in vivo model showed dephosphorylation of the EGFR without affecting downstream signal transductors, similar to the human glioblastoma. Taken together, gefitinib reaches high concentrations in the tumor tissue and efficiently dephosphorylates its target. However, regulation of downstream signal transducers in the EGFR pathway seems to be dominated by regulatory circuits independent of EGFR phosphorylation.

Integrative molecular bioinformatics study of human adrenocortical tumors : microRNA, tissue-specific target prediction, and pathway analysis.

MicroRNAs (miRs) are involved in the pathogenesis of several neoplasms; however, there are no data on their expression patterns and possible roles in adrenocortical tumors. Our objective was to study adrenocortical tumors by an integrative bioinformatics analysis involving miR and transcriptomics profiling, pathway analysis, and a novel, tissue-specific miR target prediction approach. Thirty-six tissue samples including normal adrenocortical tissues, benign adenomas, and adrenocortical carcinomas (ACC) were studied by simultaneous miR and mRNA profiling. A novel data-processing software was used to identify all predicted miR-mRNA interactions retrieved from PicTar, TargetScan, and miRBase. Tissue-specific target prediction was achieved by filtering out mRNAs with undetectable expression and searching for mRNA targets with inverse expression alterations as their regulatory miRs. Target sets and significant microarray data were subjected to Ingenuity Pathway Analysis. Six miRs with significantly different expression were found. miR-184 and miR-503 showed significantly higher, whereas miR-511 and miR-214 showed significantly lower expression in ACCs than in other groups. Expression of miR-210 was significantly lower in cortisol-secreting adenomas than in ACCs. By calculating the difference between dCT(miR-511) and dCT(miR-503) (delta cycle threshold), ACCs could be distinguished from benign adenomas with high sensitivity and specificity. Pathway analysis revealed the possible involvement of G2/M checkpoint damage in ACC pathogenesis. To our knowledge, this is the first report describing miR expression patterns and pathway analysis in sporadic adrenocortical tumors. miR biomarkers may be helpful for the diagnosis of adrenocortical malignancy. This tissue-specific target prediction approach may be used in other tumors too.

Meta-analysis of estrogen response in MCF-7 distinguishes early target genes involved in signaling and cell proliferation from later target genes involved in cell cycle and DNA repair.

ABSTRACT: BACKGROUND: Many studies have been published outlining the global effects of 17 beta-estradiol (E2) on gene expression in human epithelial breast cancer derived MCF-7 cells. These studies show large variation in results, reporting between ~100 and ~1500 genes regulated by E2, with poor overlap. RESULTS: We performed a meta-analysis of these expression studies, using the Rank product method to obtain a more accurate and stable list of the differentially expressed genes, and of pathways regulated by E2. We analyzed 9 time-series data sets, concentrating on response at 3-4 hrs (early) and at 24 hrs (late). We found >1000 statistically significant probe sets after correction for multiple testing at 3-4 hrs, and >2000 significant probe sets at 24 hrs. Differentially expressed genes were examined by pathway analysis. This revealed 15 early response pathways, mostly related to cell signaling and proliferation, and 20 late response pathways, mostly related to breast cancer, cell division, DNA repair and recombination. CONCLUSIONS: Our results show that meta-analysis identified more differentially expressed genes than the individual studies, and that these genes act together in networks. These results provide new insight into E2 regulated mechanisms, especially in the context of breast cancer.

Genome-wide enrichment analysis between endometriosis and obesity-related traits reveals novel susceptibility loci.

Endometriosis is a chronic inflammatory condition in women that results in pelvic pain and subfertility, and has been associated with decreased body mass index (BMI). Genetic variants contributing to the heritable component have started to emerge from genome-wide association studies (GWAS), although the majority remain unknown. Unexpectedly, we observed an intergenic locus on 7p15.2 that was genome-wide significantly associated with both endometriosis and fat distribution (waist-to-hip ratio adjusted for BMI; WHRadjBMI) in an independent meta-GWAS of European ancestry individuals. This led us to investigate the potential overlap in genetic variants underlying the aetiology of endometriosis, WHRadjBMI and BMI using GWAS data. Our analyses demonstrated significant enrichment of common variants between fat distribution and endometriosis (P = 3.7 × 10(-3)), which was stronger when we restricted the investigation to more severe (Stage B) cases (P = 4.5 × 10(-4)). However, no genetic enrichment was observed between endometriosis and BMI (P = 0.79). In addition to 7p15.2, we identify four more variants with statistically significant evidence of involvement in both endometriosis and WHRadjBMI (in/near KIFAP3, CAB39L, WNT4, GRB14); two of these, KIFAP3 and CAB39L, are novel associations for both traits. KIFAP3, WNT4 and 7p15.2 are associated with the WNT signalling pathway; formal pathway analysis confirmed a statistically significant (P = 6.41 × 10(-4)) overrepresentation of shared associations in developmental processes/WNT signalling between the two traits. Our results demonstrate an example of potential biological pleiotropy that was hitherto unknown, and represent an opportunity for functional follow-up of loci and further cross-phenotype comparisons to assess how fat distribution and endometriosis pathogenesis research fields can inform each other.

The circadian clock modulates renal sodium handling.

The circadian clock contributes to the control of BP, but the underlying mechanisms remain unclear. We analyzed circadian rhythms in kidneys of wild-type mice and mice lacking the circadian transcriptional activator clock gene. Mice deficient in clock exhibited dramatic changes in the circadian rhythm of renal sodium excretion. In parallel, these mice lost the normal circadian rhythm of plasma aldosterone levels. Analysis of renal circadian transcriptomes demonstrated changes in multiple mechanisms involved in maintaining sodium balance. Pathway analysis revealed the strongest effect on the enzymatic system involved in the formation of 20-HETE, a powerful regulator of renal sodium excretion, renal vascular tone, and BP. This correlated with a significant decrease in the renal and urinary content of 20-HETE in clock-deficient mice. In summary, this study demonstrates that the circadian clock modulates renal function and identifies the 20-HETE synthesis pathway as one of its principal renal targets. It also suggests that the circadian clock affects BP, at least in part, by exerting dynamic control over renal sodium handling.

Evidence for polygenic adaptation to pathogens in the human genome.

Most approaches aiming at finding genes involved in adaptive events have focused on the detection of outlier loci, which resulted in the discovery of individually "significant" genes with strong effects. However, a collection of small effect mutations could have a large effect on a given biological pathway that includes many genes, and such a polygenic mode of adaptation has not been systematically investigated in humans. We propose here to evidence polygenic selection by detecting signals of adaptation at the pathway or gene set level instead of analyzing single independent genes. Using a gene-set enrichment test to identify genome-wide signals of adaptation among human populations, we find that most pathways globally enriched for signals of positive selection are either directly or indirectly involved in immune response. We also find evidence for long-distance genotypic linkage disequilibrium, suggesting functional epistatic interactions between members of the same pathway. Our results show that past interactions with pathogens have elicited widespread and coordinated genomic responses, and suggest that adaptation to pathogens can be considered as a primary example of polygenic selection.

Molecular role of Cx37 in advanced atherosclerosis: a micro-array study.

Recently, we showed that connexin37 (Cx37) protects against early atherosclerotic lesion development by regulating monocyte adhesion. The expression of this gap junction protein is altered in mouse and human atherosclerotic lesions; it is increased in macrophages newly recruited to the lesions and disappears from the endothelium of advanced plaques. To obtain more insight into the molecular role of Cx37 in advanced atherosclerosis, we used micro-array analysis for gene expression profiling in aortas of ApoE(-/-) and Cx37(-/-)ApoE(-/-) mice before and after 18 weeks of cholesterol-rich diet. Out of >15,000 genes, 106 genes were significantly differentially expressed in young mice before diet (P-value of <0.05, fold change of >0.7 or <-0.7, and intensity value >2.2 times background). Ingenuity pathway analysis (IPA) revealed differences in genes involved in cell-to-cell signaling and interaction, cellular compromise and nutritional disease. In addition, we identified 100 genes that were significantly perturbed after the cholesterol-rich diet. Similar to the analysis on 10-week-old mice, IPA revealed differences in genes involved in cell-to-cell signaling and interaction as well as to immuno-inflammatory disease. Furthermore, we found important changes in genes involved in vascular calcification and matrix degradation, some of which were confirmed at protein level by (immuno-)histochemistry. In conclusion, we suggest that Cx37 deficiency alters the global differential gene expression profiles in young mice towards a pro-inflammatory phenotype, which are then further influenced in advanced atherosclerosis. The results provide new insights into the significance of Cx37 in plaque calcification.

A Novel Volume-Age-KPS (VAK) Glioblastoma Classification Identifies a Prognostic Cognate microRNA-Gene Signature.

BACKGROUND: Several studies have established Glioblastoma Multiforme (GBM) prognostic and predictive models based on age and Karnofsky Performance Status (KPS), while very few studies evaluated the prognostic and predictive significance of preoperative MR-imaging. However, to date, there is no simple preoperative GBM classification that also correlates with a highly prognostic genomic signature. Thus, we present for the first time a biologically relevant, and clinically applicable tumor Volume, patient Age, and KPS (VAK) GBM classification that can easily and non-invasively be determined upon patient admission. METHODS: We quantitatively analyzed the volumes of 78 GBM patient MRIs present in The Cancer Imaging Archive (TCIA) corresponding to patients in The Cancer Genome Atlas (TCGA) with VAK annotation. The variables were then combined using a simple 3-point scoring system to form the VAK classification. A validation set (N = 64) from both the TCGA and Rembrandt databases was used to confirm the classification. Transcription factor and genomic correlations were performed using the gene pattern suite and Ingenuity Pathway Analysis. RESULTS: VAK-A and VAK-B classes showed significant median survival differences in discovery (P = 0.007) and validation sets (P = 0.008). VAK-A is significantly associated with P53 activation, while VAK-B shows significant P53 inhibition. Furthermore, a molecular gene signature comprised of a total of 25 genes and microRNAs was significantly associated with the classes and predicted survival in an independent validation set (P = 0.001). A favorable MGMT promoter methylation status resulted in a 10.5 months additional survival benefit for VAK-A compared to VAK-B patients. CONCLUSIONS: The non-invasively determined VAK classification with its implication of VAK-specific molecular regulatory networks, can serve as a very robust initial prognostic tool, clinical trial selection criteria, and important step toward the refinement of genomics-based personalized therapy for GBM patients.

Physiology of GLP-1--lessons from glucoincretin receptor knockout mice.

GLP-1 has both peripheral and central actions, as this hormone is secreted by gut endocrine cells and brainstem neurons projecting into the hypothalamus and other brain regions. GLP-1 has multiple regulatory functions participating in the control of glucose homeostasis, beta-cell proliferation and differentiation, food intake, heart rate and even learning. GLP-1 action depends on binding to a specific G-coupled receptor linked to activation of the adenylyl cyclase pathway. Analysis of mice with inactivation of the GLP-1 receptor gene has provided evidence that absence of GLP-1 action in the mouse, despite this hormone potent physiological effects when administered in vivo, only leads to mild abnormalities in glucose homeostasis without any change in body weight. However, a critical role for this hormone and its receptor was demonstrated in the function of the hepatoportal vein glucose sensor, in contrast to that of the pancreatic beta-cells, although absence of both GLP-1 and GIP receptors leads to a more severe phenotype characterized by a beta-cell-autonomous defect in glucose-stimulated insulin secretion. Together, the studies of these glucoincretin receptor knockout mice provide evidence that these hormones are part of complex regulatory systems where multiple redundant signals are involved.

Methylation signature of lymph node metastases in breast cancer patients.

BACKGROUND: Invasion and metastasis are two important hallmarks of malignant tumors caused by complex genetic and epigenetic alterations. The present study investigated the contribution of aberrant methylation profiles of cancer related genes, APC, BIN1, BMP6, BRCA1, CST6, ESR-b, GSTP1, P14 (ARF), P16 (CDKN2A), P21 (CDKN1A), PTEN, and TIMP3, in the matched axillary lymph node metastasis in comparison to the primary tumor tissue and the adjacent normal tissue from the same breast cancer patients to identify the potential of candidate genes methylation as metastatic markers. METHODS: The quantitative methylation analysis was performed using the SEQUENOM's EpiTYPER? assay which relies on matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS). RESULTS: The quantitative DNA methylation analysis of the candidate genes showed higher methylation proportion in the primary tumor tissue than that of the matched normal tissue and the differences were significant for the APC, BIN1, BMP6, BRCA1, CST6, ESR-b, P16, PTEN and TIMP3 promoter regions (P<0.05). Among those candidate methylated genes, APC, BMP6, BRCA1 and P16 displayed higher methylation proportion in the matched lymph node metastasis than that found in the normal tissue (P<0.05). The pathway analysis revealed that BMP6, BRCA1 and P16 have a role in prevention of neoplasm metastasis. CONCLUSIONS: The results of the present study showed methylation heterogeneity between primary tumors and metastatic lesion. The contribution of aberrant methylation alterations of BMP6, BRCA1 and P16 genes in lymph node metastasis might provide a further clue to establish useful biomarkers for screening metastasis.

Modulation of genetic associations with serum urate levels by body-mass-index in humans.

We tested for interactions between body mass index (BMI) and common genetic variants affecting serum urate levels, genome-wide, in up to 42569 participants. Both stratified genome-wide association (GWAS) analyses, in lean, overweight and obese individuals, and regression-type analyses in a non BMI-stratified overall sample were performed. The former did not uncover any novel locus with a major main effect, but supported modulation of effects for some known and potentially new urate loci. The latter highlighted a SNP at RBFOX3 reaching genome-wide significant level (effect size 0.014, 95% CI 0.008-0.02, Pinter= 2.6 x 10-8). Two top loci in interaction term analyses, RBFOX3 and ERO1LB-EDARADD, also displayed suggestive differences in main effect size between the lean and obese strata. All top ranking loci for urate effect differences between BMI categories were novel and most had small magnitude but opposite direction effects between strata. They include the locus RBMS1-TANK (men, Pdifflean-overweight= 4.7 x 10-8), a region that has been associated with several obesity related traits, and TSPYL5 (men, Pdifflean-overweight= 9.1 x 10-8), regulating adipocytes-produced estradiol. The top-ranking known urate loci was ABCG2, the strongest known gout risk locus, with an effect halved in obese compared to lean men (Pdifflean-obese= 2 x 10-4). Finally, pathway analysis suggested a role for N-glycan biosynthesis as a prominent urate-associated pathway in the lean stratum. These results illustrate a potentially powerful way to monitor changes occurring in obesogenic environment.

Inference of Evolutionary Forces Acting on Human Biological Pathways.

Because natural selection is likely to act on multiple genes underlying a given phenotypic trait, we study here the potential effect of ongoing and past selection on the genetic diversity of human biological pathways. We first show that genes included in gene sets are generally under stronger selective constraints than other genes and that their evolutionary response is correlated. We then introduce a new procedure to detect selection at the pathway level based on a decomposition of the classical McDonald-Kreitman test extended to multiple genes. This new test, called 2DNS, detects outlier gene sets and takes into account past demographic effects and evolutionary constraints specific to gene sets. Selective forces acting on gene sets can be easily identified by a mere visual inspection of the position of the gene sets relative to their two-dimensional null distribution. We thus find several outlier gene sets that show signals of positive, balancing, or purifying selection but also others showing an ancient relaxation of selective constraints. The principle of the 2DNS test can also be applied to other genomic contrasts. For instance, the comparison of patterns of polymorphisms private to African and non-African populations reveals that most pathways show a higher proportion of nonsynonymous mutations in non-Africans than in Africans, potentially due to different demographic histories and selective pressures.

The role of MAR elements in DNA recombination

The untargeted integration of foreign DNA into the mammalian cell genome, extensively used in gene therapy and biotechnology, remains an incompletely understood process. It is believed to be based on cellular DNA double strand break (DSB) repair machinery and to involve two major steps: i) the formation of long gene arrays (concatemers), and ii) recombination of the resulting concatemer with the genome. The main DSB repair pathways in eukaryotes include non-homologous end-joining (NHEJ), homologous recombination (HR), and microhomology-mediated end-joining (MMEJ). However, it is still not clear, which of these pathways are responsible for transgene integration. Here, we show that NHEJ is not the primary pathway used by mammalian cells in the transgene integration process, while the components of the HR pathway seem to be important for genomic integration but not concatemerization. Instead, concatemer formation appears to be mediated by a subset of the MMEJ pathway, termed synthesis-dependent MMEJ (SD-MMEJ). This mechanism also seems to be preferentially used for plasmid integration into the genome, as confirmed by the analysis of plasmid-to-genome junction sequences, which were found to display an SD-MMEJ pattern. Therefore, we propose the existence of two distinct SD-MMEJ subpathways, relying on different subsets of enzymes. One of these mechanisms appears to be responsible for concatemerization, while the other mechanism, partially dependent in HR enzymes, seems to mediate recombination with the genome. Previous studies performed by our group suggested that matrix attachment regions (MARs), which are epigenetic regulatory DNA elements that participate in the formation of chromatin boundaries and augment transcription, may mediate increased plasmid integration into the genome of CHO cells by stimulating DNA recombination. In the present work, we demonstrate that MAR-mediated plasmid integration results from the enhanced SD-MMEJ pathway. Analysis of transgene integration loci and junction DNA sequences validated the prevalent use of this pathway by the MAR elements to target plasmid DNA into gene-rich areas of the CHO genome. We propose that this finding should in the future help to engineer cells for improved recombinant protein production. In addition to investigating the process of transgene integration, we designed recombination assays to better characterize the components of the MMEJ and SD-MMEJ pathways. We also used CHO cells expressing cycle-sensitive reporter genes to demonstrate a potential role of HR proteins in the cell cycle regulation.

Further molecular profiling of tumors harboring therapeutic targets within non-small cell lung cancer

Background: Treatment of NSCLC has been revolutionized in recent years with the introduction of several targeted therapies for selected genetically altered subtypes of NSCLC. A better understanding of molecular characteristics of NSCLC, which features common drug targets, may identify new therapeutic options. Methods: Over 6,700 non-small cell lung cancer cases referred to Caris Life Sciences between 2009 and 2014. Diagnoses and history were collected from referring physicians. Specific testing was performed per physician request and included a combination of sequencing (Sanger, NGS or pyrosequencing), protein expression (IHC), gene amplification/rearrangement (CISH or FISH), and/or RNA fragment analysis. Results: Tumors profiles from patients with hormone receptor positive disease (HER2, ER, PR, or AR positive by IHC) (n=629), HER2 mutations (n=8) ALK rearrangements (n=55), ROS1 rearrangement (n=17), cMET amplification or mutation (n=126), and cKIT mutation (n=11) were included in this analysis and compared to the whole cohort. Tumors with ALK rearrangement overexpressed AR in 18% of cases, and 7% presented with concomitant KRAS mutation. Lower rates of PTEN loss, as assessed by IHC, were observed in ALK positive (20%), ROS1 positive (9%) and cKIT mutated tumors (25%) compared to the overall NSCLC population (58%). cMET was overexpressed in 66% of ROS1 translocated and 57% of HER2 mutated tumors. cKIT mutations were found co-existing with APC (20%) and EGFR (20%) mutations. Pathway analysis revealed that hormone receptor positive disease carried more mutations in the ERK pathway (32%) compared to 9% in the mTOR pathway. 25% of patients with HER2 mutations harbored a co-existing mutation in the mTOR pathway. Conclusions: Pathway profiling reveals that NSCLC tumors present more often than reported with several concomitant alterations affecting the ERK or AKT pathway. Additionally, they are also characterized by the expression of potential biological modifiers of the cell cycle like hormonal receptors, representing a rationale for dual inhibition strategies in selected patients. Further refining of the understanding of NSCLC biomarker profile will optimize research for new treatment strategies.

Molecular profiling of non-small cell lung cancer by histologic subtype

Background: A substantial proportion of NSCLC has been shown to harbour specific molecular alterations affecting tumour proliferation and resulting in sensitivity to inhibition of the corresponding activated oncogenic pathway by targeted therapies. Comprehensive tumor profiling can diagnose such alterations and may identify new alterations opening additional treatment options for all distinct NSCLC subtypes. Methods: Over 6,700 non-small cell lung cancer cases referred to Caris Life Sciences between 2009 and 2014 were evaluated; clinical diagnoses and detailed tumor pathology were collected from referring physicians. Specific profiling was performed per physician request and included a combination of sequencing (Sanger, NGS or pyrosequencing), protein expression (IHC), gene amplification/rearrangement (CISH or FISH), and/or RNA fragment analysis within potential cancer-related genes and pathways. Results: Patients were grouped into cohorts according to histological subtype - adenocarcinoma (AD) (n=4,286), squamous cell carcinoma (SCC) (n=1,280), large cell carcinoma (LCC) (n=153) and bronchioalveolar carcinoma (BAC) (n=94). Protein overexpression of cMET (>2+ in >50% cells) was higher in AD (35.9%) compared to other subgroups (12-20%) while RRM1 and TOP2A levels were lower in AD. ALK or ROS1 were rearranged in 5.3% of patients with AD compared to 3.7% of patients with LCC and 1.2% of patients with SCC. EGFR mutations were found at low prevalence in both the LCC (0%) and SCC cohorts (2.8%) compared to 21% in AD. Similar lower rates of BRAF mutations were observed in the LCC and SCC cohorts compared to AD (0%, 1.1% and 5.1%). Pathway analysis showed activating mutations in the ERK pathway in 40% of patients with AD. Only 10-12% of patients with LCC or SCC had activating mutations in the ERK pathway. Conclusions: Despite the limitations of this retrospective series, we report comprehensive profiling of the largest cohort of NSCLC. Tumor profiling reveals that ADs may be more addicted to the ERK pathway than other histological subtypes. Drugs which target cMET may also have most utility in AD. Full analysis by histological subtype and additional correlative data on protein expression, gene copy number and mutations will be presented.