Unknown,Transcriptomics,Genomics,Proteomics

Dataset Information

Mutational activation of a mechanosensitive megachannel drives metastatic cell survival.

ABSTRACT: Next-generation sequencing has revolutionized cancer biology by accelerating the unbiased discovery of novel mutations across human cancers. Transforming such discoveries into a conceptual framework of cancer progression requires narrowing the vast number of mutations down to the driver elements, and further reducing these to mutations that govern cancer progression as opposed to tumor initiation. By integrating next-generation RNA-sequencing (RNA-seq) with in vivo selection, we devise an approach that identifies a series of novel recurrent non-synonymous amino acid mutations that are enriched in metastatic breast cancer cells and predicted to significantly alter protein function. These mutations, found in PANX1, RBFA, REST, KRIT1 and ZSWIM6, are detected at higher frequencies in the transcriptomes of two patientsâ highly metastatic sub-lines relative to their poorly metastatic parental lines. We functionally characterize the cellular and molecular roles of one of these mutationsâa nonsense alteration that yields a truncated pannexin-1 (PANX11-89) plasma membrane megachannel subunitâin metastatic progression. PANX11-89 interacts with full-length PANX1 and augments PANX1 channel activity to promote the survival of cancer cells as they are mechanically deformed. Protection from deformation-induced cell death requires PANX1 channels to release ATP, which acts as a cell autonomous survival signal during mechanical stress. Functional characterization of additional nonsense and missense PANX1 mutations detected in epithelial cancers of the colon, lung, and prostate reveals that these mutants also enhance PANX1-mediated ATP release. In vivo testing of one such truncating mutation detected in a metastatic colorectal tumor also enhances early survival, dissemination and liver metastatic colonization by human colon cancer cells. Finally, pharmacological inhibition of PANX1 inhibits breast cancer metastasis, implicating PANX1 as a novel therapeutic target in cancer. Our findings reveal that ATP release through mechanosensitive PANX1 channels enables cancer cells to overcome a major metastasis suppressive barrier—deformation-induced death in the microvasculature. To systematically identify base-pair mutations present in metastatic cells that may drive cancer progression, we performed whole-transcriptome RNA-sequencing (RNA-seq) of in vivo-selected, highly metastatic human breast cancer cell sub-lines, CN-LM1A and MDA-LM2, as well as the CN34 and MDA-MB-231 parental lines from which they were derived. To minimize the false positive rate and allow for subsequent statistical analysis, we sequenced biological replicates of each cell line. Mutations conferring enhanced metastatic capacity should be enriched in the transcriptomes of highly metastatic cells relative to their less metastatic parental populations.

ORGANISM(S): Homo sapiens

SUBMITTER: David Soong

PROVIDER: E-GEOD-45162 | biostudies-arrayexpress |

REPOSITORIES: biostudies-arrayexpress

ACCESS DATA

Similar Datasets

Project description:Next-generation sequencing has revolutionized cancer biology by accelerating the unbiased discovery of novel mutations across human cancers. Transforming such discoveries into a conceptual framework of cancer progression requires narrowing the vast number of mutations down to the driver elements, and further reducing these to mutations that govern cancer progression as opposed to tumor initiation. By integrating next-generation RNA-sequencing (RNA-seq) with in vivo selection, we devise an approach that identifies a series of novel recurrent non-synonymous amino acid mutations that are enriched in metastatic breast cancer cells and predicted to significantly alter protein function. These mutations, found in PANX1, RBFA, REST, KRIT1 and ZSWIM6, are detected at higher frequencies in the transcriptomes of two patients’ highly metastatic sub-lines relative to their poorly metastatic parental lines. We functionally characterize the cellular and molecular roles of one of these mutations—a nonsense alteration that yields a truncated pannexin-1 (PANX11-89) plasma membrane megachannel subunit—in metastatic progression. PANX11-89 interacts with full-length PANX1 and augments PANX1 channel activity to promote the survival of cancer cells as they are mechanically deformed. Protection from deformation-induced cell death requires PANX1 channels to release ATP, which acts as a cell autonomous survival signal during mechanical stress. Functional characterization of additional nonsense and missense PANX1 mutations detected in epithelial cancers of the colon, lung, and prostate reveals that these mutants also enhance PANX1-mediated ATP release. In vivo testing of one such truncating mutation detected in a metastatic colorectal tumor also enhances early survival, dissemination and liver metastatic colonization by human colon cancer cells. Finally, pharmacological inhibition of PANX1 inhibits breast cancer metastasis, implicating PANX1 as a novel therapeutic target in cancer. Our findings reveal that ATP release through mechanosensitive PANX1 channels enables cancer cells to overcome a major metastasis suppressive barrier—deformation-induced death in the microvasculature.

Project description:In-depth information regarding the DFCI OncoPanel sequencing panel has been described previously. Briefly, sequencing is performed using an Illumina HiSeq 2500 system (RRID:SCR_016383) with 2×100 paired-end reads. Samples must meet an average 50X coverage and minimum of 30X coverage for 80% of targets for analysis. For all samples, a board-certified molecular pathologist evaluates alterations and writes interpretations of the results. The three generations of OncoPanel use Mutect (RRID:SCR_000559) and GATK (RRID:SCR_001876) to interrogate possible alterations in the complete exonic DNA sequences of 275, 309, and 447 cancer-related genes, respectively, including substitutions, insertions, and deletions. Because tumor tissues are tested without a paired normal from individual patients, additional informatics steps are taken to identify common single nucleotide polymorphisms (SNPs). Any SNP present at >0.1% in Exome Variant Server, NHLBI GO Exome Sequencing Project (ESP) (RRID:SCR_012761), or present in dbSNP (RRID:SCR_002338) is filtered; however, variants also present in at least twice in COSMIC (RRID:SCR_002260) are rescued for manual review. Variants that appear two or more times in a panel of 150 normal samples sequenced in-house and are not present in COSMIC are also filtered. For copy number alteration (CNA) analysis, the workflow employs an in-house algorithm (RobustCNV) to cover exonic regions of targeted genes as well as select intronic regions. The depth of coverage of these regions is determined by counting the number of reads aligning within defined genomic intervals and then normalized against a panel of normals and corrected for GC bias. The copy number for a segmented genomic interval is calculated as a log2 ratio of the depth of coverage of this sample compared to a panel of normal (non-cancer) samples that are run on the same plate. In parallel (samples not included), additional cases were previously sequenced via the MSK-IMPACT platform. To ensure maximum overlap between the platforms, only single nucleotide variants (SNVs), dinucleotide/oligonucleotide variants (DNVs/ONVs), insertions/deletions (indels), high copy gains (amplifications), and homozygous deletions from 182 genes shared between all versions of both platforms were evaluated in our analyses involving samples from both platforms. To filter pathogenic and passenger mutations, missense, truncating (nonsense, frameshift, and splice site), and in-frame indel mutations from the 182 genes overlapping between the DFCI OncoPanel and MSK-IMPACT panels were subjected to analysis in the Cancer Genome Interpreter (CGI) suite (RRID:SCR_023752) with cancer type set to “Skin Cutaneous Melanoma (SKCM)”. Known pathogenic mutations, defined as those listed as pathogenic in ClinVar (RRID:SCR_006169), OncoKB (RRID:SCR_014782), and/or CGI databases, were included in downstream analyses. Mutations which had not been previously annotated but were predicted to be pathogenic mutations by CGI were also included. Missense SNV mutations called passengers in CGI were subjected to another round of analysis in the ChasmPlus suite with cancer type set to SKCM. Missense mutations called pathogenic via this approach (p < 0.05; FDR Q < 0.3) were salvaged and included in further analyses. All other mutations were excluded. TERT promoter mutations were evaluated via FATHMM-MKL. Mutations called pathogenic via this tool were included in downstream analyses. Clinical data for these cases is available in a supplementary file from the associated publication. OF NOTE, THE RAW SEQUENCING DATA FROM THESE SAMPLES CANNOT BE MADE PUBLICLY AVAILABLE BECAUSE THE RESEARCH PARTICIPANT CONSENT DOES NOT INCLUDE AUTHORIZATION TO SHARE IDENTIFIABLE DATA.

Project description:U87MG is a commonly studied grade IV glioma cell line that has been analyzed in at least 1,700 publications over four decades. In order to comprehensively characterize the genome of this cell line and to serve as a model of broad cancer genome sequencing, we have generated greater than 30x genomic sequence coverage using a novel 50-base mate paired strategy with a 1.4kb mean insert library. A total of 1,014,984,286 mate-end and 120,691,623 single-end two-base encoded reads were generated from five slides. All data were aligned using a custom designed tool called BFAST, allowing optimal color space read alignment and accurate identification of DNA variants. The aligned sequence reads and mate pair information identified 35 interchromosomal translocation events, 1,315 structural variations (>100bp), 191,743 small (<21bp) insertions and deletions (indels), and 2,384,470 single nucleotide variations (SNVs). Among these observations, the known homozygous mutation in PTEN was robustly identified, and genes involved in cell adhesion were overrepresented in the mutated gene list. Data were compared to 219,187 heterozygous single nucleotide polymorphisms assayed by Illumina 1M Duo genotyping array to assess accuracy: 93.83% of all SNPs were reliably detected at filtering thresholds that yield greater than 99.99% sequence accuracy. Protein coding sequences were disrupted predominantly in this cancer cell line due to small indels, large deletions and translocations. In total, 512 genes were homozygously mutated, including 154 by SNVs, 178 by small indels, 145 by large microdeletions and 35 by interchromosomal translocations to reveal a highly mutated cell line genome. Of the small homozygously mutated variants, 8 SNVs and 99 indels were novel events not present in dbSNP. These data demonstrate that routine generation of broad cancer genome sequence is possible outside of genome centers. The sequence analysis of U87MG provides an unparalleled level of mutational resolution compared to any cell line to date. Whole genome sequencing of the U87MG brain cancer cell line using the AB SOLiD3 sequencer and genotyping using the Illumina Human1M-Duov3 DNA Analysis BeadChip

Dataset Information

Mutational activation of a mechanosensitive megachannel drives metastatic cell survival.

Similar Datasets

OmicsDI is part of the ELIXIR infrastructure

Tweets