AbsCN-seq: a statistical method to estimate tumor purity, ploidy and absolute copy numbers from next-generation sequencing data

A deep-learning model for quantifying circulating tumour DNA from the density distribution of DNA-fragment lengths

The quantification of circulating tumour DNA (ctDNA) in blood enables non-invasive surveillance of cancer progression. Here we show that a deep-learning model can accurately quantify ctDNA from the density distribution of cell-free DNA-fragment lengths. We validated the model, which we named 'Fragle', by using low-pass whole-genome-sequencing data from multiple cancer types and healthy control cohorts. In independent cohorts, Fragle outperformed tumour-naive methods, achieving higher accuracy and lower detection limits. We also show that Fragle is compatible with targeted sequencing data. In plasma samples from patients with colorectal cancer, longitudinal analysis with Fragle revealed strong concordance between ctDNA dynamics and treatment responses. In patients with resected lung cancer, Fragle outperformed a tumour-naive gene panel in the prediction of minimal residual disease for risk stratification. The method's versatility, speed and accuracy for ctDNA quantification suggest that it may have broad clinical utility.

Cancer specific up-regulated lactate genes associated with immunotherapy resistance in a pan-cancer analysis

Background

Although the lactate pathway has been reported to lead to immune escape through the inhibition of effector T cells, the cancer-intrinsic lactate signature has not been identified, and the immunotherapeutic efficacy and potential mechanism of the lactate signature are still unclear.

Methods

We defined a pan-cancer up-lactate score by comparing malignant tissues and normal tissues in the TCGA cohort. The immunotherapeutic efficacy was evaluated in non-small cell lung cancer (NSCLC), metastatic renal cancer (mRCC), bladder cancer (BLCA) and melanoma cohorts. The cancer cell-intrinsic mechanism to immune checkpoint inhibitors (ICIs) resistance was measured using single cell sequencing (scRNA-seq) data. Pathway activation was evaluated in the TCGA cohort and CPTAC cohort with transcriptomics and proteomics. The co-occurrence of up-lactate signature and mTOR signaling was determined by spatial transcriptomics of the tissue samples. Immunotherapy resistance and pathway regulation were validated in the in-house NSCLC cohort.

Results

Patients with the high up-lactate scores had significantly short overall survival (OS) than those with the low up-lactate scores (p < 0.001) across multiple types of cancers. The up-regulated lactate signature exhibited higher expression in the malignant cells compared with stromal cells and immune cells in multiple scRNA-seq datasets. A high up-lactate score was associated with poor OS in NSCLC, mRCC, BLCA and melanoma patients who received anti-PD(L)1 antibody. The up-lactate score was higher in the responders of cancer cells, but not in immune cells and stromal cells compared with the non-responders (p < 0.05). Moreover, up-lactate score was positively correlated with mTOR signaling across multiple cancers. In patients with NSCLC who received anti-PD-1 antibody, higher up-lactate scores were associated with significantly shorter PFS compared to lower up-lactate scores (p < 0.001). Additionally, the up-lactate score was associated with cold tumor, and was positively correlated with mTOR signaling.

Conclusion

Collectively, we defined a pan-cancer up-lactate signature, which is a feature of malignant cells and is associated with ICIs resistance. This reveals a coherent program with prognostic and predictive value that may be therapeutically targeted.

Unravelling the Genetic Landscape of Hemiplegic Migraine: Exploring Innovative Strategies and Emerging Approaches

Migraine is a severe, debilitating neurovascular disorder. Hemiplegic migraine (HM) is a rare and debilitating neurological condition with a strong genetic basis. Sequencing technologies have improved the diagnosis and our understanding of the molecular pathophysiology of HM. Linkage analysis and sequencing studies in HM families have identified pathogenic variants in ion channels and related genes, including CACNA1A, ATP1A2, and SCN1A, that cause HM. However, approximately 75% of HM patients are negative for these mutations, indicating there are other genes involved in disease causation. In this review, we explored our current understanding of the genetics of HM. The evidence presented herein summarises the current knowledge of the genetics of HM, which can be expanded further to explain the remaining heritability of this debilitating condition. Innovative bioinformatics and computational strategies to cover the entire genetic spectrum of HM are also discussed in this review.

T proliferating cells derived autophagy signature associated with prognosis and immunotherapy resistance in a pan-cancer analysis

Despite autophagy modulating tumor immunity in the tumor microenvironment (TME), the immunotherapeutic efficacy and potential mechanism of autophagy signature was not explicit. We manually curated an autophagy gene set and defined a pan-cancer autophagy signature by comparing malignant tissues and normal tissues in The Cancer Genome Atlas (TCGA) cohort. The pan-cancer autophagy signature was derived from T proliferating cells as demonstrated in multiple single-cell RNA sequencing (scRNA-seq) datasets. The pan-cancer autophagy signature could influence the cell-cell interactions in the TME and predict the responsiveness of immune checkpoint inhibitors (ICIs) in the metastatic renal cell carcinoma, non-small cell lung cancer, bladder cancer, and melanoma cohorts. Metabolism inactivation accompanied with dysregulation of autophagy was investigated with transcriptomic and proteomic data. The immunotherapeutic predictive role and mechanism regulation of the autophagy signature was validated in an in-house cohort. Our study provides valuable insights into the mechanisms of ICI resistance.

Classification of HER2-negative breast cancers by ERBB2 copy number alteration status reveals molecular differences associated with chromosome 17 gene aberrations

Background

Recently, HER2-negative breast cancers have been reclassified by protein expression into 'HER2-low' and 'HER2-zero' subgroups, but the consideration of HER2-low breast cancer as a distinct biological subtype with differing prognoses remains controversial. By contrast, non-neutral ERBB2 copy number alteration (CNA) status is associated with inferior survival outcomes compared to ERBB2 CNA-neutral breast cancer, providing an alternative approach to classification.

Methods

Here, we investigated the molecular landscape of non-metastatic HER2-negative BCs in relation to ERBB2 CNA status to elucidate biological differences. Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) and The Cancer Genome Atlas (TCGA) TCGA-BRCA datasets (n = 1875) were analyzed.

Results

Nearly two-fifths of the cohort harbored ERBB2 CNAs (39.4%), which were significantly enriched within hormone receptor-negative (56.1%) than within hormone receptor-positive BCs (35.5%; p < 0.0001). Globally, CNAs across the genome were significantly higher in ERBB2 non-neutral compared to neutral cohorts (p < 0.0001). Notably, genetic aberrations on chromosome 17 - BRCA1, NF1, TP53, MAP2K4, and NCOR1 - were widespread in the ERBB2 non-neutral cases. While chromosome 17q arm-level alterations were largely in tandem with ERBB2 CNA status, arm-level loss in chromosome 17p was prevalent regardless of ERBB2 gain, amplification, or loss. Differential gene expression analysis demonstrated that pathways involved in the cell cycle, proteasome, and DNA replication were upregulated in ERBB2 non-neutral cases.

Conclusion

Classification of HER2-negative BCs according to ERBB2 CNA status reveals differences in the genomic landscape. The implications of concurrent aberrations in other genes on chromosome 17 merit further research in ERBB2 non-neutral BCs.

PUREE: accurate pan-cancer tumor purity estimation from gene expression data

Tumors are complex masses composed of malignant and non-malignant cells. Variation in tumor purity (proportion of cancer cells in a sample) can both confound integrative analysis and enable studies of tumor heterogeneity. Here we developed PUREE, which uses a weakly supervised learning approach to infer tumor purity from a tumor gene expression profile. PUREE was trained on gene expression data and genomic consensus purity estimates from 7864 solid tumor samples. PUREE predicted purity with high accuracy across distinct solid tumor types and generalized to tumor samples from unseen tumor types and cohorts. Gene features of PUREE were further validated using single-cell RNA-seq data from distinct tumor types. In a comprehensive benchmark, PUREE outperformed existing transcriptome-based purity estimation approaches. Overall, PUREE is a highly accurate and versatile method for estimating tumor purity and interrogating tumor heterogeneity from bulk tumor gene expression data, which can complement genomics-based approaches or be used in settings where genomic data is unavailable.

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost, a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.

clevRvis: visualization techniques for clonal evolution

BACKGROUND

A thorough analysis of clonal evolution commonly requires integration of diverse sources of data (e.g., karyotyping, next-generation sequencing, and clinical information). Subsequent to actual reconstruction of clonal evolution, detailed analysis and interpretation of the results are essential. Often, however, only few tumor samples per patient are available. Thus, information on clonal development and therapy effect may be incomplete. Furthermore, analysis of biallelic events-considered of high relevance with respect to disease course-can commonly only be realized by time-consuming analysis of the raw results and even raw sequencing data.

RESULTS

We developed clevRvis, an R/Bioconductor package providing an extensive set of visualization techniques for clonal evolution. In addition to common approaches for visualization, clevRvis offers a unique option for allele-aware representation: plaice plots. Biallelic events may be visualized and inspected at a glance. Analyzing 4 public datasets, we show that plaice plots help to gain new insights into tumor development and investigate hypotheses on disease progression and therapy resistance. In addition to a graphical user interface, automatic phylogeny-aware color coding of the plots, and an approach to explore alternative trees, clevRvis provides 2 algorithms for fully automatic time point interpolation and therapy effect estimation. Analyzing 2 public datasets, we show that both approaches allow for valid approximation of a tumor's development in between measured time points.

CONCLUSIONS

clevRvis represents a novel option for user-friendly analysis of clonal evolution, contributing to gaining new insights into tumor development.

Single-cell and bulk transcriptome sequencing identifies two epithelial tumor cell states and refines the consensus molecular classification of colorectal cancer

The consensus molecular subtype (CMS) classification of colorectal cancer is based on bulk transcriptomics. The underlying epithelial cell diversity remains unclear. We analyzed 373,058 single-cell transcriptomes from 63 patients, focusing on 49,155 epithelial cells. We identified a pervasive genetic and transcriptomic dichotomy of malignant cells, based on distinct gene expression, DNA copy number and gene regulatory network. We recapitulated these subtypes in bulk transcriptomes from 3,614 patients. The two intrinsic subtypes, iCMS2 and iCMS3, refine CMS. iCMS3 comprises microsatellite unstable (MSI-H) cancers and one-third of microsatellite-stable (MSS) tumors. iCMS3 MSS cancers are transcriptomically more similar to MSI-H cancers than to other MSS cancers. CMS4 cancers had either iCMS2 or iCMS3 epithelium; the latter had the worst prognosis. We defined the intrinsic epithelial axis of colorectal cancer and propose a refined ‘IMF’ classification with five subtypes, combining intrinsic epithelial subtype (I), microsatellite instability status (M) and fibrosis (F).

A single-cell transcriptomic analysis of 63 patients with colorectal cancer classifies tumor cells into two epithelial subtypes. An improved tumor classification based on epithelial subtype, microsatellite stability and fibrosis reveals differences in pathway activation and metastasis.

Computational analysis of cancer genome sequencing data

Distilling biologically meaningful information from cancer genome sequencing data requires comprehensive identification of somatic alterations using rigorous computational methods. As the amount and complexity of sequencing data have increased, so has the number of tools for analysing them. Here, we describe the main steps involved in the bioinformatic analysis of cancer genomes, review key algorithmic developments and highlight popular tools and emerging technologies. These tools include those that identify point mutations, copy number alterations, structural variations and mutational signatures in cancer genomes. We also discuss issues in experimental design, the strengths and limitations of sequencing modalities and methodological challenges for the future.

Copy number alteration signatures as biomarkers in cancer: a review

Within certain cancers, extensive copy number alterations (CNAs) contribute to a complex and heterogenic genomic profile. This makes it difficult to understand and unravel the distinct molecular dynamics shaping the disease while preventing clinically effective patient stratification. CNA signature analysis represents a novel genomic stratification tool for probing this complexity, offering an intricate framework for deriving CNA patterns at the molecular level. This allows the underlying genomic mechanisms of specific cancers to be revealed, leading to the potential identification of therapeutic targets and prognostic associations. This review outlines the molecular and methodological basis of CNA signatures and focuses on recent advances highlighting their clinical utility, limitations and prospective future as novel diagnostic and prognostic cancer biomarkers.

Genomic characteristics of invasive mucinous adenocarcinoma of the lung with multiple pulmonary sites of involvement

Invasive mucinous adenocarcinoma (IMA) of the lung frequently presents with diffuse pneumonic-type features or multifocal lesions, which are regarded as a pattern of intrapulmonary metastases. However, the genomics of multifocal IMAs have not been well studied. We performed whole exome sequencing on samples taken from 2 to 5 regions in seven patients with synchronous multifocal IMAs of the lung (24 regions total). Early initiating driver events, such as KRAS, NKX2-1, TP53, or ARID1A mutations, are clonal mutations and were present in all multifocal IMAs in each patient. The tumor mutational burden of multifocal IMAs was low (mean: 1.13/mega base), but further analyses suggested intra-tumor heterogeneity. The mutational signature analysis found that IMAs were predominantly associated with endogenous mutational process (signature 1), APOBEC activity (signatures 2 and 13), and defective DNA mismatch repair (signature 6), but not related to smoking signature. IMAs synchronously located in the bilateral lower lobes of two patients with background usual interstitial pneumonia had different mutation types, suggesting that they were double primaries. In conclusion, genomic evidence found in this study indicated the clonal intrapulmonary spread of diffuse pneumonic-type or multifocal IMAs, although they can occur in multicentric origins in the background of usual interstitial pneumonia. IMAs exhibited a heterogeneous genomic landscape despite the low somatic mutation burden. Further studies are warranted to determine the clinical significance of the genomic characteristics of IMAs in expanded cohorts.

STIC: Predicting Single Nucleotide Variants and Tumor Purity in Cancer Genome

Single nucleotide variant (SNV) plays an important role in cellular proliferation and tumorigenesis in various types of human cancer. Next-generation sequencing (NGS) has provided high-throughput data at an unprecedented resolution to predict SNVs. Currently, there exist many computational methods for either germline or somatic SNV discovery from NGS data, but very few of them are versatile enough to adapt to any situations. In the absence of matched normal samples, the prediction of somatic SNVs from single-tumor samples becomes considerably challenging, especially when the tumor purity is unknown. Here, we propose a new approach, STIC, to predict somatic SNVs and estimate tumor purity from NGS data without matched normal samples. The main features of STIC include: (1) extracting a set of SNV-relevant features on each site and training the BP neural network algorithm on the features to predict SNVs; (2) creating an iterative process to distinguish somatic SNVs from germline ones by disturbing allele frequency; and (3) establishing a reasonable relationship between tumor purity and allele frequencies of somatic SNVs to accurately estimate the purity. We quantitatively evaluate the performance of STIC on both simulation and real sequencing datasets, the results of which indicate that STIC outperforms competing methods.

The Analysis of Gene Expression Data Incorporating Tumor Purity Information

The tumor microenvironment is composed of tumor cells, stroma cells, immune cells, blood vessels, and other associated non-cancerous cells. Gene expression measurements on tumor samples are an average over cells in the microenvironment. However, research questions often seek answers about tumor cells rather than the surrounding non-tumor tissue. Previous studies have suggested that the tumor purity (TP)-the proportion of tumor cells in a solid tumor sample-has a confounding effect on differential expression (DE) analysis of high vs. low survival groups. We investigate three ways incorporating the TP information in the two statistical methods used for analyzing gene expression data, namely, differential network (DN) analysis and DE analysis. Analysis 1 ignores the TP information completely, Analysis 2 uses a truncated sample by removing the low TP samples, and Analysis 3 uses TP as a covariate in the underlying statistical models. We use three gene expression data sets related to three different cancers from the Cancer Genome Atlas (TCGA) for our investigation. The networks from Analysis 2 have greater amount of differential connectivity in the two networks than that from Analysis 1 in all three cancer datasets. Similarly, Analysis 1 identified more differentially expressed genes than Analysis 2. Results of DN and DE analyses using Analysis 3 were mostly consistent with those of Analysis 1 across three cancers. However, Analysis 3 identified additional cancer-related genes in both DN and DE analyses. Our findings suggest that using TP as a covariate in a linear model is appropriate for DE analysis, but a more robust model is needed for DN analysis. However, because true DN or DE patterns are not known for the empirical datasets, simulated datasets can be used to study the statistical properties of these methods in future studies.

Tissue-specific cell-free DNA degradation quantifies circulating tumor DNA burden

Profiling of circulating tumor DNA (ctDNA) may offer a non-invasive approach to monitor disease progression. Here, we develop a quantitative method, exploiting local tissue-specific cell-free DNA (cfDNA) degradation patterns, that accurately estimates ctDNA burden independent of genomic aberrations. Nucleosome-dependent cfDNA degradation at promoters and first exon-intron junctions is strongly associated with differential transcriptional activity in tumors and blood. A quantitative model, based on just 6 regulatory regions, could accurately predict ctDNA levels in colorectal cancer patients. Strikingly, a model restricted to blood-specific regulatory regions could predict ctDNA levels across both colorectal and breast cancer patients. Using compact targeted sequencing (<25 kb) of predictive regions, we demonstrate how the approach could enable quantitative low-cost tracking of ctDNA dynamics and disease progression.

Circulating tumour DNA (ctDNA) represents a non-invasive option to monitor cancer progression. Here, the authors perform deep sequencing of plasma cell-free DNA, and find that nucleosome-dependent cfDNA degradation at 6 specific regulatory regions is predictive of ctDNA burden.

HCMMCNVs: hierarchical clustering mixture model of copy number variants detection using whole exome sequencing technology

SUMMARY

In this article, we introduce a hierarchical clustering and Gaussian mixture model with expectation-maximization (EM) algorithm for detecting copy number variants (CNVs) using whole exome sequencing (WES) data. The R shiny package 'HCMMCNVs' is also developed for processing user-provided bam files, running CNVs detection algorithm and conducting visualization. Through applying our approach to 325 cancer cell lines in 22 tumor types from Cancer Cell Line Encyclopedia (CCLE), we show that our algorithm is competitive with other existing methods and feasible in using multiple cancer cell lines for CNVs estimation. In addition, by applying our approach to WES data of 120 oral squamous cell carcinoma (OSCC) samples, our algorithm, using the tumor sample only, exhibits more power in detecting CNVs as compared with the methods using both tumors and matched normal counterparts.

AVAILABILITY AND IMPLEMENTATION

HCMMCNVs R shiny software is freely available at github repository https://github.com/lunching/HCMM_CNVs.and Zenodo https://doi.org/10.5281/zenodo.4593371.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Pan-Cancer Analysis of Ligand-Receptor Cross-talk in the Tumor Microenvironment

Signaling between cancer and nonmalignant (stromal) cells in the tumor microenvironment (TME) is a key to tumor progression. Here, we deconvoluted bulk tumor transcriptomes to infer cross-talk between ligands and receptors on cancer and stromal cells in the TME of 20 solid tumor types. This approach recovered known transcriptional hallmarks of cancer and stromal cells and was concordant with single-cell, in situ hybridization and IHC data. Inferred autocrine cancer cell interactions varied between tissues but often converged on Ephrin, BMP, and FGFR-signaling pathways. Analysis of immune checkpoints nominated interactions with high levels of cancer-to-immune cross-talk across distinct tumor types. Strikingly, PD-L1 was found to be highly expressed in stromal rather than cancer cells. Overall, our study presents a new resource for hypothesis generation and exploration of cross-talk in the TME. SIGNIFICANCE: This study provides deconvoluted bulk tumor transcriptomes across multiple cancer types to infer cross-talk in the tumor microenvironment.

Towards multi-omics characterization of tumor heterogeneity: a comprehensive review of statistical and machine learning approaches

The multi-omics molecular characterization of cancer opened a new horizon for our understanding of cancer biology and therapeutic strategies. However, a tumor biopsy comprises diverse types of cells limited not only to cancerous cells but also to tumor microenvironmental cells and adjacent normal cells. This heterogeneity is a major confounding factor that hampers a robust and reproducible bioinformatic analysis for biomarker identification using multi-omics profiles. Besides, the heterogeneity itself has been recognized over the years for its significant prognostic values in some cancer types, thus offering another promising avenue for therapeutic intervention. A number of computational approaches to unravel such heterogeneity from high-throughput molecular profiles of a tumor sample have been proposed, but most of them rely on the data from an individual omics layer. Since the heterogeneity of cells is widely distributed across multi-omics layers, methods based on an individual layer can only partially characterize the heterogeneous admixture of cells. To help facilitate further development of the methodologies that synchronously account for several multi-omics profiles, we wrote a comprehensive review of diverse approaches to characterize tumor heterogeneity based on three different omics layers: genome, epigenome and transcriptome. As a result, this review can be useful for the analysis of multi-omics profiles produced by many large-scale consortia. Contact:sunkim.bioinfo@snu.ac.kr.

PESM: A novel approach of tumor purity estimation based on sample specific methylation sites

Background: Tumor purity is of great significance for the study of tumor genotyping and the prediction of recurrence, which is significantly affected by tumor heterogeneity. Tumor heterogeneity is the basis of drug resistance in various cancer treatments, and DNA methylation plays a core role in the generation of tumor heterogeneity. Almost all types of cancer cells are associated with abnormal DNA methylation in certain regions of the genome. The selection of tumor-related differential methylation sites, which can be used as an indicator of tumor purity, has important implications for purity assessment. At present, the selection of information sites mostly focuses on inter-tumor heterogeneity and ignores the heterogeneity of tumor growth space that is sample specificity. Results: Considering the specificity of tumor samples and the information gain of individual tumor sample relative to the normal samples, we present an approach, PESM, to evaluate the tumor purity through the specificity difference methylation sites of tumor samples. Applied to more than 200 tumor samples of Prostate adenocarcinoma (PRAD) and Kidney renal clear cell carcinoma (KIRC), it shows that the tumor purity estimated by PESM is highly consistent with other existing methods. In addition, PESM performs better than the method that uses the integrated signal of methylation sites to estimate purity. Therefore, different information sites selection methods have an important impact on the estimation of tumor purity, and the selection of sample specific information sites has a certain significance for accurate identification of tumor purity of samples.

Minimum error calibration and normalization for genomic copy number analysis

BACKGROUND

Copy number variations (CNV) are regional deviations from the normal autosomal bi-allelic DNA content. While germline CNVs are a major contributor to genomic syndromes and inherited diseases, the majority of cancers accumulate extensive "somatic" CNV (sCNV or CNA) during the process of oncogenetic transformation and progression. While specific sCNV have closely been associated with tumorigenesis, intriguingly many neoplasias exhibit recurrent sCNV patterns beyond the involvement of a few cancer driver genes. Currently, CNV profiles of tumor samples are generated using genomic micro-arrays or high-throughput DNA sequencing. Regardless of the underlying technology, genomic copy number data is derived from the relative assessment and integration of multiple signals, with the data generation process being prone to contamination from several sources. Estimated copy number values have no absolute or strictly linear correlation to their corresponding DNA levels, and the extent of deviation differs between sample profiles, which poses a great challenge for data integration and comparison in large scale genome analysis.

RESULTS

In this study, we present a novel method named "Minimum Error Calibration and Normalization for Copy Numbers Analysis" (Mecan4CNA). It only requires CNV segmentation files as input, is platform independent, and has a high performance with limited hardware requirements. For a given multi-sample copy number dataset, Mecan4CNA can batch-normalize all samples to the corresponding true copy number levels of the main tumor clones. Experiments of Mecan4CNA on simulated data showed an overall accuracy of 93% and 91% in determining the normal level and single copy alteration (i.e. duplication or loss of one allele), respectively. Comparison of estimated normal levels and single copy alternations with existing methods and karyotyping data on the NCI-60 tumor cell line produced coherent results. To estimate the method's impact on downstream analyses, we performed GISTIC analyses on the original and Mecan4CNA normalized data from the Cancer Genome Atlas (TCGA) where the normalized data showed prominent improvements of both sensitivity and specificity in detecting focal regions.

CONCLUSIONS

Mecan4CNA provides an advanced method for CNA data normalization, especially in meta-analyses involving large profile numbers and heterogeneous source data quality. With its informative output and visualization options, Mecan4CNA also can improve the interpretation of individual CNA profiles. Mecan4CNA is freely available as a Python package and through its code repository on Github.

Ploidy- and Purity-Adjusted Allele-Specific DNA Analysis Using CLONETv2

High‐throughput DNA sequencing technology provides base‐level and statistically rich information about the genomic content of a sample. In the contexts of cancer research and precision oncology, thousands of genomes from paired tumor and matched normal samples are profiled and processed to determine somatic copy‐number changes and single‐nucleotide variations. Higher‐order informative analyses, in the form of allele‐specific copy‐number assessments or subclonality quantification, require reliable estimates of tumor DNA ploidy and tumor cellularity. CLONETv2 provides a complete set of functions to process matched normal and tumor pairs using patient‐specific genotype data, is independent of low‐level tools (e.g., aligner, segmentation algorithm, mutation caller) and offers high‐level functions to compute allele‐specific copy number from segmented data and to identify subclonal population in the input sample. CLONETv2 is applicable to whole‐genome, whole‐exome and targeted sequencing data generated either from tissue or from liquid biopsy samples. © 2019 The Authors.

Detection of differentially methylated CpG sites between tumor samples with uneven tumor purities

MOTIVATION

Inference of differentially methylated (DM) CpG sites between two groups of tumor samples with different geno- or pheno-types is a critical step to uncover the epigenetic mechanism of tumorigenesis, and identify biomarkers for cancer subtyping. However, as a major source of confounding factor, uneven distributions of tumor purity between two groups of tumor samples will lead to biased discovery of DM sites if not properly accounted for.

RESULTS

We here propose InfiniumDM, a generalized least square model to adjust tumor purity effect for differential methylation analysis. Our method is applicable to a variety of experimental designs including with or without normal controls, different sources of normal tissue contaminations. We compared our method with conventional methods including minfi, limma and limma corrected by tumor purity using simulated datasets. Our method shows significantly better performance at different levels of differential methylation thresholds, sample sizes, mean purity deviations and so on. We also applied the proposed method to breast cancer samples from TCGA database to further evaluate its performance. Overall, both simulation and real data analyses demonstrate favorable performance over existing methods serving similar purpose.

AVAILABILITY AND IMPLEMENTATION

InfiniumDM is a part of R package InfiniumPurify, which is freely available from GitHub (https://github.com/Xiaoqizheng/InfiniumPurify).

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Using R and Bioconductor in Clinical Genomics and Transcriptomics

Bioinformatics pipelines are essential in the analysis of genomic and transcriptomic data generated by next-generation sequencing (NGS). Recent guidelines emphasize the need for rigorous validation and assessment of robustness, reproducibility, and quality of NGS analytic pipelines intended for clinical use. Software tools written in the R statistical language and, in particular, the set of tools available in the Bioconductor repository are widely used in research bioinformatics; and these frameworks offer several advantages for use in clinical bioinformatics, including the breath of available tools, modular nature of software packages, ease of installation, enforcement of interoperability, version control, and short learning curve. This review provides an introduction to R and Bioconductor software, its advantages and limitations for clinical bioinformatics, and illustrative examples of tools that can be used in various steps of NGS analysis.

Comparison of in situ and extraction-based methods for the detection of MET amplifications in solid tumors

In EGFR-treatment naive NSCLC patients, high-level MET amplification is detected in approximately 2-3% and is considered as adverse prognostic factor. Currently, clinical trials with two different inhibitors, capmatinib and tepotinib, are under way both defining different inclusion criteria regarding MET amplification from proven amplification only to defining an exact MET copy number. Here, 45 patient samples, including 10 samples without MET amplification, 5 samples showing a low-level MET amplification, 10 samples with an intermediate-level MET amplification, 10 samples having a high-level MET amplification by a MET/CEN7 ratio ≥2.0 and 10 samples showing a high-level MET amplification with GCN ≥6, were evaluated by MET FISH, MET IHC, a ddPCR copy number assay, a NanoString nCounter copy number assay and an amplicon-based parallel sequencing. The MET IHC had the best concordance with MET FISH followed by the NanoString copy number assay, the ddPCR copy number assay and the custom amplicon-based parallel sequencing assays. The concordance was higher in the high-level amplified cohorts than in the low- and intermediate-level amplified cohorts. In summary, currently extraction-based methods cannot replace the MET FISH for the detection of low-level, intermediate-level and high-level MET amplifications, as the number of false negative results is very high. Only for the detection of high-level amplified samples with a gene copy number ≥6 extraction-based methods are a reliable alternative.

Gliosarcoma Is Driven by Alterations in PI3K/Akt, RAS/MAPK Pathways and Characterized by Collagen Gene Expression Signature

Gliosarcoma is a very rare brain tumor reported to be a variant of glioblastoma (GBM), IDH-wildtype. While differences in molecular and histological features between gliosarcoma and GBM were reported, detailed information on the genetic background of this tumor is lacking. We intend to fill in this knowledge gap by the complex analysis of somatic mutations, indels, copy number variations, translocations and gene expression patterns in gliosarcomas. Using next generation sequencing, we determined somatic mutations, copy number variations (CNVs) and translocations in 10 gliosarcomas. Six tumors have been further subjected to RNA sequencing analysis and gene expression patterns have been compared to those of GBMs. We demonstrate that gliosarcoma bears somatic alterations in gene coding for PI3K/Akt (PTEN, PI3K) and RAS/MAPK (NF1, BRAF) signaling pathways that are crucial for tumor growth. Interestingly, the frequency of PTEN alterations in gliosarcomas was much higher than in GBMs. Aberrations of PTEN were the most frequent and occurred in 70% of samples. We identified genes differentially expressed in gliosarcoma compared to GBM (including collagen signature) and confirmed a difference in the protein level by immunohistochemistry. We found several novel translocations (including translocations in the RABGEF1 gene) creating potentially unfavorable combinations. Collected results on genetic alterations and transcriptomic profiles offer new insights into gliosarcoma pathobiology, highlight differences in gliosarcoma and GBM genetic backgrounds and point out to distinct molecular cues for targeted treatment.

Microenvironmental niche divergence shapes BRCA1-dysregulated ovarian cancer morphological plasticity

How tumor microenvironmental forces shape plasticity of cancer cell morphology is poorly understood. Here, we conduct automated histology image and spatial statistical analyses in 514 high grade serous ovarian samples to define cancer morphological diversification within the spatial context of the microenvironment. Tumor spatial zones, where cancer cell nuclei diversify in shape, are mapped in each tumor. Integration of this spatially explicit analysis with omics and clinical data reveals a relationship between morphological diversification and the dysregulation of DNA repair, loss of nuclear integrity, and increased disease mortality. Within the Immunoreactive subtype, spatial analysis further reveals significantly lower lymphocytic infiltration within diversified zones compared with other tumor zones, suggesting that even immune-hot tumors contain cells capable of immune escape. Our findings support a model whereby a subpopulation of morphologically plastic cancer cells with dysregulated DNA repair promotes ovarian cancer progression through positive selection by immune evasion.

Microenvironmental niche divergence shapes BRCA1-dysregulated ovarian cancer morphological plasticity

How tumor microenvironmental forces shape plasticity of cancer cell morphology is poorly understood. Here, we conduct automated histology image and spatial statistical analyses in 514 high grade serous ovarian samples to define cancer morphological diversification within the spatial context of the microenvironment. Tumor spatial zones, where cancer cell nuclei diversify in shape, are mapped in each tumor. Integration of this spatially explicit analysis with omics and clinical data reveals a relationship between morphological diversification and the dysregulation of DNA repair, loss of nuclear integrity, and increased disease mortality. Within the Immunoreactive subtype, spatial analysis further reveals significantly lower lymphocytic infiltration within diversified zones compared with other tumor zones, suggesting that even immune-hot tumors contain cells capable of immune escape. Our findings support a model whereby a subpopulation of morphologically plastic cancer cells with dysregulated DNA repair promotes ovarian cancer progression through positive selection by immune evasion.

Microenvironmental niche divergence shapes BRCA1-dysregulated ovarian cancer morphological plasticity

How tumor microenvironmental forces shape plasticity of cancer cell morphology is poorly understood. Here, we conduct automated histology image and spatial statistical analyses in 514 high grade serous ovarian samples to define cancer morphological diversification within the spatial context of the microenvironment. Tumor spatial zones, where cancer cell nuclei diversify in shape, are mapped in each tumor. Integration of this spatially explicit analysis with omics and clinical data reveals a relationship between morphological diversification and the dysregulation of DNA repair, loss of nuclear integrity, and increased disease mortality. Within the Immunoreactive subtype, spatial analysis further reveals significantly lower lymphocytic infiltration within diversified zones compared with other tumor zones, suggesting that even immune-hot tumors contain cells capable of immune escape. Our findings support a model whereby a subpopulation of morphologically plastic cancer cells with dysregulated DNA repair promotes ovarian cancer progression through positive selection by immune evasion.

Forward and reverse mutations in stages of cancer development

Background

Massive occurrences of interstitial loss of heterozygosity (LOH) likely resulting from gene conversions were found by us in different cancers as a type of single-nucleotide variations (SNVs), comparable in abundance to the commonly investigated gain of heterozygosity (GOH) type of SNVs, raising the question of the relationships between these two opposing types of cancer mutations.

Methods

In the present study, SNVs in 12 tetra sample and 17 trio sample sets from four cancer types along with copy number variations (CNVs) were analyzed by AluScan sequencing, comparing tumor with white blood cells as well as tissues vicinal to the tumor. Four published “nontumor”-tumor metastasis trios and 246 pan-cancer pairs analyzed by whole-genome sequencing (WGS) and 67 trios by whole-exome sequencing (WES) were also examined.

Results

Widespread GOHs enriched with CG-to-TG changes and associated with nearby CNVs and LOHs enriched with TG-to-CG changes were observed. Occurrences of GOH were 1.9-fold higher than LOH in “nontumor” tissues more than 2 cm away from the tumors, and a majority of these GOHs and LOHs were reversed in “paratumor” tissues within 2 cm of the tumors, forming forward-reverse mutation cycles where the revertant LOHs displayed strong lineage effects that pointed to a sequential instead of parallel development from “nontumor” to “paratumor” and onto tumor cells, which was also supported by the relative frequencies of 26 distinct classes of CNVs between these three types of cell populations.

Conclusions

These findings suggest that developing cancer cells undergo sequential changes that enable the “nontumor” cells to acquire a wide range of forward mutations including ones that are essential for oncogenicity, followed by revertant mutations in the “paratumor” cells to avoid growth retardation by excessive mutation load. Such utilization of forward-reverse mutation cycles as an adaptive mechanism was also observed in cultured HeLa cells upon successive replatings. An understanding of forward-reverse mutation cycles in cancer development could provide a genomic basis for improved early diagnosis, staging, and treatment of cancers.

Electronic supplementary material

The online version of this article (10.1186/s40246-018-0170-6) contains supplementary material, which is available to authorized users.

MHC class II restricted neoantigen peptides predicted by clonal mutation analysis in lung adenocarcinoma patients: implications on prognostic immunological biomarker and vaccine design

Background

Mutant peptides presented by MHC (major histocompatibility complex) Class II in cancer are important targets for cancer immunotherapy. Both animal studies and clinical trials in cancer patients showed that CD4 T cells specific to tumor-derived mutant peptides are essential for the efficacy of immune checkpoint blockade therapy by PD1 antibody.

Results

In this study, we analyzed the next generation sequencing data of 147 lung adenocarcinoma patients from The Cancer Genome Atlas and predicted neoantigens presented by MHC Class I and Class II molecules. We found 18,175 expressed clonal somatic mutations, with an average of 124 per patient. The presentation of mutant peptides by an HLA(human leukocyte antigen) Class II molecule, HLA DRB1, were predicted by NetMHCIIpan3.1. 8804 neo-peptides, including 375 strong binders and 8429 weak binders were found. For HLA DRB1*01:01, 54 strong binders and 896 weak binders were found. The most commonly mutated genes with predicted neo-antigens are KRAS, TTN, RYR2, MUC16, TP53, USH2A, ZFHX4, KEAP1, STK11, FAT3, NAV3 and EGFR.

Conclusions

Our results support the feasibility of discovering individualized HLA Class II presented mutant peptides as candidates for immunodiagnosis and immunotherapy of lung adenocarcinoma.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4958-5) contains supplementary material, which is available to authorized users.

Accounting for tumor purity improves cancer subtype classification from DNA methylation data

Motivation

Tumor sample classification has long been an important task in cancer research. Classifying tumors into different subtypes greatly benefits therapeutic development and facilitates application of precision medicine on patients. In practice, solid tumor tissue samples obtained from clinical settings are always mixtures of cancer and normal cells. Thus, the data obtained from these samples are mixed signals. The 'tumor purity', or the percentage of cancer cells in cancer tissue sample, will bias the clustering results if not properly accounted for.

Results

In this article, we developed a model-based clustering method and an R function which uses DNA methylation microarray data to infer tumor subtypes with the consideration of tumor purity. Simulation studies and the analyses of The Cancer Genome Atlas data demonstrate improved results compared with existing methods.

Availability and implementation

InfiniumClust is part of R package InfiniumPurify , which is freely available from CRAN ( https://cran.r-project.org/web/packages/InfiniumPurify/index.html ).

Contact

hao.wu@emory.edu or xqzheng@shnu.edu.cn.

Supplementary information

Supplementary data are available at Bioinformatics online.

ploidyNGS: visually exploring ploidy with Next Generation Sequencing data

Summary

ploidyNGS is a model-free, open source tool to visualize and explore ploidy levels in a newly sequenced genome, exploiting short read data. We tested ploidyNGS using both simulated and real NGS data of the model yeast Saccharomyces cerevisiae. ploidyNGS allows the identification of the ploidy level of a newly sequenced genome in a visual way.

Availability and Implementation

ploidyNGS is available under the GNU General Public License (GPL) at https://github.com/diriano/ploidyNGS. ploidyNGS is implemented in Python and R.

Contact

diriano@gmail.com.

Universal informative CpG sites for inferring tumor purity from DNA methylation microarray data

Tumor purity is an intrinsic property of tumor samples and potentially has severe impact on many types of data analysis. We have previously developed a statistical method, InfiniumPurify, which could infer purity of a tumor sample given its tumor type (available in TCGA) or a set of informative CpG (iDMC) sites. However, in many clinical practices, researchers may focus on a specific type of tumor samples that is not included in TCGA, and samples which are too few to identify reliable iDMCs. This greatly restricts the application of InfiniumPurify in cancer research. In this paper, we proposed an updated version of InfiniumPurify (termed as uiInfiniumPurify) through identifying a universal set of iDMCs (uiDMCs) and redesigning the algorithm to determine hyper- and hypo-methylation status of each uiDMC. Through the application, we estimated tumor purities of 8830 tumor samples from TCGA. Result shows that our estimates are highly consistent with those by other available methods. Consequently, the updated uiInfiniumPurify, can be applied to a single sample (or a few samples) of interest whose tumor type is not included in TCGA. This characteristic will greatly broaden the application of uiInfiniumPurify in cancer research.

Cell subpopulation deconvolution reveals breast cancer heterogeneity based on DNA methylation signature

Tumour heterogeneity describes the coexistence of divergent tumour cell clones within tumours, which is often caused by underlying epigenetic changes. DNA methylation is commonly regarded as a significant regulator that differs across cells and tissues. In this study, we comprehensively reviewed research progress on estimating of tumour heterogeneity. Bioinformatics-based analysis of DNA methylation has revealed the evolutionary relationships between breast cancer cell lines and tissues. Further analysis of the DNA methylation profiles in 33 breast cancer-related cell lines identified cell line-specific methylation patterns. Next, we reviewed the computational methods in inferring clonal evolution of tumours from different perspectives and then proposed a deconvolution strategy for modelling cell subclonal populations dynamics in breast cancer tissues based on DNA methylation. Further analysis of simulated cancer tissues and real cell lines revealed that this approach exhibits satisfactory performance and relative stability in estimating the composition and proportions of cellular subpopulations. The application of this strategy to breast cancer individuals of the Cancer Genome Atlas's identified different cellular subpopulations with distinct molecular phenotypes. Moreover, the current and potential future applications of this deconvolution strategy to clinical breast cancer research are discussed, and emphasis was placed on the DNA methylation-based recognition of intra-tumour heterogeneity. The wide use of these methods for estimating heterogeneity to further clinical cohorts will improve our understanding of neoplastic progression and the design of therapeutic interventions for treating breast cancer and other malignancies.

Quantification of Multiple Tumor Clones Using Gene Array and Sequencing Data

Cancer development is driven by genomic alterations, including copy number aberrations. The detection of copy number aberrations in tumor cells is often complicated by possible contamination of normal stromal cells in tumor samples and intratumor heterogeneity, namely the presence of multiple clones of tumor cells. In order to correctly quantify copy number aberrations, it is critical to successfully de-convolute the complex structure of the genetic information from tumor samples. In this article, we propose a general Bayesian method for estimating copy number aberrations when there are normal cells and potentially more than one tumor clones. Our method provides posterior probabilities for the proportions of tumor clones and normal cells. We incorporate prior information on the distribution of the copy numbers to prioritize biologically more plausible solutions and alleviate possible identifiability issues that have been observed by many researchers. Our model is flexible and can work for both SNP array and next-generation sequencing data. We compare our method to existing ones and illustrate the advantage of our approach in multiple datasets.

Tumor purity and differential methylation in cancer epigenomics

DNA methylation is an epigenetic modification of DNA molecule that plays a vital role in gene expression regulation. It is not only involved in many basic biological processes, but also considered an important factor for tumorigenesis and other human diseases. Study of DNA methylation has been an active field in cancer epigenomics research. With the advances of high-throughput technologies and the accumulation of enormous amount of data, method development for analyzing these data has gained tremendous interests in the fields of computational biology and bioinformatics. In this review, we systematically summarize the recent developments of computational methods and software tools in high-throughput methylation data analysis with focus on two aspects: differential methylation analysis and tumor purity estimation in cancer studies.

The Clinical Significance and Molecular Features of the Spatial Tumor Shapes in Breast Cancers

Each breast cancer has its unique spatial shape, but the clinical importance and the underlying mechanism for the three-dimensional tumor shapes are mostly unknown. We collected the data on the three-dimensional tumor size and tumor volume data of invasive breast cancers from 2,250 patients who underwent surgery between Jan 2000 and Jul 2007. The degree of tumor eccentricity was estimated by using the difference between the spheroid tumor volume and ellipsoid tumor volume (spheroid-ellipsoid discrepancy, SED). In 41 patients, transcriptome and exome sequencing data obtained. Estimation of more accurate tumor burden by calculating ellipsoid tumor volumes did not improve the outcome prediction when compared to the traditional longest diameter measurement. However, the spatial tumor eccentricity, which was measured by SED, showed significant variation between the molecular subtypes of breast cancer. Additionally, the degree of tumor eccentricity was associated with well-known prognostic factors of breast cancer such as tumor size and lymph node metastasis. Transcriptome data from 41 patients showed significant association between MMP13 and spatial tumor shapes. Network analysis and analysis of TCGA gene expression data suggest that MMP13 is regulated by ERBB2 and S100A7A. The present study validates the usefulness of the current tumor size method in determining tumor stages. Furthermore, we show that the tumors with high eccentricity are more likely to have aggressive tumor characteristics. Genes involved in the extracellular matrix remodeling can be candidate regulators of the spatial tumor shapes in breast cancer.

Identifying recurrent mutations in cancer reveals widespread lineage diversity and mutational specificity

Mutational hotspots indicate selective pressure across a population of tumor samples, but their prevalence within and across cancer types is incompletely characterized. An approach to detect significantly mutated residues, rather than methods that identify recurrently mutated genes, may uncover new biologically and therapeutically relevant driver mutations. Here we developed a statistical algorithm to identify recurrently mutated residues in tumour samples. We applied the algorithm to 11,119 human tumors, spanning 41 cancer types, and identified 470 hotspot somatic substitutions in 275 genes. We find that half of all human tumors possess one or more mutational hotspots with widespread lineage-, position-, and mutant allele-specific differences, many of which are likely functional. In total, 243 hotspots were novel and appeared to affect a broad spectrum of molecular function, including hotspots at paralogous residues of Ras-related small GTPases RAC1 and RRAS2. Redefining hotspots at mutant amino acid resolution will help elucidate the allele-specific differences in their function and could have important therapeutic implications.

BubbleTree: an intuitive visualization to elucidate tumoral aneuploidy and clonality using next generation sequencing data

Tumors are characterized by properties of genetic instability, heterogeneity, and significant oligoclonality. Elucidating this intratumoral heterogeneity is challenging but important. In this study, we propose a framework, BubbleTree, to characterize the tumor clonality using next generation sequencing (NGS) data. BubbleTree simultaneously elucidates the complexity of a tumor biopsy, estimating cancerous cell purity, tumor ploidy, allele-specific copy number, and clonality and represents this in an intuitive graph. We further developed a three-step heuristic method to automate the interpretation of the BubbleTree graph, using a divide-and-conquer strategy. In this study, we demonstrated the performance of BubbleTree with comparisons to similar commonly used tools such as THetA2, ABSOLUTE, AbsCN-seq and ASCAT, using both simulated and patient-derived data. BubbleTree outperformed these tools, particularly in identifying tumor subclonal populations and polyploidy. We further demonstrated BubbleTree's utility in tracking clonality changes from patients' primary to metastatic tumor and dating somatic single nucleotide and copy number variants along the tumor clonal evolution. Overall, the BubbleTree graph and corresponding model is a powerful approach to provide a comprehensive spectrum of the heterogeneous tumor karyotype in human tumors. BubbleTree is R-based and freely available to the research community (https://www.bioconductor.org/packages/release/bioc/html/BubbleTree.html).

Global copy number profiling of cancer genomes

In this article, we introduce a robust and efficient strategy for deriving global and allele-specific copy number alternations (CNA) from cancer whole exome sequencing data based on Log R ratios and B-allele frequencies. Applying the approach to the analysis of over 200 skin cancer samples, we demonstrate its utility for discovering distinct CNA events and for deriving ancillary information such as tumor purity.

AVAILABILITY AND IMPLEMENTATION

https://github.com/xfwang/CLOSE CONTACT: xuefeng.wang@stonybrook.edu or michael.krauthammer@yale.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mutational Profiling Can Establish Clonal or Independent Origin in Synchronous Bilateral Breast and Other Tumors

Background

Synchronous tumors can be independent primary tumors or a primary-metastatic (clonal) pair, which may have clinical implications. Mutational profiling of tumor DNA is increasingly common in the clinic. We investigated whether mutational profiling can distinguish independent from clonal tumors in breast and other cancers, using a carefully defined test based on the Clonal Likelihood Score (CLS = 100 x # shared high confidence (HC) mutations/ # total HC mutations).

Methods

Statistical properties of a formal test using the CLS were investigated. A high CLS is evidence in favor of clonality; the test is implemented as a one-sided binomial test of proportions. Test parameters were empirically determined using 16,422 independent breast tumor pairs and 15 primary-metastatic tumor pairs from 10 cancer types using The Cancer Genome Atlas.

Results

We validated performance of the test with its established parameters, using five published data sets comprising 15,758 known independent tumor pairs (maximum CLS = 4.1%, minimum p-value = 0.48) and 283 known tumor clonal pairs (minimum CLS 13%, maximum p-value <0.01), across renal cell, testicular, and colorectal cancer. The CLS test correctly classified all validation samples but one, which it appears may have been incorrectly classified in the published data. As proof-of-concept we then applied the CLS test to two new cases of invasive synchronous bilateral breast cancer at our institution, each with one hormone receptor positive (ER+/PR+/HER2-) lobular and one triple negative ductal carcinoma. High confidence mutations were identified by exome sequencing and results were validated using deep targeted sequencing. The first tumor pair had CLS of 81% (p-value < 10–15), supporting clonality. In the second pair, no common mutations of 184 variants were validated (p-value >0.99), supporting independence. A plausible molecular mechanism for the shift from hormone receptor positive to triple negative was identified in the clonal pair.

Conclusion

We have developed the statistical properties of a carefully defined Clonal Likelihood Score test from mutational profiling of tumor DNA. Under identified conditions, the test appears to reliably distinguish between synchronous tumors of clonal and of independent origin in several cancer types. This approach may have scientific and clinical utility.

BACOM2.0 facilitates absolute normalization and quantification of somatic copy number alterations in heterogeneous tumor

Most published copy number datasets on solid tumors were obtained from specimens comprised of mixed cell populations, for which the varying tumor-stroma proportions are unknown or unreported. The inability to correct for signal mixing represents a major limitation on the use of these datasets for subsequent analyses, such as discerning deletion types or detecting driver aberrations. We describe the BACOM2.0 method with enhanced accuracy and functionality to normalize copy number signals, detect deletion types, estimate tumor purity, quantify true copy numbers, and calculate average-ploidy value. While BACOM has been validated and used with promising results, subsequent BACOM analysis of the TCGA ovarian cancer dataset found that the estimated average tumor purity was lower than expected. In this report, we first show that this lowered estimate of tumor purity is the combined result of imprecise signal normalization and parameter estimation. Then, we describe effective allele-specific absolute normalization and quantification methods that can enhance BACOM applications in many biological contexts while in the presence of various confounders. Finally, we discuss the advantages of BACOM in relation to alternative approaches. Here we detail this revised computational approach, BACOM2.0, and validate its performance in real and simulated datasets.

The promise of omics-based approaches to cancer prevention

Cancer is a complex category of diseases caused in large part by genetic or genomic, transcriptomic, and epigenetic or epigenomic alterations in affected cells and the surrounding microenvironment. Carcinogenesis reflects the clonal expansion of cells that progressively acquire these genetic and epigenetic alterations-changes that, in turn, lead to modifications at the RNA level. Gradually advancing technology and most recently, the advent of next-generation sequencing (NGS), combined with bioinformatics analytic tools, have revolutionized our ability to interrogate cancer cells. The ultimate goal is to apply these high-throughput technologies to the various aspects of clinical cancer care: cancer-risk assessment, diagnosis, as well as target identification for treatment and prevention. In this article, we emphasize how the knowledge gained through large-scale omics-oriented approaches, with a focus on variations at the level of nucleic acids, can inform the field of chemoprevention.

ENVE: a novel computational framework characterizes copy-number mutational landscapes in colorectal cancers from African American patients

Reliable detection of somatic copy-number alterations (sCNAs) in tumors using whole-exome sequencing (WES) remains challenging owing to technical (inherent noise) and sample-associated variability in WES data. We present a novel computational framework, ENVE, which models inherent noise in any WES dataset, enabling robust detection of sCNAs across WES platforms. ENVE achieved high concordance with orthogonal sCNA assessments across two colorectal cancer (CRC) WES datasets, and consistently outperformed a best-in-class algorithm, Control-FREEC. We subsequently used ENVE to characterize global sCNA landscapes in African American CRCs, identifying genomic aberrations potentially associated with CRC pathogenesis in this population. ENVE is downloadable at https://github.com/ENVE-Tools/ENVE.

Predicting tumor purity from methylation microarray data

MOTIVATION

In cancer genomics research, one important problem is that the solid tissue sample obtained from clinical settings is always a mixture of cancer and normal cells. The sample mixture brings complication in data analysis and results in biased findings if not correctly accounted for. Estimating tumor purity is of great interest, and a number of methods have been developed using gene expression, copy number variation or point mutation data.

RESULTS

We discover that in cancer samples, the distributions of data from Illumina Infinium 450 k methylation microarray are highly correlated with tumor purities. We develop a simple but effective method to estimate purities from the microarray data. Analyses of the Cancer Genome Atlas lung cancer data demonstrate favorable performance of the proposed method.

AVAILABILITY AND IMPLEMENTATION

The method is implemented in InfiniumPurify, which is freely available at https://bitbucket.org/zhengxiaoqi/infiniumpurify.

CONTACT

xqzheng@shnu.edu.cn or hao.wu@emory.edu

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Single-cell mRNA sequencing identifies subclonal heterogeneity in anti-cancer drug responses of lung adenocarcinoma cells

Background

Intra-tumoral genetic and functional heterogeneity correlates with cancer clinical prognoses. However, the mechanisms by which intra-tumoral heterogeneity impacts therapeutic outcome remain poorly understood. RNA sequencing (RNA-seq) of single tumor cells can provide comprehensive information about gene expression and single-nucleotide variations in individual tumor cells, which may allow for the translation of heterogeneous tumor cell functional responses into customized anti-cancer treatments.

Results

We isolated 34 patient-derived xenograft (PDX) tumor cells from a lung adenocarcinoma patient tumor xenograft. Individual tumor cells were subjected to single cell RNA-seq for gene expression profiling and expressed mutation profiling. Fifty tumor-specific single-nucleotide variations, including KRAS^G12D, were observed to be heterogeneous in individual PDX cells. Semi-supervised clustering, based on KRAS^G12D mutant expression and a risk score representing expression of 69 lung adenocarcinoma-prognostic genes, classified PDX cells into four groups. PDX cells that survived in vitro anti-cancer drug treatment displayed transcriptome signatures consistent with the group characterized by KRAS^G12D and low risk score.

Conclusions

Single-cell RNA-seq on viable PDX cells identified a candidate tumor cell subgroup associated with anti-cancer drug resistance. Thus, single-cell RNA-seq is a powerful approach for identifying unique tumor cell-specific gene expression profiles which could facilitate the development of optimized clinical anti-cancer strategies.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0692-3) contains supplementary material, which is available to authorized users.