The eSNV-detect: a computational system to identify expressed single nucleotide variants from transcriptome sequencing data.

Rapid development of next generation sequencing technology has enabled the identification of genomic alterations from short sequencing reads. There are a number of software pipelines available for calling single nucleotide variants from genomic DNA but, no comprehensive pipelines to identify, annotate and prioritize expressed SNVs (eSNVs) from non-directional paired-end RNA-Seq data. We have developed the eSNV-Detect, a novel computational system, which utilizes data from multiple aligners to call, even at low read depths, and rank variants from RNA-Seq. Multi-platform comparisons with the eSNV-Detect variant candidates were performed. The method was first applied to RNA-Seq from a lymphoblastoid cell-line, achieving 99.7% precision and 91.0% sensitivity in the expressed SNPs for the matching HumanOmni2.5 BeadChip data. Comparison of RNA-Seq eSNV candidates from 25 ER+ breast tumors from The Cancer Genome Atlas (TCGA) project with whole exome coding data showed 90.6-96.8% precision and 91.6-95.7% sensitivity. Contrasting single-cell mRNA-Seq variants with matching traditional multicellular RNA-Seq data for the MD-MB231 breast cancer cell-line delineated variant heterogeneity among the single-cells. Further, Sanger sequencing validation was performed for an ER+ breast tumor with paired normal adjacent tissue validating 29 out of 31 candidate eSNVs. The source code and user manuals of the eSNV-Detect pipeline for Sun Grid Engine and virtual machine are available at http://bioinformaticstools.mayo.edu/research/esnv-detect/.

DIAPH3 governs the cellular transition to the amoeboid tumour phenotype.

Therapies for most malignancies are generally ineffective once metastasis occurs. While tumour cells migrate through tissues using diverse strategies, the signalling networks controlling such behaviours in human tumours are poorly understood. Here we define a role for the Diaphanous-related formin-3 (DIAPH3) as a non-canonical regulator of metastasis that restrains conversion to amoeboid cell behaviour in multiple cancer types. The DIAPH3 locus is close to RB1, within a narrow consensus region of deletion on chromosome 13q in prostate, breast and hepatocellular carcinomas. DIAPH3 silencing in human carcinoma cells destabilized microtubules and induced defective endocytic trafficking, endosomal accumulation of EGFR, and hyperactivation of EGFR/MEK/ERK signalling. Silencing also evoked amoeboid properties, increased invasion and promoted metastasis in mice. In human tumours, DIAPH3 down-regulation was associated with aggressive or metastatic disease. DIAPH3-silenced cells were sensitive to MEK inhibition, but showed reduced sensitivity to EGFR inhibition. These findings have implications for understanding mechanisms of metastasis, and suggest that identifying patients with chromosomal deletions at DIAPH3 may have prognostic value.

Comparative genomic hybridization to detect variation in the copy number of large DNA segments.

Array comparative genomic hybridization (CGH) is an excellent tool to scan the genome for copy number variations (CNVs) when used conscientiously. This article is intended to provide an understanding of the basic principles of array CGH and the different options available to the user to design their array CGH experiments. Specifically, the six subsections discuss the different array platforms available, test and reference DNA preparation, reference DNA choice, the basics of hybridization, data processing, and our current understanding of CNVs in the human genome.

Genomic instability in breast cancer: pathogenesis and clinical implications.

Breast cancer is a heterogeneous disease, appreciable by molecular markers, gene-expression profiles, and most recently, patterns of genomic alteration. In particular, genomic profiling has revealed three distinct patterns of DNA copy-number alteration: a "simple" type with few gains or losses of whole chromosome arms, an "amplifier" type with focal high-level DNA amplifications, and a "complex" type marked by numerous low-amplitude changes and copy-number transitions. The three patterns are associated with distinct gene-expression subtypes, and preferentially target different loci in the genome (implicating distinct cancer genes). Moreover, the different patterns of alteration imply distinct underlying mechanisms of genomic instability. The amplifier pattern may arise from transient telomere dysfunction, although new data suggest ongoing "amplifier" instability. The complex pattern shows similarity to breast cancers with germline BRCA1 mutation, which also exhibit "basal-like" expression profiles and complex-pattern genomes, implicating a possible defect in BRCA1-associated repair of DNA double-strand breaks. As such, targeting presumptive DNA repair defects represents a promising area of clinical investigation. Future studies should clarify the pathogenesis of breast cancers with amplifier and complex-pattern genomes, and will likely identify new therapeutic opportunities.

Comparative analyses of chromosome alterations in soft-tissue metastases within and across patients with castration-resistant prostate cancer.

Androgen deprivation is the mainstay of therapy for progressive prostate cancer. Despite initial and dramatic tumor inhibition, most men eventually fail therapy and die of metastatic castration-resistant (CR) disease. Here, we characterize the profound degree of genomic alteration found in CR tumors using array comparative genomic hybridization (array CGH), gene expression arrays, and fluorescence in situ hybridization (FISH). Bycluster analysis, we show that the similarity of the genomic profiles from primary and metastatic tumors is driven by the patient. Using data adjusted for this similarity, we identify numerous high-frequency alterations in the CR tumors, such as 8p loss and chromosome 7 and 8q gain. By integrating array CGH and expression array data, we reveal genes whose correlated values suggest they are relevant to prostate cancer biology. We find alterations that are significantly associated with the metastases of specific organ sites, and others with CR tumors versus the tumors of patients with localized prostate cancer not treated with androgen deprivation. Within the high-frequency sites of loss in CR metastases, we find an overrepresentation of genes involved in cellular lipid metabolism, including PTEN. Finally, using FISH, we verify the presence of a gene fusion between TMPRSS2 and ERG suggested by chromosome 21 deletions detected by array CGH. We find the fusion in 54% of our CR tumors, and 81% of the fusion-positive tumors contain cells with multiple copies of the fusion. Our investigation lays the foundation for a better understanding of and possible therapeutic targets for CR disease, the poorly responsive and final stage of prostate cancer.

Genomic alterations indicate tumor origin and varied metastatic potential of disseminated cells from prostate cancer patients.

Disseminated epithelial cells can be isolated from the bone marrow of a far greater fraction of prostate-cancer patients than the fraction of patients who progress to metastatic disease. To provide a better understanding of these cells, we have characterized their genomic alterations. We first present an array comparative genomic hybridization method capable of detecting genomic changes in the small number of disseminated cells (10-20) that can typically be obtained from bone marrow aspirates of prostate-cancer patients. We show multiple regions of copy-number change, including alterations common in prostate cancer, such as 8p loss, 8q gain, and gain encompassing the androgen-receptor gene on Xq, in the disseminated cell pools from 11 metastatic patients. We found fewer and less striking genomic alterations in the 48 pools of disseminated cells from patients with organ-confined disease. However, we identify changes shared by these samples with their corresponding primary tumors and prostate-cancer alterations reported in the literature, evidence that these cells, like those in advanced disease, are disseminated tumor cells (DTC). We also show that DTCs from patients with advanced and localized disease share several abnormalities, including losses containing cell-adhesion genes and alterations reported to associate with progressive disease. These shared alterations might confer the capability to disseminate or establish secondary disease. Overall, the spectrum of genomic deviations is evidence for metastatic capacity in advanced-disease DTCs and for variation in that capacity in DTCs from localized disease. Our analysis lays the foundation for elucidation of the relationship between DTC genomic alterations and progressive prostate cancer.

Array comparative genomic hybridization analysis of genomic alterations in breast cancer subtypes.

In this study, we performed high-resolution array comparative genomic hybridization with an array of 4153 bacterial artificial chromosome clones to assess copy number changes in 44 archival breast cancers. The tumors were flow sorted to exclude non-tumor DNA and increase our ability to detect gene copy number changes. In these tumors, losses were more frequent than gains, and gains in 1q and loss in 16q were the most frequent alterations. We compared gene copy number changes in the tumors based on histologic subtype and estrogen receptor (ER) status, i.e., ER-negative infiltrating ductal carcinoma, ER-positive infiltrating ductal carcinoma, and ER-positive infiltrating lobular carcinoma. We observed a consistent association between loss in regions of 5q and ER-negative infiltrating ductal carcinoma, as well as more frequent loss in 4p16, 8p23, 8p21, 10q25, and 17p11.2 in ER-negative infiltrating ductal carcinoma compared with ER-positive infiltrating ductal carcinoma (adjusted P values < or = 0.05). We also observed high-level amplifications in ER-negative infiltrating ductal carcinoma in regions of 8q24 and 17q12 encompassing the c-myc and c-erbB-2 genes and apparent homozygous deletions in 3p21, 5q33, 8p23, 8p21, 9q34, 16q24, and 19q13. ER-positive infiltrating ductal carcinoma showed a higher frequency of gain in 16p13 and loss in 16q21 than ER-negative infiltrating ductal carcinoma. Correlation analysis highlighted regions of change commonly seen together in ER-negative infiltrating ductal carcinoma. ER-positive infiltrating lobular carcinoma differed from ER-positive infiltrating ductal carcinoma in the frequency of gain in 1q and loss in 11q and showed high-level amplifications in 1q32, 8p23, 11q13, and 11q14. These results indicate that array comparative genomic hybridization can identify significant differences in the genomic alterations between subtypes of breast cancer.