TSpred: a robust prediction framework for TCR-epitope interactions using paired chain TCR sequence data.

MOTIVATION: Prediction of T-cell receptor (TCR)-epitope interactions is important for many applications in biomedical research, such as cancer immunotherapy and vaccine design. The prediction of TCR-epitope interactions remains challenging especially for novel epitopes, due to the scarcity of available data. RESULTS: We propose TSpred, a new deep learning approach for the pan-specific prediction of TCR binding specificity based on paired chain TCR data. We develop a robust model that generalizes well to unseen epitopes by combining the predictive power of CNN and the attention mechanism. In particular, we design a reciprocal attention mechanism which focuses on extracting the patterns underlying TCR-epitope interactions. Upon a comprehensive evaluation of our model, we find that TSpred achieves state-of-the-art performances in both seen and unseen epitope specificity prediction tasks. Also, compared to other predictors, TSpred is more robust to bias related to peptide imbalance in the dataset. In addition, the reciprocal attention component of our model allows for model interpretability by capturing structurally important binding regions. Results indicate that TSpred is a robust and reliable method for the task of TCR-epitope binding prediction. AVAILABILITY: Source code is available at https://github.com/ha01994/TSpred. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Biomarker-driven phase 2 umbrella trial: Clinical efficacy of olaparib monotherapy and combination with ceralasertib (AZD6738) in small cell lung cancer.

BACKGROUND: Based on a high incidence of genomic alteration in the cell cycle and DNA damage and response (DDR)-related pathways in small cell lung cancer (SCLC), the clinical efficacy of the DDR-targeting agent olaparib (PARP inhibitor) as monotherapy and in combination with ceralasertib (ATR inhibitor) in relapsed or refractory SCLC was evaluated. METHODS: As part of a phase 2 biomarker driven umbrella study, patients with SCLC and predefined DDR gene alterations who failed to benefit from prior platinum-based regimens were allocated to the olaparib monotherapy arm and nonbiomarker-selected patients were allocated to the olaparib and ceralasertib combination arm. RESULTS: In the olaparib monotherapy arm (n = 15), the objective response rate was 6.7% (one partial response), and the disease control rate was 33.3%, including three patients with stable disease. The median progression-free survival was 1.3 months (95% CI, 1.2-NA). In the combination arm (n = 26), the objective response rate and disease control rate were 3.8% and 42.3%, respectively, with one partial response and 10 patients with stable disease. The median progression-free survival was 2.8 months (95% CI, 1.8-5.4). Treatment was generally well tolerated except for one fatal case of neutropenic fever in the combination arm. CONCLUSIONS: Targeting DDR pathways with olaparib as a single agent or in combination with ceralasertib did not meet the predefined efficacy end point. However, disease stabilization was more evident in the combination arm. Further investigation of the combination of olaparib in SCLC should be performed with diverse combinations and patient selection strategies to maximize efficacy.

Polygenic risk score-based phenome-wide association study of head and neck cancer across two large biobanks.

BACKGROUND: Numerous observational studies have highlighted associations of genetic predisposition of head and neck squamous cell carcinoma (HNSCC) with diverse risk factors, but these findings are constrained by design limitations of observational studies. In this study, we utilized a phenome-wide association study (PheWAS) approach, incorporating a polygenic risk score (PRS) derived from a wide array of genomic variants, to systematically investigate phenotypes associated with genetic predisposition to HNSCC. Furthermore, we validated our findings across heterogeneous cohorts, enhancing the robustness and generalizability of our results. METHODS: We derived PRSs for HNSCC and its subgroups, oropharyngeal cancer and oral cancer, using large-scale genome-wide association study summary statistics from the Genetic Associations and Mechanisms in Oncology Network. We conducted a comprehensive investigation, leveraging genotyping data and electronic health records from 308,492 individuals in the UK Biobank and 38,401 individuals in the Penn Medicine Biobank (PMBB), and subsequently performed PheWAS to elucidate the associations between PRS and a wide spectrum of phenotypes. RESULTS: We revealed the HNSCC PRS showed significant association with phenotypes related to tobacco use disorder (OR, 1.06; 95% CI, 1.05-1.08; P = 3.50 × 10-15), alcoholism (OR, 1.06; 95% CI, 1.04-1.09; P = 6.14 × 10-9), alcohol-related disorders (OR, 1.08; 95% CI, 1.05-1.11; P = 1.09 × 10-8), emphysema (OR, 1.11; 95% CI, 1.06-1.16; P = 5.48 × 10-6), chronic airway obstruction (OR, 1.05; 95% CI, 1.03-1.07; P = 2.64 × 10-5), and cancer of bronchus (OR, 1.08; 95% CI, 1.04-1.13; P = 4.68 × 10-5). These findings were replicated in the PMBB cohort, and sensitivity analyses, including the exclusion of HNSCC cases and the major histocompatibility complex locus, confirmed the robustness of these associations. Additionally, we identified significant associations between HNSCC PRS and lifestyle factors related to smoking and alcohol consumption. CONCLUSIONS: The study demonstrated the potential of PRS-based PheWAS in revealing associations between genetic risk factors for HNSCC and various phenotypic traits. The findings emphasized the importance of considering genetic susceptibility in understanding HNSCC and highlighted shared genetic bases between HNSCC and other health conditions and lifestyles.

Association between genetic risk and adherence to healthy lifestyle for developing age-related hearing loss.

BACKGROUND: Previous studies have shown that lifestyle/environmental factors could accelerate the development of age-related hearing loss (ARHL). However, there has not yet been a study investigating the joint association among genetics, lifestyle/environmental factors, and adherence to healthy lifestyle for risk of ARHL. We aimed to assess the association between ARHL genetic variants, lifestyle/environmental factors, and adherence to healthy lifestyle as pertains to risk of ARHL. METHODS: This case-control study included 376,464 European individuals aged 40 to 69 years, enrolled between 2006 and 2010 in the UK Biobank (UKBB). As a replication set, we also included a total of 26,523 individuals considered of European ancestry and 9834 individuals considered of African-American ancestry through the Penn Medicine Biobank (PMBB). The polygenic risk score (PRS) for ARHL was derived from a sensorineural hearing loss genome-wide association study from the FinnGen Consortium and categorized as low, intermediate, high, and very high. We selected lifestyle/environmental factors that have been previously studied in association with hearing loss. A composite healthy lifestyle score was determined using seven selected lifestyle behaviors and one environmental factor. RESULTS: Of the 376,464 participants, 87,066 (23.1%) cases belonged to the ARHL group, and 289,398 (76.9%) individuals comprised the control group in the UKBB. A very high PRS for ARHL had a 49% higher risk of ARHL than those with low PRS (adjusted OR, 1.49; 95% CI, 1.36-1.62; P < .001), which was replicated in the PMBB cohort. A very poor lifestyle was also associated with risk of ARHL (adjusted OR, 3.03; 95% CI, 2.75-3.35; P < .001). These risk factors showed joint effects with the risk of ARHL. Conversely, adherence to healthy lifestyle in relation to hearing mostly attenuated the risk of ARHL even in individuals with very high PRS (adjusted OR, 0.21; 95% CI, 0.09-0.52; P < .001). CONCLUSIONS: Our findings of this study demonstrated a significant joint association between genetic and lifestyle factors regarding ARHL. In addition, our analysis suggested that lifestyle adherence in individuals with high genetic risk could reduce the risk of ARHL.

Discovering comorbid diseases using an inter-disease interactivity network based on biobank-scale PheWAS data.

MOTIVATION: Understanding comorbidity is essential for disease prevention, treatment and prognosis. In particular, insight into which pairs of diseases are likely or unlikely to co-occur may help elucidate the potential relationships between complex diseases. Here, we introduce the use of an inter-disease interactivity network to discover/prioritize comorbidities. Specifically, we determine disease associations by accounting for the direction of effects of genetic components shared between diseases, and categorize those associations as synergistic or antagonistic. We further develop a comorbidity scoring algorithm to predict whether diseases are more or less likely to co-occur in the presence of a given index disease. This algorithm can handle networks that incorporate relationships with opposite signs. RESULTS: We finally investigate inter-disease associations among 427 phenotypes in UK Biobank PheWAS data and predict the priority of comorbid diseases. The predicted comorbidities were verified using the UK Biobank inpatient electronic health records. Our findings demonstrate that considering the interaction of phenotype associations might be helpful in better predicting comorbidity. AVAILABILITY AND IMPLEMENTATION: The source code and data of this study are available at https://github.com/dokyoonkimlab/DiseaseInteractiveNetwork. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

A nonsense TMEM43 variant leads to disruption of connexin-linked function and autosomal dominant auditory neuropathy spectrum disorder.

Genes that are primarily expressed in cochlear glia-like supporting cells (GLSs) have not been clearly associated with progressive deafness. Herein, we present a deafness locus mapped to chromosome 3p25.1 and an auditory neuropathy spectrum disorder (ANSD) gene, TMEM43, mainly expressed in GLSs. We identify p.(Arg372Ter) of TMEM43 by linkage analysis and exome sequencing in two large Asian families segregating ANSD, which is characterized by inability to discriminate speech despite preserved sensitivity to sound. The knock-in mouse with the p.(Arg372Ter) variant recapitulates a progressive hearing loss with histological abnormalities in GLSs. Mechanistically, TMEM43 interacts with the Connexin26 and Connexin30 gap junction channels, disrupting the passive conductance current in GLSs in a dominant-negative fashion when the p.(Arg372Ter) variant is introduced. Based on these mechanistic insights, cochlear implant was performed on three subjects, and speech discrimination was successfully restored. Our study highlights a pathological role of cochlear GLSs by identifying a deafness gene and its causal relationship with ANSD.

Potential Role of Macrophage Polarization in the Progression of Hunner-Type Interstitial Cystitis.

BACKGROUND: Hunner-type interstitial cystitis (HIC) is a chronic inflammatory condition of the bladder. However, it remains unclear whether there is a causal relationship between the presence of Hunner lesions and seemingly normal-appearing areas in the bladder (non-Hunner lesions). This study aimed to investigate the fundamental aspects of HIC by examining potential genetic differences between Hunner and non-Hunner lesions and elucidate their role as potential markers in the progression and suppression of the disease. METHODS: This cross-sectional study enrolled patients with HIC (n = 10) who underwent supratrigonal cystectomy along with augmentation cystoplasty. Full-thickness bladder tissue was collected from Hunner and non-Hunner lesions in the same patient. Normal bladder tissue biopsies were also obtained as controls. Whole transcriptome analysis was performed to analyze the gene expression patterns and immune cell populations. RESULTS: The mucosal layers of patients exhibited similar pathway dysregulation across Hunner and non-Hunner lesions, with immunerelated pathways being prominently affected. In the mucosal layer, genes related to anti-inflammatory and immune suppression were downregulated in Hunner lesions compared to non-Hunner lesions. Moreover, in Hunner lesions, genes related to macrophage differentiation and polarization, such as VSIG4, CD68, MAFB, and LIRB4, were downregulated. The cell fraction of M2 macrophages was found to decrease in Hunner lesions. Immunohistochemical staining revealed an elevated fraction of M1 macrophages and a reduced fraction of M2 macrophages in Hunner lesions compared to those in non-Hunner lesions. In the muscular layer, transcriptomic evidence of muscle thickness was observed in both Hunner and non-Hunner lesions; however, the difference was not significant. CONCLUSION: Hunner lesions showed a reduced expression of anti-inflammatory and immunosuppressive factors compared to non-Hunner lesions, along with alterations in immune cell populations. This study suggests the possibility that macrophage polarization is related to the progression from non-Hunner lesions to Hunner lesions, suggesting its relevance to the characteristics of autoimmune diseases.

Migratory Tumor Cells Cooperate with Cancer Associated Fibroblasts in Hormone Receptor-Positive and HER2-Negative Breast Cancer.

Hormone receptor-positive and HER2-negative breast cancer (HR+/HER2-BC) is the most common type with a favorable prognosis under endocrine therapy. However, it still demonstrates unpredictable progression and recurrences influenced by high tumoral diversity and microenvironmental status. To address these heterogeneous molecular characteristics of HR+/HER2-BC, we aimed to simultaneously characterize its transcriptomic landscape and genetic architecture at the same resolution. Using advanced single-cell RNA and DNA sequencing techniques together, we defined four distinct tumor subtypes. Notably, the migratory tumor subtype was closely linked to genomic alterations of EGFR, related to the tumor-promoting behavior of IL6-positive inflammatory tumor-associated fibroblast, and contributing to poor prognosis. Our study comprehensively utilizes integrated analysis to uncover the complex dynamics of this breast cancer subtype, highlighting the pivotal role of the migratory tumor subtype in influencing surrounding cells. This sheds light on potential therapeutic targets by offering enhanced insights for HR+/HER2-BC treatment.

Genome-wide association study of occupational attainment as a proxy for cognitive reserve.

Occupational attainment, which represents middle-age cognitive activities, is a known proxy marker of cognitive reserve for Alzheimer's disease. Previous genome-wide association studies have identified numerous genetic variants and revealed the genetic architecture of educational attainment, another marker of cognitive reserve. However, the genetic architecture and heritability for occupational attainment remain elusive. We performed a large-scale genome-wide association study of occupational attainment with 248 847 European individuals from the UK Biobank using the proportional odds logistic mixed model method. In this analysis, we defined occupational attainment using the classified job levels formulated in the UK Standard Occupational Classification system considering the individual professional skill and academic level. We identified 30 significant loci (P < 5 × 10-8); 12 were novel variants, not associated with other traits. Among them, four lead variants were associated with genes expressed in brain tissues by expression quantitative trait loci mapping from 10 brain regions: rs13002946, rs3741368, rs11654986 and rs1627527. The single nucleotide polymorphism-based heritability was estimated to be 8.5% (standard error of the mean = 0.004) and partitioned heritability was enriched in the CNS and brain tissues. Genetic correlation analysis showed shared genetic backgrounds between occupational attainment and multiple traits, including education, intelligence, leisure activities, life satisfaction and neuropsychiatric disorders. In two-sample Mendelian randomization analysis, we demonstrated that high occupation levels were associated with reduced risk for Alzheimer's disease [odds ratio (OR) = 0.78, 95% confidence interval (CI) = 0.65-0.92 in inverse variance weighted method; OR = 0.73, 95% CI = 0.57-0.92 in the weighted median method]. This causal relationship between occupational attainment and Alzheimer's disease was robust in additional sensitivity analysis that excluded potentially pleiotropic single nucleotide polymorphisms (OR = 0.72, 95% CI = 0.57-0.91 in the inverse variance weighted method; OR = 0.72, 95% CI = 0.53-0.97 in the weighted median method). Multivariable Mendelian randomization confirmed that occupational attainment had an independent effect on the risk for Alzheimer's disease even after taking educational attainment into account (OR = 0.72, 95% CI = 0.54-0.95 in the inverse variance weighted method; OR = 0.68, 95% CI = 0.48-0.97 in the weighted median method). Overall, our analyses provide insights into the genetic architecture of occupational attainment and demonstrate that occupational attainment is a potential causal protective factor for Alzheimer's disease as a proxy marker of cognitive reserve.

ICD2Vec: Mathematical representation of diseases.

BACKGROUND: The International Classification of Diseases (ICD) codes represent the global standard for reporting disease conditions. The current ICD codes connote direct human-defined relationships among diseases in a hierarchical tree structure. Representing the ICD codes as mathematical vectors helps to capture nonlinear relationships in medical ontologies across diseases. METHODS: We propose a universally applicable framework called "ICD2Vec" designed to provide mathematical representations of diseases by encoding corresponding information. First, we present the arithmetical and semantic relationships between diseases by mapping composite vectors for symptoms or diseases to the most similar ICD codes. Second, we investigated the validity of ICD2Vec by comparing the biological relationships and cosine similarities among the vectorized ICD codes. Third, we propose a new risk score called IRIS, derived from ICD2Vec, and demonstrate its clinical utility with large cohorts from the UK and South Korea. RESULTS: Semantic compositionality was qualitatively confirmed between descriptions of symptoms and ICD2Vec. For example, the diseases most similar to COVID-19 were found to be the common cold (ICD-10: J00), unspecified viral hemorrhagic fever (ICD-10: A99), and smallpox (ICD-10: B03). We show the significant associations between the cosine similarities derived from ICD2Vec and the biological relationships using disease-to-disease pairs. Furthermore, we observed significant adjusted hazard ratios (HR) and area under the receiver operating characteristics (AUROC) between IRIS and risks for eight diseases. For instance, the higher IRIS for coronary artery disease (CAD) can be the higher probability for the incidence of CAD (HR: 2.15 [95% CI 2.02-2.28] and AUROC: 0.587 [95% CI 0.583-0.591]). We identified individuals at substantially increased risk of CAD using IRIS and 10-year atherosclerotic cardiovascular disease risk (adjusted HR: 4.26 [95% CI 3.59-5.05]). CONCLUSIONS: ICD2Vec, a proposed universal framework for converting qualitatively measured ICD codes into quantitative vectors containing semantic relationships between diseases, exhibited a significant correlation with actual biological significance. In addition, the IRIS was a significant predictor of major diseases in a prospective study using two large-scale datasets. Based on this clinical validity and utility evidence, we suggest that publicly available ICD2Vec can be used in diverse research and clinical practices and has important clinical implications.

Genetic network structure of 13 psychiatric disorders in the general population.

Psychiatric disorders frequently co-occur and share common symptoms and genetic backgrounds. Previous research has used genome-wide association studies to identify the interrelationships among psychiatric disorders and identify clusters of disorders; however, these methods have limitations in terms of their ability to examine the relationships among disorders as a network structure and their generalizability to the general population. In this study, we explored the network structure of the polygenic risk score (PRS) for 13 psychiatric disorders in a general population (276,249 participants of European ancestry from the UK Biobank) and identified communities and the centrality of the network. In this network, the nodes represented a PRS for each psychiatric disorder and the edges represented the connections between nodes. The psychiatric disorders comprised four robust communities. The first community included attention-deficit hyperactivity disorder, autism spectrum disorder, major depressive disorder, and anxiety disorder. The second community consisted of bipolar I and II disorders, schizophrenia, and anorexia nervosa. The third group included Tourette's syndrome and obsessive-compulsive disorder. Cannabis use disorder, alcohol use disorder, and post-traumatic stress disorder make up the fourth community. The PRS of schizophrenia had the highest values for the three metrics (strength, betweenness, and closeness) in the network. Our findings provide a comprehensive genetic network of psychiatric disorders and biological evidence for the classification of psychiatric disorders.

Anti-Inflammatory Actions of Soluble Ninjurin-1 Ameliorate Atherosclerosis.

BACKGROUND: Macrophages produce many inflammation-associated molecules, released by matrix metalloproteinases, such as adhesion molecules, and cytokines, as well, which play a crucial role in atherosclerosis. In this context, we investigated the relationship between Ninjurin-1 (Ninj1 [nerve injury-induced protein]), a novel matrix metalloproteinase 9 substrate, expression, and atherosclerosis progression. METHODS: Ninj1 expression and atherosclerosis progression were assessed in atherosclerotic aortic tissue and serum samples from patients with coronary artery disease and healthy controls, and atheroprone apolipoprotein e-deficient (Apoe-/-) and wild-type mice, as well. Apoe-/- mice lacking systemic Ninj1 expression (Ninj1-/-Apoe-/-) were generated to assess the functional effects of Ninj1. Bone marrow transplantation was also used to generate low-density lipoprotein receptor-deficient (Ldlr-/-) mice that lack Ninj1 specifically in bone marrow-derived cells. Mice were fed a Western diet for 5 to 23 weeks, and atherosclerotic lesions were investigated. The anti-inflammatory role of Ninj1 was verified by treating macrophages and mice with the peptides Ninj11-56 (ML56) and Ninj126-37 (PN12), which mimic the soluble form of Ninj1 (sNinj1). RESULTS: Our in vivo results conclusively showed a correlation between Ninj1 expression in aortic macrophages and the extent of human and mouse atherosclerotic lesions. Ninj1-deficient macrophages promoted proinflammatory gene expression by activating mitogen-activated protein kinase and inhibiting the phosphoinositide 3-kinase/Akt signaling pathway. Whole-body and bone marrow-specific Ninj1 deficiencies significantly increased monocyte recruitment and macrophage accumulation in atherosclerotic lesions through elevated macrophage-mediated inflammation. Macrophage Ninj1 was directly cleaved by matrix metalloproteinase 9 to generate a soluble form that exhibited antiatherosclerotic effects, as assessed in vitro and in vivo. Treatment with the sNinj1-mimetic peptides, ML56 and PN12, reduced proinflammatory gene expression in human and mouse classically activated macrophages, thereby attenuating monocyte transendothelial migration. Moreover, continuous administration of mPN12 alleviated atherosclerosis by inhibiting the enhanced monocyte recruitment and inflammation characteristics of this disorder in mice, regardless of the presence of Ninj1. CONCLUSIONS: Ninj1 is a novel matrix metalloproteinase 9 substrate in macrophages, and sNinj1 is a secreted atheroprotective protein that regulates macrophage inflammation and monocyte recruitment in atherosclerosis. Moreover, sNinj1-mediated anti-inflammatory effects are conserved in human macrophages and likely contribute to human atherosclerosis.

Renal Cell Carcinoma-Infiltrating CD3low Vγ9Vδ1 T Cells Represent Potentially Novel Anti-Tumor Immune Players.

Due to the highly immunogenic nature of renal cell carcinoma (RCC), the tumor microenvironment (TME) is enriched with various innate and adaptive immune subsets. In particular, gamma-delta (γδ) T cells can act as potent attractive mediators of adoptive cell transfer immunotherapy because of their unique properties such as non-reliance on major histocompatibility complex expression, their ability to infiltrate human tumors and recognize tumor antigens, relative insensitivity to immune checkpoint molecules, and broad tumor cytotoxicity. Therefore, it is now critical to better characterize human γδ T-cell subsets and their mechanisms in RCCs, especially the stage of differentiation. In this study, we aimed to identify γδ T cells that might have adaptive responses against RCC progression. We characterized γδ T cells in peripheral blood and tumor-infiltrating lymphocytes (TILs) in freshly resected tumor specimens from 20 RCC patients. Furthermore, we performed a gene set enrichment analysis on RNA-sequencing data from The Cancer Genome Atlas (TCGA) derived from normal kidneys and RCC tumors to ascertain the association between γδ T-cell infiltration and anti-cancer immune activity. Notably, RCC-infiltrating CD3low Vγ9Vδ1 T cells with a terminally differentiated effector memory phenotype with up-regulated activation/exhaustion molecules were newly detected as predominant TILs, and the cytotoxic activity of these cells against RCC was confirmed in vitro. In an additional analysis of the TCGA RCC dataset, γδ T-cell enrichment scores correlated strongly with those for CTLs, Th1 cells, "exhausted" T cells, and M1 macrophages, suggesting active involvement of γδ T cells in anti-tumor rather than pro-tumor activity, and Vδ1 cells were more abundant than Vδ2 or Vδ3 cells in RCC tumor samples. Thus, we posit that Vγ9Vδ1 T cells may represent an excellent candidate for adoptive immunotherapy in RCC patients with a high risk of relapse after surgery.

Linking transcriptional and genetic tumor heterogeneity through allele analysis of single-cell RNA-seq data.

Characterization of intratumoral heterogeneity is critical to cancer therapy, as the presence of phenotypically diverse cell populations commonly fuels relapse and resistance to treatment. Although genetic variation is a well-studied source of intratumoral heterogeneity, the functional impact of most genetic alterations remains unclear. Even less understood is the relative importance of other factors influencing heterogeneity, such as epigenetic state or tumor microenvironment. To investigate the relationship between genetic and transcriptional heterogeneity in a context of cancer progression, we devised a computational approach called HoneyBADGER to identify copy number variation and loss of heterozygosity in individual cells from single-cell RNA-sequencing data. By integrating allele and normalized expression information, HoneyBADGER is able to identify and infer the presence of subclone-specific alterations in individual cells and reconstruct the underlying subclonal architecture. By examining several tumor types, we show that HoneyBADGER is effective at identifying deletions, amplifications, and copy-neutral loss-of-heterozygosity events and is capable of robustly identifying subclonal focal alterations as small as 10 megabases. We further apply HoneyBADGER to analyze single cells from a progressive multiple myeloma patient to identify major genetic subclones that exhibit distinct transcriptional signatures relevant to cancer progression. Other prominent transcriptional subpopulations within these tumors did not line up with the genetic subclonal structure and were likely driven by alternative, nonclonal mechanisms. These results highlight the need for integrative analysis to understand the molecular and phenotypic heterogeneity in cancer.

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.

Dynamics of circulating tumor DNA during postoperative radiotherapy in patients with residual triple-negative breast cancer following neoadjuvant chemotherapy: a prospective observational study.

BACKGROUND: This study was performed to evaluate circulating tumor DNA (ctDNA) kinetics during postoperative radiotherapy (PORT) in patients with residual triple-negative breast cancer (TNBC) at surgery following neoadjuvant chemotherapy (NAC). METHODS: Stage II/III patients with post-NAC residual TNBC who required PORT were prospectively included in this study between March 2019 and July 2020. For 11 TNBC patients, next-generation sequencing targeting 38 genes was conducted in 55 samples, including tumor tissue, three plasma samples, and leukocytes from each patient. The plasma samples were collected at three-time points; pre-PORT (T0), after 3 weeks of PORT (T1), and 1 month after PORT (T2). Serial changes in ctDNA variant allele frequency (VAF) were analyzed. RESULTS: Somatic variants were found in the tumor specimens in 9 out of 11 (81.8%) patients. Mutated genes included TP53 (n = 7); PIK3CA (n = 2); and AKT1, APC, CSMD3, MYC, PTEN, and RB1 (n = 1). These tumor mutations were not found in plasma samples. Plasma ctDNA variants were detected in three (27.3%) patients at T0. Mutations in EGFR (n = 1), CTNNB1 (n = 1), and MAP2K (n = 1) was identified with ctDNA analysis. In two (18.2%) patients, the ctDNA VAF decreased through T1 and T2 while increasing at T2 in one (9.1%) patient. After a median follow-up of 22 months, no patient showed cancer recurrence. CONCLUSION: Among patients with post-NAC residual TNBC, more than a quarter exhibited a detectable amount of ctDNA after curative surgery. The ctDNA VAF changed variably during the course of PORT. Therefore, ctDNA kinetics can serve as a biomarker for optimizing adjuvant treatment.

Computational modeling of malignant ascites reveals CCL5-SDC4 interaction in the immune microenvironment of ovarian cancer.

Fluid accumulation in the abdominal cavity is commonly found in advanced-stage ovarian cancer patients, which creates a specialized tumor microenvironment for cancer progression. Using single-cell RNA sequencing (scRNA-seq) of ascites cells from five patients with ovarian cancer, we identified seven cell types, including heterogeneous macrophages and ovarian cancer cells. We resolved a distinct polarization state of macrophages by MacSpectrum analysis and observed subtype-specific enrichment of pathways associated with their functions. The communication between immune and cancer cells was predicted through a putative ligand-receptor pair analysis using NicheNet. We found that CCL5, a chemotactic ligand, is enriched in immune cells (T cells and NK cells) and mediates ovarian cancer cell survival in the ascites, possibly through SDC4. Moreover, SDC4 expression correlated with poor overall survival in ovarian cancer patients. Our study highlights the potential role of T cells and NK cells in long-term survival patients with ovarian cancer, indicating SDC4 as a potential prognostic marker in ovarian cancer patients.

The role of PDGFRA as a therapeutic target in young colorectal cancer patients.

BACKGROUND: Young patients with colorectal cancer (CRC) exhibit poor prognoses compared to older patients due to the difficulty in early diagnosis and treatment. However, the underlying molecular characteristics are still unclear. METHODS: We conducted a comprehensive analysis of 49 CRC patients without hereditary CRC using the whole-exome and RNA sequencing with tumor and matched normal samples. A total of 594 TCGA samples and 4 patient-derived cells were utilized for validation. RESULTS: Consensus molecular subtype 4 (CMS4) (53.85%) and CMS2 (38.46%) were enriched in the young (≤ 40 years) and old (> 60 years) age groups, respectively. A CMS4-associated gene, platelet-derived growth factor receptor α (PDGFRA), was significantly upregulated in young patients with CRC (FC = 3.21, p = 0.0001) and was negatively correlated with age (p = 0.0001, R = - 0.526). Moreover, PDGFRA showed a positive co-expression with metastasis-related genes in young CRC patients. In vitro validation confirmed that young patient-derived cells (PDCs) showed an enriched expression of PDGFRA compared to old PDCs and a reduced proliferation rate by knockdown of PDGFRA. Furthermore, young CRC patients were more sensitive to regorafenib, a PDGFRA-targeting drug, than old CRC patients. CONCLUSIONS: Our study suggests that CRC in young patients is associated with CMS4 and PDGFRA. In addition, PDGFRA may serve potential of novel therapeutic strategies and represent a predictive biomarker of response to regorafenib for young CRC patients.

Neutrophils expressing lysyl oxidase-like 4 protein are present in colorectal cancer liver metastases resistant to anti-angiogenic therapy.

Colorectal cancer liver metastases (CRCLM) that present with a replacement histopathological growth pattern (HGP) are resistant to neoadjuvant anti-angiogenic therapy. Surrogate biomarkers are not available to preoperatively identify patients with these tumors. Here we identify differentially expressed genes between CRCLM with a replacement HGP and those with a desmoplastic HGP using RNA sequencing. We demonstrate that LOXL4 is transcriptionally upregulated in replacement HGP CRCLM compared with desmoplastic HGP CRCLM and the adjacent normal liver. Interestingly, lysyl oxidase-like 4 (LOXL4) protein was expressed by neutrophils present in the tumor microenvironment in replacement HGP CRCLM. We further demonstrate that LOXL4 expression is higher in circulating neutrophils of cancer patients compared with healthy control patients and its expression can be induced by stimulation with lipopolysaccharide and TNF-α. Our study is the first to show the expression of LOXL4 in neutrophils and reveals the potential for LOXL4-expressing neutrophils to support the replacement HGP phenotype and to serve as a surrogate biomarker for this subtype of CRCLM. © 2020 Pathological Society of Great Britain and Ireland. Published by John Wiley & Sons, Ltd.

Single Cell Genomics for Tumor Heterogeneity.

Single cell genomics became a universal and powerful tool to study cellular diversity at genomic levels in normal and disease conditions. Cancer is a disease of genomic instability which instigates clonal evolution and intra-tumoral heterogeneity. Cancer progression also accompanies gross alterations in the microenvironment, and the stromal or immune cell types comprising the tumor microenvironment can be explored by single cell genomics. So far, breast cancer has been analyzed by single cell genomic tools for the clonal evolution, inter- and intra-tumoral heterogeneity in molecular signatures, and tumor microenvironment. We will briefly go over those studies and discuss the potential application of single cell genomics for the diagnostics and management of cancer.