High-throughput single-cell sequencing with reduced amplification bias

Abstract

A method is provided for preparing a sequencing library containing nucleic acids from a plurality of cells. [Solution] The method includes the steps of: (a) providing isolated nuclei from a plurality of cells; (b) distributing the isolated nuclei into a first plurality of compartments, each compartment containing a subset of the isolated nuclei, and contacting each subset with a transposome complex, the transposome complex in each compartment containing a transposase, a first index sequence that is different from first index sequences in other compartments, and at least one ubiquitous sequence; and (c) fragmenting nucleic acid in the subset of nucleosome-depleted nuclei into a plurality of nucleic acid fragments and incorporating the first index sequence into at least one strand of the nucleic acid fragment to generate indexed nuclei containing indexed nucleic acid fragments.

Description

[Background technology]

[0001] CROSS-REFERENCE TO RELATED APPLICATIONS This application claims the benefit of U.S. Provisional Patent Application No. 62 / 673,023, filed May 17, 2018, and U.S. Provisional Patent Application No. 62 / 821,864, filed March 21, 2019, each of which is incorporated by reference in its entirety. government funding

[0002] This invention was made with Government support under Grant No. DP1 HG007811 awarded by the National Institutes of Health. The Government has certain rights in this invention. [Technical field]

[0003] The embodiments of the present disclosure relate to nucleic acid sequencing.In particular, the embodiments of the methods and compositions provided herein relate to generating indexed single-cell sequencing libraries and obtaining sequence data therefrom to characterize rare events, including crossover and chromosome missegregation events.In some embodiments, the methods relate to revealing cancer heterogeneity at single-cell level.

[0004] Background technology Modern single-cell genome sequencing techniques have two major limitations. First, most methods require segmenting individual cells, which can limit throughput. Second, most amplification methods are PCR-based and therefore suffer from exponential amplification bias. To solve the first problem, we and colleagues developed single-cell combinatorial indexing ("sci-"), which performs several rounds of split-pool molecular barcoding to specifically tag the nucleic acid content of single cells, thereby allowing an exponential increase in throughput with each successive round of indexing. Sci-methods have been successfully developed to profile chromatin accessibility (sci-ATAC-seq), transcriptome (sci-RNA-seq), genome (sci-DNA-seq), methylome (sci-MET), and chromosomal conformation (sci-Hi-C) in large numbers of single cells (Cao et al., 2017, Science 357:661-667; Cusanovich et al., 2015, Science, 348:910-914; Mulqueen et al., 2018, Nat. Biotechnol. 36:428-431; Ramani et al., 2017, Nat. Methods 14:263-266; Vitak et al., 2017, Nat. Methods 14:302-308). To solve the second problem, linear amplification by T7-based transcription offers a potential solution that has been previously utilized in the context of single-cell assays (Eberwine et al., 1992; Proceedings of the National Academy of Sciences 89:3010-3014; Hashimshony et al., 2012, Cell Rep. 2:666-673; Sos et al., 2016, Genome Biolol., 17:20). For example, Chen et al. recently developed linear amplification by transposon insertion ("LIANTI"), which uses Tn5 transposons to fragment genomes and simultaneously insert T7 RNA promoters for in vitro transcription (IVT).The RNA copies generated from the DNA template cannot serve as templates for further amplification. Therefore, all copies are derived directly from the original DNA template. By avoiding exponential amplification, LIANTI maintains uniformity and minimizes sequence errors. However, this method is low-throughput because it requires continuous library preparation from each single cell (Chen et al., 2017, Science 356:189-194). Summary of the Invention

[0005] A method is described herein that incorporates single-cell combinatorial indexing and linear amplification to enable exponential increase in throughput while minimizing amplification bias. By multiple rounds of molecular barcoding, the method improves throughput to at least thousands and potentially millions of cells per experiment while retaining the advantages of linear amplification. The inventors demonstrate the generalizability of the method through proof-of-concept of single-cell whole genome sequencing ("sci-L3-WGS"), targeted genome sequencing ("sci-L3-targeted-seq") and genome and transcriptome co-assay ("sci-L3-RNA / DNA"). As a further demonstration, single-cell whole genome sequencing is applied to map an unprecedented number of meiotic crossovers and rare chromosome missegregation events in immature and mature male germ cells from infertile interspecific (B6 x Spretus) F1 male mice and fertile intraspecific (B6 x Cast) F1 male mice.

[0006] definition It is understood that terms used herein take their ordinary meaning in the relevant art unless otherwise specified. Some terms used herein and their meanings are set forth herein.

[0007] As used herein, the terms "organism" and "subject" are used interchangeably and refer to microorganisms (e.g., prokaryotes or eukaryotes), animals, and plants. An example of an animal is a mammal, such as a human.

[0008] As used herein, the term "cell type" is intended to identify cells based on morphology, phenotype, developmental origin or other known or recognizable characteristic cellular properties. A variety of different cell types can be obtained from a single organism (or from organisms of the same species). Exemplary cell types include, but are not limited to, gametes (female gametes, e.g., eggs or egg cells, and male gametes, e.g., sperm), epithelial ovary, ovarian fibroblasts, testes, urinary bladder, immune cells, B cells, T cells, natural killer cells, dendritic cells, cancer cells, eukaryotic cells, stem cells, blood cells, muscle cells, adipocytes, skin cells, neural cells, bone cells, pancreatic cells, endothelial cells, pancreatic epithelium, pancreatic alpha, pancreatic beta, pancreatic endothelium, bone marrow lymphoblasts, bone marrow B lymphoblasts, bone marrow macrophages, bone marrow erythroblasts, bone marrow dendritic, bone marrow adipocytes, bone marrow osteocytes, bone marrow chondrocytes, promyeloblasts, bone marrow megakaryoblasts, urinary bladder, brain B lymphocytes, brain glia, neurons, brain astrocytes, neuroectoderm, brain macrophages, brain microglia, brain epithelium, cortical neurons, brain fibroblast, breast epithelium, colon epithelium, colon B lymphocyte, breast epithelium, breast myoepithelium, breast fibroblast, colon enterocyte, cervical epithelium, breast ductal epithelium, tongue epithelium, tonsillar dendritic, tonsillar B lymphocyte, peripheral blood lymphoblast, peripheral blood T lymphoblast, peripheral blood skin T lymphocyte, peripheral blood natural killer, peripheral blood B lymphoblast, peripheral blood monocyte, peripheral blood myeloblast, peripheral blood monoblast, peripheral blood promyeloblast, peripheral blood macrophage, peripheral blood basophil, liver endothelium, liver mast, liver epithelium, liver B lymphocyte, splenic endothelium, splenic epithelium, splenic B lymphocyte, hepatocyte, liver, fibroblast, lung epithelium, bronchial epithelium, lung fibroblast, lung B lymphocyte, lung Schwann, lung squamous epithelium, lung macrophage, lung osteoblast, neuroendocrine, lung alveolar, gastric epithelium and gastric fibroblast.

[0009] As used herein, the term "tissue" refers to a collection or aggregation of cells that act together to perform one or more specific functions in an organism. The cells may be morphologically similar. Exemplary tissues include, but are not limited to, epididymis, eye, muscle, skin, tendon, vein, artery, blood, heart, spleen, lymph node, bone, bone marrow, lung, bronchus, trachea, intestine, small intestine, large intestine, colon, rectum, salivary gland, tongue, gallbladder, appendix, liver, pancreas, brain, stomach, skin, kidney, ureter, bladder, urethra, gonads, testes, ovaries, uterus, fallopian tubes, thymus, pituitary gland, thyroid, adrenal gland, or parathyroid gland. The tissues can be derived from any of a variety of organs of a human or other organism. The tissues can be healthy or unhealthy tissues. Examples of unhealthy tissues include, but are not limited to, malignant tumors in reproductive tissues, lung, breast, colorectal, prostate, nasopharynx, stomach, testes, skin, nervous system, bone, ovaries, liver, blood tissue, pancreas, uterus, kidney, lymphatic tissue, etc. Malignant tumors can be of various histological subtypes, e.g., carcinoma, adenocarcinoma, sarcoma, fibroadenocarcinoma, neuroendocrine, or undifferentiated.

[0010] As used herein, the term "nucleosome" refers to the basic repeating unit of chromatin. The human genome consists of several meters of DNA packed into the nucleus of a cell, with an average diameter of about 10 μm. In eukaryotic nuclei, DNA is packaged into a nuclear protein complex known as chromatin. A nucleosome (basic repeating unit of chromatin) typically contains about 146 base pairs of DNA wrapped approximately 1.7 times around a core histone octamer. The histone octamer consists of two copies each of histones H2A, H2B, H3 and H4. Nucleosomes are regularly spaced along the DNA in a beads-on-a-string fashion.

[0011] As used herein, the term "compartment" refers to an area or volume that separates or isolates one thing from another. Exemplary compartments include, but are not limited to, vials, tubes, wells, droplets, boluses, beads, containers, surface characteristics, or areas or volumes that are separated by physical forces, such as liquid flow, magnetism, electric current, and the like. In one embodiment, the compartment is a well of a multi-well plate, such as a 96- or 384-well plate. As used herein, a droplet can include hydrogel beads, which are beads for encapsulating one or more nuclei or cells, including hydrogel compositions or droplet-based microfluidics. In some embodiments, the droplet is a uniform droplet of hydrogel material or a hollow droplet with a polymer hydrogel shell. The droplet, whether uniform or hollow, may be capable of encapsulating one or more nuclei or cells.

[0012] As used herein, "transposome complex" refers to an integration enzyme and a nucleic acid that includes an integration recognition site. A "transposome complex" is a functional complex formed by a transposase and a transposase recognition site that can catalyze a transposition reaction (see, for example, Gunderson et al., WO2016 / 130704). Examples of integration enzymes include, but are not limited to, integrases or transposases. Examples of integration recognition sites include, but are not limited to, transposase recognition sites.

[0013] As used herein, the term "nucleic acid" is consistent with its use in the art and includes naturally occurring nucleic acids or functional analogs thereof. Particularly useful functional analogs are capable of sequence-specific hybridization to nucleic acids or can be used as templates for replicating specific nucleotide sequences. Naturally occurring nucleic acids generally have a backbone containing phosphodiester bonds. Analog structures can have alternating backbone linkages, including any of a variety known in the art. Naturally occurring nucleic acids generally have a deoxyribose sugar (e.g., found in deoxyribonucleic acid (DNA)) or a ribose sugar (e.g., found in ribonucleic acid (RNA)). Nucleic acids can contain any of a variety of analogs of these sugar moieties known in the art. Nucleic acids can include natural or unnatural bases. In this regard, natural deoxyribonucleic acids can have one or more bases selected from the group consisting of adenine, thymine, cytosine or guanine, and ribonucleic acids can have one or more bases selected from the group consisting of adenine, uracil, cytosine or guanine. Useful unnatural bases that may be included in nucleic acids are known in the art. Examples of unnatural bases include locked nucleic acids (LNA), bridged nucleic acids (BNA) and pseudo-complementary bases (Trilink Biotechnologies, San Diego, CA). LNA and BNA bases can be incorporated into DNA oligonucleotides to increase the strength and specificity of hybridization of the oligonucleotide. LNA and BNA bases and the use of such bases are known and routine to those skilled in the art.

[0014] As used herein, the term "target", when used in reference to a nucleic acid, is intended to be a semantic identifier for the nucleic acid in the context of the methods or compositions set forth herein and does not necessarily limit the structure or function of the nucleic acid beyond that otherwise expressly set forth. A target nucleic acid may be virtually any nucleic acid of known or unknown sequence. It may be, for example, a fragment of genomic DNA (e.g., chromosomal DNA), extrachromosomal DNA, such as a plasmid, cell-free DNA, RNA (e.g., RNA or non-coding RNA), protein (e.g., intracellular or cell surface protein), or cDNA. Sequencing can result in the determination of the sequence of the entire or part of the target molecule. Targets can be derived from a primary nucleic acid sample, such as a nucleus. In one embodiment, targets can be processed into a template suitable for amplification by the installation of a universal sequence at one or both ends of each target fragment. Targets can also be obtained from an initial RNA sample by reverse transcription into cDNA. In one embodiment, targets are used in reference to a subset of DNA, RNA, or proteins present in a cell. Targeted sequencing generally uses PCR amplification (e.g., region-specific primers) or hybridization-based capture methods (e.g., using capture probes) or antibodies to select and isolate genes or regions or proteins of interest. Targeted enrichment can occur at various stages of the method. For example, targeted RNA representation can be obtained using target-specific primers in the reverse transcription step or hybridization-based enrichment of subsets from more complex libraries. Examples are exome sequencing or L1000 assays (Subramanian et al., 2017, Cell, 171;1437-1452). Targeted sequencing can include any enrichment process known to those skilled in the art.

[0015] As used herein, the term "ubiquitous" when used to describe a nucleotide sequence refers to a region of sequence that is common to two or more nucleic acid molecules, where the molecules also have regions of sequence that are different from each other. A ubiquitous sequence present in different members of a collection of molecules can allow the capture of multiple different nucleic acids using a population of ubiquitous capture nucleic acids, such as a collection of capture oligonucleotides that are complementary to a portion of the ubiquitous sequence, such as a ubiquitous capture sequence. Non-limiting examples of ubiquitous capture sequences include sequences that are identical to or complementary to P5 and P7 primers. Similarly, a ubiquitous sequence present in different members of a collection of molecules can allow the replication (sequencing) or amplification of multiple different nucleic acids using a population of ubiquitous primers that are complementary to a portion of the ubiquitous sequence, such as a ubiquitous anchor sequence. Non-limiting examples of ubiquitous anchor sequences include sequences that are identical to or complementary to spacer sequences, such as sp1 and sp2. In one embodiment, the ubiquitous anchor sequence is used as the site to which the ubiquitous primer (e.g., the sequencing primer for read 1 or read 2) anneals for sequencing. Thus, the capture oligonucleotide or ubiquitous primer comprises a sequence that can specifically hybridize to the ubiquitous sequence.

[0016] The terms "P5" and "P7" can be used when referring to a ubiquitous capture sequence or capture oligonucleotide. The terms "P5'" (P5 prime) and "P7'" (P7 prime) refer to the complements of P5 and P7, respectively. It will be understood that any suitable ubiquitous capture sequence or capture oligonucleotide can be used in the methods presented herein, and the use of P5 and P7 is only an exemplary embodiment. The use of capture oligonucleotides such as P5 and P7 or their complements on flow cells is known in the art, as exemplified by the disclosures of WO2007 / 010251, WO2006 / 064199, WO2005 / 065814, WO2015 / 106941, WO1998 / 044151 and WO2000 / 018957. For example, any suitable forward amplification primer, whether immobilized or in solution state, can be useful in the method presented herein for hybridization with complementary sequence and amplification of sequence.Similarly, any suitable reverse amplification primer, whether immobilized or in solution state, can be useful in the method presented herein for hybridization with complementary sequence and amplification of sequence.Those skilled in the art will understand how to design and use suitable primer sequence for capturing and / or amplifying nucleic acid as presented herein.

[0017] As used herein, the term "primer" and its derivatives generally refer to any nucleic acid that can hybridize to a target sequence of interest. Generally, a primer serves as a substrate onto which nucleotides can be polymerized by a polymerase or to which a nucleotide sequence, such as an index, can be ligated. However, in some embodiments, a primer can be incorporated into a synthesized nucleic acid strand and provide a site to which another primer can hybridize to prime the synthesis of a new strand complementary to the synthesized nucleic acid molecule. A primer can contain any combination of nucleotides or analogs thereof. In some embodiments, a primer is a single-stranded oligonucleotide or polynucleotide. The terms "polynucleotide" and "oligonucleotide" are used interchangeably herein to refer to a polymeric form of nucleotides of any length and can include ribonucleotides, deoxyribonucleotides, analogs thereof, or mixtures thereof. The term should be understood to include, as equivalents, antibody-oligoconjugates made from DNA, RNA, cDNA analogs, or nucleotide analogs, as well as being applicable to single-stranded (such as sense or antisense) and double-stranded polynucleotides. The term as used herein also encompasses cDNA, which is a complementary or copy DNA generated from an RNA template, for example, by the action of reverse transcriptase.The term refers only to the primary structure of the molecule.Thus, the term includes triple-stranded, double-stranded and single-stranded deoxyribonucleic acid ("DNA"), as well as triple-stranded, double-stranded and single-stranded ribonucleic acid ("RNA").

[0018] As used herein, the term "adapter" and its derivatives, such as ubiquitous adapters, generally refer to any linear oligonucleotide that can be ligated to a nucleic acid molecule of the present disclosure. In some embodiments, the adapter is substantially non-complementary to either the 3' or 5' end of any target sequence present in the sample. In some embodiments, suitable adapter lengths are in the range of about 10-100 nucleotides, about 12-60 nucleotides, or about 15-50 nucleotides. In general, the adapter can comprise any combination of nucleotides and / or nucleic acids. In some aspects, the adapter can comprise one or more cleavable groups at one or more positions. In another aspect, the adapter can comprise a sequence that is substantially identical to or substantially complementary to at least a portion of a primer, such as a ubiquitous primer. In some embodiments, the adapter can comprise a barcode (also referred to herein as a tag or index) to aid in downstream error correction, identification, or sequencing. The terms "adapter" and "adaptor" are used interchangeably.

[0019] As used herein, the term "each," when used in reference to a collection of items, unless the context clearly dictates otherwise, identifies each individual item in the collection, but does not necessarily refer to every item in the set.

[0020] As used herein, the term "transport" refers to the movement of molecules through a fluid. The term can include passive transport, such as the movement of molecules along their concentration gradient (e.g., passive diffusion). The term can also include active transport, in which molecules can move along or against their concentration gradient. Thus, transport can include the application of energy to move one or more molecules in a desired direction or to a desired location, such as an amplification site.

[0021] As used herein, "amplify", "amplifying" or "amplification reaction" and their derivatives generally refer to any act or process in which at least a portion of a nucleic acid molecule is replicated or copied into at least one additional nucleic acid molecule. The additional nucleic acid molecule may comprise a sequence that is substantially identical or substantially complementary to at least a portion of the template nucleic acid molecule. The template nucleic acid molecule may be single-stranded or double-stranded, and the additional nucleic acid molecule may be independently single-stranded or double-stranded. The amplification may comprise linear or exponential replication of the nucleic acid molecule. In some embodiments, such amplification may be performed using isothermal conditions. In other embodiments, such amplification may comprise thermocycling. In some embodiments, the amplification is a multiplex amplification that comprises simultaneous amplification of multiple target sequences in a single amplification reaction. In some embodiments, "amplification" comprises amplification of at least a portion of DNA and RNA-based nucleic acids, alone or in combination. The amplification reaction may comprise any of the amplification processes known to those skilled in the art. In some embodiments, the amplification reaction comprises polymerase chain reaction (PCR).

[0022] As used herein, "amplification conditions" and its derivatives generally refer to conditions suitable for amplifying one or more nucleic acid sequences. Such amplification may be linear or exponential. In some embodiments, amplification conditions may include isothermal conditions, or may include thermocycling conditions or a combination of isothermal and thermocycling conditions. In some embodiments, conditions suitable for amplifying one or more nucleic acid sequences include polymerase chain reaction (PCR) conditions. Generally, amplification conditions refer to a reaction mixture that is sufficient to amplify a nucleic acid, such as one or more target sequences flanked by ubiquitous sequences, or to amplify an amplified target sequence that is ligated to one or more adapters. Generally, amplification conditions include a catalyst for amplification or nucleic acid synthesis, such as a polymerase; a primer that has some complementarity to the nucleic acid to be amplified; and nucleotides, such as deoxyribonucleotide triphosphates (dNTPs), which facilitate the extension of a primer hybridized to the nucleic acid. Amplification conditions may require hybridization or annealing of the primer to the nucleic acid, extension of the primer, and a denaturation step in which the extended primer is separated from the nucleic acid sequence being amplified. Typically, but not necessarily, amplification conditions can include thermocycling, and in some embodiments, amplification conditions include multiple cycles in which the steps of annealing, extension and separation are repeated. Typically, amplification conditions include, but are not limited to, thermocycling, and in some embodiments, amplification conditions include multiple cycles in which the steps of annealing, extension and separation are repeated. 2+ or Mn 2+ and may also include various regulators of ionic strength.

[0023] As used herein, "reamplification" and its derivatives generally refer to any process in which at least a portion of an amplified nucleic acid molecule is further amplified through any suitable amplification process (in some embodiments, referred to as a "secondary" amplification), thereby generating a re-amplified nucleic acid molecule. The secondary amplification need not be identical to the original amplification process that generated the amplified nucleic acid molecule; the re-amplified nucleic acid molecule need not be completely identical or completely complementary to the amplified nucleic acid molecule; all that is required is that the re-amplified nucleic acid molecule contains at least a portion of the amplified nucleic acid molecule or its complement. For example, reamplification can include different amplification conditions and / or the use of different primers, including different target-specific primers than the primary amplification.

[0024] As used herein, the term "polymerase chain reaction" ("PCR") refers to the method of Mullis, U.S. Patent Nos. 4,683,195 and 4,683,202, which describes a method for increasing the concentration of a segment of a polynucleotide of interest in a mixture of genomic DNA without cloning or purification. This process for amplifying a polynucleotide of interest consists of introducing a large excess of two oligonucleotide primers to a DNA mixture containing the desired polynucleotide of interest, followed by a series of thermal cycles in the presence of a DNA polymerase. The two primers are complementary to their respective strands of the double-stranded polynucleotide of interest. The mixture is first denatured at a higher temperature, and the primers are then annealed to complementary sequences within the polynucleotide molecule of interest. After annealing, the primers are extended with a polymerase to form a new pair of complementary strands. The steps of denaturation, primer annealing and polymerase extension can be repeated many times (called thermocycling) to obtain a high concentration of the amplified segment of the desired polynucleotide of interest. The length of the amplified segment (amplicon) of the desired polynucleotide of interest is determined by the relative positions of the primers with respect to each other, and therefore this length is a controllable parameter. Because this process is repeated, the method is called PCR. Because the desired amplified segments of the polynucleotide of interest become the predominant nucleic acid sequences (in terms of concentration) in the mixture, they are said to be "PCR amplified". In a modification of the above method, the target nucleic acid molecule can be PCR amplified using multiple different primers, in some cases using one or more primer pairs per target nucleic acid molecule of interest, thereby forming a multiplex PCR reaction.

[0025] As defined herein, "multiplex amplification" refers to the selective, non-random amplification of two or more target sequences in a sample using at least one target-specific primer. In some embodiments, the multiplex amplification is performed such that some or all of the target sequences are amplified in a single reaction vessel. The "plexi" or "plex" of a given multiplex amplification generally refers to the number of different target-specific sequences amplified during that single multiplex amplification. In some embodiments, the plexi may be about 12-plex, 24-plex, 48-plex, 96-plex, 192-plex, 384-plex, 768-plex, 1536-plex, 3072-plex, 6144-plex or more. The amplified target sequences can also be detected by several different methods (e.g., gel electrophoresis followed by densitometry, quantification by bioanalyzer or quantitative PCR, hybridization with a labeled probe; incorporation of biotinylated primers followed by avidin-enzyme conjugate detection; detection of the amplified target sequences ... 32 P-labeled deoxynucleotide triphosphate incorporation).

[0026] As used herein, "amplified target sequence" and its derivatives generally refer to a nucleic acid sequence generated by amplifying a target sequence using the target-specific primers and methods provided herein. The amplified target sequence may be either of the same sense (i.e., the positive strand) or of the antisense (i.e., the negative strand) relative to the target sequence.

[0027] As used herein, the terms "ligate", "ligation" and their derivatives generally refer to a process for covalently linking two or more molecules, e.g., covalently linking two or more nucleic acid molecules to one another. In some embodiments, ligation involves joining nicks between adjacent nucleotides of a nucleic acid. In some embodiments, ligation involves forming a covalent bond between an end of a first nucleic acid molecule and an end of a second nucleic acid molecule. In some embodiments, ligation can involve forming a covalent bond between a 5' phosphate group of one nucleic acid and a 3' hydroxyl group of a second nucleic acid, thereby forming a ligated nucleic acid molecule. Generally, for purposes of this disclosure, an amplified target sequence can be ligated to an adaptor to generate an adaptor-ligated amplified target sequence.

[0028] As used herein, "ligase" and its derivatives generally refer to any agent capable of catalyzing the ligation of two substrate molecules. In some embodiments, a ligase includes an enzyme capable of catalyzing the joining of nicks between adjacent nucleotides of a nucleic acid. In some embodiments, a ligase includes an enzyme capable of catalyzing the formation of a covalent bond between the 5' phosphate of one nucleic acid molecule and the 3' hydroxyl of another nucleic acid molecule, thereby forming a ligated nucleic acid molecule. Suitable ligases can include, but are not limited to, T4 DNA ligase, T4 RNA ligase, and E. coli DNA ligase.

[0029] As used herein, "ligation conditions" and derivatives thereof generally refer to conditions suitable for ligating two molecules to each other. In some embodiments, the ligation conditions are suitable for sealing a nick or gap between nucleic acids. As used herein, the term nick or gap is consistent with the use of the term in the art. Generally, a nick or gap can be ligated in the presence of an enzyme such as a ligase at an appropriate temperature and pH. In some embodiments, T4 DNA ligase can ligate a nick between nucleic acids at a temperature of about 70-72°C.

[0030] As used herein, the term "flow cell" refers to a chamber that includes a solid surface through which one or more liquid reagents can flow.Examples of flow cell and related liquid systems and detection platforms that can be readily used in the method of the present disclosure are described in, for example, Bentley et al., Nature 456:53-59 (2008), WO 04 / 018497; U.S. Patent No. 7,057,026; WO 91 / 06678; WO 07 / 123744; U.S. Patent No. 7,329,492; U.S. Patent No. 7,211,414; U.S. Patent No. 7,315,019; U.S. Patent No. 7,405,281 and U.S. Patent Application Publication No. 2008 / 0108082.

[0031] As used herein, the term "amplicon" when used in reference to a nucleic acid refers to a copy product of a nucleic acid, where the product has a nucleotide sequence that is the same as or complementary to at least a portion of the nucleotide sequence of the nucleic acid. An amplicon can be generated by any of a variety of amplification methods using a nucleic acid or its amplicon as a template, including, for example, polymerase extension, polymerase chain reaction (PCR), rolling circle amplification (RCA), ligation extension or ligation chain reaction. An amplicon can be a nucleic acid molecule that has a single copy of a particular nucleotide sequence (e.g., a PCR product) or multiple copies of a nucleotide sequence (e.g., a concatemer product of RCA). A first amplicon of a target nucleic acid is generally a complementary copy. Subsequent amplicons are copies made from the target nucleic acid or from the first amplicon after the generation of the first amplicon. Subsequent amplicons can have a sequence that is substantially complementary to the target nucleic acid or substantially identical to the target nucleic acid.

[0032] As used herein, the term "amplification site" refers to a site in or on an array at which one or more amplicons can be generated. An amplification site can be further configured to contain, hold, or attach at least one amplicon generated at the site.

[0033] As used herein, the term "array" refers to a collection of sites that can be distinguished from one another by their relative positions. Different molecules at different sites of an array can be distinguished from one another by the site's position in the array. An individual site of an array can contain one or more molecules of a particular type. For example, a site can contain a single target nucleic acid molecule with a particular sequence, or a site can contain several nucleic acid molecules with the same sequence (and / or its complementary sequence). The sites of an array can be different features located on the same substrate. Exemplary features include, but are not limited to, wells in a substrate, beads (or other particles) in or on a substrate, protrusions from a substrate, bumps on a substrate, or channels in a substrate. The sites of an array can be separate substrates each with different molecules. The different molecules attached to the separate substrates can be identified by the substrate's position on a surface with which it is associated, or by the substrate's position in a liquid or gel. Exemplary arrays with separate substrates located on a surface include, but are not limited to, those with beads in wells.

[0034] As used herein, the term "capacity" refers to the maximum amount of nucleic acid material that can occupy a site when used in relation to a site and nucleic acid material.For example, the term can refer to the total number of nucleic acid molecules that can occupy a site under a particular condition.Other measures can also be used, including, for example, the total mass of nucleic acid material that can occupy a site under a particular condition or the total number of copies of a particular nucleotide sequence.In general, the capacity of a site for a target nucleic acid is substantially equivalent to the capacity of the site for an amplicon of the target nucleic acid.

[0035] As used herein, the term "capture agent" refers to a material, chemical, molecule or portion thereof that is capable of attaching to, holding, or binding to a target molecule (e.g., a target nucleic acid). Exemplary capture agents include, but are not limited to, a capture nucleic acid (also referred to herein as a capture oligonucleotide) that is complementary to at least a portion of a target nucleic acid, a member of a receptor-ligand binding pair (e.g., avidin, streptavidin, biotin, lectins, carbohydrates, nucleic acid binding proteins, epitopes, antibodies, etc.) that is capable of binding to a target nucleic acid (or a linking moiety attached thereto), or a chemical agent that is capable of forming a covalent bond with a target nucleic acid (or a linking moiety attached thereto).

[0036] As used herein, the term "reporter moiety" can refer to any identifiable tag, label, index, barcode, or group that allows for determining the composition, identity, and / or source of the analyte being investigated. In some embodiments, the reporter moiety can comprise an antibody that specifically binds to a protein. In some embodiments, the antibody can comprise a detectable label. In some embodiments, the reporter can comprise an antibody or affinity reagent labeled with a nucleic acid tag. The nucleic acid tag can be detectable, for example, through a proximity ligation assay (PLA) or a proximity extension assay (PEA) or a sequencing-based readout (Shahi et al. Scientific Reports volume 7, Article number: 44447, 2017), or CITE-seq (Stoeckius et al. Nature Methods 14:865-868, 2017).

[0037] As used herein, the term "clonal population" refers to a population of nucleic acids that are homogeneous with respect to a particular nucleotide sequence. A homogeneous sequence is generally at least 10 nucleotides in length, but may be even longer, including, for example, at least 50, 100, 250, 500 or 1000 nucleotides in length. A clonal population can be derived from a single target or template nucleic acid. Generally, all of the nucleic acids in a clonal population have the same nucleotide sequence. It will be understood that a small number of mutations (e.g., due to amplification artifacts) can occur in a clonal population without departing from clonality.

[0038] As used herein, the term "unique molecular identifier" or "UMI" refers to a molecular tag that can be attached to a nucleic acid molecule, whether random, non-random, or semi-random. When incorporated into a nucleic acid molecule, the UMI can be used to correct for subsequent amplification bias by directly counting the unique molecular identifiers (UMIs) that are sequenced after amplification.

[0039] As used herein, "providing" in the context of a composition, article, nucleic acid, or nucleus means making the composition, article, nucleic acid, or nucleus, purchasing the composition, article, nucleic acid, or nucleus, or otherwise obtaining the compound, composition, article, or nucleus.

[0040] The term "and / or" means one or all of the listed elements or a combination of any two or more of the listed elements.

[0041] The words "preferred" and "preferably" refer to embodiments of the present disclosure that may provide certain benefits, under certain circumstances. However, other embodiments may be preferred, under the same or other circumstances. Moreover, the recitation of one or more preferred embodiments does not imply that other embodiments are not beneficial and does not exclude other embodiments from the scope of the present disclosure.

[0042] The terms "comprises" and variations thereof do not have a limiting meaning where these terms appear in the description and claims.

[0043] Whenever embodiments are described herein in conjunction with words such as "comprises," "including," or "comprising," it is understood that otherwise similar embodiments described in terms of "consisting of" and / or "consisting essentially of" are also provided.

[0044] Unless otherwise stated, "a," "an," "the," and "at least one" are used interchangeably to mean one or more.

[0045] Further herein, the recitations of numerical ranges by endpoints include all numbers subsumed within that range (eg, 1 to 5 includes 1, 1.5, 2, 2.75, 3, 3.80, 4, 5, etc.).

[0046] For any method disclosed herein that includes distinct steps, the steps can be performed in any practicable order. Further, any combination of two or more steps can be performed simultaneously, where appropriate.

[0047] Throughout this specification, references to "one embodiment," "an embodiment," "a particular embodiment," or "some embodiments," etc., mean that the particular features, configurations, compositions, or characteristics described in connection with that embodiment are included in at least one embodiment of the disclosure. Thus, the appearances of such phrases in various places throughout this specification do not necessarily refer to the same embodiment of the disclosure. Furthermore, the particular features, configurations, compositions, or characteristics may be combined in any suitable manner in one or more embodiments.

[0048] The following detailed description of the exemplary embodiments of the present disclosure can be best understood when read in conjunction with the following figures. [Brief description of the drawings]

[0049] [Figure 1] 1A-B show a general block diagram of a general exemplary method for single-cell combinatorial indexing according to the present disclosure. [Diagram 2] FIG. 2 shows a general block diagram of a general exemplary method for single-cell combinatorial indexing according to the present disclosure. [Diagram 3]Figure 3A-F shows that sci-L3-WGS enables high-throughput single-cell linear whole genome amplification. (A) Schematic of the sci-L3-WGS workflow with three levels of indexing. (B) Top: Barcode structures of the resulting amplified DNA duplexes, compatible with various library preparation methods. bc, barcode; sp, spacer; gDNA, genomic DNA. Middle: Example of library structure for sci-L3-WGS. P5 and P7 sequencing adaptors are added by A-tailing and ligation. Note that it is equally possible to have P7 at the UMI end and P5 at the gDNA end due to the left-right symmetry of the ligation. Bottom: Example of library structure for sci-L3-targeted-seq. P5 and P7 sequencing adaptors are added by A-tailing and ligation, respectively, priming from spacer 2 (sp2) and the targeted locus of interest in the genome. Note that a new third barcode bc3' corresponding to each bc3 in the WGS library is also added by PCR, and a new UMI' is added outside of bc3'. (C) Scatter plot of the number of unique Tn5 insertion sites from human and mouse cells at low sequencing depth, 24bc1 × 64bc2 × 6bc3 sci-L3-WGS, sorting 100-300 cells per well. Blue, inferred mouse cell numbers (percentage of mouse reads >95%, median 98.7%, n=315); red, inferred human cell numbers (percentage of human reads >95%, median 99.8%, n=719); grey, inferred collisions (n=48,4%). (D) Box diagram showing the number of unique Tn5 insertion sites per cell at an average of 2.4M raw reads per cell and 1.78× depth. Depth is defined as the ratio between the number of unique IVT transcripts and the number of unique Tn5 insertion sites. Thick horizontal line, median; upper and lower edges of the box, first and third quartiles, respectively; whiskers, 1.5 times the interquartile range; circles, outliers). For characterization of libraries generated with an improved version of the protocol, see also Figure 5 and Example 2, section "Methods and molecular design of sci-L3-WGS and sci-L3-targeted-seq". (E) Example of chromosomal CNV diagram for individual cells.Top, HEK293T cells, 2.6M raw reads, 2.4M unique molecules, 1.3M unique Tn5 insertion sites with MAPQ>1. Bottom, 3T3 cells, 2.7M raw reads, 2.4M unique molecules, 1.2M unique Tn5 insertion sites with MAPQ>1. (F) Box diagram for copy number variation across 822 293T cells or 1,453 HAP1 cells. The Y-axis shows the percentage of reads per chromosome normalized by chromosome length, such that euploid chromosomes with no segmental copy gain or loss are expected to have a value of 1. [Figure 4]Figure 4A-F shows the molecular structure of sci-LIANTI in each step. Dashed line: RNA, solid line: DNA. (A) Tn5 adapters are phosphorylated at both 5' ends, one required for insertion and the other for ligation. The overhang of the annealed transposon contains the first-round barcode ("bc1") and spacer ("sp1") for ligation. (B) The ligation molecule is preannealed as a hairpin loop, which reduces intermolecular ligation from three molecules to two; the hairpin structure also helps improve RT efficiency in downstream steps. The hairpin contains 1) an overhang that anneals to "sp1" for ligation, 2) a second-round barcode ("bc2") and spacer ("sp2") that serve as a priming site in the stem for SSS in downstream steps, and 3) a T7 promoter in the loop for IVT. (C) Gap extension converts the looped T7 promoter into a double strand. Note that if ligation is successful at both ends, the T7 promoter is present on both sides; however, if ligation is successful at one end, the box portion is missing. Nevertheless, both can be reverse transcribed in downstream steps by different RT primers. (D) IVT generates single-stranded RNA amplicons downstream of the T7 promoter. (E) If ligation is successful at both ends, RT is preferably primed by self-looping RT primers, which are inherited from the looped ligation molecule; if ligation is successful at only one end, RT is primed by an additional RNA RT primer that is added in excess. Then, to avoid interference with the subsequent SSS reaction, the excess RNA primer is removed before SSS. (F) Double-stranded DNA molecules are generated by SSS to de-prime "sp2" to simultaneously add a third barcode and UMI tag to each transcript. A more detailed explanation is provided in Example 2, "Methods and molecular design of sci-L3-WGS and sci-L3-targeted-seq" section. [Diagram 5]5A-G show the number of reads in different sci-L3-WGS experiments with different Tn5 transposome concentrations. Box plot showing the number of unique Tn5 insertion sites per cell at the indicated depth. Depth is defined as the ratio between the number of unique IVT transcripts and the number of unique Tn5 insertion sites. Thick horizontal line, median; upper and lower edges of the box, first and third quartiles, respectively; whiskers, 1.5 times the interquartile range; circles, outliers). Concentrated Tn5 transposomes: 0.2 µM, diluted Tn5 transposomes: 0.1 µM. (A) Human vs. mouse specific reads for yi128 (median depth: 1.19x) with enriched Tn5 (median human specific reads: 215k, n=115 cells; median mouse specific reads: 169k, n=44); human specific reads for enriched Tn5 (median specific reads: 215k) vs. diluted Tn5 (median specific reads: 46k). (B) human specific reads for yi129 (median depth: 1.78x) with enriched Tn5 (median specific reads: 635k) vs. diluted Tn5 (median specific reads: 183k). Mouse specific reads are provided in Figure 3D. (C) yi140 and yi141 (median depth: 1.37x; median human specific reads: 660k) with enriched Tn5. See also Table 2 and Example 2. (D) yi144 and yi145 with enriched Tn5 (median depth: 1.05x; median human specific reads: 97.3k). See also Table 2. Note that yi140, yi141, yi144 and yi145 are libraries from the optimized protocol discussed in Example 2. (E) Human / mouse specific reads of yi174 (median depth: 1.06x) with enriched Tn5 (median human specific reads: 100k, n=103; median mouse specific reads: 23k, n=35); human specific reads with enriched Tn5 (median specific reads: 100k) and diluted Tn5 (median specific reads: 54k). (F) Mouse embryonic stem cell libraries: yi186, yi187, yi188 prepared with diluted Tn5; yi190, yi192, yi193 prepared with concentrated Tn5. (G) Number of unique Tn5 insertion sites as a function of sequencing depth.Blue and red lines indicate sci-L3-WGS with and without RNA RT primers, respectively (Example 2). Human unique insertions (median unique insertions: 635k) enriched for yi129 (as in panel B, median depth: 1.78x). When estimated at 5x and 10x depth, the expected number of unique insertions is 1.9M and 2.6M, respectively. yi140 and yi141 combined have a median depth of 1.37x and a median unique insertion of 660k. When estimated at 1.78x, 5x and 10x depth, the expected number of unique insertions is 1.5M, 4.2M and 6.0M, respectively. [Figure 6]Figure 6A-E show that Sci-L3-based RNA / DNA co-assay jointly enables high-throughput linear amplification for genome and transcriptome from the same single cell. (A) Schematic of sci-L3-RNA / DNA co-assay workflow with three levels of indexing. Note that both Tn5 transposon and cDNA synthesis primer contain the same phosphorylated ligation landing pad (pink) on the outer 5' overhang of the first-round barcode. (B) Resulting amplified duplex barcode structures corresponding to genome and transcriptome (left and right, respectively) compatible with various library preparation methods. bc, barcode; sp, spacer; gDNA, genomic DNA. (C) Scatter plot of the number of unique Tn5 insertion sites from human and mouse cells at low and high sequencing depth plotted together, 24bc1 × 64bc2 × 6bc3 sci-L3-RNA / DNA co-assay, 100–300 cells sorted per well. Blue, inferred mouse cell counts (percentage of mouse reads >95%, median 99.5%, n=2002); red, inferred human cell counts (percentage of human reads >95%, median 99.8%, n=2419); grey, inferred collisions (n=149, 6.6% low and high depth combined; 5 / 270, 3.7% high depth). (D) Same as in (C) for RNA. Blue, inferred mouse cells (median purity of mouse reads 95.1%); red, inferred human cells (median purity of human reads 91.5%); grey, inferred collisions (n=272, 12% low and high depth combined; 7 / 270, 5.2% high depth). (E) Seurat with RNA-seq signal shows clear clusters corresponding to BJ-5ta human dermal fibroblasts (male) and HEK293T (female) cells. Based on the presence or absence of a Y chromosome, 988 / 1024 cells (96.5%) were correctly assigned. [Figure 7]Figure 7A-E show mitotic / equal and meiotic / reductive chromosome segregation with and without crossover. Each vertical segment represents one chromatid (DNA strand not shown). Black and blue colors represent homologs. Ovals represent centromeres. It should be noted that mouse chromosomes are acrocentric. Grey crossovers indicate crossover sites after DNA replication in phase 4C. Red boxes indicate mitotic daughter cells that are heterozygous, and black and blue boxes indicate meiosis I (MI) daughter cells that are homozygous for the respective lineage background at the centromere-proximal region. LOH regions in daughter cells are characterized by curly brackets. (A) Mitotic / equal segregation without crossover. Both daughter cells retain heterozygosity. (B) Mitotic / equal segregation with crossover between homologs. Recombined chromatids separate, resulting in LOH centromeric distal to the crossover. (C) Mitotic / equal segregation with crossover between homologs. The recombined chromatids segregate together such that both daughter cells retain heterozygosity but one daughter cell has the linkage switch. (D) Meiosis / reductive segregation with crossing over, which results in LOH centromeric proximal to the crossing over, unlike (B). (E) Meiosis / reductive segregation without crossing over, which results in reciprocal parthenogenetic disomy (UPD) in the daughter cells. Note that MI with equal chromosome segregation resembles (B) and (C). In the main text, because our study focuses primarily on MI, we refer to the expected meiosis / reductive segregation during MI where sister chromatids segregate together as "reductive segregation" and the unexpected mitosis-like / equal segregation during MI where sister chromatids segregate apart as "equal segregation." [Figure 8]Figure 8A-G shows sperm and sperm precursors and their ploidy by FACS. (A) Visualization of B6 sperm. (B) Visualization of (B6 x Spret) F1 sperm. A small number of round germ cells of unknown ploidy and very few morphologically mature sperm (arrows) are observed. (C) (B6 x Spret) F1 sperm and sperm precursors isolated from the epididymis unexpectedly contain a large proportion of 2C cells. DAPI voltage of 375. (D) HEK293 / Patski mix, DAPI voltage of 350. The Patski peak (2C) is slightly shifted to the left with respect to the 2C peak in (C) due to the lower DAPI voltage. (E) (B6 x Cast) F1 sperm isolated from the epididymis are almost entirely composed of 1C cells. DAPI voltage of 375. (F) (B6×Cast)F1 sperm precursors, pre-sorting for 2C cells from dissociated testes; numerous 1C cells are still present. DAPI voltage at 375. (G) (B6×Cast)F1 sperm and sperm precursors at the FACS step (after two rounds of barcoding) during sci-L3-WGS, still consisting mostly of 1C cells. Based on the percentage of contaminating 1C nuclei among the pre-sorted 2C nuclei from (F), the percentage of tagmented nuclei that are 2C is estimated to be 18%, a 7.2-fold enrichment from the 2.5% of 2C nuclei in homogenized testes. We sorted from the 2C population (approximately 15.4% of total cells, similar to the 18% estimated in the tagmentation step). DAPI voltage at 375. [Figure 9]Figure 9A-F shows that sci-L3-WGS of interspecific hybrid mouse male germlines reveals numerous examples of non-independent equal segregation in MI. In (A), (B) and (C), the red lines indicate the crossover transpositions fitted by HMM. The centromere is located at the left-most position of each chromosome picture. (A) Example of a crossover diagram for 1C cells. The grey dots have a value of 1 for the Spret allele and a value of 0 for the B6 allele. In (B) and (C), the grey dots indicate the allele frequency of Spret at an average of 40 SNP sites. (B) Example of a LOH diagram for M2 cells with reduced segregation (see also Figure 7D). LOH is present in the centromere-proximal region of the crossover site. (C) Example of a LOH diagram for M2 cells with equal segregation (see also Figure 7B). LOH is present in the centromere-distal region of the crossover site, unlike (B). (D-F) Number of reductionally (red, pink, black) and evenly (blue, green) segregated chromosomes for each of the M2 cells. Each column represents one single M2 cell (19 chromosomes per cell, distributed as indicated by the colors). (D) Expected distribution of reduced vs even segregation based on the binomial distribution and assuming that the probability of reduced segregation, p, is equal to MLE 0.76 from the observed data. (E) Observed data in M2 cells. In rare cases (27 / 5,548 chromosomes), it was not possible to distinguish reduced vs even segregation due to sparse SNP coverage (white space at the top of the panel). Black bars indicate MI nondisjunctions (NDJ, 40 chromosomes in total) where 0 or 4 copies of chromatids were observed. It should be noted that NDJ is considered reduced segregation because sister chromatids segregate together. (F) Same as (E), but further subdivided by the number of chromosomes with and without crossovers (abbreviated as "CO"). Cells are sorted first by the number of evenly segregated chromosomes (light green and blue, descending order) and then by the number of evenly segregated chromosomes observed without crossovers (blue, descending order). [Figure 10]Figure 10A-G shows meiotic crossover and parthenogenetic chromosome distributions at the chromosome scale. (A) After normalizing for chromosome size, the number of haploid cells with at least one crossover over each chromosome is negatively correlated with chromosome size (r=-0.87, p=2e-6). (B6xSpret) crossovers are shown. For (B6xCast) crossovers, see Figure 14C. (B) Same as (A) for M2 cells (r=-0.91, p=8e-8). For (B6xCast) crossovers, see Figure 14D. (C) Distribution of the number of crossovers (CO) per chromosome per haploid cell (mean=0.62 for (B6xSpret), mean=0.58 for (B6xCast)). (D) Same as (C) for M2 cells (mean=0.92 for (B6xSpret), mean=1.03 for (B6xCast)). (E) Crossover distances for all chromosomes for chromosomes with at least two crossovers. The expected number distribution is generated by randomly placing two crossovers per chromosome. (B6 x Spret) crossovers are shown. For (B6 x Cast) crossovers, see FIG. 14E. (F) Number (top) and chromosome distribution (bottom) of UPD and LOH events in Patski cells. (G) Disaggregated mitochondrial copy numbers (normalized) for M2 cells with reductionist vs. even segregation of the majority of chromosomes. (B6 x Spret) crossovers. [Figure 11]Figure 11A-E shows that sci-L3-WGS of interspecific hybrid mouse male germlines also reveals numerous examples of non-independent equal segregation. (A-B) Number of reductively (red) and evenly (blue) segregated chromosomes for artificial "2C" cells from barcoded group 1, derived from a doublet of two random 1C cells. Each column represents one single 2C cell (19 chromosomes per cell, distributed as indicated by colors). (A) Expected distribution of reduced vs. equal segregation based on a binomial distribution and assuming the probability of equal segregation p equal to 0.5. (B) Observed data in 2C cells, consistent with the expected distribution shown in (A). (C-E) Number of reductively (red, pink, black) and evenly (blue, green) segregated chromosomes for non-1C cells from barcoded group 2, a mixture of both artificial doublets of two random 1C nuclei and real 2C secondary spermatocytes. Each column represents one single non-1C cell (19 chromosomes per cell, distributed as indicated by colors). (C) All non-1C cells from barcoded group 2. (D) Non-1C cells with only biased chromosome segregation, i.e., with at least 15 chromosomes with even or reduced segregation. Black bars indicate meiosis I nondisjunctions (NDJ, 2 out of a total of 2,185 chromosomes) where 0 or 4 copies of chromatids were observed. (E) Same as (D), but further subdivided by the number of chromosomes with and without crossovers (abbreviated as "CO"). Cells are sorted first by the number of evenly segregated chromosomes (light green and blue, descending order) and then by the number of evenly segregated chromosomes observed without crossovers (blue, descending order). [Figure 12] Figures 12A-C show fitted finite mixture models with three binomial distributions (top) compared to observed data from sci-L3-WGS of mouse male germline (bottom). For details of mixture modeling, see Example 2. (A) Mixture modeling of non-1C cells from barcode group 1 in (B6 x Cast) hybrids. (B) Mixture modeling of non-1C cells from barcode group 2 in (B6 x Cast) hybrids. (C) Mixture modeling of 2C cells from (B6 x Spret) crosses. [Figure 13]Figures 13A-I show meiotic crossovers and parthenogenetic chromosome distributions on the chromosome scale. (A) The number of crossovers (cM / Mb) normalized by chromosome size is negatively correlated with chromosome size in haploid cells (r=-0.66, p=0.002). (B6×Spret) crossovers are shown. For (B6×Cast) crossovers, see Figure 14A. (B) Same as (A) for M2 cells (r=-0.83, p=1e-5). (B6×Spret) crossovers are shown. For (B6×Cast) crossovers, see Figure 14B. (C) Distribution of crossover (CO) frequencies per chromosome per haploid cell. For number distribution, see Figure 10C. (D) Same as (C) for M2 cells. For number distribution, see Figure 10D. (E) Distances (Mb) between crossovers for chromosomes 1, 2, 12 and 13 for chromosomes with at least two crossovers. For all chromosomes, see FIG. 10E. (B6×Spret) crosses are shown. For (B6×Cast) crosses, see FIG. 14E. The expected distribution of numbers is generated by randomly placing two crossovers per chromosome. The box plot shows that (B6×Cast) crosses have stronger crossover interference than (B6×Spret) crosses (p=5e-91). (F) Histogram of parthenogenetic chromosome numbers per haploid (median=8, mean=8.1), M2 cell (median=1, mean=1.1) or other diploid / 4C (median=0, mean=0.4) cells. (B6×Spret) crosses are shown. For (B6×Cast) crosses, see FIG. 14F. (G) Parthenogenetic chromosome distribution for haploid (r=-0.87, p=2e-6), M2 cells (r=-0.75, p=2e-4) and other diploid / 4C (r=-0.68, p=0.001) cells. (B6×Spret) crosses are shown. For (B6×Cast) crosses see FIG. 14G. (H) Chromosomal distribution of reverse segregation events in (B6×Spret) (left) and (B6×Cast) (right) crosses. (I) Number of mitochondrial reads per cell normalized by read depth for haploid, M2 cells and other diploid / 4C diploid cells. (B6×Spret) crosses. [Figure 14]14A-G show chromosome distribution for meiotic crossovers and UPD (B6×Cast). (A) The number of crossovers (cM / Mb) normalized by chromosome size is negatively correlated with chromosome size in haploid cells (r=-0.65, p=0.003). (B6×Cast) crossovers. (B) Same as (A) for M2 cells (r=-0.9, p=2e-7). (B6×Cast) crossovers. (C) After normalizing for chromosome size, the number of haploid cells with at least one crossover on each chromosome is negatively correlated with chromosome size (r=-0.85, p=5e-6). (B6×Cast) crossovers. (D) Same as (C) for M2 cells (r=-0.94, p=3e-9). (B6×Cast) crossovers. (E) Crossover distance for all chromosomes for chromosomes with at least two crossovers. Expected number distributions are generated by randomly placing two crossovers per chromosome. (B6 x Cast) cross. (F) Parthenogenetic chromosome numbers per cell in haploid (median = 8, mean = 8.9) and M2 cells (median = 0, mean = 0.54). (B6 x Cast) cross. (G) Parthenogenetic chromosome distribution (correlation with chromosome size shown in brackets), haploid (r = -0.8, p = 4e-5) and M2 cells (r = -0.45, p = 0.05). (B6 x Cast) cross. [Figure 15]Figure 15A-C shows crossover breakpoint pileup profiles. (A) From top to bottom: meiotic DSB hotspots by SSDS map for B6, Cast and (B6 x Cast) F1 hybrids, crossover maps in (B6 x Spret) and (B6 x Cast) created in this study. For haploid vs. M2 cell breakdown and Spol1-oligo maps see (B) and (C). (B) From top to bottom: 1) meiotic DSB hotspots by SSDS map for (B6 x Cast) F1 hybrids, 2) haploid crossover map of (B6 x Cast), and 3) M2 cell crossover map of (B6 x Cast). (C) From top to bottom: 1) meiotic DSB hotspots by Spo11-oligo map with "symmetric" hotspots, 2) meiotic DSB hotspots by Spo11-oligo map with all hotspots: PRDM9 motifs are not considered. 3) haploid cross map in (B6 x Spret), and 4) M2 cell cross map in (B6 x Spret). [Figure 16]Figure 16A-F shows meiotic crossover hotlines and explanatory genomic features. (A) Borderline inclusion probability for features associated with crossover hotlines by BMA. The x-axis ranks models by posterior probability, where grey boxes indicate features not included in each model (vertical lines, indicating top 20 models), and the orange scale indicates the posterior probability of the models. A combined dataset from (B6 x Spret) and (B6 x Cast) crosses is shown here. For the two crosses analyzed separately, see Figure 15. (B) Distribution of sizes for breakpoint resolution (log normal distribution). Left: (B6 x Spret), median 150 kb. Right: (B6 x Cast), median 250 kb. (C-D) Location of rightmost crossovers for each chromosome. Chromosome length is indicated by SNPs (black bars) farther to the right than the red line range. (C) M2 cells. Crossovers in the (B6 x Cast) (left) crossovers favor the centromere-distal ends of chromosomes, whereas crossovers in the (B6 x Spret) crossovers (right) favor the central regions of each chromosome arm. After accounting for variability between chromosomes, we estimate that crossovers in the (B6 x Spret) crossovers are on average 5.5 Mb more proximal to the centromere. Similar except for 1C cells, see Figure 20A. (D) Comparing 1C and M2 cells, the (B6 x Spret) crossovers. After accounting for variability between chromosomes, we estimate that crossovers in the M2 cells (right) are on average 9.4 Mb more proximal to the centromere than 1C (left) in the (B6 x Spret) crossovers. The same trend is observed to a lesser extent in the (B6 x Cast) crossover (see Figure 20B). (E) An AUC of 0.73 quantifies the expected accuracy in predicting whether a region drawn from the mouse genome comes from a B6×Spret crossover zone or an equal number of randomly sampled zones. Left: 76 total features. Right: subset of 25 features from BMA with MIP>0.5. (F) An AUC of 0.85 quantifies the expected accuracy in predicting whether a region drawn from the mouse genome comes from a B6×Cast crossover zone or an equal number of randomly sampled zones. Left: 69 total features. Right: subset of 25 features from BMA with MIP>0.5. [Figure 17]17A-B show the borderline inclusion probability for features associated with cross hot varnish by BMA. The x-axis ranks the models by posterior probability. (A) (B6 x Cast) cross. (B) (B6 x Spret) cross. [Figure 18] Figure 18 shows the correlation matrix for both crossover events and genomic features in the (B6 x Cast) cross. Here, we show all possible pairwise correlations between various crossover pileup tracks and genomic features, calculated in 100 kb windows. The crossover pileup tracks are the first five columns or rows ("event" prefix; red text), while the rest are the same genomic features used in the modeling (blue text). The crossover pileup tracks with the suffixes "hp_m2", "hp", "m2", "mt" and "me" are from haploid and M2 cells, haploid, M2 cells, M2 cells with biased even segregation and M2 cells with biased reduced segregation, respectively. Blue boxes indicate positive correlations and red boxes indicate negative correlations. Features are ordered by hierarchical clustering. Open ovals highlight the features "telomere" and "quantile_75_100", which show different trends in the two crosses as described in the text. [Figure 19] Figure 19 shows the correlation matrix for both crossover events and genomic features in the (B6 x Spret) cross. Same format as described in Figure 18 legend. [Figure 20]Figures 20A-E show the location of the crossovers on the right edge of each chromosome. (A) Haploid cells. In both crossovers, crossovers favor the centromere-distal ends of the chromosomes. (B) Comparing haploid and M2 cells (B6 x Cast cross). After accounting for variability between chromosomes, we estimate that crossovers in M2 cells are, on average, 5.2 Mb proximal to the centromere than haploids in the (B6 x Cast) cross. (C) Comparing M2 cells with biased chromosome segregation. After accounting for variability between chromosomes, we estimate that crossovers in M2 cells with biased even segregation are, on average, 13.7 Mb distal to the centromere than those of M2 cells with biased reduced segregation in the (B6 x Cast) cross. (D) Same as (C) for the (B6 x Spret) cross. Crossovers are, on average, 8.7 Mb more distal to the centromere. (E) Model of the effect of crossover location on proper chromosome segregation. Crossovers closer to the centromere (in the middle two quartiles rather than the last quartile) can promote reduction segregation by having stronger arm cohesion; however, crossovers near the ends of the chromosome arms can promote MII segregation by having stronger CEN cohesion. [Figure 21] Figure 21 shows a principal component analysis of features distinguishing crossover hotspots in B6 x Spret crosses. Note that "chr3_bp (breakpoint)" and "chr1_upc (parthenogenetic chromosome)" are representative of features included for all chromosomes. 44 out of 115 total features are shown. 36 other chromosomal breakpoints and 35 other features other than the UPC feature that were omitted are not shown due to lack of obvious trends. [Figure 22] Figure 22 shows a principal component analysis of features distinguishing crossover hotspots of B6 x Cast crosses. Note that "chr3_bp (breakpoint)" and "chr1_upc (parthenogenetic chromosome)" are representative of features included for all chromosomes. 19 of 108 total features are shown. 53 other features are not shown due to lack of obvious trends, except for the omitted 36 other chromosomal breakpoint and UPC features. [Figure 23] 23 shows a model of the relationship between meiotic crossover and chromosome missegregation. "MI": meiosis I, "CEN": centromere (oval or round circle), "IH": interhomolog. The following detailed description of exemplary embodiments of the present disclosure can be best understood when read in conjunction with the following figures. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0050] The schematic drawings are not necessarily to scale. Like numbers used in the figures refer to like components, steps, etc. However, it will be understood that the use of a number to refer to a component in a given figure is not intended to limit the component in another figure labeled with the same number. Moreover, the use of different numbers to refer to a component is not intended to indicate that the differently numbered component is not necessarily the same as or similar to the other numbered component. For example, the disclosure provides the following: [Section 1] 1. A method for preparing a sequencing library comprising nucleic acids from a plurality of single nuclei or cells, comprising: Providing a plurality of isolated nuclei or cells in a first plurality of compartments, Each compartment contains an isolated nucleus or a subset of cells, the nucleus or cell contains the nucleic acid fragment; introducing a linear amplification agent into said cell or nucleus; amplifying the nucleic acid fragment by linear amplification; processing each subset of nuclei or cells to generate indexed nuclei or cells; said processing comprising adding a first compartment-specific index sequence to a nucleic acid fragment present in said isolated nucleus or cell to provide an indexed nucleic acid present in the isolated nucleus or cell; The processing includes a step of ligation, primer extension, hybridization, amplification or transposition; combining said indexed nuclei or cells to generate pooled indexed nuclei or cells, thereby generating a sequencing library from a plurality of said nuclei or cells; The method includes: [Section 2] 2. The method of claim 1, wherein the amplification occurs prior to the treatment. [Section 3] 2. The method of claim 1, wherein the treatment occurs before the amplification. [Section 4] 1. A method for preparing a sequencing library comprising nucleic acids from a plurality of single nuclei or cells, comprising: Providing a plurality of isolated nuclei or cells, the nucleus or cell contains the nucleic acid fragment; introducing a linear amplification agent into said isolated nuclei or cells; Distributing the isolated nuclei or cells into a first plurality of compartments, each compartment containing an isolated nucleus or a subset of cells; amplifying the nucleic acid fragment by linear amplification; Processing each subset of isolated nuclei or cells to generate indexed nuclei or cells, said processing comprising adding a first compartment-specific index sequence to a nucleic acid fragment present in said isolated nucleus or cell to provide an indexed nucleic acid present in the isolated nucleus or cell; The processing includes a step of ligation, primer extension, amplification or transposition; combining said indexed nuclei to generate pooled indexed nuclei or cells, thereby generating a sequencing library from a plurality of said nuclei or cells; The method includes: [Section 5] 1. A method for preparing a sequencing library comprising nucleic acids from a plurality of single nuclei or cells, comprising: Providing a plurality of isolated nuclei or cells in a first plurality of compartments, Each compartment contains an isolated nucleus or a subset of cells, the nucleus or cell contains the nucleic acid fragment; processing each subset of nuclei or cells to generate indexed nuclei or cells; The processing includes adding to the nucleic acid fragments present in the isolated nuclei or cells (i) a first compartment-specific index sequence to provide an indexed nucleic acid present in the isolated nuclei or cells, and (ii) adding a nucleotide sequence recognized by a linear amplification agent; The processing includes a step of ligation, primer extension, hybridization, amplification or transposition; introducing a linear amplification agent into said cell or nucleus; amplifying the nucleic acid fragment by linear amplification; combining said indexed nuclei or cells to generate pooled indexed nuclei or cells, thereby generating a sequencing library from a plurality of said nuclei or cells; The method includes: [Section 6] 6. The method of any one of paragraphs 1, 4 or 5, wherein the linear amplification agent comprises a phage RNA polymerase or a linear amplification primer. [Section 7] 7. The method of claim 6, wherein the nucleic acid fragment comprises a T7 promoter and the phage RNA polymerase comprises T7 RNA polymerase. [Section 8] 7. The method of claim 6, wherein introducing the linear amplification agent comprises adding the linear amplification agent to the nucleic acid fragment present in the isolated nucleus or cell. [Section 9] The method of paragraph 1 or 5, further comprising exposing a plurality of the isolated nuclei or cells of each compartment to a predetermined condition. [Section 10] 10. The method of paragraph 9, further comprising isolating nuclei from a plurality of said cells after said exposing. [Section 11] The method of claim 4, further comprising exposing a plurality of the isolated nuclei or cells to a predetermined condition. [Section 12] 6. The method of any one of paragraphs 1, 4 or 5, further comprising subjecting the isolated nuclei to conditions that generate nucleosome-depleted nuclei while maintaining the integrity of the isolated nuclei. [Section 13] The process, contacting each subset with a transposome complex, the transposome complexes in each compartment comprising a first index sequence that is different from first index sequences in other compartments; fragmenting the nucleic acids in said subset into a plurality of nucleic acids and incorporating said first index sequence into at least one strand of said nucleic acids to generate said indexed nuclei or cells comprising said indexed nucleic acids; The method according to any one of items 1, 4 or 5, comprising: [Section 14] The process, contacting each subset with a reverse transcriptase and a primer that anneals to RNA molecules in said isolated nuclei, said primer in each compartment comprising a first index sequence that differs from a first index sequence in other compartments, to generate indexed said nuclei or cells comprising said indexed nucleic acids. The method according to any one of items 1, 4 or 5, comprising: [Section 15] 15. The method of claim 14, wherein said contacting further comprises a target-specific primer that anneals to a specific nucleotide sequence. [Section 16] 6. The method of any one of paragraphs 1, 4, or 5, wherein the process of adding the first compartment-specific index sequence comprises a two-step process of adding a nucleotide sequence comprising a ubiquitous sequence to the nucleic acid fragments and then adding the first compartment-specific index sequence to the nucleic acid fragments. [Section 17] 17. The method of claim 16, wherein said adding comprises a transposome complex comprising said ubiquitous sequence. [Section 18] 6. The method of any one of paragraphs 1, 4 or 5, wherein the processing comprises adding a first index to DNA nucleic acid present in the isolated nucleus or cell, adding a first index to RNA nucleic acid present in the isolated nucleus or cell, or a combination thereof. [Section 19] adding said first index sequence to an RNA nucleic acid, contacting each subset with a reverse transcriptase and a primer that anneals to RNA molecules in said isolated nuclei or cells; the primers in each compartment contain a first compartment-specific index sequence, generating indexed nuclei or cells that contain the indexed nucleic acid. Item 19. The method according to item 18, comprising: [Section 20] adding said first index sequence to a DNA nucleic acid, contacting each subset with a transposome complex, the transposome complex in each compartment comprises the first compartment-specific index sequence; fragmenting the nucleic acids in said subset into a plurality of nucleic acids and incorporating said first compartment-specific index sequence into at least one strand of said nucleic acids to generate said indexed nuclei or cells comprising said indexed nucleic acids; Item 19. The method according to item 18, comprising: [Section 21] 20. The method of claim 19, wherein the first index sequence added to the DNA nucleic acid and the first index sequence added to the RNA nucleic acid in each compartment are identical. [Section 22] 20. The method of claim 19, wherein the first index sequence added to the DNA nucleic acid and the first index sequence added to the RNA nucleic acid in each compartment are not identical. [Section 23] 6. The method of any one of paragraphs 1, 4 or 5, further comprising exponential amplification of said nucleic acid fragments, said exponential amplification comprising target specific primers that anneal to specific nucleotide sequences. [Section 24] After said combining, distributing said pooled indexed subsets of nuclei or cells into a second plurality of compartments; introducing a second compartment-specific index sequence into the indexed nucleic acid to generate a doubly-indexed nucleus or cell comprising the doubly-indexed nucleic acid, The introducing step may include ligation, primer extension, amplification, or transposition. The method of any one of paragraphs 1, 4 or 5, further comprising: [Section 25] combining said doubly indexed nuclei to generate pooled doubly indexed nuclei or cells; distributing said pooled dual-indexed subsets of nuclei or cells into a third plurality of compartments; introducing a third compartment-specific index sequence into the indexed nucleic acid to generate a triply indexed nucleus or cell comprising the triply indexed nucleic acid, The introducing step may include ligation, primer extension, amplification, or transposition. 25. The method of claim 24, further comprising: [Section 26] 6. The method of any one of paragraphs 1, 4 or 5, further comprising processing the indexed nuclei or cells for methylation analysis to generate nucleic acid fragments suitable for methylation analysis. [Section 27] 6. The method of any one of paragraphs 1, 4 or 5, further comprising subjecting the indexed nuclei or cells to proximity ligation to generate nucleic acid fragments suitable for chromatin conformation analysis. [Section 28] Item 6. The method of any one of items 1, 4 or 5, further comprising amplifying the nucleic acid fragments of the sequencing library to generate DNA nanoballs. [Section 29] 29. The method of any one of paragraphs 1 to 28, wherein the compartment comprises a well or a droplet. [Section 30] 30. The method according to any one of paragraphs 1 to 29, wherein each compartment of the first plurality of compartments comprises 50 to 100,000,000 nuclei or cells. [Section 31] 30. The method of any one of paragraphs 1 to 29, wherein each compartment of the second plurality of compartments comprises 50 to 100,000,000 nuclei or cells. [Section 32] Providing a surface comprising a plurality of amplification sites, the amplification site comprising at least two populations of attached single-stranded capture oligonucleotides having free 3' ends; contacting said surface containing amplification sites with said indexed fragments under conditions suitable to generate a plurality of amplification sites, each comprising a clonal population of amplicons from an individual fragment containing a plurality of indexes; Item 32. The method according to any one of items 1 to 31, further comprising: [Section 33] 1. A method for preparing a sequencing library comprising nucleic acids from a plurality of single cells, comprising: (a) providing isolated nuclei from a plurality of cells; (b) subjecting the isolated nuclei to a chemical treatment that maintains the integrity of the isolated nuclei while generating nucleosome-depleted nuclei; (c) distributing a subset of the nucleosome-depleted nuclei into a first plurality of compartments and contacting each subset with a transposome complex, wherein the transposome complex in each compartment comprises a transposase and a first index sequence that is different from the first index sequence in the other compartments; (d) fragmenting nucleic acid in said subset of nucleosome-depleted nuclei into a plurality of nucleic acid fragments and incorporating said first index sequence into at least one strand of said nucleic acid fragments to generate indexed nuclei comprising indexed nucleic acid fragments, wherein said indexed nucleic acid fragments remain attached to said transposase; (d) combining the indexed nuclei to generate pooled indexed nuclei; (e) distributing a subset of the pooled indexed nuclei into a second plurality of compartments and contacting each subset with a hairpin ligation duplex under conditions suitable for ligation of the hairpin ligation duplex to one or both ends of an indexed nucleic acid fragment to provide a doubly-indexed nucleic acid fragment, wherein the hairpin ligation duplex comprises a second index sequence that differs from the second index sequences in the other compartments; (f) combining the doubly indexed nuclei to generate pooled indexed nuclei; (g) distributing the pooled doubly indexed subset of nuclei into a third plurality of compartments; (h) lysing the dual-indexed nuclei; (i) processing the dual-indexed nuclear fragments to include a third index sequence that is different from a third index sequence in another compartment; (j) combining the triple index fragments, thereby generating a sequencing library comprising the entire genome nucleic acid from a plurality of the single cells; The method includes:

[0051] The methods provided herein can be used to generate single-cell combinatorial indexing (sci) sequencing libraries of multiple single cells or nuclei, including, for example, whole genome (sci-WGS), transcriptome (sci-RNA), genome and transcriptome (sci-DNA / RNA) and / or methylome (sci-MET) co-assays. In one embodiment, the method can be used for targeted sequencing of one or more specific regions of interest. For example, primers hybridizing to specific regions (e.g., coding regions, non-coding regions, etc.), guide RNAs or nucleotide sequences inserted by guide RNAs can be used to selectively enrich for targeted sequences. In one embodiment, information for individual gene edits, DNA edits or markers for editing, gene signatures, perturbations and / or functional readouts (RNA, DNA, protein or combinations) from cells or nuclei can be collected and analyzed (Perturb-seq). In other embodiments, the methods can be used to assess chromatin accessibility (sci-ATAC), chromatin higher order structure (Hi-C) and other single-cell combinatorial indexing methods.

[0052] The method includes providing isolated nuclei or cells, distributing a subset of nuclei or cells into compartments, treating the nuclei or cells so that they contain nucleic acid fragments, adding compartment-specific indexes to the nucleic acid fragments, and amplifying the nucleic acid fragments by linear amplification. These steps can occur in different orders and can be combined in different ways. Three embodiments are shown in Figures 1A and 1B. In one embodiment, the method includes providing a distributed subset of isolated nuclei or cells that contain nucleic acid fragments (Figure 1A, Block 1, and Figure 1B, Block 1). As shown in Figure 1AB, following amplification of the nucleic acid fragments by linear amplification (Figure 1A, Block 2), an index is added to the amplified nucleic acid fragments (Figure 1A, Block 3). As shown in Figure 1B, the nucleic acid fragments in the distributed nuclei or cells contain an index and the nucleic acid fragments are amplified by linear amplification (Figure 1B, Block 2). Described herein are steps of providing isolated nuclei or cells, distributing a subset of the isolated nuclei or cells, processing the isolated nuclei or cells to contain nucleic acid fragments, adding compartment-specific indexes, and amplifying the nucleic acid fragments by linear amplification.

[0053] Providing isolated nuclei or cells The methods provided herein include providing a cell or a nucleus isolated from a plurality of cells. The cells and nuclei can be from any sample, such as any organism(s), and from any cell type or any tissue of the organism(s). In one embodiment, the cells can be reproductive cells, such as sperm cells or egg cells. In one embodiment, the tissue can be reproductive tissue, such as epididymis. In one embodiment, the cells or nuclei can be from cancer or diseased tissue. The method can further include dissociating the cells and / or isolating the nuclei. Methods for isolating nuclei from cells are known and routine to those skilled in the art. The number of nuclei or cells can be at least two. The upper limit depends on the practical limitations of the equipment (e.g., multi-well plates) used in other steps of the methods described herein. The number of nuclei or cells that can be used is not limiting and can number in the billions. For example, in one embodiment, the number of nuclei or cells may be 100,000,000 or less, 10,000,000 or less, 1,000,000,000 or less, 100,000,000 or less, 10,000,000 or less, 1,000,000 or less, 100,000 or less, 10,000 or less, 1,000 or less, 500 or less, or 50 or less. One or more samples can be provided. For example, a sample can be one cell type or tissue from one organism. To identify and then combine samples, the indexing method described herein can be used to separately index multiple samples, such as different cell types from one organism, one cell type or tissue from two or more organisms, or different cell types or tissues from two or more organisms, with a first index. Those skilled in the art will appreciate that in some embodiments, the nucleic acid molecule within each nucleus is a genomic DNA molecule that represents the entire genetic complement of the organism (also referred to as the entire genome of the organism) and includes both intron and exon sequences, as well as non-coding regulatory sequences such as promoter and enhancer sequences.

[0054] Isolation of nuclei can be accomplished by incubating cells in cell lysis buffer for at least 1-20 minutes, e.g., 5, 10 or 15 minutes. Optionally, cells can be exposed to an external force to aid in lysis, such as movement with a pipette. An example of a cell lysis buffer includes 10 mM Tris-HCl, pH 7.4, 10 mM NaCl, 3 mM MgCl2, 0.1% IGEPAL CA-630, and 1% SUPERase In RNase Inhibitor. Those skilled in the art will recognize that these levels of components can be modified slightly without reducing the usefulness of the cell lysis buffer to isolate nuclei. Those skilled in the art will recognize that RNase inhibitors, BSA, and / or detergents may be beneficial in buffers used for isolation of nuclei, and other additives can be added to the buffer for other downstream single-cell combinatorial indexing applications.

[0055] In one embodiment, nuclei are isolated from individual cells that are adherent or in suspension. Methods for isolating nuclei from individual cells are known to those skilled in the art. In one embodiment, nuclei are isolated from cells present in tissue. The method for obtaining isolated nuclei generally includes preparing tissue and isolating nuclei from the prepared tissue. In one embodiment, all steps are carried out on ice.

[0056] Preparation of tissue may involve flash freezing the tissue in liquid nitrogen, followed by mincing or application of blunt force to reduce the size of the tissue to fragments of 1 mm or less in diameter. Optionally, cold proteases and / or other enzymes to disrupt intercellular communication may be used. Mincing may be accomplished with a blade to cut the tissue into small pieces. Application of blunt force may be accomplished by crushing the tissue with a hammer or similar object, and the resulting composition of the crushed tissue is called a powder.

[0057] Traditional tissue nucleus extraction techniques typically incubate tissues with tissue-specific enzymes (e.g., trypsin) at high temperatures (e.g., 37°C) for 30 minutes to several hours, and then lyse the cells with cell lysis buffer for nuclei extraction. The nucleus isolation method described herein and in U.S. Provisional Patent Application No. 62 / 680,259 has several advantages: (1) no artificial enzymes are introduced, and all steps are performed on ice. This reduces potential perturbations to the cell state (e.g., transcriptome state, chromatin state, or methylation state). (2) It has been validated across most tissue types, including disease samples such as brain, lung, kidney, spleen, heart, cerebellum, and tumor tissues. Compared to traditional tissue nucleus extraction techniques that use different enzymes for different tissue types, the new technique can potentially reduce bias when comparing cell states from different tissues. (3) The method also reduces costs and increases efficiency by removing the enzyme treatment step. (4) Compared to other nuclear extraction techniques (e.g., Dounce tissue grinder), the present technique is more robust for different tissue types (e.g., the Dounce method requires optimizing the Dounce cycle for different tissues) and allows for high-throughput processing of large sample pieces (e.g., the Dounce method is limited by the size of the grinder).

[0058] The isolated nuclei or cells may contain nucleosomes, may be nucleosome-free, or may be subjected to conditions that deplete the nucleosomes to produce nucleosome-depleted nuclei. Nucleosome-depleted nuclei are useful in methods for determining the DNA sequence of the entire genome of a cell, or a fraction thereof.

[0059] In one embodiment, the conditions used for nucleosome depletion maintain the integrity of isolated nuclei. Generally, the nucleosome depletion method is used on a pellet or suspension of single cells, thus in embodiments where adherent cell cultures or tissues are used as a source of cells, the source is processed to obtain a pellet or suspension of single cells.

[0060] Methods for depleting nucleosomes are known and routine, including, but not limited to, enzymatic and chemical treatments. In one embodiment, the conditions for depleting nucleosomes include chemical treatment with a chaotropic agent capable of disrupting nucleic acid-protein interactions. An example of a useful chaotropic agent includes, but is not limited to, 3,5-lithium diiodosalicylic acid. The conditions for using 3,5-lithium diiodosalicylic acid include adding it to a pellet of cells and incubating on ice.

[0061] In a preferred embodiment, the conditions include chemical treatment with a detergent capable of disrupting nucleic acid-protein interactions. Examples of useful detergents include, but are not limited to, sodium dodecyl sulfate (SDS). Conditions for using SDS include adding it to a pellet of cells and incubating at an elevated temperature, such as 42°C, followed by adding a non-ionic detergent, such as Triton™ X-100, and incubating at an elevated temperature, such as 42°C.

[0062] In some embodiments, when using detergents such as SDS, the nuclei are exposed to a crosslinking agent prior to nucleosome depletion (WO2018 / 018008). In one embodiment, the nuclei are exposed to the crosslinking agent while in the cells, and in another embodiment, the isolated nuclei are exposed to the crosslinking agent. Useful examples of crosslinking agents include, but are not limited to, formaldehyde (Hoffman et al., 2015, J. Biol. Chem., 290:26404-26411). Treating the cells with formaldehyde can include adding formaldehyde to the cell suspension and incubating at room temperature. In one embodiment, after formaldehyde treatment, the nuclei can be exposed to a non-ionic, non-denaturing detergent such as glycine and Igepal®.

[0063] During the process of depleting nucleosomes in the isolated nuclei, the integrity of the isolated nuclei is maintained. Whether the nuclei remain intact after exposure to conditions for depleting nucleosomes can be determined by visualizing the state of the nuclei by routine methods such as phase contrast imaging. In one embodiment, the number of nuclei that are intact after nucleosome depletion can be 1 to 1,000, 1,000 to 10,000, 10,000 to 100,000, 100,000 to 1,000,000, 1,000,000 to 10,000,000, or 10,000,000 to 100,000,000.

[0064] The nucleus or cell manipulation, including the providing, pooling and distributing steps described herein, can include the use of nucleus buffer. An example of nucleus buffer includes 10 mM Tris-HCl, pH 7.4, 10 mM NaCl, 3 mM MgCl2, 1% SUPERase In RNase Inhibitor (20 U / μL, Ambion) and 1% BSA (20 mg / ml, NEB). Those skilled in the art will recognize that these levels of components can be slightly changed without reducing the usefulness of the nucleus buffer to suspend nuclei. Those skilled in the art will also recognize that various components can be substituted without reducing the usefulness of the nucleus buffer to suspend nuclei.

[0065] In one embodiment, the cells (including the cells from which the nuclei are isolated) are exposed to different predetermined conditions. For example, a subset of cells can be exposed to different predetermined conditions. The different conditions can include, for example, different culture conditions (e.g., different media, different environmental conditions), different doses of drugs, different drugs or drug combinations. Drugs are described herein. The nuclei or cells of each subset of cells and / or one or more samples are indexed with one or more index sequences, pooled, and analyzed by a massively multiplexed single nucleus or single cell sequencing method. Single nucleus transcriptome sequencing (US Provisional Patent Application No. 62 / 680,259 and Gunderson et al. (WO2016 / 130704)), single nucleus whole genome sequencing (US Patent Application Publication No. 2018 / 0023119), or single nucleus sequencing of transposon accessible chromatin (US Pat. No. 10,059,989), sci-HiC (Ramani et al., Nature Methods, 2017, 14:263-266), DRUG-seq (Ye et al., Nature Commun., 9, article number 4307), Perturb-seq (Dixit et al., Cell, 2016, 167(7):1853-1866.e17), or any combination of analytes from DNA, RNA and protein, such as sci-CAR (Cao et al., Science, Virtually any single nucleus or single cell sequencing method can be used, including but not limited to, methods such as those described in J. Immunol. 2018, 361(6409):1380-1385. Droplet-based single cell analysis can also be applied after initial split-and-pool indexing (examples include the 10X Genome Research Chromium™ system or Biorad ddseq system), including the use of the index as a sample index. Nuclear hashing is used to demultiplex and identify individual cells or nuclei from different conditions.

[0066] In one embodiment, each subset of cells is exposed to a drug or perturbation. The drug can be virtually anything that causes a change in the cell. For example, the drug can change the transcriptome of the cell, change the chromatin structure of the cell, change the activity of a protein in the cell, change the DNA of the cell, change the DNA editing of the cell, or cause other changes. Examples of drugs include, but are not limited to, compounds, such as proteins (including antibodies), non-ribosomal proteins, polyketides, organic molecules (including organic molecules of 900 daltons or less), inorganic molecules, RNA or RNAi molecules, carbohydrates, glycoproteins, nucleic acids, or combinations thereof. In one embodiment, the drug causes a genetic perturbation, such as a DNA editing protein and / or a guide RNA, such as CRISPR or Talen. In one embodiment, the drug is a therapeutic agent. In one embodiment, the cell can be a wild-type cell, and in another embodiment, the cell can be genetically modified to include a genetic perturbation, such as a gene knock-in or gene knock-out (Szlachta et al., Nat Commun., 2018, 9:4275). Subsets of cells can be exposed to the same drug, but different variables can be varied across the compartments of the multi-well device, allowing multiple variables to be tested in a single experiment. For example, different dosages, different exposure times and different cell types can be tested in a single plate. In one embodiment, the cells can express a protein with a known activity, and the effect of the drug on the activity can be evaluated in different conditions. The use of index sequences to label nucleic acid fragments allows for subsequent identification of nucleic acids originating from a specific subset of nuclei or cells, for example, one well of a multi-well plate.

[0067] Distributing the subset The methods provided herein include distributing a subset of nuclei, e.g., nucleosome-depleted nuclei or cells, into a plurality of compartments. The methods can include a plurality of partitioning steps, in which a population of isolated nuclei or cells (also referred to herein as a pool) is divided into subsets. Generally, partitioning of a subset of isolated nuclei or cells from a pool into a plurality of compartments occurs prior to the addition of indexes to the nucleic acid fragments present in the subset of isolated nuclei or cells. Thus, the methods include at least one "split and pool" step of taking the pooled isolated nuclei or cells and distributing them, where the number of "split and pool" steps can depend on the number of different indexes added to the nucleic acid fragments. After indexing, the subsets can be pooled, divided into subsets, indexed, and pooled again as necessary until a sufficient number of indexes have been added to the nucleic acid fragments.

[0068] The number of nuclei or cells present in a subset, and thus in each compartment, may be at least 1. In one embodiment, the number of nuclei or cells present in a subset is 100,000,000 or less, 10,000,000 or less, 1,000,000 or less, 100,000 or less, 10,000 or less, 4,000 or less, 3,000 or less, 2,000 or less, or 1,000 or less, 500 or less, or 50 or less. In one embodiment, the number of nuclei or cells present in a subset may be 1 to 1,000, 1,000 to 10,000, 10,000 to 100,000, 100,000 to 1,000,000, 1,000,000 to 10,000,000, or 10,000,000 to 100,000,000. In one embodiment, the number of nuclei or cells present in each subset is approximately equal. The number of nuclei present in each subset and therefore in each compartment is based in part on the need to reduce index collisions, which is the presence of two nuclei with the same transposase index present in the same compartment in this step of the method. Methods for distributing nuclei or cells into subsets are known and routine to those skilled in the art. Examples include, but are not limited to, fluorescence activated cell sorting (FACS) cytometry and simple dilution. Optionally, nuclei of different ploidy can be gated and enriched by staining, for example DAPI (4',6-diamidino-2-phenylindole) staining.

[0069] The number of compartments in the partitioning step (and subsequent addition of index) can depend on the format used. For example, the number of compartments can be 2 to 96 compartments (when a 96-well plate is used), 2 to 384 compartments (when a 384-well plate is used), or 2 to 1536 compartments (when a 1536-well plate is used). In one embodiment, each compartment can be a droplet. When the type of compartment used is a droplet containing two or more nuclei or cells, any number of droplets can be used, for example at least 10,000, at least 100,000, at least 1,000,000 or at least 10,000,000 droplets. In one embodiment, the number of compartments is 24.

[0070] Processing to generate nucleic acid fragments In one embodiment, processing of isolated nuclei or cells can be used to fragment DNA nucleic acids, such as chromosomes and / or plasmids, in isolated nuclei or cells into nucleic acid fragments. Processing is generally necessary when the target nucleic acid to be sequenced is derived from DNA present in the nuclei or cells. However, in some embodiments, processing is appropriate when the target nucleic acid to be sequenced is derived from RNA (e.g., mRNA and / or non-coding RNA) present in the nuclei or cells, since RNA molecules often do not need to be fragmented. Processing the nucleic acid in the nuclei or cells typically adds a nucleotide sequence to one or both ends of the nucleic acid fragments generated by the processing, and the nucleotide sequence can and typically does include one or more ubiquitous sequences. The ubiquitous sequences can be used as "landing pads" in subsequent steps to anneal nucleotide sequences that can be used as primers for the addition of another nucleotide sequence, such as an index, to the nucleic acid fragments by subsequent steps of ligation, primer extension, or amplification. The nucleotide sequence of such a primer may include an index sequence. Processing the nucleic acid in the nuclei or cells can add one or more unique molecular identifiers to one or both ends of the nucleic acid fragments generated by the processing.

[0071] There are several points in the method where processing of nucleic acids into nucleic acid fragments can occur. For example, in one embodiment, isolated nuclei or cells can be processed prior to partitioning of isolated nuclei or cell subsets. In such embodiments, processing generally includes the addition of ubiquitous sequences and / or ubiquitous molecular identifiers to the nucleic acid fragments, rather than compartment-specific indexes, since adding a compartment-specific index generally serves no purpose when all isolated nuclei or cells are combined. In another embodiment, isolated nuclei or cells can be processed after partitioning of subsets into different compartments (see, e.g., Figures 1A and 1B). In one aspect of this embodiment, processing does not add an index (Figure 1A, block 1), and in another aspect of this embodiment, processing can include the addition of a compartment-specific index (Figure 1B, block 1). Processing at any point in the method can include the addition of ubiquitous sequences and / or ubiquitous molecular identifiers to one or both ends of the nucleic acid fragments.

[0072] Various methods are known for processing nucleic acids in a nucleus or cell into nucleic acid fragments. Examples include CRISPR and Talen-like enzymes, as well as enzymes that unwind DNA (e.g., helicases) that can create single-stranded regions to which DNA fragments can hybridize and initiate extension or amplification. For example, helicase-based amplification can be used (Vincent et al., 2004, EMBO Rep., 5(8):795-800). In one embodiment, extension or amplification is initiated with random primers. In one embodiment, a transposome complex is used. A transposome complex is a transposase that is bound to a transposase recognition site and can insert the transposase recognition site into a target nucleic acid in a nucleus in a manner sometimes called "tagmentation". In some such insertion events, one strand of the transposase recognition site can be transferred to the target nucleic acid. Such a strand is called the "transfer strand". In one embodiment, the transposome complex includes a dimeric transposase with two subunits and two discontinuous transposon sequences. In another embodiment, the transposase comprises a dimeric transposase having two subunits and a contiguous transposon sequence. In one embodiment, the 5' end of one or both strands of the transposase recognition site can be phosphorylated.

[0073] Some embodiments may include the use of a hyperactive Tn5 transposase and a Tn5-type transposase recognition site (Goryshin and Reznikoff, J. Biol. Chem., 273:7367 (1998)), or a MuA transposase and a Mu transposase recognition site including R1 and R2 end sequences (Mizuuchi, K., Cell, 35: 785, 1983; Savilahti, H, et al., EMBO J., 14: 4893, 1995). A Tn5 Mosaic End (ME) sequence may also be used, as optimized by one of skill in the art.

[0074] More examples of transposition systems that can be used in certain embodiments of the compositions and methods provided herein include Staphylococcus aureus Tn552 (Colegio et al., J. Bacteriol., 183: 2384-8, 2001; Kirby C et al., Mol. Microbiol., 43: 173-86, 2002), Ty1 (Devine & Boeke, Nucleic Acids Res., 22: 3765-72, 1994 and WO 95 / 23875), transposon Tn7 (Craig, NL, Science. 271: 1512, 1996; Craig, NL, Review in: Curr Top Microbiol Immunol., 204:27-48, 1996), Tn / O and IS10 (Kleckner N, et al., Curr Top Microbiol Immunol., 204:49-82, 1996), Mariner transposase (Lampe DJ, et al., EMBO J., 15: 5470-9, 1996), Tc1 (Plasterk RH, Curr. Topics Microbiol. Immunol., 204: 125-43, 1996), P elements (Gloor, GB, Methods Mol. Biol., 260: 97-114, 2004), Tn3 (Ichikawa & Ohtsubo, J Biol. Chem. 265:18829-32, 1990), bacterial insertion sequences (Ohtsubo & Sekine, Curr. Top. Microbiol. Immunol. 204:1-26, 1996), retroviruses (Brown, et al., Proc Natl Acad Sci USA, 86:2525-9, 1989), and yeast retrotransposons (Boeke & Corces, Annu Rev Microbiol. 43:403-34, 1989).More examples include IS5, Tn10, Tn903, IS911 and engineered versions of the transposase family of enzymes (Zhang et al., (2009) PLoS Genet. 5:e1000689. Epub 2009 Oct 16; Wilson C. et al (2007) J. Microbiol. Methods 71:332-5).

[0075] Other examples of integrases that can be used in the methods and compositions provided herein include retroviral integrases and integrase recognition sequences for such retroviral integrases, such as integrases from HIV-1, HIV-2, SIV, PFV-1, RSV.

[0076] Transposon sequences useful in the methods and compositions described herein are provided in US Patent Application Publication No. 2012 / 0208705, US Patent Application Publication No. 2012 / 0208724 and International Application No. WO2012 / 061832. In some embodiments, the transposon sequence comprises a first transposase recognition site and a second transposase recognition site. In embodiments in which a transposome complex is used to introduce an index sequence, the index sequence can be present between the transposase recognition sites or within the transposon.

[0077] Some transposome complexes useful herein include a transposase with two transposon sequences.In some such embodiments, the two transposon sequences are not linked to each other, in other words, the transposon sequences are not contiguous with each other.Examples of such transposomes are known in the art (see, for example, US Patent Application Publication No. 2010 / 0120098).

[0078] In some embodiments, the transposome complex comprises a transposon sequence nucleic acid that binds to two transposase subunits to form a "looped complex" or "looped transposome". In one example, the transposome comprises a dimeric transposase and a transposon sequence. The looped complex can ensure that the transposon is inserted into the target DNA while maintaining the instruction information of the original target DNA and without fragmenting the target DNA. As will be appreciated, the looped structure can insert a desired nucleic acid sequence, such as an index, into the target nucleic acid while maintaining the physical continuity of the target nucleic acid. In some embodiments, the transposon sequence of the looped transposome complex can include a fragmentation site so that the transposon sequence can be fragmented to create a transposome complex containing two transposon sequences. Such a transposome complex is beneficial to ensure that the neighboring target DNA fragments into which the transposon is inserted receive a code combination that can be unambiguously assembled at a later stage of the assay.

[0079] In one embodiment, fragmenting the nucleic acid is accomplished using a fragmentation site present in the nucleic acid. Generally, the fragmentation site is introduced into the target nucleic acid using a transposome complex. In one embodiment, after the nucleic acid is fragmented, the transposase remains attached to the nucleic acid fragments, so that nucleic acid fragments originating from the same genomic DNA molecule remain physically linked (Adey et al., 2014, Genome Res., 24:2041-2049). For example, a looped transposome complex can include a fragmentation site. The fragmentation site can be used to sever the physical association, but not the informational association, between index sequences inserted into the target nucleic acid. The cleavage can be biochemical, chemical, or by other means. In some embodiments, the fragmentation site can include nucleotides or nucleotide sequences that can be fragmented by various means. Examples of fragmentation sites include, but are not limited to, restriction endonuclease sites, at least one ribonucleotide cleavable by an RNase, a nucleotide analogue cleavable in the presence of certain chemical agents, a diol bond cleavable by treatment with periodate, a disulfide group cleavable with a chemical reducing agent, a cleavable moiety amenable to photochemical cleavage, and a peptide cleavable by a peptidase enzyme or other suitable means (see, e.g., U.S. Patent Application Publication No. 2012 / 0208705, U.S. Patent Application Publication No. 2012 / 0208724, and WO2012 / 061832).

[0080] The transposome complex may contain at least one index sequence, which may be referred to as a transposase index. The index sequence is present as part of the transposon sequence. In one embodiment, the index sequence may be present on the transposed strand, the strand of the transposase recognition site that has been transferred to the target nucleic acid.

[0081] The transposome complex may include at least one nucleotide sequence that can be used by a linear amplification mediator. Examples of such nucleotide sequences include, but are not limited to, an RNA polymerase when the nucleic acid fragment includes a phage promoter, such as T7 RNA polymerase for a T7 promoter, and a linear amplification primer. Examples of linear amplification primers include a single primer for use in PCR-type amplification or a linear amplification mediator. Other embodiments of nucleotide sequences that can be used by a linear amplification mediator are sequences that are recognized by strand-displacing polymerases. The mediator may contain a nicking site to initiate replication. In some cases, the nicking site is regenerated for further amplification.

[0082] Adding Compartment-Specific Indexes The index sequence, also called tag or barcode, is useful as a marker characteristic of the compartment in which a particular nucleic acid was present.The index is thus a nucleic acid sequence tag attached to each of the target nucleic acids present in a particular compartment, the presence of which indicates or is used to identify the compartment in which the isolated nuclei or cell population was present at a particular stage of the method.The addition of the index to the nucleic acid fragment is achieved in a subset of the isolated nuclei or cells that are distributed in different compartments.

[0083] Index sequences can be any suitable number of nucleotides in length, e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 or more nucleotides. Four-nucleotide tags offer the possibility of multiplexing 256 samples on the same array, while six-base tags allow processing 4096 samples on the same array.

[0084] In one embodiment, the addition of the index is accomplished during processing of the nucleic acid to the nucleic acid fragment. For example, a transposome complex containing the index can be used. In other embodiments, the index is added after processing produces nucleic acid fragments containing nucleotide sequences at one or both ends. Methods for adding the index include, but are not limited to, ligation, extension (including extension using reverse transcriptase), hybridization, adsorption, specific or non-specific interaction of primers, or amplification. The nucleotide sequences added to one or both ends of the nucleic acid fragment can also include one or more ubiquitous sequences and / or unique molecular identifiers. The ubiquitous sequences can be used as "landing pads" in subsequent steps to anneal nucleotide sequences that can be used as primers for the addition of another nucleotide sequence, such as another index and / or another ubiquitous sequence, to the nucleic acid fragment.

[0085] For example, in embodiments involving the use of nucleic acid fragments derived from mRNA, various methods can be used to add an index to the mRNA in one or two steps. For example, the index can be added using the type of methods used to generate cDNA. A primer with a 3'-terminal poly-T sequence can be annealed to the mRNA molecule and extended using reverse transcriptase. Exposing isolated nuclei or cells to these components under conditions suitable for reverse transcription results in one-step addition of the index, resulting in a population of indexed nuclei or cells, where each nucleus or cell contains an indexed nucleic acid fragment. Alternatively, the primer with the poly-T sequence contains a ubiquitous sequence instead of the index, and the index is added by subsequent steps of ligation, primer extension, and amplification. The indexed nucleic acid fragment can, and typically does, contain an index sequence that indicates a specific compartment on the synthesized strand.

[0086] In embodiments involving the use of nucleic acid fragments derived from non-coding RNA, various methods can be used to add an index to the non-coding RNA in one or two steps. For example, a first primer and a template switch primer containing a random sequence can be used to add an index, where either primer can contain an index. A reverse transcriptase with terminal transferase activity can be used that results in the addition of non-template nucleotides to the 3' end of the synthesized strand, and the template switch primer contains nucleotides that anneal with the non-template nucleotides added by the reverse transcriptase. One example of a useful reverse transcriptase is Moloney Murine Leukemia Virus Reverse Transcriptase. In certain embodiments, SMARTer™ Reagents (Cat. No. 634926) available from Takara Bio USA, Inc. are used for the use of template switching to add an index to non-coding RNA, and optionally, mRNA. Alternatively, the first primer and / or the template switch primer can contain a ubiquitous sequence instead of an index, and the index is added by a subsequent step of ligation, primer extension, or amplification. Indexed nucleic acid fragments can, and typically do, include an index sequence that indicates a particular compartment on the synthesized strand. Other embodiments include 5' or 3' profiling of RNA or full-length RNA profiling.

[0087] Other methods can be used for the addition of indexes to nucleic acid fragments, and how the index is added is not limiting. For example, in one embodiment, the incorporation of an index sequence includes ligating a primer to one or both ends of a nucleic acid fragment. Ligation of a ligation primer can be aided by the presence of a ubiquitous sequence at the end of a nucleic acid fragment. A non-limiting example of a primer is a hairpin ligation duplex. A ligation duplex can be ligated to one or preferably both ends of a nucleic acid fragment. In one embodiment, a primer such as a hairpin ligation duplex can contain a nucleotide sequence that is recognized by a linear amplification agent. An example of a hairpin adapter containing such nucleotides is described in Example 1, FIG. 2. An assay scheme such as that described in Example 1, which introduces an amplification agent that only requires successful ligation at one of the two ends of a barcoded molecule to generate an amplification product of that molecule, is desirable because it has the advantage of increasing the efficiency of template conversion. For example, if a single ligation event has an efficiency of 50%, this modification gives a 75% success rate in the ligation step of amplifying the molecule instead of 25% (Example 1, Figure 2).

[0088] In another embodiment, the incorporation of the index sequence involves the use of a single-stranded nucleic acid fragment and the synthesis of a second DNA strand. In one embodiment, the second DNA strand is generated using a primer that includes a sequence complementary to the nucleotides present at the end of the single-stranded nucleic acid fragment.

[0089] In another embodiment, incorporation of the index occurs in one, two, three or more rounds of split and pool barcoding resulting in single-, double-, triple- or multiple-indexed single-cell libraries.

[0090] In another embodiment, the incorporation of index and nucleotide sequences that the amplification agent can use is designed unidirectionally, allowing the preparation of targeted single-cell sequencing libraries (see Example 1, FIG. 3b).

[0091] Linear Amplification of Nucleic Acid Fragments The methods provided herein include linear amplification of nucleic acid fragments. Most amplification methods are PCR-based and therefore suffer from exponential amplification bias. Linear amplification used herein can reduce or eliminate exponential amplification bias, leading to better uniformity and reduced sequence errors. In all single-cell genomic research methods utilizing whole genome amplification, the amplification products are placed in a compartment (e.g., a well or a droplet) and a barcode is directly or indirectly attached to the amplification product. In this way, there is only a single cell per compartment, limiting throughput and increasing costs. A unique aspect of the present invention is that multiple single-cell libraries can be amplified without exponential amplification bias in a single compartment. Libraries from single cells can be assigned based on the unique barcode or barcodes for each unique single cell.

[0092] In one embodiment, linear amplification is achieved by adding a phage promoter to one or both ends of the nucleic acid fragment. When placed upstream of the nucleic acid fragment, the phage promoter can be used to promote transcription using the corresponding phage RNA polymerase by in vitro transcription to generate single-stranded RNA. The RNA copies generated from the DNA template cannot serve as templates for further amplification. Thus, all copies are directly derived from the original DNA template, and exponential amplification is avoided. In one embodiment, subsequent steps can then include reverse transcription of the RNA copies to obtain single-stranded DNA, and second strand synthesis to convert the single-stranded DNA copies into double-stranded molecules. Second strand synthesis typically requires the use of a primer, which can be used to introduce one or more of the index, ubiquitous sequence and / or ubiquitous molecular identifier.

[0093] Other methods of linear amplification can be used. For example, PCR amplification can be used with one primer or with two primers, one extra. In some embodiments, linear PCR can be used for amplification of flanking sequences adjacent to the transposon insertion site (Xianbo et al. AMB Express, 2017, 7:195). Concatenated linear amplification (Reyes et al., Clin. Chem., 2001, 47(1):31-40), linear extension and linear extension and ligation, strand displacement amplification (SDA) (Walker et al., Nucl. Acids Res., 1992, 20(7): 1691-1696) and rolling circle amplification (Ali et al., Chem. Soc. Rev., 2014, 43:3324-3341) can also be used in some embodiments. In some embodiments, indexes, ubiquitous sequences and / or unique molecular identifiers can be added to the nucleic acid fragments during linear amplification.

[0094] Generally, linear amplification involves introducing a linear amplification agent into isolated nuclei or cells. Examples of linear amplification agents include RNA polymerase when the nucleic acid fragment contains a phage promoter, such as T7 RNA polymerase for a T7 promoter, and linear amplification primers. Examples of linear amplification primers include a single primer or linear amplification agent for use in PCR-type amplification. Another embodiment of an amplification agent is a strand-displacing polymerase that recognizes a nucleotide sequence. The agent can contain a nicking site to initiate replication. In some cases, the nicking site is regenerated for further amplification. The agent can contain a unique barcode or primer that allows for replication of the barcode during amplification or labeling of the amplification product.

[0095] Addition of ubiquitous sequences for immobilization In one embodiment, the addition of nucleotides during the processing and / or indexing step adds ubiquitous sequences that are beneficial in immobilizing and sequencing the fragments. In another embodiment, the indexed nucleic acid fragments can be further processed to add ubiquitous sequences that are beneficial in immobilizing and sequencing the nucleic acid fragments. Those skilled in the art will recognize that in embodiments in which the compartments are droplets, the sequences for immobilizing the nucleic acid fragments are appropriate. In one embodiment, the incorporation of ubiquitous sequences that are beneficial in immobilizing and sequencing the fragments includes ligating identical ubiquitous adapters (also called "mismatch adapters," the general features of which are described in Gormley et al., U.S. Pat. No. 7,741,463, and Bignell et al., U.S. Pat. No. 8,053,192) to the 5' and 3' ends of the indexed nucleic acid fragments. In one embodiment, the ubiquitous adapters include all sequences required for sequencing, including sequences for immobilizing the indexed nucleic acid fragments on the array.

[0096] In one embodiment, blunt-end ligation can be used. In another embodiment, the nucleic acid fragments are prepared with a single overhanging nucleotide, for example, by the activity of certain types of DNA polymerases, such as Taq polymerase or Klenow exo minus polymerase, that have non-template-dependent terminal transferase activity that adds a single deoxynucleotide, such as deoxyadenosine (A), to the 3' end of the indexed nucleic acid fragment. In some cases, the overhanging nucleotide is multi-base. Such enzymes can be used to add a single nucleotide "A" to the blunt 3' end of each strand of the nucleic acid fragment. Thus, an "A" can be added to the 3' end of each strand of the double-stranded target fragment by reaction with Taq or Klenow exo minus polymerase, and the additional sequence added to each end of the nucleic acid fragment can include a matching "T" overhang present on the 3' end of each region of the double-stranded nucleic acid that is added. This end modification also prevents self-ligation of the nucleic acid, so there is a bias toward the formation of indexed nucleic acid fragments flanked by the sequences added in this embodiment.

[0097] In another embodiment, when the ubiquitous adapters ligated to the indexed nucleic acid fragments do not contain all the sequences required for sequencing, an amplification step such as PCR can be used to further modify the ubiquitous adapters present in each indexed nucleic acid fragment before immobilization and sequencing.For example, an initial primer extension reaction can be performed using a ubiquitous anchor sequence that is complementary to the ubiquitous sequence present in the indexed nucleic acid fragments, where an extension product is formed that is complementary to both strands of each individual indexed nucleic acid fragment.Generally, PCR adds additional ubiquitous sequences, such as a ubiquitous capture sequence.

[0098] After ubiquitous adapters are added, either by a single-step method of ligating ubiquitous adapters that contain all the sequences required for sequencing, or by a two-step method of ligating and then amplifying the ubiquitous adapters to further modify them, the final index fragments contain a ubiquitous capture sequence and an anchor sequence. The result of adding a ubiquitous adapter to each end is a plurality of indexed nucleic acid fragments, or a library thereof.

[0099] The resulting indexed nucleic acid fragments collectively provide a library of nucleic acids that can be immobilized and subsequently sequenced. The term library, also referred to herein as a sequencing library, refers to a collection of nucleic acid fragments from a single nucleus or cell that contain known ubiquitous sequences at their 3' and 5' ends.

[0100] The indexed nucleic acid fragments can be subjected to conditions that select for a given size range, for example, 150 to 400 nucleotides in length, for example, 150 to 300 nucleotides. The resulting indexed nucleic acid fragments can be pooled and subjected to a clean-up process to enhance the purity to DNA molecules by removing at least a portion of the unincorporated ubiquitous adaptors or primers. Any suitable clean-up process can be used, for example, electrophoresis, size-exclusion chromatography, etc. In some embodiments, solid-phase reversibly immobilized paramagnetic beads can be used to separate the desired DNA molecules from the unattached ubiquitous adaptors or primers and select the nucleic acids based on size. Solid-phase reversibly immobilized paramagnetic beads are commercially available from Beckman Coulter (Agencourt AMPure XP), Thermofisher (MagJet), Omega Biotek (Mag-Bind), Promega Beads (Promega) and Kapa Biosystems (Kapa Pure Beads).

[0101] A non-limiting exemplary embodiment of the present disclosure is shown in Figure 2 and described in Example 1. In this embodiment, the method includes providing isolated nuclei from a plurality of cells (Figure 2, block 22). The isolated nuclei can be free of nucleosomes or can be subjected to conditions that deplete the nucleosomes of the nuclei to generate nucleosome-depleted nuclei (Figure 2, block 23).

[0102] In this embodiment, the method includes distributing a subset of nucleosome-depleted nuclei into a first plurality of compartments (FIG. 2, block 24). The number of compartments in the first distribution step (FIG. 2, block 24) can depend on the format used. In one embodiment, the number of compartments is 24.

[0103] Each compartment contains a transposome complex. The transposome complex can be added to each compartment before, after, or at the same time as the subset of nuclei is added to the compartment. The transposome complex contains at least one index sequence and at least one ubiquitous sequence. The ubiquitous sequence present as part of the transposome complex can be referred to as a spacer sequence. The spacer sequence is present as part of the transposon sequence. In one embodiment, the spacer sequence can be present on the transferred strand, the strand of the transposase recognition site transferred to the target nucleic acid. The spacer sequence is useful as a site for annealing with a complementary sequence. For example, the spacer sequence can be a ubiquitous primer or the complement of a ubiquitous primer. The spacer sequence of the transposome complex can be the same in each compartment. In one embodiment, the index ("bc1") and spacer ("sp1") are as shown in Figure 1 of Example 1. 4 They are arranged in the orientation shown in A and are present in the protrusions.

[0104] The method includes generating indexed nuclei (FIG. 2, block 25). In one embodiment, generating indexed nuclei includes processing nucleic acids present in a subset of nucleosome-depleted nuclei (e.g., nucleic acids present in each compartment) into a plurality of nucleic acid fragments. In one embodiment, after the nucleic acids are fragmented, the transposase remains attached to the nucleic acid fragments, such that nucleic acid fragments originating from the same genomic DNA molecule remain physically linked (Adey et al., 2014, Genome Res., 24:2041-2049). The result of fragmentation is a population of indexed nuclei, where each nucleus contains an indexed nucleic acid fragment. The index sequence of the transposome complex is different in each compartment, and thus the indexed nucleic acid fragments can, and typically do, include an index sequence on at least one strand that indicates a particular compartment. An example of indexed nucleic acid fragments is shown in FIG. 1 of Example 1. 4 Shown in box A.

[0105] The indexed nuclei from the multiple compartments can be combined (FIG. 2, block 26). A subset of these combined indexed nuclei are then distributed to a second multiple compartment. The number of nuclei present in the subsets, and therefore in each compartment, is based in part on the desire to reduce index collisions, that is, the presence of two nuclei with the same transposase index present in the same compartment at this step of the method. In one embodiment, the number of nuclei present in each subset is approximately equal.

[0106] Partitioning of the nuclei into subsets is followed by incorporation of a second index sequence into the indexed nucleic acid fragments in each compartment to generate dual index fragments. This results in further indexing of the indexed nucleic acid fragments (FIG. 2, block 27). In embodiments in which the cells are crosslinked by a crosslinking agent, the transposase attached to the indexed nucleic acid fragments can be dissociated from the indexed nucleic acid fragments. A detergent can be used to dissociate the transposase, and in one embodiment, the detergent is sodium dodecyl sulfate (SDS).

[0107] In one embodiment, the incorporation of the second index sequence comprises ligating a hairpin ligation duplex to the indexed nucleic acid fragment in each compartment. The ligation duplex can be ligated to one or preferably both ends of the double indexed nucleic acid fragment. In one embodiment, the ligation duplex comprises five elements: 1) the reverse complement of the first spacer sequence (e.g., see Figure 1 of Example 1); 4 1) a second strand synthesis (SSS) primer region that starts with GGG to enhance T7 transcription (see FIG. 1B); 2) the reverse complement of the second round barcode; 3) the reverse complement of the second strand synthesis (SSS) primer; 4) a T7 promoter, which is preferably the hairpin loop region; 5) a second strand synthesis (SSS) primer region that starts with GGG to enhance T7 transcription (see FIG. 1B). 4 5) the second spacer sequence in B, "sp2"; and 6) the second barcode, second index sequence (see FIG. 1). 4 The second index sequence is unique to each compartment, where the indexed nuclei to be distributed were placed after the first index was added by tagmentation (Figure 2, block 27).

[0108] The indexed nuclei from the multiple compartments can be combined (FIG. 2, block 28). A subset of these combined indexed nuclei are then distributed to a third multiple compartment. The number of nuclei present in the subsets, and therefore in each compartment, is based in part on the desire to reduce index collisions, that is, the presence of two nuclei with the same transposase index present in the same compartment at this step of the method. In one embodiment, between 100 and 300 cells are distributed to each well. In one embodiment, up to 300 cells are distributed to each well. In one embodiment, the number of nuclei present in each subset is approximately equal.

[0109] Partitioning of the doubly indexed nuclei into subsets is followed by lysis and further manipulation (Figure 2, block 29). Methods for lysis of nuclei are known and routine to those skilled in the art. Further manipulations include, but are not limited to, gap extension, in vitro transcription (IVT) and reverse transcription.

[0110] Gap extension converts the hairpin T7 promoter structure into a duplex (Figure 1 in Example 1). 4 C) For gap extension, polymerases with strand displacement activity are generally used. Polymerases with this activity, such as Bst polymerase, are available.

[0111] IVT generates linear amplified single-stranded RNA molecules downstream of a T7 promoter (Figure 1 in Example 1). 4 D). Methods for IVT are known and routine.

[0112] Reverse transcription can occur by one of two pathways (see Figure 1). 4E). The ligation reaction described herein results in two types of nucleic acid fragments: a nucleic acid fragment with a ligation duplex at both ends and a nucleic acid fragment with a ligation duplex at one end. If ligation is successful at both ends, reverse transcription can be primed by the self-looping reverse transcription primers, which are inherited from the looped ligation duplex. If ligation is successful at only one end, reverse transcription is primed by an additional RNA reverse transcription primer added in excess.

[0113] Lysis of the nuclei and processing of the nucleic acid fragments is followed by incorporation of a third index sequence into the doubly indexed nucleic acid fragments in each compartment to generate triply indexed fragments, where the third index sequence in each compartment is different from the first and second index sequences in the other compartments, and where the third index sequence in each compartment is different from the third index sequence in the other compartments. This results in further indexing of the indexed nucleic acid fragments prior to fixation and sequencing (Figure 2, block 30; Figure 2 of Example 1). 4 F). The third index can be incorporated by synthesis of the second DNA strand. In one embodiment, the second DNA strand is generated using a primer that contains a sequence complementary to the nucleotides present at the end of the doubly indexed nucleic acid fragment. For example, the primer can contain a second spacer sequence (sp2) that anneals with the reverse complement of the second spacer sequence (Figure 1 in Example 1). 4 F). The primers were the third index (Figure 1). 4 F) and other unique molecular identifiers (UMIs). The resulting double-stranded DNA can be purified using routine methods.

[0114] A plurality of triple-indexed fragments can be prepared for sequencing. After pooling the triple-indexed fragments, they are enriched, typically by immobilization and / or amplification, prior to sequencing (Figure 2, block 31).

[0115] Preparation of immobilized samples for sequencing A plurality of indexed fragments can be prepared for sequencing. For example, in an embodiment in which a library of triple-indexed fragments is generated, the triple-indexed fragments are typically enriched by immobilization and / or amplification prior to sequencing (FIG. 2, block 21). Methods for attaching indexed fragments from one or more sources to a substrate are known in the art. In one embodiment, the indexed fragments are enriched using a plurality of capture oligonucleotides with specificity for the indexed fragments, and the capture oligonucleotides can be immobilized to the surface of a solid substrate. For example, the capture oligonucleotides can include a first member of a ubiquitous binding pair, where a second member of the binding pair is immobilized to the surface of a solid substrate. Similarly, methods for amplifying immobilized dual-indexed fragments include, but are not limited to, cross-linking amplification and kinetic exclusion. Methods for immobilization and amplification prior to sequencing are described, for example, in Bignell et al. (U.S. Patent No. 8,053,192), Gunderson et al. (WO2016 / 130704), Shen et al. (U.S. Patent No. 8,895,249) and Pipenburg et al. (U.S. Patent No. 9,309,502).

[0116] The pooled samples can be immobilized in preparation for sequencing. Sequencing can be performed as an array of single molecules or can be amplified prior to sequencing. Amplification can be performed using one or more immobilized primers. The immobilized primer(s) can be, for example, on a flat surface or on a lawn on the pool of beads. The pool of beads can be isolated in an emulsion with a single bead in each "compartment" of the emulsion. At a concentration of only one template per "compartment", only a single template is amplified on each bead.

[0117] The term "solid-phase amplification" as used herein refers to any nucleic acid amplification reaction that is carried out on or in conjunction with a solid support, so that all or part of the amplification products are immobilized on the solid support as they are formed. In particular, the term encompasses solid-phase polymerase chain reaction (solid-phase PCR) and solid-phase isothermal amplification, which are reactions similar to standard solution-phase amplification, except that one or both of the forward and reverse amplification primers are immobilized on a solid support. Solid-phase PCR encompasses systems such as colony formation in emulsions where one primer is anchored to a bead and the other is in free solution, and solid-phase gel matrices where one primer is anchored to a surface and one is in free solution.

[0118] In some embodiments, the solid support comprises a patterned surface. "Patterned surface" refers to an arrangement of different regions in or on an exposed layer of a solid support. For example, one or more of the regions may be features in which one or more amplification primers are present. The features may be separated by intermediate regions in which no amplification primers are present. In some embodiments, the pattern may be an xy format of features in rows and columns. In some embodiments, the pattern may be a repeating arrangement of features and / or intermediate regions. In some embodiments, the pattern may be a random arrangement of features and / or intermediate regions. Exemplary patterned surfaces that can be used in the methods and compositions provided herein are described in U.S. Pat. Nos. 8,778,848, 8,778,849 and 9,079,148, and U.S. Patent Application Publication No. 2014 / 0243224.

[0119] In some embodiments, the solid support comprises an array of wells or surface depressions. This can be fabricated using a variety of techniques, including but not limited to photolithography, stamping techniques, molding techniques and microetching techniques, as are generally known in the art. As one of ordinary skill in the art will appreciate, the technique used will depend on the composition and shape of the array substrate.

[0120] The features of the patterned surface may be wells in an array of wells (e.g., microwells or nanowells) on glass, silicon, plastic, or other suitable solid support with a patterned covalently attached gel, such as poly(N-(5-azidoacetamidylpentyl)acrylamide-co-acrylamide) (PAZAM, see, e.g., U.S. Patent Application Publication No. 2013 / 184796, WO2016 / 066586, and WO2015 / 002813). This process creates a gel pad used for sequencing that can be stable over multiple cycles of sequencing runs. Covalent attachment of the polymer to the wells is effective in maintaining the gel in the structured features over the life of the structured substrate during various uses. However, in many embodiments, the gel need not be covalently attached to the wells. For example, in some conditions, silane-free acrylamide (SFA, see, e.g., U.S. Patent No. 8,563,477), which is not covalently attached to any part of the structured substrate, can be used as the gel material.

[0121] In certain embodiments, a structured substrate can be made by patterning a solid support with wells (e.g., microwells or nanowells), coating the patterned support with a gel material (e.g., PAZAM, SFA or chemically modified variants thereof, such as the azidolylated version of SFA (azido-SFA)), and polishing the gel-coated support, e.g., by chemical or mechanical polishing, thereby retaining the gel in the wells but removing or inactivating substantially all of the gel from the intermediate regions of the surface of the structured substrate between the wells. Primer nucleic acid can be attached to the gel material. A solution of indexed fragments can then be contacted with the polished substrate such that individual indexed fragments seed individual wells through interaction with primers attached to the gel material. However, due to the absence or inactivity of the gel material, the target nucleic acid does not occupy the intermediate regions. Amplification of the indexed fragments is restricted to the wells, as the absence or inactivity of the gel in the intermediate regions prevents outward migration of the growing nucleic acid colonies. This process can be conveniently fabricated and is scalable, utilizing conventional micro- or nanofabrication methods.

[0122] Although the present disclosure encompasses "solid-phase" amplification methods in which only one amplification primer is immobilized (the other primer is usually in free solution), in one embodiment, it is preferred that a solid support is provided with immobilized forward and reverse primers. In practice, since the amplification process requires an excess of primers to sustain amplification, there are "multiple" identical forward primers and / or "multiple" identical reverse primers immobilized on the solid support. Reference herein to forward and reverse primers should be construed as encompassing "multiple" such primers, unless the context indicates otherwise.

[0123] As the skilled reader will understand, any given amplification reaction requires at least one type of forward primer and at least one type of reverse primer that are specific to the template to be amplified. However, in certain embodiments, the forward and reverse primers can contain template-specific portions of the same sequence and can have completely identical nucleotide sequences and structures (including any non-nucleotide modifications). In other words, it is possible to carry out solid-phase amplification using only one type of primer, and such a single primer method is encompassed within the scope of the present disclosure. Other embodiments can use forward and reverse primers that contain the same template-specific sequence but differ in some other structural features. For example, one type of primer can contain non-nucleotide modifications that are not present in the other.

[0124] In all embodiments of the present disclosure, primers for solid-phase amplification are preferably immobilized by single-point covalent attachment to the solid support at or near the 5' end of the primer, leaving the template-specific portion of the primer free to anneal to its cognate template and free for primer extension of the 3' hydroxyl group. Any suitable covalent attachment means known in the art can be used for this purpose. The attachment chemistry selected will depend on the nature of the solid support and any derivatization or functionalization applied to it. The primers themselves can contain moieties that may be non-nucleotide chemical modifications to facilitate attachment. In certain embodiments, the primers can contain a sulfur-containing nucleophile at the 5' end, such as a phosphorothioate or thiophosphate. In the case of solid-supported polyacrylamide hydrogels, this nucleophile binds to the bromoacetamide groups present in the hydrogel. A more specific means of attaching primers and templates to a solid support is by 5' phosphorothioate attachment to a hydrogel composed of polymerized acrylamide and N-(5-bromoacetamidylpentyl)acrylamide (BRAPA), as described in WO 05 / 065814.

[0125] Certain embodiments of the present disclosure may utilize solid supports that include an inert substrate or matrix (e.g., glass slides, polymeric beads, etc.) that are "functionalized" by application of a layer or coating of an intermediate material that includes reactive groups that allow for covalent attachment to a biomolecule, such as a polynucleotide. Examples of such supports include, but are not limited to, polyacrylamide hydrogels supported on an inert substrate such as glass. In such embodiments, the biomolecule (e.g., polynucleotide) may be covalently attached directly to the intermediate material (e.g., hydrogel), although the intermediate material itself may be non-covalently attached to the substrate or matrix (e.g., glass substrate). The term "covalent attachment to a solid support" should be interpreted to encompass this type of arrangement.

[0126] The pooled samples can be amplified on beads, each containing forward and reverse amplification primers. In a particular embodiment, the library of indexed fragments is used to prepare clustered arrays of nucleic acid colonies, similar to those described in US Patent Application Publication No. 2005 / 0100900, US Patent No. 7,115,400, WO 00 / 18957 and WO 98 / 44151, by solid-phase amplification, more particularly solid-phase isothermal amplification. The terms "cluster" and "colony" are used interchangeably herein to refer to individual sites on a solid support that contain a plurality of identical immobilized nucleic acid strands and a plurality of identical immobilized complementary nucleic acid strands. The term "clustered array" refers to an array formed from such clusters or colonies. In this context, the term "array" should not be understood to require an orderly arrangement of clusters.

[0127] The terms "solid phase" or "surface" are used to mean a flat array in which primers are attached to a flat surface, such as a glass, silica or plastic microscope slide or similar flow cell device; beads in which one or two primers are attached to the beads and the beads are amplified; or an array of beads on a surface after the beads have been amplified.

[0128] Clustered arrays can be prepared using the process of thermocycling described in WO 98 / 44151, or a process in which the temperature is kept constant and cycles of extension and denaturation are carried out by changing reagents. Such isothermal amplification methods are described in patent application WO 02 / 46456 and US Patent Application Publication No. 2008 / 0009420. Due to the lower temperatures beneficial in isothermal processes, this is particularly preferred in some embodiments.

[0129] It is understood that any of the amplification methods described herein or generally known in the art can be used with ubiquitous or target-specific primers to amplify immobilized DNA fragments.As described in U.S. Pat. No. 8,003,354, suitable methods for amplification include, but are not limited to, polymerase chain reaction (PCR), strand displacement amplification (SDA), transcription-mediated amplification (TMA) and nucleic acid sequence-based amplification (NASBA).The above amplification methods can be used to amplify one or more nucleic acids of interest.For example, PCR, including multiplex PCR, SDA, TMA, NASBA, etc., can be utilized to amplify immobilized DNA fragments.In some embodiments, primers specifically directed to the polynucleotide of interest are included in the amplification reaction.

[0130] Other suitable methods for amplifying polynucleotides can include oligonucleotide extension and ligation, rolling circle amplification (RCA) (Lizardi et al., Nat. Genet. 19:225-232 (1998)) and oligonucleotide ligation assay (OLA) (see generally U.S. Pat. Nos. 7,582,420, 5,185,243, 5,679,524 and 5,573,907; EP 0320308; EP 0336731; EP 0439182; WO 90 / 01069; WO 89 / 12696; and WO 89 / 09835). It is understood that these amplification methods can be designed to amplify immobilized DNA fragments. For example, in some embodiments, the amplification method can include a ligation probe amplification or oligonucleotide ligation assay (OLA) reaction containing a primer specifically directed to the nucleic acid of interest. In some embodiments, the amplification method can include a primer extension-ligation reaction containing a primer specifically directed to the nucleic acid of interest. As a non-limiting example of primer extension and ligation primers that can be specifically designed to amplify the nucleic acid of interest, the amplification can include the primers used for the GoldenGate assay (Illumina, Inc., San Diego, CA) exemplified in U.S. Patent Nos. 7,582,420 and 7,611,869.

[0131] DNA nanoballs can also be used in combination with the methods and compositions described herein.The method of making and using DNA nanoballs for genome sequencing can be found in, for example, U.S. Patent and Application Publication Nos. 7,910,354, 2009 / 0264299, 2009 / 0011943, 2009 / 0005252, 2009 / 0155781, 2009 / 0118488, and as described in, for example, Drmanac et al., 2010, Science 327(5961): 78-81. Briefly, following genomic library DNA fragmentation, adapters are ligated to the fragments, the adapter-ligated fragments are circularized by ligation with circular ligase, and rolling circle amplification is performed (as described in Lizardi et al., 1998. Nat. Genet. 19:225-232 and US Patent Application Publication No. 2007 / 0099208). The extended concatemeric structure of the amplicons promotes coiling, thereby creating compact DNA nanoballs. The DNA nanoballs can be captured on a substrate to create an array that is preferably ordered or patterned, such that the distance between each nanoball is maintained, thereby allowing sequencing of individual DNA nanoballs. In some embodiments, such as those used by Complete Genome Research (Mountain View, Calif.), successive rounds of adapter ligation, amplification, and digestion are performed prior to circularization to generate head-to-tail constructs with several genomic DNA fragments separated by adapter sequences.

[0132] Exemplary isothermal amplification methods that can be used in the methods of the present disclosure include, but are not limited to, multiple displacement amplification (MDA), as exemplified by, for example, Dean et al., Proc. Natl. Acad. Sci. USA 99:5261-66 (2002), or isothermal strand displacement nucleic acid amplification, as exemplified by, for example, U.S. Patent No. 6,214,587.Other non-PCR-based methods that can be used in the present disclosure include strand displacement amplification (SDA), as described, for example, in Walker et al., Molecular Methods for Virus Detection, Academic Press, Inc., 1995; U.S. Patent Nos. 5,455,166 and 5,130,238, and Walker et al., Nucl. Acids Res. 20:1691-96 (1992), or hyperbranched strand displacement amplification, as described, for example, in Lage et al., Genome Res. 13:294-307 (2003). Isothermal amplification methods can be used, for example, with strand-displacing Phi 29 polymerase or large fragments of Bst DNA polymerase, 5'->3' exo-, for random-primed amplification of genomic DNA. The use of these polymerases takes advantage of their high processivity and strand-displacing activity. High processivity allows the polymerase to generate fragments that are 10-20 kb in length. As shown above, polymerases with low processivity and strand-displacing activity, such as Klenow polymerase, can be used to generate smaller fragments in isothermal conditions. Further description of the amplification reaction, conditions and components are set forth in detail in the disclosure of U.S. Pat. No. 7,670,810.

[0133] Another polynucleotide amplification method that is useful in the present disclosure is tagged PCR, which uses a population of two-domain primers with a constant 5' region followed by a random 3' region, as described in, for example, Grothues et al. Nucleic Acids Res. 21(5):1321-2 (1993).The first few rounds of amplification are carried out to allow multiple initiations on heat-denatured DNA based on individual hybridization from randomly synthesized 3' regions.Due to the nature of the 3' region, initiation sites are considered to be random throughout the genome.Unbound primers can then be removed, and further replication can occur using primers that are complementary to the constant 5' region.

[0134] In some embodiments, isothermal amplification can be performed using kinetic exclusion amplification (KEA), also known as exclusion amplification (ExAmp). The nucleic acid library of the present disclosure can be made using a method that includes reacting amplification reagents to generate a plurality of amplification sites, each containing a substantially clonal population of amplicons from an individual target nucleic acid seeded at the site. In some embodiments, the amplification reaction proceeds until a sufficient number of amplicons are generated to fill the capacity of each amplification site. In this manner, filling a previously seeded site to capacity inhibits target nucleic acids from landing and amplifying at the site, thereby generating a clonal population of amplicons at the site. In some embodiments, apparent clonality can be achieved even if an amplification site is not filled to capacity before a second target nucleic acid arrives at the site. Under some conditions, amplification of a first target nucleic acid can proceed to a point where a sufficient number of copies are generated to effectively outpace or overwhelm the generation of copies from a second target nucleic acid transported to the site. For example, in one embodiment using a bridge amplification process on circular features less than 500 nm in diameter, it was determined that after 14 cycles of exponential amplification for a first target nucleic acid, contamination with a second target nucleic acid at the same site produces an insufficient number of contaminating amplicons to adversely affect sequence-by-synthesis analysis on an Illumina sequencing platform.

[0135] In some embodiments, the amplification sites in the array may, but need not, be completely clonal. Rather, for some applications, individual amplification sites may be predominantly populated with amplicons from a first indexed fragment, and may also have low levels of contaminating amplicons from a second target nucleic acid. An array may have one or more amplification sites with low levels of contaminating amplicons, so long as that level of contamination does not unacceptably affect subsequent use of the array. For example, when the array is used in a detection application, an acceptable level of contamination would be a level that does not unacceptably affect the signal-to-noise or resolution of the detection technique. Thus, apparent clonality is generally relevant to a particular use or application of an array made by the methods set forth herein. Exemplary levels of contamination that may be tolerated in individual amplification sites for a particular application include, but are not limited to, at most 0.1%, 0.5%, 1%, 5%, 10% or 25% contaminating amplicons. An array may include one or more amplification sites with these exemplary levels of contaminating amplicons. For example, up to 5%, 10%, 25%, 50%, 75% or even 100% of the amplification sites in an array may have some contaminating amplicons. It will be understood that in arrays or other collections of sites, at least 50%, 75%, 80%, 85%, 90%, 95% or 99% or more of the sites may be clonal or clearly clonal.

[0136] In some embodiments, kinetic exclusion can occur when a process occurs at a rate that is fast enough to effectively exclude the occurrence of another event or process.For example, take the creation of a nucleic acid array, where the array sites are randomly seeded with triple-indexed fragments from a solution, and copies of the triple-indexed fragments are generated in an amplification process that fills each of the seeded sites to capacity.The kinetic exclusion method of the present disclosure allows the seeding and amplification processes to proceed simultaneously under conditions where the amplification rate exceeds the seeding rate.In this way, the relatively fast rate at which copies are created at the site seeded by the first target nucleic acid effectively excludes the second nucleic acid from seeding at that site for amplification.The kinetic exclusion amplification method can be carried out as described in detail in the disclosure of US Patent Application Publication No. 2013 / 0338042.

[0137] Kinetic exclusion can take advantage of a relatively slow rate to initiate amplification (e.g., a slow rate of making the first copy of the indexed fragment) relative to a relatively fast rate to make subsequent copies of the triply indexed fragment (or of the first copy of the indexed fragment). In the example of the previous paragraph, kinetic exclusion occurs due to a relatively slow rate of seeding of the indexed fragment (e.g., relatively slow diffusion or transport) relative to a relatively fast rate at which amplification occurs to fill the site with seeded copies of the indexed fragment. In another exemplary embodiment, kinetic exclusion can occur due to a delay (e.g., delayed or slow activation) in the formation of the first copy of the indexed fragment seeded at the site relative to a relatively fast rate at which subsequent copies are made to fill the site. In this example, an individual site may be seeded with several different indexed fragments (e.g., several indexed fragments may be present at each site prior to amplification). However, the first copy formation for any given indexed fragment may be activated randomly such that the average rate of first copy formation is relatively slow compared to the rate at which subsequent copies are generated. In this case, although an individual site may be seeded with several different indexed fragments, kinetic exclusion allows only one of the indexed fragments to be amplified. More specifically, when a first indexed fragment is activated for amplification, the site is rapidly filled to capacity with its copy, thereby preventing copies of a second indexed fragment from being made at that site.

[0138] In one embodiment, the method is carried out to (i) transport the indexed fragments to the amplification site at an average transport rate, and (ii) simultaneously amplify the indexed fragments at the amplification site at an average amplification rate, where the average amplification rate exceeds the average transport rate (US Pat. No. 9,169,513). Thus, kinetic exclusion can be achieved in such an embodiment using a relatively slow rate of transport. For example, a sufficiently low concentration of the indexed fragments can be selected to achieve a desired average transport rate, with lower concentrations resulting in a slower average rate of transport. Alternatively, or in addition, a high viscosity solution and / or the presence of a molecular crowding agent in the solution can be used to reduce the transport rate. Examples of useful molecular crowding agents include, but are not limited to, polyethylene glycol (PEG), ficoll, dextran, or polyvinyl alcohol. Exemplary molecular crowding agents and formulations are shown in US Pat. No. 7,399,590, which is incorporated herein by reference. Another factor that can be adjusted to achieve a desired transport rate is the average size of the target nucleic acid.

[0139] Amplification reagents can include additional components that facilitate amplicon formation and in some cases increase the rate of amplicon formation. An example is recombinase. Recombinase can facilitate amplicon formation by allowing repeated invasion / extension. More specifically, recombinase can facilitate the invasion of the indexed fragment by polymerase and the extension of the primer by polymerase, using the indexed fragment as a template for amplicon formation. This process can be repeated as a chain reaction, with the amplicons generated from each round of invasion / extension serving as templates in subsequent rounds. This process can occur more quickly than standard PCR, since no denaturation cycles (e.g., through heating or chemical denaturation) are required. Thus, recombinase-facilitated amplification can be carried out isothermally. To facilitate amplification, it is generally desirable for recombinase-facilitated amplification reagents to include ATP or other nucleotides (or analogs thereof that are hydrolysis-resistant in some cases). A mixture of recombinase and single-stranded binding (SSB) proteins is particularly beneficial, since SSB can further facilitate amplification. Exemplary formulations for recombinase-facilitated amplification include those marketed as TwistAmp kits by TwistDx (Cambridge, UK). Useful components and reaction conditions for recombinase-facilitated amplification reagents are set forth in U.S. Patent No. 5,223,414 and U.S. Patent No. 7,399,590.

[0140] Another example of a component that can be included in an amplification reagent that promotes amplicon formation and in some cases increases the rate of amplicon formation is a helicase. Helicase can promote amplicon formation by allowing a chain reaction of amplicon formation. This process can occur more quickly than standard PCR because no denaturation cycles (e.g., through heating or chemical denaturation) are required. Thus, helicase-promoted amplification can be carried out isothermally. A mixture of helicase and single-stranded binding (SSB) proteins is particularly useful, since SSB can further promote amplification. Exemplary formulations for helicase-promoted amplification include those commercially available as IsoAmp kits from Biohelix (Beverly, MA). Further examples of useful formulations that include helicase proteins are described in U.S. Patent No. 7,399,590 and U.S. Patent No. 7,829,284.

[0141] Yet another example of a component that can be included in an amplification reagent that facilitates amplicon formation, and in some cases increases the rate of amplicon formation, is an origin binding protein.

[0142] Use in sequencing / sequencing methods After the indexed fragment is attached to the surface, the sequence of the indexed fragment that is immobilized and amplified is determined.Sequencing can be carried out using any suitable sequencing technology, and the method of determining the sequence of the indexed fragment that is immobilized and amplified, including strand resynthesis, is known in the art, and is described, for example, in Bignell et al. (US Pat. No. 8,053,192), Gunderson et al. (WO2016 / 130704), Shen et al. (US Pat. No. 8,895,249) and Pipenburg et al. (US Pat. No. 9,309,502).

[0143] The methods described herein can be used with a variety of nucleic acid sequencing techniques. Particularly applicable techniques are those in which the nucleic acids are attached to fixed positions on an array so that their relative positions do not change, and the array is imaged repeatedly. Particularly applicable are embodiments in which images are obtained in different color channels that correspond to different labels used, for example, to distinguish one nucleotide base type from another. In some embodiments, the process of determining the nucleotide sequence of the indexed fragments can be an automated process. A preferred embodiment includes sequencing-by-synthesis ("SBS") techniques.

[0144] SBS technology generally involves the enzymatic extension of nascent nucleic acid chain through repeated addition of nucleotides to template chain.In the traditional method of SBS, a single nucleotide monomer can be provided to target nucleic acid in each delivery in the presence of polymerase.However, in the method described herein, multiple types of nucleotide monomers can be provided to target nucleic acid in each delivery in the presence of polymerase.

[0145] In one embodiment, the nucleotide monomer comprises a locked nucleic acid (LNA) or a bridged nucleic acid (BNA). The use of an LNA or BNA in the nucleotide monomer increases the hybridization strength between the nucleotide monomer and the sequencing primer sequence present on the immobilized indexed fragment.

[0146] SBS can use nucleotide monomers with terminator moieties or those that lack any terminator moieties. As shown in more detail herein, methods that use nucleotide monomers that lack terminators include, for example, pyrosequencing and sequencing using γ-phosphate-labeled nucleotides. In methods that use nucleotide monomers that lack terminators, the number of nucleotides added in each cycle is generally variable and depends on the template sequence and nucleotide delivery mode. For SBS techniques that use nucleotide monomers with terminator moieties, the terminator can be effectively irreversible under the sequencing conditions used, as in the case of traditional Sanger sequencing that uses dideoxynucleotides, or the terminator can be reversible, as in the case of the sequencing method developed by Solexa (now Illumina, Inc.).

[0147] SBS techniques can use nucleotide monomers with or without a label moiety. Thus, incorporation events can be detected based on the properties of the label, such as the fluorescence of the label; the properties of the nucleotide monomer, such as molecular weight or charge; by-products of nucleotide incorporation, such as the release of pyrophosphate; and the like. In embodiments in which two or more different nucleotides are present in the sequencing reagent, the different nucleotides can be distinguished from each other, or the two or more different labels can be indistinguishable under the detection technique used. For example, the different nucleotides present in the sequencing reagent can have different labels, which can be distinguished using appropriate optics, as exemplified by the sequencing method developed by Solexa (now Illumina, Inc.).

[0148] A preferred embodiment involves pyrosequencing technology. Pyrosequencing detects the release of inorganic pyrophosphate (PPi) as specific nucleotides are incorporated into the nascent strand (Ronaghi, M., Karamohamed, S., Pettersson, B., Uhlen, M. and Nyren, P. (1996) "Real-time DNA sequencing using detection of pyrophosphate release." Analytical Biochemistry 242(1), 84-9; Ronaghi, M. (2001) "Pyrosequencing sheds light on DNA sequencing." Genome Res. 11(1), 3-11; Ronaghi, M., Uhlen, M. and Nyren, P. (1998) "A sequencing method based on real-time pyrophosphate." Science 281(5375), 363; U.S. Patent Nos. 6,210,891; 6,258,568 and 6,274,320). In pyrosequencing, the released PPi can be detected by its immediate conversion to adenosine triphosphate (ATP) by ATP sulfurase, and the level of ATP produced is detected through photons generated by luciferase. The nucleic acid to be sequenced can be attached to features in the array, and the array can be imaged to capture the chemiluminescent signal generated by the incorporation of nucleotides in the features of the array. Images can be obtained after treating the array with a particular nucleotide type (e.g., A, T, C, or G). The images obtained after the addition of each nucleotide type differ in which features in the array are detected. These differences in the images reflect the different sequence content of the features in the array. However, the relative position of each feature remains unchanged in the image. The images can be stored, processed, and analyzed using the methods presented herein.For example, the images obtained after processing of the array with each different nucleotide type can be treated in the same manner as illustrated herein for images obtained from different detection channels for reversible terminator-based sequencing methods.

[0149] In another exemplary type of SBS, cycle sequencing is achieved by stepwise addition of reversible terminator nucleotides, which contain cleavable or photobleachable dye labels, such as those described in WO 04 / 018497 and U.S. Patent No. 7,057,026.This approach has been commercialized by Solexa (now Illumina Inc.) and is also described in WO 91 / 06678 and WO 07 / 123744.The availability of fluorescently labeled terminators that can both reverse the termination and cleave the fluorescent label facilitates efficient cyclic reversible termination (CRT) sequencing.Polymerase can also be simultaneously engineered to efficiently incorporate and extend these modified nucleotides.

[0150] In some reversible terminator-based sequencing embodiments, the label does not substantially inhibit extension under SBS reaction conditions. However, the detection label may be removable, for example, by cleavage or degradation. Following incorporation of the label into the arrayed nucleic acid features, an image can be captured. In certain embodiments, each cycle includes simultaneous delivery of four different nucleotide types to the array, each nucleotide type having a spectrally distinguishable label. Four images can then be obtained, each using a detection channel that is selective for one of the four different labels. Alternatively, different nucleotide types can be added sequentially, and images of the array can be obtained between each addition step. In such embodiments, each image shows nucleic acid features that have incorporated a particular type of nucleotide. Different features are present or absent in different images due to the different sequence content of each feature. However, the relative positions of the features remain unchanged within the images. As shown herein, images obtained from such reversible terminator-SBS methods can be stored, processed, and analyzed. After the image capture step, the label can be removed, and the reversible terminator moiety can be removed for subsequent cycles of nucleotide addition and detection. Removal of the label after they are detected in a particular cycle and before subsequent cycles can provide the advantage of reducing background signals and crosstalk during cycles. Examples of useful labeling and removal methods are provided herein.

[0151] In certain embodiments, some or all of the nucleotide monomers can include a reversible terminator. In such embodiments, the reversible terminator / cleavable fluorophore can include a fluorophore linked to a ribose moiety through a 3' ester linkage (Metzker, Genome Res. 15:1767-1776 (2005)). Other approaches have separated the terminator chemistry from the cleavage of the fluorescent label (Ruparel et al., Proc Natl Acad Sci USA 102: 5932-7 (2005)). Ruparel et al. described the development of a reversible terminator that uses a small 3' allyl group to block extension, but can be easily unblocked by a short treatment with palladium catalyst. The fluorophore was attached to the base through a photocleavable linker that can be easily cleaved by 30 seconds of exposure to long wavelength UV light. Thus, disulfide reduction or photocleavage can be used as a cleavable linker. Another approach to reversible termination is the use of natural termination that occurs after the placement of a bulky dye on the dNTP. The presence of a charged bulky dye on the dNTP can act as an effective terminator through steric and / or electrostatic hindrance. The presence of one incorporation event blocks further incorporation unless the dye is removed. Cleavage of the dye removes the fluorophore, effectively reversing the termination. Examples of modified nucleotides are also described in U.S. Pat. Nos. 7,427,673 and 7,057,026.

[0152] Further exemplary SBS systems and methods that can be utilized in the methods and systems described herein are described in U.S. Patent Application Publication Nos. 2007 / 0166705, 2006 / 0188901, 2006 / 0240439, 2006 / 0281109, 2012 / 0270305 and 2013 / 0260372, U.S. Patent No. 7,057,026, WO 05 / 065814, U.S. Patent Application Publication No. 2005 / 0100900, WO 06 / 064199 and WO 07 / 010251.

[0153] Some embodiments may use detection of four different nucleotides using less than four different labels. For example, SBS may be performed using the methods and systems described in the incorporated materials of US Patent Application Publication No. 2013 / 0079232. As a first example, a set of nucleotide types may be detected at the same wavelength, but may be differentiated based on a difference in the intensity of one member of the set compared to the others, or based on a change to one member of the set (e.g., through chemical, photochemical, or physical modification) that causes a clear signal to appear or disappear compared to the signal detected for the other members of the set. As a second example, three of the four different nucleotide types may be detected under certain conditions, but the fourth nucleotide type may lack a label that is detectable under those conditions, or may be minimally detected under those conditions (e.g., minimal detection due to background fluorescence, etc.). The incorporation of the first three nucleotide types into the nucleic acid may be determined based on the presence of their respective signals, and the incorporation of the fourth nucleotide type into the nucleic acid may be determined based on the absence or minimal detection of any signal. As a third example, one nucleotide type may include label(s) that are detected in two different channels, while the other nucleotide type is detected in only one of the channels. The three exemplary configurations above are not considered to be mutually exclusive and may be used in various combinations.An exemplary embodiment combining all three examples is a fluorescence-based SBS method that uses a first nucleotide type detected in a first channel (e.g., dATP having a label that is detected in the first channel when excited by a first excitation wavelength), a second nucleotide type detected in a second channel (e.g., dCTP having a label that is detected in the second channel when excited by a second excitation wavelength), a third nucleotide type detected in the first and second channels (e.g., dTTP having at least one label that is detected in both channels when excited by the first and / or second excitation wavelengths), and a fourth nucleotide type that is not detected or minimally detected in any channel and lacks a label (e.g., dGTP having no label).

[0154] Furthermore, as described in the incorporated materials of U.S. Patent Application Publication No. 2013 / 0079232, sequencing data can be obtained using a single channel. In such a so-called one-dye sequencing approach, a first nucleotide type is labeled, but the label is removed after the first image is generated, a second nucleotide type is labeled only after the first image is generated, a third nucleotide type retains its label in both the first and second images, and a fourth nucleotide type remains unlabeled in both images.

[0155] Some embodiments may use sequencing by ligation techniques. Such techniques use DNA ligase to incorporate oligonucleotides and identify the incorporation of such oligonucleotides. The oligonucleotides generally have different labels that correlate with the identity of the specific nucleotide in the sequence to which the oligonucleotide hybridizes. As with other SBS methods, images can be obtained after treatment of an array of nucleic acid features with labeled sequencing reagents. Each image shows nucleic acid features that incorporate a specific type of label. Due to the different sequence content of each feature, different features are present or absent in different images, but the relative positions of the features remain unchanged in the images. As shown herein, images obtained from ligation-based sequencing methods can be stored, processed, and analyzed. Exemplary SBS systems and methods that can be utilized in the methods and systems described herein are described in U.S. Patent Nos. 6,969,488, 6,172,218, and 6,306,597.

[0156] Some embodiments may use nanopore sequencing (Deamer, DW & Akeson, M. "Nanopores and nucleic acids: prospects for ultrarapid sequencing." Trends Biotechnol. 18, 147-151 (2000); Deamer, D. and D. Branton, "Characterization of nucleic acids by nanopore analysis", Acc. Chem. Res. 35:817-825 (2002); Li, J., M. Gershow, D. Stein, E. Brandin, and JA Golovchenko, "DNA molecules and configurations in a solid-state nanopore microscope" Nat. Mater. 2:611-615 (2003)). In such embodiments, the indexed fragment passes through a nanopore. The nanopore may be a synthetic pore or a biological membrane protein, such as α-hemolysin. As the indexed fragment passes through the nanopore, each base pair can be identified by measuring the fluctuation in the electrical conductance of the pore.(U.S. Pat. No. 7,001,792; Soni, GV & Meller, "A. Progress toward ultrafast DNA sequencing using solid-state nanopores." Clin. Chem. 53, 1996-2001 (2007); Healy, K. "Nanopore-based single-molecule DNA analysis." Nanomed. 2, 459-481 (2007); Cockroft, SL, Chu, J., Amorin, M. & Ghadiri, MR "A single-molecule nanopore device detects DNA polymerase activity with single-nucleotide resolution." J. Am. Chem. Soc. 130, 818-820 (2008)). As provided herein, data obtained from nanopore sequencing can be stored, processed, and analyzed. In particular, the data can be processed as images by the exemplary processing of optical images and other images provided herein.

[0157] Some embodiments can use methods that include real-time monitoring of DNA polymerase activity.Nucleotide incorporation can be detected through fluorescence resonance energy transfer (FRET) interactions between fluorophore-containing polymerases and γ-phosphate-labeled nucleotides, for example, as described in U.S. Patent Nos. 7,329,492 and 7,211,414, or nucleotide incorporation can be detected with zero-mode waveguides, for example, as described in U.S. Patent No. 7,315,019, and using fluorescent nucleotide analogs and engineered polymerases, for example, as described in U.S. Patent No. 7,405,281 and U.S. Patent Application Publication No. 2008 / 0108082. Illumination intensity can be restricted to a zeptoliter-scale volume around the surface-linked polymerase so that incorporation of fluorescently labeled nucleotides can be observed with low background (Levene, MJ et al. "Zero-mode waveguides for single-molecule analysis at high concentrations." Science 299, 682-686 (2003); Lundquist, PM et al. "Parallel confocal detection of single molecules in real time." Opt. Lett. 33, 1026-1028 (2008); Korlach, J. et al. "Selective aluminum passivation for targeted immobilization of single DNA polymerase molecules in zero-mode waveguide nano structures." Proc. Natl. Acad. Sci. USA 105, 1176-1181 (2008)). Images obtained from such methods can be stored, processed and analyzed as described herein.

[0158] Some SBS embodiments include detection of protons released after incorporation of nucleotides into extension products.For example, sequencing based on detection of released protons can use electric detectors and related technologies available from Ion Torrent (Guilford, CT, a subsidiary of Life Technologies), or the sequencing methods and systems described in US Patent Application Publication Nos. 2009 / 0026082, 2009 / 0127589, 2010 / 0137143 and 2010 / 0282617.The method described herein for amplifying target nucleic acid using kinetic exclusion can be easily applied to the substrate used for detection of protons.More specifically, the method described herein can be used to generate a clonal population of amplicons used to detect protons.

[0159] The SBS method described above can be conveniently carried out in a multiplex format to simultaneously manipulate multiple different indexed fragments. In certain embodiments, the different indexed fragments can be processed in a common reaction vessel or on the surface of a specific substrate. This allows for convenient delivery of sequencing reagents, removal of unreacted reagents, and multiplexed detection of incorporation events. In embodiments using surface-bound target nucleic acids, the indexed fragments can be in the form of an array format. In an array format, the indexed fragments can be generally bound to the surface in a spatially distinguishable manner. The indexed fragments can be bound by direct covalent attachment, attachment to beads or other particles, or binding to a polymerase or other molecule attached to the surface. The array can contain a single copy of the indexed fragment at each site (also called a feature), or multiple copies with the same sequence can be present at each site or feature. Multiple copies can be generated by amplification methods such as bridge amplification or emulsion PCR, as described in more detail herein.

[0160] The methods presented herein can be used to fabricate, for example, at least about 10 features / cm 2 , 100 features / cm 2 , 500 features / cm 2 , 1,000 features / cm 2 , 5,000 features / cm 2 , 10,000 features / cm 2 , 50,000 features / cm 2 , 100,000 features / cm 2 , 1,000,000 features / cm 2 , 5,000,000 features / cm 2 Arrays having any of a variety of densities of features, including 100 sq. m or higher, can be used.

[0161] The advantage of the methods presented herein is that they can be used in multiple 2The object of the present invention is to provide rapid and efficient detection of nucleic acids in parallel. Thus, the present disclosure provides an integrated system capable of preparing and detecting nucleic acids using techniques known in the art, such as those exemplified herein. Thus, the integrated system of the present disclosure can include fluidic components capable of delivering amplification and / or sequencing reagents to one or more immobilized indexed fragments, and the system includes components such as pumps, valves, reservoirs, fluidic systems, etc. A flow cell can be configured and / or used in the integrated system for detection of target nucleic acids. Exemplary flow cells are described, for example, in U.S. Patent Application Publication No. 2010 / 0111768 and U.S. Patent Application No. 13 / 273,666. As exemplified for the flow cell, one or more of the fluidic components of the integrated system can be used for amplification and detection methods. Taking a nucleic acid sequencing embodiment as an example, one or more of the fluidic components of the integrated system can be used for the amplification method shown herein and for the delivery of sequencing reagents in a sequencing method such as the one exemplified above. Alternatively, the integrated system can include separate fluidic systems for performing the amplification method and for performing the detection method. Examples of integrated sequencing systems that can also generate amplified nucleic acids and sequence nucleic acids include, but are not limited to, the MiSeq™ platform (Illumina, Inc., San Diego, CA) and the device described in US Patent Application No. 13 / 273,666.

[0162] Compositions are also provided herein.During the implementation of the methods described herein, various compositions can be produced.For example, compositions comprising triple-indexed nucleic acid fragments can be produced.Multi-well plates are also provided, wherein the wells of the multi-well plates comprise triple-indexed nucleic acid fragments.

[0163] Kits are also provided herein. In one embodiment, the kit is for preparing a sequencing library. The kit includes the transposomes and / or linear amplification mediators described herein in a suitable packaging material in an amount sufficient for at least one assay or use. Other components may be included, such as one or more nucleic acids including primers, indexes, ubiquitous sequences, or combinations thereof, as appropriate. Other components that may be included are reagents and solutions, such as buffers. Instructions for use of the packaged components are also typically included. As used herein, the phrase "packaging material" refers to one or more physical structures used to contain the contents of the kit. The packaging material is constructed by routine methods to generally provide a sterile, contaminant-free environment. The packaging material can have a label indicating that the components can be used to generate a sequencing library. Additionally, the packaging material contains instructions indicating how the materials in the kit are used. As used herein, the term "package" refers to a container, such as glass, plastic, paper, foil, etc., that is capable of holding the components of the kit within fixed limits. "Instructions for use" typically include specific language describing at least one assay method parameter, such as reagent concentrations or relative amounts of reagents and sample to be mixed, maintenance time for the reagent / sample mixture, temperature, buffer conditions, etc.

[0164] Exemplary embodiments Embodiment 1. A method for preparing a sequencing library comprising nucleic acids from a plurality of single nuclei or cells, comprising: providing a plurality of isolated nuclei or cells in a first plurality of compartments, each compartment comprising a subset of the isolated nuclei or cells, the nuclei or cells comprising nucleic acid fragments; introducing a linear amplification agent into a cell or nucleus; amplifying the nucleic acid fragment by linear amplification; processing each subset of nuclei or cells to generate indexed nuclei or cells, the processing comprising adding a first compartment-specific index sequence to a nucleic acid fragment present in the isolated nuclei or cells to provide an indexed nucleic acid present in the isolated nuclei or cells, the processing comprising ligation, primer extension, hybridization, amplification or transposition; and The method includes combining the indexed nuclei or cells to generate pooled indexed nuclei or cells, thereby generating a sequencing library from the plurality of nuclei or cells.

[0165] Embodiment 2. The method of embodiment 1, wherein amplification occurs prior to treatment.

[0166] Embodiment 3. The method of embodiment 1, wherein the treatment occurs before amplification.

[0167] Embodiment 4. A method for preparing a sequencing library comprising nucleic acids from a plurality of single nuclei or cells, comprising: providing a plurality of isolated nuclei or cells containing the nucleic acid fragment; introducing a linear amplification agent into the isolated nucleus or cell; distributing the isolated nuclei or cells into a first plurality of compartments, each compartment containing a subset of the isolated nuclei or cells; amplifying the nucleic acid fragment by linear amplification; processing each subset of isolated nuclei or cells to generate indexed nuclei or cells, the processing comprising adding a first compartment-specific index sequence to a nucleic acid fragment present in the isolated nuclei or cells to provide an indexed nucleic acid present in the isolated nuclei or cells, the processing comprising ligation, primer extension, amplification or transposition; The method includes combining the indexed nuclei to generate pooled indexed nuclei or cells, thereby generating a sequencing library from the plurality of nuclei or cells.

[0168] Embodiment 5. A method for preparing a sequencing library comprising nucleic acids from a plurality of single nuclei or cells, comprising: providing a plurality of isolated nuclei or cells in a first plurality of compartments, each compartment comprising a subset of the isolated nuclei or cells, the nuclei or cells comprising the nucleic acid fragment; processing each subset of nuclei or cells to generate indexed nuclei or cells, the processing comprising adding to the nucleic acid fragments present in the isolated nuclei or cells (i) a first compartment-specific index sequence to provide indexed nucleic acids present in the isolated nuclei or cells, and (ii) adding a nucleotide sequence recognized by a linear amplification agent, the processing comprising ligation, primer extension, hybridization, amplification or transposition; introducing a linear amplification agent into a cell or nucleus; amplifying the nucleic acid fragment by linear amplification; and The method includes combining the indexed nuclei or cells to generate pooled indexed nuclei or cells, thereby generating a sequencing library from the plurality of nuclei or cells.

[0169] Embodiment 6. The method of any one of embodiments 1-5, wherein the linear amplification agent comprises a phage RNA polymerase or a linear amplification primer.

[0170] Embodiment 7. The method of any one of embodiments 1-6, wherein the nucleic acid fragment comprises a T7 promoter and the phage RNA polymerase comprises T7 RNA polymerase.

[0171] Embodiment 8. The method of any one of embodiments 1-7, wherein introducing a linear amplification agent comprises adding a linear amplification agent to the nucleic acid fragment present in the isolated nucleus or cell.

[0172] Embodiment 9. The method of any one of embodiments 1-8, further comprising exposing a plurality of the isolated nuclei or cells of each compartment to a predetermined condition.

[0173] Embodiment 10. The method of any one of embodiments 1-9, further comprising isolating nuclei from the plurality of cells after exposing.

[0174] Embodiment 11. The method of any one of embodiments 1-10, further comprising exposing the plurality of isolated nuclei or cells to a predetermined condition.

[0175] Embodiment 12. The method of any one of embodiments 1 to 11, further comprising subjecting the isolated nuclei to conditions that generate nucleosome-depleted nuclei while maintaining the integrity of the isolated nuclei.

[0176] Embodiment 13. The process comprises: contacting each subset with a transposome complex, wherein the transposome complex in each compartment comprises a first index sequence that is different from the first index sequences in the other compartments; and 13. The method of any one of embodiments 1-12, comprising fragmenting the nucleic acids in the subset into a plurality of nucleic acids and incorporating a first index sequence into at least one strand of the nucleic acids to generate an indexed nucleus or cell comprising the indexed nucleic acid.

[0177] Embodiment 14. The method of any one of embodiments 1-13, wherein the treating comprises contacting each subset with a reverse transcriptase and a primer that anneals to RNA molecules in the isolated nuclei, wherein the primer in each compartment comprises a first index sequence that differs from the first index sequence in the other compartment, to generate indexed nuclei or cells that comprise the indexed nucleic acid.

[0178] Embodiment 15. The method of any one of embodiments 1 to 14, wherein the contacting further comprises a target-specific primer that anneals to a specific nucleotide sequence.

[0179]

[0023] Embodiment 16. The method of any one of embodiments 1 to 15, wherein the process of adding a first compartment-specific index sequence comprises a two-step process of adding a nucleotide sequence comprising a ubiquitous sequence to the nucleic acid fragments and then adding a first compartment-specific index sequence to the nucleic acid fragments.

[0180] Embodiment 17. The method of any one of embodiments 1 to 16, wherein the adding comprises a transposome complex comprising a ubiquitous sequence.

[0181] Embodiment 18. The method of any one of embodiments 1 to 17, wherein the processing comprises adding a first index to a DNA nucleic acid present in the isolated nucleus or cell, adding a first index to an RNA nucleic acid present in the isolated nucleus or cell, or a combination thereof.

[0182] Embodiment 19. The method of any one of embodiments 1-18, wherein adding a first index sequence to the RNA nucleic acids comprises contacting each subset with a reverse transcriptase and a primer that anneals to an RNA molecule in the isolated nuclei or cells, wherein the primer in each compartment comprises a first compartment-specific index sequence, to generate indexed nuclei or cells comprising the indexed nucleic acids.

[0183]

[0023] Embodiment 20. Adding a first index sequence to a DNA nucleic acid comprises: contacting each subset with a transposome complex, wherein the transposome complex in each compartment comprises a first compartment-specific index sequence; and 20. The method of any one of embodiments 1-19, comprising fragmenting the nucleic acids in the subset into a plurality of nucleic acids and incorporating a first compartment-specific index sequence into at least one strand of the nucleic acids to generate an indexed nucleus or cell comprising the indexed nucleic acid.

[0184] Embodiment 21. The method of any one of embodiments 1 to 20, wherein the first index sequence added to the DNA nucleic acid and the first index sequence added to the RNA nucleic acid in each compartment are identical.

[0185]

[0023] Embodiment 22. The method of any one of embodiments 1 to 21, wherein the first index sequence added to the DNA nucleic acid and the first index sequence added to the RNA nucleic acid in each compartment are not identical.

[0186] Embodiment 23. The method of any one of embodiments 1 to 22, further comprising exponential amplification of the nucleic acid fragment, the exponential amplification comprising a target-specific primer that anneals to a specific nucleotide sequence.

[0187] Embodiment 24. After combining, distributing the pooled subset of indexed nuclei or cells into a second plurality of compartments; and 24. The method of any one of embodiments 1-23, further comprising introducing a second compartment-specific index sequence into the indexed nucleic acid to generate a dual-indexed nucleus or cell comprising the dual-indexed nucleic acid, wherein the introducing comprises ligation, primer extension, amplification or transposition.

[0188]

[0023] Embodiment 25. Combining the dual-indexed nuclei to generate pooled dual-indexed nuclei or cells. distributing the pooled subset of doubly indexed nuclei or cells into a third plurality of compartments; and 25. The method of any one of embodiments 1-24, further comprising introducing a third compartment-specific index sequence into the indexed nucleic acid to generate a triply-indexed nucleus or cell comprising the triply-indexed nucleic acid, wherein the introducing comprises ligation, primer extension, amplification or transposition.

[0189]

[0031] Embodiment 26. The method of any one of embodiments 1 to 25, further comprising processing the indexed nuclei or cells for methylation analysis to generate nucleic acid fragments suitable for methylation analysis.

[0190]

[0031] Embodiment 27. The method of any one of embodiments 1 to 26, further comprising subjecting the indexed nuclei or cells to proximity ligation to generate nucleic acid fragments suitable for analysis of chromatin conformation.

[0191]

[0023] Embodiment 28. The method of any one of embodiments 1 to 27, further comprising amplifying the nucleic acid fragments of the sequencing library to generate DNA nanoballs.

[0192] Embodiment 29. The method of any one of embodiments 1 to 28, wherein the compartment comprises a well or a droplet.

[0193]

[0036] Embodiment 30. The method of any one of embodiments 1-29, wherein each compartment of the first plurality of compartments comprises between 50 and 100,000,000 nuclei or cells.

[0194]

[0036] Embodiment 31. The method of any one of embodiments 1-29, wherein each compartment of the second plurality of compartments comprises between 50 and 100,000,000 nuclei or cells.

[0195] Embodiment 32. A method comprising the steps of: providing a surface comprising a plurality of amplification sites, the amplification sites comprising at least two populations of attached single-stranded capture oligonucleotides having free 3' ends; and 32. The method of any one of embodiments 1-31, further comprising contacting a surface comprising the amplification sites with the indexed fragments under conditions suitable to generate a plurality of amplification sites, each comprising a clonal population of amplicons from an individual fragment comprising a plurality of indexes.

[0196] Embodiment 33. A method for preparing a sequencing library comprising nucleic acids from a plurality of single cells, comprising: (a) providing isolated nuclei from a plurality of cells; (b) subjecting the isolated nuclei to a chemical treatment that produces nucleosome-depleted nuclei while maintaining the integrity of the isolated nuclei; (c) distributing a subset of the nucleosome-depleted nuclei into a first plurality of compartments and contacting each subset with a transposome complex, wherein the transposome complex in each compartment comprises a transposase and a first index sequence that is different from the first index sequences in the other compartments; (d) fragmenting the nucleic acid in the subset of nucleosome-depleted nuclei into a plurality of nucleic acid fragments and incorporating a first index sequence into at least one strand of the nucleic acid fragment to generate indexed nuclei comprising indexed nucleic acid fragments, wherein the indexed nucleic acid fragments remain attached to the transposase; (d) combining the indexed nuclei to generate pooled indexed nuclei; (e) distributing a subset of the pooled indexed nuclei into a second plurality of compartments and contacting each subset with a hairpin ligation duplex under conditions suitable for ligation of the hairpin ligation duplex to one or both ends of the indexed nucleic acid fragments to provide doubly-indexed nucleic acid fragments, wherein the hairpin ligation duplex comprises a second index sequence that differs from the second index sequences in the other compartments; (f) combining the dual-indexed nuclei to generate pooled indexed nuclei; (g) distributing a subset of the pooled doubly indexed nuclei into a third plurality of compartments; (h) lysing the doubly indexed nuclei; (i) processing the dual-indexed nuclear fragments to contain a third index sequence that is different from the third index sequence in the other compartment; and (j) combining the triplet index fragments, thereby generating a sequencing library comprising the entire genome nucleic acid from a plurality of single cells. The method includes: EXAMPLES

[0197] The present disclosure is illustrated by the following examples, it being understood that the specific examples, materials, amounts and procedures are to be interpreted broadly in accordance with the scope and spirit of the disclosure set forth herein. EXAMPLES

[0198] High-throughput single-cell sequencing by linear amplification Conventional methods for single-cell genome sequencing are limited in terms of uniformity and throughput. Here, we describe "sci-L3", a high-throughput, high-coverage single-cell sequencing method that combines single-cell combinatorial indexing ("sci") and linear ("L") amplification. The sci-L3 method employs a unidirectional three-level ("3") indexing scheme that minimizes amplification bias while enabling an exponential increase in throughput. We demonstrate the generalizability of the sci-L3 framework through proof-of-concept studies of single-cell whole genome sequencing ("sci-L3-WGS"), targeted genome sequencing ("sci-L3-targeted-seq") and genome and transcriptome co-assays ("sci-L3-RNA / DNA"). We apply sci-L3-WGS to profile the genomes of >10,000 sperm and sperm precursors from F1 hybrid male mice, map 86,786 crossovers, and characterize rare chromosome missegregation events, including cases of genome-wide equal chromosome segregation in male meiosis. We expect that the sci-L3 assay can be applied to link CRISPR perturbations and measurements of genome stability, to fully characterize the recombination landscape, and for other goals requiring high-throughput, high-coverage single-cell genome sequencing.

[0199] Introduction Modern single-cell genome sequencing techniques have two major limitations. First, most methods require segmenting individual cells, which limits throughput. Second, most amplification methods are PCR-based and therefore suffer from exponential amplification bias. To solve the first problem, we and colleagues developed single-cell combinatorial indexing ("sci-"), where we perform several rounds of split-pool molecular barcoding to specifically tag the nucleic acid content of single cells, thereby allowing an exponential increase in throughput with each successive round of indexing. Sci-methods have been successfully developed to profile chromatin accessibility (sci-ATAC-seq), transcriptome (sci-RNA-seq), genome (sci-DNA-seq), methylome (sci-MET), and chromosomal conformation (sci-Hi-C) in large numbers of single cells (Cao et al., 2017;Cusanovich et al., 2015;Mulqueen et al., 2018;Ramani et al., 2017;Vitak et al., 2017). To solve the second problem, linear amplification by T7-based transcription offers a potential solution that has been previously utilized in the context of single-cell assays (Eberwine et al., 1992;Hashimshony et al., 2012;Sos et al., 2016). For example, Chen et al. recently developed Linear Amplification by Transposon Insertion (LIANTI), which uses a Tn5 transposon to fragment the genome and simultaneously insert a T7 RNA promoter for in vitro transcription (IVT). The RNA copies generated from the DNA template cannot serve as templates for further amplification. Thus, all copies are derived directly from the original DNA template. By avoiding exponential amplification, LIANTI maintains uniformity and minimizes sequence errors. However, this method is low-throughput, as it requires sequential library preparation from each single cell (Chen et al., 2017).

[0200] We developed sci-L3, which incorporates single-cell combinatorial indexing and linear amplification to enable an exponential increase in throughput while minimizing amplification bias. Through triple molecular barcoding, sci-L3 increases the throughput of LIANTI to at least thousands and potentially millions of cells per experiment while retaining the advantages of linear amplification. We demonstrate the generalizability of the sci-L3 framework through proof-of-concept for single-cell whole genome sequencing ("sci-L3-WGS"), targeted genome sequencing ("sci-L3-targeted-seq") and genome and transcriptome co-assays ("sci-L3-RNA / DNA"). As further demonstration, we apply sci-L3-WGS to map an unprecedented number of meiotic crossovers and rare chromosomal missegregation events in immature and mature male germ cells from infertile interspecific (B6 x Spretus) F1 male mice and fertile intraspecific (B6 x Cast) F1 male mice.

[0201] design A potential technological route to minimize amplification bias while increasing throughput would be to simply combine the "sci" and "LIANTI" methods. However, the molecular structure of LIANTI, in which the T7 promoter is inserted through a Tn5 transposon, offers the opportunity for cell barcoding only twice, which would limit the throughput to a few thousand single cells per experiment. It is further limited to profiling genomic DNA (Chen et al., 2017; Sos et al., 2016). In developing sci-L3, we integrated three-times cell barcoding ("3-level") by introducing the T7 promoter through single-cell combinatorial indexing, linear amplification and ligation (Figure 3A). Compared to simply combining "sci" and "LIANTI", the sci-L3 approach has several major advantages. First, the potential throughput is exponentially increased to over 1,000,000 cells per experiment at a much reduced cost by three-level indexing (Cao et al., 2019). Second, the unidirectional nature of single-cell barcoding allows sci-L3 to be readily converted to targeted sequencing ("targeted-seq") in addition to whole genome sequencing ("WGS"), enabling linking CRISPR perturbations and resulting genomic instability, as well as other applications where it is desirable to sequence specific genomic loci across a large number of single cells. Third, as a generalizable linear amplification and high-throughput cell barcoding scheme, sci-L3 offers the flexibility to be adapted to other single-cell assays and co-assays with minor modifications of the protocol, as demonstrated here by our proof-of-concept of a sci-L3-based single-cell RNA / DNA co-assay.

[0202] result Proof of concept for sci-L3-WGS and sci-L3-targeted-seq The three-level combinatorial indexing and amplification scheme for sci-L3-WGS and sci-L3-targeted-seq is shown in Figure 3A: (i) cells are fixed with formaldehyde and nucleosomes are depleted by SDS (Vitak et al., 2017). The resulting nuclei are then evenly distributed into 24 wells. (ii) A first round of barcodes is added into each of the 24 wells by indexed Tn5 insertion ("tagmentation"). Unlike LIANTI, where the Tn5 transposon contains a barcode-free T7 promoter, a spacer sequence is included 5' of the barcode, which serves as a "landing pad" for the subsequent ligation step (for details of the Tn5 transposon design, see Figure 4 and Example 2, "Methods and molecular design of sci-L3-WGS and sci-L3-targeted-seq" section). (iii) All of the nuclei are pooled and redistributed uniformly into 64 new wells; a second barcode is added by ligation, which contains a T7 promoter sequence located outside both barcodes. (iv) All of the nuclei are pooled together again, sorted by fluorescence activated cell sorting (FACS) cytometry, and distributed into final round wells at up to 300 cells per well. Note that nuclei of different ploidy can be gated and enriched by DAPI (4',6-diamidino-2-phenylindole) staining. Additionally, simple dilution is an alternative to FACS that can reduce loss rates. (v) Sorted nuclei are lysed and subjected to in situ gap extension to form double-stranded T7 promoters. This is followed by IVT, reverse transcription (RT) and second strand synthesis (SSS) to linearly amplify the genome. A third barcode is added during the SSS step together with a unique molecular identifier (UMI) to tag individual IVT transcripts.(vi) Double-stranded DNA molecules, each containing three barcodes that define their cell of origin (Figure 3B, top), are compatible with traditional library construction methods if the goal is single-cell WGS (e.g., adding sequence adapters by ligation (Figure 3B, middle) or tagmentation) or with slightly modified methods if the goal is single-cell targeted DNA-seq (e.g., adding a PCR step in which one of the primers is target specific (Figure 3B, bottom)).

[0203] As an initial proof of concept, we performed sci-L3-WGS on a mix of mouse and human cells. For >95% of the resulting single-cell genomes, the majority of reads mapped to either the mouse or human genome. Occasional "collisions" arise from the fortuitous use of the same combination of barcodes by two or more cells (Figure 3C). The performance of sci-L3-WGS is compared to LIANTI as well as our previous PCR-based sci-DNA-seq method in Table 1. We highlight several advantages of sci-L3-WGS: 1) 90% of sorted cells are typically recovered, compared to 60% recovery by PCR-based sci-DNA-seq (Vitak et al., 2017); 2) with 40% fewer raw reads (329M by sci-L3-WGS vs. 549M by sci-DNA-seq), sci-L3-WGS yielded sequence coverage of approximately 97,000 unique Tn5 insertions per cell, compared to approximately 30,000 unique insertions by sci-DNA-seq (>3-fold improvement). While sequencing fewer cells to a higher depth and maintaining a higher library complexity than sci-DNA-seq, approximately 660,000 unique Tn5 insertions per cell were observed, suggesting a further improvement of >20-fold; 3) the rate of mappable reads improves from 61% with LIANTI to 86% with sci-L3-WGS, likely due to the fact that LIANTI is fully in-tube and therefore difficult to remove artifactual sequences (e.g., secondary to Tn5 self-insertion), whereas in sci-L3-WGS, nuclei are pelleted several times to remove excess free DNA; 4) Unlike PCR-based amplification, where double-repeat reads provide no information for SNP calling, the "double-repeat" reads of sci-L3-WGS mostly come from independent IVT transcripts polymerized from the original template and are therefore informative for the discovery of novel SNVs or for genotyping of known SNPs.

[0204] Table 1. Performance comparison of sci-DNA-seq vs. sci-L3-WGS vs. LIANTI. sci-DNA-seq data from xSDS method in (Vitak et al., 2017). LIANTI from in-tube method in (Chen et al., 2017). For sci-L3-WGS, the results of libraries yi140 and yi141 (high depth sequencing) and yi144 and yi145 (low depth sequencing) are shown. These four libraries use an optimized protocol with enriched Tn5 transposomes (0.2 μM) and improved RT reaction with additional RNA primers (for details, see Figure 5 and Example 2, "Method and molecular design of sci-L3-WGS and sci-L3-targeted-seq" section). The same color indicates the comparison of interest. Green: percentage of single cells recovered from sorting is improved by 1.9-fold with sci-L3-WGS compared to sci-DNA-seq; Pink: mapping rate of raw reads is improved by 1.4-fold with sci-L3-WGS compared to LIANTI; Yellow: unique insertion sites at various sequencing depths; Rows 1 and 2 are compared with similar number of raw reads, with 3.3-fold improvement with 40% less raw reads with sci-L3-WGS compared to sci-DNA-seq, and rows 1 and 3 are compared with similar library complexity, with 22.4-fold improvement with 20% better Tn5 insertion complexity with sci-L3-WGS compared to LIANTI; Blue: intermediate library complexity showing that methods including both LIANTI and sci-L3-WGS have minimal PCR2 repeats; Orange: number of cells with >50k unique reads recovered is improved by 1.8-fold with sci-L3-WGS compared to LIANTI;

[0205] [Table 1]

[0206] In sci-L3-WGS, Tn5 insertions are on average every 0.5-1.5 kb of the human genome, giving IVTs about 1,000 transcripts. This corresponds to 2,000,000-6,000,000 unique Tn5 insertions and therefore 2,000,000-6,000,000,000 unique genome-derived IVT transcripts per single cell. Sequencing the resulting libraries to saturation with respect to the number of unique IVT transcripts is clearly not currently practical. Here, we define the "sequencing depth" for each library as the ratio between the number of unique transcripts sequenced and the number of unique Tn5 insertion sites mapped. In this study, most libraries were sequenced to a depth of 1.1× to 2×, resulting in a coverage of 0.5%-5% of each cell's genome. The distribution of unique Tn5 insertion sites per cell in the human / mouse proof-of-concept experiment is shown in Figure 3D, and for other experiments in Figure 5. The estimated relative chromosome copy numbers for a representative single cell are shown in Figure 3E, and their distribution across all cells is shown in Figure 3F. To extrapolate the expected genome coverage per single cell at higher sequencing depths, the number of unique insertion sites is fitted as a function of sequencing depth (Figure 5G). We expect to observe 4.2M and 6.0M unique insertions per cell at 5x and 10x sequencing depths, respectively, corresponding to 16% and 22% coverage of the genome of individual cells.

[0207] As described above, the double-stranded amplicons generated by sci-L3 (Figure 3B, top) are compatible not only with single-cell WGS (sci-L3-WGS; Figure 3B, middle) but also with single-cell targeted DNA sequencing ("sci-L3-targeted-seq"). Specifically, for targeted sequencing, sequencing adapters can be added by PCR after second strand synthesis, with one primer carrying the third cell barcode while the other primer targets a specific region of the genome (Figure 3B, bottom). To quantify the efficiency of recovery by sci-L3-targeted-seq, lentiviral CRISPR libraries were incorporated at low MOI (for details, see Example 2, "Methods and molecular design of sci-L3-WGS and sci-L3-targeted-seq" section) and DNA sequences corresponding to sgRNA spacers were recovered by sci-L3-targeted-seq. For 97 out of 1003 single cells, a single integrated sgRNA can be successfully recovered. This 10% efficiency per haplotype is broadly consistent with the 22% genome coverage estimated above by estimating the sequencing depth (Figure 5G).

[0208] At the molecular level, we note that we modified the "sci" and "LIANTI" methods in several ways. In summary: 1) we changed the design of the Tn5 transposon to accommodate ligation, thus allowing more than two rounds of "sci", an approach that potentially generalizes to other single-cell assays; 2) we added a loop structure in the T7 promoter to promote intramolecular ligation; and 3) we modified the RT scheme to require only successful ligation at one of the two ends of the first round of barcoded molecules. Assuming that a single ligation event has an efficiency of 50%, this modification gives a 75% success rate at the ligation step instead of 25% (a comparison is shown in Figure 5). The structure of the molecules after each barcoding step is shown in Figure 4, and the rationale for these designs is discussed in Example 2, "Methods and molecular design of sci-L3-WGS and sci-L3-targeted-seq" section. Scalability and cost are also discussed in Example 2 and Table 2. For libraries of 100, 1000, 10,000, and 1,000,000 single cells, we estimate that the cost of sci-L3-WGS is 14%, 1.5%, 0.26%, and 0.014% of processing the same number of cells with LIANTI. To increase throughput (e.g., the cost of building a library for 1,000,000 cells with 3-level sci-L3-WGS at 5% collision rate is about $8,000) or to reduce collision rate (e.g., the cost of building a library for 10,000 cells with 1% collision rate with 3-level sci-L3-WGS is about $1,500), we can leverage the use of three levels of combinatorial indexing instead of two.

[0209] Table 2. Cost calculations for sci-L3-WGS. The current method involves three-level indexing, which not only increases throughput and reduces barcode collisions, but also significantly reduces the cost per cell of library preparation. This is for two reasons: 1) with two-level indexing, one needs to start with many more Tn5 transposome complexes to profile a similar number of cells, which substantially increases the cost; 2) with two-level indexing, one is limited to sorting a significantly smaller number of nuclei per well before IVT, RT and column purification, which also substantially increases the cost. For processing of ~10k and ~1M cells, we estimate that three-level sci-L3-WGS is nearly 400-fold and 7,000-fold less expensive than LIANTI per cell.

[0210] [Table 2-1] [Table 2-2] [Table 2-3] [Table 2-4]

[0211] Leveraging sci-L3-WGS for single-cell RNA / DNA co-assays We recognized that the sci-L3-WGS scheme could potentially be adapted to other aspects of molecular biology with minor modifications to the protocol. To demonstrate this, we performed a proof-of-concept experiment of the sci-L3-RNA / DNA co-assay. Briefly, the first round of DNA barcoding is performed by Tn5 insertion as in sci-L3-WGS, but the first round of RNA barcoding is performed simultaneously to tag mRNAs by barcode- and UMI-containing poly-T primers through reverse transcription (Figure 6A). Both the Tn5 insertion and the RT primer have overhangs that can mediate the second round of barcoding as well as ligation of the T7 promoter, effectively enabling three-level indexing and subsequent IVT-based linear amplification in a manner almost identical to sci-L3-WGS (Figures 6A-6B, Example 2, "Method and molecular design of sci-L3-RNA / DNA co-assay" section). As a proof-of-concept, mouse cells were mixed together with cells from two human cell lines and a sci-L3-RNA / DNA co-assay was performed. For the vast majority of cells, read data mapped to the mouse or human genome for both RNA (5.2% collision rate) and DNA (6.6% collision rate) (Figure 6C-6D). Moreover, consistent with good co-assay, 100% of cells were assigned the same species label by their RNA and DNA profiles. As a further check, we visualized human cells in the t-SNE space based on their RNA profiles. As expected, they separated into two clusters. Labeling of individual cells based on the presence or absence of a Y chromosome consistently identified clusters as corresponding to BJ cells (male) or HEK293T cells (female) with 96.5% accuracy (Figure 6E).

[0212] Single-cell DNA profiling of mouse germ cells by sci-L3-WGS During normal mitotic cell division, diploid chromosomes undergo replication to generate four copies of DNA, and sister chromatids segregate into reciprocal daughter cells. Daughter cells each receive one copy of DNA sequences inherited from the mother and one from the father, and almost always maintain heterozygosity at sequences proximal to the centromere (Figure 7A). Rarely, chromosomes undergo mitotic crossing over between chromosome homologs, which can sometimes result in diploid cells with loss of heterozygosity (LOH) at sequences distal to the centromere to the crossing over, if the two recombined chromatids segregate into different daughter cells (Figure 7B-C).

[0213] In meiosis, sister chromatids first co-segregate into identical daughter cells, while homologs segregate into reciprocal daughter cells during meiosis I ("MI"), also known as "reductional segregation," resulting in 2C cells (DNA content of non-replicating diploid cells) with loss of heterozygosity (LOH) at centromere-proximal sequences (Figure 7D-E). For successful reductional segregation of chromosomes in MI (Figure 7D), crossovers initiated by double-strand breaks (DSBs) catalyzed by Spo11 (Baudat et al., 2000; Keeney et al., 1997; Romanienko and Camerini-Otero, 2000) provide linkage and necessary tension (Hong et al., 2013) between chromosome homologs. Rarely, chromosomes segregate meiotically without any interhomolog crossovers, resulting in parthenogenetic disomy (UPD). After MI, these 2C cells then undergo mitosis-like chromosome segregation, also called "equal segregation," during meiosis II ("MII"), so that sister chromatids separate to form 1C gametes (Figure 7E). Hereafter, because our study focuses primarily on MI, we refer to meiosis / reductive segregation during MI in which sister chromatids segregate together as "reductive segregation," and mitosis-like / equal segregation during MI in which sister chromatids segregate apart as "equal segregation."

[0214] So far, most studies on the relationship between crossover locations and chromosome segregation have been performed by imaging (Wang et al., 2017a, 2017b), which cannot fully characterize the underlying genomic sequences that are prone to meiotic crossovers. Although some assays allow fine mapping of meiotic DSB hotspots (Lange et al., 2016;Smagulova et al., 2011, 2016), these assays do not directly map meiotic crossovers. Assays that dissect crossovers vs. non-crossovers at a fine scale are limited to a few hotspots (Cole et al., 2014). As a result, we know considerably less about the relationship between crossovers and chromosome-scale features, such as replication domains, than we know about meiotic DSB hotspots (Baudat et al., 2013;Choi and Henderson, 2015;Yamada et al., 2017). Genome-wide meiotic crossover maps have been generated by mapping tetrads in yeast ( Mancera et al., 2008 ; Zhang et al., 2017 ), single human sperm and complete human female meiosis ( Hou et al., 2013 ; Lu et al., 2012 ; Ottolini et al., 2015 ; Wang et al., 2012 ). With the exception of studies of human female meiosis that analyzed 87 complete meiotic divisions in total (Hou et al., 2013; Ottolini et al., 2015), most crossover maps are limited in at least three respects: 1) mature 1C gametes, whose cells have completed both rounds of meiosis, are analyzed, which precludes direct observation of the more informative intermediate 2C cells to assess whether and how often chromosomes undergo reduction vs. equal segregation during MI (Figure 7); 2) aberrant cells are counterselected because they are unable to progress to the mature gamete state; 3) analysis by sequencing of single sperm or oocytes is throughput-limited and limited to at most a few hundred cells, thus potentially missing rare events.Even for fertile crosses, the number of offspring that can be reasonably generated and genotyped is severely limited ( Liu et al., 2014 ).

[0215] To address all of these limitations at once, we applied sci-L3-WGS to the infertile progeny of an interspecific cross (female Mus musculus domesticus C57BL / 6 ("B6") × male Mus spretus SPRET / Ei (later "Spret")) as well as the fertile progeny of an intraspecific hybrid (female B6 × male Mus musculus castaneous CAST / Ei ("Cast")). By sequencing sperm with a highly scalable technology, we can map an unprecedented number of crossover events in fertile and infertile hybrids for a mammalian system. Furthermore, by leveraging the throughput of sci-L3-WGS to recover profiles from rare 2C secondary spermatocytes, we can also simultaneously investigate crossover and chromosome missegregation from the same single cell.

[0216] Unlike inbred males as well as (B6 × Cast) F1 males, whose epididymides store millions of mature sperm (Berletch et al., 2015), the epididymis of (B6 × Spret) F1 males contains very few morphologically mature sperm and a limited number of round germ cells of unknown ploidy (Figure 8A-B). Interestingly, during FACS, we observed a significantly higher percentage of 2C cells than would be expected in a "normal" epididymis dominated by 1C sperm (Figure 8C-D). The number of cells recovered and their estimated ploidy are listed in Table 3. In contrast, and as expected, the epididymis of (B6 × Cast) F1 males contained almost entirely 1C sperm (Figure 8E). For this cross, we therefore sorted 1C and 2C cells from dissociated testes (Figure 8F). ×

[0217] Table 3. Number of cells recovered and cell ploidy, (B6 x Spret) epididymis. Note that we did not generate sequencing libraries for all cells sorted; for example, the 2C library in Exp1 only contains a subset of cells. We also gated broadly for 1C cells (up to 58 DAPI signals in one particular well) and are only able to enrich for 1C cells to approximately 51-55% due to the abundance of 2C cells in this cross.

[0218] [Table 3]

[0219] For cells from F1 males from both (B6×Spret) and (B6×Cast) crosses, we proceeded with linear amplification, synthesis of the second strand to add the third barcode, library preparation and sequencing (details are in Example 2, section "Setup of sci-L3-WGS experiments in (B6×Spret) and (B6×Cast) crosses"). An important point is that although 1C and 2C cells can be distinguished informatically, their relative abundance still influences our analysis. Specifically, in (B6×Spret) crosses, 1C cells are rare, and therefore any "doublets" (e.g., two 1C cells that are anchored to each other or that incidentally receive the same barcode) do not contribute substantially to the 2C population. In contrast, in (B6×Cast) crosses, despite enrichment, the majority of cells are 1C (approximately 85%, Figure 8G), and therefore there may be many 1C doublets that mimic 2C cells. We consider how to informatively distinguish 1C doublets from genuine 2C cells in a later section.

[0220] M2 cells exhibit clustered reduced or even chromosome segregation Chromosome segregation in M2 cells from sterile (B6 x Spret) crosses. We first sought to analyze meiosis in cells from the epididymis of infertile (B6 x Spretus) F1 males, obtained as described above. Over the two sci-L3-WGS experiments, we profiled the genomes of 2,689 (92% of 2,919 sorted cells with >10k raw reads) and 4,239 (94% of 4,497 sorted cells with >30k raw reads) single cells. The number of unique mapping reads is shown in Figure 5F. At a sequencing depth of 1.6x and 1.4x for the two libraries (details in Figure 5), we obtained a median of ∼70k and ∼144k unique Tn5 sites per cell, corresponding to an intermediate genome coverage of 0.7% and 1.4%, respectively.

[0221] To identify crossover breakpoints, we performed a hidden Markov model (HMM) relying on high-quality read data that could be unambiguously assigned to B6 vs. Spret (see Example 2, "Methods for Bioinformatics and Statistical Analysis" section). We characterized crossovers in 1,663 1C cells, a representative example of which is shown in Figure 9A. Furthermore, we searched for approximately 5,200 2C cells for crossover events. While most of these 5,200 may simply be somatic cells, we were surprised to identify 292 2C cells with a significant number of crossovers, which we refer to as "M2 cells" (Figures 9B and 9C). Even more surprising, a significant proportion of these cells showed even segregation rather than reduction.

[0222] If the chromosomes separate in a reduced fashion after a crossover occurs between two chromosome homologs, the region between the centromere and the crossover location becomes homozygous, but heterozygosity is maintained downstream of the crossover (Figure 7D). However, if the chromosomes separate in an even fashion, LOH is observed centromeric distal to the crossover if the recombined chromatids separate (Figure 7B). We show one example of an M2 cell whose chromosome undergoes the expected reductional segregation in Figure 9B (note the consistent homozygosity between the centromere and the crossover), and one example of an M2 cell whose chromosome unexpectedly undergoes even segregation in Figure 9C (note the consistent heterozygosity between the centromere and the crossover). In total, across 292 M2 cells, we observed 4,162 examples of chromosomes that underwent reductional segregation, of which 3,740 had crossovers (90%), and 1,310 examples of chromosomes that underwent even segregation, of which 636 had crossovers (49%). However, it should be noted that the number of crossover events in evenly segregated chromosomes may be higher since we are unable to identify a subset of crossover outcomes (Figure 7C), whereas all crossovers in reduced segregated chromosomes can be detected.

[0223] Although we observe many examples of cells in which some chromosomes show reduced segregation and others show even segregation, the segregation patterns of individual chromosomes within M2 cells do not appear to be independent. If the chromosomes in each cell independently choose reduced vs. even segregation, we would expect a binomial distribution of reduced and evenly segregated chromosomes centered on the maximum likelihood estimate (MLE) of the probability of reduced segregation, p (p=0.76 from data, 4162 / 5472), with roughly 3 / 4 of chromosomes segregating in a reduced manner and 1 / 4 segregating evenly (Figure 9D). However, among the 292 M2 cells profiled, we observe 202 cells with at least 15 chromosomes segregated in a reduced manner and 38 cells with at least 15 chromosomes segregated in an even manner (Figure 9E; this contrasts with the 148 and 0 cells, respectively, expected when assuming independence; p=4e-23, Fisher's exact test). The bias of individual M2 cells to predominantly undergo reduction or equal segregation suggests a possible cell-autonomous global sensing mechanism for cells to decide whether to continue through meiosis or revert to mitosis.

[0224] The cells can be further classified according to whether the chromosomes of M2 cells have crossovers (Figure 9F). Reduced segregation chromosomes (pink in Figure 9F) appear to have more crossovers than evenly segregated chromosomes (green in Figure 9F). However, unlike in reduced segregation chromosomes, where all crossovers can be detected as centromeric LOH, evenly segregated chromosomes only have LOH when two recombined chromatids segregate into reciprocal daughter cells (Figure 7B). If the recombined chromatids instead cosegregate, heterozygosity is maintained across the chromosome despite an undetectable linkage switch (Figure 7C). In Figure 9F, the ratio of evenly segregated chromosomes with (shown in green) versus without (shown in blue) observable LOH is roughly 1:1. This could mean that evenly segregated chromosomes have a 50% chance of segregating the recombinant chromatids together, if those fully heterozygous chromosomes (shown in blue) do in fact harbor the linkage switch, or alternatively, that evenly segregated chromosomes will always segregate the recombinant chromatids apart, with the crossover frequency being halved compared to reduced segregation chromosomes.

[0225] Fragmentary or whole-chromosome LOH is known to be rare in mammalian mitotic cells. Nevertheless, to exclude a mitotic origin of such events, we investigated such events in the Patski cell line, a spontaneously immortalized cell line derived from a female (B6 x Spret) F1 mouse. We analyzed 1,107 single cells from Patski by sci-L3-WGS, among which we found an average of 0.36 UPD chromosomes and 0.098 fragmentary LOH events per cell, a rate that is considerably lower compared to M2 cells. We also note that these events are not necessarily independent. For example, UPDs that appear early in the passage of a cell line can be shared by a large proportion of progeny cells, so the rate of independent LOH events is probably even lower. The distribution of these events (relatively uniform in Spretus-derived chromosomes and not uniform in B6-derived chromosomes) is plotted in Figure 10F.

[0226] Taken together, the contrast between the low rate of mitotic LOH (expected) and the relatively high rate of 2C cells showing even segregation (unexpected), both measured by the same technique, confirms that the latter is highly unlikely to correspond to somatic cells. In the next section, by analyzing fertile (B6 × Cast) crosses, we further show that 1) the genome-wide even segregation events observed here are not an artifact of a doublet of two 1C cells, and 2) such segregation events also occur, albeit at a low rate, in fertile intraspecific hybrids.

[0227] Chromosome segregation in M2 cells from fertile (B6 x Cast) crosses. We wondered whether equal segregation also occurs during MI in fertile progeny of intraspecific (B6 × Cast) F1 males. As shown above, the epididymis from this cross consists almost entirely of 1C mature spermatozoa. Therefore, we enriched for 2C secondary spermatocytes from the whole testis. We then performed sci-L3-WGS on cells from both the epididymis and testis.

[0228] In the first sci-L3-WGS experiment with this cross, which was mainly performed for quality control to investigate the recovery and barcode collision rates, 1C round spermatozoa were distributed uniformly and only sorted for 1C cells after two rounds of barcoding. Doublets, identified by the fact that they are non-1C, allow the rate of barcode collisions to be quantified. Among 2,400 sorted cells (200 / well), 2,127 (89%) were recovered with >7,000 reads per cell. 2,008 of these were 1C with meiotic crossovers, indicating a barcode collision rate of 5.5%. At a sequencing depth of 1.06×, we obtained a median of about 60k unique Tn5 insertions per cell, which corresponded to an intermediate genome coverage of about 0.6%.

[0229] In the second sci-L3-WGS experiment on this cross, 1C round spermatozoa from the testis ("barcode group 1"), 2C cells from the testis ("barcode group 2"; contaminated with a large number of 1C spermatozoa as shown in Figure 8F), and 1C mature spermatozoa from the epididymis ("barcode group 3", Example 2, section "Setup of sci-L3-WGS experiments in (B6 x Spret) and (B6 x Cast) crosses") were tagged in separate wells during the first round of barcoding. As a further enrichment, for a subset of wells during the FACS step of sci-L3-WGS, we specifically gated for 2C cells (15.5% of all cells, Figure 8G). At a sequencing depth of 1.09x, we obtained a median of approximately 94k unique Tn5 insertions per cell, which corresponded to an intermediate genome coverage of approximately 0.9%.

[0230] In total, we recovered 3,539 1C and 1,477 non-1C cells from this second sci-L3-WGS experiment. Interestingly, >97% of the 1C cells were derived from barcode group 1 (n=1,853) and group 2 (n=1,598) rather than group 3 (n=88), indicating that mature sperm from the epididymis are not well recovered by sci-L3-WGS. This suggests that the 1C cells recovered from the above (B6 x Spret) cross are also probably not from mature sperm but from round sperm, consistent with the low number of sperm with mature morphology in Figure 8B.

[0231] The 1,477 non-1C cells were derived from both barcode 1 (n=1,104; likely doublets of 1C round spermatozoa) and barcode 2 (n=373; likely a mixture of bona fide M2 ​​and 1C doublets). To identify the signature of 1C doublets, we profiled non-1C cells from barcode 1 (which had been specifically presorted for 1C content and were therefore unlikely to contain bona fide M2 ​​cells). SNPs proximal to the centromeres of 1C cells that completed both rounds of meiosis should be of B6 or Cast origin. In 1C doublets, these regions have equal chances of appearing as heterozygous or homozygous. Thus, within any given 1C doublet, the number of chromosomes that appear to be evenly segregated, as well as the number that appear to be reduced segregated, should follow a binomial distribution with n=19 and p=0.5. Indeed, this is what is observed for 1C doublets from barcode group 1 (p = 0.53 for the distribution of the proportion of evenly segregating chromosomes that deviated from the binomial (19, 0.5), chi-square test; Figure 11A-B). In fact, there are only 11 1C doublet cells with at least 15 chromosomes that appear to segregate in a consistent manner, whether evenly or reduced.

[0232] In contrast, non-1C cells from barcode group 2 show a very different partitioning. Among such 373 cells, 258 resemble 1C doublets from barcode group 1 in that they have numbers of chromosomes similar to even or reduced segregation patterns. The remaining 115 cells are biased, with at least 15 chromosomes segregating in a consistent manner whether even or reduced (Figure 11C-E; 115 / 373 in barcode group 2 vs. 11 / 1,104 in barcode group 1; p=3e-70, chi-squared test), and many show a completely even (n=6) or completely reduced (n=91) pattern.

[0233] Finite mixture models for fitting three populations of non-1C cells To take this into account more formally, data from each experiment were fitted to a Bayesian finite mixture of three binomial distributions. Details are provided in Example 2, "Finite mixture model for fitting three populations of non-1C cells" section and in Figure 12, and key conclusions are summarized here. First, non-1C cells from intraspecific (B6 x Cast) F1 male testes (i.e., barcode group 2) are estimated to contain a cell subset with reduced (28%) vs. even (2%) segregation, as well as likely 1C doublets (69%) (Figure 12B). This proportion is different for M2 cells from interspecific (B6 x Spret) F1 males, which are estimated to contain a cell subset with reduced (66%) vs. even (14%) segregation, as well as likely 1C doublets (20%) (Figure 12C). These analyses support the conclusion that sterile (B6 x Spret) crosses have a significantly higher proportion of cells biased towards even segregation rather than reduced segregation.

[0234] Distribution of meiotic crossovers at the chromosome level Next, we aimed to investigate the genomic correlation of crossover events. Overall, we analyzed 1,663 1C cells with 19,601 crossover breakpoints and 240 M2 cells with 4,184 crossover breakpoints from a (B6 x Spret) cross, as well as 5,547 1C cells with 60,755 crossover breakpoints and 115 M2 cells with 2,246 crossover breakpoints from a (B6 x Cast) cross. To our knowledge, this is an unprecedented dataset in terms of the number of crossover events identified in relation to mammalian meiosis.

[0235] The high-throughput nature of sci-L3-WGS allowed us to analyze a large number of immature germ cells to identify rare cell populations that have completed MI but not MII, as well as to observe meiotic crossover and chromosomal missegregation events in the same cells. In comparing infertile interspecific (B6×Spret) hybrids with fertile intraspecific (B6×Cast) hybrids at the chromosomal level, the following disadvantages are observed for MI: 1) the proportion of M2 cells with at least one crossover on all 19 autosomes is reduced from about 2 / 3 in (B6×Cast) to about 1 / 2 in (B6×Spret); 2) the average number of crossovers per M2 cell is lower in (B6×Spret), but the average number of crossovers per 1C cell is higher; 3) crossover interference is weaker in (B6×Spret), with the median distance between adjacent crossovers being 97 Mb or less. 82 Mb; 4) in (B6×Spret) M2 cells, crossovers tend to occur in the central half of each chromosome arm, in contrast to 1C of both crosses as well as (B6×Cast) M2 cells, which prefer the quartile most distal to the centromere; 5) among M2 cells with biased even or reduced chromosome segregation, (B6×Spret) shows a significantly higher proportion (38 / 240) of whole-genome even segregation than (B6×Cast) (8 / 115); 6) among the whole-genome reduced segregation of M2 cells at MI, the average number of sporadic even segregations (also called reverse segregation (Ottolini et al., 2015)) increases from 0.2 to 1.1. These findings suggest mechanisms that may contribute to or reflect underlying factors that contribute to (B6 x Spret)F1 male infertility, including defects in the formation and placement of crossovers, impaired mechanisms to ensure at least one crossover per chromosome, and increased sporadic and genome-wide equal segregation. Details of these analyses are presented in Figures 10, 13, and 14 and in Example 2, "Distribution of Meiotic Crossovers at the Chromosome Level" section.

[0236] Distribution of meiotic crossover events in relation to genomic context Genomic features regulating cross-hot nis To assess the distribution of crossovers at a finer scale, we collapsed all crossover events to generate "hot-nise maps" along each mouse chromosome. We first compared these maps to single-stranded DNA sequencing (SSDS) maps (Brick et al., 2018; Smagulova et al., 2011, 2016) and Spo11 oligonucleotide complex maps (Lange et al., 2016), which identify meiotic DSB hotspots at the highest resolution (Figure S15A). DSB maps in the B6 strain from these two mapping methods strongly correlate with each other along 100 kb windows (rho = 0.87, p < 2e-308). Although our 1C and M2 cell cross pileups correlate with each other (rho = 0.67 for (B6 × Spret) cross and rho = 0.55 for (B6 × Cast) cross, p < 2e-308 for both, Figure 15B-C), both deviate from the DSB map. Relatedly, the PRDM9 gene, a key player for hotspot specification, has evolved to bind distinct motifs between divergent mouse lineages and even between mouse subspecies (Davies et al., 2016; Gregorova et al., 2018). We consider its potential impact on the differences between the two crosses in Example 2, section "Impact of PRDM9 on cross hotspots."

[0237] Only 10% of meiosis-specific DSBs are repaired as crossovers. We next investigated what factors besides Spo11 breaks contribute to crossover formation by constructing linear models with Bayesian model averaging (BMA) (Clyde et al., 2011). As applied here, BMA takes a weighted average of over 15,000 variable selection models examined and weights them by the posterior probability of each model, which, unlike some other variable selection techniques such as Lasso regression, accounts for the uncertainty in model selection. We quantified marginal inclusion probabilities (MIPs) for about 80 potential explanatory variables. Features known to be associated with meiotic crossovers, such as Spo11 break sites, GC content, etc., are included with high probability in almost all models (Figure S16A, Figure S17). For example, regions with high GC content are hotter for crossover formation. We also found a few more features not previously linked to meiotic crossovers, such as repeat sequences and specific families of chromatin marks, especially early replication domains. Correlation matrices between crossover hotnis and all features are plotted in Figures 18-19 for each crossover. A summary of the features used as well as simple linear models and BMA are included. The breakpoint resolution (median values ​​~150 kb for (B6 x Spret) and ~250 kb for (B6 x Cast); Figure 16B) is comparable to previous efforts to map meiotic crossovers by single-cell sequencing (150-500 kb) (Lu et al., 2012; Ottolini et al., 2015; Wang et al., 2012). However, the greater library complexity afforded by sci-L3-WGS allowed us to achieve this at a significantly lower sequencing depth.

[0238] Many of the features that correlate with crossover formation are consistent between the (B6 x Spret) and (B6 x Cast) crosses, but some are not. For example, the positional bias of crossover formation appears to be different. In 1C cells from both crosses, as well as in M2 cells from the (B6 x Cast) cross, crossovers are underrepresented within 10 Mb of the centromere and tend to occur near the telomere rather than in the rightmost positional "quartile" (Figure 18). However, in M2 cells from the (B6 x Spret) cross, crossovers are underrepresented near the centromere as well as near the telomere and tend to occur in the middle quartile (Figure 19). This trend holds true in a linear model in which we account for contributions from all other features.

[0239] The location of crossovers can greatly affect the amount of tension exerted between chromosome homologs, which in turn promotes proper chromosome segregation. Therefore, we explored this in more detail by taking only the right-most crossover for each chromosome in each cell and examining its location along the chromosome arm in each crossover (de Boer et al., 2015). Taking into account interchromosomal variability by linear mixed effects model, we estimate that the location of right-most crossovers in (B6×Spret) crossovers is on average more proximal to the 1.6 Mb centromere than those in (B6×Cast) crossovers in 1C cells (Figure 20A, p=1e-13, F-test), but more proximal to the 5.5 Mb centromere in M2 cells (Figure 16C, p=2.2e-15). We note that the rightmost crossovers in M2 cells are more proximal to the centromere than those in 1C cells in both crosses, but tend to be more so in the (B6×Spret) cross (Figure S16D) than in the (B6×Cast) cross (Figure S16B). These differences suggest that a subset of M2 cells in the (B6×Spret) cross, whose crossovers occur closer to the centromere than necessary, may be unable to mature into 1C cells, possibly due to a defect in MII segregation. Similarly, although the number of events is limited, we also compared the location of crossovers in M2 cells in which chromosome segregation is biased, and found that in both crosses, crossovers in cells with biased even segregation are more distal to the centromere than those in cells with biased reduced segregation, a difference of 13.7 Mb (p=4e-15) in the (B6×Cast) cross and 8.7 Mb (p=6e-14) in the (B6×Spret) cross (Figure S16C-D). This suggests a possible MI segregation defect in cells with crossovers too close to the telomere. We propose a tentative model to explain this observation in Figure 20E.

[0240] Cellular heterogeneity regarding crossover breakpoints Although 1C and M2 cells appear broadly similar in cross pileups (Figure 15), we wondered whether there is any structure in the features that influence cross distribution in the subset of single cells. To explore this, we collected cross-related information for each single cell for each of the 78 features (Example 2, Bioinformatics and Statistical Analysis Methods section). We then used principal component analysis (PCA) on a matrix with each row being one single cell and each column being one summarized feature value. For the (B6 x Spret) cross, the first two principal components (PCs) capture 26% of the variance, and for the (B6 x Cast) cross, PC1 and PC3 capture 17% of the variance. In both crosses, 1C and M2 cells are separated into two clusters by these PCs. In Figures 21 and 22, we plot each feature projected onto these PCs. Chromosome segregation of the crossovers, the location of the crossovers on the parthenogenetic chromosomes and chromosome tetrads are features that appear to facilitate the separation of 1C and M2 cells.

[0241] Predicting crossover zones from genomic features Finally, we aimed to exploit the large number of events observed here to build a predictive model of crossover location. Specifically, we built a binary response linear model, where 1 is crossover zone taken from genomes from the same zone length distribution and 0 is random zone (details are in Example 2, "Bioinformatics and Statistical Analysis Methods" section). Using the same 76 features as in the BMA analysis, we are able to predict crossover zones with an average receiver operator curve (ROC) area under the curve (AUC) of 0.73 for the held-out data for the (B6 x Spret) crossover. A similar average AUC of 0.72 is achieved with a subset of 25 variables with high inclusion probability (MIP>0.5) identified by BMA (Figure 16E). Similarly, for the (B6 x Cast) crossover, we achieve an average AUC of 0.85 when using all features or a subset of 25 features with MIP>0.5 (Figure 16F).

[0242] Consideration Here, we describe a framework, sci-L3, that combines three levels of single-cell combinatorial indexing and linear amplification. We demonstrate that sci-L3 is applicable to single-cell whole genome sequencing (sci-L3-WGS), single-cell targeted DNA sequencing (sci-L3-targeted-seq) and single-cell co-assays of genomes and transcriptomes (sci-L3-RNA / DNA). With sci-L3-WGS, at least tens of thousands and potentially millions of genomes of single cells can be processed in a two-day experiment with a library construction cost of $0.14 per cell for 10k cells and $0.008 per cell for 1M cells. The throughput of sci-L3-WGS is several orders of magnitude higher than alternative single-cell WGS methods based on linear amplification such as “in-tube” LIANTI (Chen et al., 2017). It further improves the number of specific molecules recovered from each single cell from low thousands (Pellegrino et al., 2018) or low tens of thousands (Vitak et al., 2017) to hundreds of thousands.

[0243] We applied sci-L3-WGS to study male mouse meiosis and identified an unexpected population of M2 cells. The single-cell nature of the data also allowed us to simultaneously characterize meiotic crossovers and chromosome missegregation. Reverse segregation events have been previously observed in a full analysis of human female meiosis (Ottolini et al., 2015), and similar events are observed here in the context of mouse male meiosis (i.e., equal segregation of one or several chromosomes). Among the 292 M2 cells we analyzed from a (B6 × Spret) cross, individual cells were biased toward equal or reduced chromosome segregation, suggesting a global sensing mechanism that determines whether the cell should continue meiosis or revert to mitotic segregation of its chromosomes. Moreover, to our knowledge, for the first time in mammalian meiosis, we observed multiple examples of whole-genome equal segregation during MI, suggesting a cell-autonomous rather than chromosome-autonomous manner of equal segregation. Although more rare in the fertile (B6 × Cast) cross, we identified such events in both crosses.

[0244] Chromosome-autonomous mechanisms, especially those that would be expected in interspecific (B6 × Spret) crosses (2 -19The high incidence of genome-wide reverse segregation compared to the rate of crossing over raises many questions that remain to be answered. We present a model in Figure 23 to highlight some of the open questions. In normal MI, centromere cohesion is maintained at reduction segregation, and sister chromatids centromeric proximal to the crossing over do not separate until MII (pattern 1 in Figure 23D). Equal segregation at MI indicates premature centromere cohesin separation (patterns 2 and / or 3 in Figure 23D). Previous studies have also shown that homolog pairing may be defective in these F1 crosses due to erosion of PRDM9 binding sites (Davies et al., 2016; Gregorova et al., 2018; Smagulova et al., 2016), and pairing problems are likely more severe in interspecific crosses. In Example 2, section "Speculations regarding the causes and consequences of reverse segregation," we speculate on: 1) what might cause premature kinetochore cohesin segregation, 2) whether one crossover is sufficient for proper reduction segregation, and 3) what the consequences of equal segregation at MI might be.

[0245] The improved genome coverage enabled high-resolution mapping of crossover breakpoints compared to other single-cell sequencing methods, and the throughput for mapping a total of about 87,000 crossovers allowed us to better characterize the genomic and epigenomic features associated with crossover hot niches by pile-up data. How the continuous band of crossover hot niches is formed by many factors is discussed in Example 2, "Crossover hot niches and associated (epi)genomic factors".

[0246] One important difference in the development of sci-L3 from simply combining a high-throughput single-cell combinatorial indexing ("sci") scheme with linear amplification through transposon insertion (LIANTI) is the introduction of a T7 promoter by ligation, which not only allows cell barcoding more than twice and further increases throughput at a significantly reduced cost, but also provides the flexibility to generalize this method to other single-cell assays by small tweaks of the protocol. As a first example, we demonstrate that sci-L3-WGS can be easily adapted to sci-L3-targeted-seq. Single-cell targeted sequencing has been reported by the 10X Genome Research platform, but to the best of our knowledge, it is of RNA transcripts rather than DNA loci. Although the current 10% "recovery rate" per haplotype may not be ideal for targeted sequencing, it is mitigated by the large number of cells that can be analyzed. As a second example, we demonstrate that sci-L3-WGS can also be adapted to sci-L3-RNA / DNA co-assays. We anticipate that it may further be possible to adapt sci-L3 to ATAC-seq, bisulfite-seq and Hi-C for single-cell profiling of chromatin accessibility, methylome and chromatin conformation, respectively, which may have advantages over published sci-methods for these goals in terms of throughput and amplification uniformity (Cusanovich et al., 2015; Mulqueen et al., 2018; Ramani et al., 2017).

[0247] In summary, sci-L3-WGS, sci-L3-targeted-seq and sci-L3-RNA / DNA co-assays expand the toolset for single-cell sequencing. In this study, we further demonstrate how sci-L3-WGS can provide a systematic and quantitative picture of meiotic recombination, shedding light on rare genome-wide chromosome missegregation events with an unprecedented combination of throughput. We anticipate that sci-L3 methods will be highly informative in other contexts where single-cell genome sequencing proves transformative, for example, to study rare interhomolog mitotic crossovers and to dissect the genetic heterogeneity and development of cancer.

[0248] References [Table A-1] [Table A-2] [Table A-3] [Table A-4] [Table A-5] EXAMPLES

[0249] Finite mixture models for fitting three populations of non-1C cells Non-1C cells recovered from (B6 × Cast) hybrids from barcode group 2 include 1C doublets, cells that appear to be biased toward even segregation, and cells that appear to be biased toward reduced segregation. To quantify their relative proportions, we fit the data to a mixture of three binomial distributions with probabilities of evenly segregating chromosomes of 0.01, 0.48, and 0.95, and mixture proportions of 0.28, 0.69, and 0.02 (Figure 12A). In contrast, when we attempt to similarly fit non-1C cells from barcode group 1 to a mixture of three binomial distributions, we obtain probabilities of evenly segregating chromosomes of 0.46, 0.5, and 0.53 (all close to 0.5), and mixture proportions of 0.24, 0.44, and 0.31 (Figure 12B).

[0250] Toward asking whether the proportion of M2 cells biased toward even versus reduced segregation differs between fertile and infertile crosses, we were able to similarly fit chromosome data from (B6 × Spret) crosses (Figure 9E), which gave probabilities of evenly segregated chromosomes of 0.05, 0.39, and 0.91, and admixture proportions of 0.66, 0.2, and 0.14 (Figure 12C). These proportions suggest that infertile (B6 × Spret) crosses have a higher proportion of cells biased toward even rather than reduced segregation.

[0251] Distribution of meiotic crossovers at the chromosome level Based on 1,663 1C cells with 19,601 crossover breakpoints and 240 M2 cells with 4,184 crossover breakpoints from (B6 x Spret) crosses, and 5,547 1C cells with 60,755 crossover breakpoints and 115 M2 cells with 2,246 crossover breakpoints from (B6 x Cast) crosses, we first considered the distribution of meiotic crossovers across chromosomes. Crossover density is defined here as the average number of crossovers per cell per division per Mb multiplied by 2 (1C cells) or 1 (M2 cells). In (B6 x Spret) crosses, a strong negative correlation was observed between chromosome size and crossover density in 1C cells (Figure 13A, r = -0.66, p = 0.002). Consistent with previous findings (Lange et al., 2016), this correlation is only partially explained by Spo11 oligonucleotide complex density (r=-0.46, p<0.05), suggesting that smaller chromosomes maintain more DSBs and that those DSBs are more likely to generate crossovers. This negative correlation is even stronger in M2 cells (Figure 13B, r=-0.83, p=1e-5). In Figure 10A-B, we consider the cases of multiple crossovers per chromosome per cell as single events, which further strengthens the negative correlation (r=-0.87, p=2e-6 in 1C cells; r=-0.91, p=8e-8 in M2 cells). These observations suggest that smaller chromosomes are hotter for crossovers, especially for having at least one crossover per cell division. In (B6 × Cast) crosses, the same trend is observed (Figure 14A-D). 1C cells had an average of 0.62 and 0.58 crossovers per chromosome per cell for inter- and intra-species crosses, respectively, whereas M2 cells had an average of 0.92 and 1.03 crossovers per chromosome per cell (Figures S13C-D, S10C-D). The crossover rate in inter-species M2 cells is only 9% lower than the number of crossovers measured by Mlh1 foci in 4C spermatocytes of B6 inbred mice (Froenicke et al., 2002), despite a 2% sequence variance.The crossover rate in 1C cells is 45% lower than observed in single human sperm sequencing (Lu et al., 2012; Wang et al., 2012). The latter difference may be due in large part to the acrocentric nature of mouse chromosomes. Although interspecific (B6×Spret) crosses have a higher average number of crossovers detected in 1C compared to (B6×Cast) crosses (p=7e-26, Mann-Whitney test), the average number of crossovers in M2 cells is lower (p=2e-10). We note that the percentage of M2 cells with reduced segregation of all 19 autosomes with crossovers to any chromosome is higher in (B6×Cast) crosses (60 / 91 or 66%) than in (B6×Spret) crosses (41 / 80 or 51%) (p=0.06, Fisher's exact test), which may contribute to the infertility of the latter.

[0252] To investigate crossover interference, we took chromosomes with at least two crossovers, plotted the distance between adjacent crossovers, and compared this distribution to the expectation based on random simulations (Figure S13E, S10E, S14E). The median observed distance between crossovers was 82 Mb for (B6 x Spret) and 97 Mb for (B6 x Cast). Both are significantly larger than the expectation of 39 and 42 Mb (p=1e-267 and p<2e-308, respectively, Mann-Whitney test). This is consistent with repulsion of very close crossovers. Crossover interference is stronger in (B6 x Cast) than in (B6 x Spret) crosses, and the distance between adjacent crossovers is longer (p=5e-91).

[0253] We also analyzed the distribution of parthenogenetic chromosomes in each single cell (i.e., no crossing over is observed) (Figure S13F), and for each chromosome in the (B6 x Spret) cross (Figure S13G) (the same trend is valid in the (B6 x Cast) cross, as shown in Figure S14F-G). When normalized by length, shorter chromosomes show higher crossing over rates, but the rate of parthenogenetic chromosomes (collapsed across all classes of cells) still negatively correlated with chromosome size (Figure S13G; r = -0.91, p = 4.6e-8).

[0254] Although we showed that M2 cells are strongly biased towards either even or reduced segregation of their chromosomes, we also observed hundreds of sporadic even segregation events among cells with at least 15 chromosomes with reduced segregation. This phenomenon has been observed previously and named "reverse segregation" (Ottolini et al., 2015). In Figure 13H, we show the chromosome distribution of these reverse segregation events. Note that the rate of reverse segregation is significantly higher in (B6 x Spret) crosses (mean = 1.1) than in (B6 x Cast) crosses (mean = 0.2, p = 2e-14, Mann-Whitney test), but chromosomes 7 and 11 have the highest reverse segregation rates in both crosses.

[0255] The normalized proportion of reads per cell that map to the mitochondrial genome was next examined (Figure 13I, Figure 10G). 1C cells show a bimodal distribution of mitochondrial DNA "copy numbers", an observation for which we lack a satisfactory explanation. We observed a mild negative correlation between mitochondrial read proportion and the number of crossovers (rho=-0.11, p=3e-6). Interestingly, a limited number of M2 cells, which had even vs. reduced segregation of at least 15 of their chromosomes, had a very different distribution of mitochondrial read proportions (Figure 10G). Consistent with this, mitochondrial read proportions were positively correlated with the number of reduced segregated chromosomes in M2 cells (r=0.18, p=0.005). It should be noted that this cannot be assessed in (B6 x Cast) crosses, as >90% of the sequenced single cells did not have any reads mapped to the mitochondrial genome. The different methods used for nuclear isolation from testis (B6 × Cast) versus epididymis (B6 × Spret), coupled with presorting of nuclei from testis, may have fractionated mitochondria from bulk nuclei.

[0256] Effect of PRDM9 on cross-hot varnish Based on the cross-hot niss map by stacking cross-over breakpoints along chromosomes across the genome (Figure S15), we found that in intraspecific (B6 x Cast) crosses, cross-hot niss correlated better with DSB hot domains mapped in Cast males than in B6 males (rho=0.28 and 0.12, p<2e-308 and p=1e-83, respectively), likely as a result of the Cast PRDM9 allele being semi-dominant in the F1 hybrids. The correlation was stronger with DSB hot domains mapped in (B6 x Cast) F1 animals (rho=0.3, p<2e-308). In the (B6 x Spret) cross, erosion of the PRDM9 consensus binding site results in four types of DSB hotspots defined by the Spo11 oligonucleotide complex map: those conserved between B6 and Spret, termed "symmetric" hotspots, those present only in B6 or Spret, termed "asymmetric" hotspots, and those that do not contain PRDM9 binding sites in either species. All four types of DSB hot domains do not correlate well with crossovers from the (B6 x Spret) cross (rho = 0.13, p = 4e-87 when using all Spo11 hotspots mapped in B6; rho = 0.11, p = 3e-63 when using only "symmetric hotspots"). One possibility is that the DSB sites in the (B6 x Spret) cross are strongly dominated by the SpretPRDM9 allele, and therefore crossover sites are not predicted from DSB hotspots mapped in the B6 phylogenetic background.

[0257] Speculation on the causes and consequences of reverse segregation We observed a high incidence of reverse segregation, especially in interspecific (B6 × Spret) crosses, leading us to speculate about: 1) what might cause premature centromeric cohesin segregation, 2) whether one crossing over is sufficient for proper reduction segregation, and 3) what the consequences of equal segregation in MI might be.

[0258] First, due to insufficient homolog pairing between the B6 and Spret chromosomes, DSBs that should have been successfully repaired from the homolog during meiosis may instead be frequently repaired using sister chromatids as templates, which may cause cohesin disruption ( Storlazzi et al., 2008 ) and lead to premature centromeric cohesin segregation.

[0259] Second, the current model suggests that one interhomolog crossover and proper sister chromatid cohesion are sufficient for chiasma formation despite initial insufficient homolog pairing in interspecific crosses (Figure 23). If the crossover is successfully formed, chromosome segregation should not be impaired. In our study, at the individual chromosome level, many of the evenly segregating chromosomes observed had normal crossovers, as evidenced by centromeric distal LOH, which may indicate that defects in initial homolog pairing affect the final outcome. However, at the genome level, we cannot confidently investigate whether cells with biased even segregation have a similar number of crossovers to their reduction-biased counterparts, because we can detect all crossovers for reduction-segregating chromosomes, whereas crossovers in evenly segregating chromosomes can only be detected when the two recombined chromatids segregate apart (Figure 5B-C and Figure S16D, patterns 2 and 3). Assuming that recombined chromatids are equally likely to segregate together or separately, the number of crossovers is not smaller than in the case of their genome-level equal segregation, but we cannot exclude the possibility that segregation is biased away from 50 / 50 due to unresolved recombination intermediates (Figure 23, pattern 3).

[0260] Third, what is the outcome of these evenly segregated chromosomes? Do they return to mitosis with extensive LOH or do they proceed to MII and, if so, contribute to the formation of 1C gametes? In yeast, a phenomenon called "growth regression" has been characterized in which cells initiating a meiotic program can, in the presence of appropriate nutrients, return to normal mitosis resulting in numerous LOH events (Dayani et al., 2011). In human female meiosis, chromosomes with reverse segregation progress to MII, leading to one euploid oocyte and one euploid polar body 2, consistent with normal MII segregation. The authors suggest that unresolved recombination intermediates may have caused reverse segregation at MI by joining otherwise unrelated homologous chromatids, facilitating proper MII segregation (Figure 23, pattern 3) (Ottolini et al., 2015). Mlh1 is important for both mismatch repair (MMR) and the resolution of Holliday junction intermediates in meiosis. Considering the 2% sequence variance between B6 and Spret, it is possible that Mlh1 is limiting due to convergent MMR and there may not be enough Mlh1 to resolve recombination intermediates. However, we emphasize that this would not lead to LOH if recombined homolog chromatids cosegregate (Figure 5C). Therefore, M2 cells with LOH and equal segregation cannot be explained by cosegregation of unresolved intermediates.

[0261] Finally, in Figure 23, one of the patterns (Pattern 4) is not distinguishable from cells with crossovers but co-segregating recombined chromatids (Pattern 3), indicating a possible contribution to gamete formation from chromosomes without any interhomolog crossovers, probably due to insufficient homolog pairing. However, if these cells without crossovers contribute significantly to 1C cells, more chromosomes without crossovers should be observed among 1C cells. Among the 1C cells observed in both crossovers, the number of chromosomes with and without crossovers is roughly 50-50, indicating that they are mainly derived from some combination of patterns 1-3 in Figure 23, and that 2C cells without interhomolog crossovers (Patterns 4 and 5) do not substantially contribute to 1C cells that successfully complete MII.

[0262] Crossover hotnis and associated (epi)genomic factors Crossover hot spots are continuous bands and are shaped by many factors. Crossovers in the (B6 × Cast) cross correlate more strongly with meiotic DSB hot spots mapped in the F1 cross than in the individual maps for the two parental lines, which is expected based on previous findings that novel meiotic hot spots can form in F1 hybrids (Smagulova et al., 2016). In the (B6 × Spret) cross, crossovers correlate weakly but positively with Spo11 breaks. It should be noted that the Spo11 map only accounts for PRDM9 sites bound by the PRDM9 protein in the B6 allele, and it is possible that the Spret copy of PRDM9 binds to different sites to form new meiotic DSB hot spots that are not accounted for in our analysis. Genomic features observed to be positively correlated with meiotic crossing over include GC-rich regions (as in yeast meiosis (Petes, 2001; Petes and Merker, 2002)), lineages (Lilue et al., 2018), gene bodies, pseudogene transcripts, CTCF binding sites, replication domains (Marchal et al., 2018), DNA transposons, satellite DNA and an increase in CNVs among a subset of histone modifications including H3K4me1, H3K27me3 and H3K36me3 (Mu et al., 2017). Interestingly, binding sites for Dmrt6, involved in regulating the mitotic to meiotic transition in male germ cells (Zhang et al., 2014), are strongly correlated with meiotic crossing over hot niches. Notably, genomic features negatively correlated with meiotic crossing over include 3'UTRs, LINEs and low-complexity DNA. Unlike in yeast, where rDNA is highly resistant to meiotic crossing over (Petes and Botstein, 1977), mouse rDNA does not appear to suppress crossing over.With these genomic features, true meiotic crossover initiation sites can be distinguished from randomly picked zones in the mouse genome with accuracies of 0.73 and 0.85 for (B6 x Spret) and (B6 x Cast), respectively, with a predicted accuracy of 0.85 in (B6 x Cast) crosses applied to a subset of 25 genomic features. We emphasize that although the various features behave fairly consistently between modeling approaches, no causality can be assigned without further experimentation.

[0263] method Methods and molecular design of sci-L3-WGS and sci-L3-targeted-seq Single cell preparation and nucleosome depletion Cell suspensions are prepared by trypsinization from petri dishes or homogenization from tissues. Male F1 mice were euthanized by CO2 and subsequent cervical dislocation according to a protocol approved by the University of Washington IACUC. For isolation of male germ cells, epididymides were dissected by slicing the duct and incubating the tissue for 15 min at room temperature in 1 ml of 1x PBS supplemented with 10% FBS. After incubation, cell suspensions were collected by pipetting. Cells isolated from epididymides were used for (B6 x Spret) cross-breeding experiments and also as a source of mature sperm ("barcode group 3") for (B6 x Cast) cross-breeding. For isolation of nuclei from whole testes as an enrichment method for 2C cells for (B6 x Cast) cross-breeding, testicular cells were first cross-linked with 1% formaldehyde and nuclei were extracted using a hypotonic buffer. 1C and 2C nuclei were then FACS sorted primarily by DNA content based on the DAPI signal. Cultured human and mouse cells are pelleted at 550 g for 5 min at 4° C. and male germ cells are pelleted at 2400 g for 10 min at 4° C.

[0264] Nucleosome depletion mostly follows the xSDS method in sci-DNA-seq (Vitak et al., 2017), except that the lysis buffer is modified to fit the downstream LIANTI protocol (Chen et al., 2017). Cells are crosslinked by 406 μL of 37% formaldehyde (final concentration 1.5%) in 10 mL of DMEM complete medium for 10 min at room temperature (gently inverting the tube). Then, 800 μL of 2.5 M glycine is added and incubated on ice for 5 min. Cells are pelleted and washed with 1 mL of lysis buffer (60 mM Tris-Ac, pH 8.3, 2 mM EDTA, pH 8.0, 15 mM DTT). Resuspend the pellet in 1 mL of lysis buffer with 0.1% IGEPAL (I8896, SIGMA) and incubate on ice for 20 min. Nuclei are then pelleted, washed in 1xNEBuffer2.1, and resuspended in 800 μL of 1xNEBuffer2.1 with 0.3% SDS for nucleosome depletion at 42° C. (vigorously shaking at 500 rpm for 30 min). Then, 180 μL of 10% Triton-X is added and vigorously shaking (500 rpm) for 30 min at 42° C. Permeabilized nuclei are then washed twice in 1 mL of lysis buffer and resuspended in lysis buffer at 20,000 nuclei per μL.

[0265] Transposome design and assembly Transposon DNA oligos are synthesized with both 5' of the two strands phosphorylated, one required for Tn5 insertion as in LIANTI and Nextera (5' / phosphate / CTGTCTCTTATACACATCT, IDT, PAGE purified (SEQ ID NO:1)) and the other required for ligation (5' / phosphate / GTCTTG XXXXXXXX [1st barcode] AGATGTGTATAAGAGACAG, IDT, standard desalted (SEQ ID NO:2)). After 1:1 annealing in annealing buffer (10 mM Tris-HCl pH 8.0, 50 mM NaCl, 1 mM EDTA, pH 8.0) with stepwise cooling (95°C for 5 min, -0.1°C / cycle, 9 sec / cycle, 700 cycles to 25°C), the Tn5 duplex with 5' overhang is diluted to 1.5 μM. Next, 7.2 μL of storage buffer (1xTE with 50% glycerol) is added to 12 μL of 1 μM Tn5 transposase (Lucigen, TNP92110) and 0.79 μL of diluted transposase is incubated with 0.4 μL of 1.5 μM Tn5 duplex for 30 min at room temperature. The transposomes dimerize to a final concentration of 0.2 μM. The transposome complexes can be stably stored at -20 °C for up to one year. 24 reactions were set up for barcoding 24 wells in the first round, but more wells may be desirable depending on the application. For each new organism application, the transposomes are first further diluted to 0.1 μM for testing experiments. Although the number of unique reads and library complexity is less optimal (Figure 5), it can be used for mapping at lower resolution.

[0266] In Figure 7, the molecular structure of sci-L3-WGS at each step is shown. Commercial Nextera library preparation loses at least half of the sequenceable DNA material due to: 1) Tn5 insertion introduces symmetric transposon sequences at the two ends of the fragmented genomic DNA, which can result in the formation of hairpin loops when denatured, blocking PCR amplification; and 2) if the two ends are tagmented by both i5 or i7 with a 50% chance, the molecule cannot be sequenced. One key advantage of LIANTI over Nextera-based library preparation is that the looped Tn5 design breaks the symmetry introduced by the transposome dimer and promotes reverse transcription (RT) using an intramolecular RT primer, which is also a feature of looped transposons. However, looped transposons are not compatible with barcoding more than twice, which limits throughput and increases library costs considerably (see Table 2 for a comparison). With the modifications made to sci-L3-WGS, the advantages provided by looped Tn5 during the ligation step are maintained.

[0267] Tagmentation (1st barcode) and ligation (2nd barcode) Next, distribute 1.5 μL of nuclei into each well of a lo-bind 96-well plate at a concentration of 20,000 / μL and add 6.5 μL of H2O and 0.7 μL of 50 mM MgCl2 (EDTA in lysis buffer is at a final concentration of 3.24 mM). Add 1.2 μL of transposomes prepared above to each well, then incubate the plate at 55°C for 20 minutes (a thermomixer is recommended but not required). Then add 5 μL of stop solution (40 mM EDTA and 1 mM spermidine) and pool the nuclei into a trough. Add an additional 1 mL of lysis buffer to the nuclei suspension before pelleting. After carefully removing the supernatant, resuspend the nuclei in 312 μL of resuspension buffer (24 μL of 10 mM dNTPs, 48 ​​μL of 10× tagmentation buffer [50 mM MgCl2, 100 mM Tris-HCl pH 8.0], 96 μL of H2O, 144 μL of lysis buffer) and distribute 4.7 μL of the nuclei mix into each well of a new lo-bind 96-well plate. The hairpin ligation duplex (1. CAAGAC 2. Y'Y'Y'Y'Y'Y'Y' [reverse complement of 2nd round barcode] 3. CAGGAGCGAGCTGCATCCC 4. AATTTAATCGACTCACTATA 5. GGGATGCAGCTCGCTCCTG 6. YYYYYYY [2nd round barcode] (SEQ ID NO: 3)) is preannealed similarly to the Tn5 transposon duplex and diluted to 1.5 μM. Note that the ligation duplex contains five elements: 1) the reverse complement of the ligation adapter on Tn5; 2) the reverse complement of the second round barcode; 3) the reverse complement of the second strand synthesis (SSS) primer; 4) the T7 promoter, note that this is the hairpin loop region; 5) the second strand synthesis (SSS) primer region starting with GGG to enhance T7 transcription ("sp2" in Figure 4B); 6) the second round barcode ("bc2" in Figure 4B). Add 0.8 μL of these duplexes to each of the 64 wells with the nuclei suspension and add 1.18 μL of ligation mix (0.6 μL 10x NEB T4 ligase buffer, 0.48 μL PEG-4000, 0.1 μL T4 DNA ligase [Thermo EL0011]) to each well and incubate at 20 °C for 30 min.Note that after ligation, the looped structure resembles that of LIANTI, facilitating efficiency in the RT step (discussed below), and that barcodes from both rounds are present 3' of the T7 promoter and therefore included in the amplified molecule. The ligation reaction is stopped by adding 4 µL of stop solution. Cells are then pooled into a new trough (~630 µL), stained with DAPI at a final concentration of 5 µg / mL, and sorted 100-300 into each new well with 3 µL of lysis buffer added before cell sorting. Note that each sorting event by FACS involves ~3-5 nL of FACS buffer depending on the size of the nozzle; we recommend keeping the total amount of liquid added to each well <1 µL to keep the salt concentration low.

[0268] Linear amplification by cell lysis, gap extension and in vitro transcription Next, for cell lysis, a total of 3.5-4 μL of sorted nuclei in each well is taken up by incubating at 75 °C for 45 min, cooling to 4 °C and treating with freshly diluted Qiagen protease (final concentration 2 mg / mL) at 55 °C for 8 h. The protease is then heat-inactivated by incubating at 75 °C for 30 min. The cell lysate can be stored at -80 °C. It is recommended to process no more than 32 wells of samples (approximately 9600 single cells) for each experiment, as the subsequent amplification steps involve RNA and are time-sensitive. For gap extension (Figure 4C), a polymerase with strand displacement activity is used by adding a mixture of 2 μL HO, 0.7 μL 10x tagmentation buffer, 0.35 μL 10 mM dNTPs and 0.35 μL Bst WarmStart2.0 polymerase with strand displacement activity and incubate at 68 °C for 5 min. Note that if ligation is successful at both ends, the duplex is symmetric with T7 promoters on both sides, but if ligation is successful at only one end, the region in the dashed box is missing on one side. Ligation between molecules is generally inefficient. We included a pre-annealed hairpin loop to minimize the need for ligation between molecules, but two molecules (instead of three without the hairpin loop) still need to find each other. If the ligation efficiency is 50%, the rate of having ligation at both ends is 25%, but the rate of having ligation at either end is 75%. We show that in the second half of the RT step, successful ligation is only required at one end. After gap extension, assemble 20 μL of the T7 in vitro transcription system by adding 2 μL of H2O, 2 μL of T7 Pol mix, and 10 μL of rNMP mix (NEB, HiScribe™ T7 Quick High Yield RNA Synthesis kit). Incubate the mixture at 37 °C for 10-16 h.

[0269] RNA purification, RT and SSS (or targeted sequencing) The transcription is terminated by adding 2.2 μL of 0.5 M EDTA. The amplified RNA molecules are then purified with RCC-5 (Zymo Research, R1016) and eluted with 18 μL of 0.1×TE. A 30 μL RT system is assembled by first adding 0.6 μL of RNA RT primer (rArGrArUrGrUrGrUrArUrArArGrArGrArCrArG, IDT (SEQ ID NO: 4)), 2 μL of 10 mM dNTPs and 0.5 μL of SUPERase*In™ RNase inhibitor (20 U / μL, Thermo Fisher AM2696). It is then incubated at 70°C for 1 min and 90°C for 20 s for denaturation and removal of secondary structures, and then snap-cooled on ice. SuperScript™ IV Reverse Transcriptase (SSIV, Thermo Fisher 18090050) is used for RT with 6 μL 5×RT buffer, 1.5 μL 0.1 M DTT, 1 μL SUPERase*In™ and 1 μL SSIV. The RT reaction is incubated at 55° C. for 15 min, 60° C. for 10 min, 65° C. for 12 min, 70° C. for 8 min, 75° C. for 5 min and 80° C. for 10 min. The reaction is allowed to cool to room temperature before adding 0.5 μL RNase H (NEB) and 0.3 μL RNase A (Life Technologies, AM2270) and incubating at 37° C. for 30 min. Note that Figure 4E shows two scenarios during the RT step: 1) if ligation of both ends is successful, RT may be primed by the foldback loop as in LIANTI; 2) if ligation is successful at only one end, RT may be primed by the RNA RT primer added before the denaturation step. Excess RNA primer and RNA transcript are degraded after cDNA synthesis. Finally, synthesize the second strand with Q5 DNA polymerase by adding 27 μL HO, 20 μL 5x Q5 buffer, 20 μL Q5GC enhancer, 1 μL Q5 polymerase and 1 μL SSS primer (NNNN[UMI]ZZZZZZ[3rd barcode]GGGATGCAGCTCGCTCCTG, IDT, standard desalted (SEQ ID NO:5).The resulting double-stranded DNA can be purified with DCC-5 (ZymoResearch, D4014) and proceeded with a library preparation kit such as NEBNext Ultra II with a minimal three cycle PCR to add sequencing adapters.

[0270] It is worth noting that the SSS process can be easily modified to use a single cell barcode primer with a P5 end (AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGAC GCTCTTCCGATCT NNNNNNN ZZZZZZ [3rd barcode] GGGATGCAGCTCGCTCCTG (SEQ ID NO: 6)) together with a targeted primer for one region in the genome to enable targeted sequencing (Figure 3B). For example, in an application incorporating a lentivirus-based CRISPR library (Shalem et al., 2014), the guide RNA sequence in each single cell can be read using a P7 end with the lentivirus-incorporated CRISPR library primer, CAAGCAGAAGACGGCATACGAGATTCGCCTTG [index 1] GTGACTGGAGTTCAGACGTGTGCTCTTCCGATCTCCGACTCGGTGCCACTTTTTCAA (SEQ ID NO: 7)), thus avoiding the need to sequence the entire genome to enrich for a specific region of interest. In this case, the library preparation step can be omitted and replaced with gel or bead purification to remove primer dimers.

[0271] Method and molecular design of sci-L3-RNA / DNA co-assay Single cell preparation and nucleosome depletion Cell suspensions are prepared using the same protocol as in sci-L3-WGS, with the differences noted below. HEK293T, BJ-5ta and 3T3 cells were trypsinized from petri dishes and fixed with 2% PFA in 1x PBS at RT for 10 min at a cell concentration of 1 M / mL. Subsequent quenching (with glycine), washing, nuclei isolation (with 0.1% IGEPAL) and nucleosome depletion (xSDS method) steps are identical to sci-L3-WGS, except that 1% Superase-In is added to all lysis buffers and 1xNEBuffer2.1. Nuclei are resuspended in lysis buffer containing 1% Superase-In at 20,000 nuclei per μL.

[0272] Transposome and reverse transcription (RT) primer design For the single-cell genome amplification component, the design and assembly of the transposome is identical to sci-L3-WGS.

[0273] For the single-cell transcriptome profiling component, the reverse transcription primer shares a similar structure to sci-RNA-seq in (Cao et al., 2017; Cusanovich et al., 2015; Mulqueen et al., 2018; Ramani et al., 2017; Vitak et al., 2017) for the reverse transcription aspect, i.e., the poly-T priming portion of the oligo, but contains a different barcode structure and landing pad for the subsequent ligation step ( / 5Phos / GTCTTG[same landing pad sequence as for sci-L3-WGS]NNNNNN[UMI1 for tagging specific transcripts]X'X'X'X'X'X'X'X'X'X'X'X'[1st round barcode for transcriptome, different sequence from Tn5 transposon barcode]TTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTTT

[0274] RT and tagmentation (1st barcode), ligation (2nd barcode), FACS and cell lysis Next, to prepare for the RT step, distribute 1.5 μL of nuclei at a concentration of 20,000 / μL into each well of a lo-bind 96-well plate and add 0.2 μL of HO, 0.3 μL of 50 mM MgCl2 (to neutralize the EDTA in the lysis buffer), 0.25 μL of 10 mM dNTPs, and 1 μL of the 25 μM RT primer described above. The nuclei mixture is then incubated at 55 °C for 5 min to remove secondary structures and quickly inactivated on ice. Next, add 1 μL of 5x RT buffer, 0.03 μL of 100 mM DTT (note there is DTT from the lysis buffer, final concentration 5 mM), 0.25 μL of SSIV, 0.25 μL of RNaseOUT (Thermo Fisher catalog number 10777019) and incubate for 1 min at 25°C, 1 min at 37°C, 1 min at 42°C, 1 min at 50°C, and 15 min at 55°C for the RT reaction. Next, add 0.4 μL of MgCl2 and 3.52 μL of H2O to each well, and 1.2 μL of transposomes prepared above. All subsequent steps until after cell lysis are identical to sci-L3-WGS.

[0275] Linear amplification by gap extension and in vitro transcription Partial NEBNext read data 1 primer was used as 5' overhang (CACGACGCTCTTCCGATCT NNNNNNN (SEQ ID NO: 9)). For gap extension, random heptamers are used. Add 1 μL of 20 μM oligos, incubate at 95 °C for 3 min to denature the DNA, and gradually cool to room temperature for the oligos to anneal (~5 min). Then add 2 μL of H2O, 0.8 μL of 10x NEBuffer2, 0.4 μL of 10 mM dNTPs, 0.4 μL Klenow fragment (3'→5' exo-, NEB M0212S) and incubate at 30 °C for 8 min and 75 °C for 10 min. After gap extension, 20 μL of T7 in vitro transcription system is assembled by the same sci-L3-WGS protocol.

[0276] RNA purification, RT and SSS Except for the different oligo sequences, all steps are identical to sci-L3-WGS. Instead of using 0.6 μL of RNA RT primer in the RT step after IVT, 0.6 μL of NEBNext Read Data 1 primer (AATGATACGGCGACCACCG AGATCTACACTCTTTCCCTACACGACGCTCTTCCGATCT, P5 end for Illumina sequencing, IDT (SEQ ID NO: 10)) is used. For the SSS primer, AAGCAGAAGACGGCATACGAGAT[P7 end]NNNN[UMI2]Z'Z'Z'Z'Z'Z'[3rd barcode]CGTCTCTAC GGGATGCAGCTCGCTCCTG (SEQ ID NO: 11) is used to add the sequencing adaptor. It should be noted that the resulting double stranded DNA now contains both P5 and P7 ends for Illumina sequencing and can be purified with 1.1× AmpureXP beads and proceeded to sequencing. The library preparation step to add sequencing adapters and the minimal three-cycle PCR in sci-L3-WGS are not required for the co-assay.

[0277] Setting up the sci-L3-WGS experiment in the (B6 × Spret) and (B6 × Cast) crosses (B6×Spret) hybrid In two separate experiments, cells isolated from six and three epididymides from 70- and 88-day-old (B6 x Spret) F1 males, respectively, were pooled and fixed in 1% formaldehyde. In each experiment, after nucleosome depletion, 30,000 cells were distributed per well and in situ indexed Tn5 insertion was performed across 24 wells to add the first barcode. All cells were then pooled and these were redistributed into 64 wells to add the second barcode and T7 promoter by ligation. After pooling all cells again, the cell mixture was split 1:6 and the majority of cells (6 / 7) were FACS sorted and the remainder (1 / 7) were diluted. The resulting wells contained 100-360 cells per well with an estimated collision rate of 4-11%.

[0278] (B6×Cast) cross. From six testes, we recovered approximately 12M 1C round spermatozoa and approximately 0.5M 2C cells. However, due to the >20-fold higher number of 1C cells, many 1C cells were still found in the population sorted for 2C cells (Figure 8F). In one of the sci-L3-WGS experiments that attempted to enrich for 2C cells, we estimated that we tagmented approximately 160k spermatozoa, approximately 160k 1C round spermatozoa, and approximately 70k 2C cells from the epididymis, which were further enriched for 2C cells during the FACS step of sci-L3-WGS (Figure 8G). However, despite two rounds of enrichment, 1C cells still predominated.

[0279] Table 4. Oligos for sci-L3. [Table 4-1] [Table 4-2] [Table 4-3] [Table 4-4] [Table 4-5] [Table 4-6] [Table 4-7] [Table 4-8] [Table 4-9]

[0280] Bioinformatics and statistical analysis methods Read data processing, alignment and SNV calling Base calls were converted to fastq files by bcl2fastq, which allowed one mismatch for errors in the index. Next, to demultiplex, we used a customized shell script "sci_lianti_v2.sh" (Python script and R Markdown files are uploaded separately as "sci_lianti_inst.tar.gz"; the R package containing intermediate data files for generating all main and supplementary figures can be downloaded and installed through the following link: https: / / drive.google.com / file / d / 19NFubouHrahZ8WoblL-tcDrrTlIZEpJh / view?usp=sharing), which calls Python scripts or NGS tools for the following steps: 1) Order the read pairs such that all single-cell combinatorial barcodes are in read1 (R1); 2) Demultiplex the third (SSS, 6nt, no errors allowed) barcodes, name the reads with both barcodes and UMIs for the transcripts, and split the libraries by the third barcodes. It should be noted that for the third barcode split of individual libraries containing 100-300 single cells, all subsequent steps are performed in parallel; 3) use cutadapt to split the first (Tn5, 8 nt, one error allowed) and second (ligation, 7 nt, one error allowed) barcodes in R1, errors are calculated by Levenshtein distance, and the barcodes of both rounds are named reads. This step is performed in paired ends mode, i.e., if R1 does not have the correct barcode and spacer structure, paired reads 2 (R2) are rejected; 4) use cutadapt to clean up R2; 5) align to the hg19 or mm10 genome with bwa mem in paired ends mode (Li and Durbin, 2009).For experiments investigating barcode collisions, use the concatenated references of hg19 and mm10 and use the differentially aligned reads to determine the relative mapping rate to the human or mouse genome; 6) Split the bam files into single-cell bam files using the first and second round barcodes attached to the read names; 7) Convert the bam files to bed files by bedtool (Quinlan and Hall, 2010) and determine unique insertion sites when R1 or R2 share the same endpoint. Unique Tn5 insertion sites are defined as fragments where both ends of a read pair must be different; 8) Call variants in allele-aware mode using the “pileup” function in the “lianti” package (https: / / github.com / lh3 / lianti / blob / master / pileup.c) (Chen et al., 2017). Note that in this step, we only need that the depth threshold for each SNP position is exceeded in the bulk file for the SNP call to be included in the final vcf file, and therefore include the bulk bam file (generated by samtool merging of all ~6900 single cells (Chen et al., 2017; Li and Durbin, 2009), >30×) together with each single cell bam file so that the raw values ​​of the REF and ALT alleles are included in the single cell columns as long as the variant is present as a heterozygous SNP in the bulk file. This avoids the problem of high false negative rates due to low depth sequencing in single cells by converting the novel SNP calling question into a genotyping question;9) Annotate the SNVs called in each single cell for SNP quality by the reference SNP vcf file for Spret (SPRET_EiJ.mgp.v5.snps.dbSNP142.vcf.gz downloaded from the Mouse Genome Project). The annotated SNP file is then used as input for subsequent crossover breakpoint analysis.

[0281] HMM for calling breakpoints The genotype at a given SNP site is determined by comparing the number of reads supporting the reference and alternative alleles. In 1C cells, the crossover position is determined by fitting a hidden Markov model of three states, reference, alternative and heterozygote.

[0282] The transition matrix is ​​specified in Table 5.

[0283] Table 5. Transition matrix. [Table 5]

[0284] Parameters were manually selected based on a visual assessment of how well the HMM captures the apparent structure in the data and that the results do not change much when changing the primary parameters by two orders of magnitude. transprob takes a very small number [1e-10 / (total number of SNPs on a given chromosome) in this case] to reflect the belief that state transitions at any individual SNP site should be very rare events. Further decomposition of transprob by ratios of 0.3 and 0.7 aims to suppress fast successive transitions of the form reference-alternate-reference or alternative-reference-alternate.

[0285] The release matrix is ​​specified in Table 6.

[0286] Table 6. Release matrix. [Table 6]

[0287] After the hidden states are called for each individual SNP, consecutive long state blocks are called by removing state blocks shorter than 50 kb. Crossover positions are then determined by where the long state block switches to a different state, where the start position of the breakpoint zone is the last SNP position of the previous state block, and the zone end position is the first SNP position of the following state block.

[0288] In M2 cells, the mean allele frequency is first obtained by averaging over alleles in a window of 40 SNPs. The sorted allele frequencies are then used to infer the underlying chromosomal state from a hidden Markov model with a single Gaussian probability distribution.

[0289] The transition matrix is ​​specified in Table 7.

[0290] Table 7. Transition matrix. [Table 7]

[0291] The release matrix is ​​specified in Table 8.

[0292] Table 8. Release matrix. [Table 8]

[0293] Long continuous state blocks are called by removing state blocks shorter than 50 kb, and then approximate breakpoint locations are determined by where the long state blocks switch to different states. The approximate breakpoint locations are then refined by a likelihood ratio test that aims to find possible breakpoints among the 20 upstream and 20 downstream SNPs surrounding the approximate breakpoint. For each SNP, the probability of observing the observed genotype is specified in Table 9.

[0294] Table 9. Probability of observing observed genotypes. [Table 9]

[0295] The error-prob is defined as 1e-3, which reflects the probability that a SNP is called incorrectly. For each SNP around the approximate breakpoint, the likelihood that it is the actual breakpoint is calculated by the distribution above. All SNPs with a likelihood greater than 0.01*maximum likelihood are considered to be within the breakpoint. The start of the break zone is determined as the leftmost SNP among these SNPs, and the end of the break zone is determined as the rightmost SNP. As in the 1C case, all M2 cell breakpoint zones are further manually inspected to remove artifacts, e.g., where two immediately adjacent switches are within 50 kb. We also performed the same breakpoint calling in mitotically dividing Patski cells. For M2 and Patski cells, the breakpoint zones were further manually inspected by comparing bin sizes of 10 and 40 SNPs for cells with sparse genome coverage.

[0296] This step generates crossover breakpoints, which are then post-processed to add chromosome segregation information based on whether the centromeric regions, i.e., the start regions of each chromosome, are heterozygous ("mt", mitotic segregation) or homozygous ("me", meiotic segregation).

[0297] Parthenogenetic chromosome analysis This step takes the rds file from the HMM output and generates parthenogenetic chromosome calls.

[0298] Analysis of meiotic crossover and chromosome segregation at the chromosome level This process produces the chromosome-level characterization of the meiotic crossovers shown in FIGS.

[0299] Fitting a finite mixture model to barcoded population 2C cells of the (B6 × Cast) cross We fit the data to a mixture of three binomial distributions, parameterized by p1, p2, and p3, respectively, which represent their probabilities of chromosomes segregating evenly. The relative contribution of these three binomial distributions is represented by the length-3 vector theta. We estimate p1, p2, p3 as well as θ by drawing samples from their posterior distributions using the R package rstan (http: / / mc-stan.org / users / interfaces / rstan) with a uniform Dirichlet leading for θ: θ~Dir(K=3, α=1) and a beta leading for p: p~beta(a=5, b=5). See the Stan file mt_mixture_model.stan for further details on the model specification.

[0300] Preprocessing of datasets from other genomic studies to build linear models of cross-hotline and cell clustering. We processed datasets from previous genome studies, as well as downloaded mouse annotation files in gff3 format and RepeatMasker (https: / / genome.ucsc.edu / cgi-bin / hgTables) from the UCSC genome browser for various genomic elements. Datasets based on mm9 are first raised to mm10. These datasets are roughly divided into two categories: numerical data in bed format or signals of various genetic or epigenetic marks in bedGraph format. In cell clustering and predictive modeling, the crossing zones have different lengths. For cell clustering analysis, we normalize the numerical data by dividing the total amount of sequence summed from all crossings in each single cell, and for each crossing zone or randomly picked zone, we normalize by dividing the zone length plus 1 kb so that very short zones are not overweighted. It should be noted that the median zone length is 150 kb so that adding 1 kb does not include a lot of extra sequence. For datasets with continuous signals of various marks, we take the average signal of the marks across the crossing or random zone. For the cross pile-up dataset, uniformly sized 100 kb windows were used, so zone lengths were not normalized when using the numerical data.

[0301] In addition to the datasets noted in the Discussion section, where features have a statistically significant association with the presence of crossovers, we also used the following datasets: 1) sequence variance (Lilue et al., 2018); 2) ATAC-seq and H3K27ac mapped from purified pachytene spermatocytes (Maezawa et al., 2018); 3) bisulfite sequencing from spermatogonia (Inoue et al., 2017); 4) MNase-based nucleosome positioning in spermatocytes (Barral et al., 2017); 5) butyrylation and acetylation of H4K5 and H4K8 in spermatocytes (Goudarzi et al., 2016); 6) H2A ubiquitination in spermatocytes (Hasegawa et al., 2015); 7) the binding site of CTCFL, a testis-specific paralog of the binding site of CTCF (Sleutels et al., 2012); 8) 5-hmC maps in pachytene spermatocytes (Gan et al., 2013); 9) terminal-seq after etoposide treatment and CTCF and RAD21 ChIP-seq in activated B cells, TOP2A and TOP2B ChIP-seq in MEFs (Canela et al., 2017); 10) ATAC-seq data of the Patski allele (Bonora et al., 2018).

[0302] PCA for cell clustering, BMA for linear models of cross-sectional hot spots, and Random Forest for predictive models of cross-sectional and random zones. We use principal component analysis to visualize in 2D the separation of 1C and M2 cells based on their breakpoint features. For each single cell, we collected the crossover-related information of a total of 78 features corresponding to three types. As the first type, we simply calculated the number of crossover or whole-chromosome LOH events for each chromosome in each cell. As the second type, we calculated the median values ​​of crossover breakpoints in each cell for features such as GC content, sequence variance, and chromatin mark intensity. As the third type, we calculated the normalized values ​​of genomic elements, such as gene bodies, long terminal repeats (LTRs), and LINE elements overlapping with crossover breakpoints in each cell.

[0303] To construct a linear model predicting the cross hot varnish, we used Bayesian model averaging using the “bas” package (Clyde et al., 2011) (function bas.lm with default settings uses 2 14 A random forest is trained to distinguish true crossover zones from zones randomly drawn from the genome that resemble a "null" distribution. Model accuracy is determined by fully nested 5-fold cross-validation, external 5-fold and 5-fold within each training set (see the section called "Model" in sci-L3-WGS-figures.Rmd for R code and annotations).

[0304] To estimate the effect of lineage (or cell type) on positioning the rightmost crossover along the chromosomes, we use a linear mixed-effects model with a fixed effect on lineage (or cell type) and a random break on chromosomes to account for interchromosomal variability (see the section called "Karyotype plots" in sci-L3-WGS-figures.Rmd for R code and annotations).

[0305] References [Table B-1] [Table B-2] [Table B-3] [Table B-4] [Table B-5]

[0306] The complete disclosures of all patents, patent applications and publications cited herein, as well as electronically available materials (including, for example, nucleotide sequence submissions in GenBank and RefSeq, and amino acid sequence submissions in, for example, SwissProt, PIR, PRF, PDB, and translations from annotated coding regions in GenBank and RefSeq) are fully incorporated by reference. Supplementary materials referenced in publications (such as supplementary tables, supplementary figures, supplementary materials and methods, and / or supplementary experimental data) are likewise fully incorporated by reference. In the event of any inconsistency between the disclosure of this application and the disclosure(s) of any document incorporated herein by reference, the disclosure of this application shall prevail. The foregoing detailed description and examples are given solely for clarity of understanding. No unnecessary limitations should be understood therefrom. The disclosure is not limited to the exact details shown and described, since variations obvious to one skilled in the art are included within the disclosure as defined by the claims.

[0307] Unless otherwise indicated, all numbers expressing amounts of components, molecular weights, and the like used in the specification and claims should be understood to be modified in all instances by the term "about." Accordingly, unless otherwise indicated to the contrary, the numerical parameters set forth in the specification and claims are approximations that may vary depending on the desired properties sought to be obtained by the present disclosure. At the very least, and without intending to limit the doctrine of equivalents to the scope of the claims, each numerical parameter should at least be construed in light of the number of reported significant digits and by applying ordinary rounding techniques.

[0308] Notwithstanding that the numerical ranges and parameters setting forth the broad scope of the disclosure are approximations, the numerical values ​​set forth in the specific examples are reported as precisely as possible, however, all numerical values ​​inherently have certain ranges necessarily resulting from the standard deviations found in their respective testing measurements.

[0309] Unless otherwise specified, all headings are for the convenience of the reader and should not be used to limit the meaning of the text that follows the heading.

Claims

1. 1. A method for preparing a sequencing library comprising nucleic acids from a plurality of cells, the method comprising: (a) providing isolated nuclei from a plurality of cells; (b) distributing the isolated nuclei into a first plurality of compartments, wherein each compartment comprises a subset of the isolated nuclei, and contacting each subset with a transposome complex, wherein the transposome complex in each compartment comprises a transposase, a first indexing sequence that is different from first indexing sequences in other compartments, and at least one ubiquitous sequence; (c) fragmenting nucleic acid in said subset of nucleosome-depleted nuclei into a plurality of nucleic acid fragments and incorporating said first index sequence into at least one strand of said nucleic acid fragments to generate indexed nuclei comprising indexed nucleic acid fragments; (d) combining the indexed nuclei to generate a pooled indexed nuclei; (e) distributing the pooled indexed subset of nuclei into a second plurality of compartments; (f) incorporating a second index sequence into the indexed nucleic acid fragments in each compartment to generate doubly indexed nucleic acid fragments, wherein the second index sequence is unique to each compartment; (g) combining the doubly indexed nuclei to generate pooled doubly indexed nuclei; (h) distributing a subset of the pooled doubly-indexed nuclei into a third plurality of compartments, lysing the doubly-indexed nuclei, and amplifying the doubly-indexed nucleic acid fragments by linear amplification; (i) incorporating a third index sequence into the doubly indexed nucleic acid fragments in each compartment to generate triply indexed nucleic acid fragments, wherein the third index sequence in each compartment is different from the first and second index sequences in other compartments, and the third index sequence in each compartment is different from the third index sequence in other compartments; (j) combining said triply indexed nucleic acid fragments to generate pooled triply indexed nucleic acid fragments, thereby producing a sequencing library from a plurality of said nuclei; and

2. The method of claim 1 , wherein the isolated nuclei are free of nucleosomes.

3. 2. The method of claim 1, further comprising subjecting the isolated nuclei prior to step (b) to conditions that produce nucleosome-depleted nuclei while maintaining the integrity of the isolated nuclei.

4. The method of claim 1 , wherein the ubiquitous sequence comprises a ubiquitous primer.

5. 2. The method of claim 1, wherein the indexed nucleic acid fragments of step (c) remain attached to the transposase such that nucleic acid fragments derived from the same genomic DNA molecule remain physically linked.

6. The method of claim 5, further comprising dissociating the transposase from the indexed nucleic acid fragment.

7. 2. The method of claim 1, wherein step (f) comprises contacting each subset with a hairpin ligation duplex under conditions suitable for ligation of the hairpin ligation duplex to one or both ends of the indexed nucleic acid fragment to provide a doubly-indexed nucleic acid fragment, wherein the hairpin ligation duplex comprises the second index sequence.

8. 2. The method of claim 1, further comprising manipulating the doubly indexed nucleic acid fragment after step (h) and before step (i), wherein the manipulating comprises gap extension, in vitro transcription, reverse transcription, or a combination thereof.

9. 2. The method of claim 1, wherein the incorporating of step (i) comprises second strand synthesis.

Images