Normalization of nucleic acid libraries

Abstract

This disclosure provides methods and compositions for removing one or more high abundance species from a plurality of nucleic acid molecules. In some embodiments, the methods and compositions can be used for normalizing nucleic acid libraries. In some embodiments, molecular labels are used in conjunction with the methods and compositions disclosed herein to improve sequencing efficiency.

Description

BACKGROUND

The level of gene expression in a biological sample can vary greatly. For examples, it has been described that gene expression level follows 3 broad categories: 1) ‘high expressers,’ which are comprised of 5-10 genes that dominate ˜20% of cellular mRNAs; 2) ‘intermediate expressers’ that are comprised of 50-200 genes that occupy 40-60% of cellular mRNAs; and 3) ‘moderate expressers’ that are comprised of 10,000-20,000 genes that occupy the rest of the cellular mRNA fraction. One challenge in molecular biology and molecular genetics is to be able to capture this highly dynamic gene expression profile efficiently and effectively in order to distinguish different cell types and phenotypes in the sample.

In recent years, next generation sequencing (NGS) has provided a high throughput method in assessing gene expression profiles. During library preparation for NGS, a sample with heterogeneous cDNA species is amplified by PCR to obtain adequate sample amount and to attach NGS-compatible adapters. The sequencing process captures the number of reads for each gene from the PCR-amplified library sample to interpret the gene expression level. However, since different genes are expressed at a large range of levels, PCR amplification can skew the native gene expression. For example, a gene has 1 molecule of cDNA would require 40 cycles of PCR to achieve the same representative amount as a gene with 1000 molecules of cDNA in 30 cycles. In a heterogeneous cDNA sample, PCR is usually performed in excess cycles to adequately amplify low expressers; in those scenarios, the native gene expression profile is usually skewed by the dominating high expresser PCR products. A method to correct for such bias in PCR product is Molecular Indexing; however, high expressers such as ribosomal protein mRNAs, mitochondrial mRNAs, or housekeeping genes often dominate the sequencing run with little contribution to the experimental interpretation, rendering the sequencing cost for Molecular Index counting to be expensive.

SUMMARY

Some embodiments disclosed herein provide methods of removing high abundance species from a plurality of nucleic acid molecules, comprising: hybridizing a plurality of first oligonucleotides comprising an affinity moiety with a first plurality of nucleic acid molecules, wherein the first plurality of nucleic acid molecules comprises at least one high abundance species; extending the plurality of oligonucleotides to generate a plurality of complementary strands of the first plurality of nucleic acid molecules comprising the affinity moiety; denaturing a plurality of double-stranded nucleic acid molecules comprising the plurality of complementary strands of the first plurality of nucleic acid molecules; partially reannealing the plurality of complementary strands of the first plurality of nucleic acid molecules; and removing the reannealed complementary strands of the first plurality of nucleic acid molecules by a capture molecule immobilized on one or more solid support to generate a second plurality of nucleic acid molecules, wherein the capture molecules specifically bind to the affinity moiety, whereby the at least one high abundance species is reduced in the second plurality of nucleic acid molecules. In some embodiments, the affinity moiety is a functional group selected from the group consisting of biotin, streptavidin, heparin, an aptamer, a click-chemistry moiety, digoxigenin, primary amine(s), carboxyl(s), hydroxyl(s), aldehyde(s), ketone(s), and any combination thereof. In some embodiments, the affinity moiety is biotin. In some embodiments, the capture molecule is streptavidin. In some embodiments, the methods further comprise synthesizing a second strand for each of the plurality of complementary strands of the first plurality of nucleic acid molecules to generate the plurality of double-stranded nucleic acid molecules comprising the plurality of complementary strands of the first plurality of nucleic acid molecules. In some embodiments, the synthesizing comprises hybridizing a plurality of second oligonucleotides to the plurality of complementary strands of the first plurality of nucleic acid molecules and extending the plurality of second oligonucleotide. In some embodiments, the plurality of first oligonucleotides or the plurality of second oligonucleotides comprises a universal primer binding site. In some embodiments, the methods further comprise amplifying the plurality of double-stranded nucleic acid molecules. In some embodiments, the first plurality of nucleic acid molecules comprises a plurality of high abundance species. In some embodiments, the plurality of high abundance species represents at least 50% of the first plurality of nucleic acid molecules. In some embodiments, the plurality of high abundance species represents at least 60% of the first plurality of nucleic acid molecules. In some embodiments, the plurality of high abundance species represents at least 70% of the first plurality of nucleic acid molecules. In some embodiments, at least 80% of the plurality of high abundance species is removed. In some embodiments, at least 90% of the plurality of high abundance species is removed. In some embodiments, at least 95% of the plurality of high abundance species is removed. In some embodiments, at least 99% of the plurality of high abundance species is removed. In some embodiments, the second plurality of nucleic acid molecules comprises the plurality of high abundance species. In some embodiments, the plurality of high abundance species in the second plurality of nucleic acid molecules represents less than 50% of the second plurality of nucleic acid molecules. In some embodiments, the plurality of high abundance species in the second plurality of nucleic acid molecules represents less than 40% of the second plurality of nucleic acid molecules. In some embodiments, the plurality of high abundance species in the second plurality of nucleic acid molecules represents less than 30% of the second plurality of nucleic acid molecules. In some embodiments, the first plurality of nucleic acid molecules comprises a plurality of low abundance species. In some embodiments, the plurality of low abundance species represents less than 10% of the first plurality of nucleic acid molecules. In some embodiments, the plurality of low abundance species represents less than 5% of the first plurality of nucleic acid molecules. In some embodiments, the plurality of low abundance species represents less than 1% of the first plurality of nucleic acid molecules. In some embodiments, the second plurality of nucleic acid molecules comprises the plurality of low abundance species. In some embodiments, the plurality of low abundance species in the second plurality of nucleic acid molecules represents at least 5% of the second plurality of nucleic acid molecules. In some embodiments, the plurality of low abundance species in the second plurality of nucleic acid molecules represents at least 10% of the second plurality of nucleic acid molecules. In some embodiments, the plurality of low abundance species in the second plurality of nucleic acid molecules represents at least 20% of the second plurality of nucleic acid molecules. In some embodiments, each of the first plurality of nucleic acid molecules or each of the second plurality of nucleic acid molecules comprises a stochastic barcode. In some embodiments, the methods further comprise sequencing the second plurality of nucleic acid molecules to generate a plurality of sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance species is less than 50% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance species is less than 40% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance species is less than 30% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance species is at least 5% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance species is at least 10% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance species is at least 20% of the total sequencing reads. In some embodiments, each of the plurality of sequencing reads comprises a molecular label. In some embodiments, the ratio of sequencing reads to molecular label for the plurality of high abundance species is less than 15. In some embodiments, the ratio of sequencing reads to molecular label for the plurality of high abundance species is less than 10. In some embodiments, the methods further comprise selectively depleting a high abundance species. In some embodiments, the selectively depleting comprises hybridizing a target-specific oligonucleotide that specifically binds to the high abundance species. In some embodiments, the selectively depleting comprises treating the first plurality of nucleic acid molecules with a Cas9 protein complexed with a guide oligonucleotide that specifically binds to the high abundance species.

Some embodiments disclosed herein provide methods of generating a normalized nucleic acid library, comprising: hybridizing a plurality of first oligonucleotides comprising an affinity moiety with a plurality of nucleic acid targets in a sample; extending the plurality of oligonucleotides to generate a plurality of complementary strands of the plurality of nucleic acid targets comprising the affinity moiety; denaturing a plurality of double-stranded nucleic acid molecules comprising the plurality of complementary strands of the plurality of nucleic acid targets; partially reannealing the plurality of complementary strands of the plurality of nucleic acid targets; and removing the reannealed complementary strands of the plurality of nucleic acid targets by a capture molecule immobilized on one or more solid support, wherein the capture molecules specifically bind to the affinity moiety, whereby a normalized nucleic acid library of the plurality of nucleic acid targets is generated. In some embodiments, the affinity moiety is a functional group selected from the group consisting of biotin, streptavidin, heparin, an aptamer, a click-chemistry moiety, digoxigenin, primary amine(s), carboxyl(s), hydroxyl(s), aldehyde(s), ketone(s), and any combination thereof. In some embodiments, the affinity moiety is a functional group selected from the group consisting of biotin, streptavidin, heparin, an aptamer, a click-chemistry moiety, digoxigenin, primary amine(s), carboxyl(s), hydroxyl(s), aldehyde(s), ketone(s), and any combination thereof. In some embodiments, the affinity moiety is biotin. In some embodiments, the capture molecule is streptavidin. In some embodiments, the methods further comprise synthesizing a second strand for each of the plurality of complementary strands of the plurality of nucleic acid targets to generate the plurality of double-stranded nucleic acid molecules comprising the plurality of complementary strands of the plurality of nucleic acid targets. In some embodiments, the synthesizing comprises hybridizing a plurality of second oligonucleotides to the plurality of complementary strands of the plurality of nucleic acid targets and extending the plurality of second oligonucleotide. In some embodiments, the plurality of first oligonucleotides or the plurality of second oligonucleotides comprises a universal primer binding site. In some embodiments, the methods further comprise amplifying the plurality of double-stranded nucleic acid molecules. In some embodiments, the plurality of nucleic acid targets comprises a plurality of low abundance nucleic acid targets. In some embodiments, the plurality of low abundance nucleic acid targets represents less than 10% of the plurality of nucleic acid targets. In some embodiments, the plurality of low abundance nucleic acid targets represents less than 5% of the plurality of nucleic acid targets. In some embodiments, the plurality of low abundance nucleic acid targets represents less than 1% of the plurality of nucleic acid targets. In some embodiments, the normalized nucleic acid library of the plurality of nucleic acid targets comprises the plurality of low abundance nucleic acid targets. In some embodiments, the plurality of low abundance nucleic acid targets in the normalized nucleic acid library represents at least 5% of the plurality of nucleic acid targets in the normalized nucleic acid library. In some embodiments, the plurality of low abundance nucleic acid targets in the normalized nucleic acid library represents at least 10% of the plurality of nucleic acid targets in the normalized nucleic acid library. In some embodiments, the plurality of low abundance nucleic acid targets in the normalized nucleic acid library represents at least 20% of the plurality of nucleic acid targets in the normalized nucleic acid library. In some embodiments, the plurality of nucleic acid targets comprises a plurality of high abundance nucleic acid targets. In some embodiments, the plurality of high abundance nucleic acid targets represents at least 50% of the plurality of nucleic acid target. In some embodiments, the plurality of high abundance nucleic acid targets represents at least 60% of the plurality of nucleic acid target. In some embodiments, the plurality of high abundance nucleic acid targets represents at least 70% of the plurality of nucleic acid target. In some embodiments, at least 80% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, at least 90% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, at least 95% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, at least 99% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, the normalized nucleic acid library of the plurality of nucleic acid targets comprises the plurality of high abundance nucleic acid targets. In some embodiments, the plurality of high abundance nucleic acid targets in the normalized nucleic acid library represents less than 50% of the plurality of nucleic acid targets in the normalized nucleic acid library. In some embodiments, the plurality of high abundance nucleic acid targets in the normalized nucleic acid library represents less than 40% of the plurality of nucleic acid targets in the normalized nucleic acid library. In some embodiments, the plurality of high abundance nucleic acid targets in the normalized nucleic acid library represents less than 30% of the plurality of nucleic acid targets in the normalized nucleic acid library. In some embodiments, each of the plurality of first oligonucleotides or each of the plurality of second oligonucleotides comprises a stochastic barcode. In some embodiments, the methods further comprise sequencing the normalized nucleic acid library to generate a plurality of sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 50% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 40% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 30% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 5% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 10% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 20% of the total sequencing reads. In some embodiments, each of the plurality of sequencing reads comprises a molecular label. In some embodiments, the ratio of sequencing reads to molecular label for the plurality of high abundance nucleic acid targets is less than 15. In some embodiments, the ratio of sequencing reads to molecular label for the plurality of high abundance nucleic acid targets is less than 10. In some embodiments, the plurality of nucleic acid targets comprises mRNA. In some embodiments, the plurality of nucleic acid targets comprises mitochondrial mRNA. In some embodiments, the plurality of nucleic acid targets comprises ribosomal protein mRNA. In some embodiments, the low abundance nucleic acid targets comprise 7,000 genes with the lowest number of transcripts. In some embodiments, the low abundance nucleic acid targets comprise 4,000 genes with the lowest number of transcripts. In some embodiments, the low abundance nucleic acid targets comprise 2,000 genes with the lowest number of transcripts. In some embodiments, the plurality of first oligonucleotides comprise target-specific primers. In some embodiments, the plurality of first oligonucleotides comprise non-target-specific primers. In some embodiments, the plurality of nucleic acid targets comprises cDNA. In some embodiments, the plurality of nucleic acid targets comprises genomic DNA. In some embodiments, the high abundance nucleic acid targets comprise short tandem repeat sequences. In some embodiments, the high abundance nucleic acid targets comprise telomeric sequences. In some embodiments, the high abundance nucleic acid targets comprise centrometic sequences. In some embodiments, the sample is a single cell. In some embodiments, the methods further comprise selectively depleting a high abundance species. In some embodiments, the selectively depleting comprises hybridizing a target-specific oligonucleotide that specifically binds to the high abundance species. In some embodiments, the selectively depleting comprises treating the first plurality of nucleic acid targets with a Cas9 protein complexed with a guide oligonucleotide that specifically binds to the high abundance species.

Some embodiments disclosed herein provide methods of generating a normalized nucleic acid library, comprising: hybridizing a plurality of first oligonucleotides comprising an affinity moiety with a plurality of nucleic acid targets in an unnormalized nucleic acid library; extending the plurality of oligonucleotides to generate a plurality of complementary strands of the plurality of nucleic acid targets comprising the affinity moiety; denaturing a plurality of double-stranded nucleic acid molecules comprising the plurality of complementary strands of the plurality of nucleic acid targets; partially reannealing the plurality of complementary strands of the plurality of nucleic acid targets; and removing the reannealed complementary strands of the plurality of nucleic acid targets, whereby a normalized nucleic acid library of the plurality of nucleic acid targets is generated. In some embodiments, the unnormalized nucleic acid library comprises one or more high abundance nucleic acid targets and one or more low abundance nucleic acid target. In some embodiments, the one or more high abundance nucleic acid targets represents at least 50% of the unnormalized nucleic acid library. In some embodiments, the one or more high abundance nucleic acid targets represents at least 60% of the unnormalized nucleic acid library. In some embodiments, the one or more high abundance nucleic acid targets represents at least 70% of the unnormalized nucleic acid library. In some embodiments, at least 80% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, at least 90% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, at least 95% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, at least 99% of the plurality of high abundance nucleic acid targets is removed. In some embodiments, the one or more low abundance nucleic acid targets represents less than 10% of the unnormalized nucleic acid library. In some embodiments, the one or more low abundance nucleic acid targets represents less than 5% of the unnormalized nucleic acid library. In some embodiments, the one or more low abundance nucleic acid targets represents less than 1% of the unnormalized nucleic acid library. In some embodiments, the one or more low abundance nucleic acid targets represents at least 5% of the normalized nucleic acid library. In some embodiments, the one or more low abundance nucleic acid targets represents at least 10% of the normalized nucleic acid library. In some embodiments, the one or more low abundance nucleic acid targets represents at least 20% of the normalized nucleic acid library. In some embodiments, the one or more high abundance nucleic acid targets represents less than 50% of the normalized nucleic acid library. In some embodiments, the one or more high abundance nucleic acid targets represents less than 40% of the normalized nucleic acid library. In some embodiments, the one or more high abundance nucleic acid targets represents less than 30% of the normalized nucleic acid library. In some embodiments, the unnormalized nucleic acid library is a cDNA library. In some embodiments, the unnormalized nucleic acid library is a genomic library. In some embodiments, the unnormalized nucleic acid library is a single-cell library. In some embodiments, the methods further comprise selectively depleting a high abundance nucleic acid target. In some embodiments, the selectively depleting comprises hybridizing a target-specific oligonucleotide that specifically binds to the high abundance nucleic acid target. In some embodiments, the selectively depleting comprises treating the first plurality of nucleic acid targets with a Cas9 protein complexed with a guide oligonucleotide that specifically binds to the high abundance nucleic acid target.

BRIEF DESCRIPTION OF THE DRAWINGS

FIGS. 1A and 1B show a schematic illustration of exemplary library normalization by biotin immobilization.

FIG. 2 shows simulated normalization results on sequencing reads using sequencing results from expression data of targeted TCR genes.

FIG. 3 shows simulated normalization results on average sequencing reads per MI using sequencing results from expression data of targeted TCR genes.

FIG. 4 shows simulated normalization results on sequencing reads using sequencing results from whole transcriptome amplification data.

FIG. 5 shows simulated normalization results on sequencing read distribution using sequencing results from whole transcriptome amplification data.

FIG. 6 shows simulated normalization results on sequencing reads per MI using sequencing results from whole transcriptome amplification data.

FIG. 7 shows simulated normalization results on sequencing reads using shallow sequencing results from single T cell whole transcriptome amplification data.

FIG. 8 shows simulated normalization results on sequencing read distribution using shallow sequencing results from single T cell whole transcriptome amplification data.

FIG. 9 shows simulated normalization results on sequencing reads per MI using shallow sequencing results from single T cell from whole transcriptome amplification data.

FIG. 10 shows simulated normalization results on sequencing reads using deep sequencing results from single tumor cell whole transcriptome amplification data.

FIG. 11 shows simulated normalization results on sequencing read distribution using deep sequencing results from single tumor cell whole transcriptome amplification data.

FIG. 12 shows simulated normalization results on sequencing reads per MI using deep sequencing results from single tumor cell from whole transcriptome amplification data.

DETAILED DESCRIPTION

Definitions

Unless otherwise defined, all technical terms used herein have the same meaning as commonly understood by one of ordinary skill in the art in the field to which this disclosure belongs. As used in this specification and the appended claims, the singular forms “a,” “an,” and “the” include plural references unless the context clearly dictates otherwise. Any reference to “or” herein is intended to encompass “and/or” unless otherwise stated.

As used herein the term “associated” or “associated with” can mean that two or more species are identifiable as being co-located at a point in time. An association can mean that two or more species are or were within a similar container. An association can be an informatics association, where for example digital information regarding two or more species is stored and can be used to determine that one or more of the species were co-located at a point in time. An association can also be a physical association. In some instances two or more associated species are “tethered”, “attached”, or “immobilized” to one another or to a common solid or semisolid surface. An association may refer to covalent or non-covalent means for attaching labels to solid or semi-solid supports such as beads. An association may comprise hybridization between a target and a label.

As used herein, the term “complementary” can refer to the capacity for precise pairing between two nucleotides. For example, if a nucleotide at a given position of a nucleic acid is capable of hydrogen bonding with a nucleotide of another nucleic acid, then the two nucleic acids are considered to be complementary to one another at that position. Complementarity between two single-stranded nucleic acid molecules may be “partial,” in which only some of the nucleotides bind, or it may be complete when total complementarity exists between the single-stranded molecules. A first nucleotide sequence can be said to be the “complement” of a second sequence if the first nucleotide sequence is complementary to the second nucleotide sequence. A first nucleotide sequence can be said to be the “reverse complement” of a second sequence, if the first nucleotide sequence is complementary to a sequence that is the reverse (i.e., the order of the nucleotides is reversed) of the second sequence. As used herein, the terms “complement”, “complementary”, and “reverse complement” can be used interchangeably. It is understood from the disclosure that if a molecule can hybridize to another molecule it may be the complement of the molecule that is hybridizing.

As used herein, the term “digital counting” can refer to a method for estimating a number of target molecules in a sample. Digital counting can include the step of determining a number of unique labels that have been associated with targets in a sample. This stochastic methodology transforms the problem of counting molecules from one of locating and identifying identical molecules to a series of yes/no digital questions regarding detection of a set of predefined labels.

As used herein, the term “label” or “labels” can refer to nucleic acid codes associated with a target within a sample. A label can be, for example, a nucleic acid label. A label can be an entirely or partially amplifiable label. A label can be entirely or partially sequencable label. A label can be a portion of a native nucleic acid that is identifiable as distinct. A label can be a known sequence. A label can comprise a junction of nucleic acid sequences, for example a junction of a native and non-native sequence. As used herein, the term “label” can be used interchangeably with the terms, “index”, “tag,” or “label-tag.” Labels can convey information. For example, in various embodiments, labels can be used to determine an identity of a sample, a source of a sample, an identity of a cell, and/or a target.

As used herein, a “nucleic acid” can generally refer to a polynucleotide sequence, or fragment thereof. A nucleic acid can comprise nucleotides. A nucleic acid can be exogenous or endogenous to a cell. A nucleic acid can exist in a cell-free environment. A nucleic acid can be a gene or fragment thereof. A nucleic acid can be DNA. A nucleic acid can be RNA. A nucleic acid can comprise one or more analogs (e.g. altered backgone, sugar, or nucleobase). Some non-limiting examples of analogs include: 5-bromouracil, peptide nucleic acid, xeno nucleic acid, morpholinos, locked nucleic acids, glycol nucleic acids, threose nucleic acids, dideoxynucleotides, cordycepin, 7-deaza-GTP, florophores (e.g. rhodamine or flurescein linked to the sugar), thiol containing nucleotides, biotin linked nucleotides, fluorescent base analogs, CpG islands, methyl-7-guanosine, methylated nucleotides, inosine, thiouridine, pseudourdine, dihydrouridine, queuosine, and wyosine. “Nucleic acid”, “polynucleotide, “target polynucleotide”, and “target nucleic acid” can be used interchangeably.

A nucleic acid can comprise one or more modifications (e.g., a base modification, a backbone modification), to provide the nucleic acid with a new or enhanced feature (e.g., improved stability). A nucleic acid can comprise a nucleic acid affinity tag. A nucleoside can be a base-sugar combination. The base portion of the nucleoside can be a heterocyclic base. The two most common classes of such heterocyclic bases are the purines and the pyrimidines. Nucleotides can be nucleosides that further include a phosphate group covalently linked to the sugar portion of the nucleoside. For those nucleosides that include a pentofuranosyl sugar, the phosphate group can be linked to the 2′, the 3′, or the 5′ hydroxyl moiety of the sugar. In forming nucleic acids, the phosphate groups can covalently link adjacent nucleosides to one another to form a linear polymeric compound. In turn, the respective ends of this linear polymeric compound can be further joined to form a circular compound; however, linear compounds are generally suitable. In addition, linear compounds may have internal nucleotide base complementarity and may therefore fold in a manner as to produce a fully or partially double-stranded compound. Within nucleic acids, the phosphate groups can commonly be referred to as forming the internucleoside backbone of the nucleic acid. The linkage or backbone of the nucleic acid can be a 3′ to 5′ phosphodiester linkage.

A nucleic acid can comprise a modified backbone and/or modified internucleoside linkages. Modified backbones can include those that retain a phosphorus atom in the backbone and those that do not have a phosphorus atom in the backbone. Suitable modified nucleic acid backbones containing a phosphorus atom therein can include, for example, phosphorothioates, chiral phosphorothioates, phosphorodithioates, phosphotriesters, aminoalkylphosphotriesters, methyl and other alkyl phosphonates such as 3′-alkylene phosphonates, 5′-alkylene phosphonates, chiral phosphonates, phosphinates, phosphoramidates including 3′-amino phosphoramidate and aminoalkylphosphoramidates, phosphorodiamidates, thionophosphoramidates, thionoalkylphosphonates, thionoalkylphosphotriesters, selenophosphates, and boranophosphates having normal 3′-5′ linkages, 2′-5′ linked analogs, and those having inverted polarity wherein one or more internucleotide linkages is a 3′ to 3′, a 5′ to 5′ or a 2′ to 2′ linkage.

A nucleic acid can comprise polynucleotide backbones that are formed by short chain alkyl or cycloalkyl internucleoside linkages, mixed heteroatom and alkyl or cycloalkyl internucleoside linkages, or one or more short chain heteroatomic or heterocyclic internucleoside linkages. These can include those having morpholino linkages (formed in part from the sugar portion of a nucleoside); siloxane backbones; sulfide, sulfoxide and sulfone backbones; formacetyl and thioformacetyl backbones; methylene formacetyl and thioformacetyl backbones; riboacetyl backbones; alkene containing backbones; sulfamate backbones; methyleneimino and methylenehydrazino backbones; sulfonate and sulfonamide backbones; amide backbones; and others having mixed N, O, S and CH2 component parts.

A nucleic acid can comprise a nucleic acid mimetic. The term “mimetic” can be intended to include polynucleotides wherein only the furanose ring or both the furanose ring and the internucleotide linkage are replaced with non-furanose groups, replacement of only the furanose ring can also be referred as being a sugar surrogate. The heterocyclic base moiety or a modified heterocyclic base moiety can be maintained for hybridization with an appropriate target nucleic acid. One such nucleic acid can be a peptide nucleic acid (PNA). In a PNA, the sugar-backbone of a polynucleotide can be replaced with an amide containing backbone, in particular an aminoethylglycine backbone. The nucleotides can be retained and are bound directly or indirectly to aza nitrogen atoms of the amide portion of the backbone. The backbone in PNA compounds can comprise two or more linked aminoethylglycine units which gives PNA an amide containing backbone. The heterocyclic base moieties can be bound directly or indirectly to aza nitrogen atoms of the amide portion of the backbone.

A nucleic acid can comprise a morpholino backbone structure. For example, a nucleic acid can comprise a 6-membered morpholino ring in place of a ribose ring. In some of these embodiments, a phosphorodiamidate or other non-phosphodiester internucleoside linkage can replace a phosphodiester linkage.

A nucleic acid can comprise linked morpholino units (i.e. morpholino nucleic acid) having heterocyclic bases attached to the morpholino ring. Linking groups can link the morpholino monomeric units in a morpholino nucleic acid. Non-ionic morpholino-based oligomeric compounds can have less undesired interactions with cellular proteins. Morpholino-based polynucleotides can be nonionic mimics of nucleic acids. A variety of compounds within the morpholino class can be joined using different linking groups. A further class of polynucleotide mimetic can be referred to as cyclohexenyl nucleic acids (CeNA). The furanose ring normally present in a nucleic acid molecule can be replaced with a cyclohexenyl ring. CeNA DMT protected phosphoramidite monomers can be prepared and used for oligomeric compound synthesis using phosphoramidite chemistry. The incorporation of CeNA monomers into a nucleic acid chain can increase the stability of a DNA/RNA hybrid. CeNA oligoadenylates can form complexes with nucleic acid complements with similar stability to the native complexes. A further modification can include Locked Nucleic Acids (LNAs) in which the 2′-hydroxyl group is linked to the 4′ carbon atom of the sugar ring thereby forming a 2′-C,4′-C-oxymethylene linkage thereby forming a bicyclic sugar moiety. The linkage can be a methylene (—CH2-), group bridging the 2′ oxygen atom and the 4′ carbon atom wherein n is 1 or 2. LNA and LNA analogs can display very high duplex thermal stabilities with complementary nucleic acid (Tm=+3 to +10° C.), stability towards 3′-exonucleolytic degradation and good solubility properties.

A nucleic acid may also include nucleobase (often referred to simply as “base”) modifications or substitutions. As used herein, “unmodified” or “natural” nucleobases can include the purine bases, (e.g. adenine (A) and guanine (G)), and the pyrimidine bases, (e.g. thymine (T), cytosine (C) and uracil (U)). Modified nucleobases can include other synthetic and natural nucleobases such as 5-methylcytosine (5-me-C), 5-hydroxymethyl cytosine, xanthine, hypoxanthine, 2-aminoadenine, 6-methyl and other alkyl derivatives of adenine and guanine, 2-propyl and other alkyl derivatives of adenine and guanine, 2-thiouracil, 2-thiothymine and 2-thiocytosine, 5-halouracil and cytosine, 5-propynyl (—C═C—CH3) uracil and cytosine and other alkynyl derivatives of pyrimidine bases, 6-azo uracil, cytosine and thymine, 5-uracil (pseudouracil), 4-thiouracil, 8-halo, 8-amino, 8-thiol, 8-thioalkyl, 8-hydroxyl and other 8-substituted adenines and guanines, 5-halo particularly 5-bromo, 5-trifluoromethyl and other 5-substituted uracils and cytosines, 7-methylguanine and 7-methyladenine, 2-F-adenine, 2-aminoadenine, 8-azaguanine and 8-azaadenine, 7-deazaguanine and 7-deazaadenine and 3-deazaguanine and 3-deazaadenine. Modified nucleobases can include tricyclic pyrimidines such as phenoxazine cytidine(1H-pyrimido(5,4-b)(1,4)benzoxazin-2(3H)-one), phenothiazine cytidine (1H-pyrimido(5,4-b)(1,4)benzothiazin-2(3H)-one), G-clamps such as a substituted phenoxazine cytidine (e.g. 9-(2-aminoethoxy)-H-pyrimido(5,4-(b) (1,4)benzoxazin-2(3H)-one), carbazole cytidine (2H-pyrimido(4,5-b)indol-2-one), pyridoindole cytidine (Hpyrido(3′,′:4,5)pyrrolo[2,3-d]pyrimidin-2-one).

As used herein, the term “sample” can refer to a composition comprising targets. Suitable samples for analysis by the disclosed methods, devices, and systems include cells, single cells, tissues, organs, or organisms.

As used herein, the term “sampling device” or “device” can refer to a device which may take a section of a sample and/or place the section on a substrate. A sample device can refer to, for example, a fluorescence activated cell sorting (FACS) machine, a cell sorter machine, a biopsy needle, a biopsy device, a tissue sectioning device, a microfluidic device, a blade grid, and/or a microtome.

As used herein, the term “solid support” can refer to discrete solid or semi-solid surfaces to which a plurality of stochastic barcodes may be attached. A solid support may encompass any type of solid, porous, or hollow sphere, ball, bearing, cylinder, or other similar configuration composed of plastic, ceramic, metal, or polymeric material (e.g., hydrogel) onto which a nucleic acid may be immobilized (e.g., covalently or non-covalently). A solid support may comprise a discrete particle that may be spherical (e.g., microspheres) or have a non-spherical or irregular shape, such as cubic, cuboid, pyramidal, cylindrical, conical, oblong, or disc-shaped, and the like. A plurality of solid supports spaced in an array may not comprise a substrate. A solid support may be used interchangeably with the term “bead.” As used herein, “solid support” and “substrate” can be used interchangeably.

As used here, the term “target” can refer to a composition which can be associated with a stochastic barcode. Exemplary suitable targets for analysis by the disclosed methods, devices, and systems include oligonucleotides, DNA, RNA, mRNA, microRNA, tRNA, and the like. Targets can be single or double stranded. In some embodiments targets can be proteins. In some embodiments targets are lipids. As used herein, “target” can be used interchangeably with “species”.

The term “reverse transcriptases” can refer to a group of enzymes having reverse transcriptase activity (i.e., that catalyze synthesis of DNA from an RNA template). In general, such enzymes include, but are not limited to, retroviral reverse transcriptase, retrotransposon reverse transcriptase, retroplasmid reverse transcriptases, retron reverse transcriptases, bacterial reverse transcriptases, group II intron-derived reverse transcriptase, and mutants, variants or derivatives thereof. Non-retroviral reverse transcriptases include non-LTR retrotransposon reverse transcriptases, retroplasmid reverse transcriptases, retron reverse transciptases, and group II intron reverse transcriptases. Examples of group II intron reverse transcriptases include the Lactococc s lactis Ll.LtrB intron reverse transcriptase, the Thermosynechococcus elongatus TeI4c intron reverse transcriptase, or the Geobacillus stearothermophilus GsI-IIC intron reverse transcriptase. Other classes of reverse transcriptases can include many classes of non-retroviral reverse transcriptases (i.e., retrons, group II introns, and diversity-generating retroelements among others).

Methods of Removing High Abundance Species

Some embodiments disclosed herein provide methods of removing high abundance species from a plurality of nucleic acid molecules. In some embodiment, the methods disclosed herein can remove high abundance species from a plurality of nucleic acid molecules without significantly removing the low abundance species or the intermediate abundance species from the plurality of nucleic acid molecules. As used herein, “significantly removing” refers to removing at least 10%, at least 20%, at least 30%, at least 40%, at least 50% or more of a low abundance species or intermediate abundance species from the plurality of nucleic acid molecules. In some embodiments, the methods disclosed herein can remove high abundance species and the intermediate abundance species from a plurality of nucleic acid molecules without significantly removing the low abundance species from the plurality of nucleic acid molecules.

As used herein, a “species” refers to the polynucleotides (for example, single-stranded polynucleotides) in the plurality of nucleic acid molecules that are the same or the complement of one another, or are capable of hybridize to one another, or are transcripts from the same genetic locus, or encode the same protein or fragment thereof, etc. In some embodiments, members of a species are at least 80%, at least 90%, at least 95%, at least 98%, at least 99%, or 100% homologous to one another, or complement thereof. In some embodiments, members of a species can hybridize to one another under high stringent hybridization conditions. In some embodiments, members of a species can hybridize to one another under moderate stringent hybridization conditions. In some embodiments, members of a species can hybridize to one another under low stringent hybridization conditions. In some embodiments, members of a species are transcripts from the same genetic locus and the transcripts can be of the same or different length. The species is, in some embodiments, cDNA or mRNA.

As used herein, a “high abundance species” refers to a species that is present in high amount in the plurality of nucleic acids, for example the species can represent at least 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 20%, 30%, 40%, 50%, or more of the plurality of nucleic acid molecules. In some embodiments, the plurality of nucleic acid molecules can comprise at least 1, at least 2, at least 3, at least 4, at least 5, at least 10, at least 20, at least 50, at least 100, at least 200, at least 500, at least 1,000, or more, high abundance species. In some embodiments, the total of all the high abundance species represent at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, or more of the plurality of nucleic acid molecules. In some embodiments, high abundance species can comprise polynucleotides encoding one or more ribosomal proteins. In some embodiments, high abundance species can comprise polynucleotides encoding one or more mitochondrial proteins. In some embodiments, high abundance species can comprise polynucleotides encoding one or more housekeeping proteins.

As used herein, “intermediate abundance species” refers to a species that is present in an amount in the plurality of nucleic acid that is lower than at least one species in the plurality of nucleic acid and is higher than at least one other species in the plurality of nucleic acid. In some embodiments, an intermediate abundance species can represent about 10%, 5%, 4%, 3%, 2%, 1%, 0.1%, 0.01%, or a range between any two of the above values, of the plurality of nucleic acid molecules. In some embodiments, the plurality of nucleic acid molecules can comprise at least 1, at least 2, at least 3, at least 4, at least 5, at least 10, at least 20, at least 50, at least 100, at least 200, at least 500, at least 1,000, or more, intermediate abundance species. In some embodiments, the total of all the intermediate abundance species represent about 1%, about 2%, about 3%, about 4%, about 5%, about 10%, about 20%, about 30%, or a range between any two of the above values, of the plurality of nucleic acid molecules.

As used herein, “low abundance species” refers to a species that is present in low amount in the plurality of nucleic acids, for example the species can represent less than 1%, 0.1%, 0.01%, 0.001%, 0.0001%, or less of the plurality of nucleic acid molecules. In some embodiments, the plurality of nucleic acid molecules can comprise at least 1, at least 2, at least 3, at least 4, at least 5, at least 10, at least 20, at least 50, at least 100, at least 200, at least 500, at least 1,000, or more, low abundance species. In some embodiments, the total of all the low abundance species represent less than 20%, less than 10%, less than 5%, less than 4%, less than 3%, less than 2%, less than 1%, less than 0.1%, or less of the plurality of nucleic acid molecules. In some embodiments, low abundance species can comprise polynucleotides encoding one or more transcription factors. In some embodiments, high abundance species can comprise polynucleotides encoding one or more T cell receptors. In some embodiments, high abundance species can comprise polynucleotides encoding one or more antibodies.

In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules. For example, the methods and compositions disclosed herein can remove at least 1, at least 2, at least 3, at least 4, at least 5, at least 10, at least 20, at least 50, at least 100, at least 200, at least 500, at least 1,000, or more, high abundance species. In some embodiments, the methods and compositions disclosed herein can remove at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of each of the one or more high abundance species from the plurality of nucleic acid molecules. In some embodiments, the methods and compositions disclosed herein can remove at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of at least one of the one or more high abundance species from the plurality of nucleic acid molecules. In some embodiments, the methods and compositions disclosed herein can remove at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the total of high abundance species from the plurality of nucleic acid molecules.

In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules without significantly removing the low abundance species or the intermediate abundance species from the plurality of nucleic acid molecules. In some embodiments, the methods and compositions disclosed herein can remove at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of each of the one or more high abundance species from the plurality of nucleic acid molecules without significantly removing the low abundance species or the intermediate abundance species from the plurality of nucleic acid molecules. In some embodiments, the methods and compositions disclosed herein can remove at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the total of high abundance species from the plurality of nucleic acid molecules without significantly removing the low abundance species or the intermediate abundance species from the plurality of nucleic acid molecules. In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules while keeping at least at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of each of the one or more low abundance species. In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules while keeping at least at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of at least one of the one or more of low abundance species. In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules while keeping at least at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the total of low abundance species. In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules while keeping at least at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of at least one of the one or more of intermediate abundance species. In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules while keeping at least at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the total of intermediate abundance species. In some embodiments, the methods and compositions disclosed herein can remove one or more high abundance species from the plurality of nucleic acid molecules while keeping at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of each of the intermediate abundance species from the plurality of nucleic acid molecules.

Plurality of Nucleic Acid Molecules

It would be appreciated by one of ordinary skill in the art that the plurality of nucleic acid molecules can comprise a variety of nucleic acid molecules. In some embodiments, the plurality of nucleic acid molecules can comprise, DNA molecules, RNA molecules, genomic DNA molecules, cDNA molecules, mRNA molecules, rRNA molecules, siRNA molecules, or a combination thereof, and can be double-stranded or single-stranded. In some embodiments, the plurality of nucleic acid molecules comprise at least 100, at least 1,000, at least 10,000, at least 20,000, at least 30,000, at least 40,000, at least 50,000, at least 100,000, at least 1,000,000, or more species. In some embodiments, the plurality of nucleic acid molecules can be from a sample, such as a single cell, or a plurality of cells. In some embodiments, the plurality of nucleic acid molecules can be pooled from a plurality of samples, such as a plurality of single cells.

In some embodiments, the plurality of nucleic acid molecules comprises an unnormalized nucleic acid library, a partially normalized nucleic acid library, or a nucleic acid library that has been normalized by other methods, such as a cDNA library, a genomic DNA library, or the like. In some embodiments, the plurality of nucleic acid molecules can comprise a pooled unnormalized nucleic acid library, such as a pooled unnormalized nucleic acid library constructed from a plurality of unnormalized nucleic acid libraries each representing a single cell. In some embodiments, the unnormalized nucleic acid library is a cDNA library. In some embodiments, the unnormalized nucleic acid library is a genomic library. In some embodiments, the unnormalized nucleic acid library is a single-cell nucleic acid library.

In some embodiments, the plurality of nucleic acid molecules can be subjected to amplification before removing the high abundance species. For example, the plurality of nucleic acid molecules can comprise an amplified nucleic acid library. In some embodiments, the plurality of nucleic acid molecules can comprise at least 2, at least 4, at least 8, at least 16, at least 100, at least 1,000 or more copies of each nucleic acid molecules.

Molecular Labels

It is contemplated that the methods and compositions disclosed herein can be used in conjunction of molecular labels, for example, nucleic acid molecules that comprise molecular labels. Accordingly, the species of nucleic acid molecules as disclosed herein can include polynucleotides in the plurality of nucleic acid molecules that are the same or the complement of one another, or are capable of hybridize to one another, or are transcripts from the same genetic locus, or encode the same protein or fragment thereof, etc., but that are associated with different molecular labels. It would be appreciated that molecular labels can be used to identify occurrences of a nucleic acid species, such as a high abundance species, a low abundance a species, and/or an intermediate abundance species. In some embodiments, a high abundance species can comprise nucleic acid molecules that comprise at least 100, at least 1,000, at least 10,000 or more different molecular labels. In some embodiments, a low abundance species can comprise nucleic acid molecules that comprise less than 100, less than 50, less than 20, less than 10, less than 5 or less different molecular labels. In some embodiments, an intermediate abundance species can comprise nucleic acid molecules that comprise about 10, about 20, about 50, about 100, or a range between any two of the above values, different molecular labels.

A molecular label may comprise a nucleic acid sequence that provides identifying information for the specific nucleic acid. A molecular label may comprise a nucleic acid sequence that provides a counter for the specific occurrence of the target nucleic acid. A molecular label may be at least about 1, 2, 3, 4, 5, 10, 15, 20, 25, 30, 35, 40, 45, 50 or more nucleotides in length. A molecular label may be at most about 300, 200, 100, 90, 80, 70, 60, 50, 40, 30, 20, 15, 12, 10, 9, 8, 7, 6, 5, 4 or fewer nucleotides in length.

It would be appreciated that in some embodiments, the methods and compositions disclosed herein may remove a high abundance species without significantly reducing the number of different molecular labels associated with the high abundance species. For example, the methods and compositions disclosed herein can remove a high abundance species while retaining at least at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the different molecular labels associated with the high abundance species. In some embodiments, the methods and compositions disclosed herein can remove at least at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or at least 99% of a high abundance species while retaining at least at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the different molecular labels associated with the high abundance species. In some embodiments, the removing of a high abundance species does not remove the number of different molecular labels associated with the high abundance species.

Accordingly, in some embodiments, the methods and compositions disclosed herein may remove a high abundance species without significantly reducing the number of different molecular labels associated with each of the intermediate or low abundance species. For example, the methods and compositions disclosed herein can remove a high abundance species while retaining at least at least at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the different molecular labels associated with an intermediate or low abundance species. In some embodiments, the methods and compositions disclosed herein can remove at least at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or at least 99% of a high abundance species while retaining at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the different molecular labels associated with an intermediate or low abundance species. In some embodiments, the removing of a high abundance species does not remove the number of different molecular labels associated with the high abundance species.

Affinity Moiety

In some embodiments, the methods disclosed herein comprise hybridizing a plurality of first oligonucleotides comprising an affinity moiety with the plurality of nuclei acid molecules. A variety of affinity moieties can be used for the methods and compositions disclosed herein. For example, an affinity moiety can be part of a binding pair. In some embodiments, the affinity moiety can be a functional group added to the oligonucleotides. In some embodiments, the affinity moiety can be biotin, streptavidin, heparin, an aptamer, a click-chemistry moiety, digoxigenin, primary amine(s), carboxyl(s), hydroxyl(s), aldehyde(s), ketone(s), or any combination thereof.

The affinity moieties as disclosed herein are capable of bind to a capture moiety such as a capture molecule. In some embodiments, the affinity moiety and capture molecule can be members of a binding pair, for example, biotin/streptavidin. The capture molecule can be immobilized on a solid support, such as a bead.

In some embodiments, the first oligonucleotides can be extended to generate a plurality of complementary strands of the plurality of nucleic acid targets comprising the affinity moiety. In some embodiments, a second strand can be synthesized using a primer that binds to a binding site on the complementary strands to produce double stranded nucleic acid molecules.

Remove High Abundance Species by Denaturation/Partial Reannealing

In some embodiments, removing of high abundance species can comprise denaturation followed by partial reannealing of the double stranded nucleic acid molecules, followed by removing the reannealed complementary strands of the plurality of nucleic acid targets by a capture molecule immobilized on one or more solid support, wherein the capture molecules specifically bind to the affinity moiety.

Denaturation can be performed by a variety of methods including heating the double stranded nucleic acid molecules, treating the double stranded nucleic acid molecules with organic solvents (e.g., DMS or formamide), changing the salt concentration of the double stranded nucleic acid molecules, and/or changing the pH of the double stranded nucleic acid molecules.

After denaturation, the single-stranded nucleic acid molecules can be partially reannealed. Partial reannealing can be performed by any method, for example, rapid cooling on ice, changing the salt concentration (e.g., reversing the salt concentration from the amount used in denaturation), and/or changing the pH (e.g., reversing the pH from the level used in denaturation), and the like.

It would be appreciated that the extent of reannealing can be adjusted according to the desired percentage of high abundance species to be removed, and/or the percentage of intermediate or low abundance species to be retained. Without being bound by any particular theory, it is believed that more abundant species (e.g., high abundance species) anneals faster than the species with lower abundance (e.g., intermediate and low abundance species) under the same anneal conditions. For example, by changing the temperature, salt concentration, pH, and/or duration of the reannealing step, the percentage of high abundance species to be removed, and/or the percentage of intermediate or low abundance species to be retained can be adjusted. In some embodiments, the temperature, salt concentration, pH, and/or duration of the reannealing step can be adjusted so that at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or at least 99% of a high abundance species reanneal. In some embodiments, the temperature, salt concentration, pH, and/or duration of the reannealing step can be adjusted so that at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the intermediate or low abundance species remain single-stranded.

In some embodiments, the molecular labels associated with the high abundance species, intermediate or low abundance species can be used as an indicator for the adjustment of the reannealing conditions. For example, the temperature, salt concentration, pH, and/or duration of the reannealing step can be adjusted so that at least 10%, at least 20%, at least 30%, at least 40%, at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the different molecular labels associated with one or more high abundance species are removed. In some embodiments, the temperature, salt concentration, pH, and/or duration of the reannealing step can be adjusted so that at least 50%, at least 60%, at least 70%, at least 80%, at least 90%, at least 95%, or 100% of the different molecular labels associated with one or more intermediate or low abundance species are retained.

Selective Depletion of High Abundance Species

In some embodiments, the methods disclosed herein comprise selective depletion of one or more high abundance species using target specific oligonucleotides comprising an affinity moiety. For example, the high abundance species can be removed by hybridizing an oligonucleotide that specifically binds to a high abundance species comprising an affinity moiety, followed by removing the high abundance species using a capture molecule immobilized on a solid support, wherein the capture molecule specifically binds to the affinity moiety. Exemplary high abundance species can comprise nucleic acid molecules that encode ribosomal proteins, mitochondrial proteins, housekeeping proteins, etc.

In some embodiments, Cas9 can be used to remove a high abundance species as described in Gu et al. Genome Biology (2016) 17:41, hereby incorporated by reference in its entirety for discussion of Cas nucleases-based removal of high abundance species.

Sequencing

In some embodiments, the normalized libraries disclosed herein may be used for sequencing. In some embodiments, the normalized libraries disclosed herein may be amplified prior to sequencing. Any suitable sequencing method known in the art can be used, preferably high-throughput approaches. For example, cyclic array sequencing using platforms such as Roche 454, Illumina Solexa, ABI-SOLiD, ION Torrent, Complete Genomics, Pacific Bioscience, Helicos, or the Polonator platform, may also be utilized. Sequencing may comprise MiSeq sequencing. Sequencing may comprise HiSeq sequencing.

It would be appreciated the normalized libraries can increase the efficiency of sequencing, by increasing the sequencing reads for intermediate or low abundance species in the normalized library, and/or decreasing the sequencing reads for high abundance species in the normalized library.

In some embodiments, after removing the high abundance species using the methods described herein, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 50%, less than 40%, less than 30%, less than 20%, less than 10%, less than 5%, or less, of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 40% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 30% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 20% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of high abundance nucleic acid targets is less than 10% of the total sequencing reads. The methods described herein for removing high abundance species can also improve the sequencing reads for intermediate and/or low abundance nucleic acid targets in a plurality of nucleic acids (e.g., a nucleic acid library). For example, the sequencing reads for the plurality of intermediate abundance nucleic acid targets can be at least 30%, at least 20%, at least 10%, at least 5%, of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets can be at least 50%, at least 40%, at least 30%, at least 20%, at least 10%, at least 5%, at least 4%, at least 3%, at least 2%, at least 1%, of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 5% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 10% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 20% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 30% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 40% of the total sequencing reads. In some embodiments, the sequencing reads for the plurality of low abundance nucleic acid targets is at least 50% of the total sequencing reads.

The methods and compositions disclosed herein can improve the sequencing efficiency of the normalized library by decreasing the sequencing reads:molecular label ratio of a high abundance species and/or increasing the sequencing reads:molecular label ratio of a low or intermediate abundance species. In some embodiments, the ratio of sequencing reads to molecular label for the plurality of high abundance species is less than 15, less than 10, less than 9, less than 8, less than 7, less than 6, less than 5, less than 4, less than 3, less than 2, less than 1. In some embodiments, the ratio of sequencing reads to molecular label for a high abundance species is less than 15, less than 10, less than 9, less than 8, less than 7, less than 6, less than 5, less than 4, less than 3, less than 2, less than 1. In some embodiments, the ratio of sequencing reads to molecular label for the plurality of low abundance species is at least 1, at least 2, at least 3, at least 4, at least 5, at least 6, at least 7, at least 8, at least 9, at least 10, or more. In some embodiments, the ratio of sequencing reads to molecular label for a low abundance species is at least 1, at least 2, at least 3, at least 4, at least 5, at least 6, at least 7, at least 8, at least 9, at least 10, or more. In some embodiments, the ratio of sequencing reads to molecular label for the plurality of intermediate abundance species is at least 1, at least 2, at least 3, at least 4, at least 5, at least 6, at least 7, at least 8, at least 9, at least 10, or more. In some embodiments, the ratio of sequencing reads to molecular label for an intermediate abundance species is at least 1, at least 2, at least 3, at least 4, at least 5, at least 6, at least 7, at least 8, at least 9, at least 10, or more.

Methods for Library Normalization

Methods and compositions described herein address the challenges of physical and enzymatic separation of ssDNA and dsDNA fractions during library normalization. As shown in FIG. 1, an unnormalized library contains a high abundance species 105 and a low abundance species 110. During initial library preparation, the library, or a fraction of the library, is asymmetrically labeled on one end with an affinity capture moiety 115 (such as biotin, etc.). The labeled double stranded nucleic acid molecules 120 are denatured 130 to generate single-stranded nucleic acid molecules including a high abundance species 135 and a low abundance species 140. The single-stranded nucleic acid molecules are partially reannealed 150 to form double-stranded molecules of the high abundance species 155 whereas the low abundance species remain single-stranded 160. After denaturation and partial reannealing, all of the labeled strands are captured on a support matrix 165, such as paramagnetic streptavidin beads, and the bound and unbound fractions are separated 170. Highly abundant sequences will have predominantly rehybridized, and both strands will be removed in the bound fraction. However, low abundance sequences will not have reannealed, so the complement of the labeled strand will be present in the unbound fraction. The unbound fraction including single-stranded high abundance species 175 and single-stranded low abundance species 180 would represent a normalized library 185, and could either be used directly or further amplified for downstream applications.

Library normalization strategies can be based on one of two principles. The first is hybridizing the library to another set of nucleic acids where the sequences are uniformly represented, such as the genomic DNA from the source organism, and retaining the hybridized fraction. The other approach relies on the concentration dependence of solution hybridization. When a set of dsDNA molecules are denatured, they will rehybridize at a rate proportional to the square of their original concentrations, as represented by the following equation:

$\frac{\partial C}{\partial t}={\mathrm{kC}}^2$

Wherein k=association constant, and C=concentration of a particular DNA species. Integration and plugging in ssDNA and ddDNA concentration:

$\left[\mathrm{dsDNA}\right]=\frac{\left[\mathrm{ssDNA}\right]}{1+\left[\mathrm{ssDNA}\right]\mathrm{kt}}$

Without experimentally verifying what k constant is with our library and annealing condition, a k value that was presented in the literature (using optimal conditions) was used:

$k=\frac{1}{5\times {10}^4\times \sqrt{\mathrm{length}}}$

Where length is according to WTA Bioanalyzer, or N1 amplicon length.

Researcher exploit this property for library normalization by denaturing a mixture and only allowing it to partially reanneal; proportionally, much more of the high concentration species will have rehybridized to dsDNA while less abundant species will still be predominantly single stranded. The dsDNA and ssDNA are then separated by physical or enzymatic means. Physical separation typically requires prohibitively large amounts of starting material, and the nuclease most typically used is of variable quality purified from a non-recombinant source and is prone to digesting ssDNA in regions of local secondary structure.

Computer modeling of the cost-reducing effect and rare-transcript identification effect was performed by normalization on existing targeted and WTA sequencing experiments on the Precise and Resolve platforms.

EXAMPLES

Some aspects of the embodiments discussed above are disclosed in further detail in the following examples, which are not in any way intended to limit the scope of the present disclosure.

Example 1: Simulated from ‘ES32 TCR Panel’ Sequencing Results

FIGS. 2-3 show computer simulated normalization results using ES32 TCR panel sequencing reads, which indicate % reads dedicated for low expressers or rare transcripts can be increased 100%+, potentially allowing higher sensitivity and precision in rare transcript counting.

With a smaller dynamic range in read ratio between high expressers and low expressers, targeted panel design can allow for high dynamic range transcripts.

The simulation results indicate normalization can reduce sequencing depth of high expressers and not the low expressers, yet retain the molecular label counting in oversequenced libraries.

Example 2: Simulated from ‘CBRepro’ Sequencing Results

FIGS. 4-6 show computer simulated normalization results using CBRepro sequencing reads, which indicate normalization can reduce sequencing cost by at least 25% for the same output data by removing high expressers such as mitochondrial RNA and ribosomal protein RNA.

FIG. 6 shows sequencing depth primarily affects highly abundant transcripts.

Example 3: Simulated from Actual T Cell Resolve WTA Data, Pretty Shallow Sequencing

FIGS. 7-9 show computer simulated normalization results using actual T Cell Resolve WTA sequencing reads. Because sequencing depth for this experiment was low, this simulation cannot predict the loss of MI during normalization.

FIG. 9 shows sequencing depth changes after normalization.

Example 4: Simulated from Actual T Cell Resolve WTA Data, Next-Generation Sequencing

FIGS. 10-12 show computer simulated normalization results using actual T Cell Resolve WTA data, next-generation reads. Normalization reduced mitochondrial/ribo RNA reads by 50%, increased rare transcript reads by >10-fold.

While this sequencing run is still pretty shallow, simulation shows no loss in MI counts in intermediate and low expresser genes.

While various aspects and embodiments have been disclosed herein, other aspects and embodiments will be apparent to those skilled in the art. The various aspects and embodiments disclosed herein are for purposes of illustration and are not intended to be limiting, with the true scope and spirit being indicated by the following claims.

One skilled in the art will appreciate that, for this and other processes and methods disclosed herein, the functions performed in the processes and methods can be implemented in differing order. Furthermore, the outlined steps and operations are only provided as examples, and some of the steps and operations can be optional, combined into fewer steps and operations, or expanded into additional steps and operations without detracting from the essence of the disclosed embodiments.

With respect to the use of substantially any plural and/or singular terms herein, those having skill in the art can translate from the plural to the singular and/or from the singular to the plural as is appropriate to the context and/or application. The various singular/plural permutations may be expressly set forth herein for sake of clarity.

It will be understood by those within the art that, in general, terms used herein, and especially in the appended claims (e.g., bodies of the appended claims) are generally intended as “open” terms (e.g., the term “including” should be interpreted as “including but not limited to,” the term “having” should be interpreted as “having at least,” the term “includes” should be interpreted as “includes but is not limited to,” etc.). It will be further understood by those within the art that if a specific number of an introduced claim recitation is intended, such an intent will be explicitly recited in the claim, and in the absence of such recitation no such intent is present. For example, as an aid to understanding, the following appended claims may contain usage of the introductory phrases “at least one” and “one or more” to introduce claim recitations. However, the use of such phrases should not be construed to imply that the introduction of a claim recitation by the indefinite articles “a” or “an” limits any particular claim containing such introduced claim recitation to embodiments containing only one such recitation, even when the same claim includes the introductory phrases “one or more” or “at least one” and indefinite articles such as “a” or “an” (e.g., “a” and/or “an” should be interpreted to mean “at least one” or “one or more”); the same holds true for the use of definite articles used to introduce claim recitations. In addition, even if a specific number of an introduced claim recitation is explicitly recited, those skilled in the art will recognize that such recitation should be interpreted to mean at least the recited number (e.g., the bare recitation of “two recitations,” without other modifiers, means at least two recitations, or two or more recitations). Furthermore, in those instances where a convention analogous to “at least one of A, B, and C, etc.” is used, in general such a construction is intended in the sense one having skill in the art would understand the convention (e.g., “a system having at least one of A, B, and C” would include but not be limited to systems that have A alone, B alone, C alone, A and B together, A and C together, B and C together, and/or A, B, and C together, etc.). In those instances where a convention analogous to “at least one of A, B, or C, etc.” is used, in general such a construction is intended in the sense one having skill in the art would understand the convention (e.g., “a system having at least one of A, B, or C” would include but not be limited to systems that have A alone, B alone, C alone, A and B together, A and C together, B and C together, and/or A, B, and C together, etc.). It will be further understood by those within the art that virtually any disjunctive word and/or phrase presenting two or more alternative terms, whether in the description, claims, or drawings, should be understood to contemplate the possibilities of including one of the terms, either of the terms, or both terms. For example, the phrase “A or B” will be understood to include the possibilities of “A” or “B” or “A and B.”

In addition, where features or aspects of the disclosure are described in terms of Markush groups, those skilled in the art will recognize that the disclosure is also thereby described in terms of any individual member or subgroup of members of the Markush group.

As will be understood by one skilled in the art, for any and all purposes, such as in terms of providing a written description, all ranges disclosed herein also encompass any and all possible subranges and combinations of subranges thereof. Any listed range can be easily recognized as sufficiently describing and enabling the same range being broken down into at least equal halves, thirds, quarters, fifths, tenths, etc. As a non-limiting example, each range discussed herein can be readily broken down into a lower third, middle third and upper third, etc. As will also be understood by one skilled in the art all language such as “up to,” “at least,” and the like include the number recited and refer to ranges which can be subsequently broken down into subranges as discussed above. Finally, as will be understood by one skilled in the art, a range includes each individual member. Thus, for example, a group having 1-3 cells refers to groups having 1, 2, or 3 cells. Similarly, a group having 1-5 cells refers to groups having 1, 2, 3, 4, or 5 cells, and so forth.

From the foregoing, it will be appreciated that various embodiments of the present disclosure have been described herein for purposes of illustration, and that various modifications may be made without departing from the scope and spirit of the present disclosure. Accordingly, the various embodiments disclosed herein are not intended to be limiting, with the true scope and spirit being indicated by the following claims.

Claims

1. A composition, comprising: a plurality of first oligonucleotides comprising an affinity moiety, wherein each of the plurality of first oligonucleotides is a non-target-specific primer comprising a universal primer binding site;a plurality of selective depletion oligonucleotides comprising an affinity moiety, wherein the plurality of selective depletion oligonucleotides specifically bind nucleic acid molecules that encode: (i) ribosomal proteins, and(ii) mitochondrial proteins encoded by mitochondrial mRNA; anda plurality of capture molecules immobilized on one or more solid supports, wherein the capture molecules are configured to specifically bind to the affinity moiety of the plurality of first oligonucleotides and the plurality of selective depletion oligonucleotides, wherein the plurality of first oligonucleotides, the plurality of selective depletion oligonucleotides, and the plurality of capture molecules immobilized on one or more solid supports are provided separately.

2. The composition of claim 1, wherein the affinity moiety of the plurality of first oligonucleotides and the plurality of selective depletion oligonucleotides is a functional group selected from the group consisting of biotin, streptavidin, heparin, an aptamer, a click-chemistry moiety, digoxigenin, primary amine(s), carboxyl(s), hydroxyl(s), aldehyde(s), ketone(s), and any combination thereof.

3. The composition of claim 1, wherein the affinity moiety of the plurality of first oligonucleotides and the plurality of selective depletion oligonucleotides is biotin and the capture molecule is streptavidin.

4. The composition of claim 1, wherein the affinity moiety of the plurality of first oligonucleotides and the plurality of selective depletion oligonucleotides is streptavidin and the capture molecule is biotin.

5. The composition of claim 1, wherein the solid support is a bead.

6. The composition of claim 1, wherein each of the plurality of first oligonucleotides comprises a stochastic barcode.

7. The composition of claim 1, comprising one or more reagents for nucleic acid extension reactions.

8. The composition of claim 7, comprising a reverse transcriptase.

9. The composition of claim 1, comprising one or more reagents for PCR reactions.

10. The composition of claim 1, comprising: instructions for performing: (i) hybridizing the plurality of first oligonucleotides and the plurality of selective depletion oligonucleotides with a first plurality of nucleic acid molecules comprising at least one nucleic acid molecule that encodes: (i) a ribosomal protein, and/or(ii) a mitochondrial protein encoded by mitochondrial mRNA;(ii) extending the first plurality of oligonucleotides and the plurality of selective depletion oligonucleotides to generate a plurality of complementary strands of the first plurality of nucleic acid molecules comprising the affinity moiety;(iii) denaturing a plurality of double-stranded nucleic acid molecules comprising the plurality of complementary strands of the first plurality of nucleic acid molecules;(iv) partially reannealing the plurality of complementary strands of the first plurality of nucleic acid molecules; and(v) removing the partially reannealed complementary strands of the first plurality of nucleic acid molecules by the plurality of capture molecules immobilized on one or more solid supports to generate a second plurality of nucleic acid molecules.

11. The composition of claim 10, further comprising a plurality of second oligonucleotides capable of hybridizing to the plurality of complementary strands of the first plurality of nucleic acid molecules.

12. The composition of claim 11, wherein each of the plurality of second oligonucleotides comprises a universal primer binding site.

13. The composition of claim 11, wherein each of the plurality of second oligonucleotides comprises a stochastic barcode.

14. The composition of claim 1, further comprising Cas9 protein.

15. The composition of claim 14, wherein the Cas9 protein is complexed with a guide oligonucleotide that specifically binds a nucleic acid molecule encoding a ribosomal protein or a mitochondrial protein.