High-throughput single-cell DNA sequencing of acute myeloid leukemia tumors with droplet microfluidics.

To enable the characterization of genetic heterogeneity in tumor cell populations, we developed a novel microfluidic approach that barcodes amplified genomic DNA from thousands of individual cancer cells confined to droplets. The barcodes are then used to reassemble the genetic profiles of cells from next-generation sequencing data. By using this approach, we sequenced longitudinally collected acute myeloid leukemia (AML) tumor populations from two patients and genotyped up to 62 disease relevant loci across more than 16,000 individual cells. Targeted single-cell sequencing was able to sensitively identify cells harboring pathogenic mutations during complete remission and uncovered complex clonal evolution within AML tumors that was not observable with bulk sequencing. We anticipate that this approach will make feasible the routine analysis of AML heterogeneity, leading to improved stratification and therapy selection for the disease.

Microfluidic droplet enrichment for targeted sequencing.

Targeted sequence enrichment enables better identification of genetic variation by providing increased sequencing coverage for genomic regions of interest. Here, we report the development of a new target enrichment technology that is highly differentiated from other approaches currently in use. Our method, MESA (Microfluidic droplet Enrichment for Sequence Analysis), isolates genomic DNA fragments in microfluidic droplets and performs TaqMan PCR reactions to identify droplets containing a desired target sequence. The TaqMan positive droplets are subsequently recovered via dielectrophoretic sorting, and the TaqMan amplicons are removed enzymatically prior to sequencing. We demonstrated the utility of this approach by generating an average 31.6-fold sequence enrichment across 250 kb of targeted genomic DNA from five unique genomic loci. Significantly, this enrichment enabled a more comprehensive identification of genetic polymorphisms within the targeted loci. MESA requires low amounts of input DNA, minimal prior locus sequence information and enriches the target region without PCR bias or artifacts. These features make it well suited for the study of genetic variation in a number of research and diagnostic applications.

RNA-Seq following PCR-based sorting reveals rare cell transcriptional signatures.

BACKGROUND: Rare cell subtypes can profoundly impact the course of human health and disease, yet their presence within a sample is often missed with bulk molecular analysis. Single-cell analysis tools such as FACS, FISH-FC and single-cell barcode-based sequencing can investigate cellular heterogeneity; however, they have significant limitations that impede their ability to identify and transcriptionally characterize many rare cell subpopulations. RESULTS: PCR-activated cell sorting (PACS) is a novel cytometry method that uses single-cell TaqMan PCR reactions performed in microfluidic droplets to identify and isolate cell subtypes with high-throughput. Here, we extend this method and demonstrate that PACS enables high-dimensional molecular profiling on TaqMan-targeted cells. Using a random priming RNA-Seq strategy, we obtained high-fidelity transcriptome measurements following PACS sorting of prostate cancer cells from a heterogeneous population. The sequencing data revealed prostate cancer gene expression profiles that were obscured in the unsorted populations. Single-cell expression analysis with PACS was subsequently used to confirm a number of the differentially expressed genes identified with RNA sequencing. CONCLUSIONS: PACS requires minimal sample processing, uses readily available TaqMan assays and can isolate cell subtypes with high sensitivity. We have now validated a method for performing next-generation sequencing on mRNA obtained from PACS isolated cells. This capability makes PACS well suited for transcriptional profiling of rare cells from complex populations to obtain maximal biological insight into cell states and behaviors.

Identification and genetic analysis of cancer cells with PCR-activated cell sorting.

Cell sorting is a central tool in life science research for analyzing cellular heterogeneity or enriching rare cells out of large populations. Although methods like FACS and FISH-FC can characterize and isolate cells from heterogeneous populations, they are limited by their reliance on antibodies, or the requirement to chemically fix cells. We introduce a new cell sorting technology that robustly sorts based on sequence-specific analysis of cellular nucleic acids. Our approach, PCR-activated cell sorting (PACS), uses TaqMan PCR to detect nucleic acids within single cells and trigger their sorting. With this method, we identified and sorted prostate cancer cells from a heterogeneous population by performing >132 000 simultaneous single-cell TaqMan RT-PCR reactions targeting vimentin mRNA. Following vimentin-positive droplet sorting and downstream analysis of recovered nucleic acids, we found that cancer-specific genomes and transcripts were significantly enriched. Additionally, we demonstrate that PACS can be used to sort and enrich cells via TaqMan PCR reactions targeting single-copy genomic DNA. PACS provides a general new technical capability that expands the application space of cell sorting by enabling sorting based on cellular information not amenable to existing approaches.

Scrib regulates HGF-mediated epithelial morphogenesis and is stabilized by Sgt1-HSP90.

Scribble was originally identified as a Drosophila protein that regulates epithelial polarity and formation of the basolateral surface. The mammalian orthologue, Scrib, is evolutionarily conserved, but does not appear to be necessary for apical-basolateral epithelial polarity. Instead, it is implicated in the regulation of cell survival, protein trafficking, adhesion and migration. A key issue is to understand the molecular pathway by which Scrib participates in these processes. We have investigated Scrib using a three-dimensional epithelial cell culture system. We show a novel association between the leucine-rich repeat domain of Scrib and the co-chaperone Sgt1 and demonstrate that these proteins are necessary for epithelial morphogenesis and tubulogenesis following hepatocyte growth factor (HGF) stimulation. The molecular chaperone HSP90 is also required for Sgt1 association with Scrib, and both Sgt1 and HSP90 are needed to ensure proper Scrib protein levels. Furthermore, reduced Scrib stability, following inhibition of Sgt1-HSP90, lowers the cellular abundance of the Scrib-ßPix-PAK complex. Inhibition of any member of this complex, Scrib, ßPix or PAK, is sufficient to block HGF-mediated epithelial morphogenesis. The identification of Scrib as an Sgt1-HSP90 client protein required for three-dimensional cell migration suggests that chaperone-mediated regulation of polarity protein stability and homeostasis is an unappreciated mechanism underlying dynamic rearrangements during morphogenesis.

p120 catenin is required for normal renal tubulogenesis and glomerulogenesis.

Defects in the development or maintenance of tubule diameter correlate with polycystic kidney disease. Here, we report that absence of the cadherin regulator p120 catenin (p120ctn) from the renal mesenchyme prior to tubule formation leads to decreased cadherin levels with abnormal morphologies of early tubule structures and developing glomeruli. In addition, mutant mice develop cystic kidney disease, with markedly increased tubule diameter and cellular proliferation, and detached luminal cells only in proximal tubules. The p120ctn homolog Arvcf is specifically absent from embryonic proximal tubules, consistent with the specificity of the proximal tubular phenotype. p120ctn knockdown in renal epithelial cells in 3D culture results in a similar cystic phenotype with reduced levels of E-cadherin and active RhoA. We find that E-cadherin knockdown, but not RhoA inhibition, phenocopies p120ctn knockdown. Taken together, our data show that p120ctn is required for early tubule and glomerular morphogenesis, as well as control of luminal diameter, probably through regulation of cadherins.

High-throughput gene expression analysis with TempO-LINC sensitively resolves complex brain, lung and kidney heterogeneity at single-cell resolution.

We report the development and performance of a novel genomics platform, TempO-LINC, for conducting high-throughput transcriptomic analysis on single cells and nuclei. TempO-LINC works by adding cell-identifying molecular barcodes onto highly selective and high-sensitivity gene expression probes within fixed cells, without having to first generate cDNA. Using an instrument-free combinatorial-indexing approach, all probes within the same fixed cell receive an identical barcode, enabling the reconstruction of single-cell gene expression profiles across as few as several hundred cells and up to 100,000+ cells per run. The TempO-LINC approach is easily scalable based on the number of barcodes and rounds of barcoding performed; however, for the experiments reported in this study, the assay utilized over 5.3 million unique barcodes. TempO-LINC has a robust protocol for fixing and banking cells and displays high-sensitivity gene detection from multiple diverse sample types. We show that TempO-LINC has an observed multiplet rate of less than 1.1% and a cell capture rate of ~50%. Although the assay can accurately profile the whole transcriptome (19,683 human or 21,400 mouse genes), it can be targeted to measure only actionable/informative genes and molecular pathways of interest - thereby reducing sequencing requirements. In this study, we applied TempO-LINC to profile the transcriptomes of 89,722 cells across multiple sample types, including nuclei from mouse lung, kidney and brain tissues. The data demonstrated the ability to identify and annotate at least 50 unique cell populations and positively correlate expression of cell type-specific molecular markers within them. TempO-LINC is a robust new single-cell technology that is ideal for large-scale applications/studies across thousands of samples with high data quality.

Ultrahigh-throughput Mammalian single-cell reverse-transcriptase polymerase chain reaction in microfluidic drops.

The behaviors of complex biological systems are often dictated by the properties of their heterogeneous and sometimes rare cellular constituents. Correspondingly, the analysis of individual cells from a heterogeneous population can reveal information not obtainable by ensemble measurements. Reverse-transcriptase polymerase chain reaction (RT-PCR) is a widely used method that enables transcriptional profiling and sequencing analysis on bulk populations of cells. Major barriers to successfully implementing this technique for mammalian single-cell studies are the labor, cost, and low-throughput associated with current approaches. In this report, we describe a novel droplet-based microfluidic system for performing ~50000 single-cell RT-PCR reactions in a single experiment while consuming a minimal amount of reagent. Using cell type-specific staining and TaqMan RT-PCR probes, we demonstrate the identification of specific cells from a mixed human cell population. The throughput, robust detection rate and specificity of this method makes it well-suited for characterizing large, heterogeneous populations of cells at the transcriptional level.

Involvement of RhoA, ROCK I and myosin II in inverted orientation of epithelial polarity.

In multicellular epithelial tissues, the orientation of polarity of each cell must be coordinated. Previously, we reported that for Madin-Darby canine kidney cells in three-dimensional collagen gel culture, blockade of beta1-integrin by the AIIB2 antibody or expression of dominant-negative Rac1N17 led to an inversion of polarity, such that the apical surfaces of the cells were misorientated towards the extracellular matrix. Here, we show that this process results from the activation of RhoA. Knockdown of RhoA by short hairpin RNA reverses the inverted orientation of polarity, resulting in normal cysts. Inhibition of RhoA downstream effectors, Rho kinase (ROCK I) and myosin II, has similar effects. We conclude that the RhoA-ROCK I-myosin II pathway controls the inversion of orientation of epithelial polarity caused by AIIB2 or Rac1N17. These results might be relevant to the hyperactivation of RhoA and disruption of normal polarity frequently observed in human epithelial cancers.

Clonal Evolution and Changes in Two AML Patients Detected with A Novel Single-Cell DNA Sequencing Platform.

Next-generation sequencing (NGS) is used to detect gene variants in genetically complex cell populations of cancer patient samples. Traditional bulk analysis can only provide average variant allele frequencies of the targeted genes across all sampled cells. It fails to resolve mutational co-occurrences and may miss rare cancer cells. Genome analysis at the single cell level offers the opportunity to more fully resolve clonal architecture. Peripheral blood mononuclear cells were sampled from acute myeloid leukemia patients longitudinally and single-cell DNA sequencing libraries were generated with a novel droplet-based microfluidics approach. Molecular profiling of single nucleotide variants across thousands of cells revealed genetic chimerism in patients after bone marrow transplantation (BMT). Importantly, hierarchical clustering analysis of single nucleotide variants (SNVs) uncovered a distinct oncogenic clone of cells carrying mutated tumor-suppressor and/or oncogene(s). This novel single-cell DNA sequencing approach enabled precise monitoring of engraftment and revealed clonal evolution of oncogenic cells during the progression and treatment of the disease.

Clonal Selection with RAS Pathway Activation Mediates Secondary Clinical Resistance to Selective FLT3 Inhibition in Acute Myeloid Leukemia.

Gilteritinib is a potent and selective FLT3 kinase inhibitor with single-agent clinical efficacy in relapsed/refractory FLT3-mutated acute myeloid leukemia (AML). In this context, however, gilteritinib is not curative, and response duration is limited by the development of secondary resistance. To evaluate resistance mechanisms, we analyzed baseline and progression samples from patients treated on clinical trials of gilteritinib. Targeted next-generation sequencing at the time of AML progression on gilteritinib identified treatment-emergent mutations that activate RAS/MAPK pathway signaling, most commonly in NRAS or KRAS. Less frequently, secondary FLT3-F691L gatekeeper mutations or BCR-ABL1 fusions were identified at progression. Single-cell targeted DNA sequencing revealed diverse patterns of clonal selection and evolution in response to FLT3 inhibition, including the emergence of RAS mutations in FLT3-mutated subclones, the expansion of alternative wild-type FLT3 subclones, or both patterns simultaneously. These data illustrate dynamic and complex changes in clonal architecture underlying response and resistance to mutation-selective tyrosine kinase inhibitor therapy in AML. SIGNIFICANCE: Comprehensive serial genotyping of AML specimens from patients treated with the selective FLT3 inhibitor gilteritinib demonstrates that complex, heterogeneous patterns of clonal selection and evolution mediate clinical resistance to tyrosine kinase inhibition in FLT3-mutated AML. Our data support the development of combinatorial targeted therapeutic approaches for advanced AML.See related commentary by Wei and Roberts, p. 998.This article is highlighted in the In This Issue feature, p. 983.

A gain-of-function allele of cbp-1, the Caenorhabditis elegans ortholog of the mammalian CBP/p300 gene, causes an increase in histone acetyltransferase activity and antagonism of activated Ras.

An RTK-Ras-mitogen-activated protein kinase (MAPK) signaling pathway plays a key role in vulval induction in Caenorhabditis elegans. We have previously carried out screens for suppressors of activated Ras to identify factors that play critical roles in the regulation of the pathway. ku258 was isolated as a semidominant allele that suppresses the Multivulva phenotype caused by activated let-60 ras. Our genetic and molecular analyses indicate that ku258 is a gain-of-function allele resulting from two point mutations in the C. elegans homolog of the transcriptional coactivator p300/CBP, cbp-1. Genetic data also suggest that cbp-1 may act downstream of the Ras signaling pathway, but not primarily downstream of the Wnt signaling pathway, to negatively regulate vulval cell fate specification. cbp-1 may function in concert with LIN-1, an Ets transcription factor family member that is one of the targets of MAPK. In vitro histone acetylation assays have revealed that together, the two point mutations cause a sevenfold increase in the histone acetyltransferase (HAT) activity of recombinant CBP-1. To our knowledge, this is the only such HAT activity mutation isolated in a CBP/p300 family protein, and this mutation may define a negative role of the HAT activity in antagonizing Ras function in a specific developmental event.

The Caenorhabditis elegans nuclear receptor gene nhr-25 regulates epidermal cell development.

The development of the epidermis of Caenorhabditis elegans involves cell fusion, migration, and differentiation events. To understand the mechanisms underlying these processes, we characterized the roles of NHR-25, a member of the nuclear receptor family of transcription factors. The NHR-25 homologs Ftz-F1 in Drosophila and SF-1 in mammals are involved in various biological processes, including regulation of patterning during development, reproduction, metabolism, metamorphosis, and homeostasis. Impairment of nhr-25 activity leads to severe phenotypes in embryos and many postembryonic tissues. Further analysis has indicated that nhr-25 activity is required for the proper development, including cell-cell fusion, of several epidermal cell types, such as the epidermal syncytial, seam, and Pn.p cells. Our results also suggest that nhr-25 is likely to regulate cell-cell junctions and/or fusion. In a subset of Pn.p cells, called vulval precursor cells, nhr-25 acts collaboratively with the lin-39 Hox gene in regulating vulval cell differentiation. Additionally, our data suggest that nhr-25 may also function with another Hox gene, nob-1, during embryogenesis. Overall, our results indicate that nhr-25 plays an integral role in regulating cellular processes of epidermal cells.

Lgr6 is a stem cell marker in mouse skin squamous cell carcinoma.

The G-protein-coupled receptors LGR4, LGR5 and LGR6 are Wnt signaling mediators, but their functions in squamous cell carcinoma (SCC) are unclear. Using lineage tracing in Lgr5-EGFP-CreERT2/Rosa26-Tomato and Lgr6-EGFP-CreERT2/Rosa26-Tomato reporter mice, we demonstrate that Lgr6, but not Lgr5, acts as an epithelial stem cell marker in SCCs in vivo. We identify, by single-molecule in situ hybridization and cell sorting, rare cells positive for Lgr6 expression in immortalized keratinocytes and show that their frequency increases in advanced SCCs. Lgr6 expression is enriched in cells with stem cell characteristics, and Lgr6 downregulation in vivo causes increased epidermal proliferation with expanded lineage tracing from epidermal stem cells positive for Lgr6 expression. Surprisingly, mice with germline knockout of Lgr6 are predisposed to SCC development, through a mechanism that includes compensatory upregulation of Lgr5. These data provide a model for human patients with germline loss-of-function mutations in Wnt pathway genes, including RSPO1 or LGR4, who show increased susceptibility to squamous tumor development.

A molecular switch for the orientation of epithelial cell polarization.

The formation of epithelial tissues containing lumens requires not only the apical-basolateral polarization of cells, but also the coordinated orientation of this polarity such that the apical surfaces of neighboring cells all point toward the central lumen. Defects in extracellular matrix (ECM) signaling lead to inverted polarity so that the apical surfaces face the surrounding ECM. We report a molecular switch mechanism controlling polarity orientation. ECM signals through a ß1-integrin/FAK/p190RhoGAP complex to downregulate a RhoA/ROCK/Ezrin pathway at the ECM interface. PKCßII phosphorylates the apical identity-promoting Podocalyxin/NHERF1/Ezrin complex, removing Podocalyxin from the ECM-abutting cell surface and initiating its transcytosis to an apical membrane initiation site for lumen formation. Inhibition of this switch mechanism results in the retention of Podocalyxin at the ECM interface and the development instead of collective front-rear polarization and motility. Thus, ECM-derived signals control the morphogenesis of epithelial tissues by controlling the collective orientation of epithelial polarization.

Picoinjection enables digital detection of RNA with droplet rt-PCR.

The ability to add reagents to drops in a sequential fashion is necessary for numerous applications of microfluidics in biology. An important method for accomplishing this is picoinjection, a technique in which reagents are injected into aqueous drops using an electric field. While picoinjection has been shown to allow the precise addition of reagents to drops, its compatibility with biological reactions is yet to be thoroughly demonstrated. Here, we investigate the compatibility of picoinjection with digital RT-PCR Taqman assays, reactions that incorporate nucleic acids, enzymes, and other common biological reagents. We find that picoinjection is compatible with this assay and enables the detection of RNA transcripts at rates comparable to workflows not incorporating picoinjection. We also find that picoinjection results in negligible transfer of material between drops and that the drops faithfully retain their compartmentalization.

Electrode-free picoinjection of microfluidic drops.

Existing methods for injecting reagents into drops utilize metal electrodes integrated into the microfluidic chip, adding cumbersome and error-prone steps to the fabrication process. Here, we introduce a method that uses the injection fluid itself as an electrode by exploiting dissolved electrolytes in the solution. This obviates the need for metal electrodes and reduces fabrication time and complexity while also affording precise control of the injected reagent volumes.

Phosphoinositides in cell architecture.

Inositol phospholipids have been implicated in almost all aspects of cellular physiology including spatiotemporal regulation of cellular signaling, acquisition of cellular polarity, specification of membrane identity, cytoskeletal dynamics, and regulation of cellular adhesion, motility, and cytokinesis. In this review, we examine the critical role phosphoinositides play in these processes to execute the establishment and maintenance of cellular architecture. Epithelial tissues perform essential barrier and transport functions in almost all major organs. Key to their development and function is the establishment of epithelial cell polarity. We place a special emphasis on highlighting recent studies demonstrating phosphoinositide regulation of epithelial cell polarity and how individual cells use phosphoinositides to further organize into epithelial tissues.

The Cdc42 GEF Intersectin 2 controls mitotic spindle orientation to form the lumen during epithelial morphogenesis.

Epithelial organs are made of tubes and cavities lined by a monolayer of polarized cells that enclose the central lumen. Lumen formation is a crucial step in the formation of epithelial organs. The Rho guanosine triphosphatase (GTPase) Cdc42, which is a master regulator of cell polarity, regulates the formation of the central lumen in epithelial morphogenesis. However, how Cdc42 is regulated during this process is still poorly understood. Guanine nucleotide exchange factors (GEFs) control the activation of small GTPases. Using the three-dimensional Madin-Darby canine kidney model, we have identified a Cdc42-specific GEF, Intersectin 2 (ITSN2), which localizes to the centrosomes and regulates Cdc42 activation during epithelial morphogenesis. Silencing of either Cdc42 or ITSN2 disrupts the correct orientation of the mitotic spindle and normal lumen formation, suggesting a direct relationship between these processes. Furthermore, we demonstrated this direct relationship using LGN, a component of the machinery for mitotic spindle positioning, whose disruption also results in lumen formation defects.