A comprehensive transcriptional portrait of human cancer cell lines.

Tumor-derived cell lines have served as vital models to advance our understanding of oncogene function and therapeutic responses. Although substantial effort has been made to define the genomic constitution of cancer cell line panels, the transcriptome remains understudied. Here we describe RNA sequencing and single-nucleotide polymorphism (SNP) array analysis of 675 human cancer cell lines. We report comprehensive analyses of transcriptome features including gene expression, mutations, gene fusions and expression of non-human sequences. Of the 2,200 gene fusions catalogued, 1,435 consist of genes not previously found in fusions, providing many leads for further investigation. We combine multiple genome and transcriptome features in a pathway-based approach to enhance prediction of response to targeted therapeutics. Our results provide a valuable resource for studies that use cancer cell lines.

Integrative analysis of two cell lines derived from a non-small-lung cancer patient--a panomics approach.

Cancer cells derived from different stages of tumor progression may exhibit distinct biological properties, as exemplified by the paired lung cancer cell lines H1993 and H2073. While H1993 was derived from chemo-naive metastasized tumor, H2073 originated from the chemo-resistant primary tumor from the same patient and exhibits strikingly different drug response profile. To understand the underlying genetic and epigenetic bases for their biological properties, we investigated these cells using a wide range of large-scale methods including whole genome sequencing, RNA sequencing, SNP array, DNA methylation array, and de novo genome assembly. We conducted an integrative analysis of both cell lines to distinguish between potential driver and passenger alterations. Although many genes are mutated in these cell lines, the combination of DNA- and RNA-based variant information strongly implicates a small number of genes including TP53 and STK11 as likely drivers. Likewise, we found a diverse set of genes differentially expressed between these cell lines, but only a fraction can be attributed to changes in DNA copy number or methylation. This set included the ABC transporter ABCC4, implicated in drug resistance, and the metastasis associated MET oncogene. While the rich data content allowed us to reduce the space of hypotheses that could explain most of the observed biological properties, we also caution there is a lack of statistical power and inherent limitations in such single patient case studies.

Comparison of endogenous and overexpressed MyoD shows enhanced binding of physiologically bound sites.

BACKGROUND: Transcription factor overexpression is common in biological experiments and transcription factor amplification is associated with many cancers, yet few studies have directly compared the DNA-binding profiles of endogenous versus overexpressed transcription factors. METHODS: We analyzed MyoD ChIP-seq data from C2C12 mouse myotubes, primary mouse myotubes, and mouse fibroblasts differentiated into muscle cells by overexpression of MyoD and compared the genome-wide binding profiles and binding site characteristics of endogenous and overexpressed MyoD. RESULTS: Overexpressed MyoD bound to the same sites occupied by endogenous MyoD and possessed the same E-box sequence preference and co-factor site enrichments, and did not bind to new sites with distinct characteristics. CONCLUSIONS: Our data demonstrate a robust fidelity of transcription factor binding sites over a range of expression levels and that increased amounts of transcription factor increase the binding at physiologically bound sites.

Genome and transcriptome sequencing of lung cancers reveal diverse mutational and splicing events.

Lung cancer is a highly heterogeneous disease in terms of both underlying genetic lesions and response to therapeutic treatments. We performed deep whole-genome sequencing and transcriptome sequencing on 19 lung cancer cell lines and three lung tumor/normal pairs. Overall, our data show that cell line models exhibit similar mutation spectra to human tumor samples. Smoker and never-smoker cancer samples exhibit distinguishable patterns of mutations. A number of epigenetic regulators, including KDM6A, ASH1L, SMARCA4, and ATAD2, are frequently altered by mutations or copy number changes. A systematic survey of splice-site mutations identified 106 splice site mutations associated with cancer specific aberrant splicing, including mutations in several known cancer-related genes. RAC1b, an isoform of the RAC1 GTPase that includes one additional exon, was found to be preferentially up-regulated in lung cancer. We further show that its expression is significantly associated with sensitivity to a MAP2K (MEK) inhibitor PD-0325901. Taken together, these data present a comprehensive genomic landscape of a large number of lung cancer samples and further demonstrate that cancer-specific alternative splicing is a widespread phenomenon that has potential utility as therapeutic biomarkers. The detailed characterizations of the lung cancer cell lines also provide genomic context to the vast amount of experimental data gathered for these lines over the decades, and represent highly valuable resources for cancer biology.

The effects of hepatitis B virus integration into the genomes of hepatocellular carcinoma patients.

Hepatitis B virus (HBV) infection is a leading risk factor for hepatocellular carcinoma (HCC). HBV integration into the host genome has been reported, but its scale, impact and contribution to HCC development is not clear. Here, we sequenced the tumor and nontumor genomes (>80× coverage) and transcriptomes of four HCC patients and identified 255 HBV integration sites. Increased sequencing to 240× coverage revealed a proportionally higher number of integration sites. Clonal expansion of HBV-integrated hepatocytes was found specifically in tumor samples. We observe a diverse collection of genomic perturbations near viral integration sites, including direct gene disruption, viral promoter-driven human transcription, viral-human transcript fusion, and DNA copy number alteration. Thus, we report the most comprehensive characterization of HBV integration in hepatocellular carcinoma patients. Such widespread random viral integration will likely increase carcinogenic opportunities in HBV-infected individuals.

Differential gene expression in benign prostate epithelium of men with and without prostate cancer: evidence for a prostate cancer field effect.

BACKGROUND: Several malignancies are known to exhibit a "field effect," whereby regions beyond tumor boundaries harbor histologic or molecular changes that are associated with cancer. We sought to determine if histologically benign prostate epithelium collected from men with prostate cancer exhibits features indicative of premalignancy or field effect. EXPERIMENTAL DESIGN: Prostate needle biopsies from 15 men with high-grade (Gleason 8-10) prostate cancer and 15 age- and body mass index-matched controls were identified from a biospecimen repository. Benign epithelia from each patient were isolated by laser capture microdissection. RNA was isolated, amplified, and used for microarray hybridization. Quantitative PCR was used to determine the expression of specific genes of interest. Alterations in protein expression were analyzed through immunohistochemistry. RESULTS: Overall patterns of gene expression in microdissected benign prostate-associated benign epithelium (BABE) and cancer-associated benign epithelium (CABE) were similar. Two genes previously associated with prostate cancer, PSMA and SSTR1, were significantly upregulated in the CABE group (false discovery rate <1%). Expression of other prostate cancer-associated genes, including ERG, HOXC4, HOXC5, and MME, were also increased in CABE by quantitative reverse transcription-PCR, although other genes commonly altered in prostate cancer were not different between the BABE and CABE samples. The expression of MME and PSMA proteins on immunohistochemistry coincided with their mRNA alterations. CONCLUSION: Gene expression profiles between benign epithelia of patients with and without prostate cancer are very similar. However, these tissues exhibit differences in the expression levels of several genes previously associated with prostate cancer development or progression. These differences may comprise a field effect and represent early events in carcinogenesis.

Data standards for flow cytometry.

Flow cytometry (FCM) is an analytical tool widely used for cancer and HIV/AIDS research, and treatment, stem cell manipulation and detecting microorganisms in environmental samples. Current data standards do not capture the full scope of FCM experiments and there is a demand for software tools that can assist in the exploration and analysis of large FCM datasets. We are implementing a standardized approach to capturing, analyzing, and disseminating FCM data that will facilitate both more complex analyses and analysis of datasets that could not previously be efficiently studied. Initial work has focused on developing a community-based guideline for recording and reporting the details of FCM experiments. Open source software tools that implement this standard are being created, with an emphasis on facilitating reproducible and extensible data analyses. As well, tools for electronic collaboration will assist the integrated access and comprehension of experiments to empower users to collaborate on FCM analyses. This coordinated, joint development of bioinformatics standards and software tools for FCM data analysis has the potential to greatly facilitate both basic and clinical research--impacting a notably diverse range of medical and environmental research areas.

Combined cDNA array comparative genomic hybridization and serial analysis of gene expression analysis of breast tumor progression.

To identify genetic changes involved in the progression of breast carcinoma, we did cDNA array comparative genomic hybridization (CGH) on a panel of breast tumors, including 10 ductal carcinoma in situ (DCIS), 18 invasive breast carcinomas, and two lymph node metastases. We identified 49 minimal commonly amplified regions (MCRs) that included known (1q, 8q24, 11q13, 17q21-q23, and 20q13) and several uncharacterized (12p13 and 16p13) regional copy number gains. With the exception of the 17q21 (ERBB2) amplicon, the overall frequency of copy number alterations was higher in invasive tumors than that in DCIS, with several of them present only in invasive cancer. Amplification of candidate loci was confirmed by quantitative PCR in breast carcinomas and cell lines. To identify putative targets of amplicons, we developed a method combining array CGH and serial analysis of gene expression (SAGE) data to correlate copy number and expression levels for each gene within MCRs. Using this approach, we were able to distinguish a few candidate targets from a set of coamplified genes. Analysis of the 12p13-p12 amplicon identified four putative targets: TEL/ETV6, H2AFJ, EPS8, and KRAS2. The amplification of all four candidates was confirmed by quantitative PCR and fluorescence in situ hybridization, but only H2AFJ and EPS8 were overexpressed in breast tumors with 12p13 amplification compared with a panel of normal mammary epithelial cells. These results show the power of combined array CGH and SAGE analysis for the identification of candidate amplicon targets and identify H2AFJ and EPS8 as novel putative oncogenes in breast cancer.