Lingsheng Dong

Transcriptome sequencing of malignant pleural mesothelioma tumors

Cancers arise by the gradual accumulation of mutations in multiple genes. We now use shotgun pyrosequencing to characterize RNA mutations and expression levels unique to malignant pleural mesotheliomas (MPMs) and not present in control tissues. On average, 266 Mb of cDNA were sequenced from each of four MPMs, from a control pulmonary adenocarcinoma (ADCA), and from normal lung tissue. Previously observed differences in MPM RNA expression levels were confirmed. Point mutations were identified by using criteria that require the presence of the mutation in at least four reads and in both cDNA strands and the absence of the mutation from sequence databases, normal adjacent tissues, and other controls. In the four MPMs, 15 nonsynonymous mutations were discovered: 7 were point mutations, 3 were deletions, 4 were exclusively expressed as a consequence of imputed epigenetic silencing, and 1 was putatively expressed as a consequence of RNA editing. Notably, each MPM had a different mutation profile, and no mutated gene was previously implicated in MPM. Of the seven point mutations, three were observed in at least one tumor from 49 other MPM patients. The mutations were in genes that could be causally related to cancer and included XRCC6, PDZK1IP1, ACTR1A, and AVEN.

Four-Gene Expression Ratio Test for Survival in Patients Undergoing Surgery for Mesothelioma

BACKGROUND Malignant pleural mesothelioma has few effective treatments, one being cytoreductive surgery. We previously developed a gene ratio test to predict outcome of malignant pleural mesothelioma patients undergoing surgery. In this study, we investigated the predictive value and technical assay performance of this test in patients with malignant pleural mesothelioma. METHODS Clinical data were obtained prospectively from 120 consecutive patients with malignant pleural mesothelioma who were scheduled for debulking surgery at one institution. Specimens were obtained at surgery or by pleural biopsy examination. Expression data for four genes were collected from tumor specimens, and three ratios of gene expression (TM4SF1/PKM2, TM4SF1/ARHGDIA, and COBLL1/ARHGDIA) were determined by quantitative reverse transcriptase-polymerase chain reaction. Patients were assigned to good or poor outcome groups by the gene ratio test. Survival was estimated by the Kaplan-Meier method and the log-rank test in univariate analyses. A multivariable Cox proportional hazards model was used to control for prognostic factors. Technical robustness was determined by using up to 30 specimens per patient, two biopsy techniques, and two performance sites. All statistical tests were two-sided. RESULTS The test predicted overall survival (P < .001) and cancer-specific survival (P = .007) in univariate analysis and overall survival in multivariable analysis (hazard ratio for death = 2.09, 95% confidence interval [CI] = 1.27 to 3.45, P = .004). The test was reproducible within patients and repeatable between two determinations for specimens with widely varying tumor cell contents. Repeatability between two determinations was 88.5% (95% CI = 84.0% to 92.2%) or, when technically unacceptable test values were excluded, 91.9% (95% CI = 87.4% to 95.1%). Reproducibility between two determinations was 96.1% (95% CI = 86.5% to 99.5%). Combining the gene ratio test and other prognostic factors allowed prospective discrimination between patients at high risk (median survival = 6.9 months, 95% CI = 2.6 to 8.9 months; 3-year survival = 0%) and low risk (median survival = 31.9 months, 95% CI = 21.9 to 41.7 months; 3-year survival = 42%). CONCLUSION The gene ratio test for survival of patients with malignant pleural mesothelioma has robust predictive value and technical assay performance.

Expression patterns of inhibitor of apoptosis proteins in malignant pleural mesothelioma

Inhibitor of apoptosis proteins (IAPs) comprise a family of structurally similar proteins, five of which are widely studied in the context of cancer: IAP‐1/MIHC/cIAP2, IAP‐2/MIHB/cIAP1, livin/ML‐IAP/KIAP, survivin, and XIAP/MIHA/hILP. IAPs are overexpressed by most neoplasms, promote tumour cell survival after a wide variety of apoptotic stimuli, and frequently have gene and/or protein expression patterns associated with a relatively poor prognosis. However, many IAPs are also expressed by normal tissues, can facilitate apoptotic cell death, and have expression patterns associated with a relatively favourable prognosis in some cases. The result is that the precise role(s) of IAPs in human tumours is not exactly known. It has been previously reported that IAP‐1 is overexpressed in malignant pleural mesothelioma (MPM) and is responsible for a large degree of the resistance of cultured MPM cells to cisplatin. Given the high homology of IAP family members, it is likely that other IAPs will be important in MPM. In the present study, the gene and protein expression patterns of IAP‐1, IAP‐2, survivin, livin, and XIAP have been determined in MPM cell lines (n = 9) and a large number of MPM tumours using high‐density oligonucleotide microarrays (n = 40) and an MPM tissue array (n = 66). Human tumours were linked to clinical data and it was found that IAP‐1 and survivin mRNA expression patterns were associated with a relatively shorter patient survival, while those of XIAP and livin were associated with a relatively longer patient survival. Abundant protein for all IAPs was also detected in MPM tumours, where they were expressed primarily in the cytoplasm. Only IAP‐1 and livin protein was expressed in the nucleus of MPM tumours. These results provide the rationale for additional study of this gene family in MPM and cancer in general. Copyright

Inhibitor of apoptosis proteins are regulated by tumour necrosis factor-α in malignant pleural mesothelioma

Inhibitor of apoptosis proteins (IAPs) are overexpressed by most neoplasms and promote tumour cell survival after a wide variety of apoptotic stimuli elicited via intrinsic (ie mitochondrial) and extrinsic (ie death receptor) pathways. It has previously been reported that one of these proteins, IAP‐1(MIHC/cIAP2), is overexpressed in malignant pleural mesothelioma (MPM) and is responsible for a large degree of the resistance of cultured MPM cells to cisplatin. Subsequent analysis in a larger number of human tumours revealed that additional IAPs (eg IAP‐2/MIHB/cIAP1, livin/ML‐IAP/KIAP, survivin, and XIAP/MIHA/hILP) are also overexpressed in MPM and, with the exception of IAP‐2, have expression patterns that correlate with prognosis. In the present study, potential regulatory mechanisms of IAP genes in MPM were investigated and it was found that tumour necrosis factor‐alpha (TNF‐α) can increase mRNA and protein levels of IAP‐1, IAP‐2, and XIAP, but not livin or survivin in MPM cell lines (n = 4). It was also found that IAP gene expression levels are increased concomitantly with translocation to the nucleus of the TNF‐responsive transcription factor NF‐κB. Co‐incubation of MPM cells with TNF‐α and pyrrolidine dithiocarbamate (PDTC), an NF‐κB inhibitor, prevented TNF‐mediated up‐regulation of IAP gene expression levels. In survival studies, TNF‐α was not toxic to MPM cells at any concentration examined. However, MPM cells exposed to TNF‐α were twice as resistant to cisplatin in dose response survival assays compared with unstimulated controls and were found to have a significantly greater fraction of surviving cells at multiple cisplatin concentrations (p < 0.0087). Finally, it was found that levels of circulating TNF‐α were statistically significantly (p = 0.031) (median 312.5 pg/ml) higher in MPM patients (n = 6) prior to surgical tumour debulking compared with those after surgery (median 0 pg/ml). These results when combined with previous observations by our laboratory and others strongly suggest that IAPs act synergistically with TNF family members to promote survival of MPM tumour cells after exposure to cisplatin and possibly other chemotherapeutic drugs. Copyright

Second generation sequencing of the mesothelioma tumor genome.

The current paradigm for elucidating the molecular etiology of cancers relies on the interrogation of small numbers of genes, which limits the scope of investigation. Emerging second-generation massively parallel DNA sequencing technologies have enabled more precise definition of the cancer genome on a global scale. We examined the genome of a human primary malignant pleural mesothelioma (MPM) tumor and matched normal tissue by using a combination of sequencing-by-synthesis and pyrosequencing methodologies to a 9.6X depth of coverage. Read density analysis uncovered significant aneuploidy and numerous rearrangements. Method-dependent informatics rules, which combined the results of different sequencing platforms, were developed to identify and validate candidate mutations of multiple types. Many more tumor-specific rearrangements than point mutations were uncovered at this depth of sequencing, resulting in novel, large-scale, inter- and intra-chromosomal deletions, inversions, and translocations. Nearly all candidate point mutations appeared to be previously unknown SNPs. Thirty tumor-specific fusions/translocations were independently validated with PCR and Sanger sequencing. Of these, 15 represented disrupted gene-encoding regions, including kinases, transcription factors, and growth factors. One large deletion in DPP10 resulted in altered transcription and expression of DPP10 transcripts in a set of 53 additional MPM tumors correlated with survival. Additionally, three point mutations were observed in the coding regions of NKX6-2, a transcription regulator, and NFRKB, a DNA-binding protein involved in modulating NFKB1. Several regions containing genes such as PCBD2 and DHFR, which are involved in growth factor signaling and nucleotide synthesis, respectively, were selectively amplified in the tumor. Second-generation sequencing uncovered all types of mutations in this MPM tumor, with DNA rearrangements representing the dominant type.

Drop-on-Demand Single Cell Isolation and Total RNA Analysis

Technologies that rapidly isolate viable single cells from heterogeneous solutions have significantly contributed to the field of medical genomics. Challenges remain both to enable efficient extraction, isolation and patterning of single cells from heterogeneous solutions as well as to keep them alive during the process due to a limited degree of control over single cell manipulation. Here, we present a microdroplet based method to isolate and pattern single cells from heterogeneous cell suspensions (10% target cell mixture), preserve viability of the extracted cells (97.0±0.8%), and obtain genomic information from isolated cells compared to the non-patterned controls. The cell encapsulation process is both experimentally and theoretically analyzed. Using the isolated cells, we identified 11 stem cell markers among 1000 genes and compare to the controls. This automated platform enabling high-throughput cell manipulation for subsequent genomic analysis employs fewer handling steps compared to existing methods.

Fine-Needle Aspiration Biopsies for Gene Expression Ratio-Based Diagnostic and Prognostic Tests in Malignant Pleural Mesothelioma

Purpose: Malignant pleural mesothelioma (MPM) is an aggressive disease associated with median survival between 9 and 12 months. The correct diagnosis of MPM is sometimes challenging and usually requires solid tissue biopsies rather than fine-needle aspiration biopsies (FNA). We postulated that the accuracy of FNA-based diagnosis might be improved by the addition of molecular tests using a gene expression ratio-based algorithm and that prognostic tests could be similarly performed. Experimental Design: Two MPM and 2 lung cancer cell lines were used to establish the minimal quantity of RNA required to perform the gene ratio test. On the basis of these results, 276 ex vivo FNA biopsies from 63 MPM patients and 250 ex vivo FNA samples from 92 lung cancer patients were analyzed using previously described diagnostic and prognostic tests based on gene expression ratios. Results: We found that the sensitivity of the diagnostic test for MPM was 100% [95% confidence interval (CI): 95%–100%] and the specificity in primary lung adenocarcinoma was 90% (95% CI: 81%–95%). The FNA-based prognostic classification was concordant among 76% (95% CI: 65%–87%) of patients with the risk assignment in a subset of the matched surgical specimens previously analyzed by the prognostic test. Conclusions: Sufficient RNA can be extracted from most FNA biopsies to perform gene expression molecular tests. In particular, we show that the gene expression ratio algorithms performed well when applied to diagnosis and prognosis in MPM. This study provides support for the development of additional RNA molecular tests that may enhance the utility of FNA in the management of other solid cancers. Clin Cancer Res; 17(2); 310–16. ©2010 AACR.

Differentially Expressed Alternatively Spliced Genes in Malignant Pleural Mesothelioma Identified Using Massively Parallel Transcriptome Sequencing

BackgroundAnalyses of Expressed Sequence Tags (ESTs) databases suggest that most human genes have multiple alternative splice variants. The alternative splicing of pre-mRNA is tightly regulated during development and in different tissue types. Changes in splicing patterns have been described in disease states. Recently, we used whole-transcriptome shotgun pryrosequencing to characterize 4 malignant pleural mesothelioma (MPM) tumors, 1 lung adenocarcinoma and 1 normal lung. We hypothesized that alternative splicing profiles might be detected in the sequencing data for the expressed genes in these samples.MethodsWe developed a software pipeline to map the transcriptome read sequences of the 4 MPM samples and 1 normal lung sample onto known exon junction sequences in the comprehensive AceView database of expressed sequences and to count how many reads map to each junction. 13,274,187 transcriptome reads generated by the Roche/454 sequencing platform for 5 samples were compared with 151,486 exon junctions from the AceView database. The exon junction expression index (EJEI) was calculated for each exon junction in each sample to measure the differential expression of alternative splicing events. Top ten exon junctions with the largest EJEI difference between the 4 mesothelioma and the normal lung sample were then examined for differential expression using Quantitative Real Time PCR (qRT-PCR) in the 5 sequenced samples. Two of the differentially expressed exon junctions (ACTG2.aAug05 and CDK4.aAug05) were further examined with qRT-PCR in additional 18 MPM and 18 normal lung specimens.ResultsWe found 70,953 exon junctions covered by at least one sequence read in at least one of the 5 samples. All 10 identified most differentially expressed exon junctions were validated as present by RT-PCR, and 8 were differentially expressed exactly as predicted by the sequence analysis. The differential expression of the AceView exon junctions for the ACTG2 and CDK4 genes were also observed to be statistically significant in an additional 18 MPM and 18 normal lung samples examined using qRT-PCR. The differential expression of these two junctions was shown to successfully classify these mesothelioma and normal lung specimens with high sensitivity (89% and 78%, respectively).ConclusionWhole-transcriptome shotgun sequencing, combined with a downstream bioinformatics pipeline, provides powerful tools for the identification of differentially expressed exon junctions resulting from alternative splice variants. The alternatively spliced genes discovered in the study could serve as useful diagnostic markers as well as potential therapeutic targets for MPM.

A microaliquoting technique for precise histological annotation and optimization of cell content in frozen tissue specimens.

Knowledge of the exact cell content of frozen tissue samples is of growing importance in genomic research. We developed a microaliquoting technique to measure and optimize the cell composition of frozen tumor specimens for molecular studies. Frozen samples of 31 mesothelioma cases were cut in alternating thin and thick sections. Thin sections were stained and evaluated visually. Thick sections, i.e., microaliquots, were annotated using bordering stained sections. A range of cellular heterogeneity was observed among and within samples. Precise annotation of samples was obtained by integration and compared to conventional single face and “front and back” section estimates of cell content. Front and back estimates were more highly correlated with block annotation by microaliquoting than were single face estimates. Both methods yielded discrepant estimates, however, and for some studies may not adequately account for the heterogeneity of mesothelioma or other malignancies with variable cellular composition. High yield and quality RNA was extracted from precision annotated, tumor-enriched subsamples prepared by combining individual microaliquots with the highest tumor cellularity estimates. Microaliquoting provides accurate cell content annotation and permits genomic analysis of enriched subpopulations of cells without fixation or amplification.

A gene expression ratio-based diagnostic test for bladder cancer

Purpose Bladder cancer is relatively common but early detection techniques such as cystoscopy and cytology are somewhat limited. We developed a broadly applicable, platform-independent and clinically relevant method based on simple ratios of gene expression to diagnose human cancers. In this study, we sought to determine whether this technique could be applied to the diagnosis of bladder cancer. Experimental design We developed a model for the diagnosis of bladder cancer using expression profiling data from 80 normal and tumor bladder tissues to identify statistically significant discriminating genes with reciprocal average expression levels in each tissue type. The expression levels of select genes were used to calculate individual gene pair expression ratios in order to assign diagnosis. The optimal model was examined in two additional published microarray data sets and using quantitative RT-PCR in a cohort of 13 frozen benign bladder urothelium samples and 13 bladder cancer samples from our institution. Results A five-ratio test utilizing six genes proved to be 100% accurate (26 of 26 samples) for distinguishing benign from malignant bladder tissue samples (P < 10−6). Conclusions : We have provided a proof of principle study for the use of gene expression ratios in the diagnosis of bladder cancer. This technique may ultimately prove to be a useful adjunct to cytopathology in screening urine specimens for bladder cancer.