<?xml version="1.0" encoding="UTF-8" standalone="no"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink" article-type="research-article">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Bioeng. Biotechnol.</journal-id>
<journal-title>Frontiers in Bioengineering and Biotechnology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Bioeng. Biotechnol.</abbrev-journal-title>
<issn pub-type="epub">2296-4185</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="doi">10.3389/fbioe.2016.00070</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Bioengineering and Biotechnology</subject>
<subj-group>
<subject>Technology Report</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>Fluidic Logic Used in a Systems Approach to Enable Integrated Single-Cell Functional Analysis</article-title>
</title-group>
<contrib-group>
<contrib contrib-type="author">
<name><surname>Ramalingam</surname> <given-names>Naveen</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x02020;</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/335538"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Fowler</surname> <given-names>Brian</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x02020;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Szpankowski</surname> <given-names>Lukasz</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="author-notes" rid="fn001"><sup>&#x02020;</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Leyrat</surname> <given-names>Anne A.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Hukari</surname> <given-names>Kyle</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/352722"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Maung</surname> <given-names>Myo Thu</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/377208"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Yorza</surname> <given-names>Wiganda</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/377216"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Norris</surname> <given-names>Michael</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Cesar</surname> <given-names>Chris</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Shuga</surname> <given-names>Joe</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Gonzales</surname> <given-names>Michael L.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/361660"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Sanada</surname> <given-names>Chad D.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Wang</surname> <given-names>Xiaohui</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Yeung</surname> <given-names>Rudy</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Hwang</surname> <given-names>Win</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Axsom</surname> <given-names>Justin</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Devaraju</surname> <given-names>Naga Sai Gopi Krishna</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Angeles</surname> <given-names>Ninez Delos</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/352777"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Greene</surname> <given-names>Cassandra</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/376384"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Zhou</surname> <given-names>Ming-Fang</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Ong</surname> <given-names>Eng-Seng</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Poh</surname> <given-names>Chang-Chee</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/376389"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Lam</surname> <given-names>Marcos</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Choi</surname> <given-names>Henry</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Htoo</surname> <given-names>Zaw</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Lee</surname> <given-names>Leo</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Chin</surname> <given-names>Chee-Sing</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Shen</surname> <given-names>Zhong-Wei</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Lu</surname> <given-names>Chong T.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Holcomb</surname> <given-names>Ilona</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/378151"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Ooi</surname> <given-names>Aik</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/374246"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Stolarczyk</surname> <given-names>Craig</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author">
<name><surname>Shuga</surname> <given-names>Tony</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/369057"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Livak</surname> <given-names>Kenneth J.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/312517"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Larsen</surname> <given-names>Cate</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<uri xlink:href="http://frontiersin.org/people/u/403352"/>
</contrib>
<contrib contrib-type="author">
<name><surname>Unger</surname> <given-names>Marc</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
</contrib>
<contrib contrib-type="author" corresp="yes">
<name><surname>West</surname> <given-names>Jay A. A.</given-names></name>
<xref ref-type="aff" rid="aff1"><sup>1</sup></xref>
<xref ref-type="corresp" rid="cor1">&#x0002A;</xref>
</contrib>
</contrib-group>
<aff id="aff1"><sup>1</sup><institution>New Technologies Research Department, Fluidigm Corporation</institution>, <addr-line>South San Francisco, CA</addr-line>, <country>USA</country></aff>
<author-notes>
<fn fn-type="edited-by"><p>Edited by: Xinghua Pan, Yale University, USA</p></fn>
<fn fn-type="edited-by"><p>Reviewed by: Senentxu Lanceros-Mendez, University of Minho, Portugal; Lin Han, Shandong University, China</p></fn>
<corresp content-type="corresp" id="cor1">&#x0002A;Correspondence: Jay A. A. West, <email>jay.west&#x00040;fluidigm.com</email></corresp>
<fn fn-type="other" id="fn001"><p><sup>&#x02020;</sup>These authors have contributed equally to this work.</p></fn>
<fn fn-type="other" id="fn002"><p>Specialty section: This article was submitted to Genomic Assay Technology, a section of the journal Frontiers in Bioengineering and Biotechnology</p></fn>
</author-notes>
<pub-date pub-type="epub">
<day>21</day>
<month>09</month>
<year>2016</year>
</pub-date>
<pub-date pub-type="collection">
<year>2016</year>
</pub-date>
<volume>4</volume>
<elocation-id>70</elocation-id>
<history>
<date date-type="received">
<day>31</day>
<month>05</month>
<year>2016</year>
</date>
<date date-type="accepted">
<day>23</day>
<month>08</month>
<year>2016</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#x000A9; 2016 Ramalingam, Fowler, Szpankowski, Leyrat, Hukari, Maung, Yorza, Norris, Cesar, Shuga, Gonzales, Sanada, Wang, Yeung, Hwang, Axsom, Devaraju, Angeles, Greene, Zhou, Ong, Poh, Lam, Choi, Htoo, Lee, Chin, Shen, Lu, Holcomb, Ooi, Stolarczyk, Shuga, Livak, Larsen, Unger and West.</copyright-statement>
<copyright-year>2016</copyright-year>
<copyright-holder>Ramalingam, Fowler, Szpankowski, Leyrat, Hukari, Maung, Yorza, Norris, Cesar, Shuga, Gonzales, Sanada, Wang, Yeung, Hwang, Axsom, Devaraju, Angeles, Greene, Zhou, Ong, Poh, Lam, Choi, Htoo, Lee, Chin, Shen, Lu, Holcomb, Ooi, Stolarczyk, Shuga, Livak, Larsen, Unger and West</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/"><p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p></license>
</permissions>
<abstract>
<p>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.</p>
</abstract>
<kwd-group>
<kwd>single-cell</kwd>
<kwd>mRNA-seq</kwd>
<kwd>functional studies</kwd>
<kwd>Fluidigm</kwd>
<kwd>Polaris</kwd>
</kwd-group>
<counts>
<fig-count count="15"/>
<table-count count="1"/>
<equation-count count="0"/>
<ref-count count="27"/>
<page-count count="19"/>
<word-count count="8877"/>
</counts>
</article-meta>
</front>
<body>
<sec id="S1" sec-type="introduction">
<title>Introduction</title>
<p>Recent single-cell transcriptomic analyses have documented the importance of cellular heterogeneity in studying cancer (Ennen et al., <xref ref-type="bibr" rid="B13">2014</xref>; Saadatpour et al., <xref ref-type="bibr" rid="B23">2014</xref>; Kim et al., <xref ref-type="bibr" rid="B18">2015</xref>), immunology (Shalek et al., <xref ref-type="bibr" rid="B24">2014</xref>), developmental biology (Briggs et al., <xref ref-type="bibr" rid="B6">2015</xref>), stem cell research (Wilson et al., <xref ref-type="bibr" rid="B26">2015</xref>), and neurobiology (Pollen et al., <xref ref-type="bibr" rid="B22">2015</xref>). It has been estimated that the human body contains 37.2 trillion cells (Bianconi et al., <xref ref-type="bibr" rid="B5">2013</xref>), excluding the complex microbiome that lives in the human body. High-throughput single-cell mRNA sequencing provides an unbiased path to classifying this vast number of cells into cell types. This endeavor has stimulated the development of methods to increase throughput (Fan et al., <xref ref-type="bibr" rid="B14">2015</xref>; Klein et al., <xref ref-type="bibr" rid="B19">2015</xref>; Macosko et al., <xref ref-type="bibr" rid="B21">2015</xref>). The classification of cell types can be thought of as a high-resolution anatomy. At the single-cell level, moving from anatomy to physiology or from description to mechanism means moving from cell type to cell function. This will require integrating transcriptional data with other cellular measurements. In this regard, progress has been made in obtaining transcriptomic and genomic information (Dey et al., <xref ref-type="bibr" rid="B12">2015</xref>; Macaulay, <xref ref-type="bibr" rid="B20">2015</xref>), transcriptomic and epigenomic information (Angermueller et al., <xref ref-type="bibr" rid="B2">2016</xref>), or transcriptomic and proteomic information (Darmanis et al., <xref ref-type="bibr" rid="B9">2016</xref>; Frei et al., <xref ref-type="bibr" rid="B15">2016</xref>) from the same single cell.</p>
<p>Moving from cell type to cell function will also require understanding how single-cell profiles change in response to perturbations. It is important to examine these effects at the single-cell level because cell-to-cell heterogeneity has been observed in a diverse set of circumstances, such as the response of macrophages to bacterial invasion (Avraham et al., <xref ref-type="bibr" rid="B3">2015</xref>), the response of hematopoietic cells to various drugs (Bendall et al., <xref ref-type="bibr" rid="B4">2011</xref>), and drug resistance in adenocarcinoma cells (Kim et al., <xref ref-type="bibr" rid="B18">2015</xref>). Progress in the long-term culture of circulating tumor cells (Gao et al., <xref ref-type="bibr" rid="B16">2014</xref>; Yu et al., <xref ref-type="bibr" rid="B27">2014</xref>; Cayrefourcq et al., <xref ref-type="bibr" rid="B7">2015</xref>; Alix-Panabi&#x000E8;res et al., <xref ref-type="bibr" rid="B1">2016</xref>) enables single-cell functional studies on this important class of cells, which should lead to improved cancer diagnosis and therapy. Performing perturbation experiments on single cells requires care in maintaining the appropriate microenvironment. Examining the effects of serum on mouse embryonic stem cells (ESCs), researchers (Guo et al., <xref ref-type="bibr" rid="B17">2016</xref>) concluded that &#x0201C;a large proportion of intracellular network variability is due to the extracellular culture environment.&#x0201D; Microfluidic-based approaches are attractive for the precise control of the microenvironment because they enable structures at a size appropriate for single cells. Microfluidic systems for high-throughput preparation of sequencing libraries, though, have cell lysis as the initial step and thus are not suitable to maintain single cells for experimentation. What is required is a system specifically designed to capture, maintain, perturb, and observe single cells and then prepare these cells for high-dimensional analysis.</p>
<p>In this paper, we report development of an integrated fluidic circuit (IFC) that uses fluidic logic to actively select and sequester desired single cells based on particular biological markers of interest. This Polaris&#x02122; IFC can sequester up to 48 single cells. If required, the cells can be cultured in appropriate medium in order to control and manipulate the microenvironment around the sequestered cells. For adherent cells, appropriate extracellular matrix (ECM) can be coated inside the culture chambers. The single cells can be perturbed with a drug or other stimuli (i.e., mRNA, cytokines, bacteria, or viruses), with the response to perturbation monitored and measured by fluorescence imaging. Subsequently, the single cells are processed for cell lysis, reverse transcription (RT), and full-length transcriptome amplification using template-switching chemistry. Following harvest from the IFC, sequencing libraries are generated using a modified Nextera<sup>&#x000AE;</sup> protocol and sequenced on any Illumina<sup>&#x000AE;</sup> platform (Figure <xref ref-type="fig" rid="F1">1</xref>).</p>
<fig id="F1" position="float">
<label>Figure 1</label>
<caption><p><bold>Typical workflow for single-cell functional studies</bold>. The input single-cell suspension can be obtained from blood, primary cells, or cell culture. In the case of rare cells or a subset of cell population, there is an option to enrich them prior to use with Polaris using either fluorescence-assisted cell sorting (FACS) or other methods. Subsequently, the cells are labeled with a universal fluorescent marker for tracking the cells on the Fluidigm<sup>&#x000AE;</sup> Polaris system. Most of the single-cell functional study steps are automated on the Polaris system. The Polaris system generates preamplified full-length cDNA, which can be further processed for library preparation and massive parallel sequencing for mRNA sequencing.</p></caption>
<graphic xlink:href="fbioe-04-00070-g001.tif"/>
</fig>
</sec>
<sec id="S2" sec-type="materials|methods">
<title>Materials and Methods</title>
<sec id="S2-1">
<title>Design and Fabrication of Logic-Based Integrated Fluidic Circuit</title>
<p>The nanoscale IFC consists of a plastic carrier and a polydimethylsiloxane (PDMS) core (Figure <xref ref-type="fig" rid="F2">2</xref>A) or fluidic circuit. The carrier contains reservoir wells for input and output of reagents and circuit control. It provides a platform to facilitate interfacing with the fluidic circuit. The fluidic circuit with the desired microfluidic control components was fabricated using multilayer soft lithography (MSL<sup>&#x000AE;</sup>) process (Unger et al., <xref ref-type="bibr" rid="B25">2000</xref>). Fluidic circuit components include flow and control channels, valves, multiplexors, and logic devices [such as serial-to-parallel shift register (SR)]. Fabrication and operational details of the fluidic logic circuits and devices were reported earlier (Devaraju and Unger, <xref ref-type="bibr" rid="B10">2012</xref>). The IFC is designed to have the capability to actively select single cells based on fluorescent markers, isolate them to a desired holding location (cell capture site), apply individual conditions (feed medium and dose reagents to cells), and finally study the functional response. Execution of all these complex functions in a routine fashion requires flexible, programmable operational control, which in turn requires many controls in a parallel manner. Traditionally, in microfluidics, a dedicated external control line is required to independently control a set of valves. This imposes a limitation on the number of practical on-chip control operations and poses a challenge for scalability by requiring more external hardware. On-chip control architecture capable of receiving and processing data by elementary computation and decision-making can integrate programmability of controls on-chip and allows an increase in the number of on-chip control lines for the same number of external chip connections (Figure <xref ref-type="fig" rid="F2">2</xref>B). We developed such a microfluidic fluidic logic and implemented it on our Polaris IFC.</p>
<fig id="F2" position="float">
<label>Figure 2</label>
<caption><p><bold>(A)</bold> Polaris mRNA-seq dosing integrated fluidic circuit (IFC). The polydimethylsiloxane (PDMS) elastomer is attached to a plastic carrier. In order to image the cells with low-fluorescence background on the Polaris system, the PDMS is backed with a black material. This black backing is removed after cell perturbation and fluorescence response measurement. Removal of black backing enables high-resolution imaging of cells on a fluorescent microscope. The CAD drawing of the microfluidic components is shown on the right. The green channels are the control lines. <bold>(B)</bold> Polaris IFC uses fluidic processor to receive serial control inputs and converts them to parallel shift register elements. Traditional microfluidic control elements use one external control for every internal control.</p></caption>
<graphic xlink:href="fbioe-04-00070-g002.tif"/>
</fig>
<p>The state-based microfluidic fluidic logic devices and circuits utilize static gain and normally closed valves (NCVs). NCVs are fabricated by filling specialized control channels with a flash curable prepolymer and curing while the valve is closed. The resulting closed valve exerts certain force against fluidic pressure to keep the valve closed. The valves are characterized by breakthrough pressure: the threshold pressure in the flow channel required to push open the valve and restore the continuity of the flow. Breakthrough pressure for an NCV can be tailored by controlling the pressure at which they are cured. Using these NCVs, we have developed static gain valves (SGV) that have the ability to control higher (or equal) fluidic pressure using a lower pressure. This type of valve is essential to create any logic/feedback structures (to account for signal strength losses), which can receive the output of the previous element/gate and use it as an input for decision making.</p>
<p>Utilizing the SGV, we next built an inverter (NOT gate), which was further used to build more complex circuits including bi-stable flip flops, clocked flip flops (latches), delay flip flops (D-FF, one bit of the SR), and complex microprocessors (SR). A SR that is capable of processing <italic>n</italic>&#x02009;&#x0002B;&#x02009;1 bits of data is formed by combining <italic>n</italic> D flip flops (bits of SR). The SR presented here uses air as the medium and receives three active high-pressure inputs: source, clock, and data (Figure <xref ref-type="fig" rid="F3">3</xref>). The pneumatic output of the SR cannot be used to control the flow of liquids in microchannels directly, due to risk of introducing bubbles. In order to address this issue, the signal medium is converted from air to liquid using an inverter.</p>
<fig id="F3" position="float">
<label>Figure 3</label>
<caption><p><bold>CAD drawing of microfluidic control components on Polaris mRNA-seq dosing IFC (A)</bold>. The shift register <bold>(B)</bold> enables active selection of single cells. The dilute single-cell suspension is loaded into a serpentine partition channel <bold>(D)</bold>. The cell suspension liquid flow is stopped, and the partition channel is imaged to identify single cells based on a particular set of fluorescence markers. The selected cells are then microfluidically moved downstream to a cell capture site <bold>(F)</bold> through a multiplexer <bold>(E)</bold>. The IFC is capable of accepting 20 reagents <bold>(G)</bold> as input. The shift register uses inverter <bold>(C)</bold> and a set of source, clock, and data <bold>(H)</bold>.</p></caption>
<graphic xlink:href="fbioe-04-00070-g003.tif"/>
</fig>
<p>The Polaris IFC microprocessor receives 28 external signals serially and processes them into 28 parallel independent controls capable of controlling individual valves or a set of valves. Five dedicated high-pressure external active signals are required for a SR. The CAD drawing of the various microfluidic components on a Polaris IFC is shown in Figure <xref ref-type="fig" rid="F3">3</xref>. The IFC can accept up to 20 independent reagents. The fluorescently labeled cells are loaded in a serpentine partition channel. Based on a desired combination of up to three fluorescent markers (refer to Section &#x0201C;<xref ref-type="sec" rid="S2-2">Polaris Instrument Design</xref>&#x0201D; for excitation and emission details), single cells are selected and sequentially isolated to the cell capture sites through a multiplexer. Up to 48 single cells can be sequestered on a single Polaris IFC. Subsequently, these 48 cells are processed through template-switching chemistry for full-length cDNA generation for mRNA-seq. In brief, the cells are lysed and reverse-transcribed, and full-length cDNA is preamplified by long and accurate PCR.</p>
</sec>
<sec id="S2-2">
<title>Polaris Instrument Design</title>
<p>The Fluidigm Polaris system (Figure <xref ref-type="fig" rid="F4">4</xref>A) consists of four major modules: (1) thermal control module; (2) imaging module; (3) pneumatic control module; and (4) environmental control (EC) module. The thermal module consists of a Peltier-based thermoelectric couple (TEC) device for heating/cooling. The TEC module can provide temperature in the range of 4&#x02013;99&#x000B0;C. Vacuum grooves on the thermal module are designed to enable tight contact with the glass-based integrated heat spreader (IHS) on the Polaris IFC. This ensures thermal uniformity across the fluidic circuit. The imaging module contains a five-color LED light engine for excitation (Ex wavelengths: 438, 475, 530, 575, and 632&#x02009;nm). The light source from the engine is collected and projected onto the fluidic circuit using fiber optics. The emitted signal from the fluidic circuit passes through an emission filter (five Em wavelengths: 488, 525, 570, 630, and 700&#x02009;nm) and is collected by CCD camera with 6-&#x003BC;m pixel resolution through a custom-designed collimator lens.</p>
<fig id="F4" position="float">
<label>Figure 4</label>
<caption><p><bold>(A)</bold> Components of Polaris instrument. The instrument consist of four major modules: (1) thermal module, enables preparative chemistry on sequestered single cells; (2) imaging module, consists of LED excitation and emission collection by a camera; (3) pneumatic module, controls the movement of reagents inside microfluidic channel by application of positive pressure on the IFC carrier; and (4) environmental control module, maintains the temperature, humidity, and blood gas flow rate for single-cell culture on-IFC. <bold>(B)</bold> Components of environmental control interface plate. The top of the interface plate contains glass coated with indium tin oxide. Internal heaters are used to maintain the temperature of the closed chamber between the interface plate and Polaris IFC. The temperature and relative humidity inside the closed chamber (after clamping with Polaris IFC) are measured by T/RH sensor. Blood or premixed gas required for cell culture is pumped through mixed gas inlet port on the interface plate. Polaris IFC is shown for reference.</p></caption>
<graphic xlink:href="fbioe-04-00070-g004.tif"/>
</fig>
<p>The pneumatic control module generates and stores air with volume up to 1&#x02009;L. The system can achieve a maximum pressure of 100&#x02009;psi. The pneumatic controller generates a vacuum on the thermal chuck, clamps the Polaris IFC against the EC interface plate (IP) to enable a closed environment around the IFC, and loads reagents from the inlets on the IFC carrier to the microchannels and reagent chambers of the fluidic circuit. The system contains segregated zones to regulate four different pressures simultaneously. The EC module provides an environment suitable for cell culture using user-desired gas composition. Environmental parameters such as temperature, relative humidity (RH), and mix gas flow rate across the fluidic circuit are monitored and controlled. The gas inside the closed chamber is heated by two heater coils. The gas inlet on the EC IP is used to regulate the flow of gas across the fluidic circuit. The EC IP (Figure <xref ref-type="fig" rid="F4">4</xref>B) contains an indium tin oxide (ITO) coated glass on the top to maintain thermal control in the EC while allowing imaging through the EC IP. During cell culture operation, the ITO glass is heated to prevent condensation. During cell culture, the environment around the fluidic circuit is maintained by blood gas (5% CO<sub>2</sub>, 5% oxygen, and 90% nitrogen) or premixed gas of choice (for example, 5% CO<sub>2</sub>, 20% oxygen, and 75% nitrogen). Before on-IFC cell culture, a rectangular sponge saturated with water is installed inside the closed chamber to provide the desired humidity through heating. The EC IP is equipped with two sensors (T/RH) to measure and maintain temperature at 37&#x000B0;C and RH at 90%.</p>
</sec>
<sec id="S2-3">
<title>K562 Cell Culture and CD59 Staining</title>
<p>K562 cells (ATCC<sup>&#x000AE;</sup> CCL-243) are cultured in T25 flasks in a volume ranging from 10 to 15&#x02009;mL in an incubator (37&#x000B0;C, 5% CO<sub>2</sub>). The culture medium contains IMDM&#x02009;&#x0002B;&#x02009;GlutaMAX&#x02122;-I&#x02009;&#x0002B;&#x02009;25&#x02009;mM HEPES&#x02009;&#x0002B;&#x02009;3.024&#x02009;g/L sodium bicarbonate (Gibco, 31980-030) and is supplemented with 10% FBS. The cells were fed every 2&#x02013;3&#x02009;days by dilution to 200,000 cells/mL. The K562 cells were stained with CellTracker&#x02122; Orange (CTO) CMRA Dye (Thermo Fisher Scientific, C34551) as universal marker and Alexa Fluor<sup>&#x000AE;</sup> 647 conjugated CD59 antibody. The recommended dyes and corresponding excitation and emission filters on the Polaris system are shown in Table <xref ref-type="table" rid="T1">1</xref>. Immediately before use, the cell staining solution was prepared by adding 0.6&#x02009;&#x003BC;L of 1&#x02009;mM CTO to 2&#x02009;mL of HBSS without calcium or magnesium (&#x02212;/&#x02212;) at a final concentration of 0.3&#x02009;&#x003BC;M. The cell staining solution was protected from light until use within 30&#x02009;min. A total of &#x0007E;1.5&#x02009;&#x000D7;&#x02009;10<sup>6</sup> cells was aliquoted in a 15&#x02009;mL non-pyrogenic conical tube. The cell suspension was centrifuged at 300&#x02009;&#x000D7;&#x02009;<italic>g</italic> for 3&#x02009;min. Following this, the medium was aspirated without disturbing the pellet, and 2&#x02009;mL of cell staining solution was added to the pellet and gently suspended by pipetting up and down three times. The cells were then incubated in the dark at 37&#x000B0;C for 20&#x02009;min with occasional inverting and flicking. Following this, the cells were washed by adding 12&#x02009;mL of HBSS to the cells in the 2&#x02009;mL of staining buffer and then centrifuged at 300&#x02009;&#x000D7;&#x02009;<italic>g</italic> for 5&#x02009;min. Supernatant was aspirated and discarded without disturbing the pellet. The pellet was then resuspended in 200&#x02009;&#x003BC;L of HBSS. The CTO-stained K562 cells were split into two tubes of 100&#x02009;&#x003BC;L each. One tube was used as negative surface-stained cell population, and the other tube was processed further to stain CD59 epitope. In order to stain the surface CD59 epitope, 10&#x02009;&#x003BC;L of CD59 biotinylated antibody (BD Biosciences, 555762, 100 tests, 2.0&#x02009;mL) was added to 100&#x02009;&#x003BC;L of CTO-stained cells. For negative surface-stained cell control, 10&#x02009;&#x003BC;L of HBSS was added. Both the tubes were incubated at room temperature for 20&#x02009;min with occasional inverting and flicking. Subsequently, 13&#x02009;mL of HBSS was added to each tube and centrifuged at 300&#x02009;&#x000D7;&#x02009;<italic>g</italic> for 5&#x02009;min. The pellet was resuspended in 100&#x02009;&#x003BC;L of HBSS. To this, 0.5&#x02009;&#x003BC;L of Streptavidin Alexa Fluor<sup>&#x000AE;</sup> 647 (Thermo Fisher Scientific, S32357, 2&#x02009;mg/mL stock) was added to positive-stain tube with CD59 biotinylated antibody in 100-&#x003BC;L cell suspension. This solution was mixed gently by pipetting up and down five times. Following this, the stain solution was incubated at room temperature for 15&#x02009;min with occasional flicking. Again, 13&#x02009;mL of HBSS was added to each tube, mixed by gently pipetting up and down, and centrifuged at 300&#x02009;&#x000D7;&#x02009;<italic>g</italic> for 5&#x02009;min. The supernatant was removed, and the pellet was resuspended in &#x0007E;100&#x02013;150&#x02009;&#x003BC;L culture medium with FBS, but without phenol red, to prevent high background fluorescence during cell selection on the Polaris system. The resuspension volume of culture medium accounts for cell losses during the staining procedure and was chosen to yield a cell concentration greater than the target concentration of 550 cells/&#x003BC;L. Typically, 10&#x02009;&#x003BC;L of cell mix is loaded into a C-Chip&#x02122; Disposable Hemocytometer (INCYTO, DHC-N01) and imaged on the Polaris system to estimate the staining intensity and purity. In order to achieve optimal buoyancy, cells in the range of 333&#x02013;550 cells/&#x003BC;L are mixed with suspension reagent (Fluidigm, 101-0434). Typically, the ratio of cells to cell suspension reagent is 3:2. However, this ratio might need optimization depending on the cell type.</p>
<table-wrap position="float" id="T1">
<label>Table 1</label>
<caption><p><bold>Recommended stains and corresponding excitation and emission filter on Polaris system</bold>.</p></caption>
<table frame="hsides" rules="groups">
<thead>
<tr>
<th valign="top" align="left">Channel name</th>
<th valign="top" align="left">Excitation<xref ref-type="table-fn" rid="tfn1"><sup>a</sup></xref></th>
<th valign="top" align="left">Emission</th>
<th valign="top" align="left">Recommended stains</th>
</tr>
</thead>
<tbody>
<tr>
<td align="left" valign="top">FAM&#x02122;</td>
<td align="left" valign="top">475/40</td>
<td align="left" valign="top">525/25</td>
<td align="left" valign="top">Alexa Fluor 488 (selection marker)</td>
</tr>
<tr>
<td align="left" valign="top">VIC<sup>&#x000AE;</sup></td>
<td align="left" valign="top">530/20</td>
<td align="left" valign="top">570/30</td>
<td align="left" valign="top">CellTracker Orange CMRA Dye (universal marker)</td>
</tr>
<tr>
<td align="left" valign="top">Cy5<sup>&#x000AE;</sup></td>
<td align="left" valign="top">632/28</td>
<td align="left" valign="top">700/30</td>
<td align="left" valign="top">Alexa Fluor 647 (selection marker)</td>
</tr>
</tbody>
</table>
<table-wrap-foot>
<fn id="tfn1"><p><italic><sup>a</sup>Excitation values are center wavelength/band pass (&#x02265;90%)</italic>.</p></fn>
</table-wrap-foot>
</table-wrap>
</sec>
<sec id="S2-4">
<title>IFC Operation</title>
<p>The Polaris IFC is first primed to fill the control lines on the fluidic circuit, load cell capture beads, and the inside of PDMS channels is blocked to prevent non-specific absorption/adsorption of proteins. In order to capture and maintain the single cells in the sites, the capture sites (48 sites) are preloaded with beads that are linked on-IFC to fabricate a tightly packed bead column during the IFC prime step. In the case of adherent cells, ECM is coated inside the cell capture chambers during prime step. After completion of the prime step, the cell mix (cells with suspension reagent) is loaded on the Polaris IFC and single CTO<sup>&#x0002B;</sup>/CD59<sup>&#x0002B;</sup> cells are selected to capture sites. We extensively tested the performance of the Polaris IFC and system at three different cell purities (3, 10, and 50%). The cell purity is defined as the ratio of CTO<sup>&#x0002B;</sup>/CD59<sup>&#x0002B;</sup> cells to CTO<sup>&#x0002B;</sup>/CD59<sup>&#x02212;</sup> cells. During the cell selection step, the suspended cells are loaded into the serpentine partition channel (Figure <xref ref-type="fig" rid="F3">3</xref>). Subsequently, the flow inside the partition channel is stopped (Figure <xref ref-type="fig" rid="F5">5</xref>A), and the cells are imaged in the partition channel for different fluorescent markers, as selected by the user. Based on automated image analyses by the system&#x02019;s software, only single cells with the desired combination of fluorescent markers are selected and isolated to the cell capture site. Any doublets or single cells with undesired fluorescent combinations are not selected by the software for further experimentation (Figure <xref ref-type="fig" rid="F5">5</xref>B). The selected single cells are moved to capture sites through a multiplexer (Figure <xref ref-type="fig" rid="F5">5</xref>C). The system takes images of the capture sites to confirm the arrival of single cells from a particular position in the partition to a particular capture site number. Figure <xref ref-type="fig" rid="F5">5</xref>C shows a typical image from the Polaris system showing K562 single cells captured inside sites packed with a column of beads. The system will then select and isolate all available single cells from a partition fill as per desired fluorescent marker combination. Once it completes selection of candidate cells, the system refills the serpentine partition to look for more candidate single cells. The system repeats this process to select and isolate single cells until it fills all 48 capture sites.</p>
<fig id="F5" position="float">
<label>Figure 5</label>
<caption><p><bold>(A)</bold> Cells in suspension are loaded into a serpentine partition channel (CAD design image). <bold>(B)</bold> Image analysis shows movement of fluorescent cells inside the serpentine channel during reagent flow. Once the reagent flow is stopped, single cells are separated from each other. The system software identifies single cells based on desired fluorescent markers. <bold>(C)</bold> After identification, single cells are isolated by moving them to the capture sites through a fluidic multiplexer. <bold>(D)</bold> Image of cultured BJ fibroblast cells. The Polaris IFC can be imaged on a microscope to obtain high-resolution micrographs.</p></caption>
<graphic xlink:href="fbioe-04-00070-g005.tif"/>
</fig>
<p>If desired, the single cells can then be cultured in the capture sites. It is possible to culture either suspension (e.g., K562) or adherent (e.g., BJ fibroblast) cells. For adherent cells, extracellular matrix can be coated inside the capture site during the IFC priming step. Figure <xref ref-type="fig" rid="F5">5</xref>D shows a Polaris image of a cultured BJ fibroblast (adhered). Based on the experimental design, it is possible to dose these single cells and on-IFC-cultured single cells with drugs or other cell stimuli. Finally, the single cells are processed through template-switching mRNA-seq chemistry for full-length cDNA generation and preamplification on-IFC.</p>
</sec>
<sec id="S2-5">
<title>Full-Length cDNA Generation</title>
<p>Preamplified full-length cDNA of selected single cells are generated on-IFC, and the amplicons are harvested through 48 different outlets. We used the SMARTer Ultra<sup>&#x000AE;</sup> Low RNA Kit for Illumina Sequencing (Clontech<sup>&#x000AE;</sup>, 634936) to generate preamplified cDNA. The selected and sequestered single cells were lysed using Polaris cell lysis mixture. The 28-&#x003BC;L cell lysis mix consists of 8.0&#x02009;&#x003BC;L of Polaris Lysis Reagent (Fluidigm, 101-1637), 9.6&#x02009;&#x003BC;L of Polaris Lysis Plus Reagent (Fluidigm, 101-1635), 9.0&#x02009;&#x003BC;L of 3&#x02032; SMART&#x02122; CDS Primer II A (12&#x02009;&#x003BC;M, Clontech, 634936), and 1.4&#x02009;&#x003BC;L of Loading Reagent (20X, Fluidigm, 101-1004). Synthetic RNA spikes can be optionally used with cell lysis mix. We typically use ArrayControl&#x02122; RNA spikes 1, 4, and 7 (Thermo Fisher Scientific, AM1780) to establish the functionality of RT and PCR on-IFC. We also use ERCC spikes at 1:50,000 dilution (final in lysis mix) for efficiency and quantification estimations. In order to implement synthetic RNA spikes, we thoroughly mix 96.5&#x02009;&#x003BC;L of loading reagent with 2.5&#x02009;&#x003BC;L of SMARTer Kit RNase Inhibitor (40&#x02009;U/&#x003BC;L; Clontech, 634936) and subsequently add 1&#x02009;&#x003BC;L of synthetic RNA spike to this spike mix. If RNA spike is used, then 1.4&#x02009;&#x003BC;L of the loading reagent is replaced with the spike mix. The thermal profile for single-cell lysis is 37&#x000B0;C for 5&#x02009;min, 72&#x000B0;C for 3&#x02009;min, 25&#x000B0;C for 1&#x02009;min, and hold at 4&#x000B0;C.</p>
<p>The 48-&#x003BC;L preparation volume for RT contains 1X SMARTer Kit 5X First-Strand Buffer (5X; Clontech, 634936), 2.5-mM SMARTer Kit Dithiothreitol (100&#x02009;mM; Clontech, 634936), 1-mM SMARTer Kit dNTP Mix (10&#x02009;mM each; Clontech, 634936), 1.2-&#x003BC;M SMARTer Kit SMARTer II A Oligonucleotide (12&#x02009;&#x003BC;M; Clontech, 634936), 1-U/&#x003BC;L SMARTer Kit RNase Inhibitor (40&#x02009;U/&#x003BC;L; Clontech, 634936), 10-U/&#x003BC;L SMARTScribe&#x02122; Reverse Transcriptase (100&#x02009;U/&#x003BC;L; Clontech, 634936), and 3.2&#x02009;&#x003BC;L of Polaris RT Plus Reagent (Fluidigm, 101-1366). All the concentrations correspond to those found in the RT chambers inside the Polaris IFC. The thermal protocol for RT is 42&#x000B0;C for 90&#x02009;min (RT), 70&#x000B0;C for 10&#x02009;min (enzyme inactivation), and a final hold at 4&#x000B0;C.</p>
<p>The 90-&#x003BC;L preparation volume for PCR contains 1X Advantage 2 PCR Buffer [not short amplicon (SA)] (10X, Clontech, 639206, Advantage<sup>&#x000AE;</sup> 2 PCR Kit), 0.4-mM dNTP Mix (50X/10&#x02009;mM, Clontech, 639206), 0.48-&#x003BC;M IS PCR Primer (12&#x02009;&#x003BC;M, Clontech, 639206), 2X Advantage 2 Polymerase Mix (50X, Clontech, 639206), and 1X Loading Reagent (20X, Fluidigm, 101-1004). All the concentrations correspond to those found in the PCR chambers inside the Polaris IFC. The thermal protocol for preamplification consists of 95&#x000B0;C for 1&#x02009;min (enzyme activation), five cycles (95&#x000B0;C for 20&#x02009;s, 58&#x000B0;C for 4&#x02009;min, and 68&#x000B0;C for 6&#x02009;min), nine cycles (95&#x000B0;C for 20&#x02009;s, 64&#x000B0;C for 30&#x02009;s, and 68&#x000B0;C for 6&#x02009;min), seven cycles (95&#x000B0;C for 30&#x02009;s, 64&#x000B0;C for 30&#x02009;s, and 68&#x000B0;C for 7&#x02009;min), and final extension at 72&#x000B0;C for 10&#x02009;min. The preamplified cDNAs are harvested into 48 separate outlets on the Polaris IFC carrier.</p>
</sec>
<sec id="S2-6">
<title>qPCR Analysis on Biomark&#x02122;</title>
<p>Harvested samples from Polaris IFCs were analyzed by qPCR using 96.96 Dynamic Array&#x02122; IFCs and the Biomark&#x02122; HD system from Fluidigm. Processing of the IFCs and operation of the instruments were performed according to the manufacturer&#x02019;s procedures. For detection using the RNA expression and splice variant assays, a Master Mix was prepared consisting of 360-&#x003BC;L SsoFast&#x02122; EvaGreen<sup>&#x000AE;</sup> Supermix with Low ROX (BioRad 172-5211) and 36-&#x003BC;L 20&#x02009;&#x000D7;&#x02009;DNA Binding Dye Sample Loading Reagent (Fluidigm 100-5360), and 3.3&#x02009;&#x003BC;L of this mix was dispensed to each well of a 96-well assay plate. Diluted harvest product (2.7&#x02009;&#x003BC;L) was added to each well, and the plate was briefly vortexed and centrifuged. Following priming of the IFC in the IFC Controller HX, 5&#x02009;&#x003BC;L of the sample&#x02009;&#x0002B;&#x02009;Master Mix were dispensed to each sample inlet of the 96.96 IFC. Five microliters of each 10&#x02009;&#x000D7;&#x02009;assay [5&#x02009;&#x003BC;M each primer, 1&#x02009;&#x000D7;&#x02009;assay Loading Reagent (Fluidigm 100-5359)] were dispensed to each Detector Inlet of the 96.96 IFC. After loading the assays and samples into the IFC in the IFC Controller HX, the IFC was transferred to the Biomark HD, and PCR was performed using the thermal protocol GE Fast 96&#x02009;&#x000D7;&#x02009;96 PCR&#x02009;&#x0002B;&#x02009;Melt v2.pcl. This protocol consists of a thermal mix of 70&#x000B0;C, 40&#x02009;min; 60&#x000B0;C, 30&#x02009;s, hot start at 95&#x000B0;C, 1&#x02009;min, PCR cycle of 30 cycles of 96&#x000B0;C, 5&#x02009;s; 60&#x000B0;C, 20&#x02009;s, and melting using a ramp from 60 to 95&#x000B0;C at 1&#x000B0;C/3&#x02009;s. Data were analyzed using Fluidigm Real-Time PCR Analysis software using the Linear (Derivative) Baseline Correction Method and the Auto (Global) Ct Threshold Method. The data are exported as a.csv file into an Excel<sup>&#x000AE;</sup> macro to compile and compare the data against in-house specifications.</p>
</sec>
</sec>
<sec id="S3">
<title>Results</title>
<sec id="S3-1">
<title>Performance Evaluation of Polaris IFC</title>
<p>In order to statistically evaluate the performance of the Polaris IFC, we designed and developed two performance tests: (1) total-RNA-based performance test (RNA PT) and (2) single-cell-based key performance test (KPT). Since single cells are heterogeneous, it would be difficult to evaluate the performance uniformity across 48 capture sites using a cell-based test method. Hence, we developed a 20-cell-equivalent total-RNA PT to evaluate and improve the performance of the Polaris IFC during the initial phase of the IFC development process.</p>
<sec id="S3-1-1">
<title>Total RNA-Based Performance Test</title>
<p>The primary objective of this test is to statistically validate a workflow that is very close to the cell-based experiments on the Polaris system and yet collects critical information about uniformity of cDNA synthesis across IFC, reaction line cross-talk (on-IFC), and IFC-to-IFC correlation. To achieve this objective, we simulated steps such as loading of cell capture beads and the thermal step for cell lysis in the total-RNA PT. The workflow of the total-RNA PT is shown in Figure <xref ref-type="fig" rid="F6">6</xref>A. Briefly, the RNA-PT is a two-step procedure. In the first step, the control lines on the Polaris IFC are primed, channels are blocked, and cell capture beads are back-loaded with ArrayControl RNA SPIKES (1, 4, and 7 only, Thermo Fisher Scientific, AM1780; henceforth referred to as RNAspikes 147) in eight specific capture sites (Figure <xref ref-type="fig" rid="F6">6</xref>B).</p>
<fig id="F6" position="float">
<label>Figure 6</label>
<caption><p><bold>Total-RNA-based performance test (two-step) workflow and pipetting map</bold>. <bold>(A)</bold> Workflow to prime Polaris IFC with beads, back-load them with RNAspikes 147, simulate front dosing with diluted ERCC spikes, and generate cDNA using mRNA-seq chemistry. The cDNA amplicons from Polaris IFC were analyzed on the Fluidigm M96.96 IFC using 85 qPCR assays for specific genes, 8 assays for ERCC RNA spikes, and 3 assays for RNAspikes 147. <bold>(B)</bold> Pipetting map for IFC prime step. During prime step, RNAspikes 147 are back-loaded to 8 specific inlets with cell capture beads. <bold>(C)</bold> Pipetting map for front-dosing simulation with ERCC RNA spikes, followed by mRNA-seq chemistry.</p></caption>
<graphic xlink:href="fbioe-04-00070-g006.tif"/>
</fig>
<p>The ArrayControl RNA Spikes are used to evaluate the back-dosing cross-talk using highly sensitive qPCR assay designed to detect RNA spikes 1, 4, and 7 (three total ArrayControl RNA Spikes). After the priming step, six specific capture sites are loaded with ERCC RNA Spike-In Mix (Thermo Fisher Scientific, 4456740) to estimate the cross-talk for front-loaded reagents and dosing agents. Although the ERCC RNA mix contains 92 spike-ins, only 8 ERCC spike-ins were probed using qPCR assays. The front-dosing strategy is illustrated in Figure <xref ref-type="fig" rid="F7">7</xref>A and the pipetting map is shown in Figure <xref ref-type="fig" rid="F6">6</xref>C.</p>
<fig id="F7" position="float">
<label>Figure 7</label>
<caption><p><bold>(A)</bold> Front-dosing strategy. The capture sites on Polaris IFC are segregated as six groups of eight sites. ERCC spike-in mix is loaded into the second capture site of every group (Group alt 2). For negative control, 1X Preloading Reagent was loaded into Group alt 4, 6, and 8. <bold>(B)</bold> Sequential chambers for multistep chemistry. The capture site with beads holds 2&#x02009;nL of liquid. The lysis mix is loaded into chamber 1 (9&#x02009;nL), then reverse transcription mix is loaded into chambers 2 and 3 (9&#x02009;nL each), and PCR mix is loaded into chambers 4 and 5 (135&#x02009;nL each). <bold>(C)</bold> Chemistry loading sequence and its function.</p></caption>
<graphic xlink:href="fbioe-04-00070-g007.tif"/>
</fig>
<p>For negative control, 1X Preloading Reagent (Fluidigm, 100-9942) was loaded into specific inlets and capture sites (Figure <xref ref-type="fig" rid="F7">7</xref>C). After completion of front dosing, the mRNA-seq chemistry prep is integrated with the dosing step. The lysis mixture for the RNA-PT contains Leukemia (K562) Total RNA (Thermo Fisher Scientific, AM7832) at a concentration equivalent to 20 cells of total RNA in every cell capture site (48 sites). The cell capture site is serially connected to five chambers to enable multistep reaction chemistry (Figure <xref ref-type="fig" rid="F7">7</xref>B). Cell lysis mixture is loaded into the first 9-nL chamber. Then, RT mixture is loaded in 18-nL volume (two 9-nL chambers). Finally, PCR mixture for preamplification of full-length cDNA is loaded in 270-nL volume (two 135-nL chambers).</p>
<p>The preamplified cDNA is harvested in &#x0007E;7&#x02009;&#x003BC;L volume. The harvest is further diluted by addition of 10&#x02009;&#x003BC;L of DNA Dilution RGT (Fluidigm, 100-9167). In order to evaluate IFC uniformity and other performance metrics, we designed 88 Delta Gene&#x02122; assays (Fluidigm) for K562 (85 genes covering high and low expressors) and RNAspikes 147. In addition to these 88 assays, we used 8 ERCC qPCR assays from a published work (Devonshire et al., <xref ref-type="bibr" rid="B11">2011</xref>). In total, we used 96 intercalating dye-based qPCR assays for read-out of RNA PT on an M96.96 Dynamic Array&#x02122; IFC (Fluidigm). We routinely test positive and negative tube controls for every chemistry preparation by qPCR assays on the M96.96. In order to do this, 2 out of the 48 samples from a Polaris IFC are replaced by positive and negative tube controls on M96.96. The positive tube control contains total RNA from K562, RNAspikes 147, and ERCC. The tube controls are used to validate the functionality of chemistry preparation on a particular day. Harvest products from two Polaris IFC are tested on a single M96.96 Dynamic Array IFC run. A typical qPCR Ct heat map and associated Excel macro for two Polaris RNA PTs are shown in Figures <xref ref-type="fig" rid="F8">8</xref>A,B. In order to statistically validate the performance, we tested 44 Polaris IFCs with RNA PTs. IFC and reagents from minimum of three manufacturing lots were used. Tolerance limit or interval analyses were performed on more than 40 Polaris IFC runs. The distribution of data and tolerance limit analyses for different metrics for the RNA PTs are shown in Figure <xref ref-type="fig" rid="F9">9</xref>.</p>
<fig id="F8" position="float">
<label>Figure 8</label>
<caption><p><bold>(A)</bold> Typical heat map of high-throughput qPCR assay for RNA-based performance test. The M96.96 IFC (96 samples) can accept amplicons from two Polaris IFCs (48 samples each). For every Polaris IFC, we replace two samples with positive and negative control samples. The columns are assays (85 high- and low-expressing assays; 8 ERCC spike assays; and 3 RNAspike 147 assays). The rows are diluted amplicons from Polaris IFCs. <bold>(B)</bold> Excel macro for the RNA-based performance test.</p></caption>
<graphic xlink:href="fbioe-04-00070-g008.tif"/>
</fig>
<fig id="F9" position="float">
<label>Figure 9</label>
<caption><p><bold>Tolerance limit analyses from 44 Polaris IFC runs using RNA-based performance test</bold>. <italic>n</italic> is the sample number after removal of outliers. Outliers are estimated by Tukey&#x02019;s box plot. The percentage of outliers removed ranged from 0 to 5%. <bold>(A)</bold> Average assay SD. The data fit a normal distribution with goodness of fit <italic>P</italic>-value of 0.2813. The upper tolerance interval for this metric is 0.395. <bold>(B)</bold> Distribution of drop outs for high-expressing assays only. <bold>(C)</bold> Distribution of drop outs for both high and low expressors. <bold>(D)</bold> Back-loaded dosing cross-talk. <bold>(E)</bold> Distribution of slope for IFC-to-IFC correlation. The data fit a normal distribution with tolerance interval of 1&#x02013;1.08. <bold>(F)</bold> Distribution of <italic>R</italic><sup>2</sup> for IFC-to-IFC correlation. The data fit a normal distribution with lower tolerance limit of 0.959.</p></caption>
<graphic xlink:href="fbioe-04-00070-g009.tif"/>
</fig>
</sec>
<sec id="S3-1-2">
<title>Single-Cell-Based Key Performance Test</title>
<p>Key performance test was developed and validated using one cell type each for suspension (K562) and adherent (BJ fibroblast) cells. As described in Section &#x0201C;<xref ref-type="sec" rid="S2">Materials and Methods</xref>,&#x0201D; cells are stained with the universal fluorescent marker, CTO. A subset of these cells were stained for surface marker using antibody conjugated with Alexa 647. In the case of K562, we used Anti-Human CD59-Biotin (BD Biosciences, 555762) with Streptavidin Alexa 647, and for BJ fibroblast, we used Anti-Mouse/Human CD44-Alexa 647 (BioLegend 103018). The double-stained cells (universal CTO and surface marker Alexa 647) were mixed with cells stained with CTO only to achieve three different purity percentages (3, 10, and 50%). The cells were selected for universal CTO and surface marker. For BJ fibroblasts, we tested two different workflows, one with cell selection followed by chemistry (immediate) and another to dose the BJ fibroblasts with medium every 4&#x02009;h for 24-h adherent culture, followed by chemistry (BJ dosing). In order to evaluate the cell viability prior to the cell lysis step, we used Zombie Yellow&#x02122; cell viability stain (BioLegend, 423103; &#x003BB;<sub>ex</sub>&#x02009;&#x0003D;&#x02009;396&#x02009;nm and &#x003BB;<sub>em</sub>&#x02009;&#x0003D;&#x02009;572&#x02009;nm), which stains dead cells. Performance metrics, such as number of sites occupied with single cells out of the total 48 sites (cell selection), number of cells retained after dosing and prior to cell lysis (cell retention), and number of viable cells prior to lysis, were evaluated. On average from 20 Polaris IFC runs, our cell selection was &#x0007E;95% for K562 and BJ fibroblast with different purity percentages (Figure <xref ref-type="fig" rid="F10">10</xref>A). For cell retention, &#x0003E;47/48 sites showed presence of single cells as enumerated after the cell selection step and prior to cell lysis step (Figure <xref ref-type="fig" rid="F10">10</xref>B). The average cell viability was &#x0007E;90% as estimated from 20 Polaris IFC runs (Figure <xref ref-type="fig" rid="F10">10</xref>C).</p>
<fig id="F10" position="float">
<label>Figure 10</label>
<caption><p><bold>Single-cell selection, retention, and viability performance using single-cell key performance test</bold>. <bold>(A)</bold> Percentage of single cells selected (out of 48 total capture sites) for K562 (suspension) and BJ fibroblast (adherent) at different cell purity concentrations (3, 10, and 50%). <italic>n</italic> is number of Polaris IFCs tested. <bold>(B)</bold> Single cells retained after dosing and culture. Cells are imaged and counted before the cell lysis step. K562 and BJ fibroblasts were selected and immediately processed for mRNA-seq chemistry preparation. BJ fibroblasts were also selected, cultured, and dosed with cell culture medium for 24&#x02009;h (BJ dosing). Data from three IFCs are presented here. <bold>(C)</bold> Single-cell viability as assessed by Zombie stain on-IFC.</p></caption>
<graphic xlink:href="fbioe-04-00070-g010.tif"/>
</fig>
<p>The Polaris system generates very high quality (size distribution) and quantity (yield) of preamplified cDNA from single cells. The size distribution of preamplified cDNA from single cells, as evaluated using Bioanalyzer 2100 and the DNA high-sensitivity chip (Agilent), is typically in the size range of 0.3&#x02013;7&#x02009;kb (Figure <xref ref-type="fig" rid="F11">11</xref>A). For yield, preamplified cDNA from single cells was quantified using PicoGreen-based dsDNA quantification assay (Quant-iT&#x02122; PicoGreen<sup>&#x000AE;</sup> dsDNA Assay Kit, Thermo Fisher Scientific, P7589). The average total cDNA yield per single K562 cell is 38.42&#x02009;&#x000B1;&#x02009;8.08&#x02009;ng (Figure <xref ref-type="fig" rid="F11">11</xref>B). We randomly selected &#x0007E;14 single cells from three Polaris IFC runs, barcoded them using modified Nextera library prep, and pooled them to generate a single sequencing library. A representative library profile from 42 single cells is shown in Figure <xref ref-type="fig" rid="F11">11</xref>C. The majority of the single-cell library falls in the range of 200&#x02013;2,000&#x02009;bp. For three sequencing libraries from nine Polaris IFCs tested with K562 immediate chemistry, sequencing data from three MiSeq&#x02122; runs using v2 150&#x02009;bp PE kit were compiled, and tolerance limits (90% confidence with 95% population coverage) were estimated for two key sequencing metrics (Figure <xref ref-type="fig" rid="F12">12</xref>). The average percentage of reads mapping to rRNA/total reads is 0.122%. The Box&#x02013;Cox transformed data fit a normal distribution with a Shapiro&#x02013;Wilk <italic>P</italic>-value of 0.0983. Based on the normal distribution, the upper tolerance limit for percentage of reads mapping to rRNA is 0.3% (Figure <xref ref-type="fig" rid="F12">12</xref>A). The mean number of genes detected is 6,967&#x02009;&#x000B1;&#x02009;115. The data fit a normal distribution with a lower tolerance limit of 5,919 genes as estimated from 115 single-cell datapoints (Figure <xref ref-type="fig" rid="F12">12</xref>B).</p>
<fig id="F11" position="float">
<label>Figure 11</label>
<caption><p><bold>(A)</bold> Typical cDNA distribution from single cells from three different Polaris IFC runs (A1&#x02013;A3). The peak for RNAspikes 147 is around 1,000&#x02009;bp. <bold>(B)</bold>&#x02009;Distribution of cDNA yield from 358 single K562 cells from multiple Polaris IFC runs. The average cDNA yield per cell is 38.42&#x02009;ng with a SD of 8.08&#x02009;ng. <bold>(C)</bold>&#x02009;Library profile of 42 single cells.</p></caption>
<graphic xlink:href="fbioe-04-00070-g011.tif"/>
</fig>
<fig id="F12" position="float">
<label>Figure 12</label>
<caption><p><bold>Distribution of sequencing metrics for single cells</bold>. <bold>(A)</bold> Percentage of reads mapped to rRNA/total reads. The average mapping to rRNA is 0.122% from 116 single-cell sequence datapoints. The upper tolerance limit for percentage rRNA mapping is 0.3%. <bold>(B)</bold> Number of genes detected (TPM&#x02009;&#x0003E;&#x02009;1). Distribution of genes detected from 115 single-cell sequence data. The average genes detected is 6,967, with a SD of 115. The data fit a normal distribution, with a lower tolerance limit of 5,919.</p></caption>
<graphic xlink:href="fbioe-04-00070-g012.tif"/>
</fig>
<p>We extensively analyzed our single-cell sequencing data for transcript coverage bias and possible positional bias of single cells selected across 48 capture sites on the Polaris IFC (Figure <xref ref-type="fig" rid="F13">13</xref>). We noted uniform coverage along the transcript length (Figure <xref ref-type="fig" rid="F13">13</xref>B) without any positional bias on the Polaris IFC. The plot of normalized coverage vs. normalized distance along the transcript with respect to capture sites (2, 3, 4 and 40, 41, 42) from a Polaris IFC is shown in Figure <xref ref-type="fig" rid="F13">13</xref>B. A plot of median 3&#x02032; end bias of transcript coverage with respect to capture site number indicates no positional bias across three Polaris IFC runs (Figure <xref ref-type="fig" rid="F13">13</xref>A). In order to understand if there is any possible effect of hypoxia on single cells due to spatial location of capture sites on the Polaris IFC, we analyzed the expression value of <italic>HIF1A</italic> gene across different capture sites. Up-regulation of hypoxia-induced factor 1 gene (<italic>HIF1A</italic>) is a known consequence of hypoxia (Choudhry and Mole, <xref ref-type="bibr" rid="B8">2015</xref>). Expression analyses of <italic>HIF1A</italic> did not show any positional bias with respect to the capture sites (Figure <xref ref-type="fig" rid="F13">13</xref>C). It should be noted that we recommend strictly following the Polaris workflow as described in the Polaris protocol document (Fluidigm, 101-0082). Any deviation from the validated workflow might lead to introduction of possible bias at multiple levels.</p>
<fig id="F13" position="float">
<label>Figure 13</label>
<caption><p><bold>Bias vs. capture site number</bold>. <bold>(A)</bold> Median 3&#x02032; end bias was analyzed for random single cells picked for sequencing. No 3&#x02032; bias was noted for cells from three Polaris IFC runs <bold>(B)</bold> <italic>HIF1A</italic> gene expression of cells across different capture sites. <bold>(C)</bold> Sequencing coverage as estimated by Picard for single cells. Uniform coverage noted across different capture sites from a Polaris IFC.</p></caption>
<graphic xlink:href="fbioe-04-00070-g013.tif"/>
</fig>
</sec>
</sec>
<sec id="S3-2">
<title>Sensitivity Studies Using ERCC Spike-Ins</title>
<p>An alternative way to evaluate performance of single-cell mRNA-seq on the Polaris system is to implement use of the ERCC RNA Spike-In Mix 1 in the lysis mix. The ERCC control mix consists of 92 polyadenylated transcripts with a size range of 273&#x02013;2,022 bases and six orders of magnitude range in concentration. We tested both qPCR- and sequencing-based methods for detection of ERCC spikes. Ninety-two primer pairs were designed to target the corresponding transcripts for qPCR testing. qPCR was performed on a Fluidigm 96.96 Dynamic Array IFC. Stochastic distribution of transcripts was observed when the input concentration was 25 copies per reaction or less on the Polaris IFC followed by qPCR detection on the 96.96 Dynamic Array (Figure <xref ref-type="fig" rid="F14">14</xref>A). Single-copy RNA detection is demonstrated, although intermittently, likely due to sampling at the reaction site. Transcripts at 1.6 copies per reaction were intermittently detected by qPCR on the 96.96 Dynamic Array IFC. We also evaluated the detection rate of ERCC spikes (&#x0003E;1.6 copies/reaction) using an approach based on massive parallel sequencing. There were 7 ERCC spikes (ERCC-00170; ERCC-00148; ERCC-00126; ERCC-00099; ERCC-00054; ERCC-00163; ERCC-00059), which were at a concentration of 1.6 copies per Polaris reaction chamber. One of the 7 ERCC spikes (ERCC-00054) was not detected in any of the 19 single-cell samples. If we remove this datapoint as an outlier, the average detection rate of &#x0007E;1.6 copies is 28%. Based on Poisson estimates, single-copy detection rate should be &#x0007E;33% (67% should be a failure event). The single-copy detection rate (28%) from the Polaris system is very close to expected theoretical estimates based on Poisson statistics (Figure <xref ref-type="fig" rid="F14">14</xref>B).</p>
<fig id="F14" position="float">
<label>Figure 14</label>
<caption><p><bold>(A)</bold> Heat map (qPCR assay) for ERCC spike-in with single SKBR3 cells. The columns are 92 qPCR assays designed and developed for ERCC spike-in. The rows are amplicons from single SKBR3 cells with ERCC spike-in. The ERCC assays on the columns are sorted by concentrations. <bold>(B)</bold> Detection rate of ERCC spikes with &#x0003E;1.58 copies per Polaris IFC chambers.</p></caption>
<graphic xlink:href="fbioe-04-00070-g014.tif"/>
</fig>
</sec>
<sec id="S3-3">
<title>Single-Cell Transfection of nGFP mRNA and GFP Expression Analyses</title>
<p>In order to demonstrate dosing and functional response analyses, we transfected single K562 cells with nuclear green fluorescent protein (nGFP) mRNA and cultured the transfected single K562 for 16&#x02009;h. During this culture duration, the cells translated the nGFP mRNA and expressed the GFP inside the cell. The imaging capability of Polaris enables monitoring of GFP expression. Subsequently, the cells were processed for mRNA-seq chemistry on-IFC and sequenced on MiSeq to quantify the reads mapped to GFP. K562 cells stained with CTO were selected on Polaris IFC. To carry out the single-cell transfection, 10&#x02009;&#x003BC;L of Stemfect RNA transfection reagent was mixed with 240-&#x003BC;L Stemfect transfection buffer (Stemgent<sup>&#x000AE;</sup> Stemfect&#x02122; RNA Transfection Kit, 00-0069) (Mix A). The stock nGFP mRNA (Stemgent, 05-0019) at 100&#x02009;ng/&#x003BC;L was diluted with Stemfect transfection buffer first and then further diluted with Mix A to make mRNA transfection complex. This complex was incubated at room temperature for 15&#x02009;min and further diluted with K562 cell culture medium (refer to Section &#x0201C;<xref ref-type="sec" rid="S2-3">K562 Cell Culture and CD59 Staining</xref>&#x0201D;) to achieve different final concentrations (0.5 and 1&#x02009;ng/&#x003BC;L) of nGFP mRNA. Selected single K562 cells were cultured with the nGFP mRNA transfection complex with culture medium at 37&#x000B0;C with 5% CO<sub>2</sub> on the Polaris IFC. During this cell culture incubation time, images were taken every hour to monitor the onset of GFP expression. Image analyses (Figure <xref ref-type="fig" rid="F15">15</xref>) showed that single cells picked up nGFP mRNA at 0.5 and 1&#x02009;ng/&#x003BC;L concentrations and expressed the green fluorescent proteins, thereby reinforcing the fact that single cells on Polaris IFC are healthy and are able to uptake naked mRNA and translate it to protein capable of being transfected. Figure <xref ref-type="fig" rid="F15">15</xref>A shows typical time-series images, which can be obtained from the Polaris system. It should be noted that for this particular experiment, the imaging interval was set to 1&#x02009;h, but the Polaris system is capable of taking successive images in a rapid mode. We noted onset of GFP gene expression around the 3-h time frame at single-cell resolution. The cDNA pool showed a length range from 0.3 to 9.2&#x02009;kb, with an average length &#x0007E;2&#x02009;kb. It is also noted that &#x0003E;85% of the total cDNA pool lies between 0.5 and 9.2&#x02009;kb (Figure <xref ref-type="fig" rid="F15">15</xref>B). Sequencing data show that the cells transfected with nGFP mRNA harbored the extracellular mRNA even after 16&#x02009;h of culture. As expected, the control cells without nGFP transfection did not show any mapping to GFP sequence. The transfection of nGFP did not alter the mapping to genome and transcriptome when compared to the control cells (Figure <xref ref-type="fig" rid="F15">15</xref>C). The nGFP-transfected cells showed percentage average reads mapping of 0.57, 87.99, and 47.82% to GFP, genome, and transcriptome, respectively (<italic>n</italic>&#x02009;&#x0003D;&#x02009;7), while the control K562 showed percentage average mapping of 0, 88.03, and 49.22 (<italic>n</italic>&#x02009;&#x0003D;&#x02009;7).</p>
<fig id="F15" position="float">
<label>Figure 15</label>
<caption><p><bold>(A)</bold> Time-series images from the Polaris system taken every hour to monitor onset of GFP expression. For visual clarity, only two cell capture sites are shown. The onset GFP expression is noted around 3&#x02009;h from the time of transfection. <bold>(B)</bold> cDNA yield and distribution from single K562 processed through mRNA-seq chemistry. <bold>(C)</bold> Mapping metrics of GFP-transfected single K562 cells to genome, transcriptome, and GFP sequence.</p></caption>
<graphic xlink:href="fbioe-04-00070-g015.tif"/>
</fig>
</sec>
</sec>
<sec id="S4" sec-type="discussion">
<title>Discussion</title>
<p>In this work, we report design and development of an integrated system to perform functional studies on single cells. We developed a nanoscale IFC, which employs fluidic logic to actively select single cells, and a system capable of performing multiple functionalities. The performance of the developed IFC and system was extensively tested using RNA-based and single-cell-based performance tests. These tests were specifically designed to evaluate different functionalities of the IFC and system. The functional capability of the Polaris IFC and system has been successfully demonstrated using transfection of naked nGFP mRNA, followed by monitoring of nGFP expression and finally analysis of the whole set of mRNA transcripts by massive parallel sequencing. It is noted that it is not currently possible to perform studies reported in this work on any other single-cell platforms. The limitation of the current system includes limited number of cells for functional studies (up to 48 cells). However, the requirement on the number of cells depends on the biological question, and it is possible to expand the capability of the IFC consumable to process more cells in the future.</p>
</sec>
<sec id="S5">
<title>Author Contributions</title>
<p>NR, BF, JS, and JAAW conceived and designed the RNA-based performance test; BF, LS, AAL, JS, and JAAW conceived and designed the single-cell-based performance test; NR, LS, AAL, JS, MLG, CDS, NDA, CG, CTL, IH, AO, CS, and JAAW performed experiments; BF, NSGKD, MZ, EO, and CP were involved in Polaris IFC development; KH, MTM, WY, MN, CC, ML, HC, ZH, LL, CC, and ZS were involved in Polaris system development; RY, WH, JA, and ZH were involved in Polaris software development; NR, LS, JS, CDS, XW, and JAAW analyzed the data; TS edited the manuscript; CL drafted the Polaris user guide; MU and JAAW supervised the project, helped with design and interpretation, and provided laboratory space and financial support; and NR, LS, KJL, and JAAW wrote the manuscript with input from all authors. All authors listed, have made substantial, direct and intellectual contribution to the work, and approved it for publication.</p>
</sec>
<sec id="S6">
<title>Conflict of Interest Statement</title>
<p>All authors are employees of Fluidigm Corporation.</p>
</sec>
</body>
<back>
<ref-list>
<title>References</title>
<ref id="B1"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Alix-Panabi&#x000E8;res</surname> <given-names>C.</given-names></name> <name><surname>Bartkowiak</surname> <given-names>K.</given-names></name> <name><surname>Pantel</surname> <given-names>K.</given-names></name></person-group> (<year>2016</year>). <article-title>Functional studies on circulating and disseminated tumor cells in carcinoma patients</article-title>. <source>Mol. Oncol.</source> <volume>10</volume>, <fpage>443</fpage>&#x02013;<lpage>449</lpage>.<pub-id pub-id-type="doi">10.1016/j.molonc.2016.01.004</pub-id><pub-id pub-id-type="pmid">26847851</pub-id></citation></ref>
<ref id="B2"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Angermueller</surname> <given-names>C.</given-names></name> <name><surname>Clark</surname> <given-names>S. J.</given-names></name> <name><surname>Lee</surname> <given-names>H. J.</given-names></name> <name><surname>Macaulay</surname> <given-names>I. C.</given-names></name> <name><surname>Teng</surname> <given-names>M. J.</given-names></name> <name><surname>Hu</surname> <given-names>T. X.</given-names></name> <etal/></person-group> (<year>2016</year>). <article-title>Parallel single-cell sequencing links transcriptional and epigenetic heterogeneity</article-title>. <source>Nat. Methods</source> <volume>13</volume>, <fpage>229</fpage>&#x02013;<lpage>232</lpage>.<pub-id pub-id-type="doi">10.1038/nmeth.3728</pub-id><pub-id pub-id-type="pmid">26752769</pub-id></citation></ref>
<ref id="B3"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Avraham</surname> <given-names>R.</given-names></name> <name><surname>Haseley</surname> <given-names>N.</given-names></name> <name><surname>Brown</surname> <given-names>D.</given-names></name> <name><surname>Penaranda</surname> <given-names>C.</given-names></name> <name><surname>Jijon</surname> <given-names>H. B.</given-names></name> <name><surname>Trombetta</surname> <given-names>J. J.</given-names></name> <etal/></person-group> (<year>2015</year>). <article-title>Pathogen cell-to-cell variability drives heterogeneity in host immune responses</article-title>. <source>Cell</source> <volume>162</volume>, <fpage>1309</fpage>&#x02013;<lpage>1321</lpage>.<pub-id pub-id-type="doi">10.1016/j.cell.2015.08.027</pub-id><pub-id pub-id-type="pmid">26343579</pub-id></citation></ref>
<ref id="B4"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bendall</surname> <given-names>S. C.</given-names></name> <name><surname>Simonds</surname> <given-names>E. F.</given-names></name> <name><surname>Qiu</surname> <given-names>P.</given-names></name> <name><surname>Ad</surname> <given-names>A.</given-names></name> <name><surname>Krutzik</surname> <given-names>P. O.</given-names></name> <name><surname>Finck</surname> <given-names>R.</given-names></name> <etal/></person-group> (<year>2011</year>). <article-title>Single-cell mass cytometry of differential immune and drug responses across a human hematopoietic continuum</article-title>. <source>Science</source> <volume>332</volume>, <fpage>687</fpage>&#x02013;<lpage>696</lpage>.<pub-id pub-id-type="doi">10.1126/science.1198704</pub-id><pub-id pub-id-type="pmid">21551058</pub-id></citation></ref>
<ref id="B5"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Bianconi</surname> <given-names>E.</given-names></name> <name><surname>Piovesan</surname> <given-names>A.</given-names></name> <name><surname>Facchin</surname> <given-names>F.</given-names></name> <name><surname>Beraudi</surname> <given-names>A.</given-names></name> <name><surname>Casadei</surname> <given-names>R.</given-names></name> <name><surname>Frabetti</surname> <given-names>F.</given-names></name> <etal/></person-group> (<year>2013</year>). <article-title>An estimation of the number of cells in the human body</article-title>. <source>Ann. Hum. Biol.</source> <volume>40</volume>, <fpage>463</fpage>&#x02013;<lpage>471</lpage>.<pub-id pub-id-type="doi">10.3109/03014460.2013.807878</pub-id><pub-id pub-id-type="pmid">23829164</pub-id></citation></ref>
<ref id="B6"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Briggs</surname> <given-names>S. F.</given-names></name> <name><surname>Dominguez</surname> <given-names>A. A.</given-names></name> <name><surname>Chavez</surname> <given-names>S. L.</given-names></name> <name><surname>Reijo Pera</surname> <given-names>R. A.</given-names></name></person-group> (<year>2015</year>). <article-title>Single-cell XIST expression in human preimplantation embryos and newly reprogrammed female induced pluripotent stem cells</article-title>. <source>Stem Cells</source> <volume>33</volume>, <fpage>1771</fpage>&#x02013;<lpage>1781</lpage>.<pub-id pub-id-type="doi">10.1002/stem.1992</pub-id><pub-id pub-id-type="pmid">25753947</pub-id></citation></ref>
<ref id="B7"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Cayrefourcq</surname> <given-names>L.</given-names></name> <name><surname>Mazard</surname> <given-names>T.</given-names></name> <name><surname>Joosse</surname> <given-names>S.</given-names></name> <name><surname>Solassol</surname> <given-names>J.</given-names></name> <name><surname>Ramos</surname> <given-names>J.</given-names></name> <name><surname>Assenat</surname> <given-names>E.</given-names></name> <etal/></person-group> (<year>2015</year>). <article-title>Establishment and characterization of a cell line from human circulating colon cancer cells</article-title>. <source>Cancer Res.</source> <volume>75</volume>, <fpage>892</fpage>&#x02013;<lpage>901</lpage>.<pub-id pub-id-type="doi">10.1158/0008-5472.CAN-14-2613</pub-id><pub-id pub-id-type="pmid">25592149</pub-id></citation></ref>
<ref id="B8"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Choudhry</surname> <given-names>H.</given-names></name> <name><surname>Mole</surname> <given-names>D. R.</given-names></name></person-group> (<year>2015</year>). <article-title>Hypoxic regulation of the noncoding genome and NEAT1</article-title>. <source>Brief. Funct. Genomics</source>. <volume>15</volume>, <fpage>174</fpage>&#x02013;<lpage>185</lpage>.<pub-id pub-id-type="doi">10.1093/bfgp/elv050</pub-id></citation></ref>
<ref id="B9"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Darmanis</surname> <given-names>S.</given-names></name> <name><surname>Gallant</surname> <given-names>C. J.</given-names></name> <name><surname>Marinescu</surname> <given-names>V. D.</given-names></name> <name><surname>Niklasson</surname> <given-names>M.</given-names></name> <name><surname>Segerman</surname> <given-names>A.</given-names></name> <name><surname>Flamourakis</surname> <given-names>G.</given-names></name></person-group> (<year>2016</year>). <article-title>Simultaneous multiplexed measurement of RNA and proteins in single cells</article-title>. <source>Cell Rep.</source> <volume>14</volume>, <fpage>380</fpage>&#x02013;<lpage>389</lpage>.<pub-id pub-id-type="doi">10.1016/j.celrep.2015.12.021</pub-id><pub-id pub-id-type="pmid">26748716</pub-id></citation></ref>
<ref id="B10"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Devaraju</surname> <given-names>N. S. G. K.</given-names></name> <name><surname>Unger</surname> <given-names>M. A.</given-names></name></person-group> (<year>2012</year>). <article-title>Pressure driven digital logic in PDMS based microfluidic devices fabricated by multilayer soft lithography</article-title>. <source>Lab. Chip</source> <volume>12</volume>, <fpage>4809</fpage>&#x02013;<lpage>4815</lpage>.<pub-id pub-id-type="doi">10.1039/c2lc21155f</pub-id><pub-id pub-id-type="pmid">23000861</pub-id></citation></ref>
<ref id="B11"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Devonshire</surname> <given-names>A. S.</given-names></name> <name><surname>Elaswarapu</surname> <given-names>R.</given-names></name> <name><surname>Foy</surname> <given-names>C. A.</given-names></name></person-group> (<year>2011</year>). <article-title>Applicability of RNA standards for evaluating RT-qPCR assays and platforms</article-title>. <source>BMC Genomics</source> <volume>12</volume>:<fpage>118</fpage>.<pub-id pub-id-type="doi">10.1186/1471-2164-12-118</pub-id><pub-id pub-id-type="pmid">21332979</pub-id></citation></ref>
<ref id="B12"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Dey</surname> <given-names>S. S.</given-names></name> <name><surname>Kester</surname> <given-names>L.</given-names></name> <name><surname>Spanjaard</surname> <given-names>B.</given-names></name> <name><surname>Bienko</surname> <given-names>M.</given-names></name> <name><surname>Oudenaarden</surname> <given-names>A.</given-names></name></person-group> (<year>2015</year>). <article-title>Integrated genome and transcriptome sequencing of the same cell</article-title>. <source>Nat. Biotechnol.</source> <volume>33</volume>, <fpage>285</fpage>&#x02013;<lpage>289</lpage>.<pub-id pub-id-type="doi">10.1038/nbt.3129</pub-id><pub-id pub-id-type="pmid">25599178</pub-id></citation></ref>
<ref id="B13"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Ennen</surname> <given-names>M.</given-names></name> <name><surname>Keime</surname> <given-names>C.</given-names></name> <name><surname>Kobi</surname> <given-names>D.</given-names></name> <name><surname>Mengus</surname> <given-names>G.</given-names></name> <name><surname>Lipsker</surname> <given-names>D.</given-names></name> <name><surname>Thibault-Carpentier</surname> <given-names>C.</given-names></name> <etal/></person-group> (<year>2014</year>). <article-title>Single-cell gene expression signatures reveal melanoma cell heterogeneity</article-title>. <source>Oncogene</source> <volume>34</volume>, <fpage>3251</fpage>&#x02013;<lpage>3263</lpage>.<pub-id pub-id-type="doi">10.1038/onc.2014.262</pub-id></citation></ref>
<ref id="B14"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Fan</surname> <given-names>H. C.</given-names></name> <name><surname>Fu</surname> <given-names>G. K.</given-names></name> <name><surname>Fodor</surname> <given-names>S. P. A.</given-names></name></person-group> (<year>2015</year>). <article-title>Combinatorial labeling of single cells for gene expression cytometry</article-title>. <source>Science</source> <volume>347</volume>, <fpage>1258367</fpage>.<pub-id pub-id-type="doi">10.1126/science.1258367</pub-id></citation></ref>
<ref id="B15"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Frei</surname> <given-names>A. P.</given-names></name> <name><surname>Bava</surname> <given-names>F.-A.</given-names></name> <name><surname>Zunder</surname> <given-names>E. R.</given-names></name> <name><surname>Hsieh</surname> <given-names>E. W. Y.</given-names></name> <name><surname>Chen</surname> <given-names>S.-Y.</given-names></name> <name><surname>Nolan</surname> <given-names>G. P.</given-names></name> <etal/></person-group> (<year>2016</year>). <article-title>Highly multiplexed simultaneous detection of RNAs and proteins in single cells</article-title>. <source>Nat. Methods</source> <volume>13</volume>, <fpage>269</fpage>&#x02013;<lpage>275</lpage>.<pub-id pub-id-type="doi">10.1038/nmeth.3742</pub-id><pub-id pub-id-type="pmid">26808670</pub-id></citation></ref>
<ref id="B16"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Gao</surname> <given-names>D.</given-names></name> <name><surname>Vela</surname> <given-names>I.</given-names></name> <name><surname>Sboner</surname> <given-names>A.</given-names></name> <name><surname>Iaquinta</surname> <given-names>P. J.</given-names></name> <name><surname>Karthaus</surname> <given-names>W. R.</given-names></name> <name><surname>Gopalan</surname> <given-names>A.</given-names></name> <etal/></person-group> (<year>2014</year>). <article-title>Organoid cultures derived from patients with advanced prostate cancer</article-title>. <source>Cell</source> <volume>159</volume>, <fpage>176</fpage>&#x02013;<lpage>187</lpage>.<pub-id pub-id-type="doi">10.1016/j.cell.2014.08.016</pub-id><pub-id pub-id-type="pmid">25201530</pub-id></citation></ref>
<ref id="B17"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Guo</surname> <given-names>G.</given-names></name> <name><surname>Pinello</surname> <given-names>L.</given-names></name> <name><surname>Han</surname> <given-names>X.</given-names></name> <name><surname>Lai</surname> <given-names>S.</given-names></name> <name><surname>Shen</surname> <given-names>L.</given-names></name> <name><surname>Lin</surname> <given-names>T.-W.</given-names></name> <etal/></person-group> (<year>2016</year>). <article-title>Serum-based culture conditions provoke gene expression variability in mouse embryonic stem cells as revealed by single-cell analysis</article-title>. <source>Cell Rep.</source> <volume>14</volume>, <fpage>956</fpage>&#x02013;<lpage>965</lpage>.<pub-id pub-id-type="doi">10.1016/j.celrep.2015.12.089</pub-id><pub-id pub-id-type="pmid">26804902</pub-id></citation></ref>
<ref id="B18"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Kim</surname> <given-names>K.-T.</given-names></name> <name><surname>Lee</surname> <given-names>H. W.</given-names></name> <name><surname>Lee</surname> <given-names>H.-O.</given-names></name> <name><surname>Kim</surname> <given-names>S. C.</given-names></name> <name><surname>Seo</surname> <given-names>Y. J.</given-names></name> <name><surname>Chung</surname> <given-names>W.</given-names></name> <etal/></person-group> (<year>2015</year>). <article-title>Single-cell mRNA sequencing identifies subclonal heterogeneity in anti-cancer drug responses of lung adenocarcinoma cells</article-title>. <source>Genome Biol.</source> <volume>16</volume>, <fpage>1</fpage>&#x02013;<lpage>15</lpage>.<pub-id pub-id-type="doi">10.1186/s13059-015-0692-3</pub-id><pub-id pub-id-type="pmid">26084335</pub-id></citation></ref>
<ref id="B19"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Klein</surname> <given-names>A. M.</given-names></name> <name><surname>Mazutis</surname> <given-names>L.</given-names></name> <name><surname>Akartuna</surname> <given-names>I.</given-names></name> <name><surname>Tallapragada</surname> <given-names>N.</given-names></name> <name><surname>Veres</surname> <given-names>A.</given-names></name> <name><surname>Li</surname> <given-names>V.</given-names></name> <etal/></person-group> (<year>2015</year>). <article-title>Droplet barcoding for single-cell transcriptomics applied to embryonic stem cells</article-title>. <source>Cell</source> <volume>161</volume>, <fpage>1187</fpage>&#x02013;<lpage>1201</lpage>.<pub-id pub-id-type="doi">10.1016/j.cell.2015.04.044</pub-id><pub-id pub-id-type="pmid">26000487</pub-id></citation></ref>
<ref id="B20"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Macaulay</surname> <given-names>I. C.</given-names></name></person-group> (<year>2015</year>). <article-title>G&#x00026;T-seq: parallel sequencing of single-cell genomes and transcriptomes</article-title>. <source>Nat. Methods</source> <volume>12</volume>, <fpage>519</fpage>&#x02013;<lpage>522</lpage>.<pub-id pub-id-type="doi">10.1038/nmeth.3370</pub-id></citation></ref>
<ref id="B21"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Macosko</surname> <given-names>E. Z.</given-names></name> <name><surname>Basu</surname> <given-names>A.</given-names></name> <name><surname>Satija</surname> <given-names>R.</given-names></name> <name><surname>Nemesh</surname> <given-names>J.</given-names></name> <name><surname>Shekhar</surname> <given-names>K.</given-names></name> <name><surname>Goldman</surname> <given-names>M.</given-names></name> <etal/></person-group> (<year>2015</year>). <article-title>Highly parallel genome-wide expression profiling of individual cells using nanoliter droplets</article-title>. <source>Cell</source> <volume>161</volume>, <fpage>1202</fpage>&#x02013;<lpage>1214</lpage>.<pub-id pub-id-type="doi">10.1016/j.cell.2015.05.002</pub-id><pub-id pub-id-type="pmid">26000488</pub-id></citation></ref>
<ref id="B22"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Pollen</surname> <given-names>A. A.</given-names></name> <name><surname>Nowakowski</surname> <given-names>T. J.</given-names></name> <name><surname>Chen</surname> <given-names>J.</given-names></name> <name><surname>Retallack</surname> <given-names>H.</given-names></name> <name><surname>Sandoval-Espinosa</surname> <given-names>C.</given-names></name> <name><surname>Nicholas</surname> <given-names>C. R.</given-names></name> <etal/></person-group> (<year>2015</year>). <article-title>Molecular identity of human outer radial glia during cortical development</article-title>. <source>Cell</source> <volume>163</volume>, <fpage>55</fpage>&#x02013;<lpage>67</lpage>.<pub-id pub-id-type="doi">10.1016/j.cell.2015.09.004</pub-id><pub-id pub-id-type="pmid">26406371</pub-id></citation></ref>
<ref id="B23"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Saadatpour</surname> <given-names>A.</given-names></name> <name><surname>Guo</surname> <given-names>G.</given-names></name> <name><surname>Orkin</surname> <given-names>S. H.</given-names></name> <name><surname>Yuan</surname> <given-names>G.-C.</given-names></name></person-group> (<year>2014</year>). <article-title>Characterizing heterogeneity in leukemic cells using single-cell gene expression analysis</article-title>. <source>Genome Biol.</source> <volume>15</volume>, <fpage>1</fpage>&#x02013;<lpage>13</lpage>.<pub-id pub-id-type="doi">10.1186/s13059-014-0525-9</pub-id><pub-id pub-id-type="pmid">25517911</pub-id></citation></ref>
<ref id="B24"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Shalek</surname> <given-names>A. K.</given-names></name> <name><surname>Satija</surname> <given-names>R.</given-names></name> <name><surname>Shuga</surname> <given-names>J.</given-names></name> <name><surname>Trombetta</surname> <given-names>J. J.</given-names></name> <name><surname>Gennert</surname> <given-names>D.</given-names></name> <name><surname>Lu</surname> <given-names>D.</given-names></name> <etal/></person-group> (<year>2014</year>). <article-title>Single-cell RNA-seq reveals dynamic paracrine control of cellular variation</article-title>. <source>Nature</source> <volume>510</volume>, <fpage>363</fpage>&#x02013;<lpage>369</lpage>.<pub-id pub-id-type="doi">10.1038/nature13437</pub-id><pub-id pub-id-type="pmid">24919153</pub-id></citation></ref>
<ref id="B25"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Unger</surname> <given-names>M. A.</given-names></name> <name><surname>Chou</surname> <given-names>H.-P.</given-names></name> <name><surname>Thorsen</surname> <given-names>T.</given-names></name> <name><surname>Scherer</surname> <given-names>A.</given-names></name> <name><surname>Quake</surname> <given-names>S. R.</given-names></name></person-group> (<year>2000</year>). <article-title>Monolithic microfabricated valves and pumps by multilayer soft lithography</article-title>. <source>Science</source> <volume>288</volume>, <fpage>113</fpage>&#x02013;<lpage>116</lpage>.<pub-id pub-id-type="doi">10.1126/science.288.5463.113</pub-id><pub-id pub-id-type="pmid">10753110</pub-id></citation></ref>
<ref id="B26"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Wilson</surname> <given-names>N. K.</given-names></name> <name><surname>Kent</surname> <given-names>D. G.</given-names></name> <name><surname>Buettner</surname> <given-names>F.</given-names></name> <name><surname>Shehata</surname> <given-names>M.</given-names></name> <name><surname>Macaulay</surname> <given-names>I. C.</given-names></name> <name><surname>Calero-Nieto</surname> <given-names>F. J.</given-names></name> <etal/></person-group> (<year>2015</year>). <article-title>Combined single-cell functional and gene expression analysis resolves heterogeneity within stem cell populations</article-title>. <source>Cell Stem Cell</source> <volume>16</volume>, <fpage>712</fpage>&#x02013;<lpage>724</lpage>.<pub-id pub-id-type="doi">10.1016/j.stem.2015.04.004</pub-id><pub-id pub-id-type="pmid">26004780</pub-id></citation></ref>
<ref id="B27"><citation citation-type="journal"><person-group person-group-type="author"><name><surname>Yu</surname> <given-names>M.</given-names></name> <name><surname>Bardia</surname> <given-names>A.</given-names></name> <name><surname>Aceto</surname> <given-names>N.</given-names></name> <name><surname>Bersani</surname> <given-names>F.</given-names></name> <name><surname>Madden</surname> <given-names>M. W.</given-names></name> <name><surname>Donaldson</surname> <given-names>M. C.</given-names></name> <etal/></person-group> (<year>2014</year>). <article-title>Ex vivo culture of circulating breast tumor cells for individualized testing of drug susceptibility</article-title>. <source>Science</source> <volume>345</volume>, <fpage>216</fpage>&#x02013;<lpage>220</lpage>.<pub-id pub-id-type="doi">10.1126/science.1253533</pub-id></citation></ref>
</ref-list>
</back>
</article>
