High-Throughput Multimodal Single-Cell Targeted DNA and Surface Protein Analysis Using the Mission Bio Tapestri Platform.

An important aspect of understanding cancer biology is to connect the diverse repertoire of genotype-to-phenotype displays in individual specimens and ultimately resolve disease course outcome through informative datasets. A focus of cancer genomics has strived to provide predictive capabilities using genomic information to further inform therapeutic strategies. The advent of single-cell sequencing and analysis now provides a route to decipher high-resolution genomic diversity in individual samples and facilitate detailed understanding of clonal evolution in clinical research settings. In addition to generating high-throughput single-cell genomic SNV and CNV data, this protocol describes a new analytical ability that adds a second dimension which provides for interrogation of surface protein marker expression. The first immediate application of this technology is quite suitable to heme cancer cell studies. This multimodal approach allows for correlation of diverse genomic signatures to key phenotypic biomarkers such as immunophenotypes in leukemic diseases.

Simultaneous DNA, RNA, and Protein Analysis from Single Cells Using a High-Throughput Microfluidic Workflow for Resolution of Genotype-to-Phenotype Modalities.

Understanding the genomic landscape of cancer in single cells can be valuable for the characterization of molecular events that drive evolution of tumorigenesis and fostering progress in identifying druggable regimens for patient treatment scenarios. We report a new approach to measure multiple modalities simultaneously from up to 10,000 individual cells using microfluidics paired with next-generation sequencing. Our procedure determines targeted protein levels, mRNA transcript levels, and somatic gDNA sequence variations including copy number variants. This approach can resolve over 20 proteins, 100s of targeted transcripts, and DNA amplicons.

Personalized Single-Cell Proteogenomics to Distinguish Acute Myeloid Leukemia from Non-Malignant Clonal Hematopoiesis.

Genetic mutations associated with acute myeloid leukemia (AML) also occur in age-related clonal hematopoiesis, often in the same individual. This makes confident assignment of detected variants to malignancy challenging. The issue is particularly crucial for AML post-treatment measurable residual disease monitoring, where results can be discordant between genetic sequencing and flow cytometry. We show here, that it is possible to distinguish AML from clonal hematopoiesis and to resolve the immunophenotypic identity of clonal architecture. To achieve this, we first design patient-specific DNA probes based on patient's whole-genome sequencing, and then use them for patient-personalized single-cell DNA sequencing with simultaneous single-cell antibody-oligonucleotide sequencing. Examples illustrate AML arising from DNMT3A and TET2 mutated clones as well as independently. The ability to personalize single-cell proteogenomic assessment for individual patients based on leukemia-specific genomic features has implications for ongoing AML precision medicine efforts.

Single-cell mutation analysis of clonal evolution in myeloid malignancies.

Myeloid malignancies, including acute myeloid leukaemia (AML), arise from the expansion of haematopoietic stem and progenitor cells that acquire somatic mutations. Bulk molecular profiling has suggested that mutations are acquired in a stepwise fashion: mutant genes with high variant allele frequencies appear early in leukaemogenesis, and mutations with lower variant allele frequencies are thought to be acquired later1-3. Although bulk sequencing can provide information about leukaemia biology and prognosis, it cannot distinguish which mutations occur in the same clone(s), accurately measure clonal complexity, or definitively elucidate the order of mutations. To delineate the clonal framework of myeloid malignancies, we performed single-cell mutational profiling on 146 samples from 123 patients. Here we show that AML is dominated by a small number of clones, which frequently harbour co-occurring mutations in epigenetic regulators. Conversely, mutations in signalling genes often occur more than once in distinct subclones, consistent with increasing clonal diversity. We mapped clonal trajectories for each sample and uncovered combinations of mutations that synergized to promote clonal expansion and dominance. Finally, we combined protein expression with mutational analysis to map somatic genotype and clonal architecture with immunophenotype. Our findings provide insights into the pathogenesis of myeloid transformation and how clonal complexity evolves with disease progression.

Benchmarking single-cell RNA-sequencing protocols for cell atlas projects.

Single-cell RNA sequencing (scRNA-seq) is the leading technique for characterizing the transcriptomes of individual cells in a sample. The latest protocols are scalable to thousands of cells and are being used to compile cell atlases of tissues, organs and organisms. However, the protocols differ substantially with respect to their RNA capture efficiency, bias, scale and costs, and their relative advantages for different applications are unclear. In the present study, we generated benchmark datasets to systematically evaluate protocols in terms of their power to comprehensively describe cell types and states. We performed a multicenter study comparing 13 commonly used scRNA-seq and single-nucleus RNA-seq protocols applied to a heterogeneous reference sample resource. Comparative analysis revealed marked differences in protocol performance. The protocols differed in library complexity and their ability to detect cell-type markers, impacting their predictive value and suitability for integration into reference cell atlases. These results provide guidance both for individual researchers and for consortium projects such as the Human Cell Atlas.

Single-Cell Dosing and mRNA Sequencing of Suspension and Adherent Cells Using the PolarisTM System.

Single-cell functional analysis provides a natural next step in the now widely adopted single-cell mRNA sequencing studies. Functional studies can be designed to study cellular context by using single-cell culture, perturbation, manipulation, or treatment. Here we present a method for a functional study of 48 single cells by single-cell isolation, dosing, and mRNA sequencing with an integrated fluidic circuit (IFC) on the Fluidigm® Polaris™ system. The major procedures required to execute this protocol are (1) cell preparation and staining; (2) priming, single-cell selection, cell dosing, cell staining, and cDNA generation on the Polaris IFC; and (3) preparation and sequencing of single-cell mRNA-seq libraries. The cell preparation and staining steps employ the use of a universal tracking dye to trace all cells that enter the IFC, while additional fluorescence dyes chosen by the user can be used to differentiate cell types in the overall mix. The steps on the Polaris IFC follow standard protocols, which are also described in the Fluidigm user documentation. The library preparation step adds Illumina® Nextera® XT indexes to the cDNA generated on the Polaris IFC. The resulting sequencing libraries can be sequenced on any Illumina sequencing platform.

Simultaneous Targeted Detection of Proteins and RNAs in Single Cells.

Simultaneous detection of both RNA and protein in individual single cells offers a powerful tool for genotype-to-phenotype investigations. Proximity extension assay (PEA) is a quantitative, sensitive, and multiplex protein detection system that has superb utility in single-cell omic analysis. We implemented PEA using the flexible microfluidic workflow of the Fluidigm® C1™ system followed by real-time quantitative polymerase chain reaction (RT-qPCR) on the Fluidigm Biomark™ HD system. With this workflow, targeted quantification of RNAs and proteins within individual cells is readily conducted.

Silencing the Snail-Dependent RNA Splice Regulator ESRP1 Drives Malignant Transformation of Human Pulmonary Epithelial Cells.

Epithelial-to-mesenchymal transition (EMT) is organized in cancer cells by a set of key transcription factors, but the significance of this process is still debated, including in non-small cell lung cancer (NSCLC). Here, we report increased expression of the EMT-inducing transcription factor Snail in premalignant pulmonary lesions, relative to histologically normal pulmonary epithelium. In immortalized human pulmonary epithelial cells and isogenic derivatives, we documented Snail-dependent anchorage-independent growth in vitro and primary tumor growth and metastatic behavior in vivo Snail-mediated transformation relied upon silencing of the tumor-suppressive RNA splicing regulatory protein ESRP1. In clinical specimens of NSCLC, ESRP1 loss was documented in Snail-expressing premalignant pulmonary lesions. Mechanistic investigations showed that Snail drives malignant progression in an ALDH+CD44+CD24- pulmonary stem cell subset in which ESRP1 and stemness-repressing microRNAs are inhibited. Collectively, our results show how ESRP1 loss is a critical event in lung carcinogenesis, and they identify new candidate directions for targeted therapy of NSCLC.Significance: This study defines a Snail-ESRP1 cancer axis that is crucial for human lung carcinogenesis, with implications for new intervention strategies and translational opportunities. Cancer Res; 78(8); 1986-99. ©2018 AACR.

Single-cell protein-mRNA correlation analysis enabled by multiplexed dual-analyte co-detection.

We have investigated the correlation between proteins and mRNAs in single cells employing an integrated workflow for dual-analyte co-detection. This is achieved by combining the oligo extension reaction (OER), which converts protein levels to DNA levels, with reverse transcription for mRNA detection. Unsupervised gene expression profiling analysis, including principal component analysis and hierarchical clustering, revealed different aspects of the protein-mRNA relationship. Violin plot analysis showed that some genes exhibited similar distribution patterns for proteins and mRNAs. We also demonstrate that cells can be separated into subpopulations based on their protein-mRNA expression profiles, and that different subpopulations have distinct correlation coefficient values. Our results demonstrated that integrated investigations of mRNA and protein levels in single cells allows comprehensive analysis not attainable at bulk levels.

Corrigendum: Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

[This corrects the article on p. 70 in vol. 4, PMID: 27709111.].

Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis.

The study of single cells has evolved over the past several years to include expression and genomic analysis of an increasing number of single cells. Several studies have demonstrated wide spread variation and heterogeneity within cell populations of similar phenotype. While the characterization of these populations will likely set the foundation for our understanding of genomic- and expression-based diversity, it will not be able to link the functional differences of a single cell to its underlying genomic structure and activity. Currently, it is difficult to perturb single cells in a controlled environment, monitor and measure the response due to perturbation, and link these response measurements to downstream genomic and transcriptomic analysis. In order to address this challenge, we developed a platform to integrate and miniaturize many of the experimental steps required to study single-cell function. The heart of this platform is an elastomer-based integrated fluidic circuit that uses fluidic logic to select and sequester specific single cells based on a phenotypic trait for downstream experimentation. Experiments with sequestered cells that have been performed include on-chip culture, exposure to various stimulants, and post-exposure image-based response analysis, followed by preparation of the mRNA transcriptome for massively parallel sequencing analysis. The flexible system embodies experimental design and execution that enable routine functional studies of single cells.

Simple Method To Prepare Oligonucleotide-Conjugated Antibodies and Its Application in Multiplex Protein Detection in Single Cells.

The diversity of nucleic acid sequences enables genomics studies in a highly multiplexed format. Since multiplex protein detection is still a challenge, it would be useful to use genomics tools for this purpose. This can be accomplished by conjugating specific oligonucleotides to antibodies. Upon binding of the oligonucleotide-conjugated antibodies to their targets, the protein levels can be converted to oligonucleotide levels. In this report we describe a simple method for preparing oligonucleotide-conjugated antibodies and discuss this method's application in oligonucleotide extension reaction (OER) for multiplex protein detection. Conjugation is based on strain-promoted alkyne-azide cycloaddition (the Cu-free click reaction), in which the antibody is activated with a dibenzocyclooctyne (DBCO) moiety and subsequently linked covalently with an azide-modified oligonucleotide. In the functional test, the reaction conditions and purification processes were optimized to achieve maximum yield and best performance. The OER assay employs a pair of antibody binders (two antibodies, each conjugated with its own oligonucleotide) developed for each protein target. The two oligonucleotides contain unique six-base complementary regions at their 3' prime ends to allow annealing and extension by DNA synthesis enzymes to form a DNA template. Following preamplification, the DNA template is detected by qPCR. Distinct oligonucleotide sequences are assigned to different antibody binders to enable multiplex protein detection. When tested using recombinant proteins, some antibody binders, such as those specific to CSTB, MET, EpCAM, and CASP3, had dynamic ranges of 5-6 logs. The antibody binders were also used in a multiplexed format in OER assays, and the binders successfully detected their protein targets in cell lysates, and in single cells in combination with the C1 system. This click reaction-based antibody conjugation procedure is cost-effective, needs minimal hands-on time, and is well-suited for the development of affordable multiplex protein assays, which provides the potential to accelerate proteomics research.

Molecular Pathways: Targeting Cellular Energy Metabolism in Cancer via Inhibition of SLC2A1 and LDHA.

Reprogramming of cellular energy metabolism is widely accepted to be one of the main hallmarks of cancer. The aberrant expression pattern of key regulators in the glycolysis pathway in cancer cells corroborates with the hypothesis that most cancer cells utilize aerobic glycolysis as their main ATP production method instead of mitochondrial oxidative phosphorylation. Overexpression of SLC2A1 and LDHA, both important regulators of the glycolysis pathway, was detected in the premalignant lesions and tumors of lung cancer patients, suggesting the involvement of these proteins in early carcinogenesis and tumor progression in cancer. Preclinical studies demonstrated that inhibiting SLC2A1 or LDHA led to diminished tumor growth in vitro and in vivo. SLC2A1 and LDHA inhibitors, when administered in combination with other chemotherapeutic agents, showed synergistic antitumor effects by resensitizing chemoresistant cancer cells to the chemotherapies. These results indicate that disrupting SLC2A1, LDHA, or other regulators in cancer cell energetics is a very promising approach for new targeted therapies.

Molecular profiling of premalignant lesions in lung squamous cell carcinomas identifies mechanisms involved in stepwise carcinogenesis.

Lung squamous cell carcinoma (SCC) is thought to arise from premalignant lesions in the airway epithelium; therefore, studying these lesions is critical for understanding lung carcinogenesis. Previous microarray and sequencing studies designed to discover early biomarkers and therapeutic targets for lung SCC had limited success identifying key driver events in lung carcinogenesis, mostly due to the cellular heterogeneity of patient samples examined and the interindividual variability associated with difficult to obtain airway premalignant lesions and appropriate normal control samples within the same patient. We performed RNA sequencing on laser-microdissected representative cell populations along the SCC pathologic continuum of patient-matched normal basal cells, premalignant lesions, and tumor cells. We discovered transcriptomic changes and identified genomic pathways altered with initiation and progression of SCC within individual patients. We used immunofluorescent staining to confirm gene expression changes in premalignant lesions and tumor cells, including increased expression of SLC2A1, CEACAM5, and PTBP3 at the protein level and increased activation of MYC via nuclear translocation. Cytoband enrichment analysis revealed coordinated loss and gain of expression in chromosome 3p and 3q regions, respectively, during carcinogenesis. This is the first gene expression profiling study of airway premalignant lesions with patient-matched SCC tumor samples. Our results provide much needed information about the biology of premalignant lesions and the molecular changes that occur during stepwise carcinogenesis of SCC, and it highlights a novel approach for identifying some of the earliest molecular changes associated with initiation and progression of lung carcinogenesis within individual patients.

Potential genetic anticipation in hereditary leiomyomatosis-renal cell cancer (HLRCC).

Hereditary leiomyomatosis-renal cell cancer (HLRCC) is an autosomal dominant disorder characterised by cutaneous leiomyomas, symptomatic uterine leiomyomas and aggressive type II papillary renal cell carcinoma. It is caused by heterozygous mutations in the fumarate hydratase (FH) gene on chromosome 1q43. We present evidence of genetic anticipation in HLRCC syndrome. A comprehensive literature review was performed to determine the potential for genetic anticipation in HLRCC syndrome. The normal random effects model was used to evaluate for genetic anticipation to ensure reduction in bias. A total of 11 FH kindreds with available multi-generational data were identified for analysis. The mean difference in age at diagnosis of RCC between the first and second generation was -18.6 years (95 % CI -26.6 to -10.6, p < 0.001). The mean difference in age at diagnosis of RCC between the first and third generation was -36.2 years (95 % CI -47.0 to -25.4, p < 0.001). No evidence of anticipation for uterine leiomyomas was observed (p = 0.349). We report preliminary evidence of genetic anticipation of RCC in HLRCC syndrome. Additional clinical validation is important to confirm this observation, which may have practical implications on counseling and timing of surveillance initiation. Exploration of the underlying mechanisms of anticipation in HLRCC would be of considerable biological interest.

Aldehyde dehydrogenase activity enriches for proximal airway basal stem cells and promotes their proliferation.

Both basal and submucosal gland (SMG) duct stem cells of the airway epithelium are capable of sphere formation in the in vitro sphere assay, although the efficiency at which this occurs is very low. We sought to improve this efficiency of sphere formation by identifying subpopulations of airway basal stem cells (ABSC) and SMG duct cells based on their aldehyde dehydrogenase (ALDH) activity. ALDH(hi) ABSCs and SMG duct cells were highly enriched for the population of cells that could make spheres, while the co-culture of ALDH(hi) differentiated cells with the ALDH(hi) ABSCs increased their sphere-forming efficiency. Specific ALDH agonists and antagonists were used to show that airway specific ALDH isozymes are important for ABSC proliferation. Pathway analysis of gene expression profiling of ALDH(hi) and ALDH(lo) ABSCs revealed a significant upregulation of the arachidonic acid (AA) metabolism pathway in ALDH(hi) ABSCs. We confirmed the importance of this pathway in the metabolism of proliferating ALDH(hi) ABSCs using bioenergetics studies as well as agonists and antagonists of the AA pathway. These studies could lead to the development of novel strategies for altering ABSC proliferation in the airway epithelium.

Isolation and in vitro characterization of basal and submucosal gland duct stem/progenitor cells from human proximal airways.

Basal cells and submucosal gland (SMG) duct cells have been isolated and shown to be stem/progenitor cell populations for the murine airway epithelium. However, methods for the isolation of basal and SMG duct cells from human airways have not been defined. We used an optimized two-step enzyme digestion protocol to strip the surface epithelium from tracheal specimens separate from SMG cells, and we then sorted the basal and duct stem/progenitors using fluorescence-activated cell sorting. We used nerve growth factor receptor, as well as a combination of CD166 and CD44, to sort basal cells and also used CD166 to isolate SMG duct cells. Sorted stem/progenitor cells were cultured to characterize their self-renewal and differentiation ability. Both basal and SMG duct cells grew into spheres. Immunostaining of the spheres showed mostly dense spheres with little to no central lumen. The spheres expressed cytokeratins 5 and 14, with some mucus- and serous-secreting cells. The sphere-forming efficiency and the rate of growth of the spheres varied widely between patient samples and correlated with the degree of hyperplasia of the epithelium. We found that only aldehyde dehydrogenase (ALDH)(hi) basal and duct cells were capable of sphere formation. Global inhibition of ALDH, as well as specific inhibition of the ALDH2 isoform, inhibited self-renewal of both basal and duct cells, thereby producing fewer and smaller spheres. In conclusion, we have developed methods to isolate basal and SMG duct cells from the surface epithelium and SMGs of human tracheas and have developed an in vitro model to characterize their self-renewal and differentiation.

Changes in specialized blood vessels in lymph nodes and their role in cancer metastasis.

BACKGROUND: High endothelial venules (HEV) have been recognized to play a role in metastasis by its changes seen in the cancer microenvironment of lymph nodes (LN) and solid cancers. Squamous cell carcinoma (SCC) of the tongue is a prevalent tumor of the head and neck region with high propensity for LN metastasis. The extent of LN metastasis is the most reliable adverse prognostic factor. Primary tumors can induce vasculature reorganization within sentinel LN before the arrival of tumor cells and HEV represents these remodelled vessels. This study aims to evaluate the cancer induced vascular changes in regional lymph nodes (LN) of patients by studying the morphological and functional alterations of HEV and its correlation with clinical outcome and pathological features. METHODS: This study was based on 65 patients with SCC tongue who underwent primary surgical treatment including neck dissection. The patients were categorized into 2 groups based on the presence of malignancy in their cervical lymph nodes. A review of the patients' pathological and clinical data was performed from a prospective database. Immunohistochemical staining of the tissue blocks for HEV and high-power-field image analysis were performed and analyzed with correlation to the patients' clinical and pathological features. RESULTS: The total number of HEV was found to be significantly associated to disease-free interval. There was a similar association comparing the HEV parameters to overall survival. The density of abnormal HEV was significantly higher in patients with established metastases in their lymph nodes and HEV was shown to be a better prognosis factor than conventional tumor staging. The HEV morphological metamorphosis demonstrates a spectrum that correlates well with disease progression and clinical outcome. CONCLUSIONS: The results suggest that the HEV displays a spectrum of morphological changes in the presence of cancer and LN metastasis, and that HEV is possibly involved in the process of cancer metastasis. We revealed the relationship of HEV and their metamorphosis in pre-metastatic and metastatic environments in regional lymph nodes of tongue cancer patients in relation to clinical outcomes. The significant observation of modified dilated HEV containing red blood cells in lymph nodal basin of a cancer suggests the shifting of its function from one primarily of immune response to that of a blood carrying vessel. It also demonstrated potential prognostic value. More studies are needed to elucidate its potential role in cancer immunotherapy and as a potential novel therapeutic approach to preventing metastasis by manipulating the remodelling processes of HEV.

Identification of an interleukin 13-induced epigenetic signature in allergic airway inflammation.

Epigenetic changes have been implicated in the pathogenesis of asthma. We sought to determine if IL13, a key cytokine in airway inflammation and remodeling, induced epigenetic DNA methylation and miRNAs expression changes in the airways in conjunction with its transcriptional gene regulation. Inducible expression of an IL13 transgene in the airways resulted in significant changes in DNA methylation in 177 genes, most of which were associated with the IL13 transcriptional signature in the airways. A large number of genes whose expression was induced by IL13 were found to have decreased methylation, including those involved in tissue remodeling (Olr1), leukocyte influx (Cxcl3, Cxcl5, CSFr2b), and the Th2 response (C3ar1, Chi3l4). Reciprocally, some genes whose expression was suppressed were found to have increased methylation (e.g. Itga8). In addition, miRNAs were identified with targets for lung development and Wnt signaling, amongst others. These results indicate that IL13 confers an epigenetic methylation and miRNA signature that accompanies its transcriptional program in the airways, which may play a critical role in airway inflammation and remodeling.

2011 Young Surgeon's Award Winner: high endothelial venules: a novel prognostic marker in cancer metastasis and the missing link?

INTRODUCTION: The extent of lymph nodes (LNs) metastasis is a major determinant for the staging and the most reliable adverse prognostic factor. Primary tumours can induce lymphatics and vasculature reorganisations within sentinel LN before the arrival of cancer cells and these key blood vessels are identified as high endothelial venules (HEV). The alterations of HEV in the presence of cancer, coupled with the increased proliferation rate of the endothelial cells, results in a functional shift of HEV from immune response mediator to blood flow carrier. We aim to evaluate tumour-induced vascularisation in regional LN of cancer patients by studying the morphological and functional alterations of HEV and its correlation to clinico-pathological features. MATERIALS AND METHODS: This multi-centre study with a prospective database identified 65 consecutive patients with tongue squamous cell carcinoma (SCC) who underwent primary surgical treatment from 2001 to 2005. Immunohistochemical staining for HEV and image analysis were performed and analysed with correlation to the patients' clinico-pathological features. RESULTS: The total number of HEV is significantly associated to disease-free interval when controlling for the group (P = 0.022) as well as combining both groups as one cohort (P = 0.023). There is also a similar association comparing the HEV parameters to overall survival. CONCLUSION: Our results suggest that HEV possibly plays a key role in the pathogenesis of lymphatic and subsequent distant metastases and may provide the missing link in cancer metastasis. Confirmation of this hypothesis would offer a novel therapeutic approach to preventing metastasis by blocking the remodeling processes of HEV in LN.