Identification of B cell subsets based on antigen receptor sequences using deep learning.

B cell receptors (BCRs) denote antigen specificity, while corresponding cell subsets indicate B cell functionality. Since each B cell uniquely encodes this combination, physical isolation and subsequent processing of individual B cells become indispensable to identify both attributes. However, this approach accompanies high costs and inevitable information loss, hindering high-throughput investigation of B cell populations. Here, we present BCR-SORT, a deep learning model that predicts cell subsets from their corresponding BCR sequences by leveraging B cell activation and maturation signatures encoded within BCR sequences. Subsequently, BCR-SORT is demonstrated to improve reconstruction of BCR phylogenetic trees, and reproduce results consistent with those verified using physical isolation-based methods or prior knowledge. Notably, when applied to BCR sequences from COVID-19 vaccine recipients, it revealed inter-individual heterogeneity of evolutionary trajectories towards Omicron-binding memory B cells. Overall, BCR-SORT offers great potential to improve our understanding of B cell responses.

Mapping cancer biology in space: applications and perspectives on spatial omics for oncology.

Technologies to decipher cellular biology, such as bulk sequencing technologies and single-cell sequencing technologies, have greatly assisted novel findings in tumor biology. Recent findings in tumor biology suggest that tumors construct architectures that influence the underlying cancerous mechanisms. Increasing research has reported novel techniques to map the tissue in a spatial context or targeted sampling-based characterization and has introduced such technologies to solve oncology regarding tumor heterogeneity, tumor microenvironment, and spatially located biomarkers. In this study, we address spatial technologies that can delineate the omics profile in a spatial context, novel findings discovered via spatial technologies in oncology, and suggest perspectives regarding therapeutic approaches and further technological developments.

Hema-seq reveals genomic aberrations in a rare simultaneous occurrence of hematological malignancies.

Co-occurrence of multiple myeloma and acute myelogenous leukemia is rare, with both malignancies often tracing back to multipotent hematopoietic stem cells. Cytogenetic techniques are the established baseline for diagnosis and characterization of complex hematological malignancies. In this study, we develop a workflow called Hema-seq to delineate clonal changes across various hematopoietic lineages through the integration of whole-genome sequencing, copy-number variations, cell morphology, and cytogenetic aberrations. In Hema-seq, cells are selected from Wright-stained slides and fluorescent probe-stained slides for sequencing. This technique therefore enables direct linking of whole-genome sequences to cytogenetic profiles. Through this method, we mapped sequential clonal alterations within the hematopoietic lineage, identifying critical shifts leading to myeloma and acute myeloid leukemia (AML) cell formations. By synthesizing data from each cell lineage, we provided insights into the hematopoietic tree's clonal evolution. Overall, this study highlights Hema-seq's capability in deciphering genomic heterogeneity in complex hematological malignancies, which can enable better diagnosis and treatment strategies.

Barcoded multiple displacement amplification for high coverage sequencing in spatial genomics.

Determining mutational landscapes in a spatial context is essential for understanding genetically heterogeneous cell microniches. Current approaches, such as Multiple Displacement Amplification (MDA), offer high genome coverage but limited multiplexing, which hinders large-scale spatial genomic studies. Here, we introduce barcoded MDA (bMDA), a technique that achieves high-coverage genomic analysis of low-input DNA while enhancing the multiplexing capabilities. By incorporating cell barcodes during MDA, bMDA streamlines library preparation in one pot, thereby overcoming a key bottleneck in spatial genomics. We apply bMDA to the integrative spatial analysis of triple-negative breast cancer tissues by examining copy number alterations, single nucleotide variations, structural variations, and kataegis signatures for each spatial microniche. This enables the assessment of subclonal evolutionary relationships within a spatial context. Therefore, bMDA has emerged as a scalable technology with the potential to advance the field of spatial genomics significantly.

Analysis of clinical and genomic profiles of therapy-related myeloid neoplasm in Korea.

BACKGROUND: Therapy-related myeloid neoplasm (T-MN) rarely occurs among cancer survivors, and was characterized by poor prognosis. T-MN has germline predisposition in a considerable proportion. Here, clinical characteristics and germline/somatic variant profiles in T-MN patients were investigated, and the findings were compared with those of previous studies. METHODS: A review of medical records, cytogenetic study, targeted sequencing by next-generation sequencing, and survival analysis were performed on 53 patients with T-MN at a single institution in Korea. RESULTS: The patients were relatively younger compared to T-MN patients in other studies. Our T-MN patients showed a high frequency of complex karyotypes, -5/del(5q), and -7/del(7q), which was similar to the Japanese study group but higher than the Australian study group. The most common primary disease was non-Hodgkin lymphoma, followed by breast cancer. The detailed distributions of primary diseases were different across study groups. Seven patients (13.2%) harbored deleterious presumed/potential germline variants in cancer predisposition genes (CPG) such as BRIP1, CEBPA, DDX41, FANCM, NBN, NF1, and RUNX1. In the somatic variant profile, TP53 was the most frequently mutated gene, which was consistent with the previous studies about T-MN. However, the somatic variant frequency in our study group was lower than in other studies. Adverse factors for overall survival were male sex, older age, history of previous radiotherapy, previous longer cytotoxic therapy, and -5/del(5q). CONCLUSION: The findings of our study corroborate important information about T-MN patients. As well as a considerable predisposition to CPG, the clinical characteristics and somatic variant profile showed distinctive patterns. Germline variant testing should be recommended for T-MN patients. If the T-MN patients harbor pathogenic germline variants, the family members for stem cell donation should be screened for carrier status through germline variant testing to avoid donor-derived myeloid neoplasm. For the prediction of the prognosis in T-MN patients, sex, age, past treatment history, and cytogenetic findings can be considered.

Telomere integrated scoring system of myelodysplastic syndrome.

INTRODUCTION: Recently, multigene target sequencing is widely performed for the purpose of prognostic prediction and application of targeted therapy. Here, we proposed a new scoring system that encompasses gene variations, telomere length, and Revised International Prognostic Scoring System (IPSS-R) together in Asian myelodysplastic syndrome. METHODS: We developed a new scoring model of these variables: age ≥ 65 years + IPSS-R score + ASXL1 mutation + TP53 mutation + Telomere length (<5.37). According to this new scoring system, patients were divided into four groups: very good score cutoff (≤3.0), good (3.0-4.5), poor (4.5-7.0), and very poor (>7.0). RESULTS: The median OS was 170.1, 100.4, 46.0, and 12.0 months for very good, good, poor, and very poor, retrospectively (p < 0.001). Meanwhile, according to the conventional IPSS-R scoring system, the median OS was 141.3, 50.2, 93.0, 36.0, and 16.2 months for very low, low, intermediate, high, and very high, retrospectively (p < 0.001). CONCLUSIONS: The newly developed model incorporating molecular variations and TL yielded more clear separations of the survival curves. By adding the presence of gene mutation and telomere length to the existing IPSS-R, its predictive ability can be further improved in myelodysplastic syndrome.

Novel THPO variant in hereditary thrombocytopenia: A potential candidate variant for predisposition to myeloid neoplasm.

Hereditary thrombocytopenia is a heterogeneous group of congenital disorders with a wide range of symptoms depending on the severity of platelet dysfunction or thrombocytopenia. Because of its clinical phenotypes and the bone marrow morphology associated with this condition, hereditary thrombocytopenia can be misdiagnosed as primary immune thrombocytopenia and myelodysplastic syndrome. Therefore, genetic evidence is necessary for the accurate diagnosis of hereditary thrombocytopenia. Refractory cytopenia of childhood is a subgroup of myelodysplastic syndrome that was added to the World Health Organization classification in 2008. To investigate the germline and somatic variants associated with refractory cytopenia of childhood, we performed targeted multigene sequencing in three patients with refractory cytopenia of childhood. Of the three patients, one progressed from megakaryocytic hypoplasia with thrombocytopenia, and targeted multigene sequencing revealed THPO variants in this patient and his sister. We propose that the monoallelic deletion of THPO is a potential candidate for germline predisposition to myeloid malignancy.

Accurate Diagnosis of COVID-19 from Self-Collectable Biospecimens Using Synthetic Apolipoprotein H Peptide-Coated Nanoparticle Assay.

A high-throughput, accurate screening is crucial for the prevention and control of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Current methods, which involve sampling from the nasopharyngeal (NP) area by medical staffs, constitute a fundamental bottleneck in expanding the testing capacity. To meet the scales required for population-level surveillance, self-collectable specimens can be used; however, its low viral load has hindered its clinical adoption. Here, we describe a magnetic nanoparticle functionalized with synthetic apolipoprotein H (ApoH) peptides to capture, concentrate, and purify viruses. The ApoH assay demonstrates a viral enrichment efficiency of >90% for both SARS-CoV-2 and its variants, leading to an order of magnitude improvement in analytical sensitivity. For validation, we apply the assay to a total of 84 clinical specimens including nasal, oral, and mouth gargles obtained from COVID-19 patients. As a result, a 100% positivity rate is achieved from the patient-collected nasal and gargle samples, which exceeds that of the traditional NP swab method. The simple 12 min pre-enrichment assay enabling the use of self-collectable samples will be a practical solution to overcome the overwhelming diagnostic capacity. Furthermore, the methodology can easily be built on various clinical protocols, allowing its broad applicability to various disease diagnoses.

The importance of critically short telomere in myelodysplastic syndrome.

A few critically short telomeres trigger genomic instability regardless of average telomere length (TL). Recently, the telomere shortest length assay (TeSLA) was developed to detect critically short telomeres and measure absolute telomeres. Using TeSLA with the internally labeled biotin probe, we measured the TL of bone marrow (BM) aspirates from 52 patients with myelodysplastic syndrome (MDS). A percentage of shortest telomeres (< 1.0 kb (ShTL1.0)) were calculated. ShTL1.0 was correlated to IPSS-R risk (spearman's rho = 0.35 and p = 0.0196), and ShTL1.0 and BM blast (2.61% in < 5% blast, 4.15% in 5-10% blast, and 6.80% in 10-20% blast, respectively, p = 0.0332). Interestingly, MDS patients with a shortest TL ≥ 0.787 kb at the time of diagnosis showed better overall survival (OS) and progression-free survival (PFS) than patients with a shortest TL < 0.787 kb in the multivariate analyses (HR = 0.13 and 0.30, p = 0.011 and 0.048 for OS and PFS, respectively). Our results clearly show the presence and abundance of critically short telomeres in MDS patients. These pathologic telomeres are associated with IPSS-R which is a validated prognostic scoring system in MDS. Furthermore, they are independent prognostic factors for OS in MDS patients. Future prospective studies are needed to validate our results.

Derivation of prognostic contextual histopathological features from whole-slide images of tumours via graph deep learning.

Methods of computational pathology applied to the analysis of whole-slide images (WSIs) do not typically consider histopathological features from the tumour microenvironment. Here, we show that a graph deep neural network that considers such contextual features in gigapixel-sized WSIs in a semi-supervised manner can provide interpretable prognostic biomarkers. We designed a neural-network model that leverages attention techniques to learn features of the heterogeneous tumour microenvironment from memory-efficient representations of aggregates of highly correlated image patches. We trained the model with WSIs of kidney, breast, lung and uterine cancers and validated it by predicting the prognosis of 3,950 patients with these four different types of cancer. We also show that the model provides interpretable contextual features of clear cell renal cell carcinoma that allowed for the risk-based retrospective stratification of 1,333 patients. Deep graph neural networks that derive contextual histopathological features from WSIs may aid diagnostic and prognostic tasks.

Spatial epitranscriptomics reveals A-to-I editome specific to cancer stem cell microniches.

Epitranscriptomic features, such as single-base RNA editing, are sources of transcript diversity in cancer, but little is understood in terms of their spatial context in the tumour microenvironment. Here, we introduce spatial-histopathological examination-linked epitranscriptomics converged to transcriptomics with sequencing (Select-seq), which isolates regions of interest from immunofluorescence-stained tissue and obtains transcriptomic and epitranscriptomic data. With Select-seq, we analyse the cancer stem cell-like microniches in relation to the tumour microenvironment of triple-negative breast cancer patients. We identify alternative splice variants, perform complementarity-determining region analysis of infiltrating T cells and B cells, and assess adenosine-to-inosine base editing in tumour tissue sections. Especially, in triple-negative breast cancer microniches, adenosine-to-inosine editome specific to different microniche groups is identified.

Advances in Tumor Sampling and Sequencing in Breast Cancer and their Application in Precision Diagnostics and Therapeutics.

Intra- and Inter-tumoral heterogeneity is one of the main hurdles in diagnosing and treating breast cancer. Selecting, sampling, and sequencing the samples appropriately provide unique opportunities in realizing precision medicine. This chapter reviews some of the past landmarks, state-of-the-art technologies, and future directions of translational research in terms of tumor sampling technologies and sequencing in breast cancer. In the state-of-the-art technologies section, the technologies are categorized in terms of scientific, precision diagnostic, and precision therapeutic tools. Finally, limitations and future directions regarding various translational research for clinical applications using these technologies will be discussed.

Phenotype-based single cell sequencing identifies diverse genetic subclones in CD133 positive cancer stem cells.

Tumor heterogeneity is one of the ongoing huddles in the field of colon cancer therapy. It is evident that there are countless clones which exhibit different phenotypes and therefore, single cell analysis is inevitable. Cancer stem cells (CSCs) are rare cell population within tumor which is known to function in cancer metastasis and recurrence. Although there have been trials to prove intra-tumoral heterogeneity using single cell sequencing, that of CSCs has not been clearly elucidated. Here, we articulate the presence of heterogeneous subclones within CD133 positive cancer stem cells through single cell sequencing. As a proof of principle, we performed phenotype-based high-throughput laser isolation and single cell sequencing (PHLI-seq) of CD133 positive cells in a frozen tumor tissue obtained from a patient with colorectal cancer. The result proved that CD133 positive cells were shown to be heterogeneous both in copy number and mutational profiles. Single cancer stem cell specific mutations such as RNF144A, PAK2, PARP4, ADAM21, HYDIN, KRT38 and CELSR1 could be also detected in liver metastatic tumor of the same patient. Collectively, these data suggest that single cell analysis used to spot subclones with genetic variation within rare population, will lead to new strategies to tackle colon cancer metastasis.

OPENchip: an on-chip in situ molecular profiling platform for gene expression analysis and oncogenic mutation detection in single circulating tumour cells.

Liquid biopsy holds promise towards practical implementation of personalized theranostics of cancer. In particular, circulating tumour cells (CTCs) can provide clinically actionable information that can be directly linked to prognosis or therapy decisions. In this study, gene expression patterns and genetic mutations in single CTCs are simultaneously analysed by strategically combining microfluidic technology and in situ molecular profiling technique. Towards this, the development and demonstration of the OPENchip (On-chip Post-processing ENabling chip) platform for single CTC analysis by epithelial CTC enrichment and subsequent in situ molecular profiling is reported. For in situ molecular profiling, padlock probes that identify specific desired targets to examine biomarkers of clinical relevance in cancer diagnostics were designed and used to create libraries of rolling circle amplification products. We characterize the OPENchip in terms of its capture efficiency and capture purity, and validate the probe design using different cell lines. By integrating the obtained results, molecular analyses of CTCs from metastatic breast cancer (HER2 (ERBB2) gene expression and PIK3CA mutations) and metastatic pancreatic cancer (KRAS gene mutations) patients were demonstrated without any off-chip processes. The results substantiate the potential implementation of early molecular detection of cancer through sequencing-free liquid biopsy.

Targeted sequencing aids in identifying clonality in chronic myelomonocytic leukemia.

Chronic myelomonocytic leukemia (CMML) typically shows monocytosis in the peripheral blood (PB), which must be differentiated from reactive monocytosis. To determine the clonality of CMML, we performed molecular and cytogenetic analysis in Korean patients. To investigate whether monocytes in the PB harbored clonal mutational changes, we performed single-cell sequencing after selecting monocytes, neutrophils, and lymphocytes by morphology-aided laser microdissection. Targeted sequencing was performed in 35 patients with CMML with 41 bone marrow samples. Single-cell analysis was performed in two cases. Most (94.3%) patients harbored at least one variant, in genes considered as potential therapeutic targets, while cytogenetic aberrations occurred in only 28.6% of cases. ASXL1 (54.3%), SRSF2 (37.1%), NRAS (31.4%), and TET2 (25.7%) were frequently mutated, with lower frequencies of TET2 mutation and higher frequencies of NRAS, DNMT3A (17.1%), and NPM1 (11.4%) mutations compared to in previous studies of Caucasians. Patients with SETBP1 mutation and those with more than two variants showed poorer survival than those without mutation (P < 0.001 and P = 0.007, respectively). Most (70.8%) variants were detected at diagnosis and follow-up with no significant differences in variant allele frequency, warranting sequencing during follow-up if diagnostic samples were unavailable. Single-cell analysis revealed clonal monocytes with mutations, and the same mutations were also identified in lymphocytes and neutrophils. Targeted sequencing aided in clonality detection in most patients with CMML and single-cell sequencing facilitated identification of clonal monocytes and the co-existence of mutations in non-myeloid cells, suggesting that certain mutations are acquired by pluripotent stem cells.

Characteristics of Waldenström Macroglobulinemia in Korean Patients According to Mutational Status of MYD88 and CXCR4: Analysis Using Ultra-Deep Sequencing.

BACKGROUND: Little is known about the mutational frequency of myeloid differentiation factor 88 (MYD88) and C-X-C chemokine receptor type 4 (CXCR4) and the corresponding characteristics in Asian individuals afflicted with Waldenström macroglobulinemia (WM). We investigated the characteristics of WM according to mutational status of MYD88/CXCR4, and attempted to determine the lineage commitment among hematopoietic cells by MYD88L265P single-cell sequencing on bone marrow (BM) smear slides. MATERIALS AND METHODS: CXCR4 mutations (muts) were detected using ultra-deep sequencing using target capture. Mutational burden of MYD88 was assessed using real-time polymerase chain reaction. Single-cell sequencing for MYD88 was performed on lymphocytes, plasmacytoid lymphocytes, plasma cells, and neutrophils using laser microdissection. RESULTS: Among 31 patients, the frequencies of MYD88/CXCR4 muts were as follows: MYD88 wild type (WT) CXCR4WT (6 patients, 19.4%), MYD88L265PCXCR4WT (19 patients, 61.4%), MYD88L265PCXCR4mut (6 patients, 19.4%; 1 frameshift and 5 nonsense muts). Immunoglobulin M levels of MYD88L265CXCR4WT patients were significantly higher than those of MYD88WTCXCR4WT patients (P = .024). Tumor burden in BM was highest in patients with MYD88L265PCXCR4mut (82.0%), followed by MYD88L265PCXCR4WT (52.8%) and MYD88WTCXCR4WT (14.2%) (P < .001). The quantity of MYD88-mutated DNA tended to correlate with tumor burden in BM (correlation coefficient 0.647; P = .009). MYD88L265P was detected in plasma cells, plasmacytoid lymphocytes, and lymphocytes but not neutrophils. CONCLUSION: The frequency of MYD88/CXCR4 muts in Korean and Caucasian patients with WM was similar, however 5 of the 6 CXCR4 muts were nonsense-a proportion higher than reported frequencies in Caucasian individuals. Ultra-deep sequencing was capable of detecting CXCR4 muts not detectable using Sanger sequencing, suggesting a possible replacement of the B-cell sorting.

Whole Genome Sequencing of Single Circulating Tumor Cells Isolated by Applying a Pulsed Laser to Cell-Capturing Microstructures.

Single cell analysis of heterogeneous circulating tumor cells (CTCs), by which the genomic profiles of rare single CTCs are connected to the clinical status of cancer patients, is crucial for understanding cancer metastasis and the clinical impact on patients. However, the heterogeneity in genotypes and phenotypes and rarity of CTCs have limited extensive single CTC genome research, further hindering clinical investigation. Despite recent efforts to build platforms that separate CTCs, the investigation on CTCs is difficult due to the lack of a retrieval process at the single cell level. In this study, laser-induced isolation of microstructures on an optomechanically-transferrable-chip and sequencing (LIMO-seq) is applied for whole genome sequencing of single CTCs. Also, the whole genome sequences and the molecular profiles of the isolated single cells from the whole blood of a breast cancer patient are analyzed.

Efficient Selection of Antibodies Reactive to Homologous Epitopes on Human and Mouse Hepatocyte Growth Factors by Next-Generation Sequencing-Based Analysis of the B Cell Repertoire.

: YYB-101 is a humanized rabbit anti-human hepatocyte growth factor (HGF)-neutralizing antibody currently in clinical trial. To test the effect of HGF neutralization with antibody on anti-cancer T cell immunity, we generated surrogate antibodies that are reactive to the mouse homologue of the epitope targeted by YYB-101. First, we immunized a chicken with human HGF and monitored changes in the B cell repertoire by next-generation sequencing (NGS). We then extracted the VH gene repertoire from the NGS data, clustered it into components by sequence homology, and classified the components by the change in the number of unique VH sequences and the frequencies of the VH sequences within each component following immunization. Those changes should accompany the preferential proliferation and somatic hypermutation or gene conversion of B cells encoding HGF-reactive antibodies. One component showed significant increases in the number and frequencies of unique VH sequences and harbored genes encoding antibodies that were reactive to human HGF and competitive with YYB-101 for HGF binding. Some of the antibodies also reacted to mouse HGF. The selected VH sequences shared 98.3% identity and 98.9% amino acid similarity. It is therefore likely that the antibodies encoded by them all react to the epitope targeted by YYB-101.

Author Correction: Evaluating Tumor Evolution via Genomic Profiling of Individual Tumor Spheroids in a Malignant Ascites.

A correction to this article has been published and is linked from the HTML and PDF versions of this paper. The error has not been fixed in the paper.

PHLI-seq: constructing and visualizing cancer genomic maps in 3D by phenotype-based high-throughput laser-aided isolation and sequencing.

Spatial mapping of genomic data to tissue context in a high-throughput and high-resolution manner has been challenging due to technical limitations. Here, we describe PHLI-seq, a novel approach that enables high-throughput isolation and genome-wide sequence analysis of single cells or small numbers of cells to construct genomic maps within cancer tissue in relation to the images or phenotypes of the cells. By applying PHLI-seq, we reveal the heterogeneity of breast cancer tissues at a high resolution and map the genomic landscape of the cells to their corresponding spatial locations and phenotypes in the 3D tumor mass.