Improved methods and kits for generating DNA libraries for massively parallel sequencing - Patents.com

Abstract

A method for generating a massively parallel sequencing library, comprising: a) generating a WGA library; a primary WGA library containing fragments containing universal sequence adapters; The disclosed method includes the steps of: preparing a DNA library (pWGAlib); b) performing one PCR cycle on the pWGAlib using a first primer (1PR) comprising, in 5' to 3' order, a first sequencing adapter (1PR5SA) and a first primer 3' section (1PR3S) that hybridizes to the reverse complement of the WGA library universal sequence adapter; c) performing one PCR cycle on the product of step b) using a second primer (2PR) comprising, in 5' to 3' order, a second sequencing adapter (2PR5SA) different from 1PR5SA and a second primer 3' section (2PR3S) that hybridizes to the reverse complement of the WGA library universal sequence adapter; and d) amplifying the product of step c) by PCR using a third primer comprising 1PR5SA and a fourth primer comprising 2PR5SA. [Selected Figure] Figure 2A

Description

[Technical field]

[0001] [CROSS REFERENCE TO RELATED APPLICATIONS] This application claims priority to European Patent Application No. 17182693.6, filed July 21, 2017, the disclosure of which is incorporated herein by reference.

[0002] The present invention relates to a method and a kit for generating massively parallel sequencing libraries for whole genome sequencing from whole genome amplification products (WGA). In particular, said method can be advantageously applied to deterministic restriction site whole genome amplification (DRS-WGA) DNA products.

[0003] The libraries can be advantageously used for low-pass whole genome sequencing and genome-wide copy number profiling. [Background technology]

[0004] To simplify and / or enable performing various types of genetic analyses on single cells, including sequencing, SNP detection, etc., it is useful to perform whole genome amplification (WGA) to obtain more DNA.

[0005] WGA by LM-PCR based on definitive restriction enzyme sites (e.g., as described in patent document 1) is known from the art (hereinafter simply referred to as DRS-WGA). A commercial kit for LM-PCR-based DRS-WGA (Ampli1™ WGA kit, Menarini Silicon Biosystems) was used in non-patent document 1. In this study, copy number analysis by low-pass whole genome sequencing was performed on single-cell WGA material. However, in the standard workflow used in this paper, Illumina library preparation requires multiple steps, including i) digestion of WGA adapters, ii) DNA fragmentation, iii) EndRepair, iv) A-tailing, v) ligation of barcoded adapters, and vi) standard Illumina workflow steps such as the general steps of sample pooling for barcoded NGS libraries, and vii) sequencing. As shown in the aforementioned paper (Figure 5b), WBCs showed a small number of copy number calls that were suspected to be false positives, while typical CTCs showed more deviations.

[0006] Ampli1™ WGA is compatible with array comparative genomic hybridization (aCGH) and has indeed been shown by several groups (Non-Patent Document 2, Non-Patent Document 3) to be suitable for high-resolution copy number analysis. However, aCGH technology is expensive and laborious, so a different method such as low-pass whole genome sequencing (LPWGS) may be preferred for somatic copy number alteration (CNA) detection.

[0007] Although DRS-WGA provides the best results in terms of homogeneous and balanced amplification, current protocols based on aCGH or metaphase CGH are laborious and / or expensive. Low-pass whole genome sequencing has been presented as a high-throughput method to analyze multiple samples with higher throughput and lower cost than aCGH. However, known methods for the generation of massively parallel sequencing libraries for WGA products (such as DRS-WGA) still require protocols that include several enzymatic steps and reactions.

[0008] A more streamlined method that combines the reproducibility and quality of DRS-WGA with the ability to analyze genome-wide copy number variations (CNVs) is desirable. In addition, it is desirable to determine whole-genome copy number profiles from small amounts of cell, FFPE, or tissue biopsy samples.

[0009] WO 2005 / 023991 in the name of the applicant discloses a method for the efficient generation of massively parallel sequencing libraries, also called NGS (Next Generation Sequencing) libraries, starting from WGA products.

[0010] The method comprises amplifying a primary WGA DNA library with two primers, each carrying a different sequencing adapter at the 5' end that allows sequencing on a specific sequencing platform. Sequencing platforms that can be used are, for example, the Ion Torrent platform or the Illumina platform.

[0011] When using certain sequencing platforms, for example the Illumina platform, it is particularly advantageous to select library fragments that contain two different sequencing adapters (P5 and P7 in the case of Illumina) at opposite ends of the fragment. These fragments are referred to as "heteroadapter fragments". To this end, one of the embodiments disclosed in the above mentioned application provides that one of the two primers for amplifying the primary WGA DNA library is biotinylated at the 5' end. This embodiment is summarized in Figure 1. Once the primary WGA DNA library is amplified with the two primers, the fragments are selected with streptavidin beads. Fragments with identical sequencing adapters at both ends (hereinafter referred to as "homoadapter fragments") are either eluted (if not biotinylated) or remain bound to the streptavidin beads (if both ends are biotinylated), while ssDNA heteroadapter fragments (with different sequencing adapters at both ends) are denatured and eluted for selection.

[0012] Selection with biotinylated primers has several drawbacks. In particular, it makes single-stranded DNA libraries less quantifiable and less suitable for storage. In addition, for Illumina sequencing workflows, the use of double-stranded DNA libraries is preferred. This issue can be addressed by synthesizing the second strand or performing additional amplification cycles with P5 and P7 primers, but this of course makes the method somewhat more complicated and single-tube reactions impossible. As a result, the design of the kit becomes more complicated, for example by requiring the use of special buffers designed for the above purpose. [Prior art documents] [Patent documents]

[0013] [Patent Document 1] International Publication No. 2000 / 017390 [Patent Document 2] International Application No. PCT / EP2017 / 059075 [Non-patent literature]

[0014] [Non-Patent Document 1] Hodgkinson CL et al., Tumorigenicity and genetic profiling of circulating tumor cells in small-cell lung cancer, Nature Medicine 20, 897-903 (2014) [Non-Patent Document 2] Moehlendick B, et al. (2013) A Robust Method to Analyze Copy Number Alterations of Less than 100 kb in Single Cells Using Oligonucleotide Array CGH. PLoS ONE 8(6): e67031 [Non-Patent Document 3] Czyz ZT, et al (2014) Reliable Single Cell Array CGH for Clinical Samples. PLoS ONE 9(1): e85907 Summary of the Invention

[0015] It is therefore an object of the present invention to provide a method for the generation of massively parallel sequencing libraries starting from WGA products, which overcomes the above mentioned problems.

[0016] Another object of the present invention is to provide a method for low-pass whole genome sequencing and a method for genome-wide copy number profiling starting from a WGA product and using the library preparation method according to the present invention.

[0017] It is a further object of the present invention to provide a massively parallel sequencing library preparation kit and a low-pass whole genome sequencing kit for carrying out the above-mentioned method. [Brief description of the drawings]

[0018] [Figure 1] FIG. 1 is a diagram summarizing the steps of the method for selecting heterozygous adapter fragments disclosed in Patent Document 2 by the applicant. [Figure 2A] FIG. 1 summarizes the steps of the method according to the invention, specifically for the situation where the primary WGA DNA library is obtained by DRS-WGA, without, however, intending to limit the scope of the invention to this particular type of WGA. [Figure 2B] Figure 2B shows the structure of the primary WGA DNA library and the four primers used in the method according to the invention shown in Figure 2A. In the claims and the specification, the acronyms used for the various segments of the primary WGA DNA library and the four primers are also shown in brackets. [Figure 3A] FIG. 1 shows dual index sequencing on MiSeq, HiSeq2000 / 2500 and 1000 / 1500 using one custom sequencing primer according to the present invention. [Figure 3B] FIG. 1 shows dual index sequencing on MiniSeq, NextSeq, and HiSeq3000 / 4000 Illumina sequencing platforms using two custom sequencing primers according to the present invention. [Figure 4A] Figure 1 shows the sequencing results of low-pass whole genome sequencing performed by the method according to the invention, showing the copy number alteration (CNA) profiles of two single cells belonging to NCI-H441 and SW-480 cell lines separated by DEPArray™ (Menarini Silicon Biosystems). [Figure 4B]Figure 1 shows the sequencing results of low-pass whole genome sequencing performed by the method according to the invention. Copy number alteration (CNA) profiles of single blood cells (circulating tumor cells (CTCs) and white blood cells (WBCs)) separated by DEPArray™ (Menarini Silicon Biosystems) are shown. [Diagram 5] Figure 1 shows the sequencing results of low-pass whole genome sequencing performed by the method according to the invention, showing copy number alteration (CNA) profiles of single tumor cells and 50 stromal cell populations originating from disaggregated FFPE sections, numerically separated by DEPArray™ (Menarini Silicon Biosystems) and whole genome amplified using the Ampli1™ WGA kit. [Figure 6A] Figure 1 shows inter-method and inter-platform (Ion Torrent and Illumina) comparison results, specifically showing copy number alteration (CNA) profiles of NCI-H23 single cells obtained with a low-pass whole genome sequencing method (as shown in US Patent No. 5,999,336) for Ion Torrent and the method according to the present invention for the Illumina platform. [Figure 6B] FIG. 1 shows inter-method and inter-platform (Ion Torrent and Illumina) comparison results, specifically showing hierarchical clustering of NCI-H441 and WBC single cells based on CNA profiles. DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

[0019] definition Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this invention belongs. Although a number of methods and materials similar or equivalent to those described herein can be used in the practice or testing of the present invention, the preferred methods and materials are described below. Unless otherwise noted, the techniques described herein used in connection with the present invention are standard techniques known to those of ordinary skill in the art.

[0020] The term "original DNA" intends genomic DNA (gDNA) prior to amplification by DRS-WGA.

[0021] The term "adapter" or "WGA adapter" or "WGA PCR primer" or "WGA library universal sequence adapter" refers, in the case of DRS-WGA, to an additional oligonucleotide that is ligated to each fragment generated by the action of a restriction enzyme, or, in the case of MALBAC, to a known polynucleotide sequence present in the 5' section of each molecule of a WGA DNA library, resulting from the extension and PCR process.

[0022] The term "copy number alteration (CNA)" intends a change in the copy number of a genomic region in somatic cells, typically defined relative to the genome of the same individual.

[0023] The term "copy number variation (CNV)" refers to a germline variation in the copy number of a genomic region, typically defined relative to a reference genome. Since most arguments can be applied to both situations, CNA and CNV may be used interchangeably in the specification. Unless otherwise specified, each of these terms is intended to refer to both situations.

[0024] The term "massively parallel sequencing (MPS)" or "next-generation sequencing (NGS)" refers to a method of DNA sequencing that involves generating a library of clonally sequenced (pre-clonally amplified or not) DNA molecules separated in space and / or time. Examples include the Illumina platform (Illumina), the Ion Torrent platform (ThermoFisher Scientific), the Pacific Biosciences platform, and MinION (Oxford Nanopore Technologies).

[0025] The term "target sequence" intends a region of interest in the original DNA.

[0026] The term "primary WGA DNA library (pWGAlib)" intends a DNA library obtained from a WGA reaction.

[0027] The term "multiple annealing and looping amplification cycles (MALBAC)" refers to a quasi-linear whole genome amplification method (Zong et al., Genome-wide detection of single-nucleotide and copy-number variations of a single human cell, Science. 2012 Dec 21;338(6114):1622-6. doi: 10.1126 / science.1229164.). MALBAC primers have an 8-nucleotide 3' random sequence that hybridizes to the template, and a 27-nucleotide 5' consensus sequence (GTG AGT GAT GGT TGA GGT AGT GTG GAG). After the first extension, the semi-amplicon is used as a template for a second extension to obtain a full amplicon with complementary 5' and 3' ends. After several cycles of quasi-linear amplification, the full amplicon can be exponentially amplified by subsequent PCR cycles.

[0028] The term "DNA library purification" refers to a process of separating DNA library material from undesired reaction components such as enzymes, dNTPs, salts and / or other molecules that are not part of the desired DNA library. Examples of DNA library purification processes include purification by paramagnetic bead-based techniques (often referred to as "magnetic beads" for convenience in the scientific literature and below), such as Agencourt AMPure XP, or Solid Phase Reversible Immobilization (SPRI) beads by Beckman Coulter, or spin column purification, such as Amicon spin columns by Merck Millipore. Other examples of DNA library purification processes include purification by magnetic beads conjugated directly or via protein-protein interactions with oligonucleotide baits, such as streptavidin-coated magnetic beads that interact with biotinylated oligonucleotides.

[0029] The term "DNA library selection" intends a process in which either or both of DNA library purification or DNA library size selection are performed.

[0030] The term "sequencing adapter (SA)" contemplates one or more molecules that provide a means for sequencing a DNA insert. Each molecule may contain none, one or more of a polynucleotide sequence, a functional group, or both. In particular, it contemplates a non-information carrying polynucleotide sequence that is required to be present in a massively parallel sequencing library in order for the sequencer to generate an accurate output sequence (non-limiting examples include polynucleotide sequences that hybridize ssDNA to a flow cell in Illumina sequencing, or to an ion sphere in Ion Torrent sequencing, or polynucleotide sequences required to prime a sequencing by synthesis reaction).

[0031] The term "sequencing barcode" intends a polynucleotide sequence that, when sequenced within a sequencer read, allows the assignment of the read to a particular sample associated with that barcode.

[0032] The term "low-pass whole-genome sequencing" intends whole-genome sequencing in which the average sequencing depth is less than one.

[0033] The term "average sequencing depth" as used herein refers to the total number of bases sequenced and mapped to a reference genome per sample divided by the size of the entire reference genome. The total number of bases sequenced and mapped can be approximated by multiplying the number of mapped reads by the average read length.

[0034] By "homogenization" is intended the operation of adjusting and equalizing the concentration of one or more samples.

[0035] "Normalization" contemplates the operation of adjusting the concentrations of one or more samples to a desired ratio between them (equalization is a special case where the ratio is 1). For convenience in the specification, the terms normalization and equalization are used interchangeably, since they are clearly conceptually identical.

[0036] Detailed Description of the Invention Referring to FIG. 2A, which is an example in which the primary WGA DNA library is obtained by DRS-WGA and the sequencing platform is an Illumina platform, the massively parallel sequencing library generation method according to the present invention includes the following steps:

[0037] In step a, a primary WGA DNA library (pWGAlib) is prepared, which comprises a fragment comprising a known 5' sequence section (5SS), a central sequence section (MSS), and a known 3' sequence section (3SS) that is reverse-complementary to the known 5' sequence section. The known 5' sequence section (5SS) comprises a WGA library universal sequence adaptor. The central sequence section (MSS) comprises at least an insert section (IS) that corresponds to the DNA sequence of the original unamplified DNA prior to WGA. The central sequence section (MSS) optionally comprises, in addition to the insert section (IS), a flanking 5' middle section (F5) and / or a flanking 3' middle section (F3) (e.g., when the primary WGA DNA library is generated by MALBAC or DRS-WGA according to the teachings of WO 2015 / 118077). In Figures 2A and 2B, the 5' known sequence section (5SS) corresponds to the LIB sequence specific for DRS-WGA (SEQ ID NO:50), and the 3' known sequence section (3SS) is the reverse complement of the LIB sequence (LIBrc).

[0038] In step b, one PCR cycle is carried out on the primary WGA DNA library using at least one first primer (1PR), the first primer comprising at least a first primer 5' section (1PR5S) and a first primer 3' section (1PR3S). The first primer 5' section (1PR5S) comprises at least one first sequencing adaptor (1PR5SA). The first primer 3' section (1PR3S) hybridizes to the known 3' sequence section (3SS). The one PCR cycle of step b results in a first primer-extended WGA DNA library.

[0039] One PCR cycle includes a double-stranded DNA denaturation step, a primer annealing step, and an extension step of the annealed primers. A preferred embodiment includes a denaturation step at 95° C. for 30 seconds, an annealing step at 62° C. for 30 seconds, and an extension step at 72° C. for 3 minutes.

[0040] Preferably, the first primer (1PR) further comprises at least one lead sequencing primer sequence (1PRSEQ) at the 3' position of the first primer 5' section (1PR5S) and at the 5' position of the first primer 3' section (1PR3S).

[0041] The first primer (1PR) preferably has a sequence selected from the group consisting of SEQ ID NO:1 to SEQ ID NO:12.

[0042] After step b, it is necessary to prevent the first primer (1PR) from polymerizing other primary WGA products, i.e. the WGA DNA library fragments extended with the first primer. In a preferred embodiment, the first primer-extended WGA DNA library is purified after step b. This purification is preferably performed with SPRIselect beads (Beckman Coulter).

[0043] In step c, one PCR cycle is carried out on the first primer-extended WGA DNA library using at least one second primer (2PR), the second primer comprising a second primer 5' section (2PR5S) and a second primer 3' section (2PR3S). The second primer 5' section (2PR5S) comprises at least one second sequencing adaptor (2PR5SA) that is different from the at least one first sequencing adaptor (1PR5SA). The second primer 3' section (2PR3S) hybridizes to the known 3' sequence section (3SS).

[0044] One PCR cycle in step c results in a WGA DNA library extended with the first and second primers.

[0045] One PCR cycle includes a double-stranded DNA denaturation step, a primer annealing step, and an extension step of the annealed primers. A preferred embodiment includes a denaturation step at 95° C. for 30 seconds, an annealing step at 60° C. for 30 seconds, and an extension step at 72° C. for 3 minutes.

[0046] The second primer (2PR) preferably has a sequence selected from the group consisting of SEQ ID NO:13 to SEQ ID NO:20.

[0047] After step c, the WGA DNA library extended with the first and second primers is purified, preferably with a 2.5M NaCl PEG 20% solution.

[0048] In step d, the WGA DNA library extended with the first and second primers is amplified using at least one third primer (3PR) containing a first sequencing adaptor (1PR5SA) and at least one fourth primer (4PR) containing a second sequencing adaptor (2PR5SA). The PCR amplification in step d results in an amplified WGA DNA library extended with the first and second primers. The yield of this amplification step for the DNA library is sufficient to perform its sequencing.

[0049] After step d, the amplified WGA DNA library extended with the first and second primers is purified, preferably with a 2.5M NaCl PEG 20% solution.

[0050] The third primer (3PR) preferably has the sequence of SEQ ID NO:22 and the fourth primer (4PR) preferably has the sequence of SEQ ID NO:21.

[0051] The first primer 5' section (1PR5S) of the first primer (1PR) preferably further comprises at least one first sequencing barcode (1PR5BC) at the 3' position of at least one first sequencing adaptor (1PR5SA) and at the 5' position of the first primer 3' section (1PR3S). The second primer 5' section (2PR5S) of the second primer (2PR) preferably further comprises at least one second sequencing barcode (2PR5BC) at the 3' position of at least one second sequencing adaptor (2PR5SA) and at the 5' position of the second primer 3' section (2PR3S). This allows for higher multiplexing. Specifically, if SEQ ID NO:1 to SEQ ID NO:12 are used as first primers (1PR), each containing a different barcode, and SEQ ID NO:13 to SEQ ID NO:20 are used as second primers (2PR), each containing a different barcode, then 96 barcode combinations can be obtained and 96 libraries can be analyzed.

[0052] The WGA library universal sequence adaptor is preferably a DRS-WGA library universal sequence adaptor or a MALBAC library universal sequence adaptor, and more preferably a DRS-WGA library universal sequence adaptor.

[0053] The DRS-WGA library universal sequence adaptor preferably has SEQ ID NO:50, and the MALBAC library universal sequence adaptor preferably has SEQ ID NO:51.

[0054] The method for low-pass whole genome sequencing according to the present invention comprises the following steps.

[0055] First, prepare multiple barcoded massively parallel sequencing libraries obtained by the above-disclosed method for generating a massively parallel sequencing library, and pool the samples obtained using different sequencing barcodes (BC).Then, sequence the pooled library.

[0056] Pooling samples using different sequencing barcodes (BC) further comprises quantifying the DNA in each of the barcoded, massively parallel sequencing libraries and normalizing the amounts of the barcoded, massively parallel sequencing libraries.

[0057] The massively parallel sequencing library preparation kit according to the present invention includes at least one first primer (1PR), one second primer (2PR), one third primer (3PR) and one fourth primer (4PR), the structures of which have been disclosed above.

[0058] In a preferred embodiment, in which the primary WGA DNA library is DRS-WGA, the massively parallel sequencing library preparation kit comprises one or more primers selected from the group consisting of SEQ ID NO:1 to SEQ ID NO:12, one or more primers selected from the group consisting of SEQ ID NO:13 to SEQ ID NO:20, and primers of SEQ ID NO:21 and SEQ ID NO:22.

[0059] In an alternative preferred embodiment, where the primary WGA DNA library is DRS-WGA, the massively parallel sequencing library preparation kit comprises one or more primers selected from the group consisting of SEQ ID NO:52 to SEQ ID NO:63, one or more primers selected from the group consisting of SEQ ID NO:13 to SEQ ID NO:20, and primers of SEQ ID NO:21 and SEQ ID NO:22.

[0060] In an alternative preferred embodiment, where the primary WGA DNA library is MALBAC WGA, the massively parallel sequencing library preparation kit comprises one or more primers selected from the group consisting of SEQ ID NO:27 to SEQ ID NO:38, one or more primers selected from the group consisting of SEQ ID NO:39 to SEQ ID NO:46, and primers of SEQ ID NO:21 and SEQ ID NO:22.

[0061] In an alternative preferred embodiment, where the primary WGA DNA library is MALBAC WGA, the massively parallel sequencing library preparation kit comprises one or more primers selected from the group consisting of SEQ ID NO:64 to SEQ ID NO:75, one or more primers selected from the group consisting of SEQ ID NO:39 to SEQ ID NO:46, and primers of SEQ ID NO:21 and SEQ ID NO:22.

[0062] A low-pass whole genome sequencing kit according to one embodiment of the present invention includes at least one primer selected from the group consisting of SEQ ID NO:1 to SEQ ID NO:12, at least one primer selected from the group consisting of SEQ ID NO:13 to SEQ ID NO:20, primers of SEQ ID NO:21 and SEQ ID NO:22, and a custom sequencing primer of SEQ ID NO:23. The kit preferably further includes a primer selected from the group consisting of SEQ ID NO:24, SEQ ID NO:25, and SEQ ID NO:26. Specifically, one of the primers of SEQ ID NO:24, SEQ ID NO:25, and SEQ ID NO:26 is used with MiniSeq, NextSeq, HiSeq3000 / 4000 Illumina sequencing platforms to read index 2. Further description is provided below. Among the three primers, SEQ ID NO:24 is particularly preferred.

[0063] In an alternative embodiment, the low-pass whole genome sequencing kit includes at least one primer selected from the group consisting of SEQ ID NO:52 to SEQ ID NO:63, at least one primer selected from the group consisting of SEQ ID NO:13 to SEQ ID NO:20, primers of SEQ ID NO:21 and SEQ ID NO:22, a custom sequencing primer of SEQ ID NO:23, and a custom index 2 primer selected from the group consisting of SEQ ID NO:24, SEQ ID NO:25, and SEQ ID NO:26.

[0064] A low-pass whole genome sequencing kit according to another embodiment of the present invention comprises at least one primer selected from the group consisting of SEQ ID NO:27 to SEQ ID NO:38, at least one primer selected from the group consisting of SEQ ID NO:39 to SEQ ID NO:46, primers of SEQ ID NO:21 and SEQ ID NO:22, and a primer of SEQ ID NO:47. Preferably, the kit also comprises SEQ ID NO:48. Even more preferably, the kit comprises a primer of SEQ ID NO:49.

[0065] In an alternative embodiment, the low-pass whole genome sequencing kit comprises at least one primer selected from the group consisting of SEQ ID NO:64 to SEQ ID NO:75, at least one primer selected from the group consisting of SEQ ID NO:39 to SEQ ID NO:46, primers of SEQ ID NO:21 and SEQ ID NO:22, and primer of SEQ ID NO:47. Preferably, the kit also comprises SEQ ID NO:48. Even more preferably, the kit comprises primer of SEQ ID NO:49.

[0066] The method of genome-wide copy number profiling according to the present invention comprises: Sequencing the developed DNA library using one of the sequencing library preparation kits described above; Analyzing sequencing read depth across different regions of the genome; determining a copy number value for a region of the genome by comparing the number of reads in the region of the genome to the number of reads expected in the same region of a reference genome; Includes.

[0067] Protocol 1 for low-pass whole-genome sequencing on the Illumina platform Definitive Restriction Site Whole Genome Amplification (DRS-WGA) Single cell DNA was amplified using the Ampli1™ WGA kit (Menarini Silicon Biosystems) according to the manufacturer's instructions. 5 μL of WGA amplified DNA was diluted with 5 μL of nuclease-free water and purified using the SPRIselect beads (Beckman Coulter) system (1.8× ratio). DNA was eluted in 12.5 μL and quantified by dsDNA HS Assay in a Qubit™ 2.0 fluorometer.

[0068] P7 alone elongation A one-step PCR extension was performed using the Ampli1™ PCR kit (Menarini Silicon Biosystems) and LIB_IL_index D7xx (one of the primers SEQ ID NO:1-SEQ ID NO:12) in a volume of 15 μL. Each PCR reaction mix contained 1.5 μL Ampli1™ PCR reaction buffer (10×), 3 μL of one primer LIB_IL_index D7xx [2.5 μM] ranging from SEQ ID NO:1-SEQ ID NO:12, 0.51 μL Ampli1™ PCR dNTPs (10 mM), 0.37 μL BSA, 0.12 μL Ampli1™ PCR Taq polymerase, WGA-purified DNA (10 ng-75 ng) and Ampli1™ water to a final volume of 15 μL.

[0069] An Applied Biosystems™ 2720 thermal cycler was set as follows: 95° C. for 4 min, 1 cycle (95° C. for 30 sec, 62° C. for 30 sec, 72° C. for 3 min).

[0070] SPRIselect bead purification and P5 single extension 15 μL of Ampli1™ WGA amplicon from the previous step was purified using SPRIselect beads (Beckman Coulter) system (1.5× ratio). DNA was eluted in 15 μL of PCR reaction mix prepared as follows: 1.5 μL of Ampli1™ PCR reaction buffer (10×), 3 μL of one primer LIB_IL_index D5xx (one primer of SEQ ID NO:13 to SEQ ID NO:20) [2.5 μM], 0.51 μL of Ampli1™ PCR dNTPs (10 mM), 0.37 μL of BSA, 0.12 μL of Ampli1™ PCR Taq polymerase and 9.5 μL of Ampli1™ water. P5 single extension PCR reaction was performed in the presence of beads.

[0071] An Applied Biosystems™ 2720 thermal cycler was set as follows: 95° C. for 4 min, 1 cycle (95° C. for 30 sec, 60° C. for 30 sec, 72° C. for 3 min).

[0072] 2.5M NaCl PEG20% solution purification and library amplification 15 μL of Ampli1™ WGA amplicon from the previous step was purified using 2.5M NaCl PEG 20% solution (1.5× ratio). DNA was eluted in 15 μL of PCR reaction mix prepared as follows: 1.5 μL of Ampli1™ PCR reaction buffer (10×), 1 μL of one primer adapter P5 (SEQ ID NO:21) and 1 μL of one primer adapter P7 (SEQ ID NO:22) (7.5 μL each), 0.51 μL of Ampli1™ PCR dNTPs (10 mM), 0.37 μL of BSA, 0.12 μL of Ampli1™ PCR Taq polymerase and 10.5 μL of Ampli1™ water.

[0073] The library amplification PCR reactions were carried out in the presence of beads.

[0074] An Applied Biosystems™ 2720 thermal cycler was set as follows: 95° C. for 4 min, 1 cycle (95° C. for 30 sec, 60° C. for 30 sec, 72° C. for 2 min), 10 cycles (95° C. for 30 sec, 60° C. for 30 sec, 72° C. for 2 min (20 sec extension per cycle)) and a final extension at 72° C. for 7 min.

[0075] In a preferred embodiment, two additional cycles during library amplification can be added to increase the library concentration by at least 2 (~4) fold by extending the cycle to 12: 95°C for 30 sec, 60°C for 30 sec, 72°C for 2 min, increasing the total number of cycles from 11 to 13.

[0076] Final library purification The amplified libraries (containing Illumina sequencing adapter sequences) were finally purified using 2.5M NaCl PEG20% solution (1.5x ratio) and eluted in 15μL Ampli1™ water. The quality of the purified libraries was qualified by Agilent DNA 7500, DNA 1000 or DNA HS kits on a 2100 Bioanalyzer™ to obtain equimolar pools, and quantified by dsDNA HS Assay on a Qubit™ 2.0 fluorometer. Based on the average size of the library (typically an average length of 700 bp is experimentally observed in DRS-WGA products obtained using the Ampli1™ WGA kit (Menarini Silicon Biosystems)), the library concentration obtained in the quantification step can be converted to nM as known to those skilled in the art (e.g., 1 ng / μL=2.5 nM for an average library size of 600 bp, 1 ng / μL=2 nM for an average library size of 700 bp, and 1 ng / μL=1.9 nM for an average library size of 800 bp). See Illumina Technical Note: DNA Sequencing "Nextera™ Library Validation and Cluster Density Optimization" Pub No. 770-2013-003.

[0077] As a preferred alternative, the purified libraries were quantified by quantitative real-time PCR (qPCR), which accurately quantifies functional libraries, specifically fragments with the correct adapters at each end (heteroadapter fragments), based on a standard curve generated from control template dilutions.

[0078] Sequencing on the MiSeq Sequencing System 4 nM of the pool was denatured with 0.1 N NaOH for 5 min. The denatured sample was diluted with HT1 buffer (Illumina) to obtain a 20 pM denatured library. 600 μL of the denatured library was loaded onto a MiSeq reagent cartridge (Illumina).

[0079] Single-end 150 base reads or paired-end reads (75PE) were generated using Illumina MiSeq v3 chemistry.

[0080] The custom Read 1 sequencing primer (SEQ ID NO: 23) was then diluted in HT1 to give a final concentration of 0.5 μM. 600 μL of the diluted custom Read 1 sequencing primer was loaded into a MiSeq reagent cartridge (Illumina).

[0081] Sequencing on HiSeq 1000 / 1500 and 2000 / 2500 Systems The 4 nM pool was denatured with 0.1 N NaOH for 5 min.

[0082] Single-end 100 base reads or paired-end reads (100PE) were generated using Illumina HiSeq v2 chemistry in rapid run mode or Illumina HiSeq v4 chemistry in high output run mode.

[0083] Next, the custom Read 1 sequencing primer (SEQ ID NO:23) was diluted in HT1 to give a final concentration of 0.5 μM.

[0084] Sequencing on NextSeq, HiSeq 3000 and 4000, NovaSeq series and HiSeq X Ten systems The 4 nM pool was denatured with 0.1 N NaOH for 5 min.

[0085] Single-end 150 base reads or paired-end reads (100PE) were generated using specific chemistry on the Illumina platform.

[0086] Next, the custom Read 1 sequencing primer (SEQ ID NO: 23) and custom primer index 2A (i5) [LNA-5'] (SEQ ID NO: 24) were diluted in HT1 to give a final concentration of 0.5 μM.

[0087] The following table summarizes the primer sequences compatible with DRS-WGA for all Illumina platforms (5' to 3' sequences, 5' and 3' omitted):

[0088] [Table 1]

[0089] The following table summarizes the primer sequences compatible with MALBAC-WGA for all Illumina platforms (5' to 3' sequences, 5' and 3' omitted):

[0090] [Table 2]

[0091] As shown in Figures 2A and 2B, when using DRS-WGA, the LIB reverse complement is the target of the primers SEQ ID NO: 1 to SEQ ID NO: 20 listed in Table 1. Additionally, because the final library lacks the target sequence of the Illumina Read 1 sequencing primer, a custom Read 1 sequencing primer (SEQ ID NO: 23) was designed. The custom Read 1 sequencing primer (SEQ ID NO: 23) contains the LIB sequence and is complementary to the LIB reverse complement sequence.

[0092] Of note, in this sequencing setup, this method allows the construction of high complexity libraries based on our input Ampli1™ WGA products, thus avoiding the use of the PhiX spike-in control (Illumina).

[0093] Additionally, the sequencing run is preferably performed using the custom Read 1 sequencing primer (SEQ ID NO: 23), and since the PhiX DNA library lacks the target sequence for the custom Read 1 sequencing primer, the PhiX DNA is not sequenced.

[0094] Furthermore, the final library obtained by the method of the present invention does not have the regular Illumina sequencing adapters used on the MiniSeq, NextSeq, HiSeq 3000 and 4000 Illumina systems to read index 2 (i5).

[0095] For this reason, these platforms use custom primer index 2 (SEQ ID NO: 24 or SEQ ID NO: 25 or SEQ ID NO: 26) to enable accurate reading of index i5. It is noteworthy that custom sequencing primer index 2 contains LIB sequence. In particular, custom primer index 2A (i5) [LNA-5'] (SEQ ID NO: 24) and custom primer index 2A (i5) [LNA-3'] (SEQ ID NO: 25) have three LNA (locked nucleic acid [LNA™], Exiqon) modified nucleotides, indicated by a "+" next to them in Table 1 (e.g., "+A"). Furthermore, custom primer index 2 (i5) [RNA] (SEQ ID NO: 26) is formed of 15 RNA nucleotides, indicated by a "r" next to them (e.g., "rA").

[0096] Similar considerations as above also apply mutatis mutandis when using the MALBAC compatible primers listed in Table 2 (SEQ ID NO: 27 to SEQ ID NO: 49).

[0097] As a further alternative embodiment capable of producing libraries suitable for all Illumina platforms (and where the primary WGA library is a DRS-WGA library), the following primer combinations can be used: At least one primer selected from the group consisting of SEQ ID NO: 52 to SEQ ID NO: 63; At least one primer selected from the group consisting of SEQ ID NO: 13 to SEQ ID NO: 20; Primers of SEQ ID NO: 21 and SEQ ID NO: 22, A custom sequencing primer of SEQ ID NO:23 and a custom index 2 primer selected from the group consisting of SEQ ID NO:24, SEQ ID NO:25 and SEQ ID NO:26.

[0098] The primers of SEQ ID NO:52 to SEQ ID NO:63 contain the sequence of the 8 nucleotides required for the 5' of the first primer 3' section instead of the RDSP sequence to increase the annealing temperature of the custom sequencing primer. It is noteworthy that the second primers (SEQ ID NO:13 to SEQ ID NO:20) contain the same sequence of the 8 nucleotides 5' of the second primer 3' section. The custom sequencing primer of SEQ ID NO:23 is used on the Illumina sequencing platform to read Read 1 and / or Read 2, and the custom sequencing primer of SEQ ID NO:24 (or SEQ ID NO:25 or SEQ ID NO:26) is used on the Illumina sequencing platform to read Index 1 and Index 2.

[0099] The following table summarizes the primer sequences compatible with DRS-WGA on the Illumina platform according to this embodiment (sequences in the 5' to 3' direction, 5' and 3' omitted):

[0100] [Table 3]

[0101] As a further alternative embodiment capable of producing libraries suitable for all Illumina platforms (and where the primary WGA library is a MALBAC library), the following primer combinations can be used: At least one primer selected from the group consisting of SEQ ID NO: 64 to SEQ ID NO: 75; At least one primer selected from the group consisting of SEQ ID NO: 39 to SEQ ID NO: 46; Primers of SEQ ID NO: 21 and SEQ ID NO: 22, Custom read primer of sequence number 47.

[0102] It is preferred to also use the primer of SEQ ID NO: 48. It is even more preferred to also use the primer of SEQ ID NO: 49.

[0103] The following table summarizes the primer sequences compatible with DRS-WGA on the Illumina platform according to this embodiment (sequences in the 5' to 3' direction, 5' and 3' omitted):

[0104] [Table 4] EXAMPLES

[0105] Example 1 Sequenced reads were aligned to the hg19 human reference genome using the BWA MEM algorithm (Li H. and Durbin R., 2010).

[0106] Copy number calls were obtained without control samples using the Control-FREEC (Boeva ​​V. et al., 2011) algorithm. Read counts were corrected for GC content and mappability (uniqMatch option), and window sizes were determined by the software using a coefficient of variation = 0.06. Major ploidy parameters were set based on the ploidy of the genetic material tested, and no contamination correction was used.

[0107] Plots of CNA profiles of circulating tumor cells (CTCs) and white blood cells (WBCs), two single cells belonging to NCI-441 and SW-480 cell lines separated by DEPArray™ (Menarini Silicon Biosystems), and single blood cells separated by DEPArray™, were obtained using custom python scripts as shown in Figures 4A and 4B.

[0108] Of note from the figure is that autochromosome 22 shows large gains and losses that are manifested as absolute copy numbers in isolated single tumor cells.

[0109] Ploidy values ​​are represented on the y-axis to provide a better correspondence of profiles with the partitioned data (black line) and improve CNA calling: dots above the assessed predominant ploidy can be considered as gains and dots below the assessed predominant ploidy can be considered as losses.

[0110] On the other hand, there is no expected increase or decrease in normal WBC cells.

[0111] Example 2 Single tumor cells and 50 stromal cell populations from disaggregated FFPE sections were quantitatively isolated by DEPArray™ (Menarini Silicon Biosystems) and whole genomes were amplified using the Ampli1™ WGA kit. Figure 5 shows the results of low-pass whole genome sequencing performed by the method disclosed above. The figure shows gain and loss of copy number alteration (CNA) profiles only in single tumor cells. These high-quality CNA profiles confirm that the method is highly resistant to DNA degradation and demonstrates that it is a reliable and useful method for molecular characterization of tumor heterogeneity down to the single cell level in FFPE tissues.

[0112] Example 3 Figures 6A and 6B show the results of an inter-method and inter-platform (Ion Torrent and Illumina) comparison.

[0113] Specifically, Figure 6A shows the copy number alteration (CNA) profile of NCI-H23 single cells obtained by the low-pass whole genome sequencing method on Ion Torrent (as shown in US Patent No. 5,999,633) and the method according to the present invention on the Illumina platform. The results obtained on the two platforms are in good agreement with each other.

[0114] Furthermore, hierarchical clustering based on the CNA profiles of NCI-H441 and WBC single cells shows that samples are stratified by sample type (tumor and normal) but not by method or platform (Figure 6B).

[0115] In conclusion, the Ampli1™ low-pass method on both the Ion Torrent and Illumina sequencing platforms shows high concordance of CNA profiles.

[0116] Although the invention has been described only with respect to the methodology of Ampli1™ WGA, the techniques described clearly apply mutatis mutandis to any other type of WGA (e.g. MALBAC) involving libraries with self-complementary 5' and 3' regions, as will be apparent to the skilled artisan.

[0117] In an alternative embodiment of the invention, the first primer or the second primer (1PR and / or 2PR) may be modified to bind at its 5' end to a solid substrate in at least one step of the process.

[0118] As an example, the first primer (1PR) may contain a biotinylated 5' end. After the first PCR cycle (step b), streptavidin-coated beads are added to the reaction tube to capture the DNA library fragments extended with the first primer, while the unextended fragments not including the first primer are eluted with primers dNTPs and polymerase. In step c, a second primer (2PR) is provided (along with other PCR reagents) to hybridize with the DNA library fragments extended with the first primer, and polymerization occurs using the bead-bound DNA library fragments left in the tube as templates to generate heteroadapter DNA fragments. After washing the reaction mix of step c, the heteroadapter can be further amplified by PCR in the same tube using a third primer (3PR) and a fourth primer (4PR) (along with other PCR reagents) (step d). Alternatively, the heteroadapter can be denatured, eluted from the tube, and further amplified by PCR in another tube using a third primer (3PR) and a fourth primer (4PR) (step d).

[0119] As another example, the first primer (1PR) may be covalently bound to magnetic beads. After the first PCR cycle (step b), the DNA library fragments extended with the first primer are retained bound to the beads, while the unextended fragments containing the known 3' sequence section (3SS) but not the first primer are eluted. In step c, the second primer (2PR) hybridizes to the known 3' sequence section (3SS) of the DNA library fragments extended with the first primer, and polymerization occurs using the bead-bound DNA library fragments left in the tube as templates to produce heteroadapter DNA fragments. After washing the reaction mix of step c, the heteroadapter can be further amplified by PCR in the same tube using a third primer (3PR) and a fourth primer (4PR) (step d). Alternatively, the heteroadapter can be denatured, eluted from the tube, and further amplified by PCR in another tube using a third primer (3PR) and a fourth primer (4PR) (step d).

[0120] One advantage of these embodiments is that the size selective effect associated with the use of SPRI beads is not included in the process. This allows smaller and / or larger length fragments than normally retained by SPRI beads to be present in the final sequenceable library. One disadvantage is the increased complexity of kit reagents that is inherent to the binding of primers to beads or biotin. The presence of a wider range of fragments is particularly beneficial when seeking to achieve full genome sequencing at higher depths by low-pass WGS and complete resequencing across genome-wide copy number profiling.

[0121] As yet another example of the embodiment, after using the first primer (1PR) in the first PCR cycle (step b), SPRI purification is performed to remove the remaining first primer, and then the second primer (2PR) bound to magnetic beads (either directly by covalent bonding or indirectly by biotin modification of the second primer and biotin-streptavidin interaction with streptavidin-coated magnetic beads) is used for the above-mentioned one PCR cycle in step c, then the reaction mix and the remainder of the mono-adapter fragment are eluted, and only the resulting hetero-adapter fragment is retained in the tube by magnetic force, and then proceed to step d of PCR amplification with the third primer (3PR) and fourth primer (4PR).

[0122] advantage The massively parallel sequencing library generation method of the present invention allows rapid, efficient, reliable, and cost-effective generation of libraries containing only heterologous adapter fragments for sequencing platforms that require such libraries. In particular, the method of the present invention allows generation of sequence-ready libraries of double-stranded DNA in a streamlined, one-day, single-tube workflow.

[0123] Finally, the advantages of low-pass whole genome sequencing over array CGH (aCGH), the current leading technology for CNV investigation, should be highlighted. Array CGH (aCGH) is based on the use of differentially labeled test and reference genomic DNA samples that are simultaneously hybridized to DNA targets arrayed on glass slides or other solid platforms. However, the use of low quality / low amount of DNA has proven difficult and remains limited for certain applications. Furthermore, the assessment of chromosome copy number based on low-pass whole genome sequencing methods may bring several advantages compared to aCGH, including reduced costs of DNA sequencing, improved partial or segmental aneuploidy detection as a result of the possible increased resolution of chromosome analysis, and copy number variation calling without controls. In addition, the possibility of automating sequencing library preparation allows for the minimization of human error, reduction of hands-on time, and the realization of higher throughput and consistency.

Claims

**Claim 1** A method for generating a super-parallel sequencing library, comprising: a. preparing a primary WGA DNA library (pWGA lib) comprising fragments including a known 5' sequence section (5SS), a central sequence section (MSS), and a known 3' sequence section (3SS) that is reverse complementary to the known 5' sequence section, wherein the known 5' sequence section (5SS) includes a WGA library universal sequencing adapter, the central sequence section (MSS) includes an insert section (IS) corresponding to at least the DNA sequence of the unamplified original DNA before WGA, and the central sequence section (MSS) optionally further includes an adjacent 5' intermediate section (F5) and / or an adjacent 3' intermediate section (F3); b. performing one PCR cycle in the primary WGA DNA library using at least one first primer (1PR) including at least a first primer 5' section (1PR5S) and a first primer 3' section (1PR3S) to obtain a WGA DNA library extended by the first primer, wherein the first primer 5' section (1PR5S) includes at least one first sequencing adapter (1PR5SA), and the first primer 3' section (1PR3S) hybridizes to the known 3' sequence section (3SS); c. performing one PCR cycle in the WGA DNA library extended by the first primer using at least one second primer (2PR) including a second primer 5' section (2PR5S) and a second primer 3' section (2PR3S) to obtain a WGA DNA library extended by the first primer and the second primer, wherein the second primer 5' section (2PR5S) includes at least one second sequencing adapter (2PR5SA) different from the at least one first sequencing adapter (1PR5SA), and the second primer 3' section (2PR3S) hybridizes to the known 3' sequence section (3SS). d. Amplifying the WGA DNA library extended with the first primer and the second primer by PCR using at least one third primer (3PR) containing the first sequencing adapter (1PR5SA) and at least one fourth primer (4PR) containing the second sequencing adapter (2PR5SA) to obtain an amplified WGA DNA library extended with the first primer and the second primer; A method comprising the above. **Claim 2** The method according to claim 1, wherein the at least one first primer (1PR) further comprises at least one read sequencing primer sequence (1PRSEQ) at the 3' position of the first primer 5' section (1PR5S) and at the 5' position of the first primer 3' section (1PR3S). **Claim 3** The first primer 5' section (1PR5S) further comprises at least one first sequencing barcode (1PR5BC) at the 3' position of the at least one first sequencing adapter (1PR5SA) and at the 5' position of the first primer 3' section (1PR3S), and / or the second primer 5' section (2PR5S) further comprises at least one second sequencing barcode (2PR5BC) at the 3' position of the at least one second sequencing adapter (2PR5SA) and at the 5' position of the second primer 3' section (2PR3S). The method according to claim 1 or 2. **Claim 4** The method according to any one of claims 1 to 3, further comprising a step of purifying the WGA DNA library extended with the first primer after step b, and / or a step of purifying the WGA DNA library extended with the first primer and the second primer after step c, and / or a step of purifying the amplified WGA DNA library extended with the first primer and the second primer after step d. **Claim 5** The method according to any one of claims 1 to 4, wherein the WGA library universal sequencing adapter is a DRS-WGA library universal sequencing adapter or a MALBAC library universal sequencing adapter. **Claim 6** The method according to claim 5, wherein the WGA library universal sequence adapter is a DRS-WGA library universal sequence adapter.

7. The method according to claim 5 or 6, wherein the DRS-WGA library universal sequence adapter is SEQ ID NO: 50 and the MALBAC library universal sequence adapter is SEQ ID NO:

51.

8. A method for low-pass whole genome sequencing, comprising: preparing a plurality of barcoded ultra-parallel sequencing libraries by the method according to any one of claims 3 to 7, and pooling samples obtained using different sequencing barcodes (BC); sequencing the pooled libraries; A method comprising:

9. The step of pooling samples using different sequencing barcodes (BC) further comprises quantifying DNA in each of the barcoded ultra-parallel sequencing libraries and normalizing the amount of the barcoded ultra-parallel sequencing libraries. The method for low-pass whole genome sequencing according to claim 8.

10. At least one first primer (1PR) comprising at least a first primer 5' section (1PR5S) and a first primer 3' section (1PR3S), wherein the first primer 5' section (1PR5S) comprises at least one first sequencing adapter (1PR5SA), and the first primer 3' section (1PR3S) hybridizes to a known 3' sequence section (3SS) that is reverse complementary to a known 5' sequence section (5SS) containing the WGA library universal sequence adapter of the fragment of the primary WGA DNA library (pWGA lib), and the fragment further comprises a central sequence section (MSS) at the 3' position of the known 5' sequence section (5SS) and the 5' position of the known 3' sequence section (3SS). At least one first primer, ​ ​ At least one second primer (2PR) comprising at least a second primer 5' section (2PR5S) and a second 3' section (2PR3S), wherein the second primer 5' section (2PR5S) comprises at least one second sequencing adapter (2PR5SA) different from the at least one first sequencing adapter (1PR5SA), and the second 3' section hybridizes to the known 3' sequence section (3SS) of the fragment, and at least one second primer, At least one third primer (3PR) comprising the first sequencing adapter (1PR5SA), At least one fourth primer (4PR) comprising the second sequencing adapter (2PR5SA), A massively parallel sequencing library preparation kit comprising.

11. a) One or more primers selected from the group consisting of SEQ ID NO: 1 to SEQ ID NO: 12, one or more primers selected from the group consisting of SEQ ID NO: 13 to SEQ ID NO: 20, and the primers of SEQ ID NO: 21 and SEQ ID NO: 22, or, b) One or more primers selected from the group consisting of SEQ ID NO: 27 to SEQ ID NO: 38, one or more primers selected from the group consisting of SEQ ID NO: 39 to SEQ ID NO: 46, and the primers of SEQ ID NO: 21 and SEQ ID NO: 22, or, c) One or more primers selected from the group consisting of SEQ ID NO: 52 to SEQ ID NO: 63, one or more primers selected from the group consisting of SEQ ID NO: 13 to SEQ ID NO: 20, and the primers of SEQ ID NO: 21 and SEQ ID NO: 22, or, d) One or more primers selected from the group consisting of SEQ ID NO: 64 to SEQ ID NO: 75, one or more primers selected from the group consisting of SEQ ID NO: 39 to SEQ ID NO: 46, and the primers of SEQ ID NO: 21 and SEQ ID NO: 22, A massively parallel sequencing library preparation kit comprising.

12. a) At least one primer selected from the group consisting of SEQ ID NO: 1 to SEQ ID NO: 12, at least one primer selected from the group consisting of SEQ ID NO: 13 to SEQ ID NO: 20, the primers of SEQ ID NO: 21 and SEQ ID NO: 22, and the custom sequencing primer of SEQ ID NO: 23, or, b) at least one kind of primer selected from the group consisting of SEQ ID NO: 27 to SEQ ID NO: 38, at least one kind of primer selected from the group consisting of SEQ ID NO: 39 to SEQ ID NO: 46, the primers of SEQ ID NO: 21 and SEQ ID NO: 22, and the primer of SEQ ID NO: 47, or, c) at least one kind of primer selected from the group consisting of SEQ ID NO: 52 to SEQ ID NO: 63, at least one kind of primer selected from the group consisting of SEQ ID NO: 13 to SEQ ID NO: 20, the primers of SEQ ID NO: 21 and SEQ ID NO: 22, the custom sequencing primer of SEQ ID NO: 23, and the custom index 2 primer selected from the group consisting of SEQ ID NO: 24, SEQ ID NO: 25 and SEQ ID NO: 26, or, d) at least one kind of primer selected from the group consisting of SEQ ID NO: 64 to SEQ ID NO: 75, at least one kind of primer selected from the group consisting of SEQ ID NO: 39 to SEQ ID NO: 46, the primers of SEQ ID NO: 21 and SEQ ID NO: 22, and the primer of SEQ ID NO: 47, A low-pass whole genome sequencing kit comprising.

13. The low-pass whole genome sequencing kit according to claim 12, wherein option a) further comprises a primer selected from the group consisting of SEQ ID NO: 24, SEQ ID NO: 25 and SEQ ID NO: 26, and option b) or option d) further comprises the primer of SEQ ID NO:

48.

14. The low-pass whole genome sequencing kit according to claim 13, wherein option b) or option d) further comprises the primer of SEQ ID NO:

49.

15. A method for genome-wide copy number profiling, comprising: Preparing a DNA library by the method according to any one of claims 1 to 7 using the sequencing library preparation kit according to claim 10 or 11; Sequencing the prepared DNA library; Analyzing the sequencing read depth across different regions of the genome; Determining the copy number value of the region of the genome by comparing the number of reads in the genomic region with the expected number of reads in the same region of the reference genome; A method comprising.

Images