Philippe Hupé

Analysis of array CGH data: from signal ratio to gain and loss of DNA regions

MOTIVATION Genomic DNA regions are frequently lost or gained during tumor progression. Array Comparative Genomic Hybridization (array CGH) technology makes it possible to assess these changes in DNA in cancers, by comparison with a normal reference. The identification of systematically deleted or amplified genomic regions in a set of tumors enables biologists to identify genes involved in cancer progression because tumor suppressor genes are thought to be located in lost genomic regions and oncogenes, in gained regions. Array CGH profiles should also improve the classification of tumors. The achievement of these goals requires a methodology for detecting the breakpoints delimiting altered regions in genomic patterns and assigning a status (normal, gained or lost) to each chromosomal region. RESULTS We have developed a methodology for the automatic detection of breakpoints from array CGH profile, and the assignment of a status to each chromosomal region. The breakpoint detection step is based on the Adaptive Weights Smoothing (AWS) procedure and provides highly convincing results: our algorithm detects 97, 100 and 94% of breakpoints in simulated data, karyotyping results and manually analyzed profiles, respectively. The percentage of correctly assigned statuses ranges from 98.9 to 99.8% for simulated data and is 100% for karyotyping results. Our algorithm also outperforms other solutions on a public reference dataset. AVAILABILITY The R package GLAD (Gain and Loss Analysis of DNA) is available upon request.

miR-181a and miR-630 regulate cisplatin-induced cancer cell death.

MicroRNAs (miRNA) are noncoding RNAs that regulate multiple cellular processes, including proliferation and apoptosis. We used microarray technology to identify miRNAs that were upregulated by non-small cell lung cancer (NSCLC) A549 cells in response to cisplatin (CDDP). The corresponding synthetic miRNA precursors (pre-miRNAs) per se were not lethal when transfected into A549 cells yet affected cell death induction by CDDP, C2-ceramide, cadmium, etoposide, and mitoxantrone in an inducer-specific fashion. Whereas synthetic miRNA inhibitors (anti-miRNAs) targeting miR-181a and miR-630 failed to modulate the response of A549 to CDDP, pre-miR-181a and pre-miR-630 enhanced and reduced CDDP-triggered cell death, respectively. Pre-miR-181a and pre-miR-630 consistently modulated mitochondrial/postmitochondrial steps of the intrinsic pathway of apoptosis, including Bax oligomerization, mitochondrial transmembrane potential dissipation, and the proteolytic maturation of caspase-9 and caspase-3. In addition, pre-miR-630 blocked early manifestations of the DNA damage response, including the phosphorylation of the ataxia-telangiectasia mutated (ATM) kinase and of two ATM substrates, histone H2AX and p53. Pharmacologic and genetic inhibition of p53 corroborated the hypothesis that pre-miR-630 (but not pre-miR-181a) blocks the upstream signaling pathways that are ignited by DNA damage and converge on p53 activation. Pre-miR-630 arrested A549 cells in the G0-G1 phase of the cell cycle, correlating with increased levels of the cell cycle inhibitor p27(Kip1) as well as with reduced proliferation rates and resulting in greatly diminished sensitivity of A549 cells to the late S-G2-M cell cycle arrest mediated by CDDP. Altogether, these results identify miR-181a and miR-630 as novel modulators of the CDDP response in NSCLC.

Frequent PTEN genomic alterations and activated phosphatidylinositol 3-kinase pathway in basal-like breast cancer cells

IntroductionBasal-like carcinomas (BLCs) and human epidermal growth factor receptor 2 overexpressing (HER2+) carcinomas are the subgroups of breast cancers that have the most aggressive clinical behaviour. In contrast to HER2+ carcinomas, no targeted therapy is currently available for the treatment of patients with BLCs. In order to discover potential therapeutic targets, we aimed to discover deregulated signalling pathways in human BLCs.MethodsIn this study, we focused on the oncogenic phosphatidylinositol 3-kinase (PI3K) pathway in 13 BLCs, and compared it with a control series of 11 hormonal receptor negative- and grade III-matched HER2+ carcinomas. The two tumour populations were first characterised by immunohistochemistry and gene expression. The PI3K pathway was then investigated by gene copy-number analysis, gene expression profiling and at a proteomic level using reverse-phase protein array technology and tissue microarray. The effects of the PI3K inhibition pathway on proliferation and apoptosis was further analysed in three human basal-like cell lines.ResultsThe PI3K pathway was found to be activated in BLCs and up-regulated compared with HER2+ tumours as shown by a significantly increased activation of the downstream targets Akt and mTOR (mammalian target of rapamycin). BLCs expressed significantly lower levels of the tumour suppressor PTEN and PTEN levels were significantly negatively correlated with Akt activity within that population. PTEN protein expression correlated significantly with PTEN DNA copy number and more importantly, reduced PTEN DNA copy numbers were observed specifically in BLCs. Similar to human samples, basal-like cell lines exhibited an activation of PI3K/Akt pathway and low/lack PTEN expression. Both PI3K and mTOR inhibitors led to basal-like cell growth arrest. However, apoptosis was specifically observed after PI3K inhibition.ConclusionsThese data provide insight into the molecular pathogenesis of BLCs and implicate the PTEN-dependent activated Akt signalling pathway as a potential therapeutic target for the management of patients with poor prognosis BLCs.

VAMP: Visualization and analysis of array-CGH, transcriptome and other molecular profiles

MOTIVATION Microarray-based CGH (Comparative Genomic Hybridization), transcriptome arrays and other large-scale genomic technologies are now routinely used to generate a vast amount of genomic profiles. Exploratory analysis of this data is crucial in helping to understand the data and to help form biological hypotheses. This step requires visualization of the data in a meaningful way to visualize the results and to perform first level analyses. RESULTS We have developed a graphical user interface for visualization and first level analysis of molecular profiles. It is currently in use at the Institut Curie for cancer research projects involving CGH arrays, transcriptome arrays, SNP (single nucleotide polymorphism) arrays, loss of heterozygosity results (LOH), and Chromatin ImmunoPrecipitation arrays (ChIP chips). The interface offers the possibility of studying these different types of information in a consistent way. Several views are proposed, such as the classical CGH karyotype view or genome-wide multi-tumor comparison. Many functionalities for analyzing CGH data are provided by the interface, including looking for recurrent regions of alterations, confrontation to transcriptome data or clinical information, and clustering. Our tool consists of PHP scripts and of an applet written in Java. It can be run on public datasets at http://bioinfo.curie.fr/vamp AVAILABILITY The VAMP software (Visualization and Analysis of array-CGH,transcriptome and other Molecular Profiles) is available upon request. It can be tested on public datasets at http://bioinfo.curie.fr/vamp. The documentation is available at http://bioinfo.curie.fr/vamp/doc.

Genomic Profiling and Identification of High-Risk Uveal Melanoma by Array CGH Analysis of Primary Tumors and Liver Metastases

PURPOSE Incurable metastases develop in approximately 50% of patients with uveal melanoma (UM). The purpose of this study was to analyze genomic profiles in a large series of ocular tumors and liver metastases and design a genome-based classifier for metastatic risk assessment. METHODS A series of 86 UM tumors and 66 liver metastases were analyzed by using a BAC CGH (comparative genomic hybridization) microarray. A clustering was performed, and correlation with the metastatic status was sought among a subset of 71 patients with a minimum follow-up of 24 months. The status of chromosome 3 was further examined in the tumors, and metastases with disomy 3 were checked with an SNP microarray. A prognostic classifier was constructed using a log-linear model on minimal regions and leave-one-out cross-validation. RESULTS The clustering divides the groups of tumors with disomy 3 and monosomy 3 into two and three subgroups, respectively. Same subgroups are found in primary tumors and in metastases, but with different frequencies. Isolated monosomy 3 was present in 0% of metastatic ocular tumors and in 3% of metastases. The highest metastatic rate in ocular tumors was observed in a subgroup defined by the gain of 8q with a proximal breakpoint, and losses of 3, 8p, and 16q, also most represented in metastases. A prognostic classifier that included the status of these markers led to an 85.9% classification accuracy. CONCLUSIONS The analysis of the status of these specific chromosome regions by genome profiling on SNP microarrays should be a reliable tool for identifying high-risk patients in future adjuvant therapy protocols.

Live-Cell Chromosome Dynamics and Outcome of X Chromosome Pairing Events during ES Cell Differentiation

Random X inactivation represents a paradigm for monoallelic gene regulation during early ES cell differentiation. In mice, the choice of X chromosome to inactivate in XX cells is ensured by monoallelic regulation of Xist RNA via its antisense transcription unit Tsix/Xite. Homologous pairing events have been proposed to underlie asymmetric Tsix expression, but direct evidence has been lacking owing to their dynamic and transient nature. Here we investigate the live-cell dynamics and outcome of Tsix pairing in differentiating mouse ES cells. We find an overall increase in genome dynamics including the Xics during early differentiation. During pairing, however, Xic loci show markedly reduced movements. Upon separation, Tsix expression becomes transiently monoallelic, providing a window of opportunity for monoallelic Xist upregulation. Our findings reveal the spatiotemporal choreography of the X chromosomes during early differentiation and indicate a direct role for pairing in facilitating symmetry-breaking and monoallelic regulation of Xist during random X inactivation.

Computation of recurrent minimal genomic alterations from array-CGH data

MOTIVATION The identification of recurrent genomic alterations can provide insight into the initiation and progression of genetic diseases, such as cancer. Array-CGH can identify chromosomal regions that have been gained or lost, with a resolution of approximately 1 mb, for the cutting-edge techniques. The extraction of discrete profiles from raw array-CGH data has been studied extensively, but subsequent steps in the analysis require flexible, efficient algorithms, particularly if the number of available profiles exceeds a few tens or the number of array probes exceeds a few thousands. RESULTS We propose two algorithms for computing minimal and minimal constrained regions of gain and loss from discretized CGH profiles. The second of these algorithms can handle additional constraints describing relevant regions of copy number change. We have validated these algorithms on two public array-CGH datasets. AVAILABILITY From the authors, upon request. CONTACT celine@lri.fr SUPPLEMENTARY INFORMATION Supplementary data are available at Bioinformatics online.

Spatial normalization of array-CGH data

BackgroundArray-based comparative genomic hybridization (array-CGH) is a recently developed technique for analyzing changes in DNA copy number. As in all microarray analyses, normalization is required to correct for experimental artifacts while preserving the true biological signal. We investigated various sources of systematic variation in array-CGH data and identified two distinct types of spatial effect of no biological relevance as the predominant experimental artifacts: continuous spatial gradients and local spatial bias. Local spatial bias affects a large proportion of arrays, and has not previously been considered in array-CGH experiments.ResultsWe show that existing normalization techniques do not correct these spatial effects properly. We therefore developed an automatic method for the spatial normalization of array-CGH data. This method makes it possible to delineate and to eliminate and/or correct areas affected by spatial bias. It is based on the combination of a spatial segmentation algorithm called NEM (Neighborhood Expectation Maximization) and spatial trend estimation. We defined quality criteria for array-CGH data, demonstrating significant improvements in data quality with our method for three data sets coming from two different platforms (198, 175 and 26 BAC-arrays).ConclusionWe have designed an automatic algorithm for the spatial normalization of BAC CGH-array data, preventing the misinterpretation of experimental artifacts as biologically relevant outliers in the genomic profile. This algorithm is implemented in the R package MANOR (Micro-Array NORmalization), which is described at http://bioinfo.curie.fr/projects/manor and available from the Bioconductor site http://www.bioconductor.org. It can also be tested on the CAPweb bioinformatics platform at http://bioinfo.curie.fr/CAPweb.

Multiparametric analysis of cytokine-driven human Th17 differentiation reveals a differential regulation of IL-17 and IL-22 production

T helper 17 (Th17) cells produce IL-17 but can also make tumor necrosis factor, interleukin (IL)-6, IL-10, IL-21, and IL-22. These cytokines collectively contribute to the functional outcome of the Th response. IL-22 plays a critical role in some Th17-associated diseases, such as psoriasis, but its relationship to IL-17 remains controversial. Here, we used a systematic multiparametric analysis of Th-17-associated cytokines, which revealed the unexpected finding that the regulation pattern of IL-22 was most closely related to interferon-gamma, the prototypical Th1 cytokine, and not to IL-17. To explain this observation, we systematically tested the role of Th1- and Th17-inducing cytokines. We could show that IL-12 and IL-23 induced high levels of IL-22 but no IL-17. Conversely, transforming growth factor-beta inhibited IL-22 production but promoted IL-17. Thus, IL-17 and IL-22 are differentially regulated during cytokine-induced Th cell differentiation. This has important implications for the understanding and pharmacologic manipulation of Th17-associated pathologies.

Prognostic Impact of Vitamin B6 Metabolism in Lung Cancer

Patients with non-small cell lung cancer (NSCLC) are routinely treated with cytotoxic agents such as cisplatin. Through a genome-wide siRNA-based screen, we identified vitamin B6 metabolism as a central regulator of cisplatin responses in vitro and in vivo. By aggravating a bioenergetic catastrophe that involves the depletion of intracellular glutathione, vitamin B6 exacerbates cisplatin-mediated DNA damage, thus sensitizing a large panel of cancer cell lines to apoptosis. Moreover, vitamin B6 sensitizes cancer cells to apoptosis induction by distinct types of physical and chemical stress, including multiple chemotherapeutics. This effect requires pyridoxal kinase (PDXK), the enzyme that generates the bioactive form of vitamin B6. In line with a general role of vitamin B6 in stress responses, low PDXK expression levels were found to be associated with poor disease outcome in two independent cohorts of patients with NSCLC. These results indicate that PDXK expression levels constitute a biomarker for risk stratification among patients with NSCLC.