Identification of intracellular bacteria from multiple single-cell RNA-seq platforms using CSI-Microbes.

The study of the tumor microbiome has been garnering increased attention. We developed a computational pipeline (CSI-Microbes) for identifying microbial reads from single-cell RNA sequencing (scRNA-seq) data and for analyzing differential abundance of taxa. Using a series of controlled experiments and analyses, we performed the first systematic evaluation of the efficacy of recovering microbial unique molecular identifiers by multiple scRNA-seq technologies, which identified the newer 10x chemistries (3' v3 and 5') as the best suited approach. We analyzed patient esophageal and colorectal carcinomas and found that reads from distinct genera tend to co-occur in the same host cells, testifying to possible intracellular polymicrobial interactions. Microbial reads are disproportionately abundant within myeloid cells that up-regulate proinflammatory cytokines like IL1Β and CXCL8, while infected tumor cells up-regulate antigen processing and presentation pathways. These results show that myeloid cells with bacteria engulfed are a major source of bacterial RNA within the tumor microenvironment (TME) and may inflame the TME and influence immunotherapy response.

Distinct functions of wild-type and R273H mutant Δ133p53α differentially regulate glioblastoma aggressiveness and therapy-induced senescence.

Despite being mutated in 92% of TP53 mutant cancers, how mutations on p53 isoforms affect their activities remain largely unknown. Therefore, exploring the effect of mutations on p53 isoforms activities is a critical, albeit unexplored area in the p53 field. In this article, we report for the first time a mutant Δ133p53α-specific pathway which increases IL4I1 and IDO1 expression and activates AHR, a tumor-promoting mechanism. Accordingly, while WT Δ133p53α reduces apoptosis to promote DNA repair, mutant R273H also reduces apoptosis but fails to maintain genomic stability, increasing the risks of accumulation of mutations and tumor's deriving towards a more aggressive phenotype. Furthermore, using 2D and 3D spheroids culture, we show that WT Δ133p53α reduces cell proliferation, EMT, and invasion, while the mutant Δ133p53α R273H enhances all three processes, confirming its oncogenic potential and strongly suggesting a similar in vivo activity. Importantly, the effects on cell growth and invasion are independent of mutant full-length p53α, indicating that these activities are actively carried by mutant Δ133p53α R273H. Furthermore, both WT and mutant Δ133p53α reduce cellular senescence in a senescence inducer-dependent manner (temozolomide or radiation) because they regulate different senescence-associated target genes. Hence, WT Δ133p53α rescues temozolomide-induced but not radiation-induced senescence, while mutant Δ133p53α R273H rescues radiation-induced but not temozolomide-induced senescence. Lastly, we determined that IL4I1, IDO1, and AHR are significantly higher in GBMs compared to low-grade gliomas. Importantly, high expression of all three genes in LGG and IL4I1 in GBM is significantly associated with poorer patients' survival, confirming the clinical relevance of this pathway in glioblastomas. These data show that, compared to WT Δ133p53α, R273H mutation reorientates its activities toward carcinogenesis and activates the oncogenic IL4I1/IDO1/AHR pathway, a potential prognostic marker and therapeutic target in GBM by combining drugs specifically modulating Δ133p53α expression and IDO1/Il4I1/AHR inhibitors.

Urinary Metabolite Diagnostic and Prognostic Liquid Biopsy Biomarkers of Lung Cancer in Non-smokers and Tobacco Smokers.

PURPOSE: Non-smokers account for 10-13% of all lung cancer cases in the United States. Etiology is attributed to multiple risk factors including exposure to secondhand smoking, asbestos, environmental pollution, and radon, but these exposures are not within the current eligibility criteria for early lung cancer screening by low-dose computed tomography (LDCT). EXPERIMENTAL DESIGN: Urine samples were collected from two independent cohorts comprising 846 participants (exploratory cohort) and 505 participants (validation cohort). The cancer urinary biomarkers, creatine riboside (CR) and N-acetylneuraminic acid (NANA) were analyzed and quantified using liquid chromatography-mass spectrometry to determine if non-smoker cases can be distinguished from sex and age-matched controls in comparison to tobacco smoker cases and controls, potentially leading to more precise eligibility criteria for LDCT screening. RESULTS: Urinary levels of CR and NANA were significantly higher and comparable in non-smokers and tobacco smoker cases as compared to population controls in both cohorts. Receiver Operating Characteristics (ROC) analysis for combined CR and NANA levels in non-smokers of the exploratory cohort resulted in better predictive performance with the area under the curve (AUC) of 0.94, whereas the validation cohort non-smokers had an AUC of 0.80. Kaplan-Meier survival curves showed that high levels of CR and NANA were associated with increased cancer-specific death in non-smokers as well as tobacco smoker cases in both cohorts. CONCLUSIONS: Measuring CR and NANA in urine liquid biopsies could identify non-smokers at high risk for lung cancer as candidates for LDCT screening and warrant prospective studies of these biomarkers.

The landscape of etiological patterns of hepatocellular carcinoma and intrahepatic cholangiocarcinoma in Thailand.

Thailand is among countries with the highest global incidence and mortality rates of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA). While viral hepatitis and liver fluke infections have been associated with HCC and iCCA, respectively, other environmental risk factors, overall risk factor commonality and combinatorial roles, and effects on survival have not been systematically examined. We conducted a TIGER-LC consortium-based population study covering all high-incidence areas of both malignancies across Thailand: 837 HCC, 1474 iCCA, and 1112 controls (2011-2019) were comprehensively queried on lifelong environmental exposures, lifestyle, and medical history. Multivariate logistic regression and Cox proportional hazards analyses were used to evaluate risk factors and associated survival patterns. Our models identified shared risk factors between HCC and iCCA, such as viral hepatitis infection, liver fluke infection, and diabetes, including novel and shared associations of agricultural pesticide exposure (OR range of 1.50; 95% CI: 1.06-2.11 to 2.91; 95% CI: 1.82-4.63) along with vulnerable sources of drinking water. Most patients had multiple risk factors, magnifying their risk considerably. Patients with lower risk levels had better survival in both HCC (HR 0.78; 95% CI: 0.64-0.96) and iCCA (HR 0.84; 95% CI: 0.70-0.99). Risk factor co-exposures and their common associations with HCC and iCCA in Thailand emphasize the importance for future prevention and control measures, especially in its large agricultural sector. The observed mortality patterns suggest ways to stratify patients for anticipated survivorship and develop plans to support medical care of longer-term survivors, including behavioral changes to reduce exposures.

Acidovorax temperans skews neutrophil maturation and polarizes Th17 cells to promote lung adenocarcinoma development.

Change within the intratumoral microbiome is a common feature in lung and other cancers and may influence inflammation and immunity in the tumor microenvironment, affecting growth and metastases. We previously characterized the lung cancer microbiome in patients and identified Acidovorax temperans as enriched in tumors. Here, we instilled A. temperans in an animal model driven by mutant K-ras and Tp53. This revealed A. temperans accelerates tumor development and burden through infiltration of proinflammatory cells. Neutrophils exposed to A. temperans displayed a mature, pro-tumorigenic phenotype with increased cytokine signaling, with a global shift away from IL-1ß signaling. Neutrophil to monocyte and macrophage signaling upregulated MHC II to activate CD4+ T cells, polarizing them to an IL-17A+ phenotype detectable in CD4+ and γδ populations (T17). These T17 cells shared a common gene expression program predictive of poor survival in human LUAD. These data indicate bacterial exposure promotes tumor growth by modulating inflammation.

Enhancing chimeric antigen receptor T cell therapy by modulating the p53 signaling network with Δ133p53α.

Chimeric antigen receptor (CAR) T cell dysfunction is a major barrier to achieving lasting remission in hematologic cancers, especially in chronic lymphocytic leukemia (CLL). We have shown previously that Δ133p53α, an endogenous isoform of the human TP53 gene, decreases in expression with age in human T cells, and that reconstitution of Δ133p53α in poorly functional T cells can rescue proliferation [A. M. Mondal et al., J. Clin. Invest. 123, 5247-5257 (2013)]. Although Δ133p53α lacks a transactivation domain, it can form heterooligomers with full-length p53 and modulate the p53-mediated stress response [I. Horikawa et al., Cell Death Differ. 24, 1017-1028 (2017)]. Here, we show that constitutive expression of Δ133p53α potentiates the anti-tumor activity of CD19-directed CAR T cells and limits dysfunction under conditions of high tumor burden and metabolic stress. We demonstrate that Δ133p53α-expressing CAR T cells exhibit a robust metabolic phenotype, maintaining the ability to execute effector functions and continue proliferating under nutrient-limiting conditions, in part due to upregulation of critical biosynthetic processes and improved mitochondrial function. Importantly, we show that our strategy to constitutively express Δ133p53α improves the anti-tumor efficacy of CAR T cells generated from CLL patients that previously failed CAR T cell therapy. More broadly, our results point to the potential role of the p53-mediated stress response in limiting the prolonged antitumor functions required for complete tumor clearance in patients with high disease burden, suggesting that modulation of the p53 signaling network with Δ133p53α may represent a translationally viable strategy for improving CAR T cell therapy.

Pan-viral serology uncovers distinct virome patterns as risk predictors of hepatocellular carcinoma and intrahepatic cholangiocarcinoma.

This study evaluates the pan-serological profiles of hepatocellular carcinoma (HCC) and intrahepatic cholangiocarcinoma (iCCA) compared to several diseased and non-diseased control populations to identify risk factors and biomarkers of liver cancer. We used phage immunoprecipitation sequencing, an anti-viral antibody screening method using a synthetic-phage-displayed human virome epitope library, to screen patient serum samples for exposure to over 1,280 strains of pathogenic and non-pathogenic viruses. Using machine learning methods to develop an HCC or iCCA viral score, we discovered that both viral scores were positively associated with several liver function markers in two separate at-risk populations independent of viral hepatitis status. The HCC score predicted all-cause mortality over 8 years in patients with chronic liver disease at risk of HCC, while the viral hepatitis status was not predictive of survival. These results suggest that non-hepatitis viral infections may contribute to HCC and iCCA development and could be biomarkers in at-risk populations.

Gut dysbiosis in Thai intrahepatic cholangiocarcinoma and hepatocellular carcinoma.

Primary liver cancer (PLC), which includes intrahepatic cholangiocarcinoma (iCCA) and hepatocellular carcinoma (HCC), has the highest incidence of all cancer types in Thailand. Known etiological factors, such as viral hepatitis and chronic liver disease do not fully account for the country's unusually high incidence. However, the gut-liver axis, which contributes to carcinogenesis and disease progression, is influenced by the gut microbiome. To investigate this relationship, fecal matter from 44 Thai PLC patients and 76 healthy controls were subjected to whole-genome metagenomic shotgun sequencing and then analyzed by marker gene-based and assembly based methods. Results revealed greater gut microbiome heterogeneity in iCCA compared to HCC and healthy controls. Two Veillonella species were found to be more abundant in iCCA samples and could distinguish iCCA from HCC and healthy controls. Conversely, Ruminococcus gnavus was depleted in iCCA patients and could distinguish HCC from iCCA samples. High Veillonella genus counts in the iCCA group were associated with enriched amino acid biosynthesis and glycolysis pathways, while enriched phospholipid and thiamine metabolism pathways characterized the HCC group with high Blautia genus counts. These findings reveal distinct landscapes of gut dysbiosis among Thai iCCA and HCC patients and warrant further investigation as potential biomarkers.

Pilot study of plasma creatine riboside as a potential biomarker for cervical cancer.

This pilot study aimed primarily to evaluate plasma levels of a novel metabolite, creatine riboside, in patients with cervical cancer (discovery and validation cohorts, n = 11 for each) compared with non-cancer subjects (controls, n = 30). We found that the pre-treatment plasma creatine riboside level was significantly higher in the discovery cohort than in controls. The cut-off value determined from the discovery cohort distinguished 90.9% of the patients in the validation cohort from controls. Unbiased principal component analysis of plasma metabolites in high-creatine riboside samples demonstrated enrichment of pathways involved in arginine and creatine metabolism. These data indicate the potential utility of plasma creatine riboside as a biomarker of cervical cancer.

Matrin3 regulates mitotic spindle dynamics by controlling alternative splicing of CDC14B.

Matrin3 is an RNA-binding protein that regulates diverse RNA-related processes, including mRNA splicing. Although Matrin3 has been intensively studied in neurodegenerative diseases, its function in cancer remains unclear. Here, we report Matrin3-mediated regulation of mitotic spindle dynamics in colorectal cancer (CRC) cells. We comprehensively identified RNAs bound and regulated by Matrin3 in CRC cells and focused on CDC14B, one of the top Matrin3 targets. Matrin3 knockdown results in increased inclusion of an exon containing a premature termination codon in the CDC14B transcript and simultaneous down-regulation of the standard CDC14B transcript. Knockdown of CDC14B phenocopies the defects in mitotic spindle dynamics upon Matrin3 knockdown, and the elongated and misoriented mitotic spindle observed upon Matrin3 knockdown are rescued upon overexpression of CDC14B, suggesting that CDC14B is a key downstream effector of Matrin3. Collectively, these data reveal a role for the Matrin3/CDC14B axis in control of mitotic spindle dynamics.

Bile salt hydrolase in non-enterotoxigenic Bacteroides potentiates colorectal cancer.

Bile salt hydrolase (BSH) in Bacteroides is considered a potential drug target for obesity-related metabolic diseases, but its involvement in colon tumorigenesis has not been explored. BSH-expressing Bacteroides is found at high abundance in the stools of colorectal cancer (CRC) patients  with overweight and in the feces of a high-fat diet (HFD)-induced CRC mouse model. Colonization of B. fragilis 638R, a strain with low BSH activity, overexpressing a recombinant bsh gene from B. fragilis NCTC9343 strain, results in increased unconjugated bile acids in the colon and accelerated progression of CRC under HFD treatment. In the presence of high BSH activity, the resultant elevation of unconjugated deoxycholic acid and lithocholic acid activates the G-protein-coupled bile acid receptor, resulting in increased ß-catenin-regulated chemokine (C-C motif) ligand 28 (CCL28) expression in colon tumors. Activation of the ß-catenin/CCL28 axis leads to elevated intra-tumoral immunosuppressive CD25+FOXP3+ Treg cells. Blockade of the ß-catenin/CCL28 axis releases the immunosuppression to enhance the intra-tumoral anti-tumor response, which decreases CRC progression under HFD treatment. Pharmacological inhibition of BSH reduces HFD-accelerated CRC progression, coincident with suppression of the ß-catenin/CCL28 pathway. These findings provide insights into the pro-carcinogenetic role of Bacteroides in obesity-related CRC progression and characterize BSH as a potential target for CRC prevention and treatment.

Extracellular matrix profiles determine risk and prognosis of the squamous cell carcinoma subtype of non-small cell lung carcinoma.

BACKGROUND: Squamous cell carcinoma (SqCC) is a subtype of non-small cell lung cancer for which patient prognosis remains poor. The extracellular matrix (ECM) is critical in regulating cell behavior; however, its importance in tumor aggressiveness remains to be comprehensively characterized. METHODS: Multi-omics data of SqCC human tumor specimens was combined to characterize ECM features associated with initiation and recurrence. Penalized logistic regression was used to define a matrix risk signature for SqCC tumors and its performance across a panel of tumor types and in SqCC premalignant lesions was evaluated. Consensus clustering was used to define prognostic matreotypes for SqCC tumors. Matreotype-specific tumor biology was defined by integration of bulk RNAseq with scRNAseq data, cell type deconvolution, analysis of ligand-receptor interactions and enriched biological pathways, and through cross comparison of matreotype expression profiles with aging and idiopathic pulmonary fibrosis lung profiles. RESULTS: This analysis revealed subtype-specific ECM signatures associated with tumor initiation that were predictive of premalignant progression. We identified an ECM-enriched tumor subtype associated with the poorest prognosis. In silico analysis indicates that matrix remodeling programs differentially activate intracellular signaling in tumor and stromal cells to reinforce matrix remodeling associated with resistance and progression. The matrix subtype with the poorest prognosis resembles ECM remodeling in idiopathic pulmonary fibrosis and may represent a field of cancerization associated with elevated cancer risk. CONCLUSIONS: Collectively, this analysis defines matrix-driven features of poor prognosis to inform precision medicine prevention and treatment strategies towards improving SqCC patient outcome.

Creatine riboside is a cancer cell-derived metabolite associated with arginine auxotrophy.

The metabolic dependencies of cancer cells have substantial potential to be exploited to improve the diagnosis and treatment of cancer. Creatine riboside (CR) is identified as a urinary metabolite associated with risk and prognosis in lung and liver cancer. However, the source of high CR levels in patients with cancer as well as their implications for the treatment of these aggressive cancers remain unclear. By integrating multiomics data on lung and liver cancer, we have shown that CR is a cancer cell-derived metabolite. Global metabolomics and gene expression analysis of human tumors and matched liquid biopsies, together with functional studies, revealed that dysregulation of the mitochondrial urea cycle and a nucleotide imbalance were associated with high CR levels and indicators of a poor prognosis. This metabolic phenotype was associated with reduced immune infiltration and supported rapid cancer cell proliferation that drove aggressive tumor growth. CRhi cancer cells were auxotrophic for arginine, revealing a metabolic vulnerability that may be exploited therapeutically. This highlights the potential of CR not only as a poor-prognosis biomarker but also as a companion biomarker to inform the administration of arginine-targeted therapies in precision medicine strategies to improve survival for patients with cancer.

Δ133p53α Protects Human Astrocytes from Amyloid-beta Induced Senescence and Neurotoxicity.

Cellular senescence is an important contributor to aging and age-related diseases such as Alzheimer's disease (AD). Senescent cells are characterized by a durable cell proliferation arrest and the acquisition of a proinflammatory senescence-associated secretory phenotype (SASP), which participates in the progression of neurodegenerative disorders. Clearance of senescent glial cells in an AD mouse model prevented cognitive decline suggesting pharmacological agents targeting cellular senescence might provide novel therapeutic approaches for AD. Δ133p53α, a natural protein isoform of p53, was previously shown to be a negative regulator of cellular senescence in primary human astrocytes, with clinical implications from its diminished expression in brain tissues from AD patients. Here we show that treatment of proliferating human astrocytes in culture with amyloid-beta oligomers (Aß), an endogenous pathogenic agent of AD, results in reduced expression of Δ133p53α, as well as induces the cells to become senescent and express proinflammatory SASP cytokines such as IL-6, IL-1ß and TNFα. Our data suggest that Aß-induced astrocyte cellular senescence is associated with accelerated DNA damage, and upregulation of full-length p53 and its senescence-inducing target gene p21WAF1. We also show that exogenously enhanced expression of Δ133p53α rescues human astrocytes from Aß-induced cellular senescence and SASP through both protection from DNA damage and dominant-negative inhibition of full-length p53, leading to inhibition of Aß-induced, astrocyte-mediated neurotoxicity. The results presented here demonstrate that Δ133p53α manipulation could modulate cellular senescence in the context of AD, possibly opening new therapeutic avenues.

MPAPASS software enables stitched multiplex, multidimensional EV repertoire analysis and a standard framework for reporting bead-based assays.

Extracellular vesicles (EVs) of various types are released or shed from all cells. EVs carry proteins and contain additional protein and nucleic acid cargo that relates to their biogenesis and cell of origin. EV cargo in liquid biopsies is of widespread interest owing to its ability to provide a retrospective snapshot of cell state at the time of EV release. For the purposes of EV cargo analysis and repertoire profiling, multiplex assays are an essential tool in multiparametric analyte studies but are still being developed for high-parameter EV protein detection. Although bead-based EV multiplex analyses offer EV profiling capabilities with conventional flow cytometers, the utilization of EV multiplex assays has been limited by the lack of software analysis tools for such assays. To facilitate robust EV repertoire studies, we developed multiplex analysis post-acquisition analysis (MPAPASS) open-source software for stitched multiplex analysis, EV database-compatible reporting, and visualization of EV repertoires.

Accounting for EGFR Mutations in Epidemiologic Analyses of Non-Small Cell Lung Cancers: Examples Based on the International Lung Cancer Consortium Data.

BACKGROUND: Somatic EGFR mutations define a subset of non-small cell lung cancers (NSCLC) that have clinical impact on NSCLC risk and outcome. However, EGFR-mutation-status is often missing in epidemiologic datasets. We developed and tested pragmatic approaches to account for EGFR-mutation-status based on variables commonly included in epidemiologic datasets and evaluated the clinical utility of these approaches. METHODS: Through analysis of the International Lung Cancer Consortium (ILCCO) epidemiologic datasets, we developed a regression model for EGFR-status; we then applied a clinical-restriction approach using the optimal cut-point, and a second epidemiologic, multiple imputation approach to ILCCO survival analyses that did and did not account for EGFR-status. RESULTS: Of 35,356 ILCCO patients with NSCLC, EGFR-mutation-status was available in 4,231 patients. A model regressing known EGFR-mutation-status on clinical and demographic variables achieved a concordance index of 0.75 (95% CI, 0.74-0.77) in the training and 0.77 (95% CI, 0.74-0.79) in the testing dataset. At an optimal cut-point of probability-score = 0.335, sensitivity = 69% and specificity = 72.5% for determining EGFR-wildtype status. In both restriction-based and imputation-based regression analyses of the individual roles of BMI on overall survival of patients with NSCLC, similar results were observed between overall and EGFR-mutation-negative cohort analyses of patients of all ancestries. However, our approach identified some differences: EGFR-mutated Asian patients did not incur a survival benefit from being obese, as observed in EGFR-wildtype Asian patients. CONCLUSIONS: We introduce a pragmatic method to evaluate the potential impact of EGFR-status on epidemiological analyses of NSCLC. IMPACT: The proposed method is generalizable in the common occurrence in which EGFR-status data are missing.

Integration of adeno-associated virus (AAV) into the genomes of most Thai and Mongolian liver cancer patients does not induce oncogenesis.

BACKGROUND: Engineered versions of adeno-associated virus (AAV) are commonly used in gene therapy but evidence revealing a potential oncogenic role of natural AAV in hepatocellular carcinoma (HCC) has raised concerns. The frequency of potentially oncogenic integrations has been reported in only a few populations. AAV infection and host genome integration in another type of liver cancer, cholangiocarcinoma (CCA), has been studied only in one cohort. All reported oncogenic AAV integrations in HCC come from strains resembling the fully sequenced AAV2 and partly sequenced AAV13. When AAV integration occurs, only a fragment of the AAV genome is detectable in later DNA or RNA sequencing. The integrated fragment is typically from the 3' end of the AAV genome, and this positional bias has been only partly explained. Three research groups searched for evidence of AAV integration in HCC RNAseq samples in the Cancer Genome Atlas (TCGA) but reported conflicting results. RESULTS: We collected and analyzed whole transcriptome and viral capture DNA sequencing in paired tumor and non-tumor samples from two liver cancer Asian cohorts from Thailand (N = 147, 47 HCC and 100 intrahepatic cholangiocarcinoma (iCCA)) and Mongolia (N = 70, all HCC). We found only one HCC patient with a potentially oncogenic integration of AAV, in contrast to higher frequency reported in European patients. There were no oncogenic AAV integrations in iCCA patients. AAV genomic segments are present preferentially in the non-tumor samples of Thai patients. By analyzing the AAV genome positions of oncogenic and non-oncogenic integrated fragments, we found that almost all the putative oncogenic integrations overlap the X gene, which is present and functional only in the strain AAV2 among all fully sequenced strains. This gene content difference could explain why putative oncogenic integrations from other AAV strains have not been reported. We resolved the discrepancies in previous analyses of AAV presence in TCGA HCC samples and extended it to CCA. There are 12 TCGA samples with an AAV segment and none are in Asian patients. AAV segments are present in preferentially in TCGA non-tumor samples, like what we observed in the Thai patients. CONCLUSIONS: Our findings suggest a minimal AAV risk of hepatocarcinogenesis in Asian liver cancer patients. The partial genome presence and positional bias of AAV integrations into the human genome has complicated analysis of possible roles of AAV in liver cancer.