First somatic mutation of E2F1 in a critical DNA binding residue discovered in well-differentiated papillary mesothelioma of the peritoneum
- Willie Yu†1, 2, 3,
- Waraporn Chan-On†1, 2,
- Melissa Teo4,
- Choon Kiat Ong1, 2,
- Ioana Cutcutache5,
- George E Allen1, 2,
- Bernice Wong1, 2,
- Swe Swe Myint1, 2,
- Kiat Hon Lim6,
- P Mathijs Voorhoeve7, 8,
- Steve Rozen5,
- Khee Chee Soo4,
- Patrick Tan9, 10, 11Email author and
- Bin Tean Teh1, 2, 12Email author
© Yu et al.; licensee BioMed Central Ltd. 2011
Received: 13 March 2011
Accepted: 28 September 2011
Published: 28 September 2011
Well differentiated papillary mesothelioma of the peritoneum (WDPMP) is a rare variant of epithelial mesothelioma of low malignancy potential, usually found in women with no history of asbestos exposure. In this study, we perform the first exome sequencing of WDPMP.
WDPMP exome sequencing reveals the first somatic mutation of E2F1, R166H, to be identified in human cancer. The location is in the evolutionarily conserved DNA binding domain and computationally predicted to be mutated in the critical contact point between E2F1 and its DNA target. We show that the R166H mutation abrogates E2F1's DNA binding ability and is associated with reduced activation of E2F1 downstream target genes. Mutant E2F1 proteins are also observed in higher quantities when compared with wild-type E2F1 protein levels and the mutant protein's resistance to degradation was found to be the cause of its accumulation within mutant over-expressing cells. Cells over-expressing wild-type E2F1 show decreased proliferation compared to mutant over-expressing cells, but cell proliferation rates of mutant over-expressing cells were comparable to cells over-expressing the empty vector.
The R166H mutation in E2F1 is shown to have a deleterious effect on its DNA binding ability as well as increasing its stability and subsequent accumulation in R166H mutant cells. Based on the results, two compatible theories can be formed: R166H mutation appears to allow for protein over-expression while minimizing the apoptotic consequence and the R166H mutation may behave similarly to SV40 large T antigen, inhibiting tumor suppressive functions of retinoblastoma protein 1.
Mesothelioma is an uncommon neoplasm that develops from the mesothelium, the protective lining covering a majority of the body's internal organs, and is divided into four subtypes: pleural, peritoneum, pericardium and tunica vaginalis . While malignant peritoneal mesothelioma (MPM) is an aggressive tumor mainly afflicting asbestos-exposed males in the age range of 50 to 60 years old , well-differentiated papillary mesothelioma of the peritoneum (WDPMP), a rare subtype of epithelioid mesothelioma  with fewer than 60 cases described in the literature , is generally considered to be a tumor of low malignant potential found predominately in young women with no definitive exposure to asbestos . While much scientific research has been done on asbestos-related malignant mesothelioma [4–7], the rarity of WDPMP coupled with its good prognosis relegated its research to case reports and reviews by medical oncologists concentrating in the area of diagnosis, prognosis and treatment options.
Second generation sequencing technologies coupled with newly developed whole exome capturing technologies  allow for rapid, relatively inexpensive approaches to obtain an overview of large complex genomes by concentrating on the critical coding areas of the genome. Here, we report the first exome sequencing of a WDPMP tumor, its tumor-derived cell line and a matched control sample employing Agilent SureSelect All Exon capturing technology to selectively capture all human exons followed by Illumina massively parallel genomic sequencing. We developed methodology and informatics to obtain a compact graphical view of the exome as well as detailed analysis of single nucleotide variants (SNVs). We demonstrate that while this WDPMP tumor does not exhibit any of the chromosomal aberrations and focal deletions commonly associated with asbestos-related mesothelioma , it does exhibit the first reported somatic single nucleotide mutation of E2F1 (E2F transcription factor 1) in cancer, with the mutation affecting one of two evolutionarily conserved arginine residues responsible for motif recognition and DNA binding.
WDPMP exome sequencing: mutation landscape changes big and small
Summary of overall WDPMP exome sequencing
Reads passing filter
Percentage remaining after filter
Number of PCR duplicates
Percentage PCR duplicates
Reads within exome
Percentage reads overall in exome
Percentage exome covered at 1×
Percentage exome covered at 20×
Since the Hilbert plots showed no gross anomalies, we turned our attention to mining the exome data for somatic single nucleotide mutations. The SNV discovery pipeline, described in the Materials and methods section, was performed using the Genome Analyzer Toolkit  for the tumor, normal and cell line exomes. Filtering was set to accept candidate SNVs with a quality/depth score > 3 and were present in both the tumor and cell line but not in the normal sample. Nineteen potential somatic mutations remain and these were validated using Sanger sequencing (Additional file 2); E2F1, PPFIBP2 (liprin beta 2) and TRAF7 (TNF receptor-associated factor 7) were validated to contain true somatic mutations (Additional file 3).
The E2F1 R166H mutation affects a critical DNA binding residue
The R166H mutation is detrimental to E2F1's DNA binding ability and negatively affects downstream target gene expression
To show that the R166H mutant's reduced DNA binding affinity affected the expression of E2F1 target genes, the expression of SIRT1, APAF1 and CCNE1 (cyclin E1) was examined by real-time PCR in MSTO-211H and NCI-H28 cell lines that were transfected with the E2F1 mutant or wild type. Interestingly, over-expression of the E2F1 R166H mutant (E2F1-R166H) did not up-regulate expression of SIRT1 and APAF1 as high as over-expression of wild-type E2F1 in both cell lines (Figure 6b, c). In particular, levels of SIRT1 and APAF1 expression in MSTO-211H cells observed with E2F1-R166H were significantly lower than those with the E2F1 wild type (P = 0.032 for SIRT1 and P = 0.005 for APAF1). However, the expression of cyclin E1, a well known target of E2F1 , was minimally affected in the over-expression context, which may be indicative of a compensatory effect by other members of the E2F family.
Cells over-expressing E2F1-R166H show massive protein accumulation and increased protein stability
To monitor E2F1 protein stability, we over-expressed the E2F1 wild type and mutant in MSTO-211H cells before treating the cells with 25 μg/ml cyclohexamide to block new protein synthesis in half hour intervals. As shown in Figure 6b, the protein levels of the E2F1 mutant remained almost constant throughout the 3-hour period of the experiment while those of the wild type decreased in a time-dependent manner. This result suggests that the mutant protein is more stable and resistant to degradation than the wild type and the increased stability of E2F1-R166H is the cause of its accumulation within cells over-expressing it.
Over-expression of E2F1-R166H does not adversely affect cell proliferation
For this study we have performed the first exome sequencing of a matched pair of WDPMP samples along with a cell line derived from the tumor. Analysis of the exomes revealed none of the chromosomal aberrations or focal gene deletions commonly associated with asbestos-related malignant mesothelioma. We were able to verify somatic mutations in PPFIBP2, TRAF7 and E2F1.
TRAF7 is an E3 ubiquitin ligase  shown to be involved in MEKK3 (mitogen-activated protein kinase kinase kinase 3) signaling and apoptosis . The Y621D mutation in TRAF7 occurs in the WD40 repeat domain, which has been shown to be involved in MEKK3-induced activator protein 1 (AP1) activation . Since AP1 in turn controls a large number of cellular processes involved in differentiation, proliferation and apoptosis , this mutation in TRAF7's WD40 repeat domain may de-regulate MEKK3's control over AP1 activation, which may contribute to WDPMP transformation.
PPFIBP2 is a member of the LAR protein-tyrosine-phosphatase-interacting protein (liprin) family . While no functional studies on PPFIBP2 have been published, it was reported as a potential biomarker for endometrial carcinomas . However, the Q791H mutation in PPFIBP2 is predicted by Polyphen to be benign and the Catalogue of Somatic Mutations in Cancer (COSMIC) did not show this particular mutation to occur in other cancers; thus, this mutation is likely to be of a passenger variety.
Of particular interest is the E2F1 mutation as no reported somatic mutation has ever been reported for this protein despite its critical roles in cell cycle control , apoptosis  and DNA repair . Using various bioinformatics tools, this mutation was identified to mutate an arginine residue into a histidine residue, thus altering a critical evolutionarily conserved DNA contact point responsible for DNA binding and motif recognition.
Since computational modeling is sufficient to pinpoint the mutation's structural location but is inconclusive in showing the mutation's functional effect on DNA binding, we performed a ChIP assay that showed the R166H mutation abrogates E2F1 DNA binding. Analysis of expression of selected E2F1 target genes in an over-expression system showed the inability of the E2F1 mutant to adequately up-regulate the expression of SIRT1 and APAF1 when compared with the E2F1 wild type. Of interest is the lack of changes in the expression of the gene encoding cyclin E1, a known target of E2F1 and an important component in starting the S phase of the cell cycle. A possible explanation for this is the functional redundancy of the E2F family to ensure the cell's replication machinery is operational - for example, mice studies have shown that E2F1-/- mice can be grown to maturity [29, 30].
Our study has also shown that the R166H mutant is much more stable than its wild-type counterpart, enabling massive accumulation within the cell. A previous study showed that over-expression of E2F1 results in induction of apoptosis , which is in line with our observation of a drop in proliferation when cells were over-expressing wild-type E2F1; curiously, over-expression of the mutant E2F1 protein did not lead to any noticeable effect on cellular proliferation even though mutant protein levels were many fold higher than those of the wild type in equivalent transfection conditions. One explanation for this phenomenon is that inactivation of E2F1 decreases apoptosis and its abrogated cell cycle role is compensated for by other members of its family. E2F1-/- mice can grow to maturity and reproduce normally but display a predisposition to develop various cancers , indicating the greater importance of the tumor suppressive function of E2F1 compared to its cell cycle gene activation function.
An alternative but not mutually exclusive explanation is that stable and numerous E2F1-R166H proteins behave functionally like SV40 large T antigens, serving as competitive inhibitors by taking up the lion's share of the binding capacity of retinoblastoma protein 1 (Rb), resulting in unbound wild-type E2F1, which drives the cell cycle. While the R166H mutation cripples E2F1's DNA binding ability, its other interaction domains, including the Rb interaction domain, are still active. The mutant's stability and large quantities will favor its preferential binding to Rb due to its sheer numbers and the heterozygous nature of the mutation in the WDPMP tumor would ensure active copies of wild-type E2F1 were present to drive the cell cycle. This theory is supported by two studies: Cress et al.  created an E2F1-E132 mutant that is artificially mutated in position 132 within E2F1's DNA binding domain and that was demonstrated to have lost its DNA binding capacity, like our R166H mutant; Halaban et al.  demonstrated that expression of the E2F1-E132 mutant can induce a partially transformed phenotype by conferring growth factor-independent cell cycle progression in mice melanocytes. One possible reason why the proliferation of cells over-expressing the E2F1 mutant was not greater than that of control cells is that both mesothelial cell lines used in this study already have a homozygous deletion of the CDKN2A gene resulting in p16 null cells. A key part of the G1/S checkpoint of the cell cycle is p16 deactivation of cyclin-dependent kinase 6, which keeps Rb hypophosphoyylated, thus keeping E2F1 sequestered . A p16 null cell has already lost its G1/S checkpoint control; thus, introducing another mutation that will cause the same checkpoint loss will not cause noticeable growth differences.
Given that WDPMP is a rare sub-type of mesothelioma, it is of interest to extrapolate E2F1's role to the more prevalent MPM. Given that CDKN2A homozygous deletion is prevalent in MPM, with up to 72% of tumors affected , the G1/S checkpoint is already broken in CDKN2A deleted tumors; thus, in terms of proliferation it is unlikely that an additional E2F1 R166H mutation will be useful as the mutation will be redundant in this context. On the other hand, E2F1 also plays an important role in the activation of apoptosis pathways , and the R166H mutation, with its abrogated DNA binding, may contribute to the survival of the cancer cell harboring this mutation. It would be worth checking the remaining 28% of MPMs without CDKN2A deletion for possible mutations in E2F1 and other related genes. It is interesting to note that BAP1 (BRCA1 associated protein-1), a nuclear deubiquitinase affecting E2F and Polycomb target genes, was recently shown to be inactivated by somatic mutations in 23% of MPMs , suggesting that the genes within the E2F pathways might play an important role in mesothelioma in general.
We have performed the first exome sequencing of a WDPMP tumor and a matched control sample and a tumor-derived cell line and discovered the first somatic mutation of E2F1, R166H. This mutation is found to be the critical DNA contact point in the protein's DNA binding domain responsible for gene activation and motif recognition. Experiments confirmed that the mutation abrogates DNA binding and renders the mutated protein unable to adequately up-regulate its target genes. Large accumulation of the mutant protein is observed in over-expression studies and this is due to a great increase in protein stability as shown by a cyclohexamide chase assay. Overall, two compatible theories can explain the observed results: first, E2F1-R166H decreases apoptosis and its abrogated cell cycle role is compensated for by other members of its family; and second, heterozygous E2F1-R166H behaves like SV-40 large T antigen, interfering with the tumor suppressive role of Rb and allowing its wild-type counterpart to drive cell division.
Materials and methods
Tumor and blood samples were collected from a 41-year-old Chinese female who was diagnosed with WDPMP after a laparoscopic biopsy of the omental nodules that were found during a computerized tomography scan. The patient underwent cytoreductive surgery and hyperthermic infusion of intraperitoneal chemotherapy. She completed 5 days of early post-operative intraperitoneal chemotherapy whilst hospitalized, and recovered uneventfully without any complications. She was discharged on post-operative day 15 and remains disease-free at 8 months after her surgery. Informed consent for tissue collection was obtained from the patient by SingHealth Tissue Repository (approved reference number 10-MES-197) and this study was approved by the SingHealth Centralised Institutional Review Board (CIRB reference number 2010-282-B).
Cell line establishment
Fresh tumor sections were first minced into a paste using surgical scissors in a sterile petri dish and the minced section was transferred to a 50-ml falcon conical tube along with 10 ml of 0.1% collagenase (C5138; Sigma, St. Louis, MO, USA) and incubated for 1 hour at 37°C. RPMI1640 (40 ml) was then added to the tube and spun for 5 minutes at 500 g after which the supernatant was removed and the process repeated until the pellet had a white color. The pellet was re-suspended with 14 ml of RPMI1640 containing 10% fetal bovine serum and antibiotics and seeded onto a T-75 flask. The flask was incubated for 24 hours at 37°C in a 5% CO2 environment before being checked under the microscope for cell attachment to the flask surface and the cells were passaged every 3 days.
Extraction of DNA from patient samples and cell lines
For sample DNA extraction, approximately 15 to 20 mg of frozen tissue was measured out and the sample pulverized into a fine powder using a mortar and pestle; the powdered sample was then added to a 15-ml falcon tube containing 2 ml Master mix containing 4 μl of Rnase A, 100 μl of QIAGEN (Valencia, CA, USA) protease and 2 μl of Buffer G2 and mixed thoroughly. The mixture was incubated in a 50°C incubator for 24 hours then spun at maximum speed for 25 minutes before the supernatant was extracted.
DNA was then extracted from the supernatant using QIAGEN's Blood and Cell Culture Mini kit according to the manufacturer's instructions. In brief, the supernatant was loaded into the kit-supplied column (Genomic-Tip 20/G) and the flow-through was discarded. The column was then washed and DNA eluted into a falcon tube and isopropanol added to precipitate the DNA. The tube was then spun at maximum speed for 15 minutes before washing twice with 70% ethanol. The ethanol was discarded and the remaining DNA pellet re-suspended in TE buffer.
Exome capture and paired-end sequencing
Sample exomes were captured using Agilent SureSelect Human All Exon Kit v1.01 designed to encompass 37.8 Mb of the human exon coding region. DNA (3 μg) from the WDPMP tumor, matched blood sample and the tumor-derived cell line were sheared, end-repaired and ligated with paired-end adaptors before hybridizing with biotinylated RNA library baits for 24 hours at 65°C. The DNA-bait RNA fragments were captured using streptavidin-coated magnetic beads and the captured fragments were RNA digested, with the remaining DNA fragments PCR amplified to generate the exon captured sequencing library.
A 15 picomolar concentration of the exome library was used in cluster generation in accordance with Illumina's v3 paired-end cluster-generation protocol. The cluster-generated flow cell was then loaded into the GAIIx sequencer to generate the 76 bp of the first read. After first read completion, the paired end module of GAIIx was used to regenerate the clusters within the flow cell for another 76 bp sequencing of the second read. All raw sequencing data generated are available at the NCBI Sequence Read Archive  [SRA:SRP007386].
Sequence mapping and filtering criteria
Illumina paired-end reads were first converted from Illumina quality scores to Sanger quality scores using the converter module of MAQ before paired-end read alignment to the NCBI hg18 build 36.1 reference genome using the short read aligner BWA (Burrows-Wheeler Aligner)  with default options. The aligned output from BWA was processed by SAMtools  in the following manner. The BWA output was first converted into a compressed BAM format before the aligned sequences were sorted according to chromosomal coordinates. The sorted sequences were then subjected to SAMtools' PCR duplicate removal module to discard sequence pairs with identical outer chromosomal coordinates. Because each sample was sequenced in duplicate, the resulting BAM files representing the duplicate lanes were merged into a single BAM file before the quality filtering step. Quality filtering involved selecting sequences that were uniquely aligned with the reference genome, had less than or equal to four mismatches to the reference genome and had a mapping quality score of at least one. The output result of this filter formed the core sequence file for further downstream analysis.
Generation of exome Hilbert plots
Using the core sequence file described above, we first discarded all intronic bases in the following manner. First, conversion was performed on Agilent's SureSelect exon coordinates file from BED format into space-delimited format specifying the chromosomal location of every exon base. SAMtools' pileup command, using the space-delimited exon coordinate file as a parameter, was used to exclusively output only bases belonging to the exome. Since the pileup command was coded to output only bases with non-zero depth to conserve storage, a quick R script was used to insert in the exome bases that are of zero depth into the initial exome pileup output. This final pileup contains every nucleotide of the exome and its associate sequencing depth sorted by chromosomal coordinates. For the visualization of the entire exome, we used the statistical program R, and in particular HilbertVis, a compact graphical representation of linear data packages . Instead of linearly plotting the sequencing depth versus the exome DNA string, a Hilbert plot computationally wraps the DNA string in a fractal manner onto a two-dimensional grid of pre-determined size and represents the coverage depth via a heat map similar to gene expression data. Red and blue color heat mapping is used to demarcate the borders of each chromosome.
Single nucleotide variant discovery
Additional file 5 shows the SNV discovery pipeline. Aligned reads were processed using Genome Analyzer Toolkit . Reads containing microindels were first locally re-aligned to obtain more accurate quality scores and alignments then quality filtered before consensus calling was performed to obtain the raw SNVs. These raw SNVs were subjected to further quality filtering before being compared against dbSNP130 and 1000 Genomes databases where common SNPs present in the exome were discarded; from this pool of remaining SNVs, only non-synonymous variations occurring in exons or splice sites were retained. This pipeline was performed for tumor, normal sample and cell line exomes and only SNVs that had a quality/depth score > 3 and were present in both the tumor and cell line and not in the normal sample were retained; SNVs in this final pool were considered to be candidate somatic mutations.
Sanger sequencing validation
Primers for sequencing validation were designed using Primer3 . Purified PCR products were sequenced in forward and reverse directions using the ABI PRISM BigDye Terminator Cycle Sequencing Ready Reaction kit (version 3) and an ABI PRISM 3730 Genetic Analyzer (Applied Biosystems, Foster City, CA, USA). Chromatograms were analyzed by SeqScape V2.5 and by manual review. The validation PCR primers were (where F and R stand for forward and reverse, respectively): E2F1_F, 5' GCAGCCACAGTGGGTATTACT 3'; E2F1_R, 5' GGGGAGAAGTCACGCTATGA 3'; TRAF7_F, 5' GCCTTGCTCAGTGTCTTTGA 3'; TRAF7_R, 5' CATGTTGTCCATACTCCAGACC 3'; PPFIBP2_F, 5' CCCTCGAGCCATTTGTATTT 3'; PPFIBP2_R, 5' CCACAGCAGAAGCTGAAAGA 3'.
Protein visualization and homology modeling
Protein modeling of the mutated and wild-type DNA binding domain of E2F1 was done using the automated mode of SWISS-MODEL , a web-based fully automated protein structure homology-modeling server. The basic input requirement from the user is the protein sequence of interest or its UniProt AC code (if available). Swiss-PDBviewer  provides an interface allowing users to visualize and manipulate multiple proteins simultaneously. Structures generated by SWISS-MODEL or experimentally determined structures archived at the RCSB Protein Data Bank  can be downloaded in a compact.pdb format that serves as the input source for this viewer.
Mesothelioma cell lines and mutant plasmid generation
The mesothelioma cell lines MSTO-211H and NCI-H28 (ATCC catalogue number CRL2081 and CRL5820, respectively) were cultured in RPMI-1640 supplemented with 10% fetal bovine serum (v/v). Total RNA extracted from the heterozygous E2F1 mutated mesothelioma sample was used for cDNA synthesis using an iScrip cDNA Synthesis Kit (Bio-Rad, Hercules, CA, USA). Full-length wild-type and mutant E2F1 were amplified using iProof DNA polymerase (Bio-Rad) and E2F1 primers. The primer sequences were: E2F1-ORF-F, 5'-AGTTAAGCTTGACCATGGCCTTGGCCGGGG-3'; E2F1-ORF-R, 5'-AGAATTCCAGAAATCCAGGGGGGTGAGGT-3'. The PCR products were subsequently cloned into pcDNA6/myc-His B (Invitrogen, Carlsbad, CA, USA) using HindIII and EcoRI. Plasmids expressing wild-type E2F1 (pcDNA6-E2F1) or mutant E2F1 (pcDNA6-E2F1/R166H) were validated by dideoxy terminator sequencing. pcDNA3-EGFP was constructed as described previously .
ChIP was carried out in MSTO-211H cells transiently transfected with wild-type E2F1 and E2F1-R166H for 48 hours. Transiently transfected cells were cross-linked with 1% formaldehyde. Chromatin solution pre-cleared with protein G sepharose 4 fast flow (GE Healthcare Life Sciences, Piscataway, NJ, USA) was used for immunoprecipitation with anti-Myc tag antibody (ab9132; Abcam, Boston, MA, USA) targeting Myc tag at the carboxyl terminus of E2F1. Co-precipitated chromatin was eluted from complexes and purified by QIAquick PCR Purification Kit (QIAGEN, Valencia, CA, USA). The presence of SIRT1 and APAF1 promoter was analyzed by semi-quantitative PCR using 2 μl from 35 μl of DNA extraction and GoTaq DNA Polymerase (Promega, Madison, WI, USA). Primer sequences used were: Apaf-1 pro-F, 5'-GGAGACCCTAGGACGACAAG-3'; Apaf-1 pro-R, 5'-CAGTGAAGCAACGAGGATGC-3'. Primers specific for the SIRT1 promoter have been described previously . PCR products were resolved on 2% agarose gel containing ethidium bromide.
Quantitative real-time PCR
Total RNA was extracted using TriPure (Roche, Indianapolis, IN, USA). Total RNA (1 μg) was subjected to cDNA synthesis using an iScrip cDNA Synthesis Kit (Bio-Rad). Expression of target genes was examined using specific primers in combination with SsoFast EvaGreen Supermix using a CFX96 Real-Time PCR Detection System (Bio-Rad). Primers used for detecting E2F1 targets were: SIRT1-F, 5'-TGGCAAAGGAGCAGATTAGTAGG-3'; SIRT1-R, 5'-TCATCCTCCATGGGTTCTTCT-3'; Cyclin E1-F, 5'-GGTTAATGGAGGTGTGTGAAGTC-3'; Cyclin E1-R, 5'-CCATCTGTCACATACGCAAACT-3'; APAF1-F, 5'-TGACATTTCTCACGATGCTACC-3'; APAF1-R, 5'-ATTGTCATCTCCCGTTGCCA-3'; GAPDH-F, 5'-GTGGACCTGACCTGCCGTCT-3'; GAPDH-R, 5'-GGAGGAGTGGGTGTCGCTGT-3'. Primers used for determining transfection efficiency were: E2F1-F, 5'-GCTGAAGGTGCAGAAGCGGC-3'; E2F1-R, 5'-TCCTGCAGCTGTCGGAGGTC-3'; EGFP-F, 5'-CTACGGCGTGCAGTGCTTCA-3'; EGFP-R, 5'- CGCCCTCGAACTTCACCTCG-3'.
Relative expression levels of transcripts were normalized with glyceraldehyde 3-phosphate dehydrogenase (GAPDH) expression level.
E2F1 plasmids were transiently transfected into MSTO-211H and NCI-H28 cells through the use of Effectene (QIAGEN) according to the manufacturer's instructions. Briefly, cells were plated at a density of 60% in a six-well plate. The next day, cells were transfected with 0.4 μg pcDNA6-E2F1, pcDNA6-E2F1/R166H or empty vector using Effectene. After a 48-hour transfection period, the cells were harvested for downstream assays. To determine transfection efficiency, 0.1 μg of pcDNA3-EGFP was co-transfected with 0.3 μg of E2F1 plasmids. Cells were collected for RNA and protein extraction after a 48-hour transfection. Expression of EGFP and E2F1 transcripts was assessed by real-time PCR.
Western blot analysis
Cells were lysed in phosphate-buffered saline containing 1% triton-X100 in the presence of protease inhibitor (Roche, Indianapolis, IN, USA). Total protein extracts (20 μg) were separated using 8% SDS-PAGE, transferred to nitrocellulose membranes and probed with antibody specific for E2F1 (KH95; Santa Cruz Biotechnology, Santa Cruz, CA, USA) and β-actin (AC-15; Sigma).
MSTO-211H cells were transfected with 4 μg of wild-type E2F1 or E2F1-R166H in a 99-mm dish. After 24 hours, cells were harvested and split into a six-well plate. After 20 hours, cells were treated with RPMI containing 25 μg/ml cycloheximide (Sigma). Cells were collected at 30 minute time points and lysed in lysis buffer containing 1% triton-X100 and protease inhibitor. The E2F1 level was then determined by western blot.
Transfected cells were seeded into a 96-well plate at a density of 2 × 103 cells after a 48-hour transfection period. Proliferation rates for cells over-expressing wild-type E2F1 and E2F1-R166H were assessed using the colorimetric 3-(4,5-dimethylthiazol-2yl)-5-(3-carboxymethoxyphenyl)-(4-sulfophenyl)-2H-tetrazoluim assay according to the manufacturer's protocol (MTS; Promega). The assay was performed in triplicate and repeated three times independently.
Statistical analyses were performed with PASW Statistics 18.0 (IBM, Endicott, NY, USA). Differences between individual groups were analyzed using ANOVA followed by post hoc analysis. P-values of < 0.05 are considered statistically significant.
activator protein 1
E2F dimerization partner
E2F transcription factor 1
enhanced green fluorescent protein
glyceraldehyde 3-phosphate dehydrogenase
mitogen-activated protein kinase kinase kinase 3
malignant peritoneal mesothelioma
open reading frame
liprin beta 2
retinoblastoma protein 1
single nucleotide polymorphism
single nucleotide variant
TNF receptor-associated factor 7
well differentiated papillary mesothelioma of the peritoneum.
We would like to thank Ler Lian Dee and Pan You Fu for their help in optimizing the ChIP study. This work is supported in part by funding from Lee Foundation, National Cancer Centre Foundation, NUS Graduate School for Integrative Sciences and Engineering, Cancer Science Institute of Singapore, Duke-NUS core grants, Singapore Ministry of Health and the Agency for Science Technology and Research. The funding agencies played no role in study design, the collection, analysis and interpretation of data, writing of the manuscript, or the decision to submit the manuscript for publication.
- Hoekstra A, Riben M, Frumovitz M, Liu J, Ramirez P: Well differentiated papillary mesothelioma of the peritoneum: a pathological analysis and review of the literature. Gynecol Oncol. 2005, 98: 161-167. 10.1016/j.ygyno.2005.03.031.PubMedView ArticleGoogle Scholar
- Bani-Hani K, Gharaibeh K: Malignant peritoneal mesothelioma. J Surg Oncol. 2005, 91: 17-25. 10.1002/jso.20266.PubMedView ArticleGoogle Scholar
- Clarke JM, Helft P: Long-term survival of a woman with well differentiated papillary mesothelioma of the peritoneum: a case report and review of the literature. J Med Case Reports. 2010, 4: 346-10.1186/1752-1947-4-346.PubMed CentralView ArticleGoogle Scholar
- Jaurand MC: Mechanisms of fiber-induced genotoxicity. Environ Health Perspect. 1997, 105 (Suppl 5): 1073-1084. 10.1289/ehp.97105s51073.PubMedPubMed CentralView ArticleGoogle Scholar
- Pisick E, Salgia R: Molecular biology of malignant mesothelioma: a review. Hematol Oncol Clin N Am. 2005, 19: 997-1023. 10.1016/j.hoc.2005.09.012.View ArticleGoogle Scholar
- Sugarbaker DJ, Richards WG, Gordon GJ, Dong L, De Rienzo A, Maulik G, Glickman JN, Chirieac LR, Hartman ML, Taillon BE, Du L, Bouffard P, Kingsmore SF, Miller NA, Farmer AD, Jensen RV, Gullans SR, Bueno R: Transcriptome sequencing of malignant pleural mesothelioma tumors. Proc Natl Acad Sci USA. 2008, 105: 3521-3526. 10.1073/pnas.0712399105.PubMedPubMed CentralView ArticleGoogle Scholar
- Yang H, Rivera Z, Jube S, Nasu M, Bertino P, Goparaju C, Franzoso G, Lotze MT, Krausz T, Pass HI, Bianchi ME, Carbone M: Programmed necrosis induced by asbestos in human mesothelial cells causes high-mobility group box 1 protein release and resultant inflammation. Proc Natl Acad Sci USA. 2010, 107: 12611-12616. 10.1073/pnas.1006542107.PubMedPubMed CentralView ArticleGoogle Scholar
- Ng SB, Turner EH, Robertson PD, Flygare SD, Bigham AW, Lee C, Shaffer T, Wong M, Bhattacharjee A, Eichler EE, Bamshad M, Nickerson DA, Shendure J: Targeted capture and massively parallel sequencing of 12 human exomes. Nature. 2009, 461: 272-276. 10.1038/nature08250.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Durbin R: Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics. 2009, 25: 1754-1760. 10.1093/bioinformatics/btp324.PubMedPubMed CentralView ArticleGoogle Scholar
- Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, Marth G, Abecasis G, Durbin R, 1000 Genomes Project Data Processing Subgroup: The Sequence alignment/map (SAM) format and SAMtools. Bioinformatics. 2009, 25: 2978-2079.Google Scholar
- Anders S: Visualization of genomic data with the Hilbert curve. Bioinformatics. 2009, 25: 1231-1235. 10.1093/bioinformatics/btp152.PubMedPubMed CentralView ArticleGoogle Scholar
- Musti M, Kettunen E, Dragonieri S, Lindholm P, Cavone D, Serio G, Knuutila S: Cytogenetic and molecular genetic changes in malignant mesothelioma. Cancer Genet Cytogenet. 2006, 170: 9-15. 10.1016/j.cancergencyto.2006.04.011.PubMedView ArticleGoogle Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altschuler D, Gabriel S, Daly M, DePristo MA: The Genome Analysis Toolkit: A MapReduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010, 20: 1297-1303. 10.1101/gr.107524.110.PubMedPubMed CentralView ArticleGoogle Scholar
- Chenna R, Sugawara H, Koike T, Lopez R, Gibson TJ, Higgins DG, Thompson JD: Multiple sequence alignment with the Clustal series of programs. Nucleic Acids Res. 2003, 31: 3497-3500. 10.1093/nar/gkg500.PubMedPubMed CentralView ArticleGoogle Scholar
- Yue P, Melamud E, Moult J: SNPs3D: Candidate gene and SNP selection for association studies. BMC Bioinformatics. 2006, 7: 166-10.1186/1471-2105-7-166.PubMedPubMed CentralView ArticleGoogle Scholar
- Guex N, Peitsch MC: SWISS-MODEL and the Swiss-Pdb Viewer: An environment for comparative protein modelling. Electrophoresis. 1997, 18: 2714-2723. 10.1002/elps.1150181505.PubMedView ArticleGoogle Scholar
- Zheng N, Fraenkel E, Pabo CO, Pavletich NP: Structural basis of DNA recognition by the heterodimeric cell cycle transcription factor E2F-DP. Genes Dev. 1999, 13: 666-674. 10.1101/gad.13.6.666.PubMedPubMed CentralView ArticleGoogle Scholar
- Arnold K, Bordoli L, Jurgen Kopp, Schwede T: The SWISS-MODEL Workspace: a web-based environment for protein structure homology modelling. Bioinformatics. 2006, 22: 195-201. 10.1093/bioinformatics/bti770.PubMedView ArticleGoogle Scholar
- Melo F, Feytmans E: Assessing protein structures with a non-local atomic interaction energy. J Mol Biol. 1998, 277: 1141-1152. 10.1006/jmbi.1998.1665.PubMedView ArticleGoogle Scholar
- Christen M, Hünenberger PH, Bakowies D, Baron R, Bürgi R, Geerke DP, Heinz TN, Kastenholz MA, Kräutler V, Oostenbrink C, Peter C, Trzesniak D, van Gunsteren WF: The GROMOS software for biomolecular simulation: GROMOS05. J Comput Chem. 2005, 26: 1719-1751. 10.1002/jcc.20303.PubMedView ArticleGoogle Scholar
- Bouwmeester T, Bauch A, Ruffner H, Angrand PO, Bergamini G, Croughton K, Cruciat C, Eberhard D, Gagneur J, Ghidelli S, Hopf C, Huhse B, Mangano R, Michon AM, Schirle M, Schlegl J, Schwab M, Stein MA, Bauer A, Casari G, Drewes G, Gavin AC, Jackson DB, Joberty G, Neubauer G, Rick J, Kuster B, Superti-Furga G: A physical and functional map of the human TNF-alpha/NF-kappa B signal transduction pathway. Nat Cell Biol. 2004, 6: 97-105. 10.1038/ncb1086.PubMedView ArticleGoogle Scholar
- Xu LG, Li LY, Shu HB: TRAF7 Potentiates MEKK3-induced AP1 and CHOP Activation and Induces Apoptosis. J Biol Chem. 2004, 274: 17278-17282.View ArticleGoogle Scholar
- Shaulian E, Karin M: AP-1 as a regulator of cell life and death. Nat Cell Biol. 2002, 4: E131-E136. 10.1038/ncb0502-e131.PubMedView ArticleGoogle Scholar
- Serra-Pages C, Medley QG, Tang M, Hart A, Streuli M: Liprins, a family of LAR transmembrane protein-tyrosine phosphatase-interacting proteins. J Biol Chem. 1998, 273: 15611-15620. 10.1074/jbc.273.25.15611.PubMedView ArticleGoogle Scholar
- Colas E, Perez C, Cabrera S, Pedrola N, Monge M, Castellvi J, Eyzaguirre F, Gregorio J, Ruiz A, Llaurado M, Rigau M, Garcia M, Ertekin T, Montes M, Lopez-Lopez R, Carreras R, Xercavins J, Ortega A, Maes T Rosell E, Doll A, Abal M, Reventos J, Gil-Moreno A: Molecular markers of endometrial carcinoma detected in uterine aspirates. Int J Cancer. 2011, 129: 2435-2444. 10.1002/ijc.25901.PubMedView ArticleGoogle Scholar
- Johnson DG, Ohtani K, Nevins JR: Autoregulatory control of E2F1 expression in response to positive and negative regulators of cell cycle progression. Genes Dev. 1994, 8: 1514-1525. 10.1101/gad.8.13.1514.PubMedView ArticleGoogle Scholar
- Stanelle J, Putzer BM: E2F1-induced apoptosis: turning killers into therapeutics. Trends Mol Med. 2006, 12: 177-185. 10.1016/j.molmed.2006.02.002.PubMedView ArticleGoogle Scholar
- Bracken AP, Ciro M, Cocito A, Helin K: E2F target genes: unraveling the biology. Trends Mol Med. 2004, 29: 409-417.Google Scholar
- Field SJ, Tsai FY, Kuo F, Zubiaga AM, Kaelin WG, Livingston DM, Orkin SH, Greenberg ME: E2F-1 Functions in mice to promote apoptosis and suppress proliferation. Cell. 1996, 85: 549-561. 10.1016/S0092-8674(00)81255-6.PubMedView ArticleGoogle Scholar
- Yamasaki L, Jacks T, Bronson R, Goillot E, Harlow E, Dyson NJ: Tumor induction and tissue atrophy in mice lacking E2F-1. Cell. 1996, 85: 537-548. 10.1016/S0092-8674(00)81254-4.PubMedView ArticleGoogle Scholar
- Bracken AP, Ciro M, Cocito A, Helin K: E2F target genes: unraveling the biology. Trends Mol Med. 2004, 29: 409-417.Google Scholar
- Cress WD, Johnson DG, Nevins JR: A genetic analysis of the E2F1 gene distinguishes regulation by Rb, p107 and adenovirus E4. Mol Cell Biol. 1993, 13: 6314-6325.PubMedPubMed CentralView ArticleGoogle Scholar
- Halaban R, Cheng E, Zhang Y, Mandigo CE, Miglarese MR: Release of cell cycle constraints in mouse melanocytes by overexpressed mutant E2F1E132 but not by deletion of p16INK4A or p21WAF1/CIP1. Oncogene. 1998, 16: 2489-2501. 10.1038/sj.onc.1201773.PubMedView ArticleGoogle Scholar
- Nevins JR: The Rb/E2F pathway and cancer. Hum Mol Genet. 2001, 10: 669-703. 10.1093/hmg/10.7.669.View ArticleGoogle Scholar
- Illei PB, Rusch VW, Zakowski MF, Ladanyi M: Homozygous deletion of CDKN2A and codeletion of the methyltheioadenosine phosphorylase gene in the majority of pleural mesotheliomas. Clin Can Res. 2003, 9: 2108-2113.Google Scholar
- Bott M, Brevet M, Taylor BS, Shimizu S, Ito T, Wang L, Creaney J, Lake RA, Zakowski MF, Reva B, Sander C, Delsite R, Powell S, Zhou Q, Shen R, Olshen A, Rusch V, Ladanyi M: The nuclear deubiquitinase BAP1 is commonly inactivated by somatic mutations and 3p21.1 losses in malignant pleural mesothelioma. Nat Genet. 2011, 43: 668-672. 10.1038/ng.855.PubMedPubMed CentralView ArticleGoogle Scholar
- NCBI Sequence Read Archive. [http://www.ncbi.nlm.nih.gov/Traces/sra]
- Rozen S, Skaletsky HJ: Primer3 on the WWW for general users and for biologist programmers. Bioinformatics Methods and Protocols: Methods in Molecular Biology. Edited by: Krawetz S, Misener S. 2000, Totowa, NJ: Humana Press, 365-386.Google Scholar
- RCSB Protein Data Bank. [http://www.pdb.org]
- Ong CK, Ng CY, Leong C, Ng CP, Foo KT, Tan PH, Huynh H: Genomic structure of human OKL38 gene and its differential expression in kidney carcinogenesis. J Biol Chem. 2004, 279: 743-754.PubMedView ArticleGoogle Scholar
- Wang CG, Chen LH, Hou XH, Li Z, Kabra N, Ma Y, Nemoto S, Finkel T, Gu W, Cress WD, Chen J: Interactions between E2F1 and SirT1 regulate apoptotic response to DNA damage. Nat Cell Biol. 2006, 8: 1025-1031. 10.1038/ncb1468.PubMedView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an open access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.