SIDR: simultaneous isolation and parallel sequencing of genomic DNA and total RNA from single cells.

Simultaneous sequencing of the genome and transcriptome at the single-cell level is a powerful tool for characterizing genomic and transcriptomic variation and revealing correlative relationships. However, it remains technically challenging to analyze both the genome and transcriptome in the same cell. Here, we report a novel method for simultaneous isolation of genomic DNA and total RNA (SIDR) from single cells, achieving high recovery rates with minimal cross-contamination, as is crucial for accurate description and integration of the single-cell genome and transcriptome. For reliable and efficient separation of genomic DNA and total RNA from single cells, the method uses hypotonic lysis to preserve nuclear lamina integrity and subsequently captures the cell lysate using antibody-conjugated magnetic microbeads. Evaluating the performance of this method using real-time PCR demonstrated that it efficiently recovered genomic DNA and total RNA. Thorough data quality assessments showed that DNA and RNA simultaneously fractionated by the SIDR method were suitable for genome and transcriptome sequencing analysis at the single-cell level. The integration of single-cell genome and transcriptome sequencing by SIDR (SIDR-seq) showed that genetic alterations, such as copy-number and single-nucleotide variations, were more accurately captured by single-cell SIDR-seq compared with conventional single-cell RNA-seq, although copy-number variations positively correlated with the corresponding gene expression levels. These results suggest that SIDR-seq is potentially a powerful tool to reveal genetic heterogeneity and phenotypic information inferred from gene expression patterns at the single-cell level.

Hypomethylation of ATP1A1 Is Associated with Poor Prognosis and Cancer Progression in Triple-Negative Breast Cancer.

Dysregulated DNA methylation in cancer is critical in the transcription machinery associated with cancer progression. Triple-negative breast cancer (TNBC) is the most aggressive breast cancer subtype, but no treatment targeting TNBC biomarkers has yet been developed. To identify specific DNA methylation patterns in TNBC, methyl-binding domain protein 2 (MBD) sequencing data were compared in TNBC and the three other major breast cancer subtypes. Integrated analysis of DNA methylation and gene expression identified a gene set showing a correlation between DNA methylation and gene expression. ATPase Na+/K+-transporting subunit alpha 1 (ATP1A1) was found to be specifically hypomethylated in the coding sequence (CDS) region and to show increased expression in TNBC. The Cancer Genome Atlas (TCGA) database also showed that hypomethylation and high expression of ATP1A1 were strongly associated with poor survival in patients with TNBC. Furthermore, ATP1A1 knockdown significantly reduced the viability and tumor-sphere formation of TNBC cells. These results suggest that the hypomethylation and overexpression of ATP1A1 could be a prognostic marker in TNBC and that the manipulation of ATP1A1 expression could be a therapeutic target in this disease.

Development of a Breast Cancer Risk Prediction Model Incorporating Polygenic Risk Scores and Nongenetic Risk Factors for Korean Women.

BACKGROUND: To develop a breast cancer prediction model for Korean women using published polygenic risk scores (PRS) combined with nongenetic risk factors (NGRF). METHODS: Thirteen PRS models generated from single or multiple combinations of the Asian and European PRSs were evaluated among 20,434 Korean women. The AUC and increase in OR per SD were compared for each PRS. The PRSs with the highest predictive power were combined with NGRFs; then, an integrated prediction model was established using the Individualized Coherent Absolute Risk Estimation (iCARE) tool. The absolute breast cancer risk was stratified for 18,142 women with available follow-up data. RESULTS: PRS38_ASN+PRS190_EB, a combination of Asian and European PRSs, had the highest AUC (0.621) among PRSs, with an OR per SD increase of 1.45 (95% confidence interval: 1.31-1.61). Compared with the average risk group (35%-65%), women in the top 5% had a 2.5-fold higher risk of breast cancer. Incorporating NGRFs yielded a modest increase in the AUC of women ages >50 years. For PRS38_ASN+PRS190_EB+NGRF, the average absolute risk was 5.06%. The lifetime absolute risk at age 80 years for women in the top 5% was 9.93%, whereas that of women in the lowest 5% was 2.22%. Women at higher risks were more sensitive to NGRF incorporation. CONCLUSIONS: Combined Asian and European PRSs were predictive of breast cancer in Korean women. Our findings support the use of these models for personalized screening and prevention of breast cancer. IMPACT: Our study provides insights into genetic susceptibility and NGRFs for predicting breast cancer in Korean women.

Tumor-promoting macrophages prevail in malignant ascites of advanced gastric cancer.

Gastric cancer (GC) patients develop malignant ascites as the disease progresses owing to peritoneal metastasis. GC patients with malignant ascites have a rapidly deteriorating clinical course with short survival following the onset of malignant ascites. Better optimized treatment strategies for this subset of patients are needed. To define the cellular characteristics of malignant ascites of GC, we used single-cell RNA sequencing to characterize tumor cells and tumor-associated macrophages (TAMs) from four samples of malignant ascites and one sample of cerebrospinal fluid. Reference transcriptomes for M1 and M2 macrophages were generated by in vitro differentiation of healthy blood-derived monocytes and applied to assess the inflammatory properties of TAMs. We analyzed 180 cells, including tumor cells, macrophages, and mesothelial cells. Dynamic exchange of tumor-promoting signals, including the CCL3-CCR1 or IL1B-IL1R2 interactions, suggests macrophage recruitment and anti-inflammatory tuning by tumor cells. By comparing these data with reference transcriptomes for M1-type and M2-type macrophages, we found noninflammatory characteristics in macrophages recovered from the malignant ascites of GC. Using public datasets, we demonstrated that the single-cell transcriptome-driven M2-specific signature was associated with poor prognosis in GC. Our data indicate that the anti-inflammatory characteristics of TAMs are controlled by tumor cells and present implications for treatment strategies for GC patients in which combination treatment targeting cancer cells and macrophages may have a reciprocal synergistic effect.

Single-cell RNA sequencing reveals the tumor microenvironment and facilitates strategic choices to circumvent treatment failure in a chemorefractory bladder cancer patient.

BACKGROUND: Tumor cell-intrinsic mechanisms and complex interactions with the tumor microenvironment contribute to therapeutic failure via tumor evolution. It may be possible to overcome treatment resistance by developing a personalized approach against relapsing cancers based on a comprehensive analysis of cell type-specific transcriptomic changes over the clinical course of the disease using single-cell RNA sequencing (scRNA-seq). METHODS: Here, we used scRNA-seq to depict the tumor landscape of a single case of chemo-resistant metastatic, muscle-invasive urothelial bladder cancer (MIUBC) addicted to an activating Harvey rat sarcoma viral oncogene homolog (HRAS) mutation. In order to analyze tumor evolution and microenvironmental changes upon treatment, we also applied scRNA-seq to the corresponding patient-derived xenograft (PDX) before and after treatment with tipifarnib, a HRAS-targeting agent under clinical evaluation. RESULTS: In the parallel analysis of the human MIUBC and the PDX, diverse stromal and immune cell populations recapitulated the cellular composition in the human and mouse tumor microenvironment. Treatment with tipifarnib showed dramatic anticancer effects but was unable to achieve a complete response. Importantly, the comparative scRNA-seq analysis between pre- and post-tipifarnib-treated PDX revealed the nature of tipifarnib-refractory tumor cells and the tumor-supporting microenvironment. Based on the upregulation of programmed death-ligand 1 (PD-L1) in surviving tumor cells, and the accumulation of multiple immune-suppressive subsets from post-tipifarnib-treated PDX, a PD-L1 inhibitor, atezolizumab, was clinically applied; this resulted in a favorable response from the patient with acquired resistance to tipifarnib. CONCLUSION: We presented a single case report demonstrating the power of scRNA-seq for visualizing the tumor microenvironment and identifying molecular and cellular therapeutic targets in a treatment-refractory cancer patient.

Lineage-dependent gene expression programs influence the immune landscape of colorectal cancer.

Immunotherapy for metastatic colorectal cancer is effective only for mismatch repair-deficient tumors with high microsatellite instability that demonstrate immune infiltration, suggesting that tumor cells can determine their immune microenvironment. To understand this cross-talk, we analyzed the transcriptome of 91,103 unsorted single cells from 23 Korean and 6 Belgian patients. Cancer cells displayed transcriptional features reminiscent of normal differentiation programs, and genetic alterations that apparently fostered immunosuppressive microenvironments directed by regulatory T cells, myofibroblasts and myeloid cells. Intercellular network reconstruction supported the association between cancer cell signatures and specific stromal or immune cell populations. Our collective view of the cellular landscape and intercellular interactions in colorectal cancer provide mechanistic information for the design of efficient immuno-oncology treatment strategies.

Alternative polyadenylation of single cells delineates cell types and serves as a prognostic marker in early stage breast cancer.

Alternative polyadenylation (APA) in 3' untranslated regions (3' UTR) plays an important role in regulating transcript abundance, localization, and interaction with microRNAs. Length-variation of 3'UTRs by APA contributes to efficient proliferation of cancer cells. In this study, we investigated APA in single cancer cells and tumor microenvironment cells to understand the physiological implication of APA in different cell types. We analyzed APA patterns and the expression level of genes from the 515 single-cell RNA sequencing (scRNA-seq) dataset from 11 breast cancer patients. Although the overall 3'UTR length of individual genes was distributed equally in tumor and non-tumor cells, we found a differential pattern of polyadenylation in gene sets between tumor and non-tumor cells. In addition, we found a differential pattern of APA across tumor types using scRNA-seq data from 3 glioblastoma patients and 1 renal cell carcinoma patients. In detail, 1,176 gene sets and 53 genes showed the distinct pattern of 3'UTR shortening and over-expression as signatures for five cell types including B lymphocytes, T lymphocytes, myeloid cells, stromal cells, and breast cancer cells. Functional categories of gene sets for cellular proliferation demonstrated concordant regulation of APA and gene expression specific to cell types. The expression of APA genes in breast cancer was significantly correlated with the clinical outcome of earlier stage breast cancer patients. We identified cell type-specific APA in single cells, which allows the identification of cell types based on 3'UTR length variation in combination with gene expression. Specifically, an immune-specific APA signature in breast cancer could be utilized as a prognostic marker of early stage breast cancer.

Multi-omics profiling of younger Asian breast cancers reveals distinctive molecular signatures.

Breast cancer (BC) in the Asia Pacific regions is enriched in younger patients and rapidly rising in incidence yet its molecular bases remain poorly characterized. Here we analyze the whole exomes and transcriptomes of 187 primary tumors from a Korean BC cohort (SMC) enriched in pre-menopausal patients and perform systematic comparison with a primarily Caucasian and post-menopausal BC cohort (TCGA). SMC harbors higher proportions of HER2+ and Luminal B subtypes, lower proportion of Luminal A with decreased ESR1 expression compared to TCGA. We also observe increased mutation prevalence affecting BRCA1, BRCA2, and TP53 in SMC with an enrichment of a mutation signature linked to homologous recombination repair deficiency in TNBC. Finally, virtual microdissection and multivariate analyses reveal that Korean BC status is independently associated with increased TIL and decreased TGF-ß signaling expression signatures, suggesting that younger Asian BCs harbor more immune-active microenvironment than western BCs.

Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer.

Single-cell transcriptome profiling of tumour tissue isolates allows the characterization of heterogeneous tumour cells along with neighbouring stromal and immune cells. Here we adopt this powerful approach to breast cancer and analyse 515 cells from 11 patients. Inferred copy number variations from the single-cell RNA-seq data separate carcinoma cells from non-cancer cells. At a single-cell resolution, carcinoma cells display common signatures within the tumour as well as intratumoral heterogeneity regarding breast cancer subtype and crucial cancer-related pathways. Most of the non-cancer cells are immune cells, with three distinct clusters of T lymphocytes, B lymphocytes and macrophages. T lymphocytes and macrophages both display immunosuppressive characteristics: T cells with a regulatory or an exhausted phenotype and macrophages with an M2 phenotype. These results illustrate that the breast cancer transcriptome has a wide range of intratumoral heterogeneity, which is shaped by the tumour cells and immune cells in the surrounding microenvironment.

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.