APOBEC3-Mediated RNA Editing in Breast Cancer is Associated with Heightened Immune Activity and Improved Survival

The cytidine deaminase APOBEC3A regulates nucleolar function to promote cell growth and ribosome biogenesis

Cancer initiates as a consequence of genomic mutations and its subsequent progression relies in part on increased production of ribosomes to maintain high levels of protein synthesis for unchecked cell growth. Recently, cytidine deaminases have been uncovered as sources of mutagenesis in cancer. In an attempt to form a connection between these 2 cancer driving processes, we interrogated the cytidine deaminase family of proteins for potential roles in human ribosome biogenesis. We identified and validated APOBEC3A and APOBEC4 as novel ribosome biogenesis factors through our laboratory's established screening platform for the discovery of regulators of nucleolar function in MCF10A cells. Through siRNA depletion experiments, we highlight APOBEC3A's requirement in making ribosomes and specific role within the processing and maturation steps that form the large subunit 5.8S and 28S ribosomal (r)RNAs. We demonstrate that a subset of APOBEC3A resides within the nucleolus and associates with critical ribosome biogenesis factors. Mechanistic insight was revealed by transient overexpression of both wild-type and a catalytically dead mutated APOBEC3A, which both increase cell growth and protein synthesis. Through an innovative nuclear RNA sequencing methodology, we identify only modest predicted APOBEC3A C-to-U target sites on the pre-rRNA and pre-mRNAs. Our work reveals a potential direct role for APOBEC3A in ribosome biogenesis likely independent of its editing function. More broadly, we found an additional function of APOBEC3A in cancer pathology through its function in ribosome biogenesis, expanding its relevance as a target for cancer therapeutics.

Determinants of gastric cancer immune escape identified from non-coding immune-landscape quantitative trait loci

The landscape of non-coding mutations in cancer progression and immune evasion is largely unexplored. Here, we identify transcrptome-wide somatic and germline 3' untranslated region (3'-UTR) variants from 375 gastric cancer patients from The Cancer Genome Atlas. By performing gene expression quantitative trait loci (eQTL) and immune landscape QTL (ilQTL) analysis, we discover 3'-UTR variants with cis effects on expression and immune landscape phenotypes, such as immune cell infiltration and T cell receptor diversity. Using a massively parallel reporter assay, we distinguish between causal and correlative effects of 3'-UTR eQTLs in immune-related genes. Our approach identifies numerous 3'-UTR eQTLs and ilQTLs, providing a unique resource for the identification of immunotherapeutic targets and biomarkers. A prioritized ilQTL variant signature predicts response to immunotherapy better than standard-of-care PD-L1 expression in independent patient cohorts, showcasing the untapped potential of non-coding mutations in cancer.

The implications of APOBEC3-mediated C-to-U RNA editing for human disease

Intra-organism biodiversity is thought to arise from epigenetic modification of constituent genes and post-translational modifications of translated proteins. Here, we show that post-transcriptional modifications, like RNA editing, may also contribute. RNA editing enzymes APOBEC3A and APOBEC3G catalyze the deamination of cytosine to uracil. RNAsee (RNA site editing evaluation) is a computational tool developed to predict the cytosines edited by these enzymes. We find that 4.5% of non-synonymous DNA single nucleotide polymorphisms that result in cytosine to uracil changes in RNA are probable sites for APOBEC3A/G RNA editing; the variant proteins created by such polymorphisms may also result from transient RNA editing. These polymorphisms are associated with over 20% of Medical Subject Headings across ten categories of disease, including nutritional and metabolic, neoplastic, cardiovascular, and nervous system diseases. Because RNA editing is transient and not organism-wide, future work is necessary to confirm the extent and effects of such editing in humans.

SARS-CoV-2 and innate immunity: the good, the bad, and the "goldilocks"

An ancient conflict between hosts and pathogens has driven the innate and adaptive arms of immunity. Knowledge about this interplay can not only help us identify biological mechanisms but also reveal pathogen vulnerabilities that can be leveraged therapeutically. The humoral response to SARS-CoV-2 infection has been the focus of intense research, and the role of the innate immune system has received significantly less attention. Here, we review current knowledge of the innate immune response to SARS-CoV-2 infection and the various means SARS-CoV-2 employs to evade innate defense systems. We also consider the role of innate immunity in SARS-CoV-2 vaccines and in the phenomenon of long COVID.

Computational drug discovery of an inhibitor of APOBEC3B as a treatment for epithelial cancers

Cancer is one of the leading causes of death in the U.S., and tumorous cancers such as cervical, lung, breast, and ovarian cancers are the most common types. APOBEC3B is a nonessential cytidine deaminase found in humans and theorized to defend against viral infection. However, overexpression of APOBEC3B is linked to cancer in humans, which makes APOBEC3B a potential cancer treatment target through competitive inhibition for several tumorous cancers. Computational studies can help reveal a small molecule inhibitor using high-throughput virtual screening of millions of candidates with relatively little cost. This study aims to narrow the field of potential APOBEC3B inhibition candidates for future in vitro assays and provide an effective scaffold for drug design studies. Another goal of this project is to provide critical amino acid targets in the active site for future drug design studies. This study simulated 7.8 million drug candidates using high-throughput virtual screening and further processed the top scoring 241 molecules from AutoDock Vina, DOCK 6, and de novo design. Using virtual screening, de novo design, and molecular dynamics simulations, a competitive inhibitor candidate was discovered with an average binding free energy score of -46.03 kcal/mol, more than 10 kcal/mol better than the substrate control (dCMP). These results indicate that this molecule (or a structural derivative) may be an effective inhibitor of APOBEC3B and prevent host genome mutagenesis resulting from protein overexpression. Another important finding is the confirmation of essential amino acid targets, such as Tyr250 and Gln213 within the active site of APOBEC3B. Therefore, study used novel computational methods to provide a theoretical scaffold for future drug design studies that may prove useful as a treatment for epithelial cancers.Communicated by Ramaswamy H. Sarma.

Constructing an APOBEC-related gene signature with predictive value in the overall survival and therapeutic sensitivity in lung adenocarcinoma

Background

APOBEC family play an important role in cancer mutagenesis and tumor development. The role of APOBEC family in lung adenocarcinoma (LUAD) has not been studied comprehensively.

Materials and methods

The expression data of pan-cancer as well as LUAD was obtained from public databases. The expression level of APOBEC family genes was analyzed in different normal and cancer tissues. APOBEC mutagenesis enrichment score (AMES) was utilized to evaluate the APOBEC-induced mutations and the relation of APOBEC with genomic instability. Gene set enrichment analysis was used to identify differentially enriched pathways. Univariate Cox regression and Lasso regression were applied to screen key prognostic genes. The immune cell infiltration was estimated by CIBERSORT. RT-qPCR assay, CCK-8 and Transwell assay were conducted to explore gene expression and lung cancer cell invasion.

Results

Cancer tissues had significantly altered expression of APOBEC family genes and the expression patterns of APOBEC family were different in different cancer types. APOBEC3B was the most aberrantly expressed in most cancer types. In LUAD, we observed a significantly positive correlation of AMES with intratumor heterogeneity (ITH), tumor neoantigen burden (TNB), and tumor mutation burden (TMB). High AMES group had high mutation counts of DNA damage repair pathways, and high enrichment of cell cycle and DNA repair pathways. We identified four prognostic genes (LYPD3, ANLN, MUC5B, and FOSL1) based on AMES, and constructed an AMES-related gene signature. The expressions of four genes were enhanced and accelerated the invasion ability and viability of lung cancer cells. Furthermore, we found that high group increased oxidative stress level.

Conclusions

APOBEC family was associated with genomic instability, DNA damage-related pathways, and cell cycle in LUAD. The AMES-related gene signature had a great potential to indicate the prognosis and guide immunotherapy/chemotherapy for patients suffering from LUAD.

Genome-wide perturbations of A-to-I RNA editing dysregulated circular RNAs promoting the development of cervical cancer

Cervical cancer, the second most common female malignant tumor, seriously threatens women's health and lives. Despite the availability of the HPV vaccine, effective treatment options for cervical cancer are still lacking. New research perspectives now clarify that RNA editing dysregulation and changes in circRNA expression are jointly involved in disease pathogenesis, so molecular changes associated with circRNA and RNA editing may provide clues for the development of new therapeutic strategies for cervical cancer. In this study, we designed a series of pipelines to identify and analyze dysregulated RNA editing events in circRNAs. Our findings indicate a decrease in A-to-I RNA editing levels in cervical cancer compared to normal tissues, and editing may influence the back-splicing process of circRNAs through structural modifications of Alu elements. Moreover, our research reveals that RNA editing could modulate circRNA biogenesis by influencing RNA binding protein (RBP) binding on a transcriptome-wide scale, as well as influence the expression and coding potential of circRNAs. Importantly, we identified three RNA editing sites that could serve as potential biomarkers. In summary, our study presents a comprehensive landscape of RNA editing perturbations in circRNAs, providing new insights into the complex relationship between RNA editing and circRNA dysregulation in cervical cancer.

Development of a Machine Learning-Based Prognostic Model for Hormone Receptor-Positive Breast Cancer Using Nine-Gene Expression Signature

Background

Determining the prognosis of hormone receptor positive (HR+) breast cancer (BC), which accounts for 80% of all BCs, is critical in improving survival outcomes. Stratifying individuals at high risk of BC-related mortality and improving prognosis has been the focus of research for over a decade. However, these tools are not universal as they are limited to clinical factors. We hypothesized that a new framework for predicting prognosis in HR+ BC patients can develop using artificial intelligence.

Methods

A total of 2,338 HR+ human epidermal growth factor receptor 2 negative (HER2-) BC cases were analyzed from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC), The Cancer Genome Atlas (TCGA), and Gene Expression Omnibus (GEO) cohorts. Groups were then divided into high- and low-risk categories utilizing a recurrence prediction model (RPM). An RPM was created by extracting nine prognosis-related genes from over 18,000 genes using a logistic progression model.

Results

Risk classification by RPM was significantly stratified in both the discovery cohort and validation cohort. In the time-dependent area under the curve analysis, there was some variation depending on the cohort, but accuracy was found to decline significantly after about 10 years. Cell cycle related gene sets, MYC, and PI3K-AKT-mTOR signaling were enriched in high-risk tumors by the Gene Set Enrichment Analysis. High-risk tumors were associated with high levels of immune cells from the lymphoid and myeloid lineage and immune cytolytic activity, as well as low levels of stem cells and stromal cells. High-risk tumors were also associated with poor therapeutic effects of chemotherapy and endocrine therapy.

Conclusions

This model was able to stratify prognosis in multiple cohorts. This is because the model reflects major BC therapeutic target pathways and tumor immune microenvironment and, further is supported by the therapeutic effect of chemotherapy and endocrine therapy.

Unraveling the Enzyme-Substrate Properties for APOBEC3A-Mediated RNA Editing

The APOBEC3 family of human cytidine deaminases is involved in various cellular processes, including the innate and acquired immune system, mostly through inducing C-to-U in single-stranded DNA and/or RNA mutations. Although recent studies have examined RNA editing by APOBEC3A (A3A), its intracellular target specificity are not fully characterized. To address this gap, we performed in-depth analysis of cellular RNA editing using our recently developed sensitive cell-based fluorescence assay. Our findings demonstrate that A3A and an A3A-loop1-containing APOBEC3B (A3B) chimera are capable of RNA editing. We observed that A3A prefers to edit specific RNA substrates which are not efficiently deaminated by other APOBEC members. The editing efficiency of A3A is influenced by the RNA sequence contexts and distinct stem-loop secondary structures. Based on the identified RNA specificity features, we predicted potential A3A-editing targets in the encoding region of cellular mRNAs and discovered novel RNA transcripts that are extensively edited by A3A. Furthermore, we found a trend of increased synonymous mutations at the sites for more efficient A3A-editing, indicating evolutionary adaptation to the higher editing rate by A3A. Our results shed light on the intracellular RNA editing properties of A3A and provide insights into new RNA targets and potential impact of A3A-mediated RNA editing.

The Role of C-to-U RNA Editing in Human Biodiversity

Intra-organism biodiversity is thought to arise from epigenetic modification of our constituent genes and post-translational modifications after mRNA is translated into proteins. We have found that post-transcriptional modification, also known as RNA editing, is also responsible for a significant amount of our biodiversity, substantively expanding this story. The APOBEC (apolipoprotein B mRNA editing catalytic polypeptide-like) family RNA editing enzymes APOBEC3A and APOBEC3G catalyze the deamination of cytosines to uracils (C>U) in specific stem-loop structures.1,2 We used RNAsee (RNA site editing evaluation), a tool developed to predict the locations of APOBEC3A/G RNA editing sites, to determine whether known single nucleotide polymorphisms (SNPs) in DNA could be replicated in RNA via RNA editing. About 4.5% of non-synonymous SNPs which result in C>U changes in RNA, and about 5.4% of such SNPs labelled as pathogenic, were identified as probable sites for APOBEC3A/G editing. This suggests that the variant proteins created by these DNA mutations may also be created by transient RNA editing, with the potential to affect human health. Those SNPs identified as potential APOBEC3A/G-mediated RNA editing sites were disproportionately associated with cardiovascular diseases, digestive system diseases, and musculoskeletal diseases. Future work should focus on common sites of RNA editing, any variant proteins created by these RNA editing sites, and the effects of these variants on protein diversity and human health. Classically, our biodiversity is thought to come from our constitutive genetics, epigenetic phenomenon, transcriptional differences, and post-translational modification of proteins. Here, we have shown evidence that RNA editing, often stimulated by environmental factors, could account for a significant degree of the protein biodiversity leading to human disease. In an era where worries about our changing environment are ever increasing, from the warming of our climate to the emergence of new diseases to the infiltration of microplastics and pollutants into our bodies, understanding how environmentally sensitive mechanisms like RNA editing affect our own cells is essential.

Unveiling the interplay between mutational signatures and tumor microenvironment: a pan-cancer analysis

Background

While recent studies have separately explored mutational signatures and the tumor microenvironment (TME), there is limited research on the associations of both factors in a pan-cancer context.

Materials and methods

We performed a pan-cancer analysis of over 8,000 tumor samples from The Cancer Genome Atlas (TCGA) project. Machine learning methods were employed to systematically explore the relationship between mutational signatures and TME and develop a risk score based on TME-associated mutational signatures to predict patient survival outcomes. We also constructed an interaction model to explore how mutational signatures and TME interact and influence cancer prognosis.

Results

Our analysis revealed a varied association between mutational signatures and TME, with the Clock-like signature showing the most widespread influence. Risk scores based on mutational signatures mainly induced by Clock-like and AID/APOBEC activity exhibited strong pan-cancer survival stratification ability. We also propose a novel approach to predict transcriptome decomposed infiltration levels using genome-derived mutational signatures as an alternative approach for exploring TME cell types when transcriptome data are unavailable. Our comprehensive analysis revealed that certain mutational signatures and their interaction with immune cells significantly impact clinical outcomes in particular cancer types. For instance, T cell infiltration levels only served as a prognostic biomarker in melanoma patients with high ultraviolet radiation exposure, breast cancer patients with high homologous recombination deficiency signature, and lung adenocarcinoma patients with high tobacco-associated mutational signature.

Conclusion

Our study comprehensively explains the complex interplay between mutational signatures and immune infiltration in cancer. The results highlight the importance of considering both mutational signatures and immune phenotypes in cancer research and their significant implications for developing personalized cancer treatments and more effective immunotherapy.

Tumor immune microenvironment and immunotherapy efficacy in BRAF mutation non-small-cell lung cancer

Previous small-size studies reported BRAF-mutated NSCLC patients have comparable sensitivity to immune checkpoint inhibitors (ICIs). However, how BRAF mutation affects the tumor immune microenvironment (TIME) is unknown. We performed Nanostring-panel RNA sequencing to evaluate TIME in 57 BRAF mutated and wild-type (WT) NSCLC specimens (cohort A). The efficacy of ICI monotherapy or combined therapies was determined in 417 patients with WT and BRAF mutated NSCLC (cohort B). We found that BRAF-mutant tumors had similar ratios of CD8+ T cells to Tregs, the balance of cytotoxicity gene expression signatures and immune suppressive features, and similar ICI-response-related biomarkers to WT NSCLC. A similar TIME pattern was observed between the BRAF V600E and Non-V600E subgroups of NSCLC. The further retrospective study confirmed that treatment with ICI monotherapy or combined therapies resulted in similar overall survival (OS) (HR: 0.85; 95% CI, 0.56 to 1.30; p = 0.47) and progress-free survival (PFS) (HR: 1.02; 95% CI, 0.72 to 1.44; p = 0.91) of patients with WT (n = 358) and BRAF mutant (n = 59) NSCLC. Similarly, both patients with BRAF V600E or Non-V600E NSCLC had similar responses to immunotherapy. Our findings support that BRAF mutation did not modulate TIME in NSCLC and therapeutic responses to ICIs. Patients with NSCLC harboring BRAF mutation should not be denied treatment with ICIs.

Enhanced immune response outperform aggressive cancer biology and is associated with better survival in triple-negative breast cancer

Although the value of tumor-infiltrating lymphocytes is well known, the clinical relevance of an increased immune response, specifically in breast cancer, has not been investigated across large cohorts of patients using computational algorithms. Our hypothesis stated that an enhanced immune response is associated with an improvement in outcomes. To quantify the immune response, we utilized the allograft rejection score correlated with cytolytic activity and with all the other Hallmark immune-related gene sets. The score reflected the amount of infiltrating immune cells that correlated with the immune checkpoint molecule expressions, including CD4⁺ and CD8⁺ T cells, T helper type 1 (Th1) and type 2 (Th2) cells, M1 macrophages, B cells, and plasmacytoid dendritic cells (pDC). A high score was associated with high levels of intratumor heterogeneity, homologous recombination defects, mutation rate, histological grade, advanced stage, and lymph node metastasis. Breast malignancy with a high score enriched immune-related gene sets and pro-cancer-related gene sets, including epithelial–mesenchymal transition and KRAS pathway, in ER-positive/HER2-negative and triple-negative breast cancer (TNBC) groups. TNBC had the highest score compared to other subtypes, and was associated with better survival. In conclusion, we found that breast cancer with a high immune response is associated with aggressive cancer biology, but with better survival in TNBC.

Left sided breast cancer is associated with aggressive biology and worse outcomes than right sided breast cancer

Breast cancer is more common on the left side than the right side. We aim to evaluate differences in clinicopathological and genomic characteristics based on laterality. We analyzed survival outcomes and clinical characteristics of 881,320 patients recorded by the Surveillance, Epidemiology, and End Results (SEER) program. The Cancer Genome Atlas (TCGA) was used to explore genomic and clinical features from 1,062 patients. Gene expression data was used to quantitate cytolytic activity and hallmark gene-sets were used for gene set enrichment analysis. An institutional retrospective review was conducted on 155 patients treated with neoadjuvant chemotherapy (NACT). Patient characteristics were summarized by pathological complete response (pCR). Left sided tumors were found to be more prevalent than right sided tumors. No major clinicopathological differences were noted by laterality. Left sided breast cancer demonstrated poorer outcomes versus right sided tumors (HR 1.05, 95% CI 1.01–1.08; p = 0.011). Cell proliferation gene sets, including E2F Targets, G2M Checkpoint, Mitotic spindle, and MYC Targets, were enriched on the left side compared to the right. Left sided tumors had lower pCR rates versus right sided tumors (15.4% versus 29.9%, p = 0.036). Our findings suggest that left sided breast cancer is associated with aggressive biology and worse outcomes compared to right sided breast cancer.

Emerging Roles of Long Noncoding RNAs in Breast Cancer Epigenetics and Epitranscriptomics

Breast carcinogenesis is a multistep process that involves both genetic and epigenetic changes. Epigenetics refers to reversible changes in gene expression that are not accompanied by changes in gene sequence. In breast cancer (BC), dysregulated epigenetic changes, such as DNA methylation and histone modifications, are accompanied by epitranscriptomic changes, in particular adenine to inosine modifications within RNA molecules. Factors that trigger these phenomena are largely unknown, but there is evidence for widespread participation of long noncoding RNAs (lncRNAs) that already have been linked to virtually any aspect of BC biology, making them promising biomarkers and therapeutic targets in BC patients. Here, we provide a systematic review of known and possible roles of lncRNAs in epigenetic and epitranscriptomic processes, along with methods and tools to study them, followed by a brief overview of current challenges regarding the use of lncRNAs in medical applications.

Clinical Relevance of Estrogen Reactivity in the Breast Cancer Microenvironment

Purpose

Estrogen signals play an important role in the phenotype of estrogen receptor-positive breast cancer. However, comprehensive analyses of the effect of responsiveness to estrogen signals on the tumor microenvironment and survival in large cohorts of primary breast cancer patients have been lacking. We aimed to test the hypothesis that estrogen reactivity affects gene expression and immune cell infiltration profiles in the tumor microenvironment and survival.

Methods

A total of 3,098 breast cancer cases were analyzed: 1,904 from the Molecular Taxonomy of Breast Cancer (METABRIC) cohort, 1,082 from The Cancer Genome Atlas (TCGA) cohort, and 112 from the Hokkaido University Hospital cohort. We divided the group into estrogen reactivity-high and estrogen reactivity-low groups utilizing the scores of ESTROGEN_RESPONSE_EARLY and ESTROGEN_RESPONSE_LATE in Gene Set Variation Analysis.

Results

Breast cancer with high estrogen reactivity was related to Myc targets, metabolism-related signaling, cell stress response, TGF-beta signaling, androgen response, and MTORC1 signaling gene sets in the tumor microenvironment. Low estrogen reactivity was related to immune-related proteins, IL2-STAT5 signaling, IL6-JAK-STAT3 signaling, KRAS signaling, cell cycle-related gene sets, and EMT. In addition, breast cancer with high levels of estrogen reactivity had low immune cytolytic activity and low levels of immunostimulatory cells. It also had low levels of stimulatory and inhibitory factors of the cancer immunity cycle. Patients with high estrogen reactivity were also associated with a better prognosis.

Conclusion

We demonstrated the relationship between estrogen reactivity and the profiles of immune cells and gene expression, as well as survival.

Gastric cancer with enhanced apical junction pathway has increased metastatic potential and worse clinical outcomes

Excessive intercellular connection at confluency may be limiting further cell growth or a sign of aggressive biology in the cell culture. As apical junction complex is a main component of cell-to-cell connection, we aimed to investigate gastric cancer biology using Apical Junction Pathway score that we generated using Gene set variant analysis (GSVA) of the "Hallmark Apical Junction" gene set. 1,239 gastric cancer patients from the Cancer Genome Atlas (TCGA) and two GSE cohorts were included in this study. The cohorts were dichotomized using the median of the score. Apical Junction Pathway score high gastric cancer was not consistently associated with increased cell proliferation or immune cell infiltration. On the other hand, Apical Junction Pathway score high gastric cancer was associated with significantly higher infiltration of stromal cells, such as endothelial cells; hence, increased neovascularization and angiogenesis in the tumor microenvironment (TME) were speculated. Gene set enrichment analysis (GSEA) confirmed increased expression of epithelial mesenchymal transition (EMT) and angiogenesis in the high Apical Junction Pathway score group (false discovery rate (FDR) <0.25). Lastly, the high Apical Junction Pathway score group was associated with more aggressive clinicopathological characteristics, such as significantly higher American Joint Committee on Cancer (AJCC) T-category and higher pathological stage, leading to worse disease-specific survival and overall survival (P<0.05, respectively). In conclusion, enhanced Apical Junction Pathway score gastric cancer was associated with aggressive clinical characteristics leading to shorter survival likely due to increased metastatic potential from EMT and angiogenesis.

Autoimmunity and Cancer-Two Sides of the Same Coin

Autoimmune disease results from the immune response against self-antigens, while cancer develops when the immune system does not respond to malignant cells. Thus, for years, autoimmunity and cancer have been considered as two separate fields of research that do not have a lot in common. However, the discovery of immune checkpoints and the development of anti-cancer drugs targeting PD-1 (programmed cell death receptor 1) and CTLA-4 (cytotoxic T lymphocyte antigen 4) pathways proved that studying autoimmune diseases can be extremely helpful in the development of novel anti-cancer drugs. Therefore, autoimmunity and cancer seem to be just two sides of the same coin. In the current review, we broadly discuss how various regulatory cell populations, effector molecules, genetic predisposition, and environmental factors contribute to the loss of self-tolerance in autoimmunity or tolerance induction to cancer. With the current paper, we also aim to convince the readers that the pathways involved in cancer and autoimmune disease development consist of similar molecular players working in opposite directions. Therefore, a deep understanding of the two sides of immune tolerance is crucial for the proper designing of novel and selective immunotherapies.

Functions and consequences of AID/APOBEC-mediated DNA and RNA deamination

The AID/APOBEC polynucleotide cytidine deaminases have historically been classified as either DNA mutators or RNA editors based on their first identified nucleic acid substrate preference. DNA mutators can generate functional diversity at antibody genes but also cause genomic instability in cancer. RNA editors can generate informational diversity in the transcriptome of innate immune cells, and of cancer cells. Members of both classes can act as antiviral restriction factors. Recent structural work has illuminated differences and similarities between AID/APOBEC enzymes that can catalyse DNA mutation, RNA editing or both, suggesting that the strict functional classification of members of this family should be reconsidered. As many of these enzymes have been employed for targeted genome (or transcriptome) editing, a more holistic understanding will help improve the design of therapeutically relevant programmable base editors.

In this Perspective, Pecori et al. provide an overview of the AID/APOBEC cytidine deaminase family, discussing key structural features, how they contribute to viral and tumour evolution and how they can be harnessed for (potentially therapeutic) base-editing purposes.

Immune Checkpoint Inhibitors in Triple Negative Breast Cancer: The Search for the Optimal Biomarker

Triple negative breast cancer (TNBC) is a high-risk and aggressive malignancy characterized by the absence of estrogen receptors (ER) and progesterone receptors (PR) on the surface of malignant cells, and by the lack of overexpression of human epidermal growth factor 2 (HER2). It has limited therapeutic options compared to other subtypes of breast cancer. There is now a growing body of evidence on the role of immunotherapy in TNBC, however much of the data from clinical trials is conflicting and thus, challenging for clinicians to integrate the data into clinical practice. Landmark phase III trials using immunotherapy in the early-stage neoadjuvant setting concluded that the addition of immunotherapy to chemotherapy improved the pathologic complete response (pCR) rate compared to chemotherapy with placebo while others found no significant improvement in pCR. Phase III trials have investigated the utility of immunotherapy in previously untreated metastatic TNBC, and these studies have similarly arrived at inconsistent conclusions. Some studies showed no benefit while others demonstrated a clinically significant improvement in overall survival in the PD-L1 positive population. It is not yet clear which biomarkers are most useful, and assays for these biomarkers have not been standardized. Given the often serious and severe side effects of immunotherapy, it is important and necessary to identify predictive biomarkers of response and resistance in order to enhance patient selection. In this review, we will discuss both the challenges of traditional biomarkers and the opportunities of emerging biomarkers for patient selection.

APOBEC3F expression in triple-negative breast cancer is associated with tumor microenvironment infiltration and activation of cancer immunity and improved survival

The apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like (APOBEC) causes a point mutation from cytidine to uracil in DNA and/or RNA. The role of APOBEC3A and APOBEC3B in breast cancer has been well described, whereas that of APOBEC3F remains unknown. To investigate the clinical relevance of APOBEC3F expression, we analyzed a total of 3000 breast cancer cases from multiple independent large patient cohorts including METABRIC, TCGA, GSE75688, and GSE114725. High expression of APOBEC3F was associated with improved disease-specific and overall survival in triple negative breast cancer (TNBC). APOBEC3F is not usually a reflection of cancer cell biology in TNBC or luminal breast cancer, except for homologous recombination deficiency in TNBC. In the TNBC homologous recombination deficiency group, APOBEC3F expression was not consistently associated with intratumor heterogeneity, mutation rates, or neoantigens. APOBEC3F expression did not correlate with response to any of the drugs tested in breast cancer cell lines in vitro. However, high APOBEC3F expression was associated with enrichment of several immune-related gene sets and immune activity. High APOBEC3F expression also accompanied higher infiltration of anti-cancer immune cell infiltration in TNBC. However, in luminal breast cancer, high APOBEC3F tumor significantly enriched not only immune-related gene sets, but also cell proliferation-, metastasis-, and apoptosis-related gene sets. Analysis of single-cell transcriptomes showed APOBEC3F exclusively expressed in immune cells and significantly associated with cytolytic activity of the immune cells, immune response, and immune cell proliferation. Expression of immune checkpoint genes was uniformly elevated in APOBEC3F-high tumors. We conclude that APOBEC3F is exclusively expressed in immune cells and this expression is associated with enhanced anti-cancer immune response as well as improved survival in TNBC.

Mutagenic Activity of AID/APOBEC Deaminases in Antiviral Defense and Carcinogenesis

Proteins of the AID/APOBEC family are capable of cytidine deamination in nucleic acids forming uracil. These enzymes are involved in mRNA editing, protection against viruses, the introduction of point mutations into DNA during somatic hypermutation, and antibody isotype switching. Since these deaminases, especially AID, are potent mutagens, their expression, activity, and specificity are regulated by several intracellular mechanisms. In this review, we discuss the mechanisms of impaired expression and activation of AID/APOBEC proteins in human tumors and their role in carcinogenesis and tumor progression. Also, the diagnostic and potential therapeutic value of increased expression of AID/APOBEC in different types of tumors is analyzed. We assume that in the case of solid tumors, increased expression of endogenous deaminases can serve as a marker of response to immunotherapy since multiple point mutations in host DNA could lead to amino acid substitutions in tumor proteins and thereby increase the frequency of neoepitopes.

Different members of the APOBEC3 family of DNA mutators have opposing associations with the landscape of breast cancer

APOBEC enzymes are strong mutagenic factors. In breast cancer, expression of APOBEC3B is increased and associated with mutation load and poor outcome. Other APOBEC3s can also mutate DNA but their clinical significance in breast cancer and its underpinnings have not been comprehensively studied. In our examination of 1,091 breast carcinoma cases, high expression of APOBEC3A or APOBEC3B genes was associated with greater tumor burden of mutations and other genomic aberrations. Expression of none of the five APOBEC3C-H genes had any correlation with these features, including T[C-T/G]W mutations, but their high expression levels indicated a robust anti-cancer immune response within tumors, with elevated CD8⁺ T cell abundance, T cell receptor diversity, and immune cytolytic activity. Concordantly, survival analyses of this and two other cohorts with > 3,000 patients each showed favorable prognostic benefit of high APOBEC3C-H expression for both cancer progression and mortality. A detrimental prognostic value was observed for APOBEC3A and APOBEC3B. Single-cell data revealed cancer epithelial and stromal immune cells as major sources of APOBEC3B and APOBEC3C-H expression in tumors, respectively. These observations on opposing associations with breast cancer of different APOBEC3s highlight the contrasting roles of these enzymes, promoting cancer through mutagenesis while antagonizing it through immune response.

Regulatory roles of RNA modifications in breast cancer

Collectively referred to as the epitranscriptome, RNA modifications play important roles in gene expression control regulating relevant cellular processes. In the last few decades, growing numbers of RNA modifications have been identified not only in abundant ribosomal (rRNA) and transfer RNA (tRNA) but also in messenger RNA (mRNA). In addition, many writers, erasers and readers that dynamically regulate the chemical marks have also been characterized. Correct deposition of RNA modifications is prerequisite for cellular homeostasis, and its alteration results in aberrant transcriptional programs that dictate human disease, including breast cancer, the most frequent female malignancy, and the leading cause of cancer-related death in women. In this review, we emphasize the major RNA modifications that are present in tRNA, rRNA and mRNA. We have categorized breast cancer-associated chemical marks and summarize their contribution to breast tumorigenesis. In addition, we describe less abundant tRNA modifications with related pathways implicated in breast cancer. Finally, we discuss current limitations and perspectives on epitranscriptomics for use in therapeutic strategies against breast and other cancers.

Mutational spectrum of SARS-CoV-2 during the global pandemic

Viruses accumulate mutations under the influence of natural selection and host-virus interactions. Through a systematic comparison of 351,525 full viral genome sequences collected during the recent COVID-19 pandemic, we reveal the spectrum of SARS-CoV-2 mutations. Unlike those of other viruses, the mutational spectrum of SARS-CoV-2 exhibits extreme asymmetry, with a much higher rate of C>U than U>C substitutions, as well as a higher rate of G>U than U>G substitutions. This suggests directional genome sequence evolution during transmission. The substantial asymmetry and directionality of the mutational spectrum enable pseudotemporal tracing of SARS-CoV-2 without prior information about the root sequence, collection time, and sampling region. This shows that the viral genome sequences collected in Asia are similar to the original genome sequence. Adjusted estimation of the dN/dS ratio accounting for the asymmetrical mutational spectrum also shows evidence of negative selection on viral genes, consistent with previous reports. Our findings provide deep insights into the mutational processes in SARS-CoV-2 viral infection and advance the understanding of the history and future evolution of the virus.

Octogenarians' Breast Cancer Is Associated with an Unfavorable Tumor Immune Microenvironment and Worse Disease-Free Survival

Simple Summary

In recent years, as the elderly population has grown, the number of elderly breast cancer patients has increased, but their biological characteristics are still controversial. This study investigated octogenarians’ breast cancer biology and its tumor microenvironment utilizing an in-silico translational approach to multiple large patient cohorts. We found that octogenarians’ breast cancer was associated with worse survival and an unfavorable tumor immune microenvironment such as M2 macrophage but not with aggressive cancer cell biology. Our report is important for understanding the characteristics of elderly breast cancer patients and would be critical for the development of breast cancer treatment in the future.

Abstract

Elderly patients are known to have a worse prognosis for breast cancer. This is commonly blamed on their medical comorbidities and access to care. However, in addition to these social issues, we hypothesized that the extreme elderly (octogenarians—patients over 80 years old) have biologically worse cancer with unfavorable tumor immune microenvironment. The Cancer Genomic Atlas (TCGA) and the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC) breast cancer cohorts were analyzed. The control (aged 40–65) and octogenarians numbered 668 and 53 in TCGA and 979 and 118 in METABRIC, respectively. Octogenarians had significantly worse breast cancer-specific survival in both cohorts (p < 0.01). Octogenarians had a higher ER-positive subtype rate than controls in both cohorts. Regarding PAM50 classification, luminal-A and -B subtypes were significantly higher in octogenarians, whereas basal and claudin-low subtypes were significantly lower (p < 0.05) in octogenarians. There was no difference in tumor mutation load, intratumor heterogeneity, or cytolytic activity by age. However, the octogenarian cohort was significantly associated with high infiltration of pro-cancer immune cells, M2 macrophage, and regulatory T cells in both cohorts (p < 0.05). Our results demonstrate that octogenarians’ breast cancer is associated with worse survival and with an unfavorable tumor immune microenvironment.

Enhanced Thermogenesis in Triple-Negative Breast Cancer Is Associated with Pro-Tumor Immune Microenvironment

Mild cold stress induced by housing mice with a 4T1 triple-negative breast cancer (TNBC) cell implantation model at 22 °C increases tumor growth rate with a pro-tumorigenic immune microenvironment (lower CD8 ⁺T cells, higher myeloid-derived suppressor cells (MDSCs) and regulatory T-cells (Tregs)). Since cold stress also activates thermogenesis, we hypothesized that enhanced thermogenesis is associated with more aggressive cancer biology and unfavorable tumor microenvironment (TME) in TNBC patients. A total of 6479 breast cancer patients from METABRIC, TCGA, GSE96058, GSE20194, and GSE25066 cohorts were analyzed using Kyoto Encyclopedia of Genes and Genomes (KEGG) thermogenesis score. High-thermogenesis TNBC was associated with a trend towards worse survival and with angiogenesis, adipogenesis, and fatty acid metabolism pathways. On the other hand, low-thermogenesis TNBC enriched most of the hallmark cell-proliferation-related gene sets (i.e., mitotic spindle, E2F targets, G2M checkpoint, MYC targets), as well as immune-related gene sets (i.e., IFN-α and IFN-γ response). Favorable cytotoxic T-cell-attracting chemokines CCL5, CXCL9, CXCL10, and CXCL11 were lower; while the MDSC- and Treg-attracting chemokine CXCL12 was higher. There were higher M2 but lower M1 macrophages and Tregs. In conclusion, high-thermogenesis TNBC is associated with pro-tumor immune microenvironment and may serve as biomarker for testing strategies to overcome this immunosuppression.

A mark of disease: how mRNA modifications shape genetic and acquired pathologies

RNA modifications have recently emerged as a widespread and complex facet of gene expression regulation. Counting more than 170 distinct chemical modifications with far-reaching implications for RNA fate, they are collectively referred to as the epitranscriptome. These modifications can occur in all RNA species, including messenger RNAs (mRNAs) and noncoding RNAs (ncRNAs). In mRNAs the deposition, removal, and recognition of chemical marks by writers, erasers and readers influence their structure, localization, stability, and translation. In turn, this modulates key molecular and cellular processes such as RNA metabolism, cell cycle, apoptosis, and others. Unsurprisingly, given their relevance for cellular and organismal functions, alterations of epitranscriptomic marks have been observed in a broad range of human diseases, including cancer, neurological and metabolic disorders. Here, we will review the major types of mRNA modifications and editing processes in conjunction with the enzymes involved in their metabolism and describe their impact on human diseases. We present the current knowledge in an updated catalog. We will also discuss the emerging evidence on the crosstalk of epitranscriptomic marks and what this interplay could imply for the dynamics of mRNA modifications. Understanding how this complex regulatory layer can affect the course of human pathologies will ultimately lead to its exploitation toward novel epitranscriptomic therapeutic strategies.

The effects of RNA editing in cancer tissue at different stages in carcinogenesis

RNA editing is one of the most prevalent and abundant forms of post-transcriptional RNA modification observed in normal physiological processes and often aberrant in diseases including cancer. RNA editing changes the sequences of mRNAs, making them different from the source DNA sequence. Edited mRNAs can produce editing-recoded protein isoforms that are functionally different from the corresponding genome-encoded protein isoforms. The major type of RNA editing in mammals occurs by enzymatic deamination of adenosine to inosine (A-to-I) within double-stranded RNAs (dsRNAs) or hairpins in pre-mRNA transcripts. Enzymes that catalyse these processes belong to the adenosine deaminase acting on RNA (ADAR) family. The vast majority of knowledge on the RNA editing landscape relevant to human disease has been acquired using in vitro cancer cell culture models. The limitation of such in vitro models, however, is that the physiological or disease relevance of results obtained is not necessarily obvious. In this review we focus on discussing in vivo occurring RNA editing events that have been identified in human cancer tissue using samples surgically resected or clinically retrieved from patients. We discuss how RNA editing events occurring in tumours in vivo can identify pathological signalling mechanisms relevant to human cancer physiology which is linked to the different stages of cancer progression including initiation, promotion, survival, proliferation, immune escape and metastasis.

Genome editing technologies: CRISPR, LEAPER, RESTORE, ARCUT, SATI, and RESCUE

Genome editing technologies include techniques used for desired genetic modifications and allow the insertion, modification or deletion of specific DNA fragments. Recent advances in genome biology offer unprecedented promise for interdisciplinary collaboration and applications in gene editing. New genome editing technologies enable specific and efficient genome modifications. The sources that inspire these modifications and already exist in the genome are DNA degradation enzymes and DNA repair pathways. Six of these recent technologies are the clustered regularly interspaced short palindromic repeats (CRISPR), leveraging endogenous ADAR for programmable editing of RNA (LEAPER), recruiting endogenous ADAR to specific transcripts for oligonucleotide-mediated RNA editing (RESTORE), chemistry-based artificial restriction DNA cutter (ARCUT), single homology arm donor mediated intron-targeting integration (SATI), RNA editing for specific C-to-U exchange (RESCUE). These technologies are widely used from various biomedical researches to clinics, agriculture, and allow you to rearrange genomic sequences, create cell lines and animal models to solve human diseases. This review emphasizes the characteristics, superiority, limitations, also whether each technology can be used in different biological systems and the potential application of these systems in the treatment of several human diseases.

Th2 cell infiltrations predict neoadjuvant chemotherapy response of estrogen receptor-positive breast cancer

Background

High infiltration of Th2 is linked to breast cancer progression and metastasis through the induction of cytokine release and T-cell anergy. The estrogen receptor (ER)-positive subtype, which accounts for 70% of breast cancer, is known to respond less to neoadjuvant chemotherapy (NAC) due to its low potential for proliferation. We hypothesized that Th2 high tumors are highly proliferative, and thus more likely to respond to NAC in ER-positive breast cancer.

Methods

We obtained clinicopathological data and overall survival information on 1,069 breast cancer patients from The Cancer Genome Atlas (TCGA). Computational algorithms and CIBERSORT were used to estimate immune cell infiltration. Additionally, xCell was used for validation.

Results

Th2 high tumors did not consistently associate with an unfavorable immune cell composition and tumor immune microenvironment but were found to be significantly elevated in the cancer stage. Th2 high tumors also correlated with high Nottingham pathological grade, as well as with Ki-67 and proliferation score in ER-positive subtypes. High Th2 tumors achieved a pathological complete response (pCR) significantly higher in ER-positive breast cancer.

Conclusions

In conclusion, high levels of Th2 are associated with aggressive features of breast cancer. Th2 levels may be a biomarker in patient selection for NAC in ER-positive breast cancer.

The Impact of Immunofunctional Phenotyping on the Malfunction of the Cancer Immunity Cycle in Breast Cancer

Simple Summary

The cancer-immunity cycle (CIC) is a series of self-sustaining stepwise events to fight cancer growth by the immune system. We hypothesized that immunofunctional phenotyping that represent the malfunction of the CIC is clinically relevant in breast cancer (BC) utilizing total of 2979 BC cases; 1075 from TCGA cohort, 1904 from METABRIC cohort were analyzed. The immunofunctional phenotype was classified as follows: hot T-cell infiltrated, high immune cytolytic activity (CYT), cold T-cell infiltrated, high frequency of CD8+ T cells and low CYT, and non-inflamed, low frequency of CD8+ T cells and low CYT. We demonstrated that immunofunctional phenotyping not only indicated the degree of anti-cancer immune dysfunction, but also served as a prognostic biomarker and HTI was inversely related to estrogen response.

Abstract

The cancer-immunity cycle (CIC) is a series of self-sustaining stepwise events to fight cancer growth by the immune system. We hypothesized that immunofunctional phenotyping that represent the malfunction of the CIC is clinically relevant in breast cancer (BC). Total of 2979 BC cases; 1075 from TCGA cohort, 1904 from METABRIC cohort were analyzed. The immunofunctional phenotype was classified as follows: hot T-cell infiltrated (HTI), high immune cytolytic activity (CYT), Cold T-cell infiltrated (CTI), high frequency of CD8+ T cells and low CYT, and non-inflamed, low frequency of CD8+ T cells and low CYT. The analysis of tumor immune microenvironment in the immunofunctional phenotype revealed that not only immunostimulatory factors, but also immunosuppressive factors were significantly elevated and immunosuppressive cells were significantly decreased in HTI. Patients in HTI were significantly associated with better survival in whole cohort and patients in CTI were significantly associated with worse survival in triple negative. Furthers, HTI was inversely related to estrogen responsive signaling. We demonstrated that immunofunctional phenotype not only indicated the degree of anti-cancer immune dysfunction, but also served as a prognostic biomarker and HTI was inversely related to estrogen response.

High Expression of NRF2 Is Associated with Increased Tumor-Infiltrating Lymphocytes and Cancer Immunity in ER-Positive/HER2-Negative Breast Cancer

Simple Summary

The clinical relevance of Nuclear factor erythroid 2-Related Factor 2 (NRF2) in human breast cancer remains unclear. A total of 5443 breast cancer patients with transcriptomic profile were analyzed for the clinical relevance of NRF2 expression, including cancer aggressiveness, immune cell infiltration, patient survival, and drug response. We found that tumors with high NRF2 expression were associated with better survival in ER-positive/HER2-negative breast cancer. NRF2 expression was equivalent in immune, stromal, and cancer cells in tumor microenvironment. We found that high NRF2 expression was associated with enhanced tumor-infiltrating lymphocytes in ER-positive/HER2-negative breast cancer. NRF2 expression significantly correlated with drug sensitivity in multiple ER-positive breast cancer cell lines, but not associated with pathological complete response after neoadjuvant chemotherapy in breast cancer patients regardless of subtypes.

Abstract

Nuclear factor erythroid 2-related factor 2 (NRF2) is a key modifier in breast cancer. It is unclear whether NRF2 suppresses or promotes breast cancer progression. We studied the clinical relevance of NRF2 expression by conducting in silico analyses in 5443 breast cancer patients from several large patient cohorts (METABRIC, GSE96058, GSE25066, GSE20194, and GSE75688). NRF2 expression was significantly associated with better survival, low Nottingham pathological grade, and ER-positive/HER2-negative and triple negative breast cancer (TNBC). High NRF2 ER-positive/HER2-negative breast cancer enriched inflammation- and immune-related gene sets by GSEA. NRF2 expression was elevated in immune, stromal, and cancer cells. High NRF2 tumors were associated with high infiltration of immune cells (CD8⁺, CD4⁺, and dendritic cells (DC)) and stromal cells (adipocyte, fibroblasts, and keratinocytes), and with low fraction of Th1 cells. NRF2 expression significantly correlated with area under the curve (AUC) of several drug response in multiple ER-positive breast cancer cell lines, however, there was no significant association between NRF2 and pathologic complete response (pCR) rate after neoadjuvant chemotherapy in human samples. Finally, high NRF2 breast cancer was associated with high expression of immune checkpoint molecules. In conclusion, NRF2 expression was associated with enhanced tumor-infiltrating lymphocytes in ER-positive/HER2-negative breast cancer.

Deaminase-Independent Mode of Antiretroviral Action in Human and Mouse APOBEC3 Proteins

Apolipoprotein B mRNA editing enzyme, catalytic polypeptide-like 3 (APOBEC3) proteins (APOBEC3s) are deaminases that convert cytosines to uracils predominantly on a single-stranded DNA, and function as intrinsic restriction factors in the innate immune system to suppress replication of viruses (including retroviruses) and movement of retrotransposons. Enzymatic activity is supposed to be essential for the APOBEC3 antiviral function. However, it is not the only way that APOBEC3s exert their biological function. Since the discovery of human APOBEC3G as a restriction factor for HIV-1, the deaminase-independent mode of action has been observed. At present, it is apparent that both the deaminase-dependent and -independent pathways are tightly involved not only in combating viruses but also in human tumorigenesis. Although the deaminase-dependent pathway has been extensively characterized so far, understanding of the deaminase-independent pathway remains immature. Here, we review existing knowledge regarding the deaminase-independent antiretroviral functions of APOBEC3s and their molecular mechanisms. We also discuss the possible unidentified molecular mechanism for the deaminase-independent antiretroviral function mediated by mouse APOBEC3.

Inflammation Is Associated with Worse Outcome in the Whole Cohort but with Better Outcome in Triple-Negative Subtype of Breast Cancer Patients

Inflammation has been linked with cancer, but whether it is part of the problem or part of the solution remains to be a matter of debate in breast cancer. Our group and others have demonstrated that inflammation aggravates cancer progression; however, some claim that inflammation may support immune cell infiltration and suppress cancer. We defined the gene set variation analysis of the Molecular Signatures Database Hallmark inflammatory response gene set as the inflammatory pathway score and analyzed 3632 tumors in total from 4 breast cancer cohorts (METABRIC, TCGA, GSE25066, and GSE21094). In the whole breast cancer cohort, high-score tumors were associated with aggressive clinical characteristics, such as worse disease specific survival, higher Nottingham histological grade, and younger age. Inflammatory score was significantly higher in triple-negative (TNBC) as well as basal and normal subtypes compared with the other subtypes, which suggest that the detrimental effect of high level of inflammation may be because it includes a more aggressive subtype. On the contrary, high score within TNBC was significantly associated with better survival. TNBC with high score enriched not only IFN-α, IFN-γ response, IL-2/STAT5 signaling, Allograft rejection, Complement, p53 pathway, Reactive Oxygen, and Apoptosis but also TNF-α signaling, IL6-JAK-STAT signaling, TGF-β signaling, Coagulation, Angiogenesis, EMT, KRAS signaling, and PI3K-AKT-MTOR signaling gene sets. High score was associated with mainly favorable anticancerous immune cell infiltration as well as Leukocyte fraction, TIL regional fraction, Lymphocyte infiltration, IFN-γ response, TGF-β response, and cytolytic activity scores. Although the inflammatory pathway score was not associated with neoadjuvant treatment response, it associated with expressions of immune checkpoint molecules. In conclusion, inflammation was associated with worse outcome in the whole breast cancer cohort, but with better outcome in TNBC, which was associated with favorable anticancerous immune response and immune cell infiltrations.

The interesting relationship between APOBEC3 deoxycytidine deaminases and cancer: a long road ahead

Cancer is considered a group of diseases characterized by uncontrolled growth and spread of abnormal cells and is propelled by somatic mutations. Apolipoprotein B mRNA-editing enzyme catalytic polypeptide-like 3 (APOBEC3) family of enzymes are endogenous sources of somatic mutations found in multiple human cancers. While these enzymes normally act as an intrinsic immune defence against viruses, they can also catalyse 'off-target' cytidine deamination in genomic single-stranded DNA intermediates. The deamination of cytosine forms uracil, which is promutagenic in DNA. Key factors to trigger the APOBEC 'off-target' activity are overexpression in a non-normal cell type, nuclear localization and replication stress. The resulting uracil-induced mutations contribute to genomic variation, which may result in neutral, beneficial or harmful consequences for the cancer. This review summarizes the functional and biochemical basis of the APOBEC3 enzyme activity and highlights their relationship with the most well-studied cancers in this particular context such as breast, lung, bladder, and human papillomavirus-associated cancers. We focus on APOBEC3A, APOBEC3B and APOBEC3H haplotype I because they are the leading candidates as sources of somatic mutations in these and other cancers. Also, we discuss the prognostic value of the APOBEC3 expression in drug resistance and response to therapies.

Plasmacytoid Dendritic Cell (pDC) Infiltration Correlate with Tumor Infiltrating Lymphocytes, Cancer Immunity, and Better Survival in Triple Negative Breast Cancer (TNBC) More Strongly than Conventional Dendritic Cell (cDC)

Simple Summary

Dendritic cells (DC) represent a major antigen-presenting cell in the tumor immune microenvironment (TIME) and play an essential role in cancer immunity. Clinical relevance of plasmacytoid DC (pDC) and conventional DC (cDC) infiltration were studied in large patient cohorts using transcriptomic analyses. We found that high pDC was significantly associated with better survival in triple negative breast cancer (TNBC), but not cDC. High pDC TNBC tumors enriched immune and inflammation-related gene sets including IFN-γ signaling more strongly than cDC TNBC, which also enriched metastasis-related gene sets. pDC correlated with CD8⁺ and CD4⁺ memory T cells, and had a cytolytic activity score stronger than cDC in TNBC. High pDC TNBC were associated with a high fraction of anti-cancer immune cells and high expression of all the immune checkpoint molecules examined. For the first time, we found that pDC are correlated with immune response and survival in TNBC patients more strongly than cDC.

Abstract

Dendritic cells (DC) represent a major antigen-presenting cell type in the tumor immune microenvironment (TIME) and play an essential role in cancer immunity. Conventional DC (cDC) and plasmacytoid DC (pDC) were defined by the xCell algorithm and a total of 2968 breast cancer patients (TCGA and METABRIC) were analyzed. We found that triple-negative breast cancer (TNBC) had a high fraction of cDC and pDC compared to the other subtypes. In contrast to cDC, high pDC in TNBC was significantly associated with better disease-specific and disease-free survival consistently in both cohorts. High cDC TNBC tumors enriched not only inflammation and immune-related, but also metastasis-related gene sets in Gene Set Enrichment Analysis, whereas high pDC TNBC enriched inflammation and immune -related gene sets including IFN-γ signaling more strongly than cDC. pDC TNBC correlated with CD8⁺, CD4⁺ memory, IFN-γ score, and cytolytic activity stronger than cDC TNBC. High pDC TNBC were associated with a high fraction of anti-cancer immune cells and high expression of all the immune check point molecules examined. In conclusion, pDC levels correlated with the infiltration of immune cells and patient survival in TNBC more strongly than cDC; this is the first study suggesting the clinical relevance of pDC infiltration in TNBC.

MYC Targets Scores Are Associated with Cancer Aggressiveness and Poor Survival in ER-Positive Primary and Metastatic Breast Cancer

MYC is one of the most studied oncogenes that is known to promote cell proliferation. We utilized MYC targets v1 and MYC targets v2 scores of gene set variation analysis and hypothesized that these scores correlate with tumor aggressiveness and survival outcomes. We examined a total of 3109 breast cancer patients from TCGA, METABRIC, and GSE124647 cohorts. In each cohort, the patients were divided into high- and low-score groups using the upper third value as the cut off. As expected, higher scores were related to increased cell proliferation and worse clinical and pathologic features. High MYC targets scores were associated with worse survival, specifically in primary ER-positive breast cancer, consistently in both TCGA and METABRIC cohorts. In ER-positive breast cancer, high MYC targets v1, but not v2 score, was associated with high mutation load, and high MYC targets v1 and v2 scores were both associated with increased infiltration of pro- and anti-cancerous immune cells. We found that high MYC scores were associated with worse survival in metastatic breast cancer. Our findings show that the MYC targets v1 and v2 scores are associated with tumor aggressiveness and poor prognosis in ER-positive primary tumors, as well as in metastatic breast cancer.

Abundance of Regulatory T Cell (Treg) as a Predictive Biomarker for Neoadjuvant Chemotherapy in Triple-Negative Breast Cancer

Simple Summary

Regulatory CD4⁺ T cell (Treg) is one of the suppressive immune cells, but data on its clinical relevance in large human breast cancer cohort is limited. Abundance of Tregs in 5177 breast cancer patient samples from five independent cohorts was analyzed by the xCell algorithm using tumor transcriptomics. Treg abundance was lower in metastatic tumors when compared to matched primary tumors. Treg was associated with a high mutation rate of TP53 genes and copy number mutations as well as with increased tumor infiltration of M2 macrophages and decreased infiltration of T helper type 1 cells. Interestingly, low Treg abundance was significantly associated with pathological complete response (pCR) after neoadjuvant chemotherapy (NAC) in TNBC, but not in ER-positive/Her2-negative subtype. Abundance of Treg was also associated with high expression of multiple immune checkpoint molecules. In conclusion, Treg abundance may have a potential as a predictive biomarker of pCR after NAC in TNBC.

Abstract

Regulatory CD4⁺ T cell (Treg), a subset of tumor-infiltrating lymphocytes (TILs), are known to suppress anticancer immunity but its clinical relevance in human breast cancer remains unclear. In this study, we estimated the relative abundance of Tregs in breast cancer of multiple patient cohorts by using the xCell algorithm on bulk tumor gene expression data. In total, 5177 breast cancer patients from five independent cohorts (TCGA-BRCA, GSE96058, GSE25066, GSE20194, and GSE110590) were analyzed. Treg abundance was not associated with cancer aggressiveness, patient survival, or immune activity markers, but it was lower in metastatic tumors when compared to matched primary tumors. Treg was associated with a high mutation rate of TP53 genes and copy number mutations as well as with increased tumor infiltration of M2 macrophages and decreased infiltration of T helper type 1 (Th1) cells. Pathological complete response (pCR) after neoadjuvant chemotherapy (NAC) was significantly associated with low Treg abundance in triple negative breast cancer (TNBC) but not in ER-positive/Her2-negative subtype. High Treg abundance was significantly associated with high tumor expression of multiple immune checkpoint inhibitor genes. In conclusion, Treg abundance may have potential as a predictive biomarker of pCR after NAC in TNBC.

M1 Macrophage and M1/M2 ratio defined by transcriptomic signatures resemble only part of their conventional clinical characteristics in breast cancer

Tumor associated macrophages (TAMs) play a critical role in biology of various cancers, including breast cancer. In the current study, we defined “M1” macrophage and “M1”/“M2” ratio by transcriptomic signatures using xCell. We investigated the association between high level of “M1” macrophage or “M1”/“M2” ratio and the tumor immune microenvironment by analyzing the transcriptome of publicly available cohorts, TCGA and METABRIC. We found that “M1” high tumors were not associated with prolonged survival compared with “M1” low tumors, or with the response to neoadjuvant chemotherapy. “M1” high tumors were associated with clinically aggressive features and “M1” high tumors enriched the cell proliferation and cell cycle related gene sets in GSEA. At the same time, “M1” high tumors were associated with high immune activity and favorable tumor immune microenvironment, as well as high expression of immune check point molecules. Strikingly, all these results were mirrored in “M1”/“M2” ratio high tumors. In conclusion, transcriptomically defined “M1” or “M1”/“M2” high tumors were associated with aggressive cancer biology and favorable tumor immune microenvironment but not with survival benefit, which resembled only part of their conventional clinical characteristics.

ITPKC as a Prognostic and Predictive Biomarker of Neoadjuvant Chemotherapy for Triple Negative Breast Cancer

Triple negative breast cancer (TNBC) is the most aggressive subtype of breast cancer with higher mortality than the others. Pathological complete response (pCR) to neoadjuvant chemotherapy (NAC) is considered as a surrogate to predict survival. Inositol 1,4,5-trisphosphate 3-kinase C (ITPKC) is a negative regulator of T cell activation, and reduction in ITPKC function is known to promote Kawasaki disease. Given the role of tumor infiltrating lymphocytes in NAC and since TNBC has the most abundant immune cell infiltration in breast cancer, we hypothesized that the ITPKC expression level is associated with NAC response and prognosis in TNBC. The ITPKC gene was expressed in the mammary gland, but its expression was highest in breast cancer cells among other stromal cells in a bulk tumor. ITPKC expression was highest in TNBC, associated with its survival, and was its independent prognostic factor. Although high ITPKC was not associated with immune function nor with any immune cell fraction, low ITPKC significantly enriched cell proliferation-related gene sets in TNBC. TNBC with low ITPKC achieved a significantly higher pCR rate after NAC. To the best of our knowledge, this is the first report to demonstrate that ITPKC gene expression may be useful as a prognostic and predictive biomarker in TNBC.

CD8 T Cell Score as a Prognostic Biomarker for Triple Negative Breast Cancer

CD8 T cell is an essential component of tumor-infiltrating lymphocytes (TIL) and tumor immune microenvironment (TIME). Using the xCell CD8 T cell score of whole tumor gene expression data, we estimated these cells in total of 3837 breast cancer patients from TCGA, METABRIC and various GEO cohorts. The CD8 score correlated strongly with expression of CD8 genes. The score was highest for triple-negative breast cancer (TNBC), and a high score was associated with high tumor immune cytolytic activity and better survival in TNBC but not other breast cancer subtypes. In TNBC, tumors with a high CD8 score had enriched expression of interferon (IFN)-α and IFN-γ response and allograft rejection gene sets, and greater infiltration of anti-cancerous immune cells. The score strongly correlated with CD4 memory T cells in TNBC, and tumors with both a high CD8 score and high CD4 memory T cell abundance had significantly better survival. Finally, a high CD8 score was significantly associated with high expression of multiple immune checkpoint molecules. In conclusion, a high CD8 T cell score is associated with better survival in TNBC, particularly when tumor CD4 memory T cells were elevated. Our findings also suggest a possible use of the score as a predictive biomarker for response to immune checkpoint therapy.

The E2F Pathway Score as a Predictive Biomarker of Response to Neoadjuvant Therapy in ER+/HER2- Breast Cancer

E2F transcription factors play critical roles in the cell cycle. Therefore, their activity is expected to reflect tumor aggressiveness and responsiveness to therapy. We scored 3905 tumors of nine breast cancer cohorts for this activity based on their gene expression for the Hallmark E2F targets gene set. As expected, tumors with a high score had an increased expression of cell proliferation-related genes. A high score was significantly associated with shorter patient survival, greater MKI67 expression, histological grade, stage, and genomic aberrations. Furthermore, metastatic tumors had higher E2F scores than the primary tumors from which they arose. Although tumors with a high score had greater infiltration by both pro- and anti-cancerous immune cells, they had an increased expression of immune checkpoint genes. Estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative cancer with a high E2F score achieved a significantly higher pathological complete response (pCR) rate to neoadjuvant chemotherapy. The E2F score was significantly associated with the expression of cyclin-dependent kinase (CDK)-related genes and strongly correlated with sensitivity to CDK inhibition in cell lines. In conclusion, the E2F score is a marker of breast cancer aggressiveness and predicts the responsiveness of ER-positive/HER2-negative patients to neoadjuvant chemotherapy and possibly to CDK and immune checkpoint inhibitors.

Non-Coding RNA Editing in Cancer Pathogenesis

In the last two decades, RNA post-transcriptional modifications, including RNA editing, have been the subject of increasing interest among the scientific community. The efforts of the Human Genome Project combined with the development of new sequencing technologies and dedicated bioinformatic approaches created to detect and profile RNA transcripts have served to further our understanding of RNA editing. Investigators have determined that non-coding RNA (ncRNA) A-to-I editing is often deregulated in cancer. This discovery has led to an increased number of published studies in the field. However, the eventual clinical application for these findings remains a work in progress. In this review, we provide an overview of the ncRNA editing phenomenon in cancer. We discuss the bioinformatic strategies for RNA editing detection as well as the potential roles for ncRNA A to I editing in tumor immunity and as clinical biomarkers.

Contribution of Immune Cells to Glucocorticoid Receptor Expression in Breast Cancer

Breast cancer (BC) patients experience increased stress with elevated cortisol levels, increasing risk of cancer recurrence. Cortisol binds to a cytoplasmic receptor, glucocorticoid receptor (GR) encoded by GR gene (NR3C1). We hypothesized that not only cancer cells, but even immune cells in the tumor microenvironment (TME) may contribute to GR expression in bulk tumor and influence prognosis. To test this, mRNA expression data was accessed from METABRIC and TCGA. "High" and "low" expression was based on highest and lowest quartiles of NR3C1 gene expression, respectively. Single-cell sequencing data were obtained from GSE75688 and GSE114725 cohorts. Computer algorithms CIBERSORT, Gene Set Enrichment Analysis and TIMER were used. GR-high BC has better median disease-free and disease-specific survival. Single cell sequencing data showed higher GR expression on immune cells compared to cancer and stromal cells. Positive correlation between GR-high BC and CD8⁺ T-cells was noted. In GR-high tumors, higher cytolytic activity (CYT) with decreased T-regulatory and T-follicular helper cells was observed. High GR expression was associated with lower proliferation index Ki67, enriched in IL-2_STAT5, apoptosis, KRAS, TGF-β signaling, and epithelial-to-mesenchymal transition. Immune cells significantly contribute to GR expression of bulk BC. GR-high BC has a favorable TME with higher CYT with favorable outcomes.

A Novel 4-Gene Score to Predict Survival, Distant Metastasis and Response to Neoadjuvant Therapy in Breast Cancer

We generated a 4-gene score with genes upregulated in LM2-4, a metastatic variant of MDA-MB-231 (DOK 4, HCCS, PGF, and SHCBP1) that was strongly associated with disease-free survival (DFS) in TCGA cohort (hazard ratio [HR]>1.2, p < 0.02). The 4-gene score correlated with overall survival of TCGA (HR = 1.44, p < 0.001), which was validated with DFS and disease-specific survival of METABRIC cohort. The 4-gene score was able to predict worse survival or clinically aggressive tumors, such as high Nottingham pathological grade and advanced cancer staging. High score was associated with worse survival in the hormonal receptor (HR)-positive/Her2-negative subtype. High score enriched cell proliferation-related gene sets in GSEA. The score was high in primary tumors that originated, in and metastasized to, brain and lung, and it predicted worse progression-free survival for metastatic tumors. Good tumor response to neoadjuvant chemotherapy or hormonal therapy was accompanied by score reduction. High scores were also predictive of response to neoadjuvant chemotherapy for HR-positive/Her2-negative subtype. High score tumors had increased expression of T cell exhaustion marker genes, suggesting that the score may also be a biomarker for immunotherapy response. Our novel 4-gene score with both prognostic and predictive values may, therefore, be clinically useful particularly in HR-positive breast cancer.

G2M Cell Cycle Pathway Score as a Prognostic Biomarker of Metastasis in Estrogen Receptor (ER)-Positive Breast Cancer

The vast majority of breast cancer death is a result of metastasis. Thus, accurate identification of patients who are likely to have metastasis is expected to improve survival. The G2M checkpoint plays a critical role in cell cycle. We hypothesized that breast cancer tumors with high activity of G2M pathway genes are more aggressive and likely to metastasize. To test this, we used the single-sample gene set variation analysis method to calculate the score for the Hallmark G2M checkpoint pathway using gene expression data of a total of 4626 samples from 12 human breast cancer cohorts. As expected, a high G2M pathway score correlated with enriched tumor expression of other cell proliferation-related gene sets. The score was significantly associated with clinical aggressive features of tumors and patient survival in estrogen receptor (ER)-positive/human epidermal growth factor receptor 2 (HER2)-negative breast cancer. Interestingly, a high G2M score of metastasis tumors was also significantly associated with worse survival. In primary as well as metastasis tumors with high scores, the infiltration of both pro- and anti-cancerous immune cells increased. Tumor G2M score was also associated with treatment response to systemic chemotherapy in ER-positive/HER2-negative cancer, and was predictive of response to cyclin-dependent kinase inhibition therapy.

Clinical Relevance of RNA Editing to Early Detection of Cancer in Human

DNA encodes RNA and is responsible for protein production in cells. RNA editing is the process by which genetic information is altered in the RNA molecule. RNA editing in cancer initiation, progression and development has been well documented and play an important role in tumorigenesis. Studying RNA editing and its application to change genetic information after transcription, RNA-editing technology could be an important innovation in cancer and has the potential for more effective precision treatment. Bioengineering integration approach and artificial intelligence could revolutionize the entire field of RNA editing for early detection of cancer.

KRAS signaling enriched triple negative breast cancer is associated with favorable tumor immune microenvironment and better survival

KRAS signaling is associated with cancer progression in several cancers. Upregulation of KRAS signaling is often seen in cancers that harbor high KRAS mutation rate, such as pancreatic cancer and non-small cell lung cancer (NSCLC). Less than 2% of breast cancers have KRAS mutation, however, the alteration of the effector signaling such as PI3K/AKT and MAPK pathways are well known. Mutated KRAS is known to function as immune suppressor in other cancers, but the role of KRAS signaling on tumor immune microenvironment (TIME) in breast cancer is not known. We hypothesize that the enrichment of KRAS signaling is associated with reduced patient survival as well as TIME in triple negative breast cancer (TNBC). Patient cohorts from Molecular Taxonomy of Breast Cancer International Consortium (METABRIC; n = 1903) and The Cancer Genome Atlas (TCGA; n = 982) were used. Higher expression of KRAS in breast cancer cell-lines (MCF7, BT474, and MDA-MB231) compared to MCF10A, which is a model of benign mammary cells was found. Both MEK and PI3K inhibitors suppressed MB231 cell proliferation in dose dependent manner. Gene Set Variant Analysis (GSVA) of the patient cohorts demonstrated two peaks by KRAS_SIGNALING_UP gene sets which were divided into KRAS-high and -low groups using median cutoff. There was no difference in KRAS mutation between KRAS-high and low. Despite its cell proliferation promoting role, KRAS-high patients demonstrated significantly better Disease-Free Survival and Overall Survival in triple negative breast cancer (TNBC). KRAS-high TNBC was associated with favorable tumor immune microenvironment with elevated B cells and CD8 T cells, monocytes, or M1 macrophage. It was associated with decreased CD4 central memory T-cells, but not Regulatory T-cells, or M2 macrophage detected by xCell. To elucidate the mechanism of this association, Gene Set Enrichment Analysis was performed. Inflammatory response, IL6/JAK-STAT3 signaling, and Interferon gamma response gene sets were enriched in KRAS-high TNBC patients in both METABRIC and TCGA cohorts. In agreement, cytolytic activity score, interferon gamma response score, and lymphocyte infiltrating signature score, were all significantly elevated in KRAS-high TNBC. In conclusion, we found that patients with enrichment of KRAS signaling gene sets were associated with inflammation and favorable tumor immune microenvironment as well as improved survival in TNBC.

Clinical Applications of RNA Editing Technology for the Early Detection of Cancer and Future Directions

Early detection of cancer has great clinical importance and potentially improves cure, survival rate and treatment outcome. RNA editing technology can be used as targeted and precise molecular scissors to cut and replace disease-causing genes with healthy ones. This is a post transcriptional modification that can lead to the recoding of proteins. RNA editing technology is in its infancy, but it can be used for early diagnoses and effective treatment of cancer. The full potential of precision medicine will be achieved by using the knowledge of RNA reversible-recoding to edit the protein. RNA editing technology could be used to expose chemo resistant cancer cells, dormant cancer stem cells and other malignant tumors. RNA editing generates RNA and protein diversity to accelerate and enhance the screening window for early detection of cancer. We propose that the RNA editing sites could be used as a novel tool for early detection of cancer.