CD4+ T cell activation distinguishes response to anti-PD-L1+anti-CTLA4 therapy from anti-PD-L1 monotherapy.

Cancer patients often receive a combination of antibodies targeting programmed death-ligand 1 (PD-L1) and cytotoxic T lymphocyte antigen-4 (CTLA4). We conducted a window-of-opportunity study in head and neck squamous cell carcinoma (HNSCC) to examine the contribution of anti-CTLA4 to anti-PD-L1 therapy. Single-cell profiling of on- versus pre-treatment biopsies identified T cell expansion as an early response marker. In tumors, anti-PD-L1 triggered the expansion of mostly CD8+ T cells, whereas combination therapy expanded both CD4+ and CD8+ T cells. Such CD4+ T cells exhibited an activated T helper 1 (Th1) phenotype. CD4+ and CD8+ T cells co-localized with and were surrounded by dendritic cells expressing T cell homing factors or antibody-producing plasma cells. T cell receptor tracing suggests that anti-CTLA4, but not anti-PD-L1, triggers the trafficking of CD4+ naive/central-memory T cells from tumor-draining lymph nodes (tdLNs), via blood, to the tumor wherein T cells acquire a Th1 phenotype. Thus, CD4+ T cell activation and recruitment from tdLNs are hallmarks of early response to anti-PD-L1 plus anti-CTLA4 in HNSCC.

PHGDH heterogeneity potentiates cancer cell dissemination and metastasis.

Cancer metastasis requires the transient activation of cellular programs enabling dissemination and seeding in distant organs1. Genetic, transcriptional and translational heterogeneity contributes to this dynamic process2,3. Metabolic heterogeneity has also been observed4, yet its role in cancer progression is less explored. Here we find that the loss of phosphoglycerate dehydrogenase (PHGDH) potentiates metastatic dissemination. Specifically, we find that heterogeneous or low PHGDH expression in primary tumours of patients with breast cancer is associated with decreased metastasis-free survival time. In mice, circulating tumour cells and early metastatic lesions are enriched with Phgdhlow cancer cells, and silencing Phgdh in primary tumours increases metastasis formation. Mechanistically, Phgdh interacts with the glycolytic enzyme phosphofructokinase, and the loss of this interaction activates the hexosamine-sialic acid pathway, which provides precursors for protein glycosylation. As a consequence, aberrant protein glycosylation occurs, including increased sialylation of integrin αvß3, which potentiates cell migration and invasion. Inhibition of sialylation counteracts the metastatic ability of Phgdhlow cancer cells. In conclusion, although the catalytic activity of PHGDH supports cancer cell proliferation, low PHGDH protein expression non-catalytically potentiates cancer dissemination and metastasis formation. Thus, the presence of PHDGH heterogeneity in primary tumours could be considered a sign of tumour aggressiveness.

Regeneration Defects in Yap and Taz Mutant Mouse Livers Are Caused by Bile Duct Disruption and Cholestasis.

BACKGROUND AND AIMS: The Hippo pathway and its downstream effectors YAP and TAZ (YAP/TAZ) are heralded as important regulators of organ growth and regeneration. However, different studies provided contradictory conclusions about their role during regeneration of different organs, ranging from promoting proliferation to inhibiting it. Here we resolve the function of YAP/TAZ during regeneration of the liver, where Hippo's role in growth control has been studied most intensely. METHODS: We evaluated liver regeneration after carbon tetrachloride toxic liver injury in mice with conditional deletion of Yap/Taz in hepatocytes and/or biliary epithelial cells, and measured the behavior of different cell types during regeneration by histology, RNA sequencing, and flow cytometry. RESULTS: We found that YAP/TAZ were activated in hepatocytes in response to carbon tetrachloride toxic injury. However, their targeted deletion in adult hepatocytes did not noticeably impair liver regeneration. In contrast, Yap/Taz deletion in adult bile ducts caused severe defects and delay in liver regeneration. Mechanistically, we showed that Yap/Taz mutant bile ducts degenerated, causing cholestasis, which stalled the recruitment of phagocytic macrophages and the removal of cellular corpses from injury sites. Elevated bile acids activated pregnane X receptor, which was sufficient to recapitulate the phenotype observed in mutant mice. CONCLUSIONS: Our data show that YAP/TAZ are practically dispensable in hepatocytes for liver development and regeneration. Rather, YAP/TAZ play an indirect role in liver regeneration by preserving bile duct integrity and securing immune cell recruitment and function.

c-MET/VEGFR-2 co-localisation impacts on survival following bevacizumab therapy in epithelial ovarian cancer: an exploratory biomarker study of the phase 3 ICON7 trial.

INTRODUCTION: Bevacizumab improves survival outcomes in women diagnosed with epithelial ovarian cancer (EOC). Pre-clinical data showed that the c-MET/VEGFR-2 heterocomplex negates VEGF inhibition through activation of c-MET signalling, leading to a more invasive and metastatic phenotype. We evaluated the clinical significance of c-MET and VEGFR-2 co-localisation and its association with VEGF pathway-related single nucleotide polymorphisms (SNPs) in women participating in the phase 3 trial, ICON7 (ClinicalTrials.gov identifier: NCT00262847). MATERIALS AND METHODS: Patients had FIGO stage I-IIA grade 3/poorly differentiated or clear cell carcinoma or stage IIB-IV epithelial ovarian, primary peritoneal or fallopian tube cancer. Immunofluorescence staining for co-localised c-MET and VEGFR-2 on tissue microarrays and genotyping of germline DNA from peripheral blood leukocytes for VEGFA and VEGFR-2 SNPs was performed. The significance of these biomarkers was assessed against survival. RESULTS: Tissue microarrays from 178 women underwent immunofluorescence staining. Multivariable analysis showed that greater c-MET/VEGFR-2 co-localisation predicted worse OS in patients treated with bevacizumab after adjusting for FIGO stage and debulking surgery outcome (hazard ratio [HR] 1.034, 95% confidence interval [95%CI] 1.010-1.059). Women in the c-MET/VEGFR-2HIGH group treated with bevacizumab demonstrated significantly reduced OS (39.3 versus > 60 months; HR 2.00, 95%CI 1.08-3.72). Germline DNA from 449 women underwent genotyping. In the bevacizumab group, those women with the VEGFR-2 rs2305945 G/G variant had a trend towards shorter PFS compared with G/T or T/T variants (18.3 versus 23.0 months; HR 0.74, 95%CI 0.53-1.03). CONCLUSIONS: In bevacizumab-treated women diagnosed with EOC, high c-MET/VEGFR-2 co-localisation on tumour tissue and the VEGFR-2 rs2305945 G/G variant, which may be biologically related, were associated with worse survival outcomes.

The genomic landscape of nonsmall cell lung carcinoma in never smokers.

Lung cancer is the number one cause of cancer-related death worldwide with cigarette smoking as its major risk factor. Although the incidence of lung cancer in never smokers is rising, this subgroup of patients is underrepresented in genomic studies of lung cancer. Here, we assembled a prospective cohort of 46 never-smoking, nonsmall cell lung cancer (NSCLC) patients and performed whole-exome and low-coverage whole-genome sequencing on tumors and matched germline DNA. We observed fewer somatic mutations, genomic breakpoints and a smaller fraction of the genome with chromosomal instability in lung tumors from never smokers compared to smokers. The lower number of mutations, enabled us to identify TSC22D1 as a potential driver gene in NSCLC. On the other hand, the frequency of mutations in actionable genes such as EGFR and ERBB2 and of amplifications in MET were higher, while the mutation rate of TP53, which is a negative prognostic factor, was lower in never smokers compared to smokers. Together, these observations suggest a more favorable prognosis for never smokers with NSCLC. Classification of somatic mutations into six-substitution type patterns or into 96-substitution type signatures revealed distinct clusters between smokers and never smokers. Particularly, we identified in never smokers signatures related to aging, homologous recombination damage and APOBEC/AID activity as the most important underlying processes of NSCLC. This further indicates that second-hand smoking is not driving NSCLC pathogenesis in never smokers.

Stable or seral? Fire-driven alternative states in aspen forests of western North America.

As important centres for biological diversity, aspen forests are essential to the function and aesthetics of montane ecosystems in western North America. Aspen stands are maintained by a nuanced relationship with wildfire, although in recent decades aspen mortality has increased. The need to understand the baseline environmental conditions that favour aspen is clear; however, long-term fire history reconstructions are rare due to the scarcity of natural archives in dry montane settings. Here, we analyse a high-resolution lake sediment record from southwestern, Utah, USA to quantify the compositional and burning conditions that promote stable (or seral) aspen forests. Our results show that aspen presence is negatively correlated with subalpine fir and that severe fires tend to promote persistent and diverse aspen ecosystems over centennial timescales. This information improves our understanding of aspen disturbance ecology and identifies the circumstances where critical transitions in montane forests may occur.

Uncovering the genomic heterogeneity of multifocal breast cancer.

Multifocal breast cancer (MFBC), defined as multiple synchronous unilateral lesions of invasive breast cancer, is relatively frequent and has been associated with more aggressive features than unifocal cancer. Here, we aimed to investigate the genomic heterogeneity between MFBC lesions sharing similar histopathological parameters. Characterization of different lesions from 36 patients with ductal MFBC involved the identification of non-silent coding mutations in 360 protein-coding genes (171 tumour and 36 matched normal samples). We selected only patients with lesions presenting the same grade, ER, and HER2 status. Mutations were classified as 'oncogenic' in the case of recurrent substitutions reported in COSMIC or truncating mutations affecting tumour suppressor genes. All mutations identified in a given patient were further interrogated in all samples from that patient through deep resequencing using an orthogonal platform. Whole-genome rearrangement screen was further conducted in 8/36 patients. Twenty-four patients (67%) had substitutions/indels shared by all their lesions, of which 11 carried the same mutations in all lesions, and 13 had lesions with both common and private mutations. Three-quarters of those 24 patients shared oncogenic variants. The remaining 12 patients (33%) did not share any substitution/indels, with inter-lesion heterogeneity observed for oncogenic mutation(s) in genes such as PIK3CA, TP53, GATA3, and PTEN. Genomically heterogeneous lesions tended to be further apart in the mammary gland than homogeneous lesions. Genome-wide analyses of a limited number of patients identified a common somatic background in all studied MFBCs, including those with no mutation in common between the lesions. To conclude, as the number of molecular targeted therapies increases and trials driven by genomic screening are ongoing, our findings highlight the presence of genomic inter-lesion heterogeneity in one-third, despite similar pathological features. This implies that deeper molecular characterization of all MFBC lesions is warranted for the adequate management of those cancers.

Genetic variants in VEGF pathway genes in neoadjuvant breast cancer patients receiving bevacizumab: Results from the randomized phase III GeparQuinto study.

Studies assessing the effect of bevacizumab (BEV) on breast cancer (BC) outcome have shown different effects on progression-free and overall survival, suggesting that a subgroup of patients may benefit from this treatment. Unfortunately, no biomarkers exist to identify these patients. Here, we investigate whether single nucleotide polymorphisms (SNPs) in VEGF pathway genes correlate with pathological complete response (pCR) in the neoadjuvant GeparQuinto trial. HER2-negative patients were randomized into treatment arms receiving either BEV combined with standard chemotherapy or chemotherapy alone. In a pre-planned biomarker study, DNA was collected from 729 and 724 patients, respectively from both treatment arms, and genotyped for 125 SNPs. Logistic regression assessed interaction between individual SNPs and both treatment arms to predict pCR. Five SNPs may be associated with a better response to BEV, but none of them remained significant after correction for multiple testing. The two SNPs most strongly associated, rs833058 and rs699947, were located upstream of the VEGF-A promoter. Odds ratios for the homozygous common, heterozygous and homozygous rare rs833058 genotypes were 2.36 (95% CI, 1.49-3.75), 1.20 (95% CI, 0.88-1.64) and 0.61 (95% CI, 0.34-1.12). Notably, some SNPs in VEGF-A exhibited a more pronounced effect in the triple-negative subgroup. Several SNPs in VEGF-A may be associated with improved pCR when receiving BEV in the neoadjuvant setting. Although none of the observed effects survived correction for multiple testing, our observations are consistent with previous studies on BEV efficacy in BC. Further research is warranted to clarify the predictive value of these markers.

Efflux pump ABCB1 single nucleotide polymorphisms and dose reductions in patients with metastatic renal cell carcinoma treated with sunitinib.

UNLABELLED: There is growing evidence that sunitinib plasma levels have an impact on treatment outcome in patients with metastatic renal cell carcinoma (mRCC). We studied the impact of single nucleotide polymorphisms (SNPs) in genes involved in sunitinib pharmacokinetics, and additionally, sunitinib pharmacodynamics on dose reductions of the tyrosine kinase inhibitor. METHODS: We retrospectively analyzed germ-line DNA retrieved from mRCC patients receiving sunitinib as first-line therapy. We genotyped 11 key SNPs, respectively, in ABCB1, NR1/2, NR1/3 and CYP3A5, involved in sunitinib pharmacokinetics as well as VEGFR1 and VEGFR3, which have been suggested as regulators of sunitinib pharmacodynamics. Association between these SNPs and time-to-dose-reduction (TTDR) was studied by Cox regression. RESULTS: We identified 96 patients who were treated with sunitinib and from whom germ-line DNA and data on dose reductions were available. We observed an increased TTDR in patients carrying the TT-genotype in ABCB1 rs1125803 compared to patients with CC- or CT-genotypes (19 vs. 7 cycles; p = 0.031 on univariate analysis and p = 0.012 on multivariate analysis) and an increased TTDR in patients carrying the TT/TA-variant in ABCB1 rs2032582 compared to patients with the GG- or GT/GA-variant (19 vs. 7 cycles; p = 0.046 on univariate analysis and p = 0.024 on multivariate analysis). CONCLUSION: mRCC patients carrying the rs1128503 TT-variant or the TT/TA-variant in rs2032582 in ABCB1, which encodes for an efflux pump, do require less dose reductions due to adverse events compared to patients with the wild type or heterozygote variants in these genes.

Single-cell Atlas of Penile Cancer Reveals TP53 Mutations as a Driver of an Aggressive Phenotype, Irrespective of Human Papillomavirus Status, and Provides Clues for Treatment Personalization.

BACKGROUND AND OBJECTIVE: TP53 loss-of-function (TP53LOF) mutations might be a driver of poor prognosis and chemoresistance in both human papillomavirus (HPV)-independent (HPV-) and HPV-associated (HPV+) penile squamous cell carcinoma (PSCC). Here, we aim to describe transcriptomic differences in the PSCC microenvironment stratified by TP53LOF and HPV status. METHODS: We used single-cell RNA sequencing (scRNA-seq) and T-cell receptor sequencing to obtain a comprehensive atlas of the cellular architecture of PSCC. TP53LOF and HPV status were determined by targeted next-generation sequencing and sequencing HPV-DNA reads. Six HPV+ TP53 wild type (WT), six HPV- TP53WT, and four TP53LOF PSCC samples and six controls were included. Immunohistochemistry and hematoxylin-eosin confirmed the morphological context of the observed signatures. Prognostic differences between patient groups were validated in 541 PSCC patients using Kaplan-Meier survival estimates. KEY FINDINGS AND LIMITATIONS: Patients with aberrant p53 staining fare much worse than patients with either HPV- or HPV+ tumors and WT p53 expression. Using scRNA-seq, we revealed 65 cell subtypes within 83 682 cells. TP53LOF tumors exhibit a partial epithelial-to-mesenchymal transition, immune-excluded, angiogenic, and morphologically invasive environment, underlying their aggressive phenotype. HPV- TP53WT tumors show stemness and immune exhaustion. HPV+ TP53WT tumors mirror normal epithelial maturation with upregulation of antibody-drug-conjugate targets and activation of innate immunity. Inherent to the scRNA-seq analysis, low sample size is a limitation and validation of signatures in large PSCC cohorts is needed. CONCLUSIONS AND CLINICAL IMPLICATIONS: This first scRNA-seq atlas offers unprecedented in-depth insights into PSCC biology underlying prognostic differences based on TP53 and HPV status. Our findings provide clues for testing novel biomarker-driven therapies in PSCC. PATIENT SUMMARY: Here, we analyzed tissues of penile cancer at the level of individual cells, which helps us understand why patients who harbor a deactivating mutation in the TP53 gene do much worse than patients lacking such a mutation. Such an analysis may help us tailor future therapies based on TP53 gene mutations and human papillomavirus status of these tumors.

Transcriptomic characterization of the histopathological growth patterns in breast cancer liver metastases.

Metastatic breast cancer (mBC) remains incurable and liver metastases (LM) are observed in approximately 50% of all patients with mBC. In some cases, surgical resection of breast cancer liver metastases (BCLM) is associated with prolonged survival. However, there are currently no validated marker to identify these patients. The interactions between the metastatic cancer cells and the liver microenvironment result in two main histopathological growth patterns (HGP): replacement (r-HGP), characterized by a direct contact between the cancer cells and the hepatocytes, and desmoplastic (d-HGP), in which a fibrous rim surrounds the tumor cells. In patients who underwent resection of BCLM, the r-HGP is associated with a worse postoperative prognosis than the d-HGP. Here, we aim at unraveling the biological differences between these HGP within ten patients presenting both HGP within the same metastasis. The transcriptomic analyses reveal overexpression of genes involved in cell cycle, DNA repair, vessel co-option and cell motility in r-HGP while angiogenesis, wound healing, and several immune processes were found overexpressed in d-HGP LM. Understanding the biology of the LM could open avenues to refine treatment of BC patients with LM.

PI3K/mTOR inhibition induces tumour microenvironment remodelling and sensitises pS6high uterine leiomyosarcoma to PD-1 blockade.

BACKGROUND: Uterine leiomyosarcomas (uLMS) are aggressive tumours with poor prognosis and limited treatment options. Although immune checkpoint blockade (ICB) has proven effective in some 'challenging-to-treat' cancers, clinical trials showed that uLMS do not respond to ICB. Emerging evidence suggests that aberrant PI3K/mTOR signalling can drive resistance to ICB. We therefore explored the relevance of the PI3K/mTOR pathway for ICB treatment in uLMS and explored pharmacological inhibition of this pathway to sensitise these tumours to ICB. METHODS: We performed an integrated multiomics analysis based on TCGA data to explore the correlation between PI3K/mTOR dysregulation and immune infiltration in 101 LMS. We assessed response to PI3K/mTOR inhibitors in immunodeficient and humanized uLMS patient-derived xenografts (PDXs) by evaluating tumour microenvironment modulation using multiplex immunofluorescence. We explored response to single-agent and a combination of PI3K/mTOR inhibitors with PD-1 blockade in humanized uLMS PDXs. We mapped intratumoural dynamics using single-cell RNA/TCR sequencing of serially collected biopsies. RESULTS: PI3K/mTOR over-activation (pS6high) associated with lymphocyte depletion and wound healing immune landscapes in (u)LMS, suggesting it contributes to immune evasion. In contrast, PI3K/mTOR inhibition induced profound tumour microenvironment remodelling in an ICB-resistant humanized uLMS PDX model, fostering adaptive anti-tumour immune responses. Indeed, PI3K/mTOR inhibition induced macrophage repolarisation towards an anti-tumourigenic phenotype and increased antigen presentation on dendritic and tumour cells, but also promoted infiltration of PD-1+ T cells displaying an exhausted phenotype. When combined with anti-PD-1, PI3K/mTOR inhibition led to partial or complete tumour responses, whereas no response to single-agent anti-PD-1 was observed. Combination therapy reinvigorated exhausted T cells and induced clonal hyper-expansion of a cytotoxic CD8+ T-cell population supported by a CD4+ Th1 niche. CONCLUSIONS: Our findings indicate that aberrant PI3K/mTOR pathway activation contributes to immune escape in uLMS and provides a rationale for combining PI3K/mTOR inhibition with ICB for the treatment of this patient population.

Lower respiratory tract single-cell RNA sequencing and neutrophil extracellular trap profiling of COVID-19-associated pulmonary aspergillosis: a single centre, retrospective, observational study.

BACKGROUND: COVID-19-associated pulmonary aspergillosis (CAPA) is a severe superinfection with the fungus Aspergillus affecting patients who are critically ill with COVID-19. The pathophysiology and the role of neutrophil extracellular traps (NETs) in this infection are largely unknown. We aimed to characterise the immune profile, with a focus on neutrophils and NET concentrations, of critically ill patients with COVID-19, with or without CAPA. METHODS: We conducted a single-centre, retrospective, observational study in two patient cohorts, both recruited at University Hospitals Leuven, Belgium. We included adults aged 18 years or older who were admitted to the intensive care unit because of COVID-19 between March 31, 2020, and May 18, 2021, and who were included in the previous Contagious trial (NCT04327570). We investigated the immune cellular landscape of CAPA versus COVID-19 only by performing single-cell RNA sequencing (scRNA-seq) on bronchoalveolar lavage fluid. Bronchoalveolar lavage immune cell fractions were compared between patients with CAPA and patients with COVID-19 only. Additionally, we determined lower respiratory tract NET concentrations using biochemical assays in patients aged 18 years and older who were admitted to the intensive care unit because of severe COVID-19 between March 15, 2020, and Dec 31, 2021, for whom bronchoalveolar lavage was available in the hospital biobank. Bronchoalveolar lavage NET concentrations were compared between patients with CAPA and patients with COVID-19 only and integrated with existing data on immune mediators in bronchoalveolar lavage and 90-day mortality. FINDINGS: We performed scRNA-seq of bronchoalveolar lavage on 43 samples from 39 patients, of whom 36 patients (30 male and six female; 14 with CAPA) were included in downstream analyses. We performed bronchoalveolar lavage NET analyses in 59 patients (46 male and 13 female), of whom 26 had CAPA. By scRNA-seq, patients with CAPA had significantly lower neutrophil fractions than patients with COVID-19 only (16% vs 33%; p=0·0020). The remaining neutrophils in patients with CAPA preferentially followed a hybrid maturation trajectory characterised by expression of genes linked to antigen presentation, with enhanced transcription of antifungal effector pathways. Patients with CAPA also showed depletion of mucosal-associated invariant T cells, reduced T helper 1 and T helper 17 differentiation, and transcriptional defects in specific aspects of antifungal immunity in macrophages and monocytes. We observed increased formation of NETs in patients with CAPA compared with patients with COVID-19 only (DNA complexed with citrullinated histone H3 median 15 898 ng/mL [IQR 4588-86 419] vs 7062 ng/mL [775-14 088]; p=0·042), thereby explaining decreased neutrophil fractions by scRNA-seq. Low bronchoalveolar lavage NET concentrations were associated with increased 90-day mortality in patients with CAPA. INTERPRETATION: Qualitative and quantitative disturbances in monocyte, macrophage, B-cell, and T-cell populations could predispose patients with severe COVID-19 to develop CAPA. Hybrid neutrophils form a specialised response to CAPA, and an adequate neutrophil response to CAPA is a major determinant for survival in these patients. Therefore, measuring bronchoalveolar lavage NETs could have diagnostic and prognostic value in patients with CAPA. Clinicians should be wary of aspergillosis when using immunomodulatory therapy that might inhibit NETosis to treat patients with severe COVID-19. FUNDING: Research Foundation Flanders, KU Leuven, UZ Leuven, VIB, the Fundação para a Ciência e a Tecnologia, the European Regional Development Fund, la Caixa Foundation, the Flemish Government, and Horizon 2020.

Rapid autopsies to enhance metastatic research: the UPTIDER post-mortem tissue donation program.

Research on metastatic cancer has been hampered by limited sample availability. Here we present the breast cancer post-mortem tissue donation program UPTIDER and show how it enabled sampling of a median of 31 (range: 5-90) metastases and 5-8 liquids per patient from its first 20 patients. In a dedicated experiment, we show the mild impact of increasing time after death on RNA quality, transcriptional profiles and immunohistochemical staining in tumor tissue samples. We show that this impact can be counteracted by organ cooling. We successfully generated ex vivo models from tissue and liquid biopsies from distinct histological subtypes of breast cancer. We anticipate these and future findings of UPTIDER to elucidate mechanisms of disease progression and treatment resistance and to provide tools for the exploration of precision medicine strategies in the metastatic setting.

PD-1- CD45RA+ effector-memory CD8 T cells and CXCL10+ macrophages are associated with response to atezolizumab plus bevacizumab in advanced hepatocellular carcinoma.

The combination of atezolizumab plus bevacizumab (atezo/bev) has dramatically changed the treatment landscape of advanced HCC (aHCC), achieving durable responses in some patients. Using single-cell transcriptomics, we characterize the intra-tumoural and peripheral immune context of patients with aHCC treated with atezo/bev. Tumours from patients with durable responses are enriched for PDL1+ CXCL10+ macrophages and, based on cell-cell interaction analysis, express high levels of CXCL9/10/11 and are predicted to attract peripheral CXCR3+ CD8+ effector-memory T cells (CD8 TEM) into the tumour. Based on T cell receptor sharing and pseudotime trajectory analysis, we propose that CD8 TEM preferentially differentiate into clonally-expanded PD1- CD45RA+ effector-memory CD8+ T cells (CD8 TEMRA) with pronounced cytotoxicity. In contrast, in non-responders, CD8 TEM remain frozen in their effector-memory state. Finally, in responders, CD8 TEMRA display a high degree of T cell receptor sharing with blood, consistent with their patrolling activity. These findings may help understand the possible mechanisms underlying response to atezo/bev in aHCC.

Transcriptional and spatial profiling of the kidney allograft unravels a central role for FcyRIII+ innate immune cells in rejection.

Rejection remains the main cause of premature graft loss after kidney transplantation, despite the use of potent immunosuppression. This highlights the need to better understand the composition and the cell-to-cell interactions of the alloreactive inflammatory infiltrate. Here, we performed droplet-based single-cell RNA sequencing of 35,152 transcriptomes from 16 kidney transplant biopsies with varying phenotypes and severities of rejection and without rejection, and identified cell-type specific gene expression signatures for deconvolution of bulk tissue. A specific association was identified between recipient-derived FCGR3A+ monocytes, FCGR3A+ NK cells and the severity of intragraft inflammation. Activated FCGR3A+ monocytes overexpressed CD47 and LILR genes and increased paracrine signaling pathways promoting T cell infiltration. FCGR3A+ NK cells overexpressed FCRL3, suggesting that antibody-dependent cytotoxicity is a central mechanism of NK-cell mediated graft injury. Multiplexed immunofluorescence using 38 markers on 18 independent biopsy slides confirmed this role of FcγRIII+ NK and FcγRIII+ nonclassical monocytes in antibody-mediated rejection, with specificity to the glomerular area. These results highlight the central involvement of innate immune cells in the pathogenesis of allograft rejection and identify several potential therapeutic targets that might improve allograft longevity.

Intermediate molecular phenotypes to identify genetic markers of anthracycline-induced cardiotoxicity risk.

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex disease whose polygenic component is mainly unidentified. We propose that levels of intermediate molecular phenotypes in the myocardium associated with histopathological damage could explain CDA susceptibility; so that variants of genes encoding these intermediate molecular phenotypes could identify patients susceptible to this complication. A genetically heterogeneous cohort of mice generated by backcrossing (N = 165) was treated with doxorubicin and docetaxel. Cardiac histopathological damage was measured by fibrosis and cardiomyocyte size by an Ariol slide scanner. We determine intramyocardial levels of intermediate molecular phenotypes of CDA associated with histopathological damage and quantitative trait loci (ipQTLs) linked to them. These ipQTLs seem to contribute to the missing heritability of CDA because they improve the heritability explained by QTL directly linked to CDA (cda-QTLs) through genetic models. Genes encoding these molecular subphenotypes were evaluated as genetic markers of CDA in three cancer patient cohorts (N = 517) whose cardiac damage was quantified by echocardiography or Cardiac Magnetic Resonance. Many SNPs associated with CDA were found using genetic models. LASSO multivariate regression identified two risk score models, one for pediatric cancer patients and the other for women with breast cancer. Molecular intermediate phenotypes associated with heart damage can identify genetic markers of CDA risk, thereby allowing a more personalized patient management. A similar strategy could be applied to identify genetic markers of other complex trait diseases.

Intermediate Molecular Phenotypes to Identify Genetic Markers of Anthracycline-Induced Cardiotoxicity Risk.

Cardiotoxicity due to anthracyclines (CDA) affects cancer patients, but we cannot predict who may suffer from this complication. CDA is a complex trait with a polygenic component that is mainly unidentified. We propose that levels of intermediate molecular phenotypes (IMPs) in the myocardium associated with histopathological damage could explain CDA susceptibility, so variants of genes encoding these IMPs could identify patients susceptible to this complication. Thus, a genetically heterogeneous cohort of mice (n = 165) generated by backcrossing were treated with doxorubicin and docetaxel. We quantified heart fibrosis using an Ariol slide scanner and intramyocardial levels of IMPs using multiplex bead arrays and QPCR. We identified quantitative trait loci linked to IMPs (ipQTLs) and cdaQTLs via linkage analysis. In three cancer patient cohorts, CDA was quantified using echocardiography or Cardiac Magnetic Resonance. CDA behaves as a complex trait in the mouse cohort. IMP levels in the myocardium were associated with CDA. ipQTLs integrated into genetic models with cdaQTLs account for more CDA phenotypic variation than that explained by cda-QTLs alone. Allelic forms of genes encoding IMPs associated with CDA in mice, including AKT1, MAPK14, MAPK8, STAT3, CAS3, and TP53, are genetic determinants of CDA in patients. Two genetic risk scores for pediatric patients (n = 71) and women with breast cancer (n = 420) were generated using machine-learning Least Absolute Shrinkage and Selection Operator (LASSO) regression. Thus, IMPs associated with heart damage identify genetic markers of CDA risk, thereby allowing more personalized patient management.

CD8+ T cell metabolic rewiring defined by scRNA-seq identifies a critical role of ASNS expression dynamics in T cell differentiation.

T cells dynamically rewire their metabolism during an immune response. We applied single-cell RNA sequencing to CD8+ T cells activated and differentiated in vitro in physiological medium to resolve these metabolic dynamics. We identify a differential time-dependent reliance of activating T cells on the synthesis versus uptake of various non-essential amino acids, which we corroborate with functional assays. We also identify metabolic genes that potentially dictate the outcome of T cell differentiation, by ranking them based on their expression dynamics. Among them, we find asparagine synthetase (Asns), whose expression peaks for effector T cells and decays toward memory formation. Disrupting these expression dynamics by ASNS overexpression promotes an effector phenotype, enhancing the anti-tumor response of adoptively transferred CD8+ T cells in a mouse melanoma model. We thus provide a resource of dynamic expression changes during CD8+ T cell activation and differentiation, and identify ASNS expression dynamics as a modulator of CD8+ T cell differentiation.