The sequence listing provided in the file entitled Sequence_listing_PCTCN2020105117.txt, which is an ASCII text file that was created on Jan. 25, 2022, and which comprises 84,565 bytes, is hereby incorporated by reference in its entirety.
The present invention relates to the technical field of molecular biology, and particularly to a library preparation method and application.
Generally, the current sequencing and analysis of a sequence in a target region first requires a library preparation, and methods of library preparation are nothing more than capture library preparation and amplification library preparation.
The capture library preparation is an enrichment library preparation targeted at a relatively large region of a genome, such as a tens or hundreds of gene whole exon regions, while the multiplex amplification library preparation is to perform a target capture and sequencing and analysis on specific hotspot regions, or the whole exon regions of individual genes.
The method for amplification library preparation is to design corresponding specific primers according to a target region. These primers are then used to conduct a multiplex amplification on target sequences. It should be noted that these specific primers will directly carry sequencing adapters or bridging sequences, and then sequencing adapters are added thereto by a secondary PCR amplification, which is the process of a normal amplification library preparation. There are some problems in the application of the existing amplification library preparation methods. For example, the library preparation process is relatively cumbersome, requiring at least two cycles of PCR amplification and two corresponding library purifications, calling for numerous manual operation time that impose high requirements on operators, thus not conducive to popularization. Moreover, primer design and system optimization are relatively complicated; the cost of library preparation is high; and the entire library preparation process is time-consuming.
Aiming at various problems existing in the amplification library preparation, the present invention provides the following technical solutions.
One purpose of the present invention is to provide a primer combination for preparing an amplicon library for detecting the variation of a target gene.
The primer combination provided by the present invention includes:
a forward outer primer F1, a forward inner primer F2, and a reverse primer R that are designed according to a target amplicon;
The forward outer primer F1 is sequentially composed of a sequencing adapter sequence 1, a barcode sequence for distinguishing different samples, and a universal sequence;
The forward inner primer F2 is sequentially composed of a universal sequence and a forward specific primer sequence of the target amplicon (a molecular tag is not required when detecting a tissue sample);
The reverse outer primer R is sequentially composed of a sequencing adapter 2 and a reverse specific primer sequence of the target amplicon.
In the above primer combination, optionally, a molecular tag is required when detecting low frequency mutations, and the forward inner primer F2 is sequentially composed of a universal sequence, a molecular tag sequence, and a forward specific primer sequence of the target amplicon.
The molecular tag sequence is composed of 6-30 bases, consisting random bases and 0-N(N is an integer ≥0) set(s) of specific bases; the specific bases are set in the random bases, for example, 1 set, 2 sets, 3 sets, or 4 sets; the specific bases in each set are composed of 1-5 bases, such as 1 base, 2 bases, 3 bases, 4 bases, or 5 bases.
The base sequence of each set is randomly selected, and the molecular tag sequence is used to distinguish different starting DNA template molecules. In a library preparation process, except for the fixed position and constant composition of the specific bases in the molecular tag sequence, the types of bases (A, T, C) of the random bases can be selected at will.
For example, in an embodiment of the present invention, the specific bases are set as 1 or 2 sets, with the sequence of ACT and/or TGA; for example, in the present embodiment, the molecular tag sequence is NNNNNACTNNNNTGA (SEQ ID NO: 13), where ACT and TGA are the specific bases, N is a random base of A, T, C, or G.
In the above primer combination, the sequencing adapter 1 and the sequencing adapter 2 are corresponding sequencing adapters selected according to different sequencing platforms.
In the above primer combination, the sequencing platform is an Illumina platform, the sequencing adapter 1 is 15, and the sequencing adapter 2 is 17;
or the sequencing platform is an Ion Torrent platform, the sequencing adapter 1 is A, and the sequencing adapter 2 is P;
or the sequencing platform is a BGI/MGI platform;
or, the nucleotide sequence of the universal sequence is shown in SEQ ID NO: 1.
Another purpose of the present invention is to provide a kit for preparing an amplicon library for detecting the variation of a target gene.
The kit provided by the present invention includes the above-mentioned primer combination.
The above kit further includes a polymerase chain reaction (PCR) amplification buffer and a DNA polymerase system.
Another purpose of the present invention is to provide any one of the following applications of the primer combination or the kit described above:
(1) an application in preparing the amplicon library for detecting the variation of the target gene;
(2) an application in detecting mutation sites or variations in a target region of a sample to be tested;
(3) an application in detecting a variation frequency of the target region of the sample to be tested.
Another purpose of the present invention is to provide a method of preparing an amplicon library for detecting a variation of a target gene.
The method provided by the present invention includes the following steps:
taking the DNA or cDNA of a sample to be tested as a template, carrying out a one-step PCR amplification using the above primer combination or the above kit to obtain an amplified product, i.e., the amplicon library of the target gene.
In the above method, the molar ratio of the forward outer primer F1, the forward inner primer F2, and the reverse primer R in an amplification system for the one-step PCR amplification is (5-20):(1-20):(5-20).
In the above method, the sample to be tested is a tissue sample, a frozen sample, a puncture sample, a formalin-fixed paraffin-embedded (FFPE) sample, blood, urine, cerebrospinal fluid, pleural fluid, or other body fluids.
The application of the above method in detecting mutation sites or variations of the target gene of the sample to be tested.
The application of the above method in detecting a variation frequency of the target gene of the sample to be tested.
The amplicon library prepared by the above method also falls within the protection scope of the present invention.
Another purpose of the present invention is to provide a method for detecting the variation of the target gene of the sample to be tested.
The method provided by the present invention includes the following steps:
1) preparing an amplicon library of the target gene by the above method;
2) evenly mixing the amplicon libraries of the target genes of all samples, and then diluting to obtain a sequencing DNA library;
3) sequencing the sequencing DNA library to obtain a sequencing result, and analyzing the variation of the target gene of the sample to be tested according to the sequencing result.
Another purpose of the present invention is to provide a method of detecting a variation frequency in a target region of a sample to be tested.
The method provided by the present invention includes the following steps:
1) preparing an amplicon library of the target gene by the above method;
2) evenly mixing the amplicon libraries of the target genes of all samples, and then diluting to obtain a sequencing DNA library;
3) sequencing the sequencing DNA library to obtain a sequencing result, and calculating the variation frequency of the target gene of the sample to be tested according to the sequencing result.
Variation frequency=number of mutation clusters/total number of effective clusters×100%.
In the above method, the sample to be tested is an in vitro tissue sample, a frozen sample, a puncture sample, an FFPE sample, blood, urine, cerebrospinal fluid, or pleural fluid.
In the above method, optionally,
the nucleotide sequence of the universal sequence is shown in SEQ ID NO: 1;
the nucleotide sequence of the sequencing adapter 1 is shown in SEQ ID NO: 2;
the nucleotide sequence of the sequencing adapter 2 is shown in SEQ ID NO: 17.
For example, when the target gene to be tested is EGFR, optionally, the corresponding forward specific primer sequence and reverse specific primer sequence are respectively shown in SEQ ID NO: 14 and SEQ ID NO: 18, or, SEQ ID NO: 15 and SEQ ID NO: 19, or, SEQ ID NO: 21 and SEQ ID NO: 24, or, SEQ ID NO: 22 and SEQ ID NO: 25;
When the target gene to be tested is ERBB2, optionally, the corresponding forward specific primer sequence and reverse specific primer sequence are respectively shown in SEQ ID NO: 16 and SEQ ID NO: 20, or, SEQ ID NO: 23 and SEQ ID NO: 26;
When the target gene to be tested is EML4, optionally, the corresponding forward specific primer sequence and reverse specific primer sequence are respectively shown in SEQ ID NO: 27 and SEQ ID NO: 31, or, SEQ ID NO: 28 and SEQ ID NO: 31;
When the target gene to be tested is LMNA, optionally, the corresponding forward specific primer sequence and reverse specific primer sequence are respectively shown in SEQ ID NO: 29 and SEQ ID NO: 32;
When the target gene to be tested is MYC, optionally, the corresponding forward specific primer sequence and reverse specific primer sequence are respectively shown in SEQ ID NO: 30 and SEQ ID NO: 33.
For example, the barcode sequences are all nucleotides with a length of 6-12 nt, no more than 3 consecutive bases, and a GC content of 40-60%;
The universal sequence 1 and the universal sequence 2 generally have a length of 16-25 nt, and a GC content of 35-65%, without consecutive bases or obvious secondary structure;
For example, the molecular tag sequence is a sequence containing 6-15 random bases; including but not limited to the above sequences; in the embodiment of the present invention, for example, the barcode sequences for distinguishing different samples are shown in SEQ ID NO: 3 to SEQ ID NO: 12;
The variation can be point mutation, deletion or insertion, or fragment fusion.
The experimental methods used in the following embodiments, unless otherwise specified, are all conventional methods.
The materials, reagents, etc. used in the following embodiments, unless otherwise specified, are commercially available.
I. Design of Primers for One-Step Amplicon Sequencing Library Preparation
The present invention provides an amplification library preparation method to prepare a second-generation sequencing library, and the structures of primers involved in the method are as follows (see
forward outer primer F1: 5′-sequencing adapter sequence 1+Barcode sequence+universal sequence-3′;
forward inner primer F2: 5′-universal sequence+molecular tag sequence+gene forward specific primer sequence-3′;
or forward inner primer F2: 5′-universal sequence+gene forward specific primer sequence-3′ (molecular tag is required when detecting low frequency mutations, and molecular tag is not required when detecting tissue samples);
reverse primer R: 5′-sequencing adapter sequence 2+gene reverse specific primer sequence-3′.
When detecting low-frequency mutations, the structure of the forward inner primer F2 is: 5′-universal sequence+molecular tag sequence+gene forward specific primer sequence-3′.
Among them, the barcode sequence is a nucleic acid sequence that is used to distinguish different samples; a sample to be tested corresponds to a barcode sequence. The barcode sequence is 6-12 nt in length, and has no more than 3 consecutive bases, and a GC content of 40-60%, and the primer where the Barcode sequence is introduced has no obvious secondary structure, etc.
The forward outer primer F1 is used to distinguish different samples. The same sample has the same forward outer primer F1 regardless of detection sites.
The molecular tag sequence is used to mark different starting DNA template molecules (templates of different amplicons), and a starting DNA template molecule corresponds to a molecular tag sequence.
The molecular tag sequence includes random bases and at least one set of specific bases, the specific bases are set in the random bases, for example, 1 set or 2 sets; each set of specific bases is composed of 1-5 bases, for example, 3 bases or 4 bases. In a library preparation process, except for the fixed position and constant composition of the specific bases in the molecular tag sequence, the types of bases (A, T, G, C) of the random bases are randomly selected.
The starting templates of the sequencing results are classified using the molecular tag sequences, which can eliminate amplification errors and sequencing errors. In the present embodiment, two types of specific bases are used: ACT and TGA, which can be used separately or in combination.
Gene forward specific primer sequence and gene reverse specific primer sequence are primer sequences (respectively including the required forward primers and corresponding reverse primers to amplify different target regions) used to amplify specific target regions;
The universal sequence 1 is a specific nucleic acid sequence, which can be changed according to actual needs. The universal sequence 1 has a length of 16-25 nt, and a GC content of 35-65%, without consecutive bases or obvious secondary structure.
In present embodiment, the universal sequence used is GGCACCCGAGAATTCCA (SEQ ID NO: 1), with a length of 17 nt;
The sequencing adapter sequence 1 and the sequencing adapter sequence 2 are specific sequences that need to be introduced to primers during sequencing, and can specifically correspond to Ion Torrent, Illumina, or BGISEQ/MGISEQ sequencing platforms.
If the sequencing platform is the Illumina platform, the sequencing adapter sequences 1 and 2 are I5 and I7, respectively, and the adapter sequences are complementary to the primer sequences on the chip. The adapter is introduced to link a nucleic acid fragment to a vector.
If the sequencing platform is the Ion Torrent platform, the sequencing adapter sequences 1 and 2 are A and P, respectively, the A adapter is used for sequencing and complementary to the sequencing primer, and the P adapter is complementary to the sequence on the vector, so as to link a template to the vector.
If the sequencing platform is the BIISEQ/MGISEQ platform, the sequencing adapters are required for sequencing, which are specific sequences meeting the requirements of single-strand circularization, subsequent DNB preparation, and sequencing.
When the second-generation sequencing library is used in a sequencing, multiple samples will be tested simultaneously. As such, a set of forward outer primers F1 will be designed. M forward outer primers F1 correspond to M samples, and the barcode sequence in each forward outer primer F1 is different; P forward inner primers F2 and corresponding Q (generally P=Q, but there are also situations where P does not equal Q, for example, in a detection of RNA fusion genes) reverse primers R are designed according to the number P of amplicons required for the target capture region on each sample, and the structures of the molecular tags in the P forward inner primers F2 are identical.
II. Amplification Principle of One-Step Amplicon Sequencing Library
The primers design of one-step rapid amplification library preparation technology are as described above. When amplifying a template DNA/RNA, the procedure shown in
III. Establishment of Detection Method
1. One-Step Amplification
The primers synthesized in section I above were prepared as follows:
The forward outer primer F1 was dissolved in water to a primer concentration of 100 μM, and the forward inner primers F2 were respectively dissolved in water to a primer concentration of 100 μM. Subsequently, the various primers were mixed at an equimolar ratio to form the forward outer primer MIX1. The reverse primers R were respectively dissolved in water to 100 μM, and then mixed at an equimolar ratio into the reverse primer MIX2.
The genomic DNA of multiple samples to be tested was extracted.
The reagents shown in Table 1 were successively added to a 0.2 ml eight-row tube or 96-well plate (each type of nucleic acid sample was extracted according to the instruction of the specific manufacturer's kit provided in the embodiment):
The procedure of the above PCR amplification is shown in Table 2.
After the PCR reaction was completed, the PCR product obtained was the amplicon library.
2. Magnetic Bead Purification and Qubit Quantification
After the PCR reaction was completed, the Agencourt AMPure XP Kit (Cat. No. A63880/A63881/A63882) from Beckman Coulter Inc. was used for purification. The operation steps were as follows:
1) The Agilent court AMPure XP Kit was taken out 30 min in advance, fully vortex and put aside at room temperature.
2) After the PCR reaction, the magnetic beads were fully vortexed again, 24 μl of magnetic beads were added to the system, blow repeatedly more than 5 times or vortex fully, and put aside at room temperature for 5 min.
3) The Eppendorf (EP) tubes were transferred to a magnetic stand and put aside for 5 min until the solution was clear, by using a pipette to carefully remove the supernatant without contacting the magnetic beads.
4) 100 μl of freshly prepared 80% ethanol solution was added to each tube, the EP tubes were slowly rotated for 2 cycles on the magnetic stand, followed by putting aside for 5 min and discarding the supernatant.
5) The step 4 was repeated one more time.
6) The EP tubes were opened and put aside at room temperature to allow a complete liquid volatilization until surfaces of the magnetic beads became matte. The magnetic beads should not be dried excessively.
7) The EP tubes were removed from the magnetic stand, 30 μl of PCR-grade purified water was added, followed by vortex to mix well, and putting aside at room temperature for 10 min.
8) The EP tubes in the previous step were placed on the magnetic stand for 2 min or until the solution was clear, followed by using a pipette to carefully absorb the supernatant on the side away from the magnet without contacting the magnetic beads.
A purified amplicon library was obtained.
The purified amplicon library was subjected to a DNA library concentration determination and an Agilent 2200 TapeStation Systems detection using Qubit 2.0.
3. Sequencing and Result Analysis
The purified amplicon libraries of multiple samples were mixed at an equal concentration, and then diluted to 100 PM to obtain a DNA library for amplicon sequencing. Sequencing was performed (sequenator used was Ion GeneStudio™ S5 Plus System, Thermofisher, A38195), after data processing and analysis (S5 Torrent Server), the mutations and mutation frequency of a tested sample were obtained.
The calculation method of the variation frequency of the library with molecular tags was as follows:
Since the original template was subjected to molecular marking during the library amplification process, the calculation method of the mutation frequency was as follows:
In the sequencing results, DNA molecules with the same kind of molecular tags were defined as a cluster, and DNA molecules with the same kind of molecular tags were amplified products of an initial DNA template, that is, a series of DNA molecules obtained by amplification using the same original template;
Whether mutations occurred in the cluster or not was confirmed. If the proportion of a specific type of bases in a certain position in the cluster was greater than or equal to 80%, the cluster was recorded as an effective cluster. If the number of mutant DNA molecules with molecular tags in the effective cluster accounted for greater than or equal to 80%, it was recorded as a mutation cluster;
Variation frequency=number of mutation clusters/total number of effective clusters×100%.
Notes: It is statistically significant only when the number of DNA molecules in the same cluster (a sequence sequenced) in the sequencing results is ≥2.
1. Design of Primers for One-Step Amplicon Sequencing Library Preparation
The detection region of this experiment contained three amplicons (EGFR L858R, 19del and insertion mutations of ERBB2);
The test samples included two frozen lung cancer tissue samples (sample 1, sample 2), four lung cancer FFPE (formalin fixed paraffin-embedded tissue samples) samples (sample 3, sample 4, sample 5, sample 6), and two white blood cell samples from healthy subjects (sample 7, sample 8). The mutations of the above eight samples were already known.
The primers (eight Barcode sequences were used in the present embodiment) shown in Table 3 were designed according to the three amplicons (EGFR L858R, 19del and insertion mutations of ERBB2):
The sequencing adapter is suitable for the Ion GeneStudio™ S5 Plus System sequencing platform.
II. One-Step Amplicon Sequencing Library
Nucleic acid extraction and purification kit (DNA extraction from FFPE samples: GeneRead DNA FFPE kit, Qiagen, 180134; DNA extraction from frozen tissue samples: QIAamp DNA Mini Kit 250, QIAGEN, 51306).
1. One-Step Amplification
The PCR product was obtained according to step 1 in section III of Embodiment 1.
The amplification system is shown in Table 4.
2. Magnetic Bead Purification and Qubit Quantification
Same steps were performed as those in step 2 of section III of Embodiment 1.
The PCR product was purified and recovered by the magnetic bead (Agencourt AMPure XP, Beckman Coulter, A63880), and DNA library concentration determination and Agilent 2200 TapeStation Systems detection were conducted using Qubit 2.0.
The result of the Agilent 2200 TapeStation Systems detection is shown in
3. Sequencing and Result Analysis
The PCR products of all samples were mixed at an equal concentration and diluted to 100 pM to obtain a DNA library for sequencing.
The sequencing results are shown in Table 6:
EGFR: p.E746_A750delELREA indicates a deletion of the 746th-750th amino acids ELREA (E: Glu glutamic acid; L: Leu leucine; R: Arg arginine; E: Glu glutamate; A: Ala alanine) of the EGFR gene, which is a kind of EGFR 19del;
EGFR: p.K745_E749delKELRE indicates a deletion of the 745th-749th amino acids KELRE (K: Lys lysine; E: Glu glutamic acid; L: Leu leucine; R: Arg arginine; E: Glu glutamic acid) of the EGFR gene, which is a kind of EGFR 19del;
ERBB2: p.A775_G776insYVMA indicates an insertion of YVMA (Y: Tyr Tyrosine; V: Val Valine; M: Met Methionine; A: Ala alanine) between the 775th alanine (A) and the 776th glycine (G) of the ERBB2 gene, corresponding to ERBB2 in Table 3.
The 63 gene detection product is a product of tumor liquid biopsy of Genetron Health (Beijing) Co., Ltd. It targets all solid tumor patients and applies high-throughput and high-precision second-generation sequencing technology to comprehensively detect mutations of 63 gene loci closely related to tumor-targeted therapy and occurrence and development (including mutation analysis of 58 genes, rearrangement analysis of 10 genes, and CNV detection of 7 genes), covering the target region with a sequencing depth of 20,000×, and reaching a detection sensitivity of 0.1%, which provides comprehensive and high-value reference information for precise medication, molecular typing, and curative effect and recurrence monitoring.
The above results show that the library prepared by the method of the present invention, when used for sequencing, leads to the variation information of tested tissue samples including point mutations, deletion mutations and insertion mutations consistent with that obtained by the known 63 gene detection.
The samples in this experiment were plasma samples from lung cancer patients, including plasma samples from four different patients and two healthy subjects (the variations of the samples were already known), cfDNA was extracted using the kit (MagMAX™ Cell-Free DNA Isolation Kit, Applied Biosystems™, A29319), and the library was prepared using a primer pool with molecular tags containing EGFR L858R, 19del and insertion mutations of ERBB2.
I. Design of Primers for One-Step Amplicon Sequencing Library Preparation
The primers (forward outer primers were identical, others were different, and six barcode sequences were used in the present embodiment) shown in Table 7 were designed according to three amplicons (EGR L858R, 19del and insertion mutations of ERBB2):
II. One-Step Amplicon Sequencing Library
Nucleic acid extraction and purification kit (DNA extraction from FFPE samples: GeneRead DNA FFPE kit, Qiagen, 180134; DNA extraction from frozen tissue samples: QIAamp DNA Mini Kit 250, QIAGEN, 51306).
1. One-Step Amplification
The PCR product was obtained according to step 1 in section III of Embodiment 1.
2. Magnetic Bead Purification and Qubit Quantification
Same steps were performed as those in step 2 of section III of Embodiment 1.
The PCR product was purified and recovered by the magnetic bead (Agencourt AMPure XP, Beckman Coulter, A63880), and detected by Qubit 2.0 and Agilent 2200 TapeStation Systems.
The result of the Agilent 2200 TapeStation Systems is shown in
3. Sequencing and Result Analysis
The PCR products of all samples were mixed at an equal concentration and diluted to 100 μM to obtain a DNA library for amplicon sequencing.
The sequencing results are shown in Table 10:
EGFR: p.E746_A750delELREA indicates a deletion of the 746th-750th amino acids ELREA (E: Glu glutamic acid; L: Leu leucine; R: Arg arginine; E: Glu glutamate; A: Ala alanine) of the EGFR gene, which is a kind of EGFR 19del;
ERBB2: p.A775_G776insYVMA indicates an insertion of YVMA (Y: Tyr Tyrosine; V: Val Valine; M: Met Methionine; A: Ala alanine) between the 775th alanine (A) and the 776th glycine (G) of the ERBB2 gene, corresponding to ERBB2 in Table 7.
The library prepared by the method of the present invention, when used for sequencing, leads to variation information of tested plasma cfDNA samples including point mutations, deletion mutations and insertion mutations consistent with that obtained by the known 63 gene detection. The amount of ctDNA extracted from the patient 1 sample is large. After the patient 1 sample is diluted by 5 times, the detection of L858R with a frequency of 4.6‰ is still obtained (after deduplicating the data Reads: mutation cluster=2; total cluster at the locus=4380).
The samples in this experiment were fine-needle aspiration (FNA) puncture samples of 3 thyroid cancer patients with gene fusion (gene fusion information was already known) and FNA puncture samples of 2 patients with benign thyroid nodules. RNA samples were extracted using MagMAX™ FFPE DNA/RNA Ultra Kit (Applied Biosystems™, A31881) according to the manufacturer's instruction, and then reverse transcription was conducted using SuperScript™ VILO™ MasterMix (Invitrogen™, 11755050) according to the manufacturer's kit instruction.
I. Design of Primers for One-Step Amplicon Sequencing Library Preparation
The primers (forward outer primers were identical to those in Table 3, others were different, and five barcode sequences were used in the present embodiment) shown in Table 11 were designed according to gene fusion: the primers for detecting gene fusion were designed before and after the breakpoint, and there was no fixed forward and reverse primer matching; the forward and reverse primers designed for the fusion breakpoint were shown as below, ALK_20 and ELM4_6/EML4_13 were combined separately to detect two ALK-EML4 fusion forms.
II. One-Step Amplicon Sequencing Library
1. One-Step Amplification
The PCR product was obtained according to step 1 in section III of Embodiment 1.
2. Magnetic Bead Purification and Qubit Quantification
Same steps were performed as those in step 2 of section III of Embodiment 1.
The PCR product was purified and recovered by the magnetic bead (Agencourt AMPure XP, Beckman Coulter, A63880), and detected by Qubit 2.0 and Agilent 2200 TapeStation Systems.
The result of the Agilent 2200 TapeStation Systems detection is shown in
3. Sequencing and Result Analysis
The PCR products of all samples were mixed at an equal concentration and diluted to 100 μM to obtain a DNA library for amplicon sequencing.
The sequencing results are shown in Table 14:
EML4-ALK-V3a (E6a A20) corresponds to EML4_6 and ALK_20 in Table 11;
EML4-ALK-V1 (E13 A20) corresponds to EML4_13 and ALK_20 in Table 11.
The 63 gene detection product used Agilent's customized probes to perform capture library preparation. The product has been used for detecting thousands of clinical plasma samples, and the performance of the product is stable.
The library prepared by the method of the present invention, when used for sequencing, leads to fusion mutation forms of tested samples consistent with the mutation information of samples obtained by the known 63 gene detection.
The foregoing embodiments are only used to illustrate the present invention. The structure, connection mode, and manufacturing process of each component can be changed. Any equivalent transformation and improvement based on the technical solution of the present invention should not be excluded from the protection scope of the present invention.
I. Design of Primers for One-Step Amplicon Sequencing Library Preparation
The forward outer primers were the same as those in Table 3, others were different, and the barcodes were determined according to the number of samples in the library preparation.
The forward inner primer F2: universal sequence+forward specific primer sequence
The reverse primer R: sequencing adapter 2+reverse specific primer sequence
II. One-Step Amplicon Sequencing Library
1. One-Step Amplification
The PCR product was obtained using 0.5 pg of gDNA of white blood cell in plasma from healthy subject as a starting sample according to step 1 in section III of Embodiment 1.
2. Magnetic Bead Purification and Qubit Quantification
Same steps were performed as those in step 2 of section III of Embodiment 1.
The PCR product was purified and recovered by the magnetic bead (Agencourt AMPure XP, Beckman Coulter, A63880), and detected by Qubit 2.0 and Agilent 2200 TapeStation Systems.
The result of the Agilent 2200 TapeStation Systems detection is shown in
3. Sequencing and Result Analysis
The PCR products of all samples were mixed at an equal concentration and diluted to 100 μM to obtain a DNA library for amplicon sequencing.
The sequencing results are shown in
Comparison of three primers used in the method of the comparative example and the method of the present invention and four primers in the prior art
I. Design of Primers for One-Step Amplicon Sequencing Library Preparation
The structures of the 3 primers and the 4 primers designed are shown in
The Present Invention:
The 3 primers of the present invention were designed according to the design principle of Embodiment 1:
The forward outer primers were the same as those in Table 3, specifically, there were 67 barcode sequences; the universal sequences were the same as those in Table 15, and the forward specific gene sequences were P1_B2_F1 to P1_B2_F67 in Table 15;
The sequencing adapter 2 was the same as that in Table 15, and the reverse specific primer sequences were P1_B2_R1 to P1_B2_R67 in Table 15;
Control:
The 4 primers were designed according to the following principles in the prior art:
Design Principles:
Barcode primer F1: sequencing adapter 1+barcode sequence+universal sequence 1;
Forward inner primer F2: universal sequence 1+molecular tag+specific base sequence+forward specific primer sequence;
Reverse outer primer R1: sequencing adapter 2+universal sequence 2;
Reverse inner primer R2: universal sequence 2+reverse specific primer sequence;
The sequencing adapter 1+barcode sequence described above are shown in Table 3.
The rest sequences are shown in Table 16 below:
II. One-Step Amplicon Sequencing Library
The method was the same as that in step 2 of Embodiment 2.
The sequencing results are analyzed as follows:
1. The Homogeneity Results of the Libraries Prepared by Triple-Functional Component Primer Pool and Quadruple-Functional Component Primer Pool
The homogeneity of the amplicons library is a very important indicator of the quality of the library. Good homogeneity of the library indicates a higher coverage of the target region of the library, and a better detection accuracy of the panel covering region. For this purpose, under the premise of ensuring the intact functional structure of the primer, the primer design of the amplicon is improved. The improved primer structure is optimized and simplified from the original F1+F2+R1+R2 (quadruple-functional primer components) to F1+F2+R (triple-functional primer components). This design will increase the stability of the reaction system and ensure the homogeneity of amplicons in the library.
Amplifications were respectively carried out on the primer set of the present invention and the control primer set using the same white blood cell DNA sample as a template.
The results are shown in
2. 30 ng of cfDNA was Used in a Library Preparation by One-Step Primer Pool, and the Number of Molecular Tag Types/the Number of Clusters of One of the Amplicons was Obtained after Data Analysis
Amplifications were respectively carried out on the primer set of the present invention and the control primer set using the same cfDNA sample as a template.
The results are shown in
3. Background Noise at the Level of 0.1‰-1‰ of the Libraries Prepared by Two Methods and Subjected to Sequencing (Same as 2)
Amplifications were respectively carried out on the primer set of the present invention and the control primer set using the same cfDNA sample as a template.
The results are shown in
Good amplification homogeneity, high capture efficiency of original template molecules, high-fidelity DNA polymerase, ultra-low background noise, and the introduction of molecular tags eventually facilitate the triple-functional primer component to achieve an effective detection of ultra-low frequency mutation at the level of 3‰. The primer structure of this library preparation method has been fully optimized, and the performance of this library preparation method is much better than that of the traditional library preparation method for low-frequency mutation detection.
The comparison results of the one-step rapid library preparation method of the present invention, the ordinary amplification library preparation method and the capture library preparation method are shown in Table 17 and Table 18.
After the amplicon library is prepared, there may be amplification products of target fragments, primer dimers or multimers, and fragment products of non-specific amplification in the system. A high proportion of the amplification products of target fragments becomes an extremely important indicator for evaluating the quality of the amplicon library. Table 18 shows the present method has great advantages in terms of the proportion of target fragments of the library as compared to the control method.
In order to solve the current difficulties in library preparation, the present invention has developed the one-step rapid amplification library preparation method. Compared with the traditional capture method, the amplification library preparation method has the following advantages (
The present invention has the following merits because of adopting the above technical solutions:
1. Little sample consumption and high utilization rate. The capture efficiency of the original template molecules in the sample is high, and thus a relatively low amount of starting templates is required. When performing germline mutation detection, even just a pg-level amount of starting templates is required. When performing low frequency mutation detection of cfDNA, a limited amount of starting templates can achieve a higher template capture efficiency, thereby achieving an effective capture of trace ctDNA molecules, realizing a lower detection limit and a higher sensitivity;
2. Ultra-low detection limit. The unique primer design, supporting PCR reaction system, reaction conditions, and subsequent information analysis and noise reduction system ultimately result in the lowest mutation detection limit of 3‰, making it possible to realize an accurate detection of ultra-early stage and trace amounts of ctDNA sample mutations;
3. Good homogeneity of library. The innovative primer structure design and supporting reaction system result in the optimal homogeneity of amplicons in the library. When conducting a multiplex amplification, the different structural characteristics of the sequences of various amplicons and the different amplification efficiencies of various primers will eventually result in a huge difference in the abundance of amplicons in the library. How to balance the difference in the abundance of amplicons is a key indicator to evaluate the quality of the library. The components of the triple-functional primer used in the present method have obvious advantages over the components of the quadruple-functional primer. Specifically, the cooperation of primer composition and reaction system ensures that the method can control differential amplifications of amplicons at a reduced number of cycles, and then a method like universal primer amplification is used. Since there is no competition between R1 and R2 in the following figure, a stable low differential amplification is achieved in subsequent cycles;
4. High repeatability. The components of the quadruple-functional primer will increase the uncertainty of the reaction system and reaction conditions, and are more sensitive to sample quality, reaction system and external environmental influences. While the components of the triple-functional primer have been improved in this aspect, and the simpler components result in a better system stability, and a higher repeatability and accuracy of sample detection;
5. Easy operation and time saving. The traditional capture library preparation technology has cumbersome operations and long procedures. The entire library preparation process takes nearly 48 h and imposes high requirements on operators. The ordinary amplification library preparation method requires at least two cycles of PCR and two cycles of purification, including subsequent QPCR quantification. The entire library preparation process requires at least one working day. The present invention only involves one-step PCR reaction and corresponding product purification steps, and the entire library preparation process can be completed within 1.5 h, thereby simplifying the library preparation operation process and saving time of the library preparation (the library preparation can be completed within 1.5 h, and the entire process from the library preparation to the completion of sequencing and to the completion of the bioinformatic analysis can be controlled within 22 h);
6. Able to detect multiple gene mutation types. Starting from a DNA sample, SNP, SNV, Ins/Del, methylation, gene or exon level copy number variation, and chromosome arm level copy number variation can be detected. In addition, after adding molecular tags to primers, mutations at the level of as low as 1‰ can be further detected. Starting with a RNA sample, the expression of specific genes, the fusion of specific genes, etc. can be detected;
7. Multiple sample types. The starting sample can be fresh tissue samples, frozen samples, puncture samples, FFPE samples and other tissue sample types. Meanwhile, isolated cfDNA or CTC in blood, urine, cerebrospinal fluid, and pleural fluid can also be detected. After DNA or RNA is extracted from normal samples, library preparation can be conducted by one-step rapid amplification library preparation method;
8. Effective elimination of cross-contamination between samples. The barcode sequences that distinguish different samples are added at the beginning of PCR, and the simplification of the operation process and steps effectively eliminates possible cross-contaminations during the library preparation process, especially when detecting low frequency mutations, cross-contamination between samples is extremely prone to determining as a false positive mutation;
9. Reduced cost of library preparation. Compared with the traditional capture technology, the cost required for library preparation using the present method is greatly reduced. The capture probes used in the traditional capture library preparation are expensive, and the reagents and consumables involved in the lengthy experimental process also increase the cost of capture library preparation. In contrast, the one-step library preparation process requires a greatly reduced amount of reagents and consumables, and the cost of library preparation is much lower than that of the traditional capture library preparation method. At the same time, compared with the one-step rapid amplification library preparation method, at least one cycle of additional PCR and purification and the QPCR quantitation of the library in the normal amplification library preparation method will also greatly increase the cost of library preparation. Compared with the components of the prior quadruple-functional primer, the components of the triple-functional primer lead to low consumption of total primer and each component of primer, thus having a lower cost advantage;
10. Space saving. Since this method requires only one cycle of PCR, the laboratory requires only 3 rooms (sample extraction, PCR amplification room, library purification and sequencing), which saves space as compared to the conventional library preparation where 4 rooms (sample extraction, PCR1, PCR2, and library purification and sequencing) are required.
Flexible and simple library preparation method, allowing detection of multiple mutation types, and extremely high detection sensitivity are the biggest features of the present invention.
Number | Date | Country | Kind |
---|---|---|---|
201910694844.3 | Jul 2019 | CN | national |
The present application is a U.S. National Phase of International Application Number PCT/CN2020/105117 filed Jul. 28, 2020, and claims priority to Chinese Application Number 201910694844.3 filed Jul. 30, 2019.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2020/105117 | 7/28/2020 | WO |