Search

Abstract Background Somatic copy number alterations (SCNAs) are an important class of genomic alteration in cancer. They are frequently observed in cancer samples, with studies showing that, on average, SCNAs affect 34% of a cancer cell’s genome. Furthermore, SCNAs have been shown to be major drivers of tumour development and have been associated with response to therapy and prognosis. Large-scale cancer genome studies suggest that tumours are driven by somatic copy number alterations (SCNAs) or single-nucleotide variants (SNVs). Despite the frequency of SCNAs and their clinical relevance, the use of genomics assays in the clinic is biased towards targeted gene panels, which identify SNVs but provide limited scope to detect SCNAs throughout the genome. There is a need for a comparably low-cost and simple method for high-resolution SCNA profiling. Results We present conliga, a fully probabilistic method that infers SCNA profiles from a low-cost, simple, and clinically-relevant assay (FAST-SeqS). When applied to 11 high-purity oesophageal adenocarcinoma samples, we obtain good agreement (Spearman’s rank correlation coefficient, r s =0.94) between conliga’s inferred SCNA profiles using FAST-SeqS data (approximately £14 per sample) and those inferred by ASCAT using high-coverage WGS (gold-standard). We find that conliga outperforms CNVkit (r s =0.89), also applied to FAST-SeqS data, and is comparable to QDNAseq (r s =0.96) applied to low-coverage WGS, which is approximately four-fold more expensive, more laborious and less clinically-relevant. By performing an in silico dilution series experiment, we find that conliga is particularly suited to detecting SCNAs in low tumour purity samples. At two million reads per sample, conliga is able to detect SCNAs in all nine samples at 3% tumour purity and as low as 0.5% purity in one sample. Crucially, we show that conliga’s hidden state information can be used to decide when a sample is abnormal or normal, whereas CNVkit and QDNAseq cannot provide this critical information. Conclusions We show that conliga provides high-resolution SCNA profiles using a convenient, low-cost assay. We believe conliga makes FAST-SeqS a more clinically valuable assay as well as a useful research tool, enabling inexpensive and fast copy number profiling of pre-malignant and cancer samples.

Summary: Imaging mass cytometry (IMC) is a powerful technique capable of detecting over 30 markers on a single slide. It has been increasingly used for single-cell-based spatial phenotyping in a wide range of samples. However, it only acquires a rectangle field of view (FOV) with a relatively small size and low image resolution, which hinders downstream analysis. Here, we reported a highly practical dual-modality imaging method that combines high-resolution immunofluorescence (IF) and high-dimensional IMC on the same tissue slide. Our computational pipeline uses the whole-slide image (WSI) of IF as a spatial reference and integrates small-FOV IMC into a WSI of IMC. The high-resolution IF images enable accurate single-cell segmentation to extract robust high-dimensional IMC features for downstream analysis. We applied this method in esophageal adenocarcinoma of different stages, identified the single-cell pathology landscape via reconstruction of WSI IMC images, and demonstrated the advantage of the dual-modality imaging strategy. Motivation: Highly multiplexed tissue imaging allows visualization of the spatially resolved expression of multiple proteins at the single-cell level. Although imaging mass cytometry (IMC) using metal isotope-conjugated antibodies has a significant advantage of low background signal and little autofluorescence or batch effect, it has a low resolution that hampers accurate cell segmentation and results in inaccurate feature extraction. In addition, it is also challenging for the current IMC platform to acquire large cm2-sized regions, which limits its application and efficiency when studying larger clinical samples with non-rectangle shapes. To maximize the research output of IMC, we developed a dual-modality imaging method and proposed a comprehensive computational pipeline that combines IF and IMC. The proposed stitching algorithm in this study relies on the IF WSI as a reference for registering and stitching multiple IMC images obtained from different runs. Importantly, this algorithm is not reliant on the positional accuracy of the IMC imaging stage. Considering the possibility of slide removal and reinsertion during WSI IMC data acquisition, it becomes essential to align the IMC images with the comprehensive IF WSI that covers the entire tissue. This registration process plays a critical role in achieving precise spatial alignment, going beyond sole reliance on IMC coordinates.