In situ combinatorial labeling of cellular molecules

Abstract

To provide methods of uniquely labeling RNA molecules within a plurality of cells.SOLUTION: Methods of uniquely labeling or barcoding molecules within a nucleus, a plurality of nuclei, a cell, a plurality of cells, and / or a tissue are provided. Kits for uniquely labeling or barcoding molecules within a nucleus, a plurality of nuclei, a cell, a plurality of cells, and / or a tissue are also provided. The molecules to be labeled may include, but are not limited to, RNAs and / or cDNAs.SELECTED DRAWING: Figure 20

Description

[Technical field]

[0001] CROSS-REFERENCE TO RELATED APPLICATIONS This application claims the benefit of U.S. Provisional Patent Application No. 62 / 561,806, filed Sep. 22, 2017, which is incorporated by reference in its entirety.

[0002] Technical Field The present disclosure generally relates to a method for uniquely labeling or barcoding molecules in a nucleus, a plurality of nuclei, a cell, a plurality of cells, and / or a tissue. The present disclosure also relates to a kit for uniquely labeling molecules in a nucleus, a plurality of nuclei, a cell, a plurality of cells, and / or a tissue. In particular, the method and kit may relate to labeling RNA and / or cDNA. [Background technology]

[0003] background Next-generation sequencing (NGS) can be used to identify and / or quantify individual transcripts from a sample of cells. However, such techniques may be too complex to perform on individual cells in large samples. Such methods generally involve purifying RNA transcripts from lysed cells (i.e., cells that have been broken into pieces) and then using reverse transcription to convert the RNA transcripts into complementary DNA (cDNA). The cDNA sequences can then be sequenced using NGS. In such procedures, the cDNA sequences are all mixed together before sequencing, so RNA expression of the entire sample is measured and individual sequences cannot be linked to the individual cells of origin. Summary of the Invention [Problem to be solved by the invention]

[0004] Methods for uniquely labeling or barcoding transcripts from individual cells may involve manually separating individual cells into separate reaction vessels and may require specialized equipment. An alternative approach to sequencing individual transcripts of cells is to use microscopy to identify individual fluorescent bases. However, this technique can be difficult to perform and limited to sequencing a small number of cells.

[0005] The presently disclosed embodiments will become more fully apparent from the following description and appended claims, taken in conjunction with the accompanying drawings. [Means for solving the problem]

[0006] Detailed Description The present disclosure generally relates to methods for uniquely labeling or barcoding molecules in a nucleus, multiple nuclei, a cell, multiple cells, and / or a tissue. The present disclosure also relates to kits for uniquely labeling or barcoding molecules in a nucleus, multiple nuclei, a cell, multiple cells, and / or a tissue. The molecules to be labeled can include, but are not limited to, RNA, cDNA, DNA, proteins, peptides, and / or antigens. [Brief description of the drawings]

[0007] [Figure 1] FIG. 1 shows ligation of nucleic acid tags to form labels or barcodes. [Diagram 2]Schematic diagram of cDNA formation by in situ reverse transcription. Panel A shows cells that have been fixed and permeabilized. Panel B shows the addition of a poly(T) primer that can template the reverse transcription of polyadenylated transcripts. Panel C shows the addition of random hexamers that can template the reverse transcription of virtually any transcript. Panel D shows the addition of primers that are designed to target specific transcripts so that only a subset of transcripts can be amplified. Panel E shows the cells of panel A after reverse transcription, showing that the cDNA has hybridized to RNA. [Figure 3A] FIG. 1 shows non-templated ligation of single-stranded adapters to RNA fragments. [Figure 3B] FIG. 1 shows ligation of a single-stranded adapter using a partial duplex with a random hexamer primer. [Figure 4] FIG. 1 shows primer binding. [Diagram 5] FIG. 1 shows reverse transcription after primer binding. [Figure 6] FIG. 1 shows DNA tagged antibodies used for labeling cellular proteins. [Figure 7] FIG. 1 shows aptamers used for labeling cellular proteins. [Figure 8] FIG. 1 is a schematic diagram of cell splitting, tagging, and pooling according to an embodiment of the present disclosure. As shown, cells can be split into multiple reaction vessels. One cell is highlighted to show its path through the illustrated process. [Figure 9A] FIG. 1 illustrates an exemplary workflow according to an embodiment of the present disclosure. [Figure 9B] FIG. 13 illustrates an exemplary workflow according to another embodiment of the present disclosure. [Figure 10] FIG. 1 shows a reverse transcription primer (BC_0055) according to an embodiment of the present disclosure. [Figure 11] FIG. 1 shows annealed first round barcode oligos according to an embodiment of the present disclosure. [Figure 12] FIG. 1 shows annealed second round barcode oligos according to an embodiment of the present disclosure. [Figure 13] FIG. 1 shows annealed third round barcode oligos according to an embodiment of the present disclosure. [Figure 14] FIG. 1 shows a ligation stop oligo according to an embodiment of the present disclosure. [Figure 15] FIG. 1 shows a single stranded DNA adaptor oligo (BC_0047) ligated to the 3′ end of a cDNA according to an embodiment of the disclosure. [Figure 16] FIG. 1 shows PCR products formed using primers BC_0051 and BC_0062, and the 3′ adaptor oligo (BC_0047) after ligation to the barcoded cDNA. [Figure 17] BC_0027, which contains a flow cell binding sequence and a binding site for the TRUSEQ™ read1 primer, and BC_0063, which contains a flow cell binding sequence and a TruSeq multiplex read2 and index binding sequence. Also shown in Figure 17 is the region of the sample index, which in this embodiment is GATCTG. [Figure 18] A scatter plot in which the number of reads aligned to the human genome (X-axis) and the number of reads aligned to the mouse genome (Y-axis) are plotted for each unique barcode combination. [Figure 19] FIG. 1 shows the size of cDNA before and after tagmentation. [Figure 20] Figure 1 shows an experiment in which more than 150,000 nuclei from P2 and P11 mouse brains and spinal cords were profiled in a single experiment in which more than 6 million barcode combinations were used. By recording which of the four starting samples (P2 spinal cord, P2 brain, P11 spinal cord, or P11 spinal cord) was added to each well, the first round barcode sequence can be used to identify which cells / nuclei came from which sample. [Figure 21]Figure 1 shows the number of nuclei in each well during the three rounds of barcoding. Despite manual pipetting of the cells, most wells contain approximately equal numbers of nuclei. Isolation of P2 spinal cords resulted in fewer cells than the other samples, which explains the lower number of nuclei in the corresponding first round wells. [Figure 22] FIG. 1 shows the distribution of P2 / P11 brain and spinal cord single nucleus transcriptomes within the oligodendrocyte glial lineage. The identity of each nucleus sample (P2 spinal cord, P2 brain, P11 spinal cord, or P11 spinal cord) can be determined by examining the first round barcodes derived from each cell / nucleus barcode combination. [Figure 23] FIG. 1 shows an overview of the protocol provided herein. [Figure 24] FIG. 1 shows a molecular diagram of the protocol provided herein. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0008] It will be readily understood that the embodiments generally described herein are exemplary. The following more detailed description of various embodiments is not intended to limit the scope of the present disclosure, but is merely representative of various embodiments. Furthermore, a person skilled in the art may change the order of steps or operations of the methods disclosed herein without departing from the scope of the present disclosure. In other words, the order or use of specific steps or operations may be modified, unless a specific order of steps or operations is necessary for the proper working of the embodiment.

[0009] The term "bond" is used broadly throughout this disclosure to refer to any form of attachment or coupling of two or more components, entities, or objects. For example, two or more components may be bonded to each other by chemical bonds, covalent bonds, ionic bonds, hydrogen bonds, electrostatic forces, Watson-Crick hybridization, and the like.

[0010] One aspect of the present disclosure relates to a method for labeling a nucleic acid. In some embodiments, the method may include labeling a nucleic acid in a first cell. The method may include: (a) generating complementary DNA (cDNA) in a plurality of cells, including the first cell, by reverse transcribing RNA using a reverse transcription primer that includes a 5' overhang sequence; (b) dividing the plurality of cells into (n) aliquots; (c) providing a plurality of nucleic acid tags to each of the n aliquots, where each labeling sequence of the plurality of nucleic acid tags provided in a given aliquot is the same, and each of the n aliquots is provided with a different labeling sequence; (d) binding at least one of the cDNAs of each of the n aliquots to a nucleic acid tag; (e) combining the n aliquots; and (f) repeating steps (b), (c), (d), and (e) with the combined aliquot. In various embodiments, the plurality of cells may be selected from eukaryotic cells and prokaryotic cells. In various other embodiments, the plurality of cells can be selected from, but is not limited to, at least one of mammalian cells, yeast cells, and / or bacterial cells.

[0011] In certain embodiments, each nucleic acid tag may include a first strand including a 3' hybridization sequence extending from the 3' end of the labeled sequence and a 5' hybridization sequence extending from the 5' end of the labeled sequence. Each nucleic acid tag may also include a second strand including an overhang sequence. The overhang sequence may include (i) a first portion complementary to at least one of the 5' hybridization sequence and the 5' overhang sequence, and (ii) a second portion complementary to the 3' hybridization sequence. In some embodiments, the nucleic acid tag (e.g., the final nucleic acid tag) may include a capture agent, such as, but not limited to, 5' biotin. cDNA labeled with a nucleic acid tag including 5'-biotin may enable or permit attachment or coupling of streptavidin-coated magnetic beads to the cDNA. In some other embodiments, a plurality of beads may be coated with a capture strand (i.e., a nucleic acid sequence) configured to hybridize to the final sequence overhang of the barcode. In yet some other embodiments, the cDNA may be purified or isolated by using commercially available kits (e.g., RNEASY™ kits). .

[0012] In various embodiments, step (f) (i.e., steps (b), (c), (d), and (e)) can be repeated a sufficient number of times to generate a unique set of labeled sequences for the cDNA of the first cell. Stated another way, step (f) may be repeated a number of times such that the cDNA of the first cell has a first unique set of labeled sequences, the cDNA of the second cell has a second unique set of labeled sequences, the cDNA of the third cell has a third unique set of labeled sequences, and so on. The disclosed method provides for labeling of cDNA sequences derived from a single cell with a unique barcode, which can aid in the identification or identification of the cell from which the cDNA was derived. In other words, a portion, a majority, or substantially all of the cDNA derived from a single cell may have the same barcode, which may not be repeated in the cDNA derived from one or more other cells in the sample (e.g., from a second cell, a third cell, a fourth cell, etc.).

[0013] In some embodiments, the barcoded cDNAs can be mixed together and sequenced (e.g., using NGS) to collect data on RNA expression at the single cell level. For example, certain embodiments of the disclosed methods can be useful for evaluating, analyzing, or studying the transcriptome (the various RNA species transcribed from the genome of a given cell) of one or more individual cells.

[0014] As discussed above, aliquots or groups of cells may be separated into different reaction vessels or containers, and a first set of nucleic acid tags may be added to multiple cDNA transcripts. A vessel or container may also be referred to as a container, sample, and well herein. Thus, the terms vessel, container, container, sample, and well may be used interchangeably herein. The aliquots of cells may then be regrouped, mixed, and separated again, and a second set of nucleic acid tags may be added to the first set of nucleic acid tags. In various embodiments, the same nucleic acid tag may be added to more than one aliquot of cells in a single or given round of labeling. However, after repeated rounds of separation, tagging, and repooling, the cDNA of each cell may be bound to a unique combination or sequence of nucleic acid tags that form a barcode. In some embodiments, cells in a single sample may be separated into a number of different reaction vessels. For example, the number of reaction vessels may include four 1.5 ml microcentrifuge tubes, multiple wells of a 96-well plate, or another suitable number and type of reaction vessels.

[0015] In certain embodiments, step (f) (i.e., steps (b), (c), (d), and (e)) may be repeated a number of times selected from 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, etc. In certain other embodiments, step (f) may be repeated a sufficient number of times to provide a high probability that each cell's cDNA is bound to a unique barcode. This number may be selected to provide a greater than 50% probability, a greater than 90% probability, a greater than 95% probability, a greater than 99% probability, or some other probability that each cell's cDNA is bound to a unique barcode. In yet other embodiments, step (f) may be repeated some other suitable number of times.

[0016] In some embodiments, the method of labeling nucleic acid of a first cell may include a step of fixing the plurality of cells prior to step (a). For example, components of the cells may be fixed or crosslinked such that the components are immobilized or held in place. The plurality of cells may be fixed using formaldehyde in phosphate buffered saline (PBS). For example, the plurality of cells may be fixed using formaldehyde in phosphate buffered saline (PBS). The cells may be fixed with about 1-4% formaldehyde in PBS. In various embodiments, the plurality of cells may be fixed using methanol (e.g., 100% methanol) at about -20°C or about 25°C. In various other embodiments, the plurality of cells may be fixed using methanol (e.g., 100% methanol) at about -20°C to about 25°C. In still various other embodiments, the plurality of cells may be fixed using ethanol (e.g., about 70-100% ethanol) at about -20°C or room temperature. In still various other embodiments, the plurality of cells may be fixed using ethanol (e.g., about 70-100% ethanol) at about -20°C to room temperature. In still various other embodiments, the plurality of cells may be fixed using acetic acid, for example, at about -20°C. In still various other embodiments, the plurality of cells may be fixed using acetone, for example, at about -20°C. Other suitable methods for fixing the plurality of cells are within the scope of this disclosure.

[0017] In certain embodiments, the method of labeling the nucleic acid of the first cell may include permeabilizing the plurality of cells prior to step (a). For example, holes or openings may be formed in the outer membrane of the plurality of cells. TRITON™ X-100 may be added to the plurality of cells, followed by optional addition of HCl to form one or more holes. About 0.2% TRITON™ X-100 may be added to the plurality of cells, followed by addition of, for example, about 0.1 N HCl. In certain other embodiments, the plurality of cells may be permeabilized using ethanol (e.g., about 70% ethanol), methanol (e.g., about 100% methanol), Tween 20 (e.g., about 0.2% Tween 20), and / or NP-40 (e.g., about 0.1% NP-40). In various embodiments, the method of labeling the nucleic acid of the first cell may include fixing and permeabilizing the plurality of cells prior to step (a).

[0018] In some embodiments, the cells may be adherent cells (e.g., adherent mammalian cells). Fixation, permeabilization, and / or reverse transcription may be performed or carried out on the adherent cells (e.g., on cells adhered to a plate). For example, the adherent cells may be fixed, permeabilized, and / or reverse transcribed, followed by trypsinization to detach the cells from the surface. Alternatively, the adherent cells may be detached prior to the separation and / or tagging steps. In some other embodiments, the adherent cells may be trypsinized prior to the fixation and / or permeabilization steps.

[0019] In some embodiments, the method of labeling the nucleic acid of the first cell may include ligating at least two of the nucleic acid tags attached to the cDNA. The ligation may be performed before or after the lysis and / or cDNA purification step. The ligation may include covalently linking the 5' phosphate sequence of the nucleic acid tag to the adjacent strand or 3' end of the nucleic acid tag such that the individual tags are formed into a continuous or substantially continuous barcode sequence attached to the 3' end of the cDNA sequence. In various embodiments, a double-stranded DNA or RNA ligase may be used with an additional linker strand configured to hold the nucleic acid tag with the adjacent nucleic acid in a "nicked" double-stranded conformation. The double-stranded DNA or RNA ligase may then be used to fill the "nick". In various other embodiments, a single-stranded DNA or RNA ligase may be used without an additional linker. In certain embodiments, the ligation may be performed in multiple cells.

[0020] In Fig. 1, multiple nucleic acid tags are ligated to form a substantially continuous label or barcode.For example, when multiple nucleic acid tags are added, each cDNA transcript can be bound or linked to a series of nucleic acid tags.By using ligase, the parts of nucleic acid tags can be ligated or covalently linked to form a substantially continuous label or barcode that is bound or attached to cDNA transcript.

[0021] In certain other embodiments, the method may include, for example after step (f), lysing the plurality of cells (i.e., breaking down cellular structures) to release cDNA from within the plurality of cells. In some embodiments, the plurality of cells are lysed in a lysis solution (e.g., 10 mM Tris-HCl (pH 7.9), 50 mM EDTA (pH 7.9), 0.2M The cells may be lysed in NaCl, 2.2% SDS, 0.5 mg / ml ANTI-RNase (protein ribonuclease inhibitor; AMBION®), and 1000 mg / ml Proteinase K (AMBION®), for example, at about 55° C. for about 1-3 hours with shaking (e.g., vigorous shaking). In some other embodiments, the cells may be lysed by sonication and / or by at least one passage through an 18-25 gauge needle. In yet some other embodiments, the cells may be lysed by heating to about 70-90° C. For example, the cells may be lysed by heating to about 70-90° C. for about 1 hour or longer. cDNA may then be isolated from the lysed cells. In some embodiments, RNase H may be added to the cDNA to remove RNA. The method may further include ligating at least two of the nucleic acid tags attached to the released cDNA. In some other embodiments, the method of labeling nucleic acid of a first cell may include ligating at least 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 30, 35, 40, 45, 50, etc., of the nucleic acid tags attached to the cDNA.

[0022] In various embodiments, the method of labeling the nucleic acid of a first cell may include removing one or more unbound nucleic acid tags (e.g., washing the plurality of cells). For example, the method may include removing a portion, a majority, or substantially all of the unbound nucleic acid tags. The unbound nucleic acid tags may be removed such that further rounds of the disclosed method are not contaminated with one or more unbound nucleic acid tags from a previous round of a given method. In some embodiments, the unbound nucleic acid tags may be removed by centrifugation. For example, the plurality of cells may be centrifuged such that a pellet of cells is formed at the bottom of a centrifuge tube. The supernatant (i.e., the liquid containing the unbound nucleic acid tags) may be removed from the centrifuged cells. The cells may then be resuspended in a buffer solution (e.g., a new buffer solution free or substantially free of unbound nucleic acid tags). In another example, the plurality of cells may be coupled or linked to magnetic beads that are coated with an antibody configured to bind to a cell membrane or nuclear membrane. A magnet may then be used to attract the plurality of cells to one side of a reaction vessel, causing the plurality of cells to pellet. In some other embodiments, the plurality of cells may be placed in a cell strainer (e.g., a PLURISTRAINER® cell strainer) and washed with a wash buffer. For example, the plurality of cells may be left in the cell strainer and the wash buffer is passed through the cell strainer. The wash buffer may include a surfactant, detergent, and / or about 5-60% formamide.

[0023] As discussed above, multiple cells may be repooled, and the method may be repeated any number of times to add more tags to the cDNA, generating a set of nucleic acid tags that can act as barcodes. With more rounds added, the number of paths the cells can take increases, and as a result, the number of possible barcodes that can be generated also increases. With enough rounds and splits, the number of possible barcodes will be much greater than the number of cells, increasing the likelihood that each cell will have a unique barcode. For example, if splitting is done in a 96-well plate, after four splits, 96 4 =84,934,656 possible barcodes.

[0024] In some embodiments, the reverse transcription primer may be configured to reverse transcribe all or substantially all of the RNA in a cell (e.g., random hexamers with 5' overhangs). In some other embodiments, the reverse transcription primer may be configured to reverse transcribe RNA with a poly(A) tail (e.g., a poly(dT) primer, such as a dT(15) primer with a 5' overhang). In yet some other embodiments, the reverse transcription primer may be configured to reverse transcribe a given RNA (e.g., a transcript-specific primer). For example, the reverse transcription primer may be configured to barcode a particular transcript such that fewer transcripts per cell can be profiled, but each of the transcripts can be profiled across a larger number of cells.

[0025] FIG. 2 shows the formation of cDNA by in situ reverse transcription. Panel A shows cells that have been fixed and permeabilized. Panel B shows the addition of poly(T) primers that can template the reverse transcription of polyadenylated transcripts, as discussed above. Panel C shows the addition of random hexamers that can template the reverse transcription of virtually any transcript, as discussed above. Panel D shows the addition of primers that are designed to target specific transcripts so that only a subset of transcripts can be amplified, as discussed above. Panel E shows the cells of panel A after reverse transcription, where cDNA is shown hybridized to RNA.

[0026] Reverse transcription may be performed or carried out on a plurality of cells. In certain embodiments, reverse transcription may be performed on a plurality of fixed and / or permeabilized cells. In some embodiments, a mutant of M-MuLV reverse transcriptase may be used for reverse transcription. Any suitable method of reverse transcription is within the scope of this disclosure. For example, the reverse transcription mix may include a reverse transcription primer that includes a 5' overhang, and the reverse transcription primer may be configured to initiate reverse transcription and / or act as a binding sequence for the nucleic acid tag. In some other embodiments, the portion of the reverse transcription primer that is configured to bind to the RNA and / or initiate reverse transcription may be one of the following: random hexamers, septamers, octamers, nonamers, decamers, poly(T) extensions, and and / or one or more gene-specific primers.

[0027] Another aspect of the disclosure relates to a method for uniquely labeling a molecule within a cell or a plurality of cells, in some embodiments, the method includes the steps of: (a) attaching an adaptor sequence or universal adaptor to a molecule within a plurality of cells; (b) dividing the plurality of cells into at least two primary aliquots, the at least two primary aliquots including at least a first primary aliquot and a second primary aliquot; (c) providing primary nucleic acid tags to the at least two primary aliquots, the primary nucleic acid tag provided in the first primary aliquot being different from the primary nucleic acid tag provided in the second primary aliquot; (d) attaching an adaptor sequence within each of the at least two primary aliquots to the provided primary nucleic acid tag. (e) combining the at least two primary aliquots; (f) dividing the combined primary aliquot into at least two secondary aliquots, the at least two secondary aliquots including at least a first secondary aliquot and a second secondary aliquot; (g) providing secondary nucleic acid tags to the at least two secondary aliquots, the secondary nucleic acid tag provided in the first secondary aliquot being different from the secondary nucleic acid tag provided in the second secondary aliquot; and (h) binding molecules within each of the at least two secondary aliquots to the provided secondary nucleic acid tags.

[0028] In certain embodiments, the method further comprises the step (i), i.e., using the subsequent aliquot. The method may further include repeating steps (e), (f), (g), and (h) for each of the molecules in the single cell. Step (i) may be repeated a sufficient number of times to generate a set of unique nucleic acid tags for the molecules in the single cell. In various embodiments, the number of times may be selected from 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, 100, etc. In certain other embodiments, step (i) may be repeated another suitable number of times.

[0029] In some embodiments, the molecule may be located within a cell or cells. In some other embodiments, the molecule may be coupled to a cell or cells. For example, the molecule may be a cell surface molecule. In yet some other embodiments, the molecule may be located within a cell or cells and / or coupled to a cell or cells.

[0030] As discussed above, the method may include fixing and / or permeabilizing the plurality of cells prior to step (a). In various embodiments, each of the nucleic acid tags may include a first strand. The first strand may include a barcode sequence including a 3' end and a 5' end. The first strand may further include a 3' hybridization sequence and a 5' hybridization sequence adjacent to the 3' end and the 5' end of the barcode sequence, respectively. In some embodiments, each of the nucleic acid tags may include a second strand. The second strand may include a first portion complementary to at least one of the 5' hybridization sequence and the adapter sequence, and a second portion complementary to the 3' hybridization sequence.

[0031] In certain embodiments, the molecule is a macromolecule. In various embodiments, the molecule is selected from at least one of RNA, cDNA, DNA, a protein, a peptide, and / or an antigen.

[0032] In some embodiments, the molecule is RNA and the adapter sequence may be single-stranded. Further, step (a) may include one of ligating the 5' end of the single-stranded adapter sequence to the 3' end of the RNA and / or ligating the 3' end of the single-stranded adapter sequence to the 5' end of the RNA. In some other embodiments, the molecule is RNA and step (a) may include hybridizing the adapter sequence to the RNA.

[0033] Methods relating to attaching or coupling an adaptor sequence to an RNA can be used, for example, for RNA transcriptome sequencing, ribosome profiling, small RNA sequencing, non-coding RNA sequencing, and / or RNA structural profiling. In some embodiments, the plurality of cells may be fixed and / or permeabilized. The 5' end of a single-stranded adaptor sequence may be ligated to the 3' end of the RNA (see Figures 3A and 3B). In certain embodiments, the ligation may be performed or carried out by T4 RNA ligase 1. In certain other embodiments, the ligation may be carried out by T4 RNA ligase 1 using a single-stranded adaptor sequence that includes a 5' phosphate. In various embodiments, the ligation may be carried out by THERMOSTABLE 5'APPDNA / RNA LIGASE™ (NEW ENGLAND BIOLABS®). In various other embodiments, the ligation may be carried out by THERMOSTABLE 5'APPDNA / RNA LIGASE™ using a 5'preadenylated single-stranded adaptor sequence. Other suitable ligases and adapter sequences are within the scope of this disclosure.

[0034] In some embodiments, the RNA can be labeled with an adapter sequence, for example using hybridization by Watson-Crick base pairing (see FIG. 4). The adapter sequence may be configured to initiate reverse transcription to form or generate cDNA after the labeling step and / or cell lysis, as discussed above (see FIG. 5).

[0035] The 3' end of the single stranded adaptor sequence may be ligated to the 5' end of the RNA. In certain embodiments, the ligation may be performed or carried out by T4 RNA ligase 1. In certain other embodiments, the ligation may be carried out by T4 RNA ligase 1 using an RNA that includes a 5' phosphate. In various embodiments, the ligation may be carried out by THERMOSTABLE 5'APP DNA / RNA LIGASE™ (NEW ENGLAND BIOLABS®). In various other embodiments, the ligation may be carried out by THERMOSTABLE 5'APP DNA / RNA LIGASE™ using a 5' pre-adenylated RNA. As noted above, other suitable ligases and adaptor sequences are within the scope of the present disclosure.

[0036] In some embodiments, the molecule may be a cDNA. Methods involving binding or coupling an adapter sequence to a cDNA may be used, for example, for RNA transcriptome sequencing. In certain embodiments, a plurality of cells may be fixed and / or permeabilized. Reverse transcription may be performed on a plurality of fixed and / or permeabilized cells using a primer that includes an adapter sequence at its 5' end. As discussed above, the 3' end of the primer may be a gene-specific random hexamer or a poly(T) sequence. The resulting cDNA may include an adapter sequence at its 5' end (see FIG. 5).

[0037] In some embodiments, when the molecule is DNA (e.g., genomic DNA), the method may further comprise digesting the DNA with a restriction enzyme prior to step (a), and step (a) may further comprise ligating an adapter sequence to the digested DNA.

[0038] Methods relating to binding or coupling adaptor sequences to DNA can be used, for example, for whole genome sequencing, targeted genome sequencing, DNase-Seq, ChIP sequencing, and / or ATAC-seq. In certain embodiments, DNA can be digested using one or more restriction enzymes to at least one of blunt-ended fragments and / or fragments with overhanging sequences. A single-stranded universal adaptor or a partial double-stranded sequence with an adaptor sequence, with one end overhanging, can be ligated to the digested genomic DNA. For example, a partial double-stranded sequence with a single-stranded adaptor sequence with an overhang, where the overhang is compatible with the overhang generated by one or more restriction enzymes, can be ligated to the digested genomic DNA.

[0039] In various embodiments, the adapter sequence can be integrated (e.g., directly integrated) into genomic DNA using Tn5 transposase, and the transposase can be released by adding sodium dodecyl sulfate (SDS) to expose the adapter sequence. Other transposases and methods for integrating adapter sequences into genomic DNA are within the scope of this disclosure.

[0040] In certain embodiments, the molecule is a protein, peptide, and / or antigen, and the adapter sequence is a unique identifier sequence (e.g., comprising a nucleic acid) that is coupled to an antibody. The unique identifier sequence may be configured to uniquely identify an antibody to which the unique identifier sequence is bound. Further, step (a) may include binding an antibody comprising each of the adapter sequence and the unique identifier sequence to a protein, peptide, and / or antigen. In certain other embodiments, the molecule is a protein, peptide, and / or antigen, and the adapter sequence may be incorporated into an aptamer. Further, step (a) may include binding an aptamer to a protein, peptide, and / or antigen.

[0041] Methods relating to binding or coupling adapter sequences to proteins, peptides, and / or antigens can be used, for example, for protein quantification, peptide quantification, and / or antigen quantification. In various embodiments, adapter sequences can be attached (e.g., chemically attached) to antibodies. For example, adapter sequences can be attached to antibodies using chemistry known to those skilled in the art for mediating DNA-protein binding. Antibodies against different proteins can be labeled with a nucleic acid sequence or nucleic acid strand that includes a unique identifier sequence in addition to the adapter sequence. An antibody or set of antibodies can then be used in an immunostaining experiment to label a protein or set of proteins in fixed and / or permeabilized cells or tissues (see FIG. 6). The cells can then be subjected to the labeling or barcoding procedures disclosed herein.

[0042] In some embodiments, the nucleic acid sequence (e.g., DNA molecule) attached or bound to the antibody can be released from the antibody and / or adapter sequence. A sequencing reaction can reveal the unique identifier sequence associated with a given protein as well as the label or barcode associated with one or more unique cells. In certain embodiments, such methods can reveal or identify the number and / or type of proteins present in one or more cells.

[0043] In various embodiments, DNA and / or RNA aptamers can be used instead of or in addition to nucleic acid-modified (or DNA-modified) antibodies as described above (see FIG. 7). Adapter sequences (and target protein-specific antibodies) can be incorporated (e.g., directly) into the sequence of a given aptamer.

[0044] Another aspect of the present disclosure relates to a method for barcoding nucleic acids in cells. In some embodiments, the method may include the steps of: (a) generating cDNA in a plurality of cells by reverse transcribing RNA using a reverse transcription primer comprising a 5' overhang sequence; (b) dividing the plurality of cells into at least two aliquots; (c) providing a plurality of nucleic acid tags to each of the at least two aliquots, wherein each barcode sequence of the plurality of nucleic acid tags introduced in a given aliquot is the same, and a different barcode sequence is introduced in each aliquot; (d) binding at least one of the cDNAs of each of the at least two aliquots to the nucleic acid tag; (e) combining the at least two aliquots; and (f) repeating steps (b), (c), (d), and (e) at least once with the combined aliquot.

[0045] In certain embodiments, each nucleic acid tag may include a first strand that includes a 3' hybridization sequence extending from the 3' end of the barcode sequence and a 5' hybridization sequence extending from the 5' end of the barcode sequence. Each nucleic acid tag may also include a second strand that includes an overhanging sequence, the overhanging sequence including (i) a first portion that is complementary to at least one of the 5' hybridization sequence and the 5' overhanging sequence, and (ii) a second portion that is complementary to the 3' hybridization sequence.

[0046] FIG. 8 illustrates cell splitting, tagging, and pooling according to an embodiment of the present disclosure. Reverse transcribed cells may be split into reaction vessels or wells. Four wells are shown in FIG. 8. However, as discussed above, any suitable number of reaction vessels or wells may be used. One cell is highlighted to show its path through the process. As shown, the highlighted cell first enters well "a" where a first tag is added that hybridizes to the overhang of all cDNA transcripts (shown in a box). This tag carries a unique barcode region "a" that identifies the well from which the cell originated. After hybridization, the cells are all washed to remove excess tags, regrouped, and then redistributed into an equal number of wells. The highlighted cell then enters well "c" where a second tag is added to identify the well from which the cell originated. After the second round, the cells are split into four wells, one well, one well, and one well. 2 = 16 possible pathways for a cell to follow. This process can be repeated to add more tags to the cDNA transcripts, increasing the number of possible pathways that a cell can follow. Figures 9A and 9B show two exemplary workflows according to embodiments of the present disclosure.

[0047] Another aspect of the present disclosure relates to a kit for labeling a nucleic acid in at least a first cell. In some embodiments, the kit may include at least one reverse transcription primer comprising a 5' overhang sequence. The kit may also include a plurality of first nucleic acid tags. Each first nucleic acid tag may include a first strand. The first strand may include a 3' hybridization sequence extending from the 3' end of the first labeling sequence and a 5' hybridization sequence extending from the 5' end of the first labeling sequence. Each first nucleic acid tag may further include a second strand. The second strand may include an overhang sequence, and the overhang sequence may include (i) a first portion complementary to at least one of the 5' hybridization sequence and the 5' overhang sequence of the reverse transcription primer, and (ii) a second portion complementary to the 3' hybridization sequence.

[0048] The kit may further comprise a plurality of second nucleic acid tags. Each second nucleic acid tag may comprise a first strand. The first strand may comprise a 3' hybridization sequence extending from the 3' end of the second labeling sequence and a 5' hybridization sequence extending from the 5' end of the second labeling sequence. Each second nucleic acid tag may further comprise a second strand. The second strand may comprise an overhanging sequence, and the overhanging sequence may comprise (i) a first portion complementary to at least one of the 5' hybridization sequence and the 5' overhanging sequence of the reverse transcription primer, and (ii) a second portion complementary to the 3' hybridization sequence. In some embodiments, the first labeling sequence may be different from the second labeling sequence.

[0049] In some embodiments, the kit may also include one or more additional nucleic acid tags. Each nucleic acid tag of the one or more additional nucleic acid tags may include a first strand. The first strand may include a 3' hybridization sequence extending from the 3' end of the labeling sequence and a 5' hybridization sequence extending from the 5' end of the labeling sequence. Each nucleic acid tag of the one or more additional nucleic acid tags may also include a second strand. The second strand may include an overhanging sequence, the overhanging sequence including (i) a first portion complementary to at least one of the 5' hybridization sequence and the 5' overhanging sequence of the reverse transcription primer, and (ii) a second portion complementary to the 3' hybridization sequence. In some embodiments, the labeling sequence may be different for each given additional nucleic acid tag of the plurality.

[0050] In various embodiments, the kit may further comprise at least one of a reverse transcriptase, a fixation agent, a permeabilization agent, a ligation agent, and / or a lysis agent.

[0051] Another aspect of the present disclosure relates to a kit for labeling the molecule in at least a first cell.For example, the kit disclosed above can be adapted to label one or more of RNA, cDNA, DNA, protein, peptide or antigen in at least a first cell.

[0052] Another aspect of the present disclosure relates to a method for uniquely labeling RNA molecules in a plurality of cells.The method includes the steps of: (a) fixing and permeabilizing the first plurality of cells prior to step (b), wherein the first plurality of cells is fixed and permeabilized at less than about 8° C.; (b) reverse transcribing RNA molecules in the first plurality of cells to form complementary DNA (cDNA) molecules in the first plurality of cells, wherein reverse transcribing the RNA molecules comprises coupling a primer to the RNA molecules, wherein the primer comprises at least one of a poly(T) sequence or a random sequence; (c) dividing the first plurality of cells comprising the cDNA molecules into at least two primary aliquots, wherein the at least two primary aliquots comprise a first primary aliquot and a second primary aliquot; (d) providing primary nucleic acid tags to the at least two primary aliquots, wherein the primary nucleic acid tag provided in the first primary aliquot is different from the primary nucleic acid tag provided in the second primary aliquot; (e) coupling the cDNA molecules in each of the at least two primary aliquots to the provided primary nucleic acid tags; (f) coupling the cDNA molecules in each of the at least two primary aliquots to the provided primary nucleic acid tags. (g) dividing the combined primary aliquots into at least two secondary aliquots, the at least two secondary aliquots comprising a first secondary aliquot and a second secondary aliquot; (h) providing secondary nucleic acid tags to the at least two secondary aliquots, the secondary nucleic acid tag provided in the first secondary aliquot being different from the secondary nucleic acid tag provided in the second secondary aliquot; (i) coupling cDNA molecules in each of the at least two secondary aliquots to the provided secondary nucleic acid tags; (j) repeating steps (f), (g), (h), and (i) with subsequent aliquots, the final nucleic acid tag comprising a capture agent; (k) combining the final aliquots; (l) lysing the first plurality of cells to release cDNA molecules from within the first plurality of cells to form a lysate; and / or (m) adding a protease inhibitor and / or a binding agent to the lysate, such that the cDNA molecules bind to the binding agent.

[0053] The method includes further dividing the combined final aliquot into at least two final aliquots, the at least two final aliquots including a first final aliquot and a second final aliquot. In some embodiments, the first plurality of cells may be fixed and permeabilized at less than about 8°C, less than about 7°C, less than about 6°C, less than about 5°C, less than about 4°C, less than about 4°C, less than about 3°C, less than about 2°C, less than about 1°C, or another suitable temperature. In certain embodiments, the method may include distributing the cells. For example, after the last or final round of barcoding (by ligation), the cells may be pooled and then lysed, and the cells may then be distributed into different lysate aliquots. Each lysate aliquot may contain a predetermined number of cells.

[0054] For example, with reference to step (m), the protease inhibitor may include phenylmethanesulfonyl fluoride (PMSF), 4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride (AEBSF), combinations thereof, and / or another suitable protease inhibitor. For example, with reference to steps (j), (k), (l), and / or (m), the capture agent may include biotin or another suitable capture agent. Further, the binding agent may include avidin (e.g., streptavidin) or another suitable binding agent.

[0055] In certain embodiments, the method for uniquely labeling RNA molecules in a plurality of cells (e.g., after step (m)) may further include (n) using a template switch oligonucleotide to perform template switching of the cDNA molecules bound to the binder; (o) amplifying the cDNA molecules to form an amplified cDNA molecule solution; and / or (p) introducing a solid phase reversible immobilization (SPRI) bead solution into the amplified cDNA molecule solution to remove polynucleotides less than about 200 base pairs, less than about 175 base pairs, or less than about 150 base pairs (see DeAngelis, MM et al., Nucleic Acids Research (1995) 23(22):4742). In other words, the cDNA molecules may be bound to streptavidin beads in the lysate. Template switching of the cDNA molecules attached to the beads may be performed (e.g., to add an adaptor to the 3' end of the cDNA molecules). PCR amplification of the cDNA molecules can then be performed, followed by the addition of SPRI beads to remove polynucleotides less than about 200 base pairs. The ratio of SPRI bead solution to amplified cDNA molecule solution can be about 0.9:1 to about 0.7:1, about 0.875:1 to about 0.775:1, about 0.85:1 to about 0.75:1, about 0.825:1 to about 0.725:1, about 0.8:1, or another suitable ratio. Additionally, the SPRI bead solution can contain about 1M to 4M NaCl, about 2M to 3M NaCl, about 2.25M to 2.75M NaCl, about 2.5M NaCl, or another suitable amount of NaCl. The SPRI bead solution may also contain about 15% to 25% w / v polyethylene glycol (PEG), with the molecular weight of PEG being about 7,000 g / mol to 9,000 g / mol (PEG 8000). In various embodiments, the SPRI bead solution may contain about 17% to 23% w / v PEG 8000, about 18% to 22% w / v PEG 8000, about 19% to 21% w / v PEG 8000, about 20% w / v PEG 8000, or another suitable w / v % PEG 8000.

[0056] The method for uniquely labeling RNA molecules in a plurality of cells may further comprise the step of adding a common adaptor sequence to the 3' end of the released cDNA molecules. The common adaptor sequence may be the same or substantially the same adaptor sequence for each of the cDNA molecules (i.e., within a given experiment). The addition of the common adaptor may be performed or carried out in a solution containing up to about 10 w / v % PEG, the molecular weight of which is about 7,000 g / mol to 9,000 g / mol. In certain embodiments, the common adaptor sequence may be added to the 3' end of the released cDNA molecules by performing template switching (see Picelli, S. et al., Nature Methods vol. 10, pp. 1096-1098 (2013)).

[0057] Step (j) can be repeated a sufficient number of times to generate a set of unique nucleic acid tags for the nucleic acids of a single cell, for example, the number of times can be selected from 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, and 100.

[0058] In various embodiments, the primers in step (b) may further comprise a first specific barcode. Stated differently, the first barcode added to the cDNA molecules in a particular vessel, mixture, reaction, container, sample, well, or bath may be predetermined (e.g., specific to a given vessel, mixture, reaction, container, sample, well, or bath). For example, 48 sets of different well-specific RT primers may be used (e.g., in a 48-well plate). Thus, if there are 48 samples (e.g., cells, tissues, etc.), each sample can obtain a unique well-specific barcode. However, if there are only four samples, each sample may have 12 different sets of well-specific RT primers. The user can determine which 12 correspond to each sample. can be known, and thus the user can retrieve the sample identity. Other numbers of first specific barcodes (or well-specific RT primers) are also within the scope of the present disclosure. Such configurations can enable or provide for multiplexing of the method as further described in Example 16.

[0059] The method includes the steps of: (q) reverse transcribing the RNA molecules in a second plurality of cells to form cDNA molecules in the second plurality of cells, where reverse transcribing the RNA molecules comprises coupling a specific primer to the RNA molecules, where the primer comprises a second specific barcode and at least one of a poly(T) sequence or a random sequence, where the first specific barcode is different from the second specific barcode such that the cDNA molecules derived from the first plurality of cells can be distinguished relative to the cDNA molecules derived from the second plurality of cells; (r) dividing the second plurality of cells comprising the cDNA molecules into at least two primary aliquots, where the at least two primary aliquots comprise a first primary aliquot and a second primary aliquot; (s) providing primary nucleic acid tags to the at least two primary aliquots, where the primary nucleic acid tags provided in the first primary aliquot are different from the primary nucleic acid tags provided in the second primary aliquot. The method may further include the steps of: (t) coupling the cDNA molecules in each of the at least two primary aliquots with the provided primary nucleic acid tag; (u) combining the at least two primary aliquots; (v) dividing the combined primary aliquots into at least two secondary aliquots, the at least two secondary aliquots including a first secondary aliquot and a second secondary aliquot; (w) providing secondary nucleic acid tags to the at least two secondary aliquots, the secondary nucleic acid tag provided in the first secondary aliquot being different from the secondary nucleic acid tag provided in the second secondary aliquot; (x) coupling the cDNA molecules in each of the at least two secondary aliquots with the provided secondary nucleic acid tag; and / or (y) repeating steps (u), (v), (w), and (x) with subsequent aliquots, the final nucleic acid tag comprising a capture agent. The above steps used with the first plurality of cells (e.g., steps (k), (l), and (m)) may also be adapted for use with the second plurality of cells.

[0060] In various embodiments, each of the nucleic acid tags may comprise a first strand, the first strand comprising (i) a barcode sequence comprising a 3' end and a 5' end, and (ii) a 3' hybridization sequence and a 5' hybridization sequence adjacent to the 5' end and the 3' end of the barcode sequence, respectively. Each of the nucleic acid tags may also comprise a second strand, the second strand comprising (i) a first portion complementary to at least one of the 5' hybridization sequence and the adapter sequence, and (ii) a second portion complementary to the 3' hybridization sequence.

[0061] The method of uniquely labeling RNA molecules in a plurality of cells may further comprise ligating at least two (or more) of the nucleic acid tags attached to the cDNA molecules. Ligation may be performed in the first plurality of cells.

[0062] The method may further comprise removing unbound nucleic acid tags. In some embodiments, the method may comprise ligating at least two of the nucleic acid tags bound to the released cDNA molecules. A majority of the cDNA molecules bound to nucleic acid tags from a single cell may comprise the same set of bound nucleic acid tags. In various embodiments, the plurality of cells (e.g., the first and second plurality of cells) may be selected from at least one of mammalian cells, yeast cells, and bacterial cells.

[0063] Another aspect of the present disclosure is directed to a method of labeling a nucleic acid in a first cell. In embodiments, the method includes the steps of: (a) generating cDNA molecules in a plurality of cells, including the first cell, by reverse transcribing RNA using (i) a first reverse transcription primer comprising a 5'-overhanging sequence, where the first reverse transcription primer is configured to reverse transcribe RNA having a poly(A) tail, and / or (ii) at least one second reverse transcription primer comprising a 5'-overhanging sequence and at least one of a random hexamer, a random heptamer, a random octamer, a random nonamer, and a random decamer; (b) dividing the plurality of cells into (n) aliquots; and (c) isolating a plurality of nucleic acid tags into each of the n aliquots. (d) binding at least one of the cDNA molecules of each of the n aliquots to a nucleic acid tag; (e) combining the n aliquots; (f) repeating steps (b), (c), (d), and (e) with the combined aliquots; (g) combining the final aliquots; (h) lysing the plurality of cells, including the first cell, to release the cDNA molecules from within the plurality of cells, including the first cell, to form a lysate; and / or (i) adding a protease inhibitor and / or a binding agent to the lysate, whereby the cDNA molecules bind to the binding agent.

[0064] For example, referring to step (c), each nucleic acid tag may include a first strand including (i) a 3' hybridization sequence extending from the 3' end of the labeled sequence, and (ii) a 5' hybridization sequence extending from the 5' end of the labeled sequence. Each nucleic acid tag may also include a second strand including an overhanging sequence, the overhanging sequence including (i) a first portion complementary to at least one of the 5' hybridization sequence and the 5' overhanging sequence, and (ii) a second portion complementary to the 3' hybridization sequence. In some embodiments, the labeling sequences of the multiple nucleic acid tags provided in a given aliquot may be the same, and each of the n aliquots may be provided with a different labeling sequence.

[0065] In certain embodiments, step (f) can be repeated a sufficient number of times to generate a set of unique labeled sequences for the cDNA molecules of the first cell. For example, the number of times can be selected from 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 30, 35, 40, 45, 50, 55, 60, 65, 70, 75, 80, 85, 90, 95, and 100.

[0066] The cDNA molecules may be formed or generated in an aliquot (e.g., a reaction mixture). The concentration of the first reverse transcription primer in the aliquot may be about 0.5 μM to about 10 μM, about 1 μM to about 7 μM, about 1.5 μM to about 4 μM, about 2 μM to about 3 μM, about 2.5 μM, or another suitable concentration. The concentration of the second reverse transcription primer in the aliquot may be about 0.5 μM to about 10 μM, about 1 μM to about 7 μM, about 1.5 μM to about 4 μM, about 2 μM to about 3 μM, about 2.5 μM, or another suitable concentration.

[0067] In some embodiments, the method may include a step of fixing the plurality of cells prior to step (a). The plurality of cells may be fixed at less than about 8° C., less than about 7° C., less than about 6° C., less than about 5° C., less than about 4° C., less than about 4° C., less than about 3° C., less than about 2° C., less than about 1° C., or another suitable temperature. In certain embodiments, the method may include a step of permeabilizing the plurality of cells prior to step (a). The plurality of cells may be permeabilized at less than about 8° C., less than about 7° C., less than about 6° C., less than about 5° C., less than about 4° C., less than about 4° C., less than about 3° C., less than about 2° C., less than about 1° C., or another suitable temperature.

[0068] The method of labeling nucleic acid in a first cell may also include ligating at least two of the nucleic acid tags attached to the cDNA molecule. In various embodiments, the ligation may be performed in a plurality of cells. The method may further include removing unbound nucleic acid tags. Further, at least one of the first and second reverse transcription primers may reverse transcribe the predetermined RNA or may be configured to reverse transcribe the predetermined RNA.

[0069] In various embodiments, the final nucleic acid tag may include a capture agent. Additionally, the method may include, after step (f), lysing the plurality of cells to release cDNA molecules from within the plurality of cells to form a lysate. The method may also include adding a protease inhibitor and / or a binding agent to the lysate to isolate the cDNA molecules. As discussed above, the protease inhibitor may include PMSF, AEBSF, combinations thereof, and / or another suitable protease inhibitor. The capture agent may include biotin or another suitable capture agent, and the binding agent may include avidin (e.g., streptavidin) or another suitable binding agent.

[0070] The method of labeling a nucleic acid in a first cell may also include (j) performing template switching of the cDNA molecules bound to a binder; (k) amplifying the cDNA molecules to form an amplified cDNA molecule solution; and (l) introducing a SPRI bead solution into the amplified cDNA molecule solution to remove polynucleotides less than about 200 base pairs, less than about 175 base pairs, or less than about 150 base pairs. The ratio of SPRI bead solution to amplified cDNA molecule solution may be about 0.9:1 to about 0.7:1, about 0.875:1 to about 0.775:1, about 0.85:1 to about 0.75:1, about 0.825:1 to about 0.725:1, about 0.8:1, or another suitable ratio.

[0071] Additionally, the SPRI bead solution may contain about 1M-4M NaCl, about 2M-3M NaCl, about 2.25M-2.75M NaCl, about 2.5M NaCl, or another suitable amount of NaCl. The SPRI bead solution may also contain about 15%-25% w / v PEG, with the molecular weight of the PEG being about 7,000g / mol-9,000g / mol. In various embodiments, the SPRI bead solution may contain about 17%-23% w / v PEG 8000, about 18%-22% w / v PEG 8000, about 19%-21% w / v PEG 8000, about 20% w / v PEG 8000, or another suitable w / v % PEG 8000.

[0072] The method for uniquely labeling RNA molecules in a plurality of cells may further include the step of adding a common adaptor sequence to the 3' end of the released cDNA molecules. As discussed above, the addition of the common adaptor may be performed or carried out in a solution containing up to about 10 w / v % PEG, the molecular weight of PEG being about 7,000 g / mol to 9,000 g / mol. In certain embodiments, the common adaptor sequence may be added to the 3' end of the released cDNA molecules by performing template switching.

[0073] In certain embodiments, any of the above-described methods may be adapted to label the nucleic acid molecules in nucleus or multiple nuclei.For example, the method may include the step of uniquely labeling the RNA molecules in multiple nuclei, or the step of labeling the nucleic acid in the first nucleus.

[0074] Another aspect of the present disclosure is directed to a kit for labeling a nucleic acid in a first cell. The kit may include a first reverse transcription primer comprising a 5'-overhanging sequence and may be configured to reverse transcribe RNA having a poly(A) tail. The kit may also include a second reverse transcription primer comprising a 5'-overhanging sequence and at least one of a random hexamer, a random heptamer, a random octamer, a random nonamer, and / or a random decamer. good.

[0075] In some embodiments, the kit may include a plurality of first nucleic acid tags. As discussed above, each first nucleic acid tag may include a first strand including (i) a 3' hybridization sequence extending from the 3' end of the first labeled sequence, and (ii) a 5' hybridization sequence extending from the 5' end of the first labeled sequence. Also, each first nucleic acid tag may include a second strand including an overhang sequence. The overhang sequence may include (i) a first portion that is complementary to at least one of the 5' hybridization sequence and the 5' overhang sequence of the first and second reverse transcription primers, and (ii) a second portion that is complementary to the 3' hybridization sequence.

[0076] In certain embodiments, the kit may also include a plurality of second nucleic acid tags. Each second nucleic acid tag may include a first strand including (i) a 3' hybridization sequence extending from the 3' end of the second labeled sequence, and (ii) a 5' hybridization sequence extending from the 5' end of the second labeled sequence. Each second nucleic acid tag may also include a second strand including an overhang sequence. The overhang sequence may include (i) a first portion complementary to at least one of the 5' hybridization sequence and the 5' overhang sequence of the first and second reverse transcription primers, and (ii) a second portion complementary to the 3' hybridization sequence. Furthermore, the first labeled sequence may be different from the second labeled sequence.

[0077] In various embodiments, the kit may also include a plurality of final nucleic acid tags. Each final nucleic acid tag may include a first strand including (i) a 3' hybridization sequence extending from the 3' end of the final labeled sequence, and (ii) a 5' hybridization sequence extending from the 5' end of the final labeled sequence. Each final nucleic acid tag may also include a second strand including an overhang sequence. The overhang sequence may include (i) a first portion complementary to at least one of the 5' hybridization sequence and the 5' overhang sequence of the first and second reverse transcription primers, and (ii) a second portion complementary to the 3' hybridization sequence. Each final nucleic acid tag may also include a capture agent. Furthermore, the final labeled sequence may be different from the first and second labeled sequences (and / or any other labeled sequences). In some embodiments, the kit may also include at least one of a reverse transcriptase, a fixation agent, a permeabilization agent, a ligation agent, a lysis agent, a protease inhibitor, and / or other suitable components.

[0078] As will be understood by those skilled in the art, each embodiment disclosed herein may comprise, consist essentially of, or consist of its specifically recited elements, steps, ingredients, or components. As used herein, the transitional terms "comprise" or "comprises" mean including, but not limited to, even in large amounts, and allowing for the intervening elements, steps, ingredients, or components not specifically recited. The transitional phrase "consisting of" excludes any element, step, ingredient, or component not specified. The transitional phrase "consisting essentially of" limits the scope of the embodiment to the specified elements, steps, ingredients, or components and those that do not materially affect the embodiment.

[0079] Unless otherwise indicated, all numerical values ​​expressing properties such as quantities of ingredients, molecular weights, reaction conditions, and the like used in the specification and claims are to be understood as being modified in all instances by the term "about." Accordingly, unless otherwise indicated, the numerical parameters set forth in the specification and appended claims are approximations that may vary depending upon the desired properties sought to be obtained by the present disclosure. However, and not as an attempt to limit the application of the doctrine of equivalents to the scope of the claims, each numerical parameter should, at the very least, be construed in light of the number of reported significant digits and by applying ordinary rounding techniques. In order to provide further clarity, the term "about," when used in conjunction with a listed numerical value or range, shall have the meaning reasonably recognized by one of ordinary skill in the art, i.e., to mean slightly more or slightly less than the stated value or range, and within a range of ±20% of the stated value, ±19% of the stated value, ±18% of the stated value, ±17% of the stated value, ±16% of the stated value, ±15% of the stated value, ±14% of the stated value, ±13% of the stated value, ±12% of the stated value, ±11% of the stated value, ±10% of the stated value, ±9% of the stated value, ±8% of the stated value, ±7% of the stated value, ±6% of the stated value, ±5% of the stated value, ±4% of the stated value, ±3% of the stated value, ±2% of the stated value, or ±1% of the stated value.

[0080] Notwithstanding, the numerical ranges and parameters setting forth the broad scope of the disclosure are approximations and the numerical values ​​set forth in the specific examples are reported as precisely as possible, however, any numerical value inherently contains certain errors necessarily resulting from the standard deviation found in their respective testing measurements.

[0081] The terms "a," "an," "the," and similar referents used in the context of describing this disclosure (particularly in the context of the claims below) should be construed to encompass both the singular and the plural, unless otherwise specified herein or clearly contradicted by context. The recitation of ranges of values ​​herein is merely intended to serve as a shorthand method of individually referring to each individual value falling within the range. Unless otherwise specified herein, each individual value is incorporated herein as if it were individually set forth herein. All methods described herein can be performed in any suitable order, unless otherwise specified herein or otherwise clearly contradicted by context. The use of any and all examples or exemplary language (e.g., "etc.") provided herein is intended merely to better clarify the disclosure and does not limit the scope of the disclosure as otherwise claimed. No language in this specification should be construed as indicating any non-claimed element essential to the practice of the disclosure.

[0082] Grouping of alternative elements or embodiments of the disclosure disclosed herein should not be interpreted as limiting. Each group member is referred to and claimed individually or in any combination with other members of the group or other elements found herein. It is understood that one or more members of a group may be included in or deleted from a group for reasons of convenience and / or patentability. When any such inclusion or deletion occurs, the specification is considered to include the modified group, and thus fulfills the description requirements of all Markush groups used in the appended claims.

[0083] The definitions and explanations used in this disclosure are intended and intended to govern all future constructions unless expressly and unambiguously modified in the following examples or where application of the meaning would render any construction meaningless or essentially meaningless, definitions should be taken from dictionaries known to those skilled in the art, such as Webster's Dictionary, 3rd Edition, or Oxford Dictionary of Biochemistry and Molecular Biology (edited by Anthony Smith, Oxford University Press, Oxford 2004). EXAMPLES

[0084] The following examples are illustrative of the methods and compositions of the present disclosure. In light of this disclosure, those of skill in the art will be able to Those skilled in the art will recognize that variations of these and other examples of the disclosed methods and compositions are possible without undue experimentation.

[0085] Example 1 - Immobilization and reverse transcription NIH / 3T3 (mouse) and Hela-S3 (human) cells can be grown to confluence in two separate 10 cm cell culture plates. Cells can be rinsed twice with 10 ml of 1× phosphate buffered saline (PBS), 1 ml of 0.05% trypsin can be added to each plate, and plates can be incubated at 37° C. for 5 minutes. Cells can be detached by tilting each plate at a 45° angle while pipetting trypsin across the plate, and this can be continued until all or substantially all cells are detached. Each cell line can be transferred to its own 15 ml conical centrifuge tube (FALCON™). 2 ml Dulbecco's Modified Eagle's Medium (DMEM) with 10% fetal bovine serum (FBS) can be added to each tube. The cell count for each tube can be calculated (e.g., with a hemocytometer or flow cytometer). For example, 200 μl of sample can be transferred from each tube into individual 1.7 ml microcentrifuge tubes (EPPENDORF®) and 100 μl of sample can be run on an ACCURI™ flow cytometer to calculate cell concentration.

[0086] An equal number of cells from each tube can be combined using as many cells as possible into a new single 15 ml conical centrifuge tube (FALCON™). A 5 minute centrifugation at 500 x g can be performed in the 15 ml conical centrifuge tube (FALCON™). It may be useful to use a bucket centrifuge so that the cells pellet to the bottom of the tube and not the side of the tube. The liquid can be aspirated without disturbing the cell pellet and the cells can be resuspended in 500 μl of 4% formaldehyde. The cells can then be left at room temperature (i.e., 20-25°C) for 10 minutes. 1.5 ml of 0.5% TRITON™ X-100 can be added to the tube and mixed gently with a pipette. The tube can then be centrifuged at 500 x g for 5 minutes. Again, the liquid can be aspirated without disturbing the pellet and the pellet can be washed twice with 1 ml of PBS without resuspending the pellet. If the wash disturbs the pellet, the second wash may be skipped and proceeded to the next. The pellet can then be resuspended in 1 ml of 0.1 N HCl and incubated at room temperature for 5 minutes.

[0087] 2 ml of Tris-HCl (pH 8.0) can be added to a new 15 ml conical centrifuge tube (FALCON™). The fixed cells in HCl above can be transferred to a tube with Tris-HCl to neutralize the HCl. The cell count of the tube can then be calculated as discussed above (e.g., with a hemocytometer or flow cytometer). The fixed cells in Tris-HCl can be spun down at 500×g for 5 minutes and the liquid can be aspirated without disturbing the pellet. The pellet can be washed twice with 1 ml of RNase-free molecular grade water without disturbing the pellet. The cells can then be resuspended to a concentration of 2.5 million cells / ml (the concentration calculated before the last centrifugation step can be used to do this).

[0088] A reverse transcription mix can be made (55 μl M-MuLV Reverse Transcriptase Buffer (ENZYMATICS®), 55 μl M-MuLV Reverse Transcriptase (ENZYMATICS®), 5.5 μl dNTPs (25 mM per base), 3.44 μl RNase Inhibitor (ENZYMATICS®, 40 units / μl), 210.4 μl nuclease-free water, and 2.75 μl RT primer (BC_0055, 100 μM)). Combine and pipette 300 μl of reverse transcription mix with 200 μl of fixed cells (approximately 500,000 cells) in a well of a 24-well cell culture plate. The mixture can then be incubated at room temperature for 10 minutes to allow the reverse transcription primer to anneal, after which the mixture can be incubated at 37° C. overnight (i.e., approximately 16 hours) in a humidified incubator.

[0089] A primer (BC_0055) that can be used for reverse transcription is shown in FIG. 10. This is an anchoring primer designed to bind to the beginning of the poly(A) tail of messenger RNA. A primer can be synthesized with any of the four bases at the 3' end (N) and any base other than T at the 2nd position from 3' (V). The primer can also contain 15 consecutive dTs. In some embodiments, the primer can contain more than 15 dTs. In some other embodiments, the primer can contain less than 15 dTs. In embodiments where the primer contains less than 15 dTs, the melting temperature of the primer may be reduced. The domain s0 does not have to hybridize to the messenger RNA, but can instead provide an accessible binding domain for a linker oligo. The primer also contains a 5' phosphate that allows the primer to be ligated to another oligo by T4 DNA ligase.

[0090] Example 2 - Preparation of Barcodes The barcodes were arranged in an ordered fashion in a 96-well plate at a concentration of 100 μM. Each barcode was annealed to its corresponding linker oligo (see Figures 10-12).

[0091] FIG. 11 shows the annealed first round barcode oligos. 96 first round barcode oligos with unique sequences in domain i8a were used. In the first round, the unique sequence in domain i8a is the region of the sequence that is used as the barcode. By varying 8 nucleotides, there are 65,536 possible unique sequences. In some embodiments, more than 8 nucleotides may be present in domain i8a. In some other embodiments, fewer than 8 nucleotides may be present in domain i8a. The first round barcodes were pre-annealed to the linker strand (BC_0056) via a complementary sequence in domain s1. The linker strand can include a complementary sequence to a portion of the reverse transcription primer (domain s0) that allows the 3' end of the first round barcode to approach and hybridize closest to the 5' end of the reverse transcription primer. The phosphate of the reverse transcription primer can then be ligated to the 3' end of the first round barcode by T4 DNA ligase. Domain s2 can provide an accessible binding domain for a linker oligo for use in another round of barcoding, and the first round barcode oligo can contain a 5' phosphate that can be ligated to the 3' end of another oligo by T4 DNA ligase.

[0092] FIG. 12 shows the annealed second round barcode oligos. 96 second round barcode oligos with unique sequences in domain i8b were used. In the second round, the unique sequence in domain i8b is the region of the sequence that is used as the barcode. By varying 8 nucleotides, there are 65,536 possible unique sequences. In some embodiments, more than 8 nucleotides may be present in domain i8b. In some embodiments, less than 8 nucleotides may be present in domain i8b. The second round barcodes can be pre-annealed to the linker strand (BC_0058) via a complementary sequence in domain s3. The linker strand can include a complementary sequence to a portion of the first round barcode oligo (domain s2) that can bring the 3' end of the first round barcode in close proximity and hybridize to the 5' end of the second round barcode oligo. The phosphate of the first round barcode oligo can then be ligated to the 3' end of the second round barcode by T4 DNA ligase. Domain s4 is a linker for use in another round of barcoding. The second round barcode oligo can provide an accessible binding domain for an oligo, and can contain a 5' phosphate that can be ligated to the 3' end of another oligo by T4 DNA ligase.

[0093] FIG. 13 shows the annealed third round barcode oligos. 96 third round barcode oligos with unique sequences in domain i8c were used. In the third round, the unique sequence in domain i8c is the region of the sequence that is used as the barcode. By varying 8 nucleotides, there are 65,536 possible unique sequences. In some embodiments, more than 8 nucleotides may be present in domain i8c. In some other embodiments, fewer than 8 nucleotides may be present in domain i8c. The third round barcodes can be pre-annealed to the linker strand (BC_0060) via a complementary sequence in domain s5. The linker strand can include a complementary sequence to a portion of the second round barcode oligo (domain s4) that can bring the 3' end of the second round barcode closest to and hybridize to the 5' end of the third round barcode oligo. The phosphate of the second round barcode oligo can then be ligated to the 3' end of the third round barcode by T4 DNA ligase. Third-round barcode oligos can be synthesized with unique molecular identifiers (UMIs; see Islam et al., Nature Methods, 2014) consisting of 10 random nucleotides (domain UMI: NNNNNNNNNN). Multiple sequencing reads can occur from a cDNA due to PCR amplification bias. Using UMIs, each cDNA can be counted only once. Third-round barcodes can also contain domains that correspond to parts of ILLUMINA® TruSeq adapters. Third-round barcodes can be synthesized with a biotin molecule at the 5' end to allow sufficient isolation of barcoded cDNAs with streptavidin-coated magnetic beads.

[0094] Starting with a 100 μM stock of each barcode oligo (i.e., one 96-well plate for each round), 11 μl of barcode oligos were transferred to a 96-well PCR plate. 9 μl of BC_0056 (100 μM stock) was added to each well of the plate with the round 1 barcode. 9 μl of BC_0058 (100 μM stock) was added to each well of the plate with the round 2 barcode. 9 μl of BC_0060 (100 μM stock) was added to each well of the plate with the round 3 barcode. Each plate was then placed in a thermocycler with the following program: heating to 90° C., reducing heat at 0.1° C. / sec, and stopping when the temperature reached 25° C., to allow annealing of the barcodes to the corresponding linker oligos. 2.2 μl was transferred from each well with the round 1 barcode to a new 96-well plate (referred to as plate L1). 3.8 μl was transferred from each well with a round 2 barcode to a new 96-well plate (designated plate L2). 6.1 μl was transferred from each well with a round 3 barcode to a new 96-well plate (designated plate L3).

[0095] Example 3 - Preparation of ligation stop oligos After each round of ligation, the ligation can be stopped by adding an excess of oligos complementary to the linker strand (see Figure 14). To stop each barcode ligation, an oligo strand perfectly complementary to the linker oligo can be added. Such an oligo can bind to the linker strand attached to the unligated barcode and displace the unligated barcode by a strand displacement reaction. The unligated barcode can then be made completely single stranded. The ligation reaction will stop proceeding because T4 DNA ligase cannot ligate single stranded DNA to other single stranded DNA. To ensure that all linker oligos are bound to their complementary oligos, a molar excess of complementary oligos (relative to the linker oligo) is added. To stop the first round ligation, add BC_0064 (complementary to BC_0056). To stop the second round ligation, add BC_0065 (complementary to BC_0058). To stop the third round ligation, add BC_0066 (complementary to BC_0060).

[0096] Dilutions of each stopped ligation strand (BC_0064, BC_0065, BC_0066) can be prepared as follows: 264 μl stopped ligation strand (BC_0064, BC_0065, BC_0066), 300 μl 10× T4 DNA ligase buffer, and 636 μl nuclease-free water.

[0097] Example 4 - Ligation of barcodes to cDNA 5 μl of 10% TRITON™ X-100 can be added to the 24-well plate reverse transcription reaction described above (to a final concentration of 0.1%). The reverse transcription (RT) reaction with cells can be transferred to a 15 ml conical centrifuge tube (FALCON™). The RT reaction can be centrifuged at 500×g for 10 minutes and resuspended in 2 ml of nuclease-free water. The cells can be combined with the ligase mix (600 μl of 10× T4 ligase buffer, 2040 μl of nuclease-free water, all of the resuspended cells (2000 μl), 100 μl of T4 DNA ligase (NEW ENGLAND BIOLABS®, 400,000 units / ml), and 60 μl of 10% TRITON™ X-100 in a disposable pipette reservoir (10 ml). The cells and ligase mix can be mixed by gently tilting the reservoir back and forth several times. Using a multichannel pipette, 40 µl of cells in ligase mix can be added to each well of the annealed round 1 barcodes (plate L1). Each well can be mixed by gently pipetting up and down 2-3 times. The cells in the ligase mix can be incubated at 37 °C for 60 min.

[0098] 10 μl of diluted BC_0064 can be added to each well to stop the ligation. The samples can then be incubated at 37° C. for 30 minutes. All of the cells can be collected into a new disposable pipette reservoir (10 ml). Using a 1 ml pipette, the cells can be passed through a 40 μM strainer into a new disposable pipette reservoir (10 ml). 100 μl of T4 DNA ligase (NEW ENGLAND BIOLABS®, 400,000 units / ml) can be added to the cells in the reservoir. The cells and ligase mix can be mixed by gently tilting the reservoir back and forth several times, and using a multichannel pipette, 40 μl of cells in ligase mix can be added to each well of the annealed round 2 barcodes (plate L2). Each well can be mixed by gently pipetting up and down 2-3 times, and then the samples can be incubated at 37° C. for 60 minutes.

[0099] 10 μl of diluted BC_0065 can be added to each well to stop the ligation. Samples can be incubated at 37° C. for 30 minutes, after which the cells can be collected into a new disposable pipette reservoir (10 ml). Using a 1 ml pipette, the cells can be passed through a 40 μM strainer into a new disposable pipette reservoir (10 ml). 100 μl of T4 DNA ligase (NEW ENGLAND BIOLABS®, 400,000 units / ml) can be added to the cells in the reservoir. The cells and ligase mix can be mixed by gently tilting the reservoir back and forth several times. Using a multichannel pipette, 40 μl of cells in ligase mix can be added to each well of the annealed round 3 barcodes (plate L3). Each well can then be mixed by gently pipetting up and down 2-3 times, and the samples are incubated at 37° C. for 60 minutes. The mixture can be incubated.

[0100] 10 μl of diluted BC_0066 can be added to each well to stop the ligation. Samples can be incubated at 37° C. for 30 minutes. All cells can be collected in a new disposable pipette reservoir (10 ml). Cells can be transferred to a 15 ml conical centrifuge tube (FALCON™) and the tube can be filled to 15 ml with wash buffer (nuclease-free water, 0.05% Tween 20, and 25% formamide). Samples can be incubated at room temperature for 15 minutes. Cells can then be pelleted at 500×g for 10 minutes and the liquid can be removed without disturbing the pellet. Cells in each tube can be resuspended in 100 μl of PBS and cells can be counted (e.g., with a hemocytometer or flow cytometer). In one example, 57,000 cells were retained. The number of cells to be sequenced can be selected. In one example, cells were distributed into aliquots of 25 cells, 250 cells, 2,500 cells, and 25,000 cells. 300 μl of lysis buffer (10 mM NaF, 1 mM Na3VO4, 0.5% DOC buffer, and 0.5% TRITON™ X-100) can be added to each of the cell aliquots, and each of the cell aliquots can be passed through a 25-gauge needle 8 times.

[0101] Example 5 - Binding of barcoded cDNA to streptavidin-coated beads First, the DYNABEADS® MYONE™ Streptavidin C1 beads can be resuspended. 20 μl of resuspended DYNABEADS® MYONE™ Streptavidin C1 beads (for each aliquot of cells) can be added to a 1.7 ml microcentrifuge tube (EPPENDORF®). The beads can be washed three times with 1× Phosphate Buffered Saline Tween 20 (PBST) and resuspended in 20 μl of PBST. 900 μl of PBST can be added to the cell aliquot and 20 μl of washed C1 beads can be added to the lysed cell aliquot. The samples can be placed on a gentle roller at room temperature for 15 minutes and then washed three times with 800 μl of PBST using a magnetic tube rack (EPPENDORF®). The beads can then be resuspended in 100 μl of PBS.

[0102] Example 6 - RNase treatment of beads The microcentrifuge tube (EPPENDORF®) containing the sample can be placed in a magnetic tube rack (EPPENDORF®) for 2 minutes and then the liquid can be aspirated. The beads can be resuspended in RNase reaction (3 μl RNase Mix (ROCHE™), 1 μl RNase H (NEW ENGLAND BIOLABS®), 5 μl RNase H 10× Buffer (NEW ENGLAND BIOLABS®), and 41 μl nuclease-free water). The sample can be incubated at 37° C. for 1 hour, removed from 37° C., and placed in a magnetic tube rack (EPPENDORF®) for 2 minutes. Without resuspending the beads, the sample can be washed with 750 μl nuclease-free water + 0.01% Tween 20 (H2O-T) while the tube is still in the magnetic tube rack. The liquid can then be aspirated. Without resuspending the beads, the sample can be washed with 750 μl of HO-T while the tube is still in the magnetic tube rack. The liquid can then be aspirated while the tube is still in the magnetic tube rack. The tube can then be removed from the magnetic tube rack and the sample can be resuspended in 40 μl of nuclease-free water.

[0103] Example 7 - 3' Adapter Ligation Referring to FIG. 15, a single-stranded DNA adaptor oligo (BC_0047) can be ligated to the 3' end of the cDNA to facilitate PCR amplification. Dideoxycytidine (ddC) can be included at the 3' end of the adaptor oligo to prevent concatemerization of the adaptor oligo. BC_0047 was generated with a phosphate at the 5' end and a ddC at the 3' end. Several enzymes are capable of ligating single-stranded oligos to the 3' end of single-stranded DNA. T4 RNA Ligase 1 (NEW ENGLAND BIOLABS®) was used herein. Thermostable 5' AppDNA / RNA Ligase (NEW ENGLAND BIOLABS®) can also be used with pre-adenylated adaptor oligos.

[0104] Specifically, 20 μl of RNase-treated beads can be added to a single PCR tube. 80 μl of ligase mix (5 μl of T4 RNA ligase 1 (NEW ENGLAND BIOLABS®), 10 μl of 10×T4 RNA ligase buffer, 5 μl of 50 μM BC_0047 oligo, 50 μl of 50% PEG8000, and 10 μl of 10 mM ATP) can be added to the 20 μl of beads in the PCR tube. The 50 μl of ligase mixed with the beads can be transferred to a new PCR tube to prevent too many beads from settling to the bottom of a single tube, and the sample can be incubated at 25° C. for 16 hours.

[0105] Example 8 - Generation of ILLUMINA®-compatible sequencing products The ligation reactions from both PCR tubes can be combined into a single 1.7 ml microcentrifuge tube (EPPENDORF®). 750 μl of H2O-T can be added to each sample. Each of the tubes can be placed in a magnetic tube rack (EPPENDORF®) for 2 minutes, the liquid can be aspirated, and the samples can be resuspended in 40 μl of water. The samples can be transferred to PCR tubes. 60 μl of PCR mix can be added to each tube (50 μl of 2×PHUSION® DNA Polymerase Master Mix (THERMO FISHER™ Scientific), 5 μl of BC_0051 (10 μM), and 5 μl of BC_0062 (10 μM)). 10 cycles of PCR can be performed (98° C. for 3 minutes, 10 replicates (98° C. for 10 seconds, 65° C. for 15 seconds, and 72° C. for 60 seconds), and 72° C. for 5 minutes). Figure 16 shows the PCR product. After ligating the 3' adaptor oligo (BC_0047) to the barcoded cDNA, PCR can be used to amplify the cDNA. Primers BC_0051 and BC_0062 were used, as shown in Figure 16.

[0106] The PCR samples can be sourced from the previous step and the magnetic beads can be moved to the bottom of each tube using a magnet. 90 μl of the PCR reaction can be transferred to a new 1.7 ml without transferring any of the magnetic beads. 10 μl of nuclease-free water can be added to each of the 1.7 ml tubes to bring the total volume to 100 μl. 60 μl of AMPURE™ beads can be added to the 100 μl PCR reaction (0.6×SPRI) and allowed to bind for 5 minutes. The tubes can be placed on the magnet for 2 minutes and, without resuspending the beads, the samples can be washed with 200 μl of 70% ethanol (wait 30 seconds). Without resuspending the beads, the samples can be washed again with 200 μl of 70% ethanol (wait 30 seconds) and then the samples can be air-dried for 5-10 minutes until the ethanol evaporates.

[0107] Each sample can be resuspended in 40 μl of nuclease-free water. The tubes can be placed on the magnetic rack for 2 minutes. With the microcentrifuge tubes (EPPENDORF®) still in place on the magnetic rack, 38 μl of solution can be transferred to a new 1.7 ml tube without transferring the beads. 62 μl of nuclease-free water can be added to the samples. A total volume of 100 μl can be added. Then, 60 μl of AMPURE™ beads can be added to the 100 μl PCR reaction (0.6×SPRI) and allowed to bind for 5 minutes. The tube can be placed on a magnet for 2 minutes, after which, without resuspending the beads, the sample can be washed with 200 μl of 70% ethanol (wait 30 seconds). Without resuspending the beads, the sample can be washed again with 200 μl of 70% ethanol (wait 30 seconds), after which the sample can be air-dried for 5-10 minutes until the ethanol has evaporated.

[0108] The samples can be resuspended in 40 μl of nuclease-free water and each tube can be placed on a magnetic rack for 2 minutes. With the tubes still in place on the magnetic rack, 38 μl of the solution can be transferred to a new 1.7 ml tube without transferring any of the beads. 20 μl of the 38 μl eluate can be added to an optical PCR tube. Additionally, PCR mix can be added to this tube (25 μl PHUSION® DNA Polymerase Master Mix (THERMO FISHER™ Scientific), 2.5 μl BC_0027 (10 μM), 2.5 μl BC_0063 (10 μM), and 2.5 μl 20×EVAGREEN® (BIOTIUM™)). After the PCR shown in FIG. 16, the full-length ILLUMINA® adapter sequence can be introduced by another round of PCR. As shown in Figure 17, BC_0027 contains the flow cell binding sequence and the binding site for the TRUSEQ™ Read 1 primer. BC_0063 contains the flow cell binding sequence and the TruSeq multiplex Read 2 and index binding sequence. There is also a region for the sample index, which in this example is GATCTG.

[0109] The above samples can be run on a qPCR machine with the following cycle conditions: 1) 98°C for 3 minutes, 2) 98°C for 10 seconds, 3) 65°C for 15 seconds, 4) 72°C for 60 seconds, and 5) repeat steps 2-4 (e.g., 10-40 times, depending on when the exponential increase in fluorescence stops). The tubes can be transferred to a thermocycler set at 72°C for 5 minutes. The qPCR reactions can be run on a 1.5% agarose gel for 40 minutes, and a 450-550 bp band can be excised and gel extracted (QIAQUICK® Gel Extraction Kit). The products can be sequenced using paired-end sequencing on an ILLUMINA® MISEQ™. The sequencing primers can be standard TRUSEQ™ multiplex primers. Read 1 can sequence the cDNA sequence and read 2 can cover the unique molecular identifier as well as three barcode sequences (8 nucleotides each). Index read 1 can be used to sequence the sample barcode, and thus multiple samples may be sequenced together.

[0110] Example 9 - Data Analysis Sequencing reads were grouped by cell barcode (3 barcodes of 8 nucleotides each, 96 x 96 x 96 = 884,736 total combinations). Each barcode combination should correspond to cDNA derived from a single cell. Only reads with valid barcodes were kept. Sequencing reads with each barcode combination were aligned to both the human and mouse genomes. Reads aligned to both genomes were discarded. Multiple reads with the same unique molecular identifier were counted as a single read. Reads with unique molecular identifiers with two or fewer mismatches were considered to have been generated by sequencing errors and were counted as a single read. For each unique barcode combination, the number of reads aligned to the human genome (x-axis) and the number of reads aligned to the mouse genome (y-axis) were plotted (see Figure 18). Each cell is either mouse or human and therefore should ideally contain only one type of RNA. Therefore, an ideal plot would have all points along the x-axis or y-axis. The plot in Figure 18 The fact that most of the points in the plot are near the axis indicates that the method is feasible.

[0111] Each point on the plot corresponds to a cDNA with the same combination of barcodes and should represent the cDNA derived from a single cell. For each point, the number of reads that map uniquely to the mouse genome is plotted on the y-axis, and the number of reads that map uniquely to the human genome is plotted on the x-axis. If the cDNAs with a particular combination of barcodes are derived from a single cell, then all the cDNAs with a particular combination of barcodes should map completely to the human genome or completely to the mouse genome. As mentioned above, the fact that most barcode combinations map near either the x-axis (human cells) or the y-axis (mouse cells) indicates that the method can indeed produce single-cell RNA sequencing data.

[0112] Example 10 - Method for uniquely labeling molecules in multiple cells For the protocol disclosed below, the expected experiment time is 2 (two) days. RNase inhibitors may be added to the buffers as indicated below. Thus, if any buffer contains the term "+RI", this indicates that ENZYMATICS® RNase inhibitor should be added at a final concentration of 0.1 U / μL. The centrifugation step may be performed in a swinging bucket rotor. In some embodiments, the use of a fixed angle centrifuge may result in more cell loss. Depending on the tissue type, centrifugation speeds may need to be altered to optimize cell retention (e.g., smaller cells = faster).

[0113] DNA barcoded plate generation may require: i) three IDT® 96-well plates, reverse transcription barcode primer, ligation round 1, and ligation round 2 stock DNA oligo plates (100 μM); ii) two linker oligos, BC_0215 and BC_0060 (note: these are assumed to be at 1 mM stock concentration, so modify volumes if another stock concentration is used (e.g., 100 μM stock); and iii) six 96-well PCR plates (e.g., three stock plates that will last at least 10 experiments, and three plates for the first experiment). Note that this allows 100 μL of DNA barcodes to be generated in each well. Generally, each experiment requires only 4 μL / well of reverse transcription primer solution, which can last 25 experiments. Generally, each experiment requires only 10 μL / well of barcode / linker solution, which can last a total of 10 experiments on these plates.

[0114] Round 1 reverse transcription barcode primers (final concentrations in each of the 48 wells are 12.5 μM random hexamer and 12.5 μM 15dT primer): 1) Using a multichannel pipette, add 12.5 μL of rows A-D of IDT® reverse transcription barcode primers to rows A-D of a BC stock 96-well PCR plate; 2) Using a multichannel pipette, add 12.5 μL of rows E-H of IDT® reverse transcription barcode primers to rows A-D of a BC stock 96-well PCR plate (where poly dT is mixed with the random hexamer primers); 3) Add 75 μl of water to rows A-D of a BC stock 96-well PCR plate.

[0115] Round 2 ligation round (final concentrations are 12 µM barcode, 11 µM linker-BC_0215): 1) Using a multichannel pipette, add 12 µL of IDT® Round 2 barcode to the R1 stock 96-well PCR plate; 2) Add 138.6 µl of BC_0215 (1 mM) to 10.9494 mL of water in the wells. (BC_0215_dil); and 3) using a multichannel pipette, add 88 µL of BC_0215_dil to each well of the R2 stock 96-well PCR plate.

[0116] Ligation round 3 (final concentrations are 14 µM barcode, 13 µM linker-BC_0060): 1) using a multichannel pipette, add 14 µL of round 3 barcode to an R3 stock 96-well PCR plate; 2) add 163.8 µl of BC_0060 (1 mM) to 10.6722 mL of water in the recess (BC_0060_dil); and 3) using a multichannel pipette, add 86 µL of BC_0060 to each well of an R3 stock 96-well PCR plate.

[0117] For each ligation plate (R2 and R3, not containing the reverse transcribed barcodes), anneal the barcode and linker oligos with the following thermocycling protocol: 1) heat to 95 °C for 2 (two) minutes, and 2) cool to 20 °C at a rate of 0.1 °C / sec; and 3) at 4 °C.

[0118] Aliquot 10 µL of each barcode / linker stock into three (3) new 96-well PCR plates. These are the plates to be used for DNA barcoding in the split-pool ligation step of this protocol.

[0119] I. Nuclear Extraction (Optional): 1) Prepare the following: a) Keep the Dance homogenizer at 4° C. until use, b) 15 ml 1× PBS + 37.5 SUPERASE-IN™ + 19 μl ENZYMATICS® RNase inhibitor (keep on ice), and c) pre-cool the centrifuge to 4° C.

[0120] 2) Make NIM1 buffer (Table 1):

[0121] [Table 1]

[0122] 3) Make homogenization buffer (Table 2):

[0123] [Table 2]

[0124] 4) Dance homogenizer: a) add tissue / cell sample to Dance homogenizer; for cells, resuspend in 700 μl homogenization buffer; b) add homogenization buffer to approximately 700 μl; c) perform 5 strokes with loose pestle; d) perform 10-15 strokes with tight pestle; e) add homogenization buffer to 1 ml; and f) check cell lysate with 5 μl trypan blue and 5 μl cells in a hemocytometer to see if nuclei have been released.

[0125] 5) Filter the homogenate through a 40 μm strainer into a 5 ml EPPENDORF™ tube (or a 15 mL FALCON™ tube). Tilt the filter at 45° while straining above the tube to ensure that the lysate passes through as intended. Note: This filtering process is different from the filtering process described below.

[0126] 6) Centrifuge at 600xg for 4 minutes (4°C) and remove the supernatant (you may leave about 20μL to avoid pellet aspiration). 7) Resuspend in 1ml 1xPBS+RI. 8) Add 10μl BSA. 9) Centrifuge at 600xg for 4 minutes. 10) Resuspend in 200μl 1xPBS+RI.

[0127] 11) Take 50 μl of resuspended cells from step 4 and add 150 μl of 1×PBS+RI. Count samples on a hemocytometer and / or flow cytometer. The volume of resuspended cells in step 4 can be varied based on user discretion. 12) Pass cells through a 40 μm strainer into a new 15 mL Falcon™ tube and place on ice (see notes below for step 4 of II. Fixation and Permeabilization). 13) Resuspend the desired number of nuclei (typically 2M) in 1 mL of 1×PBS+RI and proceed to step 5 of the Fixation and Permeabilization protocol below.

[0128] II. Fixation and Permeabilization: 1) Prepare the following buffers (calculated for 2 experiments): a) 1.33% formalin (360 μL of 37% formaldehyde in water (SIGMA®) + 9.66 ml of PBS) solution, stored at 4°C; b) 6 mL of 1x PBS + RI (15 μL of SUPERASE-IN™ and 7.5 μL of ENZYMATICS® RNase inhibitor); c) 2 mL of 0.5x PBS + RI (5 μL of SUPERASE-IN™ and 2.5 μL of ENZYMATICS® RNase inhibitor); d) 500 μL of 5% TRITON™ X-100 + RI (2 μL of SUPERASE-IN™); e) 500 μL of 100 mM Tris pH 8.0 + 2 μL of SUPERASE-IN™; and f) set the centrifuge to 4° C.

[0129] 2) Pellet the cells by centrifugation at 500 x g for 3 minutes at 4°C (some cells may require a faster centrifugation). 3) Resuspend the cells in 1 mL of cold PBS+RI. Keep the cells on ice between these steps. 4) Pass the cells through a 40 μm strainer into a new 15 mL Falcon™ tube and place on ice. NOTE: The cell resuspension will most likely not pass through the strainer passively, which may cause cell loss. Instead, fill a 1 mL pipette with the resuspension and press the tip end directly against the strainer to actively push the liquid out. This action should take approximately 1 (one) second. 5) Add 3 mL of cold 1.33% formaldehyde (final concentration is 1% formaldehyde). Fix the cells on ice for 10 minutes. 6) Add 160 μL of 5% TRITON™ X-100 + RI to the fixed cells and mix by gently pipetting up and down 5 times with a 1 mL pipette. Permeabilize cells for 3 minutes on ice. 7) Centrifuge cells at 500×g for 3 minutes at 4° C. 8) Carefully aspirate cells and resuspend in 500 μL of cold PBS + RI. 9) Add 500 μL of cold 100 mM Tris-HCl, pH 8.0. 10) Add 20 μL of 5% TRITON™ X-100. 11) Centrifuge cells at 500×g for 3 minutes at 4° C. 12) Aspirate cells and resuspend in 300 μL of cold 0.5× PBS + RI.

[0130] 13) Pass the cells through a 40 μM strainer into a new 1.7 mL tube (see notes above for step 4 in II. Fixation and Permeabilization). 14) Count the cells using a hemocytometer or flow cytometer and dilute the cell suspension to 1,000,000 cells / mL. Keep the cell suspension on ice while counting the cells. NOTE: This step will determine the number of cells that will proceed to the partition-pool round. It will be possible to sequence only a subset of the cells that will proceed to the partition-pool round (this can be done during sub-library generation in the lysis step). The total number of barcode combinations that will be used should be calculated to determine the maximum number of cells that can be sequenced with minimal barcode collisions. Without wishing to be bound by any one particular theory, the number of cells that will be processed should not exceed 5% of the total barcode combinations. Generally, dilutions of 500k to 1M cells / mL can be used here (equivalent to 4 to 8k cells going into each well of the reverse transcription barcoding round).

[0131] III. Reverse Transcription: 1) Aliquot 4 μL of the RT barcode stock plate into the top 4 rows (48 wells) of a new 96-well plate. Cover the plate with an adhesive plate seal until ready to use.

[0132] 2) Generate the following reverse transcription (RT) mix on ice (Table 3):

[0133] [Table 3]

[0134] 3) Add 8 μL of RT mix to each of the top 48 wells. Each well should now contain a volume of 12 μL. 4) Add 8 μL of cells in 0.5×PBS+RI to each of the top 48 wells. Each well should now contain a volume of 20 μL. 5) Add the plate to a thermocycler with the following protocol: a) 50° C. for 10 minutes; b) 3 (three) cycles of i) 8° C. for 12 seconds, ii) 15° C. for 45 seconds, iii) 20° C. for 45 seconds, and iv) 30° C. for 30 seconds, v) 42° C. for 2 minutes, vi) 50° C. for 3 minutes; c) 50° C. for 5 minutes; and d) always 4° C. thereafter.

[0135] 6) Place RT plate on ice. 7) Prepare 2 mL of 1x NEB Buffer 3.1 with 20 μL ENZYMATICS® RNase inhibitor. 8) Transfer each RT reaction to a 15 mL FALCON™ tube (again on ice). 9) Add 9.6 μL of 10% TRITON™ X-100 to give a final concentration of 0.1%. 10) Centrifuge pooled RT reactions at 500 x g for 3 minutes. 11) Aspirate supernatant and resuspend in 2 mL of 1x NEB Buffer 3.1 + 20 μL ENZYMATICS® RNase inhibitor.

[0136] IV. Ligation Barcoding: 1) Make the following ligation master mix on ice (Table 4):

[0137] [Table 4]

[0138] 2) Add 2 mL of cells in NEB Buffer 3.1 to the ligation mix. This mix should now have a volume of 4.04 mL. 3) Add the mix to the recessed areas. 4) Using a multichannel pipette, add 40 μL of ligation mix (with cells) to each well of the round 1 DNA barcode plate. 5) Cover the round 1 DNA barcode plate with an adhesive plate seal and incubate at 37 °C for 30 minutes with gentle rotation (50 rpm). 6) Make the round 1 blocking solution and add it to the new recessed areas (Table 5).

[0139] [Table 5]

[0140] 7) Remove the Round 1 DNA barcoded plate from the incubator and remove the cover. 8) Using a multichannel pipette, add 10 μL of Round 1 Blocking Solution to each of the 96 wells of the Round 1 DNA barcoded plate. 9) Cover the Round 1 DNA barcoded plate with an adhesive plate seal and incubate at 37°C for 30 minutes with gentle rotation (50 rpm). 10) Remove the Round 1 DNA barcoded plate from the incubator, remove the cover, and pool all the cells into a new well. 11) Pass all the cells from this well through a 40 μm strainer into another well (see note above on step 4 of Fixation and Permeabilization). 12) Add 100 μL of T4 DNA Ligase to this well and mix by pipetting approximately 20 times. 13) Using a multichannel pipette, add 50 μL of cell / ligation solution to each well of the Round 2 DNA barcoded plate. 14) Cover the round 2 DNA barcoded plate with an adhesive plate seal and incubate with gentle rotation (50 rpm) at 37°C for 30 minutes.

[0141] 15) Round 2 Make blocking solution and add it to the new recessed areas (Table 6).

[0142] [Table 6]

[0143] 16) Remove the Round 2 DNA barcoded plate from the incubator and remove the cover. 17) Using a multichannel pipette, add 20 μL of Round 2 Blocking and Stopping Solution to each of the 96 wells of the Round 2 DNA barcoded plate. 18) Pool all cells into a new well (do not incubate at the final blocking step). 19) Pass all cells from this well through a 40 μm strainer into a 15 mL FALCON™ tube (see note above on Fixation and Permeabilization step 4). 19) Count cells on a flow cytometer. Make sure cells are well mixed before aliquoting samples for counting.

[0144] V. Lysis: 1) Make 2x Lysis Buffer (Table 7).

[0145] [Table 7]

[0146] 2) If a white precipitate appears, warm at 37°C until the precipitate returns to solution (approximately 10 to 15 minutes).

[0147] 3) Make the following wash buffer (Table 8).

[0148] [Table 8]

[0149] 4) Add 70 μL of 10% TRITON™ X-100 to the cells (approximately 0.1% final concentration). 5) Centrifuge in a 15 ml tube at 1000×g for 5 minutes. Note: The pellet in the steps below may be very small and may not be visible. 6) Aspirate the supernatant, leaving about 30 μl to avoid removing the pellet, a) if possible, remove as much of the supernatant as possible using a 20 μL pipette. 7) Resuspend in 4 mL of wash buffer. 8) Centrifuge at 1000×g for 5 minutes. 9) Aspirate the supernatant and resuspend in 50 μl of 1×PBS+RI.

[0150] 10) Dilute 5 μl into 195 μl 1×PBS and count by flow cytometry. Alternatively, take 5 μl into 5 μl 1×PBS and count by hemocytometer (it can be difficult to distinguish debris from cells). 11) Determine how many sub-libraries to generate (number of sub-libraries = number of tubes needed) and how many cells each of those sub-libraries will have. 12) Aliquot the desired number of cells of each sub-library into new 1.7 mL tubes. Add 1×PBS to each tube to bring the final volume to 50 μL.

[0151] 13) Add 50 μL of 2x lysis buffer to each tube. 14) Add 10 μL of proteinase K (20 mg / mL) to each lysate. 15) Incubate for 2 hours at 55°C with shaking at 200 rpm. 16) Stopping point (optional): Freeze lysates at -80°C.

[0152] VI. Prepare buffer solutions. First, make the following stock solutions (Tables 9-11).

[0153] [Table 9]

[0154] [Table 10]

[0155] [Table 11]

[0156] Then make the following smaller aliquots (RNase inhibitors have been added; Tables 12-14):

[0157] [Table 12]

[0158] [Table 13]

[0159] [Table 14]

[0160] VII. Purification of cDNA NOTE: Agitation steps are performed on a vortexer at a low setting (2 / 10) using a foam 1.7 mL tube holder.

[0161] Wash MYONE™ C1 DYNABEADS®: 1) For each lysate to be processed, add 44 μL of MYONE™ C1 DYNABEADS® to a 1.5 mL tube (e.g., 1 lysate = 44 μL, 2 lysates = 88 μL, 3 lysates = 132 μL, etc.). 2) Add 800 μL of 1x B&W-T buffer. 3) Place samples on magnetic rack and wait until liquid clears (1-2 min). 4) Remove supernatant and resuspend beads in 800 μL of 1x B&W-T buffer. 5) Repeat steps 3-4 two more times for a total of 3 washes. 6) Place samples on magnetic rack and wait until liquid clears. 7) Resuspend beads in 100 μL (per sample) of 2x B&W buffer + RI.

[0162] Binding of samples to streptavidin: 1) Add 5 μL of 100 μM PMSF to each sample and leave at room temperature for 10 minutes. 2) Add 100 μl of resuspended C1 beads to each tube. 3) Stir for 60 minutes at room temperature to allow cDNA to bind to the C1 beads. 4) Place samples on magnetic rack and wait until liquid is clear (1-2 minutes). 5) Remove supernatant and resuspend beads in 250 μL of 1x B&W. 6) Stir beads for 5 minutes at room temperature. 7) Repeat steps 5 and 6. 8) Remove supernatant and resuspend beads in 250 μL of 10 mM Tris+T. 9) Stir beads for 5 minutes at room temperature. 10) Leave beads in final wash solution on ice.

[0163] Template switch. Depending on the number of samples, prepare the following mixes (Table 15).

[0164] [Table 15]

[0165] 1) Place samples on magnetic rack and wait until liquid is clear. 2) With samples still on magnetic rack, remove supernatant and wash with 250 μL water (do not resuspend beads this time). 3) Resuspend samples in 200 μl template switch mix. 4) Incubate at room temperature with stirring or rotation for 30 minutes. 5) Incubate at 42°C with stirring or rotation (shaking incubator at 100 rpm) for 90 minutes.

[0166] 6) Potential Stopping Points: If you wish to stop, do the following (otherwise skip to next section): a) place samples on magnetic rack and wait until liquid is clear, b) resuspend in 250 μL Tris-T.

[0167] VIII. cDNA Amplification. Depending on the number of samples, prepare the following PCR mix (Table 16).

[0168] [Table 16]

[0169] 1) Place samples on magnetic rack and wait until liquid clears. 2) While samples are still on magnetic rack, wash (do not resuspend) with 250 μL nuclease-free water. 3) Resuspend samples in 220 μL PCR mix and distribute equally among 4 (four) different PCR tubes.

[0170] 4) Run the following thermocycling program: a) 95 °C, 3 min; b) 98 °C, 20 s; c) 65 °C, 45 s; d) 72 °C, 3 min; e) repeat (b-d) four times (total of 5 cycles); and f) 4 °C, hold.

[0171] 5) Combine all four (4) reactions in a single 1.7 mL tube. Be sure to resuspend any beads that may be stuck to the bottom or sides of the PCR tube before combining the reactions. 6) Place samples on a magnetic rack and wait until the liquid runs clear. 7) Transfer 200 μL of the supernatant to four (4) optical grade qPCR tubes (50 μL in each tube). 8) Add 2.5 μL of 20× EVAGREEN® to each qPCR tube. 9) Run the following qPCR program (be sure to remove samples once the signal begins to go beyond log phase to prevent over-amplification): a) 95°C, 3 minutes; b) 98°C, 20 seconds; c) 67°C, 20 seconds; d) 72°C, 3 minutes; e) repeat (b-d) until the signal plateaus from exponential amplification; f) 72°C, 5 minutes; and g) 4°C, hold.

[0172] 10) Optional: Run an agarose gel or bioanalyze the resulting qPCR. There will likely be a combination of cDNA and dimers present.

[0173] SPRI Size Selection (0.8x). 1) Combine qPCR reactions in a single tube. 2) Remove 180 μL of pooled qPCR reactions and place in a new 1.7 mL tube. 3) Add 144 μL of KAPA™ Pure Beads to the tube and vortex briefly to mix. Wait 5 (five) minutes to bind DNA. 4) Place tube on magnetic rack and wait until liquid is clear. 5) Remove supernatant. 6) With tube still on magnetic rack, wash with 750 μL of 85% ethanol. Do not resuspend beads. 7) Repeat step 6.

[0174] 8) Remove the ethanol and air dry the beads (approximately 5 minutes). Do not over-dry the beads as this may cause them to crack. 9) Resuspend the beads in each tube in 20 μL of water. Once the beads are completely resuspended in water, incubate the tubes at 37°C for 10 minutes. 10) Place the tubes on a magnetic rack and wait until the liquid is clear. 11) Transfer 18.5 μL of the eluate to a new optical grade PCR tube. 12) Run a bioanalyzer trace on 10 μL of the eluate.

[0175] 13) If no dimers are present after size selection, proceed directly to the "Tagmentation and ILLUMINA® Amplicon Generation" section below. If dimers are still present, proceed to step 14 to perform a second amplification and size selection step. This may be necessary for cells with low RNA content, but should not be necessary for cells with high RNA content (e.g., HeLa-S3, NIH / 3T3, etc.).

[0176] Second qPCR (optional). 14) Depending on the number of samples, make the following PCR mix (Table 17).

[0177] [Table 17]

[0178] 15) Add 31.5 μL of qPCR master mix to each optical PCR tube with a previous PCR sample. Flick the tube gently to mix and centrifuge briefly in a tabletop centrifuge to remove any air bubbles. 16) Run the following qPCR program (be sure to remove samples once the signal begins to pass the log phase to prevent over-amplification): a) 95°C, 3 min; b) 98°C, 20 sec; c) 67°C, 20 sec; d) 72°C, 3 min; e) repeat (b-d) until the signal goes from exponential amplification to plateauing; f) 72°C, 5 min; and g) 4°C, hold.

[0179] 17) Run an agarose gel or bioanalyze the resulting qPCR. While dimers may still be present, amplified cDNA should be clearly visible from 500bp to 2500bp (see Figure 19 left for expected size distribution).

[0180] Second SPRI size selection (0.8x). 18) Combine qPCR reactions in a single tube. 19) Remove 40μL of qPCR reactions and place into new 1.7mL tubes. 20) Add 32μL of KAPA™ Pure beads to each tube and vortex briefly to mix. Wait 5 minutes to allow binding to DNA.

[0181] 21) Place the tubes on the magnetic rack and wait until the liquid runs clear. 22) Remove the supernatant. 23) With the tubes still on the magnetic rack, wash with 750 μL of 85% ethanol. 24) Repeat step 5. 25) Remove the ethanol and air dry the beads (approximately 5 minutes). Do not over-dry the beads, which may cause them to crack. 26) Resuspend the beads in each tube in 20 μL of water and wait 5 minutes. 27) Place the tubes on the magnetic rack and wait until the liquid runs clear. 28) Transfer 18.5 μL of eluate to a 1.7 mL tube. 29) Run an agarose gel or bioanalyze the resulting qPCR. There should be very little dimers present at this point. If any dimers remain, perform another round of qPCR, followed by another round of 0.8× SPRI size selection.

[0182] IX. Tagmentation and ILLUMINA® Amplicon Generation. 1) Quantify amplified cDNA in Quibit and dilute to 0.12ng / μL. 2) Preheat thermocycler to 55°C. 3) For each sample, combine 600pg purified cDNA with H2O for a total volume of 5μl. 4) Add 10μl of Nextera TD buffer and 5μl of Amplicon Tagment enzyme to each tube (total volume of reaction is now 20μl). Pipette approximately 5 times. Mix by pipetting approximately 5 times and centrifuge. 5) Incubate at 55°C for 5 minutes. 6) Add 5 μl of Neutralization Buffer. Mix by pipetting approximately 5 times and centrifuge (air bubbles are not unusual). 7) Incubate at room temperature for 5 minutes.

[0183] 8) Add each PCR to the tubes in the following order: i) 15 μl Nextera PCR mix; ii) 8 μl HO; iii) 1 μl of 10 μM (N7 indexed primer, one of BC_0076 to BC_0083); and iv) 1 μl of 10 μM Nextera (BC_0118) N501 oligo. 9) Run the following thermocycling program: i) 95 °C for 30 s; ii) 12 cycles of a) 95 °C for 10 s; b) 55 °C for 30 s; and c) 72 °C for 30 s; and iii) 72 °C for 5 min and always at 4 °C thereafter.

[0184] 10) Transfer 40 μl of the 50 μL reaction to a 1.7 mL tube. 11) Add 28 μL of KAPA™ Pure beads for 0.7× cleanup. Elute 20 μl. 12) Bioanalyze the resulting sample and quantitate in Quibit prior to sequencing (see Figure 19 right side for expected size distribution).

[0185] X. ILLUMINA® Sequencing. 1) Use paired-end sequencing performed with a 150bp kit. 2) Set read 1 to 66 nt (transcript sequence). 3) Set read 2 to 94 nt (cell-specific barcode and UMI). 4) Prepare sub-library indexes, including a read 1 index of 6 nt.

[0186] Example 11 - Keeping cells on ice for cell fixation and permeabilization on ice A unique molecular index (UMI) is a random molecular barcode that is added to each transcript / cDNA prior to amplification. The unique molecular index allows computational exclusion of PCR duplicates, since they have the same UMI sequence. Thus, even if multiple PCR replicates are sequenced, each original transcript will only be counted once.

[0187] The measure of UMI per cell indicates how many unique RNA molecules can be detected per cell. This directly relates to the degree of efficiency with which molecules are barcoded and processed to allow detection by next-generation sequencing. This measure is generally the gold standard in the field, and most researchers in the single-cell RNA sequencing space determine how successful an experiment is based on the number of UMIs per cell detected for a given amount of raw sequencing reads (e.g., 10,000 UMIs per cell for 50,000 sequence reads per cell).

[0188] For each of Examples 11-16, all conditions within the same table (i.e., Table 18, Table 19, Table 20, Table 21, and Table 22) were sequenced to the same saturation level, allowing accurate comparison across conditions. Each condition had the same number of raw sequencing reads compared to another condition. This indicated that the detected differences in detected UMIs per cell were due to changes in conditions and not sequencing depth.

[0189] In other protocols (see, e.g., Rosenberg, AB, et al., BioRxiv (2017):105163), formaldehyde fixation is performed at room temperature. Experiments were performed herein to show that there is a substantial improvement when cells are kept on ice (e.g., at 4°C) during fixation and permeabilization (see Table 18).

[0190] [Table 18]

[0191] Example 12 - Addition of protease inhibitors (PMSF) and direct binding to streptavidin beads In other protocols (see, e.g., Rosenberg, AB, et al., BioRxiv (2017):105163), nucleic acids are first isolated from the lysis solution using SPRI bead cleanup, and then the desired nucleic acid (containing the 5' biotin) is bound to streptavidin beads. Herein, by adding PMSF to the lysate and then adding streptavidin beads directly (thereby skipping the initial SPRI isolation of the nucleic acid), an improvement in the number of unique molecules detected per cell was found (see Table 19).

[0192] [Table 19]

[0193] Example 13 - 0.6x vs. 0.8x SPRI cleanup In other protocols (see, e.g., Rosenberg, AB, et al., BioRxiv (2017):105163), a 0.6x SPRI cleanup was performed after cDNA amplification. After cDNA amplification, SPRI cleanup can be used to remove short products (i.e., products less than about 200 base pairs). Herein, it was found that adding a 0.8:1 ratio of SPRI beads to PCR products resulted in more unique RNA molecules detected per cell than the previous 0.6:1 ratio of SPRI to PCR products (see Table 20).

[0194] [Table 20]

[0195] Example 14 - Concentration of Random Hexamer and Poly dT RT Primers In other protocols (see, for example, Rosenberg, AB, et al., BioRxiv (2017): 105163), reverse transcription is performed using only poly dT primers for reverse transcription. Herein, it was found that combining barcoded random hexamer primers and varying concentrations resulted in a significant increase in UMI / cell. With reference to Table 21, protocol version 2.1 uses condition number 4 (729 UMI / cell). With continued reference to Table 21, here, condition number 10 was used (1263 UMI / cell). Thus, by combining specific concentrations of random hexamer and poly dT reverse transcription primers together, the efficiency of in situ reverse transcription is increased.

[0196] [Table 21]

[0197] Example 15 - Use of first round barcodes as sample identifiers A barcoded reverse transcription primer was used to generate first round barcodes on the cDNA molecules, which were then further barcoded by adding the barcode to the 5' end by ligation as described above.

[0198] A proof of concept is described in Rosenberg, AB et al., (2018) Science, 360(6385);176-182, which is incorporated herein by reference in its entirety. As described herein, the experiment included four (4) unique samples. Which cells belonged to which sample was essentially tracked by observing the first barcode (incorporated by RT). This sample identification may require that the barcode sequence corresponding to each well and which sample was placed in which well are known. With this information, a lookup table from the first round barcode sequence to sample ID can be created.

[0199] Figures 20 and 21 show the experimental setup for multiplexing samples. In this particular case, four samples were multiplexed across 48 wells of a 96-well plate. Extending this same theory, it can be seen that up to 96 samples can be multiplexed in a single experiment.

[0200] Referring to Figure 22, the analysis of transcriptomes using T-distributed stochastic neighbor embedding (t-SNE) in an experiment with four different biological samples is shown. Each point represents the transcriptome from one individual cell. The identity of the barcode incorporated in RT (using barcoded RT primers) can be used to determine the sample identity of each transcriptome.

[0201] Example 16 - PEG incorporation for template switching of DNA / RNA bound to magnetic beads Biotinylated oligos attached to the reverse transcribed RNA (cDNA / RNA duplex) are bound to streptavidin-coated magnetic beads. Template switching is then performed on this cDNA / RNA duplex so that a single common adapter sequence can be incorporated at the 3' end of the cDNA. The efficiency of this step, as measured by transcript detection per cell, can be improved by incorporating up to 10 w / v% PEG (molecular weight 7000-9000) into the template switching reaction that contains the cDNA / RNA duplex bound to streptavidin-coated magnetic beads (Table 22).

[0202] [Table 22]

[0203] Example 18 - Exemplary methods for uniquely labeling RNA molecules Panel A of FIG. 23 shows the labeling of the transcriptome by distribution-pool barcoding. In each distribution-pool round, fixed cells or nuclei can be randomly distributed into wells and transcripts can be labeled with well-specific barcodes. Barcoded RT primers can be used for the first round. Second and third round barcodes can be added to the cDNA by ligation. A fourth barcode can be added to the cDNA molecules by PCR during sequencing library preparation. The bottom scheme shows the final barcoded cDNA molecules.

[0204] Panel B of Figure 23 shows a species mixing experiment with a library prepared from 1,758 whole cells. Human UBCs extend horizontally along the x-axis, mouse UBCs extend vertically along the y-axis, and the mixed species UBCs are distributed between the human and mouse UBCs. The estimated barcode collision frequency is 0.2%, while the species purity is >99%.

[0205] Panel C of Figure 23 shows UMI counts from mixing experiments performed on fresh frozen (stored at -80°C for 2 weeks) cells and nuclei. Median human UMI counts in fresh cells: 15,365; frozen cells: 15,078; nuclei: 12,113; frozen nuclei: 13,636.

[0206] Panel D of Figure 23 shows that gene expression measured by the methods provided herein is highly correlated between frozen and immediately processed cells (Pearson-r: 0.987). Freshly frozen cells were processed in two different SPLiT-seq experiments.

[0207] As illustrated in FIG. 24, fixed and permeabilized cells can be randomly distributed into wells, each of which contains a reverse transcription primer with a well-specific barcode. In situ reverse transcription converts the RNA to cDNA while adding a well-specific barcode. The cells can then be pooled and again randomly distributed into a second set of wells, each of which contains a unique well-specific barcode. These barcodes can be hybridized and ligated to the 5' end of a barcoded reverse transcription primer to generate a second lambda. The barcodes can be added to the barcoding of the strands. The cells can be pooled together again and a subsequent partition-ligation-pooling round can be performed. After the last round of ligation, the cDNA molecules contain a cell-specific combination of barcode, unique molecular identifier, and universal PCR handle at the 5' end. During the PCR step of library preparation, a fourth barcoding round can be performed.

[0208] The oligonucleotides listed in Table 23 were used in the above examples. "rG" is an RNA base and "+G" is a locked nucleobase.

[0209] [Table 23-1]

[0210] [Table 23-2]

[0211] Certain embodiments of the present disclosure are described herein, including the best mode known to the inventors for carrying out the disclosure. Of course, variations on these described embodiments will become apparent to those of ordinary skill in the art upon reading the above description. Applicants It is expected that those skilled in the art will use such variations as necessary, and the applicants intend that various embodiments of the present disclosure may be implemented differently from the specific description herein.Accordingly, the present disclosure includes all modifications and equivalents of the subject matter described in the claims appended hereto as permitted by applicable law.Furthermore, any combination of the above-described elements of any and all of these variations are encompassed by the present disclosure unless otherwise specified or otherwise clearly contradicted by the circumstances.

[0212] Additionally, throughout this specification, numerous references are made to patents and printed publications. Each of the above cited references and printed publications is individually incorporated herein by reference in its entirety.

[0213] It should be understood that the embodiments of the present disclosure are illustrative of the principles of the present disclosure. Other modifications that can be used are within the scope of the present disclosure. Thus, by way of example and not of limitation, alternative configurations of the present disclosure can be used in accordance with the teachings herein. Thus, the present disclosure is not limited to exactly that shown and described.

[0214] The details set forth herein are examples and are merely for illustrative purposes of the preferred embodiments of the present disclosure, presented in order to provide what is believed to be the most useful and readily understood explanation of the principles and conceptual aspects of various embodiments of the present disclosure.

[0215] It will be apparent to those skilled in the art that numerous changes may be made to the details of the above-described embodiments without departing from the underlying principles of the disclosure, and the scope of the invention should therefore be determined solely by the claims which follow.

Claims

1. 1. A method for uniquely labeling complementary DNA (cDNA) molecules in a plurality of cells, comprising: (a) fixing and permeabilizing a first plurality of cells prior to step (b), wherein the first plurality of cells is fixed and permeabilized at 4° C. or below; (b) reverse transcribing RNA molecules in the first plurality of cells to form cDNA molecules in the first plurality of cells, wherein reverse transcribing the RNA molecules comprises coupling a primer to the RNA molecules, the primer comprising at least one of a poly(T) sequence or a random sequence; (c) dividing the first plurality of cells comprising cDNA molecules into at least two primary aliquots, the at least two primary aliquots comprising a first primary aliquot and a second primary aliquot; (d) providing primary nucleic acid tags to the at least two primary aliquots, wherein the primary nucleic acid tag provided to the first primary aliquot is different from the primary nucleic acid tag provided to the second primary aliquot; (e) coupling said cDNA molecules in each of said at least two primary aliquots with said provided primary nucleic acid tags; (f) combining said at least two primary aliquots; (g) dividing the combined primary aliquot into at least two secondary aliquots, the at least two secondary aliquots comprising a first secondary aliquot and a second secondary aliquot; (h) providing secondary nucleic acid tags to the at least two secondary aliquots, wherein the secondary nucleic acid tag provided to the first secondary aliquot is different from the secondary nucleic acid tag provided to the second secondary aliquot; (i) coupling said cDNA molecules in each of said at least two secondary aliquots with said provided secondary nucleic acid tags; (j) combining the final aliquots; (k) lysing the first plurality of cells to release the cDNA molecules from within the first plurality of cells to form a lysate; and (l) adding a protease inhibitor and a binding agent to the lysate such that the cDNA molecules bind to the binding agent; The method includes:

2. 2. The method of claim 1, further comprising, prior to step (k), dividing the combined final aliquot into at least two final aliquots, wherein the at least two final aliquots include a first final aliquot and a second final aliquot.

3. 3. The method of claim 1 or 2, wherein the protease inhibitor comprises phenylmethanesulfonyl fluoride (PMSF) or 4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride (AEBSF).

4. (m) performing template switching of the cDNA molecule bound to the binding agent; (n) amplifying the cDNA molecules to form a solution of amplified cDNA molecules; and (o) introducing a solution of solid phase reversible immobilization (SPRI) beads into said solution of amplified cDNA molecules, wherein the ratio of SPRI bead solution to amplified cDNA molecule solution is between 0.9:1 and 0.7:1; The method of any one of claims 1 to 3, further comprising:

5. The method of any one of claims 1 to 4, wherein the common adaptor sequence is added to the 3' end of the released cDNA molecule by performing template switching.

6. The primers of step (b) further comprise a first specific barcode, and the method further comprises: (p) reverse transcribing RNA molecules in a second plurality of cells to form cDNA molecules in the second plurality of cells, wherein reverse transcribing the RNA molecules comprises coupling a specific primer to the RNA molecules, the primer comprising a second specific barcode and at least one of a poly(T) sequence or a random sequence, the first specific barcode being different from the second specific barcode such that the cDNA molecules derived from the first plurality of cells can be distinguished relative to the cDNA molecules derived from the second plurality of cells; (q) dividing the second plurality of cells comprising cDNA molecules into at least two primary aliquots, the at least two primary aliquots comprising a first primary aliquot and a second primary aliquot; (r) providing primary nucleic acid tags to the at least two primary aliquots, wherein the primary nucleic acid tag provided to the first primary aliquot is different from the primary nucleic acid tag provided to the second primary aliquot; (s) coupling said cDNA molecules in each of said at least two primary aliquots to said provided primary nucleic acid tags; (t) combining said at least two primary aliquots; (u) dividing the combined primary aliquot into at least two secondary aliquots, the at least two secondary aliquots comprising a first secondary aliquot and a second secondary aliquot; (v) providing secondary nucleic acid tags to the at least two secondary aliquots, the secondary nucleic acid tag provided to the first secondary aliquot being different from the secondary nucleic acid tag provided to the second secondary aliquot; and (w) coupling said cDNA molecules in each of said at least two secondary aliquots with said provided secondary nucleic acid tags. The method of any one of claims 1 to 5, further comprising:

7. Each of the nucleic acid tags comprises: a barcode sequence comprising a 3' end and a 5' end; and a 3' hybridization sequence and a 5' hybridization sequence adjacent to the 3' end and the 5' end, respectively, of the barcode sequence; a first strand comprising: a first portion complementary to at least one of the 5' hybridization sequence and the adapter sequence; and a second portion complementary to the 3' hybridization sequence A second strand comprising The method according to any one of claims 1 to 6, comprising:

8. Ligating at least two of said nucleic acid tags attached to said cDNA molecules. The method of any one of claims 1 to 7, further comprising:

9. The method of claim 8 , wherein the ligation is performed within the first plurality of cells.

10. Ligating at least two of the nucleic acid tags attached to the released cDNA molecules. The method of any one of claims 1 to 7, further comprising:

Images