This work was supported by the Institut Pasteur, INSERM, Universit&#233; Paris Diderot and by the Minist&#232;re de la Recherche (to S.C.); the Association pour la Recherche sur le Cancer (to S.C. and R.G.); the REVIVE Future Investment Program and the Agence Nationale de Recherche (ANR; grant ''Twothyme'' to A.C.); ANR grant ''Myeloten'' (to R.G.); and the Institut National du Cancer (Role of the immune microenvironment during liver carcinogenesis, to R.G.). We acknowledge the Center for Human Immunology and Cytometry platform at Institut Pasteur for their support.

All animal experiments were approved of by the Pasteur Institute Safety Committee in accordance with the French Agriculture Ministry and the EU guidelines.
1. Prepare a 96-well Single-cell Sorting Plate
<ol>
	<li>Prepare pre-amplification mix in a 1.5 ml tube by adding 5.0 &#181;l of specific retro-transcription buffer, 1.3 &#181;l of low ethylenediaminetetraacetic acid (EDTA) TE buffer (10 mM TE solution, pH 8, and 0.1 mM EDTA solution; 0.2 &#181;M filtered), and 0.2 &#181;l of Taq DNA polymerase per well (Figure 1a).</li>
	<li>On a 96-well plate, distribute 6.5 &#181;l of pre-amplification mix in each well (up to 48 wells). This plate is the 96-well single-cell sorting plate. 
	NOTE: Using an electronic pipette is recommended for this protocol. It allows for precise and reproducible volume measurements and saves time.</li>
</ol>
<img alt="Figure 1" src="/files/ftp_upload/54858/54858fig1.jpg" /> 
Figure 1: 96-well plate loading procedure. The loading of the 96-well plates should be done very carefully, as it will impact the rest of the experiment. (a) On the 96-well single-cell sorting plate, the pre-amplification mix is distributed first, followed by the 0.2x assay mix. (b) The record of each single-cell position should be kept on a spreadsheet. A well without a cell is called a &#34;no input&#34; well and can be used as a control. Two rows can be spared to control primer efficiency with cDNA dilution (sequential one-in-ten dilutions from the equivalent of 105 cells to one cell). (c) On the 96-well assay plate, the assay loading reagent is distributed first, followed by the addition of the primers. Do not forget to keep a layout of each primer position. <a href="https://cloudfront-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858fig1large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<ol start="3">
	<li>Prepare a 0.2x assay mix in a 1.5 ml tube by adding 1.4 &#181;l of each primer (primers 20x; up to 48 different probes; see Table 1). Adjust the final volume to 140 &#181;l with low EDTA TE buffer.</li>
	<li>Distribute 2.5 &#181;l of the 0.2x assay mix to the 48 wells in the 96-well single-cell sorting plate containing the pre-amplification mix. 
	NOTE: The final volume in each well of the 96-well sample plate should be 9 &#181;l.</li>
	<li>Seal the plate with a cover film. Vortex the plate and spin it down at 280 x g for 45 sec. Store the 96-well single-cell sorting plate at -20 &#176;C or use it immediately.</li>
</ol>
2. Single Cell Dissociation
<ol>
	<li>Using a CO2 delivery system, euthanize a mouse by hypoxia. Fix the mouse on a dissection plate and open the peritoneal cavity longitudinally using scissors.</li>
	<li>Recover ILCs from the liver.
	<ol>
		<li>Prior to liver isolation, flush the liver-circulating cells.
		<ol>
			<li>To bring out the portal vein, move the small and the large intestine to the left side of the mouse; the portal vein appears as the largest vein in the peritoneal cavity.</li>
			<li>With a 10 ml syringe and a 0.45 mm needle, flush the liver with phosphate-buffered saline medium (PBS) through the portal vein. To facilitate liver flushing, cut the gastric artery with scissors.</li>
			<li>To remove the gallbladder, pinch the base of the gallbladder with forceps and separate it from the liver with scissors.</li>
			<li>To isolate the liver, pinch the base of the liver with forceps, and using scissors, separate the liver from the rest of the peritoneal cavity.</li>
			<li>Transfer the liver in a potter tube with 5 ml of Roswell Park Memorial Institute medium (RPMI) 2% fetal calf serum (FCS).</li>
		</ol>
		</li>
	</ol>
	</li>
	<li>Dissociate the liver with a mechanic dissociation using a pestle. Transfer the medium to a 15 ml tube. Bring the total volume to 14 ml with RPMI 2% FCS. Leave the cell suspension on ice for 15-20 min to decant the hepatocytes.</li>
	<li>Recover the supernatant (without the hepatocytes) in a new 15 ml tube using a 1 ml pipette (12 ml of supernatant can be expected). Spin down at 120 x g for 7 min at 10 &#176;C. Discard the supernatant after centrifugation.</li>
	<li>Resuspend the pellet obtained in step 2.4 in 14 ml of density gradient medium (e.g., Percoll 40%). Spin down at 600 x g for 20 min at 20 &#176;C. Remove the supernatant by aspiration.</li>
	<li>Resuspend the pellet in 1 ml of potassium acetate (ACK). Leave at room temperature for 1 min. After 1 min, bring the total volume to 14 ml with Hank&#39;s Balanced Salted Solution (HBSS) 2% FCS. Spin down at 120 x g for 7 min at 10 &#176;C. Discard the supernatant.</li>
	<li>Stain the cells.
	<ol>
		<li>Resuspend the pellet in 300 &#181;l of biotinylated antibody mix (see the Materials Table; add all biotinylated antibodies mentioned) and transfer the cells in to 1.5 ml tube. Leave it in the dark for 20 min at 4 &#176;C.</li>
		<li>Bring the total volume to 1 ml with HBSS 2% FCS. Spin down at 120 x g for 7 min at 10 &#176;C. Resuspend the pellet in 300 &#181;l of fluorochrome-coupled antibody mix and fluorochrome-conjugated streptavidin (see the Materials Table; add all fluorochrome-coupled antibodies mentioned). Leave it in the dark for 20 min at 4 &#176;C.</li>
	</ol>
	</li>
	<li>Bring the total volume to 1 ml with HBSS 2% FCS. Spin down at 120 x g for 7 min at 10 &#176;C. Remove the supernatant using a 1 ml pipette. Resuspend the pellet in 1 ml HBSS and propidium iodide (Pi; 1:4,000). 
	NOTE: When using the widely expressed ILC cell-surface markers, 3 populations are defined: NKp46+ IL-7R&#945;-, NKp46+ IL-7R&#945;+, and NKp46- IL-7R&#945;+. All populations are sorted lineage- CD3- CD4- CD45.2+.</li>
</ol>
3. Single-cell Fluorescence-activated Cell Sorting (FACS)
<ol>
	<li>Use FACS to sort single cells13. 
	NOTE: Different cell types can be sorted on the same 96-well single-cell sorting plate (Figure 2).
	<ol>
		<li>Use FACS to sort one cell per well containing the 0.2x assay mix and the pre-amplification mix on the 96-well single-cell sorting plate. Use a freshly prepared 96-well single-cell sorting plate or thaw it if stored at -20 &#176;C.
		<ol>
			<li>Put a sealed 96-well plate on the FACS plate holder. Use an empty 96-well plate as a test. Position the plate with the A1-well on the left and toward the experimenter.</li>
			<li>Adjust the plate holder to obtain a drop in the center of the A1-well with verification beads.</li>
			<li>Repeat step 3.1.1.2 with all the wells in line A.</li>
			<li>Remove the seal from the 96-well plate and sort 100 verification beads per well. Check for drop formation in the center of the wells.</li>
			<li>When properly adjusted, place the 96-well single-cell sorting plate on the plate holder. Draw the plate layout and sort 1 cell per well. 
			NOTE: The proper positioning of each cell on the 96-well single-cell sorting plate is essential. A plate layout should be kept on a spreadsheet software. The plate layout will be used in the next several steps (Figure 1b)14. Leave one well containing 0.2x assay mix and pre-amplification mix on the 96-well sample plate without cells. This well is used as a no-input control. Leave two rows of 6 wells for a cDNA dilution. These wells are used as controls for primer efficiency (Figure 1b).</li>
			<li>Store the 96-well single-cell sorting plate at -20 &#176;C or use it immediately.</li>
		</ol>
		</li>
	</ol>
	</li>
</ol>
4. Pre-amplification
<ol>
	<li>Use a fresh or thawed 96-well single-cell sorting plate obtained in step 3.1.1.6. Vortex the plate and spin it down (280 x g for 45 sec).</li>
	<li>Place the 96-well single-cell sorting plate on the thermocycler. Perform reverse transcription and pre-amplification as per the mentioned program in Figure 3 and Table 2.</li>
</ol>
<img alt="Figure 3" src="/files/ftp_upload/54858/54858fig3.jpg" /> 
Figure 3: Pre-amplification program. In order to have enough material, pre-amplification of specific target genes on sorted single cells is required. The 96-well single-cell sorting plate is loaded on a thermocycler to follow the pre-amplification program. The pre-amplification products are then diluted with low EDTA TE buffer and can be used immediately or frozen at -20 &#176;C. <a href="https://cloudfront-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858fig3large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<ol start="3">
	<li>Dilute the pre-amplified samples by adding 36 &#181;l of low EDTA TE buffer into each well.</li>
	<li>Seal the plate with a cover film. Vortex the plate and spin it down (280 g for 45 sec). Store the pre-amplified, 96-well single-cell sorting plate at -20 &#176;C or use it immediately.</li>
</ol>
5. Prepare a 96-well Sample Plate
<ol>
	<li>Prepare 191 &#181;l of master mix by adding 175 &#181;l of qPCR master mix and 17.5 &#181;l of sample loading reagent.</li>
	<li>On a new 96-well plate, distribute 3.6 &#181;l of master mix in each well (up to 48 wells). This is the 96-well sample plate.</li>
	<li>Transfer 2.9 &#181;l of pre-amplified cDNA from the 96-well sample plate to the new 96-well plate. Keep the same position for each sample between the 96-single cell sorting plate and the 96-well sample plate.</li>
	<li>Seal the plate with a cover film. Vortex the plate and spin it down (280 x g for 45 sec). Place the new 96-well plate on ice protected from light. Store the 96-well sample plate at -20 &#176;C or use it immediately.</li>
</ol>
6. Prepare a 96-well Assay Plate
<ol>
	<li>On a new 96-well plate, distribute 3 &#181;l of assay loading reagent in each well (up to 48 wells). This plate is called the 96-well assay plate.</li>
	<li>Add 3 &#181;l of primers to each well. Keep a layout on a spreadsheet software (Figure 1c and Table 1; use all primers listed in Table 1). The proper positioning of each primer on the 96-well assay plate is essential for the next several steps.
	<ol>
		<li>If the number of samples is below 48, add water instead of primers to unused wells.</li>
	</ol>
	</li>
	<li>Seal the plate with a cover film. Vortex the plate and spin it down (280 x g for 45 sec). Store the 96-well assay plate at -20 &#176;C or use it immediately. 
	NOTE: To avoid repeated freezing and thawing of the primers, distribute the primers for the pre-amplification mix and for the 96-well assay plate at the same time.</li>
</ol>
7. Single-cell Gene-expression Chip
<ol>
	<li>Place the integrated microfluidic circuit (IFC) on the bench and check the valves using the capped syringe. Open the syringe, place it perpendicular to the valve, and press firmly. The O-ring should move. Fill the chip with control line fluid (0.3 ml of tuberculin).</li>
	<li>Repeat step 7.1 with the second valve. 
	NOTE: No liquid should be spilled over the chip.</li>
	<li>Remove the blue film from the bottom of the chip. Load the chip into the IFC controller. On the IFC controller screen, select &#34;PRIME&#34; and then &#34;RUN.&#34; The control of the microfluidic lines last 10 min.</li>
	<li>Eject the chip and re-seal the blue film on the bottom of the chip.</li>
</ol>
<img alt="Figure 4" src="/files/ftp_upload/54858/54858fig4.jpg" /> 
Figure 4: Single-cell multiplex gene-expression chip loading. These steps require great precision, especially during the transfer of the 96-well plate to the single-cell multiplex gene-expression chip. To avoid loading errors and misplacements, it is highly recommended to work sequentially. The volume taken for each transfer should be controlled during the pipetting process. Finally, it is important to avoid any bubbles and to remove them in case of formation. <a href="https://cloudfront-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858figfourlarge.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<ol start="5">
	<li>Using an 8-channel pipette, transfer 5 &#181;l from the 96-well assay plate on the A1 side (left side) of the chip. Fill the left side of the chip, as indicated in Figure 4. Be careful to always draw the same volume into the tips.
	<ol>
		<li>Do not create bubbles. If bubbles appear, remove them using 10 &#181;l tips. Change tips for each well of the chip.</li>
	</ol>
	</li>
	<li>Repeat step 7.5 on the right side of the chip using 5 &#181;l from the 96-well sample plate.</li>
	<li>Remove the blue film from the bottom of the chip. Load the chip into the IFC controller. On the IFC controller screen, select &#34;RUN SCRIPT&#34; and then &#34;RUN.&#34; The loading of the microfluidic lines lasts 45 min.</li>
	<li>Eject the chip and re-seal the blue film on the bottom of the chip.</li>
</ol>
8. Run the Chip
<ol>
	<li>On the microfluidic qPCR computer, select the &#34;Data Collection&#34; software. Once it starts, select &#34;New Run.&#34;</li>
	<li>Select &#34;Eject,&#34; remove the blue film from the bottom of the chip, and load the chip. Place the chip to get the A1 well on the left and toward the experimenter.</li>
	<li>On &#34;Project setting,&#34; select &#34;None&#34; and then &#34;Next.&#34;</li>
	<li>Select &#34;New chip run,&#34; &#34;New chip directory&#34; (create a folder for the experiment), and &#34;Next.&#34;</li>
	<li>On &#34;Gene expression&#34; select, &#34;Passive reference = ROX&#34; (carboxy-x-rhodamine) and &#34;Single probe;&#34; select the correct filter according to the primers. Select &#34;Next.&#34;</li>
	<li>Select &#34;Browse&#34; and select the correct program according to the primers used. 
	NOTE: &#34;Default-10min-Hotstart.pcl&#34; is the most commonly used program for single-cell multiplex RT-qPCR.</li>
	<li>Select &#34;Start Run&#34;. The reaction takes approximately 90 min.</li>
	<li>Once finished, select &#34;Eject-Done.&#34;</li>
</ol>
9. Data Analysis
<ol>
	<li>Open the &#34;Real-Time PCR Analysis&#34; software, select &#34;File&#34; and &#34;Open,&#34; find the experiment folder, and select &#34;ChipRun.bml file.&#34;</li>
	<li>Click on &#34;Analysis View,&#34; &#34;Results Table,&#34; and &#34;Heat Map View;&#34; boxes marked with an &#34;X&#34; are below the threshold detection level and/or had bad amplification curves.</li>
	<li>Name the samples. Go to &#34;Sample Setup&#34; and select &#34;New SBS 96.&#34; Click on &#34;Mapping&#34; and select &#34;&#8230;&#34; Copy and paste the sample layout design from the spreadsheet software. Define the pasted layout as &#34;Sample Name.&#34;</li>
	<li>Name the assays. Repeat step 9.3 with the primer names in &#34;Detector Setup.&#34; Define the pasted layout as &#34;Detector Name.&#34;</li>
	<li>Click on &#34;Analysis view&#34; and &#34;Analyze;&#34; the sample and primer names will be attributed automatically to the samples and primers (Figure 7).</li>
	<li>Calculate Ct, Delta Ct, and Fold Change for each sample. Use a housekeeping gene as an internal control (here, Gapdh): 
	&#916;Ct = Ctsample - Ctinternal control 
	Fold Change = 2-(&#916;Ctsample-&#916;Ctnegative control)
	<ol>
		<li>Click on &#34;Detector Setup&#34; and select the well with no input control (use the housekeeping gene as an internal control to normalize gene expression). Select &#34;Editor,&#34; &#34;define reference,&#34; and &#34;Update.&#34;</li>
		<li>Select all wells to which the above-defined internal control will be applied. Select &#34;Editor&#34; and &#34;define as Test.&#34; Select an appropriate reference gene and click &#34;Update.&#34;</li>
		<li>Select &#34;Sample Setup&#34; to select the negative control well. Select &#34;Editor&#34; and define &#34;Reference.&#34; Click &#34;Update.&#34;</li>
		<li>Select &#34;Analysis View&#34; and &#34;Analyze;&#34; the Delta Ct and Fold Change will automatically be calculated.</li>
	</ol>
	</li>
	<li>Export the data. Select &#34;File&#34; and &#34;Export, save as &#34;Table Results,&#34; select the destination folder, and hit &#34;Save&#34; and &#34;Exit.&#34;</li>
	<li>Repeat step 9.7, but instead of &#34;Table Results,&#34; save as &#34;HeatMap Results.&#34;</li>
</ol>

<ol><li>Sun, H., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((A+Bead-Based+Microfluidic+Approach+to+Integrated+Single-Cell+Gene+Expression+Analysis+by+Quantitative+RT-PCR%5BTitle%5D)+AND+%22RSC+Adv%22%5BJournal%5D)">A Bead-Based Microfluidic Approach to Integrated Single-Cell Gene Expression Analysis by Quantitative RT-PCR.</a> RSC Adv. 5, 4886-4893 (2015).</li>
<li>Kellogg, R. A., Berg, R. G., Leyrat, A. A., Tay, S. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((High-throughput+microfluidic+single-cell+analysis+pipeline+for+studies+of+signaling+dynamics%5BTitle%5D)+AND+%22Nat.+Protoc%22%5BJournal%5D)">High-throughput microfluidic single-cell analysis pipeline for studies of signaling dynamics.</a> Nat. Protoc. 9 (7), 1713-1726 (2014).</li>
<li>Moignard, V., Macaulay, I. C., Swiers, G., Buettner, F., Schütte, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Characterisation+of+transcriptional+networks+in+blood+stem+and+progenitor+cells+using+high-throughput+single+cell+gene+expression+analysis%5BTitle%5D)+AND+%22Nat.+Cell+Biol%22%5BJournal%5D)">Characterisation of transcriptional networks in blood stem and progenitor cells using high-throughput single cell gene expression analysis.</a> Nat. Cell Biol. 15 (4), 363-372 (2013).</li>
<li>Chea, S., Schmutz, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Single-Cell+Gene+Expression+Analyses+Reveal+Heterogeneous+Responsiveness+of+Fetal+Innate+Lymphoid+Progenitors+to+Notch+Signaling%5BTitle%5D)+AND+%22Cell+Rep%22%5BJournal%5D)">Single-Cell Gene Expression Analyses Reveal Heterogeneous Responsiveness of Fetal Innate Lymphoid Progenitors to Notch Signaling.</a> Cell Rep. 14 (6), 1500-1516 (2016).</li>
<li>Ishizuka, I. E., Chea, S., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Single-cell+analysis+defines+the+divergence+between+the+innate+lymphoid+cell+lineage+and+lymphoid+tissue+-+inducer+cell+lineage%5BTitle%5D)+AND+%22Nat.+Immunol%22%5BJournal%5D)">Single-cell analysis defines the divergence between the innate lymphoid cell lineage and lymphoid tissue - inducer cell lineage.</a> Nat. Immunol. 17, 1-9 (2016).</li>
<li>Spits, H., Artis, D., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Innate+lymphoid+cells+-+a+proposal+for+uniform+nomenclature%5BTitle%5D)+AND+%22Nat.Rev.+Immunol%22%5BJournal%5D)">Innate lymphoid cells - a proposal for uniform nomenclature.</a> Nat.Rev. Immunol. 13 (2), 145-149 (2013).</li>
<li>Walker, J. A., Barlow, J. L., McKenzie, A. N. J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Innate+lymphoid+cells-+how+did+we+miss+them%5BTitle%5D)+AND+%22Nat.+Rev.+Immunol%22%5BJournal%5D)">Innate lymphoid cells- how did we miss them.</a> Nat. Rev. Immunol. 13 (2), 75-87 (2013).</li>
<li>Vallentin, B., Barlogis, V., et al. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pubmed/26438443">Innate Lymphoid Cells in Cancer.</a> Cancer Immunol. Res. 3 (10), 1109-1114 (2015).</li>
<li>Monticelli, L. a, Artis, D. <a target="_blank" href="https://www.ncbi.nlm.nih.gov/pmc/articles/PMC3320042/">Innate lymphoid cells promote lung tissue homeostasis following acute influenza virus infection.</a> Nat. Immunol. 12 (11), 1045-1054 (2012).</li>
<li>Tang, L., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Differential+phenotypic+and+functional+properties+of+liver-resident+NK+cells+and+mucosal+ILC1s%5BTitle%5D)+AND+%22J.+Autoimmun%22%5BJournal%5D)">Differential phenotypic and functional properties of liver-resident NK cells and mucosal ILC1s.</a> J. Autoimmun. 67, 29-35 (2015).</li>
<li>Durum, K., Gilfillan, S., Colonna, M., Consortium, G. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Transcriptional+Programs+Define+Molecular+Characteristics+of+Innate+Lymphoid+Cell+Classes+and+Subsets%5BTitle%5D)+AND+%22Nat.+Immunol%22%5BJournal%5D)">Transcriptional Programs Define Molecular Characteristics of Innate Lymphoid Cell Classes and Subsets.</a> Nat. Immunol. 16 (3), 306-317 (2015).</li>
<li>Wang, H., Ramakrishnan, A., Fletcher, S., Prochownik, E. V., Genetics, M. Clonal dynamics of native haematopoiesis. Nature. 2 (2), 322-327 (2015).</li>
<li>Samsa, L. A., Fleming, N., Magness, S., Qian, L., Liu, J. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Isolation+and+Characterization+of+Single+Cells+from+Zebrafish+Embryos%5BTitle%5D)+AND+%22J.+Vis.+Exp%22%5BJournal%5D)">Isolation and Characterization of Single Cells from Zebrafish Embryos.</a> J. Vis. Exp. (109), e1-e10 (2016).</li>
<li>Klose, C., Flach, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Differentiation+of+Type+1+ILCs+from+a+Common+Progenitor+to+All+Helper-like+Innate+Lymhpoid+Cell+Lineages%5BTitle%5D)+AND+%22Cell%22%5BJournal%5D)">Differentiation of Type 1 ILCs from a Common Progenitor to All Helper-like Innate Lymhpoid Cell Lineages.</a> Cell. 157, 340-356 (2014).</li>
<li>Dicle, Y., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Bioinformatics+Approaches+to+Single-Cell+Analysis+in+Developmental+Biology%5BTitle%5D)+AND+%22MHR%22%5BJournal%5D)">Bioinformatics Approaches to Single-Cell Analysis in Developmental Biology.</a> MHR. 22, 182-192 (2016).</li>
<li>Setty, M., et al. <a target="_blank" href="http://www.ncbi.nlm.nih.gov/pubmed?term=((Wishbone+identifies+bifurcating+developmental+trajectories+from+single-cell+data%5BTitle%5D)+AND+%22Nat+Biotechnol%22%5BJournal%5D)">Wishbone identifies bifurcating developmental trajectories from single-cell data.</a> Nat Biotechnol. 34, 637-645 (2016).</li>
</ol>

The authors have nothing to disclose.

This protocol describes how to obtain exhaustive gene expression information at the clonal level. Here, we investigated liver ILC compartment heterogeneity. After single-cell sorting of different ILC populations (based on widely expressed ILC surface markers), samples were pre-amplified for specific pre-selected genes. Then, the obtained cDNA and primers were loaded onto a multiplex RT-qPCR microfluidic chip. Finally, we obtained the expression of 48 different genes from 48 single cells. Gene expression results were analyzed via an online software to investigate gene signature and cell population relationships.
The principal advantage of this technique is the ability to assess multiple gene expression simultaneously. It also allows work at the clonal level and on very rare populations. Unlike conventional RT-qPCR methods, single-cell multiplex RT-qPCR has no averaging effect, since gene expression is assessed at the clonal level. Thus, heterogeneity within a population can be detected and analyzed. Single-cell RT-qPCR expression allows work on very rare populations, since only very few cells are needed to obtain wide information on gene expression. Moreover, different populations of cells can be tested on the same plate. This enables the detection of differences in gene signature between populations and, with online data tools, the assessment of cell population relationships. This technique presents other advantages as well. Single-cell multiplex RT-qPCR has the benefit of the recent microfluidic advances. Volumes needed in the IFC chambers do not exceed the nanoliter-level. Thus, lower amounts of reagents are used, reducing the cost of experiments. Finally, the overall number of pipetting steps is reduced, thereby limiting the possibility of pipetting errors. The steps are simultaneous for all the cells and can be done with a multichannel pipette, enabling work in a fast and time-saving manner. Obtained results at the clonal level are reproducible, reliable, and can be used to perform robust data analyses.
Single-cell multiplex RT-qPCR does, however, present some limitations. First, this technique necessitates expensive and specific equipment, such as a microfluidic chip, a microfluidic chip controller, and a specific thermocycler. Unlike conventional RT-qPCR methods, this technique is time- and reagent-saving, since, with conventional RT-qPCR, numerous studies are needed to obtain statistically significant comparisons. In this protocol, we present a procedure to obtain gene expression for 48 genes. 96-96 multiplex RT-qPCR microfluidic chips are available. These chips can assess gene expression for 96 genes from 96 cells at the same time. During cell sorting, a minimum number of starting cells are needed, since precision and cell purity must be maximal. However, even with a small number of cells, it is still possible to sort for single-cell multiplex RT-qPCR gene expression (step 3). Finally, single-cell multiplex RT-qPCR is a sensitive method. Different tests must be performed before starting such experiments. For instance, for primer target specificity, primers should be tested for amplification efficiency and for the relative competition for the target.
Several steps of the protocol must be performed with caution. The manipulations of small volumes combined with a large number of samples and assays at the same time increase the risk of pipetting mistakes (steps 1, 4, 5, 6, and 7). All the reagents and mixes should be vortexed prior to pipetting in order to ensure a homogenous solution. Evaporation during pre-amplification (step 4) can also impact the final results. The plates should always be properly sealed with an appropriate cover film. Sorting strategy must be very precise (step 3) to avoid dead cells or multiple cells in wells within the plates. FACS sorting machines are now equipped with single-cell sorting index software that can record which specific cell was sorted in each well. This new tool will be of interest for the determination of whether gene signatures are correlated to cell-surface marker intensity. Sample and assay positions on chips must be organized meticulously, as this is essential for the experiment (steps 1, 5, and 6). Finally, the selection of primers (step 1 and 6) is a critical step. Every primer must be selected based on previously described results or expected results. A random choice of primers in the set of assessed genes could impair the final readout of the experiment. Furthermore, primers that correlate with the cell-surface marker used for cell sorting have to be used as controls with housekeeping genes.
Assessing multiplex gene expression at the clonal level offers a better characterization of the liver ILC heterogeneous compartment than in previous research. This technique, supplied by online software, describes more precisely the different subsets of ILC in the liver at homeostasis than studies using only cell-surface markers. Furthermore, it is possible to consider cell population relationships and to build differentiation networks. The gene signatures, based on single-cell gene expression, are also important tools to discern cell population features and to examine potential roles. Single-cell multiplex RT-qPCR is one of the best techniques for a gene expression assay to obtain reliable data. These data are easily manageable with bioinformatics software15. With respect to other techniques for gene expression studies, such as microarrays or RNA sequencing, single-cell multiplex RT-qPCR offers high sensitivity. Other approaches to single-cell analysis have been described and can be used as complementary techniques to single-cell multiplex RT-qPCR15,16. Furthermore, the technique can be performed with any primer combination, allowing users to look at personalized gene signatures. Many other applications can be used with this technique. For instance, time-course gene expression of cells throughout development can be done to assess the modification of the molecular profile during development. Drug effects on cells can also be investigated with drug dilution to determine the threshold of drug efficiency on gene expression.

Over the past few years, innate lymphoid cells (ILCs) have been increasingly investigated. Despite their lack of antigen-specific receptors, they belong to the lymphoid lineage and represent important sentinels for tissue homeostasis and inflammation. ILCs are currently divided into three groups based on their expression of specific transcription factor combinations and on their ability to produce cytokines6,7.
ILCs contribute to numerous homeostatic and pathophysiological situations in diverse organs via specific cytokine production8,9. To be able to understand the role of these cells, it is important to determine the various ILC subpopulations per organ and to identify their developmental relationships. In addition, phenomena of plasticity between the different subsets have been related. By studying the heterogeneity of the cells present in one organ, it is possible to delimit their stage of maturation and to distinguish their specific functions.
To illustrate the technique of single-cell multiplex RT-qPCR, hepatic ILCs were chosen, with a particular emphasis on their heterogeneity within the same ILC group (type 1 ILC)10. First, through the use of flow cytometry, three distinct ILC populations were characterized in the liver. Group 1 ILC represents around 80% of the innate effectors, while the two other populations are rare hepatic ILC populations (less than 5% of the innate effectors). Those populations were sorted using widely expressed cell-surface markers of ILC populations. As a result, sorted ILC populations in the liver look broadly similar one to another.
Single-cell multiplex RT-qPCR has emerged as one of the best techniques to promptly investigate the heterogeneity of these populations11. Two main characteristics are determined by taking advantage of the single-cell multiplex RT-qPCR technique. First, by looking at the clonal level, it is possible to recover cell-specific gene expression for comparison between cells that apparently display similar developmental stages. Then, by looking at a pre-selected combination of gene expression, we will determine new gene signatures based on simultaneous gene expression patterns at one time point. These aspects permit the collection of a wide variety of expression data for a large number of cells, even on rare populations, since the technique is performed at the clonal level. Thereby, ILC heterogeneity in the liver can be adequately assessed.
Next, by sorting all cells with a global ILC phenotype, a wide overview of the multiple-gene expression of the liver ILC populations is obtained, even though they represent extremely rare populations. A microfluidic-based chip allows experimentation with even a small amount of cell material. As a consequence, the gene expression profiles of rare cell populations can be obtained. Using online gene signature analysis software, cell population clusters and potential cell relationships can be investigated. Consequently, functional tests can be performed to validate the clustering data at the in vivo level.
Tens of gene expressions could be assessed concomitantly on hundreds or more single cells on the same chip3,11,12. Design of the assay is the longest and most important part of the experiment. The determination of the genes relevant to the hypothesis to be tested is paramount to obtain relevant results. Secondly, internal control (such as known surface markers used for sorting) and specific controls are needed. This is crucial to test the primer amplification specificity, the efficiency of the amplification, and the absence of primer competition. Therefore, working with single-cell multiplex RT-qPCR is a timesaving technique, as multiple-gene expression of a cell is assessed at the same time.
Using the same chip and mix of reagents for all cells limits the possible errors of manipulation and allows for reproducibility between samples. Altogether, the different aspects of single-cell multiplex RT-qPCR allows for the production of highly reliable results at the clonal level, with a great level of sensitivity for a wide variety of samples and genes. The obtained results offer powerful and robust data for biostatistical tests.
This can be achieved due to the microfluidic aspect of the method, which allows for work on very small amounts of material and leads to exhaustive results. Finally, using online software, it is possible to compare the desired populations.

<table><tbody><tr><td>Cells Direct One Step qRT-PCR kit</td><td>Applied Biosystems&amp;nbsp;</td><td>11753100</td><td>Primer probe detection kit.Contains 2x reaction mix, SSIII Platinium enzyme.</td></tr><tr><td>Low TE EDTA Buffer</td><td>Affymetrix</td><td>75793 100ML</td><td></td></tr><tr><td>96-well plates</td><td>Thermofisher Scientific</td><td>AB 1100</td><td>96-well plates adapted for cell sorting and thermocycling</td></tr><tr><td>Cover film</td><td>Dominique Dutscher</td><td>106570</td><td>aluminium cover film; avoid contamination and evaporation</td></tr><tr><td>Actb</td><td>Thermofisher Scientific</td><td>Mm00607939_s1</td><td>20x primer</td></tr><tr><td>Aes</td><td>Thermofisher Scientific</td><td>Mm01148854_s1</td><td>20x primer</td></tr><tr><td>Ahr</td><td>Thermofisher Scientific</td><td>Mm00478932_s1</td><td>20x primer</td></tr><tr><td>Bcl2</td><td>Thermofisher Scientific</td><td>Mm00477631_s1</td><td>20x primer</td></tr><tr><td>c-myc</td><td>Thermofisher Scientific</td><td>Mm00487804_s1</td><td>20x primer</td></tr><tr><td>Cbfb</td><td>Thermofisher Scientific</td><td>Mm01251026_s1</td><td>20x primer</td></tr><tr><td>Cd27</td><td>Thermofisher Scientific</td><td>Mm01185212_s1</td><td>20x primer</td></tr><tr><td>Cd49a</td><td>Thermofisher Scientific</td><td>Mm01306375_s1</td><td>20x primer</td></tr><tr><td>CD49b</td><td>Thermofisher Scientific</td><td>Mm00434371_s1</td><td>20x primer</td></tr><tr><td>Cxcr5</td><td>Thermofisher Scientific</td><td>Mm00432086_s1</td><td>20x primer</td></tr><tr><td>Cxcr6</td><td>Thermofisher Scientific</td><td>Mm02620517_s1</td><td>20x primer</td></tr><tr><td>Eomes</td><td>Thermofisher Scientific</td><td>Mm01351985_s1</td><td>20x primer</td></tr><tr><td>Ets1</td><td>Thermofisher Scientific</td><td>Mm01175819_s1</td><td>20x primer</td></tr><tr><td>Foxo1</td><td>Thermofisher Scientific</td><td>Mm00490672_s1</td><td>20x primer</td></tr><tr><td>Gapdh</td><td>Thermofisher Scientific</td><td>Mm03302249_s1</td><td>20x primer</td></tr><tr><td>Gata3</td><td>Thermofisher Scientific</td><td>Mm00484683_s1</td><td>20x primer</td></tr><tr><td>Gm-csf</td><td>Thermofisher Scientific</td><td>Mm01136644_s1</td><td>20x primer</td></tr><tr><td>Hes1</td><td>Thermofisher Scientific</td><td>Mm01342805_s1</td><td>20x primer</td></tr><tr><td>Hprt</td><td>Thermofisher Scientific</td><td>Mm00446968_s1</td><td>20x primer</td></tr><tr><td>Id2</td><td>Thermofisher Scientific</td><td>Mm01293217_s1</td><td>20x primer</td></tr><tr><td>Il-12rb2</td><td>Thermofisher Scientific</td><td>Mm00711781_s1</td><td>20x primer</td></tr><tr><td>Il-18r1</td><td>Thermofisher Scientific</td><td>Mm00515178_s1</td><td>20x primer</td></tr><tr><td>Il-1rl1</td><td>Thermofisher Scientific</td><td>Mm00434237_s1</td><td>20x primer</td></tr><tr><td>Il-22</td><td>Thermofisher Scientific</td><td>Mm001226722_s1</td><td>20x primer</td></tr><tr><td>Il-23r</td><td>Thermofisher Scientific</td><td>Mm00519943_s1</td><td>20x primer</td></tr><tr><td>Il-2ra</td><td>Thermofisher Scientific</td><td>Mm01340213_s1</td><td>20x primer</td></tr><tr><td>Il-2rb</td><td>Thermofisher Scientific</td><td>Mm01195267_s1</td><td>20x primer</td></tr><tr><td>IL-7r</td><td>Thermofisher Scientific</td><td>Mm00434295_s1</td><td>20x primer</td></tr><tr><td>Klr5</td><td>Thermofisher Scientific</td><td>Mm04207528_s1</td><td>20x primer</td></tr><tr><td>Lef1</td><td>Thermofisher Scientific</td><td>Mm00550265_s1</td><td>20x primer</td></tr><tr><td>Ncr1</td><td>Thermofisher Scientific</td><td>Mm01337324_s1</td><td>20x primer</td></tr><tr><td>Nfil3</td><td>Thermofisher Scientific</td><td>Mm01339838_s1</td><td>20x primer</td></tr><tr><td>Notch1</td><td>Thermofisher Scientific</td><td>Mm00435249_s1</td><td>20x primer</td></tr><tr><td>Notch2</td><td>Thermofisher Scientific</td><td>Mm00803069_s1</td><td>20x primer</td></tr><tr><td>Rora</td><td>Thermofisher Scientific</td><td>Mm01173766_s1</td><td>20x primer</td></tr><tr><td>Rorc</td><td>Thermofisher Scientific</td><td>Mm01261022_s1</td><td>20x primer</td></tr><tr><td>Runx3</td><td>Thermofisher Scientific</td><td>Mm00490666_s1</td><td>20x primer</td></tr><tr><td>Tbx21</td><td>Thermofisher Scientific</td><td>Mm01299453_s1</td><td>20x primer</td></tr><tr><td>Tcf3</td><td>Thermofisher Scientific</td><td>Mm01175588_s1</td><td>20x primer</td></tr><tr><td>Tcf7</td><td>Thermofisher Scientific</td><td>Mm00493445_s1</td><td>20x primer</td></tr><tr><td>Tle1</td><td>Thermofisher Scientific</td><td>Mm00495643_s1</td><td>20x primer</td></tr><tr><td>Tle3</td><td>Thermofisher Scientific</td><td>Mm00437097_s1</td><td>20x primer</td></tr><tr><td>Tsc22d3</td><td>Thermofisher Scientific</td><td>Mm01306210_s1</td><td>20x primer</td></tr><tr><td>Tnfrsf11a</td><td>Thermofisher Scientific</td><td>Mm00437132_s1</td><td>20x primer</td></tr><tr><td>Tox</td><td>Thermofisher Scientific</td><td>Mm00455231_s1</td><td>20x primer</td></tr><tr><td>Zbtb16</td><td>Thermofisher Scientific</td><td>Mm01176868_s1</td><td>20x primer</td></tr><tr><td>Zbtb7b</td><td>Thermofisher Scientific</td><td>Mm00784709_s1</td><td>20x primer</td></tr><tr><td>qPCR Master mix&amp;nbsp;</td><td>Applied BioSystems</td><td>P/N 4304437</td><td></td></tr><tr><td>2x Assay Loading Reagent</td><td>Fluidigm</td><td>P/N 85000736</td><td>Specific density medium to load assays in multiplex RT-qPCR microfluidic chip.&amp;nbsp;</td></tr><tr><td>2x Sample Loading Reagent</td><td>Fluidigm</td><td>P/N 85000735</td><td>Specific density medium to load samples in multiplex RT-qPCR microfluidic chip.&amp;nbsp;</td></tr><tr><td>48.48 mutliplex RT qPCR microfluidic chip</td><td>Fluidigm</td><td>BMK-M-48.48</td><td>48.48 Dynamic Array IFC for Gene Expression;chip for single cell multiplex RT-qPCR reaction</td></tr><tr><td>48.48 mutliplex RT qPCR microfluidic chip controller</td><td>Fluidigm</td><td>89000020</td><td>48.48 IFC Controller; control the chip internal fluidic system, load samples and assays in reaction chambers</td></tr><tr><td>mutliplex RT qPCR microfluidic thermocycler</td><td>Fluidigm</td><td>GE48.48</td><td>48.48 Dynamic Array IFC thermocycler</td></tr><tr><td>96-well plates</td><td>Thermofisher Scientific</td><td>AB 1100</td><td>96-well plates adapted for cell sorting and thermocycling</td></tr><tr><td>C57Bl/6 mice</td><td>Janvier</td><td>C57Bl/6 miceJ@RJ</td><td></td></tr><tr><td>10 ml syringe</td><td>BD Biosciences</td><td>309639</td><td></td></tr><tr><td>PBS</td><td>Life Technologies</td><td>14040174</td><td></td></tr><tr><td>HBSS</td><td>Life Technologies</td><td>24020133</td><td></td></tr><tr><td>RPMI</td><td>Life Technologies</td><td>61870044</td><td></td></tr><tr><td>FCS</td><td>CVFSVF000U</td><td>Eurobio Abcys</td><td>Standard fetal calf serum</td></tr><tr><td>Potter tube</td><td>N/A</td><td></td><td></td></tr><tr><td>15 ml tube</td><td>Corning</td><td>352097</td><td></td></tr><tr><td>1.5 ml tube</td><td>Sigma-Aldrich</td><td>T9661-1000EA</td><td></td></tr><tr><td>FACS machine</td><td>N/A</td><td></td><td></td></tr><tr><td>centrifuge&amp;nbsp;</td><td>Thermofisher Scientific</td><td>75004538</td><td></td></tr><tr><td>1,000 &amp;micro;l tips</td><td>Fisher Scientific</td><td>10313272</td><td></td></tr><tr><td>P1000&amp;nbsp;</td><td>Gilson</td><td>F123602</td><td></td></tr><tr><td>Percoll</td><td>Dominique Dutscher</td><td>17-0891-01</td><td></td></tr><tr><td>Facs tube</td><td>Falcon</td><td>352235</td><td></td></tr><tr><td>anti-CD8 Biotin mouse antibody</td><td>Sony</td><td>1103520</td><td>lineage antibody</td></tr><tr><td>anti-CD19 Biotin mouse antibody</td><td>Sony</td><td>1177520</td><td>lineage antibody</td></tr><tr><td>anti-TCRab Biotin mouse antibody</td><td>BioLegend</td><td>109204</td><td>lineage antibody</td></tr><tr><td>anti-TCRgd Biotin mouse antibody</td><td>BD Biosciences</td><td>553176</td><td>lineage antibody</td></tr><tr><td>anti-Ter119 Biotin mouse antibody</td><td>BD Biosciences</td><td>553672</td><td>lineage antibody</td></tr><tr><td>anti-Gr1 Biotin mouse antibody</td><td>BD Biosciences</td><td>553125</td><td>lineage antibody</td></tr><tr><td>anti-CD45.2 PerCPCy5.5 mouse antibody</td><td>BioLegend</td><td>109828</td><td></td></tr><tr><td>anti-IL7ra PeCy7 mouse antibody</td><td>ebioSciences</td><td>25-1271-82</td><td></td></tr><tr><td>anti-CD3 BV510 mouse antibody</td><td>BD Biosciences</td><td>563024</td><td></td></tr><tr><td>anti-CD4 BV786 mouse antibody</td><td>BD Biosciences</td><td>563727</td><td></td></tr><tr><td>anti-NKp46 PE mouse antibody</td><td>ebioSciences</td><td>12-3351-82</td><td></td></tr><tr><td>Streptavidin</td><td>Sony</td><td>2626025</td><td></td></tr><tr><td>Propidium Iodide</td><td>Sigma-Aldrich</td><td>P4864-10ML</td><td></td></tr><tr><td>Electronic pipette</td><td>Eppendorf</td><td>4986000017</td><td></td></tr><tr><td>Combitips 0.1 ml</td><td>Eppendorf</td><td>30089405</td><td></td></tr><tr><td>Multichannel pipette</td><td>Rainin</td><td>L8-10XLS+</td><td></td></tr><tr><td>Accudrop</td><td>BD Biosciences</td><td>345249</td><td>verification beads for FACS</td></tr></tbody></table>

single-cell gene expression using multiplex rt-qpcr to characterize heterogeneity of rare lymphoid populations

Gene expression heterogeneity is an interesting feature to investigate in lymphoid populations. Gene expression in these cells varies during cell activation, stress, or stimulation. Single-cell multiplex gene expression enables the simultaneous assessment of tens of genes1,2,3. At the single-cell level, multiplex gene expression determines population heterogeneity4,5. It allows for the distinction of population heterogeneity by determining both the probable mix of diverse precursor stages among mature cells and also the diversity of cell responses to stimuli.
Innate lymphoid cells (ILC) have been recently described as a population of innate effectors of the immune response6,7. In this protocol, cell heterogeneity of the ILC hepatic compartment is investigated during homeostasis.
Currently, the most widely used technique to assess gene expression is RT-qPCR. This method measures gene expression only one gene at a time. Additionally, this method cannot estimate heterogeneity of gene expression, since multiple cells are needed for one test. This leads to the measurement of the average gene expression of the population. When assessing large numbers of genes, RT-qPCR becomes a time-, reagent-, and sample-consuming method. Hence, the trade-offs limit the number of genes or cell populations that can be evaluated, increasing the risk of missing the global picture.
This manuscript describes how single-cell multiplex RT-qPCR can be used to overcome these limitations. This technique has benefited from recent microfluidics technological advances1,2. Reactions occurring in multiplex RT-qPCR chips do not exceed the nanoliter-level. Hence, single-cell gene expression, as well as simultaneous multiple gene expression, can be performed in a reagent-, sample-, and cost-effective manner. It is possible to test cell gene signature heterogeneity at the clonal level between cell subsets within a population at different developmental stages or under different conditions4,5. Working on rare populations with large numbers of conditions at the single-cell level is no longer a restriction.

Lymphoid populations display great diversity in gene expression. In this protocol, liver ILC compartment heterogeneity was investigated using single-cell multiplex RT-qPCR gene expression. Unlike other gene expression techniques, single-cell multiplex RT-qPCR gene expression allows work on several populations, even the rarest, at the same time. This specificity, coupled with a high sensitivity at the clonal level, allows for the investigation of differences in gene signatures within a population and between rare populations. As an example, the gene signatures of 3 ILC populations from the livers of 4-week-old mice were assessed. One population is frequent enough to be easily identified, whereas the two other are rare (less than 1% of lineage-negative cells) in the liver.
A 96-well single-cell sorting plate and a 96-well assay plate were prepared prior to liver dissection and dissociation (Figure 1). Each mix was prepared under a sterile hood and distributed with an electronic pipette for volume precision and reproducibility. All reagents were vortexed before pipetting to maintain homogeneity. After preparation, the 96-well cell sorting plate can be frozen until cell sorting, and the 96-well cell assay plate can be frozen until IFC loading.
Livers were dissected and dissociated into single-cell suspensions. Cells were stained and sorted on thawed 96-well single-cell sorting plates. Using FACS, 3 ILC populations were sorted based on widely expressed ILC markers (Figure 2). After cell sorting, cDNA was synthesized, and specific target genes were amplified (Figure 3 and Table 2). After the pre-amplification step, cDNA was diluted in low EDTA TE buffer. Primers and cDNA from single cells were loaded onto the IFC chip (Figures 4 and 5a). The obtained results (Figure 7a) were simplified for easier reading and analysis (Figure 7b). Examples of properly or improperly loaded chips are shown (Figure 5). A FAM figure of a properly loaded chip (Figure 6) with different amplification signals is shown.
The results show that after proper cell sorting, pre-amplification, and loading, the ILC population appears heterogeneous for gene expression in the liver of adult wild-type mice (Figure 7). Using online software, cell-specific gene expression signatures (Figures 7 and 8) and cell population relationships (Figure 8) could be identified. Each sorted population has a specific gene signature enriched in gene expression. For instance, cells sorted as NKp46- IL-7R&#945;+ have enriched expression of Rorc, the main transcription factor of group 3 ILCs; Rora, a Rorc-homologous transcription factor; and Il-23r and Cxcr6, two important receptors of group 3 ILC hepatic markers in the liver. The same observations can be made with NKp46+ IL-7R&#945;- and NKp46+ IL-7R&#945;+ populations. Each cell is checked for housekeeping gene expression (Hprt, Act, and Gapdh) to exclude invalid wells; at least two housekeeping genes must be expressed to consider the values as valid.
Interestingly, this technique allows for the definition of two subpopulations of NKp46+ IL-7R&#945;- based on gene signatures. Differences between those two signatures have to be validated by other sets of different functional experiments.
<img alt="Figure 2" src="/files/ftp_upload/54858/54858fig2.jpg" /> 
Figure 2: FACS cell strategy. Images are representative of isolated liver cells from 4-week-old mice. (a) FSC-A/SSC-A gating, (b) FSC-H/FSC-W doublet discrimination, (c, d, and e) alive, lineage-negative CD45.2+ CD4- CD3- cells gating. (f) Using widely expressed ILC cell-surface markers, we defined 3 populations: NKp46+ IL-7R&#945;-, NKp46+ IL-7R&#945;+, and NKp46- IL-7R&#945;+. <a href="https://cloudflare-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858fig2large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 5" src="/files/ftp_upload/54858/54858fig5.jpg" /> 
Figure 5: ROX figures of properly (a) and improperly (b) loaded chips. (a) Properly loaded single-cell multiplex gene-expression chip should appear with straight lines and rows. Each reaction chamber is filled and has the same dimension. (b) Improperly loaded single-cell multiplex gene-expression chip. Empty lines and rows of reaction chambers appear (green square), as well as bending lines (blue square). These features can be due to improper &#34;PRIME&#34; or &#34;LOAD&#34; steps, as well as to residual bubbles in the IFC wells. <a href="https://cloudflare-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858fig5large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 6" src="/files/ftp_upload/54858/54858fig6.jpg" /> 
Figure 6: FAM (Fluorescein amidite) figure of a properly loaded chip. After a few cycles (depending on the samples and on the primers assessed), differences in reaction chamber brightness should appear. Reaction chambers with an amplification signal should appear brighter (blue square) than reaction chambers with no or low amplification signals (green square). <a href="https://cloudflare-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858fig6large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 7" src="/files/ftp_upload/54858/54858fig7.jpg" /> 
Figure 7: Heat map obtained before (a) and after (b) modifications. (a) Direct heat map obtained without analysis or modification of assay/sample order. (b) Modified heat map obtained after sample and assay name definition. No input wells (blue square) are obtained if no amplification occurs. Inefficient primers (green square). cDNA dilution test (purple square). <a href="https://cloudflare-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858fig7large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<img alt="Figure 8" src="/files/ftp_upload/54858/54858fig8.jpg" /> 
Figure 8: Single-cell population clustering. A data analysis software is used to determine the clustering between populations. Free clustering software are available online. By assessing the whole expression profile of a single cell, the software can display gene signatures and cell population relationships. Each square represents a cell: blue are NKp46+ IL-7R&#945;-, red are NKp46+ IL-7R&#945;+, and green are NKp46- IL-7R&#945;+. Each population has a specific gene signature. Within the NKp46+ IL-7R&#945;- population, two signatures can be seen, attesting to population heterogeneity. <a href="https://cloudflare-jove-com.remotexs.ntu.edu.sg/files/ftp_upload/54858/54858fig8large.jpg" target="_blank">Please click here to view a larger version of this figure.</a>
<table border="1" fo:keep-together.within-page="1" fo:keep-with-next.within-page="always">
	<tbody>
		<tr>
			<td colspan="2">Primers</td>
		</tr>
		<tr>
			<td>Actb</td>
			<td>Il-23r</td>
		</tr>
		<tr>
			<td>Aes</td>
			<td>Il-2ra</td>
		</tr>
		<tr>
			<td>Ahr</td>
			<td>Il-2rb</td>
		</tr>
		<tr>
			<td>Bcl2</td>
			<td>II-7ra</td>
		</tr>
		<tr>
			<td>c-myc</td>
			<td>Klr5</td>
		</tr>
		<tr>
			<td>Cbfb</td>
			<td>Lef1</td>
		</tr>
		<tr>
			<td>Cd27</td>
			<td>Ncr1</td>
		</tr>
		<tr>
			<td>Cd49a</td>
			<td>Nfil3</td>
		</tr>
		<tr>
			<td>Cd49b</td>
			<td>Notch1</td>
		</tr>
		<tr>
			<td>Cxcr5</td>
			<td>Notch2</td>
		</tr>
		<tr>
			<td>Cxcr6</td>
			<td>Rora</td>
		</tr>
		<tr>
			<td>Eomes</td>
			<td>Rorc</td>
		</tr>
		<tr>
			<td>Ets1</td>
			<td>Runx3</td>
		</tr>
		<tr>
			<td>Foxo1</td>
			<td>Tbx21</td>
		</tr>
		<tr>
			<td>Gapdh</td>
			<td>Tcf3</td>
		</tr>
		<tr>
			<td>Gata3</td>
			<td>Tcf7</td>
		</tr>
		<tr>
			<td>Gm-csf</td>
			<td>Tle1</td>
		</tr>
		<tr>
			<td>Hes1</td>
			<td>Tle3</td>
		</tr>
		<tr>
			<td>Hprt</td>
			<td>Tsc22d3</td>
		</tr>
		<tr>
			<td>Id2</td>
			<td>Tnfrsf11a</td>
		</tr>
		<tr>
			<td>Il-12rb2</td>
			<td>Tox</td>
		</tr>
		<tr>
			<td>Il-18r1</td>
			<td>Zbtb16</td>
		</tr>
		<tr>
			<td>Il-1rl1</td>
			<td>Zbtb7b</td>
		</tr>
		<tr>
			<td>Il-22</td>
			<td></td>
		</tr>
	</tbody>
</table>
Table 1: Primer list.
<img alt="Table 2" src="/files/ftp_upload/54858/54858table2.jpg" /> 
Table 2: Pre-amplification program. Each step of the pre-amplification program is shown with the complete information regarding the time, temperature, and number of cycles.

Watch this Scientific Journal Video about Single-cell Gene Expression Using Multiplex RT-qPCR to Characterize Heterogeneity of Rare Lymphoid Populations at JoVE.com

Single-cell Gene Expression Using Multiplex RT-qPCR to Characterize Heterogeneity of Rare Lymphoid Populations

This protocol describes how to assess the expression of a large array of genes at the clonal level. Single-cell RT-qPCR produces highly reliable results with a strong sensitivity for hundreds of samples and genes.

Gene expression heterogeneity is an interesting feature to investigate in lymphoid populations. Gene expression in these cells varies during cell ...

single-cell-gene-expression-using-multiplex-rt-qpcr-to-characterize

Nanyang Technological University

Research

JoVE Journal

Genetics

10.9K Views.  Institut Pasteur.  Dit protocol beschrijft hoe de expressie van een groot scala van genen op het klonale niveau te beoordelen. Single-cell RT-qPCR levert zeer betrouwbare resultaten met een sterke gevoeligheid voor honderden monsters en genen. 

Video: Single-cell Gene Expression Using Multiplex RT-qPCR to Characterize Heterogeneity of Rare Lymphoid Populations

High Throughput MicroRNA Profiling: Optimized Multiplex qRT-PCR at Nanoliter Scale on the Fluidigm Dynamic ArrayTM IFCs

The broad involvement of miRNAs in critical processes underlying development, tissue homoeostasis and disease has led to a surging interest among the research and pharmaceutical communities. To study miRNAs, it is essential that the quantification of microRNA levels is accurate and robust. By comparing wild-type to small RNA deficient mouse embryonic stem cells (mESC), we revealed a lack of accuracy and robustness in previous published multiplex qRT-PCR techniques. Here, we describe an optimized method, including purifying away excessive primers from previous multiplex steps before singleplex real time detection, which dramatically increases the accuracy and robustness of the technique. Furthermore, we explain how performing the technique on a microfluidic chip at nanoliter volumes significantly reduces reagent costs and permits time effective high throughput miRNA expression profiling.

Here we describe an optimized multiplex reverse transcriptase quantitative PCR (qRT-PCR) protocol in combination with a microfluidic platform as a cost and time effective high-throughput screening tool for microRNA (miRNA) expression levels, especially when working with limited amounts of sample.

The broad involvement of miRNAs in critical processes underlying development, tissue homoeostasis and disease has led to a surging interest among the ...

Bioengineering

Single-cell Gene Expression Profiling Using FACS and qPCR with Internal Standards

Gene expression measurements from bulk populations of cells can obscure the considerable transcriptomic variation of individual cells within those populations. Single-cell gene expression measurements can help assess the role of noise in gene expression, identify correlations in the expression of pairs of genes, and reveal subpopulations of cells that respond differently to a stimulus. Here, we describe a procedure to measure the expression of up to 96 genes in single mammalian cells isolated from a population growing in tissue culture. Cells are sorted into lysis buffer by fluorescence-activated cell sorting (FACS), and the mRNA species of interest are reverse-transcribed and amplified. Gene expression is then measured using a microfluidic real-time PCR machine, which performs up to 96 qPCR assays on up to 96 samples at a time. We also describe the generation and use of PCR amplicon standards to enable the estimation of the absolute number of each transcript. Compared with other methods of measuring gene expression in single cells, this approach allows for the quantification of more distinct transcripts than RNA FISH at a lower cost than RNA-Seq.

We describe a method to sort single mammalian cells and to quantify the expression of up to 96 target genes of interest in each cell. This method includes the use of internal qPCR standards to enable the estimation of absolute transcript counts.

Gene expression measurements from bulk populations of cells can obscure the considerable transcriptomic variation of individual cells within those ...

A Combinatorial Single-cell Approach to Characterize the Molecular and Immunophenotypic Heterogeneity of Human Stem and Progenitor Populations

Immunophenotypic characterization and molecular analysis have long been used to delineate heterogeneity and define distinct cell populations. FACS is inherently a single-cell assay, however prior to molecular analysis, the target cells are often prospectively isolated in bulk, thereby losing single-cell resolution. Single-cell gene expression analysis provides a means to understand molecular differences between individual cells in heterogeneous cell populations. In bulk cell analysis an overrepresentation of a distinct cell type results in biases and occlusions of signals from rare cells with biological importance. By utilizing FACS index sorting coupled to single-cell gene expression analysis, populations can be investigated without the loss of single-cell resolution while cells with intermediate cell surface marker expression are also captured, enabling evaluation of the relevance of continuous surface marker expression. Here, we describe an approach that combines single-cell reverse transcription quantitative PCR (RT-qPCR) and FACS index sorting to simultaneously characterize the molecular and immunophenotypic heterogeneity within cell populations.
In contrast to single-cell RNA sequencing methods, the use of qPCR with specific target amplification allows for robust measurements of low-abundance transcripts with fewer dropouts, while it is not confounded by issues related to cell-to-cell variations in read depth. Moreover, by directly index-sorting single-cells into lysis buffer this method, allows for cDNA synthesis and specific target pre-amplification to be performed in one step as well as for correlation of subsequently derived molecular signatures with cell surface marker expression. The described approach has been developed to investigate hematopoietic single-cells, but have also been used successfully on other cell types.
In conclusion, the approach described herein allows for sensitive measurement of mRNA expression for a panel of pre-selected genes with the possibility to develop protocols for subsequent prospective isolation of molecularly distinct subpopulations.

Bulk gene expression measurements cloud individual cell differences in heterogeneous cell populations. Here, we describe a protocol for how single-cell gene expression analysis and index sorting by Florescence Activated Cell Sorting (FACS) can be combined to delineate heterogeneity and immunophenotypically characterize molecularly distinct cell populations.

Immunophenotypic characterization and molecular analysis have long been used to delineate heterogeneity and define distinct cell populations. FACS is ...

Combining Laser Capture Microdissection and Microfluidic qPCR to Analyze Transcriptional Profiles of Single Cells: A Systems Biology Approach to Opioid Dependence

Profound transcriptional heterogeneity in anatomically adjacent single cells suggests that robust tissue functionality may be achieved by cellular phenotype diversity. Single-cell experiments investigating the network dynamics of biological systems demonstrate cellular and tissue responses to various conditions at biologically meaningful resolution. Herein, we explain our methods for gathering single cells from anatomically specific locations and accurately measuring a subset of their gene expression profiles. We combine laser capture microdissection (LCM) with microfluidic reverse transcription quantitative polymerase chain reactions (RT-qPCR). We also use this microfluidic RT-qPCR platform to measure the microbial abundance of gut contents.

This protocol explains how to collect single neurons, microglia, and astrocytes from the central nucleus of the amygdala with high accuracy and anatomic specificity using laser capture microdissection. Additionally, we explain our use of microfluidic RT-qPCR to measure a subset of the transcriptome of these cells.

Profound transcriptional heterogeneity in anatomically adjacent single cells suggests that robust tissue functionality may be achieved by cellular ...

Multiplexed Single Cell mRNA Sequencing Analysis of Mouse Embryonic Cells

Single cell mRNA sequencing has made significant progress in the last several years and has become an important tool in the field of developmental biology. It has been successfully used to identify rare cell populations, discover novel marker genes, and decode spatial and temporal developmental information. The single cell method has also evolved from the microfluidic based Fluidigm C1 technology to the droplet-based solutions in the last two to three years. Here we used the heart as an example to demonstrate how to profile the mouse embryonic tissue cells using the droplet based scRNA-Seq method. In addition, we have integrated two strategies into the workflow to profile multiple samples in a single experiment. Using one of the integrated methods, we have simultaneously profiled more than 9,000 cells from eight heart samples. These methods will be valuable to the developmental biology field by providing a cost-effective way to simultaneously profile single cells from different genetic backgrounds, developmental stages, or anatomical locations.

Here we presented a multiplexed single cell mRNA sequencing method to profile gene expression in mouse embryonic tissues. The droplet-based single cell mRNA sequencing (scRNA-Seq) method in combination with multiplexing strategies can profile single cells from multiple samples simultaneously, which significantly reduces reagent costs and minimizes experimental batch effects.

Single cell mRNA sequencing has made significant progress in the last several years and has become an important tool in the field of developmental ...

Developmental Biology

Single Cell Transcriptional Profiling of Adult Mouse Cardiomyocytes

While numerous studies have examined gene expression changes from homogenates of heart tissue, this prevents studying the inherent stochastic variation between cells within a tissue. Isolation of pure cardiomyocyte populations through a collagenase perfusion of mouse hearts facilitates the generation of single cell microarrays for whole transcriptome gene expression, or qPCR of specific targets using nanofluidic arrays. We describe here a procedure to examine single cell gene expression profiles of cardiomyocytes isolated from the heart. This paradigm allows for the evaluation of metrics of interest which are not reliant on the mean (for example variance between cells within a tissue) which is not possible when using conventional whole tissue workflows for the evaluation of gene expression (Figure 1). We have achieved robust amplification of the single cell transcriptome yielding micrograms of double stranded cDNA that facilitates the use of microarrays on individual cells. In the procedure we describe the use of NimbleGen arrays which were selected for their ease of use and ability to customize their design. Alternatively, a reverse transcriptase - specific target amplification (RT-STA) reaction, allows for qPCR of hundreds of targets by nanofluidic PCR. Using either of these approaches, it is possible to examine the variability of expression between cells, as well as examining expression profiles of rare cell types from within a tissue. Overall, the single cell gene expression approach allows for the generation of data that can potentially identify idiosyncratic expression profiles that are typically averaged out when examining expression of millions of cells from typical homogenates generated from whole tissues.

Single cell expression profiling allows the detailed gene expression analysis of individual cells. We describe methods for the isolation of cardiomyocytes, and preparing the resulting lysates for either whole transcriptome microarray or qPCR of specific targets.

While numerous studies have examined gene expression changes from homogenates of heart tissue, this prevents studying the inherent stochastic variation ...

Biology

Single Cell Multiplex Reverse Transcription Polymerase Chain Reaction After Patch-clamp

The cerebral cortex is composed of numerous cell types exhibiting various morphological, physiological, and molecular features. This diversity hampers easy identification and characterization of these cell types, prerequisites to study their specific functions. This article describes the multiplex single cell reverse transcription polymerase chain reaction (RT-PCR) protocol, which allows, after patch-clamp recording in slices, to detect simultaneously the expression of tens of genes in a single cell. This simple method can be implemented with morphological characterization and is widely applicable to determine the phenotypic traits of various cell types and their particular cellular environment, such as in the vicinity of blood vessels. The principle of this protocol is to record a cell with the patch-clamp technique, to harvest and reverse transcribe its cytoplasmic content, and to detect qualitatively the expression of a predefined set of genes by multiplex PCR. It requires a careful design of PCR primers and intracellular patch-clamp solution compatible with RT-PCR. To ensure a selective and reliable transcript detection, this technique also requires appropriate controls from cytoplasm harvesting to amplification steps. Although precautions discussed here must be strictly followed, virtually any electrophysiological laboratory can use the multiplex single cell RT-PCR technique.

This protocol describes the critical steps and precautions required to perform single cell multiplex reverse transcription polymerase chain reaction after patch-clamp. This technique is a simple and effective method to analyze the expression profile of a predetermined set of genes from a single cell characterized by patch-clamp recordings.

The cerebral cortex is composed of numerous cell types exhibiting various morphological, physiological, and molecular features. This diversity hampers ...

Neuroscience

Determination of Immune Cell Identity and Purity Using Epigenetic-Based Quantitative PCR

Immune cell subtype population frequencies can have a large effect on the efficacy of T cell therapies. Current methods, like flow cytometry, have specific sample requirements, high sample input, are low throughput, and are difficult to standardize, all of which are detrimental to characterization of cell therapy products during their development and manufacturing.
The assays described herein accurately identify and quantify immune cell types in a heterogeneous mixture of cells using isolated genomic DNA (gDNA). DNA methylation patterns are revealed through bisulfite conversion, a process in which unmethylated cytosines are converted to uracils. Unmethylated DNA regions are detected through qPCR amplification using primers targeting converted areas. One unique locus per assay is measured and serves as an accurate identifier for a specific cell type. The assays are robust and identify CD8+, regulatory, and Th17 T cells in a high throughput manner. These optimized assays can potentially be used for in-process and product release testing for cell therapy process.

Here we describe a robust method of determining immune cell identity and purity through epigenetic signatures detected using quantitative PCR (qPCR). DNA demethylation at a specific locus serves as a unique identifier for a particular cell type and allows for identification of CD8+, regulatory, or Th17 T cells.

Immune cell subtype population frequencies can have a large effect on the efficacy of T cell therapies. Current methods, like flow cytometry, have ...

Immunology and Infection

Profiling Individual Human Embryonic Stem Cells by Quantitative RT-PCR

Heterogeneity of stem cell population hampers detailed understanding of stem cell biology, such as their differentiation propensity toward different lineages. A single cell transcriptome assay can be a new approach for dissecting individual variation. We have developed the single cell qRT-PCR method, and confirmed that this method works well in several gene expression profiles. In single cell level, each human embryonic stem cell, sorted by OCT4::EGFP positive cells, has high expression in OCT4, but a different level of NANOG expression. Our single cell gene expression assay should be useful to interrogate population heterogeneities.

Single cell gene expression assay is needed for understanding stem cell heterogeneities.

Heterogeneity of stem cell population hampers detailed understanding of stem cell biology, such as their differentiation propensity toward different ...

Multiplexed Analysis of Retinal Gene Expression and Chromatin Accessibility Using scRNA-Seq and scATAC-Seq

Powerful next generation sequencing techniques offer robust and comprehensive analysis to investigate how retinal gene regulatory networks function during development and in disease states. Single-cell RNA sequencing allows us to comprehensively profile gene expression changes observed in retinal development and disease at a cellular level, while single-cell ATAC-Seq allows analysis of chromatin accessibility and transcription factor binding to be profiled at similar resolution. Here the use of these techniques in the developing retina is described, and MULTI-Seq is demonstrated, where individual samples are labeled with a modified oligonucleotide-lipid complex, enabling researchers to both increase the scope of individual experiments and substantially reduce costs.

Here, the authors showcase the utility of MULTI-seq for phenotyping and subsequent paired scRNA-seq and scATAC-seq in characterizing the transcriptomic and chromatin accessibility profiles in retina.

Powerful next generation sequencing techniques offer robust and comprehensive analysis to investigate how retinal gene regulatory networks function during ...

Single-cell Gene Expression Using Multiplex RT-qPCR to Characterize Heterogeneity of Rare Lymphoid Populations

Samenvatting

Meer video's verkennen