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.

Isolation and Identification of a Urinary Biomarker for Lung Cancer: 27-Nor-5ß-Cholestane-3α,7α,12α,24R,25S Pentol Glucuronide and Its Deuterated Analog.

An untargeted metabolomic study identified four potential lung cancer diagnostic biomarkers in human urine. One of the potential biomarkers was an unidentified feature possessing a m/z value of 561+. "561+" was isolated from human urine and tentatively identified as 27-nor-5ß-cholestane-3α,7α,12α,24,25 pentol glucuronide with unknown C24,25 stereochemistry using 1H NMR and mass spectrometry. In a prior report, the C24,25 stereochemistry of the aglycone, 27-nor-5ß-cholestane-3α,7α,12α,24,25 pentol, was found to be 24S,25R through GC analysis of the acetonide-TMS derivative. An authentic sample was prepared and found not to have the same stereochemistry as "561+". To identify the C24,25 stereochemistry, four C24,C25 diastereoisomeric alcohols of 27-nor-5ß-cholestane-3α,7α,12α,24,25 pentol were prepared from chiral amino acids. Using an LCMS method, the C24,C25 stereochemistry of the "561+" aglycone was determined to be 24R,25S. With the correct aglycone in hand, it was coupled with glucuronic acid to complete the first reported synthesis of 27-nor-5ß-cholestane-3α,7α,12α,24R,25S pentol glucuronide. Deuterium labeled 27-nor-5ß-cholestane-3α,7α,12α,24R,25S pentol was also synthesized for use as an internal standard for MS quantitation.

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.

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

Mutations effects on p53 isoforms' activities remain largely unknown, although they are mutated in 92% of TP53 mutant cancers. 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, mutant Δ133p53α R273H increases glioblastoma cancer cells proliferation and invasion while the WT does not. Furthermore, while WT Δ133p53α reduces apoptosis to promote DNA repair, the mutant also reduces apoptosis but fails to maintain genomic stability.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, using TCGA data, we determined that IL4I1, IDO1 and AHR are significantly higher in GBMs compared to LGGs. IL4I1 expression is increased in mutant TP53 LGGs and GBMs, although only significantly in LGG. Importantly, high expression of all three genes in LGG and IL4I1 in GBM is significantly associated with poorer patients' survival. 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.

Association between duration of smoking abstinence before non-small-cell lung cancer diagnosis and survival: a retrospective, pooled analysis of cohort studies.

BACKGROUND: The association between duration of smoking abstinence before non-small-cell lung cancer (NSCLC) diagnosis and subsequent survival can influence public health messaging delivered in lung-cancer screening. We aimed to assess whether the duration of smoking abstinence before diagnosis of NSCLC is associated with improved survival. METHODS: In this retrospective, pooled analysis of cohort studies, we used 26 cohorts participating in Clinical Outcomes Studies of the International Lung Cancer Consortium (COS-ILCCO) at 23 hospitals. 16 (62%) were from North America, six (23%) were from Europe, three (12%) were from Asia, and one (4%) was from South America. Patients enrolled were diagnosed between June 1, 1983, and Dec 31, 2019. Eligible patients had smoking data before NSCLC diagnosis, epidemiological data at diagnosis (obtained largely from patient questionnaires), and clinical information (retrieved from medical records). Kaplan-Meier curves and multivariable Cox models (ie, adjusted hazard ratios [aHRs]) were generated with individual, harmonised patient data from the consortium database. We estimated overall survival for all causes, measured in years from diagnosis date until the date of the last follow-up or death due to any cause and NSCLC-specific survival. FINDINGS: Of 42 087 patients with NSCLC in the COS-ILCCO database, 21 893 (52·0%) of whom were male and 20 194 (48·0%) of whom were female, we excluded 4474 (10·6%) with missing data. Compared with current smokers (15 036 [40·0%] of 37 613), patients with 1-3 years of smoking abstinence before NSCLC diagnosis (2890 [7·7%]) had an overall survival aHR of 0·92 (95% CI 0·87-0·97), patients with 3-5 years of smoking abstinence (1114 [3·0%]) had an overall survival aHR of 0·90 (0·83-0·97), and patients with more than 5 years of smoking abstinence (10 841 [28·8%]) had an overall survival aHR of 0·90 (0·87-0·93). Improved NSCLC-specific survival was observed in 4301 (44%) of 9727 patients who had quit cigarette smoking and was significant at abstinence durations of more than 5 years (aHR 0·87, 95% CI 0·81-0·93). Results were consistent across age, sex, histology, and disease-stage distributions. INTERPRETATION: In this large, pooled analysis of cohort studies across Asia, Europe, North America, and South America, overall survival was improved in patients with NSCLC whose duration of smoking abstinence before diagnosis was as short as 1 year. These findings suggest that quitting smoking can improve overall survival, even if NSCLC is diagnosed at a later lung-cancer screening visit. These findings also support the implementation of public health smoking cessation strategies at any time. FUNDING: The Alan B Brown Chair, The Posluns Family Fund, The Lusi Wong Fund, and the Princess Margaret Cancer Foundation.

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.

Melanoma-intrinsic NR2F6 activity regulates antitumor immunity.

Nuclear receptors (NRs) are implicated in the regulation of tumors and immune cells. We identify a tumor-intrinsic function of the orphan NR, NR2F6, regulating antitumor immunity. NR2F6 was selected from 48 candidate NRs based on an expression pattern in melanoma patient specimens (i.e., IFN-γ signature) associated with positive responses to immunotherapy and favorable patient outcomes. Correspondingly, genetic ablation of NR2F6 in a mouse melanoma model conferred a more effective response to PD-1 therapy. NR2F6 loss in B16F10 and YUMM1.7 melanoma cells attenuated tumor development in immune-competent but not -incompetent mice via the increased abundance of effector and progenitor-exhausted CD8+ T cells. Inhibition of NACC1 and FKBP10, identified as NR2F6 effectors, phenocopied NR2F6 loss. Inoculation of NR2F6 KO mice with NR2F6 KD melanoma cells further decreased tumor growth compared with NR2F6 WT mice. Tumor-intrinsic NR2F6 function complements its tumor-extrinsic role and justifies the development of effective anticancer therapies.

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.