Characterizing the allele- and haplotype-specific copy number landscape of cancer genomes at single-cell resolution with CHISEL

Single-cell barcoding technologies have recently been used to perform whole-genome sequencing of thousands of individual cells in parallel. These technologies provide the opportunity to characterize genomic heterogeneity at single-cell resolution, but their extremely low sequencing coverage (<0.05X per cell) has thus far restricted their use to identification of the total copy number of large multi-megabase segments in individual cells. However, total copy numbers do not distinguish between the two homologous chromosomes in humans, and thus provide a limited view of tumor heterogeneity and evolution missing important events such as copy-neutral loss-of-heterozygosity (LOH). We introduce CHISEL, the first method to infer allele- and haplotype-specific copy numbers in single cells and subpopulations of cells by aggregating sparse signal across thousands of individual cells. We applied CHISEL to 10 single-cell sequencing datasets from 2 breast cancer patients, each dataset containing ≈2000 cells. We identified extensive allele-specific copy-number aberrations (CNAs) in these samples including copy-neutral LOH, whole-genome duplications (WGDs), and mirrored-subclonal CNAs in subpopulations of cells. These allele-specific CNAs alter the copy number of genomic regions containing well-known breast cancer genes including TP53, BRCA2, and PTEN but are invisible to total copy number analysis. We utilized CHISEL’s allele- and haplotype-specific copy numbers to derive a more refined reconstruction of tumor evolution: timing allele-specific CNAs before and after WGDs, identifying low-frequency subclones distinguished by unique CNAs, and uncovering evidence of convergent evolution. This reconstruction is supported by orthogonal analysis of somatic single-nucleotide variants (SNVs) obtained by pooling barcoded reads across clones defined by CHISEL.

Rawcopy: Improved copy number analysis with Affymetrix arrays

Rawcopy is an R package for processing of Affymetrix CytoScan HD, CytoScan 750k and SNP 6.0 microarray raw intensities (CEL files). It uses data from a large number of reference samples to produce log ratio for total copy number analysis and B-allele frequency for allele-specific copy number and heterozygosity analysis. Rawcopy achieves higher signal-to-noise ratio than commonly used free and proprietary alternatives, leading to improved identification of copy number alterations. In addition, Rawcopy visualises each microarray sample for assessment of technical quality, patient identity and genome-wide absolute copy number states.

The influence of dynein processivity control, MAPs, and microtubule ends on directional movement of a localising mRNA

Many cellular constituents travel along microtubules in association with multiple copies of motor proteins. How the activity of these motors is regulated during cargo sorting is poorly understood. In this study, we address this issue using a novel in vitro assay for the motility of localising Drosophila mRNAs bound to native dynein-dynactin complexes. High precision tracking reveals that individual RNPs within a population undergo either diffusive, or highly processive, minus end-directed movements along microtubules. RNA localisation signals stimulate the processive movements, with regulation of dynein-dynactin’s activity rather than its total copy number per RNP, responsible for this effect. Our data support a novel mechanism for multi-motor translocation based on the regulation of dynein processivity by discrete cargo-associated features. Studying the in vitro responses of RNPs to microtubule-associated proteins (MAPs) and microtubule ends provides insights into how an RNA population could navigate the cytoskeletal network and become anchored at its destination in cells.

The degree of segmental aneuploidy measured by total copy number abnormalities to predict survival and recurrence in superficial gastroesophageal adenocarcinoma.

62 Background: Prognostic biomarkers are needed for superficial gastroesophageal adenocarcinoma (EAC) to predict clinical outcomes and select therapy. Although recurrent mutations have been characterized in EAC, little is known about their clinical and prognostic significance. Aneuploidy is predictive of clinical outcome in many malignancies but has not been evaluated in superficial EAC. SNP arrays offer the opportunity to evaluate segmental aneuploidy at high resolution throughout the genome. Methods: We quantified copy number changes in 41 superficial EAC using Affymetrix SNP 6.0 arrays. We identified recurrent chromosomal gains and losses and calculated the total copy number abnormality (CNA) count for each tumor as a measure of aneuploidy. We correlated CNA count with overall survival and time to first recurrence in univariate and multivariate analyses. Results: Recurrent segmental gains and losses involved multiple genes, including: HER2, EGFR, MET, CDK6 , KRAS (recurrent gains); and FHIT, WWOX, CDKN2A/B, SMAD4, RUNX1 (recurrent losses). There was a 40-fold variation in CNA count across all cases. Tumors with the lowest and highest quartile CNA count had significantly better overall survival (p=0.032, log rank test) and time to first recurrence (p=0.010, log rank test) compared to those with intermediate CNA counts. In multivariate Cox analysis, there was a 3.4-fold (95% CI, 1.1–10.4) increased hazard of death among cases with intermediate CNA counts after adjusting for other predictors of survival (N stage, angiolymphatic invasion and tumor size). Similarly, there was a 7.3-fold (95% CI, 1.5-34) increased risk of recurrence for these patients. Conclusions: SNP arrays facilitate the assessment of recurrent chromosomal gain and loss and allow high resolution, quantitative assessment of segmental aneuploidy (total CNA count).The non-monotonic association of segmental aneuploidy with survival has been described in other tumors such as breast and ovarian carcinoma. The degree of segmental aneuploidy is a promising prognostic biomarker in a potentially curable form of EAC.

Genetic characterization of large parathyroid adenomas

In this study, we genetically characterized parathyroid adenomas with large glandular weights, for which independent observations suggest pronounced clinical manifestations. Large parathyroid adenomas (LPTAs) were defined as the 5% largest sporadic parathyroid adenomas identified among the 590 cases operated in our institution during 2005–2009. The LPTA group showed a higher relative number of male cases and significantly higher levels of total plasma and ionized serum calcium (P<0.001). Further analysis of 21 LPTAs revealed low MIB1 proliferation index (0.1–1.5%),MEN1mutations in five cases, and oneHRPT2(CDC73) mutation. Total or partial loss of parafibromin expression was observed in ten tumors, two of which also showed loss of APC expression. Using array CGH, we demonstrated recurrent copy number alterations most frequently involving loss in 1p (29%), gain in 5 (38%), and loss in 11q (33%). Totally, 21 minimal overlapping regions were defined for losses in 1p, 7q, 9p, 11, and 15q and gains in 3q, 5, 7p, 8p, 16q, 17p, and 19q. In addition, 12 tumors showed gross alterations of entire or almost entire chromosomes most frequently gain of 5 and loss of chromosome 11. While gain of 5 was the most frequent alteration observed in LPTAs, it was only detected in a small proportion (4/58 cases, 7%) of parathyroid adenomas. A significant positive correlation was observed between parathyroid hormone level and total copy number gain (r=0.48,P=0.031). These results support that LPTAs represent a group of patients with pronounced parathyroid hyperfunction and associated with specific genomic features.

Expression of a retroposon-like sequence upstream of the putative Trypanosoma brucei variant surface glycoprotein gene expression site promoter

We have cloned the region spanning the putative promoter from two variant surface glycoprotein gene expression sites that are at each end of chromosome M4 of Trypanosoma brucei IsTat 7. Both expression sites contain a retroposon-like sequence (ESR) pseudogene whose 3' end is approximately 30 bp upstream of the putative expression site promoter. The ESRs from both expression sites share considerable sequence homology and are related to LINE-like elements, especially the T. brucei ingi retroposon. Other ESRs are located on large, but not intermediate or mini-, chromosomes in the IsTaR 1 serodeme, and the total copy number is 10 to 20, similar to that estimated for variant surface glycoprotein expression sites. No DNA rearrangements in the vicinity of the ESR and putative expression site promoter were detected following antigenic switches in the IsTaR 1 serodeme. ESR transcripts are present in bloodstream, but not procyclic, forms. Variation in transcript size and sequence between bloodstream variant antigenic types implies that only the ESR from the active expression site is transcribed. This pattern of expression reflects that of sequences downstream of the putative expression site promoter, suggesting that the region of coordinately controlled expression extends upstream of this promoter.