An optimised tissue disaggregation and data processing pipeline for characterising fibroblast phenotypes using single-cell RNA sequencing

Single-cell RNA sequencing (scRNA-Seq) provides a valuable platform for characterising multicellular ecosystems. Fibroblasts are a heterogeneous cell type involved in many physiological and pathological processes, but remain poorly-characterised. Analysis of fibroblasts is challenging: these cells are difficult to isolate from tissues, and are therefore commonly under-represented in scRNA-seq datasets. Here, we describe an optimised approach for fibroblast isolation from human lung tissues. We demonstrate the potential for this procedure in characterising stromal cell phenotypes using scRNA-Seq, analyse the effect of tissue disaggregation on gene expression, and optimise data processing to improve clustering quality. We also assess the impact of in vitro culture conditions on stromal cell gene expression and proliferation, showing that altering these conditions can skew phenotypes

Single-cell analysis reveals prognostic fibroblast subpopulations linked to molecular and immunological subtypes of lung cancer.

Fibroblasts are poorly characterised cells that variably impact tumour progression. Here, we use single cell RNA-sequencing, multiplexed immunohistochemistry and digital cytometry (CIBERSORTx) to identify and characterise three major fibroblast subpopulations in human non-small cell lung cancer: adventitial, alveolar and myofibroblasts. Alveolar and adventitial fibroblasts (enriched in control tissue samples) localise to discrete spatial niches in histologically normal lung tissue and indicate improved overall survival rates when present in lung adenocarcinomas (LUAD). Trajectory inference identifies three phases of control tissue fibroblast activation, leading to myofibroblast enrichment in tumour samples: initial upregulation of inflammatory cytokines, followed by stress-response signalling and ultimately increased expression of fibrillar collagens. Myofibroblasts correlate with poor overall survival rates in LUAD, associated with loss of epithelial differentiation, TP53 mutations, proximal molecular subtypes and myeloid cell recruitment. In squamous carcinomas myofibroblasts were not prognostic despite being transcriptomically equivalent. These findings have important implications for developing fibroblast-targeting strategies for cancer therapy

Modulation of the tumour promoting functions of cancer associated fibroblasts by phosphodiesterase type 5 inhibition increases the efficacy of chemotherapy in human preclinical models of esophageal adenocarcinoma

Background and aims: Esophageal adenocarcinoma (EAC) is chemoresistant in the majority of cases. The tumor-promoting biology of cancer associated fibroblasts (CAF) make them a target for novel therapies. Phosphodiesterase type 5 inhibitors (PDE5i) have been shown to regulate the activated fibroblast phenotype in benign disease. We investigated the potential for CAF modulation in EAC using PDE5i to enhance the efficacy of chemotherapy. Methods: EAC fibroblasts were treated with PDE5i and phenotypic effects examined using immunoblotting, immunohistochemistry, gel contraction, transwell invasion, organotypics, single cell RNAseq and shotgun proteomics. The combination of PDE5i with standard-of-care chemotherapy (Epirubicin, 5-Fluorouracil and Cisplatin) was tested for safety and efficacy in validated near-patient model systems (3D tumor growth assays (3D-TGAs) and patient derived xenograft (PDX) mouse models). Results: PDE5i treatment reduced alpha-SMA expression in CAFs by 50% (p<0.05), associated with a significant reduction in the ability of CAFs to contract collagen-1 gels and induce cancer cell invasion, (p<0.05). RNAseq and proteomic analysis of CAF and EAC cell lines revealed PDE5i specific regulation of pathways related to fibroblast activation and tumor promotion. 3D-TGA assays confirmed the importance of stromal cells to chemoresistance in EAC, which could be attenuated by PDE5i. Chemotherapy+PDE5i in PDX-bearing mice was safe and significantly reduced PDX tumor volume (p<0.05). Conclusion: PDE5 is a candidate for clinical trials to alter the fibroblast phenotype in esophageal cancer. We demonstrate the specificity of PDE5i for fibroblasts to prevent transdifferentiation and revert the CAF phenotype. Finally, we confirm the efficacy of PDE5i in combination with chemotherapy in close-to-patient in vitro and in vivo PDX-based model systems

Genomic Analysis of Response to Neoadjuvant Chemotherapy in Esophageal Adenocarcinoma.

Neoadjuvant therapy followed by surgery is the standard of care for locally advanced esophageal adenocarcinoma (EAC). Unfortunately, response to neoadjuvant chemotherapy (NAC) is poor (20-37%), as is the overall survival benefit at five years (9%). The EAC genome is complex and heterogeneous between patients, and it is not yet understood whether specific mutational patterns may result in chemotherapy sensitivity or resistance. To identify associations between genomic events and response to NAC in EAC, a comparative genomic analysis was performed in 65 patients with extensive clinical and pathological annotation using whole-genome sequencing (WGS). We defined response using Mandard Tumor Regression Grade (TRG), with responders classified as TRG1-2 (n = 27) and non-responders classified as TRG4-5 (n =38). We report a higher non-synonymous mutation burden in responders (median 2.08/Mb vs. 1.70/Mb, p = 0.036) and elevated copy number variation in non-responders (282 vs. 136/patient, p < 0.001). We identified copy number variants unique to each group in our cohort, with cell cycle (CDKN2A, CCND1), c-Myc (MYC), RTK/PIK3 (KRAS, EGFR) and gastrointestinal differentiation (GATA6) pathway genes being specifically altered in non-responders. Of note, NAV3 mutations were exclusively present in the non-responder group with a frequency of 22%. Thus, lower mutation burden, higher chromosomal instability and specific copy number alterations are associated with resistance to NAC

Genetics and prognostication in splenic marginal zone lymphoma: revelations from deep sequencing

PURPOSE: Mounting evidence supports the clinical significance of gene mutations and immunogenetic features in common mature B-cell malignancies.EXPERIMENTAL DESIGN: We undertook a detailed characterization of the genetic background of splenic marginal zone lymphoma (SMZL), using targeted re-sequencing and explored potential clinical implications in a multinational cohort of 175 SMZL patients.RESULTS: We identified recurrent mutations in TP53 (16%), KLF2 (12%), NOTCH2 (10%), TNFAIP3 (7%), MLL2 (11%), MYD88 (7%) and ARID1A (6%), all genes known to be targeted by somatic mutation in SMZL. KLF2 mutations were early, clonal events, enriched in patients with del(7q) and IGHV1-2*04 B-cell receptor immunoglobulins, and were associated with a short median time-to-first-treatment (0.12 vs. 1.11 yrs; P=0.01). In multivariate analysis mutations in NOTCH2 (HR 2.12, 95%CI 1.02-4.4, P=0.044) and 100% germline IGHV gene identity (HR 2.19, 95%CI 1.05-4.55, P=0.036) were independent markers of short time-to-first-treatment, while TP53 mutations were an independent marker of short overall survival (HR 2.36, 95% CI 1.08-5.2, P=0.03).CONCLUSIONS: We identify key associations between gene mutations and clinical outcome, demonstrating for the first time that NOTCH2 and TP53 gene mutations are independent markers of reduced treatment-free and overall survival, respectively.<br/

Evaluation of high-throughput genomic assays for the Fc gamma receptor locus

Cancer immunotherapy has been revolutionised by the use of monoclonal antibodies (mAb) that function through their interaction with Fc gamma receptors (FcγRs). The low-affinity FcγR genes are highly homologous, map to a complex locus at 1p23 and harbour single nucleotide polymorphisms (SNPs) and copy number variation (CNV) that can impact on receptor function and response to therapeutic mAbs. This complexity can hinder accurate characterisation of the locus. We therefore evaluated and optimised a suite of assays for the genomic analysis of the FcγR locus amenable to peripheral blood mononuclear cells and formalin-fixed paraffin-embedded (FFPE) material that can be employed in a high-throughput manner. Assessment of TaqMan genotyping for FCGR2A-131H/R, FCGR3A-158F/V and FCGR2B-232I/T SNPs demonstrated the need for additional methods to discriminate genotypes for the FCGR3A-158F/V and FCGR2B-232I/T SNPs due to sequence homology and CNV in the region. A multiplex ligation-dependent probe amplification assay provided high quality SNP and CNV data in PBMC cases, but there was greater data variability in FFPE material in a manner that was predicted by the BIOMED-2 multiplex PCR protocol. In conclusion, we have evaluated a suite of assays for the genomic analysis of the FcγR locus that are scalable for application in large clinical trials of mAb therapy. These assays will ultimately help establish the importance of FcγR genetics in predicting response to antibody therapeutics

V-ATPase inhibition decreases mutant androgen receptor activity in castrate-resistant prostate cancer

Prostate cancer (PC) is critically dependent on androgen receptor (AR) signaling. Despite initial responsiveness to androgen deprivation, most patients with advanced PC subsequently progress to a clinically aggressive castrate-resistant (CRPC) phenotype, typically associated with expression of splice-variant or mutant AR forms. Although current evidence suggests that the vacuolar-ATPase (V-ATPase), a multi-protein complex that catalyzes proton transport across intracellular and plasma membranes, influences wild-type AR function, the effect of V-ATPase inhibition on variant AR function is unknown. Inhibition of V-ATPase reduced AR function in wild-type and mutant AR luciferase reporter models. In hormone-sensitive PC cell lines (LNCaP, DuCaP) and mutant AR CRPC cell lines (22Rv1, LNCaP-F877L/T878A), V-ATPase inhibition using bafilomycin-A1 and concanamycin-A reduced AR expression, and expression of AR target genes, at mRNA and protein levels. Furthermore, combining chemical V-ATPase inhibition with the AR antagonist enzalutamide resulted in a greater reduction in AR downstream target expression than enzalutamide alone in LNCaP cells. To investigate the role of individual subunit isoforms, siRNA and CRISPR-Cas9 were used to target the V1C1 subunit in 22Rv1 cells. Whereas transfection with ATP6V1C1-targeted siRNA significantly reduced AR protein levels and function, CRISPR-Cas9-mediated V1C1 knockout showed no substantial change in AR expression, but a compensatory increase in protein levels of the alternate V1C2 isoform. Overall, these results indicate that V-ATPase dysregulation is directly linked to both hormone responsive and CRPC via impact on AR function. In particular, V-ATPase inhibition can reduce AR signaling regardless of mutant AR expression.</p

Association analysis of brain-derived neurotrophic factor gene polymorphisms in asthmatic families

Background: A role for neuronal modulation of inflammation and airway hyper-responsiveness has been well described in asthma, and neurotrophins provide the link between inflammation and neuronal dysfunction. Brain-derived neurotrophic factor (BDNF) is an important mediator in this interaction. The aim of this study was to analyze the possible relationship between polymorphisms of the gene encoding BDNF and susceptibility to asthma. Methods: 341 families with at least 2 siblings with asthma were genotyped for 4 BDNF polymorphisms (rs6265, rs2030324, rs988748 and rs7124442). Results: Analysis by family-based association tests revealed no significant association between any polymorphisms analyzed and asthma susceptibility. Furthermore, BDNF polymorphism was not associated with asthma-related phenotypes such as FEV1 % predicted, bronchial hyper-responsiveness, IgE level, asthma and atopy severity scores. Conclusions: Our results suggest that genetic variation in the BDNF gene does not contribute significantly to asthma susceptibility or severity.<br/