Journal: bioRxiv

Article Title: WILD-seq: Clonal deconvolution of transcriptomic signatures of sensitivity and resistance to cancer therapeutics in vivo

doi: 10.1101/2021.12.09.471927

Figure Lengend Snippet: a. Lentiviral construct design. A PGK promoter drives expression of a transcript encoding zsGreen harbouring a WILD-seq barcode sequence in the 3’UTR. A spacer sequence and polyadenylation signal ensure that that the barcode is detectable as part of a standard oligo dT single cell RNA library preparation and sequencing pipeline. The barcode cassette comprises 2 distinct 12 nucleotide barcode sequences separated by a constant 20 nucleotide linker region. The library of barcode sequences was designed with Hamming distance 5 to allow for sequencing error correction. b. Schematic of WILD-seq method. Tumour cells are infected with the WILD-seq lentiviral library and an appropriate size population of zsGreen positive cells isolated, each of which will express a single unique WILD-seq barcode. This WILD-seq barcoded, heterogenous cell pool is then subjected to an intervention of interest (such as in vivo treatment of the implanted pool with a therapeutic agent) and subsequently analysed by single cell RNA sequencing using the 10X Genomics platform. An additional PCR amplification step is included that specifically enriches for the barcode sequence to increase the number of cells to which a WILD-seq barcode can be conclusively assigned. c. scRNA-seq of in vitro 4T1 WILD-seq cell pool. UMAP plot of in vitro cultured 4T1 WILD-seq cells. Cells for which a WILD-seq clonal barcode is identified are shown as dark grey or coloured spots. Cells which belong to five selected clonal lineages are highlighted. d. scRNA-seq of 4T1 WILD-seq tumours. UMAP plots of vehicle-treated 4T1 WILD-seq tumours generated by injecting the 4T1 WILD-seq pool into the mammary fatpad of BALB/c mice. Four independent experiments were performed each involving injection into 3 separate host animals. Six animals from experiments A and B received vehicle 1 (10% DMSO, 0.9%β- cyclodextrin) and six animals from experiments C and D received vehicle 2 (12.5%ethanol, 12.5% Kolliphor). e. Clonal representation. Proportion of tumour cells assigned to each clonal lineage based on the WILD-seq barcode (n = 1 for in vitro cultured cells, n = 6 for tumours from NSG mice, n = 12 for vehicle-treated tumours from BALB/c mice). Selected clones from the most abundant lineages are plotted. Data represents mean ± SEM. f. Principal component analysis of clonal transcriptomes. Pseudo-bulk analysis was performed by summing counts for all tumour cells expressing the same WILD-seq clonal barcode within an independent experiment. For in vivo tumour samples each point represents the combined cells from 3 animals. Principal component analysis of normalized pseudo-bulk count data showed separation of samples by origin with PC1 and PC2 and separation by clonality with PC3. g. Transcriptomic programs associated with principal components. The top/bottom 50 gene loadings of PC1, PC2 and PC3 were analysed using Enrichr – . h. Clonal transcriptomic signatures from vehicle-treated BALB/c tumours. An AUCell score enrichment was calculated for each clone and for each experiment by comparing cells of a specific clonal lineage of interest to all assigned tumour cells within the same experiment. All gene sets which showed consistent and statistically significant enrichment in one of the six most abundant clones across experiments are illustrated.

Article Snippet: The pHSW8 lentiviral backbone was constructed using a four-way Gibson Assembly (NEB) by inserting a reverse expression cassette, consisting of a PGK promoter, the zsGreen ORF, a cloning site for high-diversity barcode libraries and a synthetic polyA signal, into an empty pCCL-c-MNDU3-X backbone (#81071 Addgene).

Techniques: Construct, Expressing, Sequencing, RNA Library Preparation, Infection, Isolation, In Vivo, RNA Sequencing, Amplification, In Vitro, Cell Culture, Generated, Injection, Clone Assay

Journal: bioRxiv

Article Title: WILD-seq: Clonal deconvolution of transcriptomic signatures of sensitivity and resistance to cancer therapeutics in vivo

doi: 10.1101/2021.12.09.471927

Figure Lengend Snippet: 4T1 WILD-seq tumours were treated with the BET bromodomain inhibitor JQ1 or vehicle from 7 days post-implantation until endpoint (n = 4 mice per condition). Data represents mean ± SEM. b. scRNA-seq of JQ1-treated 4T1 WILD-seq tumours. UMAP plots of vehicle- or JQ1-treated 4T1 WILD-seq tumours. Combined cells from 2 independent experiments, each with 3 mice per condition are shown. Cells for which a WILD-seq clonal barcode is identified are shown as dark grey or coloured spots. Cells which belong to four selected clonal lineages are highlighted. c. JQ1-treatment results in a reduction in Cd8+ tumour-associated T-cells. Cells belonging to the T- cell compartment were computationally extracted from the single cell data and reclustered. Upper panels show combined UMAP plots from experiments A and B with Cd8a expression per cell illustrated enabling identification of the Cd8+ T cell cluster. Lower panels show neighborhood graphs of the results from differential abundance testing using Milo . Coloured nodes represent neighbourhoods with significantly different cell numbers between conditions (FDR < 0.05) and the layout of nodes is determined by the position of the neighborhood index cell in the UMAP panel above. Experiments A and B were analysed separately due to differences in cell numbers. d. Differential gene expression between JQ1- and vehicle-treated tumour cells. Single cell heatmap of expression for genes which are significantly and consistently down-regulated across clonal lineages (combined fisher p-value < 0.05 and mean logFC < −0.2 for both experiments).1600 cells are represented (400 per experiment/condition), grouped according to their clonal lineage. e. Differential gene set expression between JQ1 and vehicle-treated tumour cells. Median AUCell score per experiment/condition for selected gene sets. The 5 clonal lineages with the highest representation across experiments are shown. f. Clonal representation. Proportion of tumour cells assigned to each clonal lineage in experiment A based on the WILD-seq barcode (n = 3 tumours per condition). Clones which make up at least 2% of the assigned tumour cells under at least one condition are plotted. The most sensitive clone 473 is highlighted in blue and the most resistant clones 93, 439, 264 are highlighted in red. Data represents mean ± SEM. g. Clonal response to JQ1-treatment. Log 2 fold change in clonal proportions upon JQ1 treatment across experiments A and B. Fold change was calculated by comparing each JQ1- treated sample with the mean of the 3 corresponding vehicle-treated samples from the same experiment. p-value calculated by one-sample t-test vs a theoretical mean of 0. Data represents mean ± SEM. h. and i. Correlation of JQ1-response with baseline clonal transcriptomic signatures. Clonal gene set enrichment scores for vehicle-treated tumours were calculated by comparing cells of a specific clonal lineage of interest to all assigned tumour cells within the same experiment. Correlation between these scores and JQ1-treatment response (mean log 2 fold change clonal proportion JQ1 vs vehicle) was then calculated for each gene set. Selected gene sets with the highest positive or negative correlation values (Pearson correlation test) are shown. A positive correlation indicates a higher expression in resistant clones, whereas a negative correlation indicates a higher expression in sensitive clones. Resistant clonal lineages identified by barcodes 93, 264 and 439 were combined for the purpose of this analysis to have enough cells for analysis within the vehicle-treated samples.

Article Snippet: The pHSW8 lentiviral backbone was constructed using a four-way Gibson Assembly (NEB) by inserting a reverse expression cassette, consisting of a PGK promoter, the zsGreen ORF, a cloning site for high-diversity barcode libraries and a synthetic polyA signal, into an empty pCCL-c-MNDU3-X backbone (#81071 Addgene).

Techniques: Expressing, Gene Expression, Clone Assay

Journal: bioRxiv

Article Title: WILD-seq: Clonal deconvolution of transcriptomic signatures of sensitivity and resistance to cancer therapeutics in vivo

doi: 10.1101/2021.12.09.471927

Figure Lengend Snippet: a. Tumour growth curves with docetaxel treatment. 4T1 WILD-seq tumours were treated with docetaxel or vehicle (12.5% ethanol, 12.5% Kolliphor) from 7 days post-implantation for 2 weeks (n = 5 mice per condition). Dosing regimen was 12.5 mg/Kg docetaxel three times per week. Data represents mean ± SEM. b. scRNA-seq of docetaxel-treated 4T1 WILD- seq tumours. UMAP plots of vehicle- or docetaxel-treated 4T1 WILD-seq tumours. Combined cells from 2 independent experiments, each with 3 mice per condition are shown. Cells for which a WILD- seq clonal barcode is identified are shown as dark grey or coloured spots. Cells which belong to three selected clonal lineages are highlighted. c. Clonal representation. Proportion of tumour cells assigned to each clonal lineage in experiment C based on the WILD-seq barcode (n = 3 tumours per condition). Clones which make up at least 2% of the assigned tumour cells under at least one condition are plotted. The most sensitive clone 238 is highlighted in blue and the most resistant clone 679 is highlighted in red. Data represents mean ± SEM. d. Clonal response to docetaxel-treatment. Log 2 fold change in clonal proportions upon docetaxel treatment across experiments C and D. Fold change was calculated by comparing each docetaxel-treated sample with the mean of the 3 corresponding vehicle-treated samples from the same experiment. p-values calculated by one-sample t-test vs a theoretical mean of 0. Data represents mean ± SEM. e. and f. Correlation of docetaxel-response with baseline clonal transcriptomic signatures. Clonal gene set enrichment scores for vehicle-treated tumours were calculated by comparing cells of a specific clonal lineage of interest to all assigned tumour cells within the same experiment. Correlation between these scores and docetaxel-treatment response (mean log 2 fold change clonal proportion docetaxel vs vehicle) was then calculated for each gene set. Selected gene sets with the highest positive or negative correlation values (Pearson correlation test) are shown. A positive correlation indicates a higher expression in resistant clones, whereas a negative correlation indicates a higher expression in sensitive clones.

Article Snippet: The pHSW8 lentiviral backbone was constructed using a four-way Gibson Assembly (NEB) by inserting a reverse expression cassette, consisting of a PGK promoter, the zsGreen ORF, a cloning site for high-diversity barcode libraries and a synthetic polyA signal, into an empty pCCL-c-MNDU3-X backbone (#81071 Addgene).

Techniques: Clone Assay, Expressing

Journal: bioRxiv

Article Title: WILD-seq: Clonal deconvolution of transcriptomic signatures of sensitivity and resistance to cancer therapeutics in vivo

doi: 10.1101/2021.12.09.471927

Figure Lengend Snippet: a. D2A1 WILD-seq tumour growth curves with docetaxel treatment. D2A1 WILD-seq tumours were treated with docetaxel or vehicle from 7 days post-implantation for 2 weeks (n = 5 vehicle-treated mice, n = 4 docetaxel-treated mice). Data represents mean ± SEM. b. scRNA-seq of docetaxel-treated D2A1 WILD-seq tumours. UMAP plots of vehicle-treated D2A1 WILD-seq D2A1 tumours and reclustered barcoded-tumour cells from vehicle- and docetaxel-treated tumours. Combined cells from 3 mice per condition are shown. Cells for which a WILD-seq clonal barcode is identified are shown as dark grey or coloured spots. Cells which belong to five selected clonal lineages are highlighted. c. Comparison of EMT status of major 4T1 and D2A1 WILD-seq clones. Violin plot of AUCell scores from vehicle-treated tumour cells generated using the HOLLERN_EMT_BREAST_TUMOR_DN gene set, a set of genes that have low expression in murine mammary tumours of mesenchymal histology. 4T1 WILD-seq clones exhibit varying levels of expression of this geneset whereas D2A1 WILD-seq clones have consistently low levels of expression of these genes. d. Clonal representation. Proportion of tumour cells assigned to each clonal lineage based on the WILD-seq barcode (n = 3 tumours per condition). Clones which make up at least 2% of the assigned tumour cells under at least one condition are plotted. The most sensitive clones to docetaxel treatment 118, 2874 and 1072 are highlighted in blue and the most resistant clones 1240, 1197 and 751 are highlighted in red. Data represents mean ± SEM. e. Clonal transcriptomic signatures from vehicle-treated tumours. Heatmap of median AUCell scores per sample for each of the five most abundant clones. All gene sets which showed consistent and statistically significant enrichment (combined fisher p-value < 0.01 & mean log 2 enrichment > 0.1) in at least one of these clones are illustrated. f. Selected gene sets whose expression is associated with sensitivity to docetaxel. Median AUCell scores per sample for each of the five most abundant clones is plotted. g. Transcriptomic signatures associated with resistance to docetaxel. For vehicle-treated tumours, resistant clonal lineages identified by barcodes 1197, 751 and 1240 were combined to have enough cells for analysis. Gene sets with significantly enriched expression in these resistant clones in vehicle-treated tumours were determined (adjusted p-value < 0.01 & log 2 enrichment > 0.1). A heatmap of median AUCell scores per clone, per condition of these resistance-associated gene sets is plotted. h. Selected gene sets whose expression is enriched or depleted in resistant clones. Median AUCell scores per clone, per sample are plotted for samples with at least 20 cells per clone. Due to changes in clonal abundance with treatment some clones can only be assessed under vehicle- or docetaxel- treated conditions.

Article Snippet: The pHSW8 lentiviral backbone was constructed using a four-way Gibson Assembly (NEB) by inserting a reverse expression cassette, consisting of a PGK promoter, the zsGreen ORF, a cloning site for high-diversity barcode libraries and a synthetic polyA signal, into an empty pCCL-c-MNDU3-X backbone (#81071 Addgene).

Techniques: Comparison, Clone Assay, Generated, Expressing

Journal: Nature neuroscience

Article Title: Circuit and molecular architecture of a ventral hippocampal network

doi: 10.1038/s41593-020-0705-8

Figure Lengend Snippet: a. Experimental design. We mapped the projections from vCA1 to 7 downstream targets using MAPseq. Briefly, we injected a virus encoding random RNA sequences (barcodes) into vCA1 and allowed them to transport into the axonal processes. Then, we dissected seven known target regions and sequenced the barcodes present to determine which neurons projected to which region. The relative abundance of a barcode present in each area was used as a proxy for the relative projection strength, and an N X R matrix was generated for analysis (N-barcodes, R-brain regions). See for details on each step of the experimental pipeline. b. expression of MAPseq GFP at the injection site in vCA1, representative from n=4 mice. c. Heatmap of projection strengths (from N X R matrix) of all 2494 vCA1 neurons (from 12 mice) to 7 targets mapped with MAPseq. Each row is an individual neuron’s projection strength to each of the 7 target regions, normalized to the maximal value in that row, resulting in a projection strength scale from 0 to 1 (see ). d. Distribution of number of projection targets of vCA1 neurons mapped. e. Heat map of conditional probability for two regions (see ), the proportion of cells projecting to area B within the subset of cells that project to area A.

Article Snippet: Adult male and female mice (n=12) were anesthetized and injected with MAPseq Sindbis viral barcode library (3×10 10 GC/ml, diversity of 2×10 7 different barcode sequences, Cold Spring Harbor Laboratories) into the vCA1 region of vHPC at −3.20mm ML, −3.16mm AP, and −3.70 & −3.80mm DV from brain surface (with 96.6nL per DV site) for a total of 193.2nL virus.

Techniques: Injection, Virus, Generated, Expressing

Journal: Nature neuroscience

Article Title: Circuit and molecular architecture of a ventral hippocampal network

doi: 10.1038/s41593-020-0705-8

Figure Lengend Snippet: a. k-means clustering of MAPseq data. Clusters are plotted as a grid-structured scatter plot where each row is a cluster and the color of the dots along each column correspond to the normalized projection strength to a region b. Plots of normalized projection strength of all multiple projection barcode motifs from MAPseq data. As in each line (different color) is an individual neuron’s projection strength to each of the 7 target regions, normalized to the maximal value in that row, resulting in a projection strength scale from 0 to 1 (see ). Black line is mean projection strength for all neurons in that motif, and grey is SEM. Exact cell number for each motif is provided in the inset (observed and expected from null model, see ), and in Extended Data Table 2.

Article Snippet: Adult male and female mice (n=12) were anesthetized and injected with MAPseq Sindbis viral barcode library (3×10 10 GC/ml, diversity of 2×10 7 different barcode sequences, Cold Spring Harbor Laboratories) into the vCA1 region of vHPC at −3.20mm ML, −3.16mm AP, and −3.70 & −3.80mm DV from brain surface (with 96.6nL per DV site) for a total of 193.2nL virus.

Techniques: