Accessible high-throughput single-cell whole-genome sequencing with paired chromatin accessibility

Single-cell whole-genome sequencing (scWGS) enables the assessment of genome-level molecular differences between individual cells with particular relevance to genetically diverse systems like solid tumors. The application of scWGS was limited due to a dearth of accessible platforms capable of producing high-throughput profiles. We present a technique that leverages nucleosome disruption methodologies with the widely adopted 10× Genomics ATAC-seq workflow to produce scWGS profiles for high-throughput copy-number analysis without new equipment or custom reagents. We further demonstrate the use of commercially available indexed transposase complexes from ScaleBio for sample multiplexing, reducing the per-sample preparation costs. Finally, we demonstrate that sequential indexed tagmentation with an intervening nucleosome disruption step allows for the generation of both ATAC and WGS data from the same cell, producing comparable data to the unimodal assays. By exclusively utilizing accessible commercial reagents, we anticipate that these scWGS and scWGS+ATAC methods can be broadly adopted by the research community.

Spatial PD-L1, immune-cell microenvironment, and genomic copy-number alteration patterns and drivers of invasive-disease transition in prospective oral precancer cohort

BACKGROUND

Studies of the immune landscape led to breakthrough trials of programmed death-1 (PD-1) inhibitors for recurrent/metastatic head and neck squamous cell carcinoma therapy. This study investigated the timing, influence of somatic copy-number alterations (SCNAs), and clinical implications of PD-L1 and immune-cell patterns in oral precancer (OPC).

METHODS

The authors evaluated spatial CD3, CD3/8, and CD68 density (cells/mm2 ) and PD-L1 (membranous expression in cytokeratin-positive intraepithelial neoplastic cells and CD68) patterns by multiplex immunofluorescence in a 188-patient prospective OPC cohort, characterized by clinical, histologic, and SCNA risk factors and protocol-specified primary end point of invasive cancer. The authors used Wilcoxon rank-sum and Fisher exact tests, linear mixed effect models, mediation, and Cox regression and recursive-partitioning analyses.

RESULTS

Epithelial, but not CD68 immune-cell, PD-L1 expression was detected in 28% of OPCs, correlated with immune-cell infiltration, 9p21.3 loss of heterozygosity (LOH), and inferior oral cancer-free survival (OCFS), notably in OPCs with low CD3/8 cell density, dysplasia, and/or 9p21.3 LOH. High CD3/8 cell density in dysplastic lesions predicted better OCFS and eliminated the excess risk associated with prior oral cancer and dysplasia. PD-L1 and CD3/8 patterns revealed inferior OCFS in PD-L1 high intrinsic induction and dysplastic immune-cold subgroups.

CONCLUSION

This report provides spatial insight into the immune landscape and drivers of OPCs, and a publicly available immunogenomic data set for future precancer interrogation. The data suggest that 9p21.3 LOH triggers an immune-hot inflammatory phenotype; whereas increased 9p deletion size encompassing CD274 at 9p24.1 may contribute to CD3/8 and PD-L1 depletion during invasive transition. The inferior OCFS in PD-L1-high, immune-cold OPCs support the development of T-cell recruitment strategies.

Copy-number intratumor heterogeneity increases the risk of relapse in chemotherapy-naive stage II colon cancer

Optimal selection of high-risk patients with stage II colon cancer is crucial to ensure clinical benefit of adjuvant chemotherapy. Here, we investigated the prognostic value of genomic intratumor heterogeneity and aneuploidy for disease recurrence. We combined targeted sequencing, SNP arrays, fluorescence in situ hybridization, and immunohistochemistry on a retrospective cohort of 84 untreated stage II colon cancer patients. We assessed the clonality of copy-number alterations (CNAs) and mutations, CD8⁺ lymphocyte infiltration, and their association with time to recurrence. Prognostic factors were included in machine learning analysis to evaluate their ability to predict individual relapse risk. Tumors from recurrent patients displayed a greater proportion of CNAs compared with non-recurrent (mean 31.3% versus 23%, respectively; p = 0.014). Furthermore, patients with elevated tumor CNA load exhibited a higher risk of recurrence compared with those with low levels [p = 0.038; hazard ratio (HR) 2.46], which was confirmed in an independent cohort (p = 0.004; HR 3.82). Candidate chromosome-specific aberrations frequently observed in recurrent cases included gain of the chromosome arm 13q (p = 0.02; HR 2.67) and loss of heterozygosity at 17q22-q24.3 (p = 0.05; HR 2.69). CNA load positively correlated with intratumor heterogeneity (R = 0.52; p < 0.0001). Consistently, incremental subclonal CNAs were associated with an elevated risk of relapse (p = 0.028; HR 2.20), which we did not observe for subclonal single-nucleotide variants and small insertions and deletions. The clinico-genomic model rated an area under the curve of 0.83, achieving a 10% incremental gain compared with clinicopathological markers (p = 0.047). In conclusion, tumor aneuploidy and copy-number intratumor heterogeneity were predictive of poor outcome and improved discriminative performance in early-stage colon cancer. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

A Molecularly Integrated Grade for Meningioma

Background

Meningiomas are the most common primary intracranial tumor in adults. Clinical care is currently guided by the World Health Organization (WHO) grade assigned to meningiomas, a 3-tiered grading system based on histopathology features, as well as extent of surgical resection. Clinical behavior, however, often fails to conform to the WHO grade. Additional prognostic information is needed to optimize patient management.

Methods

We evaluated whether chromosomal copy-number data improved prediction of time-to-recurrence for patients with meningioma who were treated with surgery, relative to the WHO schema. The models were developed using Cox proportional hazards, random survival forest, and gradient boosting in a discovery cohort of 527 meningioma patients and validated in 2 independent cohorts of 172 meningioma patients characterized by orthogonal genomic platforms.

Results

We developed a 3-tiered grading scheme (Integrated Grades 1-3), which incorporated mitotic count and loss of chromosome 1p, 3p, 4, 6, 10, 14q, 18, 19, or CDKN2A. 32% of meningiomas reclassified to either a lower-risk or higher-risk Integrated Grade compared to their assigned WHO grade. The Integrated Grade more accurately identified meningioma patients at risk for recurrence, relative to the WHO grade, as determined by time-dependent area under the curve, average precision, and the Brier score.

Conclusion

We propose a molecularly integrated grading scheme for meningiomas that significantly improves upon the current WHO grading system in prediction of progression-free survival. This framework can be broadly adopted by clinicians with relative ease using widely available genomic technologies and presents an advance in the care of meningioma patients.

Single-cell analysis of copy-number alterations in serous ovarian cancer reveals substantial heterogeneity in both low- and high-grade tumors

Unusually high aneuploidy is a hallmark of epithelial serous ovarian cancer (SOC). Previous analyses have focused on aneuploidy on average across all tumor cells. With the expansion of single-cell sequencing technologies, however, an analysis of copy number heterogeneity cell-to-cell is now technically feasible. Here, we describe an analysis of single-cell RNA sequencing (scRNA-seq) data to infer arm-level aneuploidy in individual serous ovarian cancer cells. By first clustering high-quality sequenced epithelial versus non-epithelial cells, high-confidence tumor cell populations were identified. InferCNV was used to predict segmented copy-number alterations (CNAs), which were then used to determine arm-level aneuploidy at the single-cell level. Control comparisons of normal cells to normal cells showed zero arm-level aneuploidy, whereas a median of four aneuploid events were detectable in cancer cells. A heterogeneity analysis of high-grade tumor cells compared to low-grade tumor cells showed similar levels of cell-to-cell variation between cancer grades. Metastatic tumors potentially showed selection pressure with reduced cell-to-cell variation compared to cells from primary tumors. Minor cell populations with CNAs similar to metastatic cells were identified within the matched primary tumors. Taken together, these results provide a minimum estimate for single-cell aneuploidy in serous ovarian cancer and demonstrate the utility of single-cell sequencing for CNA analysis.

Quantitative analysis of somatically acquired and constitutive uniparental disomy in gastrointestinal cancers

Somatically acquired uniparental disomies (aUPDs) are frequent events in solid tumors and have been associated with cancer‐related genes. Studies assessing their functional consequences across several cancer types are therefore necessary. Here, we aimed at integrating aUPD profiles with the mutational status of cancer‐related genes in a tumor‐type specific manner. Using TCGA datasets for 1,032 gastrointestinal cancers, including colon (COAD), rectum (READ), stomach (STAD), esophageal adenocarcinoma (EAC) and esophageal squamous cell carcinoma (ESCC), we show a non‐random distribution of aUPD, suggesting the existence of a cancer‐specific landscape of aUPD events. Our analysis indicates that aUPD acts as a “second hit” in Knudson's model in order to achieve biallelic inactivation of tumor suppressor genes. In particular, APC, ARID1A and NOTCH1 were recurrently inactivated by the presence of homozygous mutation as a consequence of aUPD in COAD and READ, STAD and ESCC, respectively. Furthermore, while TP53 showed inactivation caused by aUPD at chromosome arm 17p across all tumor types, copy number losses at this genomic position were also frequent. By experimental and computationally inferring genome ploidy, we demonstrate that an increased number of aUPD events, both affecting the whole chromosome or segments of it, were present in highly aneuploid genomes compared to near‐diploid tumors. Finally, the presence of mosaic UPD was detected at a higher frequency in DNA extracted from peripheral blood lymphocytes of patients with colorectal cancer compared to healthy individuals. In summary, our study defines specific profiles of aUPD in gastrointestinal cancers and provides unequivocal evidence of their relevance in cancer.

What's new?

Somatically acquired uniparental disomies (aUPDs), in which two copies of a chromosome originate from the same parent, have been documented in various human cancers. Here, the authors examined the frequency of aUPDs in different gastrointestinal cancer types. Events involving aUPDs were found to occur at high incidence in gastrointestinal cancers and at increased frequency particularly in highly aneuploid genomes. The data also reveal a nonrandom distribution of aUPDs, with evidence of biallelic inactivation of tumor suppressor genes and activation of oncogenes in a tumor type‐specific manner. The findings suggest that aUPDs are functionally relevant in gastrointestinal malignancies.

Integrated Genomic Analysis of Hürthle Cell Cancer Reveals Oncogenic Drivers, Recurrent Mitochondrial Mutations, and Unique Chromosomal Landscapes

The molecular foundations of Hürthle cell carcinoma (HCC) are poorly understood. Here we describe a comprehensive genomic characterization of 56 primary HCC tumors that span the spectrum of tumor behavior. We elucidate the mutational profile and driver mutations and show that these tumors exhibit a wide range of recurrent mutations. Notably, we report a high number of disruptive mutations to both protein-coding and tRNA-encoding regions of the mitochondrial genome. We reveal unique chromosomal landscapes that involve whole-chromosomal duplications of chromosomes 5 and 7 and widespread loss of heterozygosity arising from haploidization and copy-number-neutral uniparental disomy. We also identify fusion genes and disrupted signaling pathways that may drive disease pathogenesis.

Differences in cancer gene copy number alterations between Epstein-Barr virus-positive and Epstein-Barr virus-negative nasopharyngeal carcinoma

BACKGROUND

Nasopharyngeal carcinoma (NPC) treatment is mainly based on clinical staging. We hypothesize that better understanding of the molecular heterogeneity of NPC can aid in better treatment decisions. Therefore, the purpose of this study was to present our exploration of cancer gene copy-number alterations (CNAs) of Epstein-Barr virus (EBV)-positive and EBV-negative NPC.

METHODS

Multiplex ligation-dependent probe amplification was applied to detect CNAs of 36 cancer genes (n = 103). Correlation between CNAs, clinicopathological features, and survival were examined.

RESULTS

The CNAs occurred significantly more in EBV-negative NPC, with PIK3CA and MCCC1 (P < .001) gain/amplification occurring more frequently. Gain/amplification of cyclin-L1 (CCNL1) and PTK2 (P < .001) predict worse disease-free survival (DFS) in EBV-positive NPC.

CONCLUSION

The EBV-positive and EBV-negative NPC show some similarities in cancer gene CNAs suggesting a common pathogenic route but also important differences possibly indicating divergence in oncogenesis. Copy number gain/amplification of CCNL1 and PTK2 are possibly good predictors of survival in EBV-positive NPC.

Differences in MEK inhibitor efficacy in molecularly characterized low-grade serous ovarian cancer cell lines

Advanced or recurrent low-grade serous ovarian cancers (LGSC) are resistant to conventional systemic treatments. LGSC carry mutations in RAS or RAF, leading to several clinical trials evaluating MEK inhibitors (MEKi). As LGSC cell lines and xenografts have been difficult to establish, little is known about the efficacy and on-target activity of MEKi treatment in this disease. We compared four different MEKi (trametinib, selumetinib, binimetinib and refametinib) in novel LGSC patient-derived cell lines. Molecular characterization of these cells included copy-number variation and hotspot mutational analysis. Proliferation, apoptosis and cell viability assays were used to study drug efficacy. MEKi on-target efficacy was measured using western blotting and isoelectric point focusing for ERK1/2 phosphorylation. Ten LGSC cell lines were derived from 7 patients with advanced/recurrent disease. Copy number variation showed significant heterogeneity among cell lines, however all samples showed deletions in chromosome 9p21.3, and frequent copy number gains in chromosomes 12 and 20. Mutations in KRAS/NRAS were identified in 4 patients (57%) and RAS mutation status was not associated with higher baseline levels of ERK phosphorylation. Different degrees of MEKi sensitivity were observed in the LGSC cell lines. Two cell lines, both with KRAS mutations, were highly sensitive to MEKi. Drug anti-proliferative efficacy correlated with the degree of inhibition of ERK phosphorylation, with trametinib being the most potent agent. Differences in MEKi efficacy were observed in LGSC cell lines. Trametinib showed the greatest anti-proliferative effects. This study serves as a basis for much needed future research on MEKi drug efficacy in LGSC.

Genes and Small RNA Transcripts Exhibit Dosage-Dependent Expression Pattern in Maize Copy-Number Alterations

Copy-number alterations are widespread in animal and plant genomes, but their immediate impact on gene expression is still unclear. In animals, copy-number alterations usually exhibit dosage effects, except for sex chromosomes which tend to be dosage compensated. In plants, genes within small duplications (<100 kb) often exhibit dosage-dependent expression, whereas large duplications (>50 Mb) are more often dosage compensated. However, little or nothing is known about expression in moderately-sized (1-50 Mb) segmental duplications, and about the response of small RNAs to dosage change. Here, we compared maize (Zea mays) plants with two, three, and four doses of a 14.6-Mb segment of chromosome 1 that contains ∼300 genes. Plants containing the duplicated segment exhibit dosage-dependent effects on ear length and flowering time. Transcriptome analyses using GeneChip and RNA-sequencing methods indicate that most expressed genes and unique small RNAs within the duplicated segments exhibit dosage-dependent transcript levels. We conclude that dosage effect is the predominant regulatory response for both genes and unique small RNA transcripts in the segmental dosage series we tested. To our knowledge this is the first analysis of small RNA expression in plant gene dosage variants. Because segmental duplications comprise a significant proportion of eukaryotic genomes, these findings provide important new insight into the regulation of genes and small RNAs in response to dosage changes.

Methylation profiling of choroid plexus tumors reveals 3 clinically distinct subgroups

BACKGROUND

Choroid plexus tumors are intraventricular neoplasms derived from the choroid plexus epithelium. A better knowledge of molecular factors involved in choroid plexus tumor biology may aid in identifying patients at risk for recurrence.

METHODS

Methylation profiles were examined in 29 choroid plexus papillomas (CPPs, WHO grade I), 32 atypical choroid plexus papillomas (aCPPs, WHO grade II), and 31 choroid plexus carcinomas (CPCs, WHO grade III) by Illumina Infinium HumanMethylation450 Bead Chip Array.

RESULTS

Unsupervised hierarchical clustering identified 3 subgroups: methylation cluster 1 (pediatric CPP and aCPP of mainly supratentorial location), methylation cluster 2 (adult CPP and aCPP of mainly infratentorial location), and methylation cluster 3 (pediatric CPP, aCPP, and CPC of supratentorial location). In methylation cluster 3, progression-free survival (PFS) accounted for a mean of 72 months (CI, 55-89 mo), whereas only 1 of 42 tumors of methylation clusters 1 and 2 progressed (P< .001). On stratification of outcome data according to WHO grade, all CPCs clustered within cluster 3 and were associated with shorter overall survival (mean, 105 mo [CI, 81-128 mo]) and PFS (mean, 55 mo [CI, 36-73 mo]). The aCPP of methylation cluster 3 also progressed frequently (mean, 69 mo [CI, 44-93 mo]), whereas no tumor progression was observed in aCPP of methylation clusters 1 and 2 (P< .05). Only 1 of 29 CPPs recurred.

CONCLUSIONS

Methylation profiling of choroid plexus tumors reveals 3 distinct subgroups (ie, pediatric low-risk choroid plexus tumors [cluster 1], adult low-risk choroid plexus tumors [cluster 2], and pediatric high-risk choroid plexus tumors [cluster 3]) and may provide useful prognostic information in addition to histopathology.

Analysis of Copy-Number Alterations in Single Cells Using Microarray-Based Comparative Genomic Hybridization (aCGH)

In this unit, we describe a workflow that enables array comparative genomic hybridization (aCGH) of single cells. The unit first describes the isolation and preparation of single peripheral mononuclear cells from blood (PBMC) to prepare a suitable reference DNA for aCGH experiments. An alternative procedure is described for the preparation of single cells of GM14667 and GM05423 cell lines to use as reference DNA and for sex-mismatched control experiments. A guide is also provided for micromanipulation of single cells. Next, the unit describes whole-genome amplification using adapter-linker PCR (Ampli1 WGA Kit) and an alternative nonlinear WGA method (PicoPLEX WGA Kit) for single-cell amplification. A protocol is also included for reamplification of Ampli1 WGA products, which can be used for aCGH as well. Finally, the use of 4 × 180k oligonucleotide microarrays to perform aCGH with single-cell WGA products is described.

Assessing the significance of chromosomal aberrations in cancer: methodology and application to glioma

Comprehensive knowledge of the genomic alterations that underlie cancer is a critical foundation for diagnostics, prognostics, and targeted therapeutics. Systematic efforts to analyze cancer genomes are underway, but the analysis is hampered by the lack of a statistical framework to distinguish meaningful events from random background aberrations. Here we describe a systematic method, called Genomic Identification of Significant Targets in Cancer (GISTIC), designed for analyzing chromosomal aberrations in cancer. We use it to study chromosomal aberrations in 141 gliomas and compare the results with two prior studies. Traditional methods highlight hundreds of altered regions with little concordance between studies. The new approach reveals a highly concordant picture involving approximately 35 significant events, including 16-18 broad events near chromosome-arm size and 16-21 focal events. Approximately half of these events correspond to known cancer-related genes, only some of which have been previously tied to glioma. We also show that superimposed broad and focal events may have different biological consequences. Specifically, gliomas with broad amplification of chromosome 7 have properties different from those with overlapping focalEGFR amplification: the broad events act in part through effects on MET and its ligand HGF and correlate with MET dependence in vitro. Our results support the feasibility and utility of systematic characterization of the cancer genome.