Patent Application
20180142293
The advent of massively parallel DNA sequencing has ushered in a new era of genomic exploration by making simultaneous genotyping of hundreds of billions of base-pairs possible at small fraction of the time and cost of traditional Sanger methods [1]. Because these technologies digitally tabulate the sequence of many individual DNA fragments, unlike conventional techniques which simply report the average genotype of an aggregate collection of molecules, they offer the unique ability to detect minor variants within heterogeneous mixtures [2].
This concept of “deep sequencing” has been implemented in a variety fields including metagenomics [3, 4], paleogenomics [5], forensics [6], and human genetics [7, 8] to disentangle subpopulations in complex biological samples. Clinical applications, such prenatal screening for fetal aneuploidy [9, 10], early detection of cancer [11] and monitoring its response to therapy [12, 13] with nucleic acid-based serum biomarkers, are rapidly being developed. Exceptional diversity within microbial [14, 15] viral [16-18] and tumor cell populations [19, 20] has been characterized through next-generation sequencing, and many low-frequency, drug-resistant variants of therapeutic importance have been so identified [12, 21, 22]. Previously unappreciated intra-organismal mosasism in both the nuclear [23] and mitochondrial [24, 25] genome has been revealed by these technologies, and such somatic heterogeneity, along with that arising within the adaptive immune system [13], may be an important factor in phenotypic variability of disease.
Deep sequencing, however, has limitations. Although, in theory, DNA subpopulations of any size should be detectable when deep sequencing a sufficient number of molecules, a practical limit of detection is imposed by errors introduced during sample preparation and sequencing. PCR amplification of heterogeneous mixtures can result in population skewing due to stoichastic and non-stoichastic amplification biases and lead to over- or under-representation of particular variants [26]. Polymerase mistakes during pre-amplification generate point mutations resulting from base mis-incorporations and rearrangements due to template switching [26, 27]. Combined with the additional errors that arise during cluster amplification, cycle sequencing and image analysis, approximately 1% of bases are incorrectly identified, depending on the specific platform and sequence context [2, 28]. This background level of artifactual heterogeneity establishes a limit below which the presence of true rare variants is obscured [29].
A variety of improvements at the level of biochemistry [30-32] and data processing [19, 21, 28, 32, 33] have been developed to improve sequencing accuracy. The ability to resolve subpopulations below 0.1%, however, has remained elusive. Although several groups have attempted to increase sensitivity of sequencing, several limitations remain. For example techniques whereby DNA fragments to be sequenced are each uniquely tagged [34, 35] prior to amplification [36-41] have been reported. Because all amplicons derived from a particular starting molecule will bear its specific tag, any variation in the sequence or copy number of identically tagged sequencing reads can be discounted as technical error. This approach has been used to improve counting accuracy of DNA [38, 39, 41] and RNA templates [37, 38, 40] and to correct base errors arising during PCR or sequencing [36, 37, 39]. Kinde et. al. reported a reduction in error frequency of approximately 20-fold with a tagging method that is based on labeling single-stranded DNA fragments with a primer containing a 14 bp degenerate sequence. This allowed for an observed mutation frequency of ˜0.001% mutations/bp in normal human genomic DNA [36]. Nevertheless, a number of highly sensitive genetic assays have indicated that the true mutation frequency in normal cells is likely to be far lower, with estimates of per-nucleotide mutation frequencies generally ranging from 10−9 to 10−11 [42]. Thus, the mutations seen in normal human genomic DNA by Kinde et al. are likely the result of significant technical artifacts.
Traditionally, next-generation sequencing platforms rely upon generation of sequence data from a single strand of DNA. As a consequence, artifactual mutations introduced during the initial rounds of PCR amplification are undetectable as errors—even with tagging techniques—if the base change is propagated to all subsequent PCR duplicates. Several types of DNA damage are highly mutagenic and may lead to this scenario. Spontaneous DNA damage arising from normal metabolic processes results in thousands of damaging events per cell per day [43]. In addition to damage from oxidative cellular processes, further DNA damage is generated ex vivo during tissue processing and DNA extraction [44]. These damage events can result in frequent copying errors by DNA polymerases: for example a common DNA lesion arising from oxidative damage, 8-oxo-guanine, has the propensity to incorrectly pair with adenine during complementary strand extension with an overall efficiency greater than that of correct pairing with cytosine, and thus can contribute a large frequency of artifactual G→T mutations [45]. Likewise, deamination of cytosine to form uracil is a particularly common event which leads to the inappropriate insertion of adenine during PCR, thus producing artifactual C→T mutations with a frequency approaching 100% [46].
It would be desirable to develop an approach for tag-based error correction, which reduces or eliminates artifactual mutations arising from DNA damage, PCR errors, and sequencing errors; allows rare variants in heterogeneous populations to be detected with unprecedented sensitivity; and which capitalizes on the redundant information stored in complexed double-stranded DNA.
In one embodiment, a single molecule identifier (SMI) adaptor molecule for use in sequencing a double-stranded target nucleic acid molecule is provided. Said SMI adaptor molecule includes a single molecule identifier (SMI) sequence which comprises a degenerate or semi-degenerate DNA sequence; and an SMI ligation adaptor that allows the SMI adaptor molecule to be ligated to the double-stranded target nucleic acid sequence. The SMI sequence may be single-stranded or double-stranded. In some embodiments, the double-stranded target nucleic acid molecule is a double-stranded DNA or RNA molecule.
In another embodiment, a method of obtaining the sequence of a double-stranded target nucleic acid is provided (also known as Duplex Consensus Sequencing or DCS) is provided. Such a method may include steps of ligating a double-stranded target nucleic acid molecule to at least one SMI adaptor molecule to form a double-stranded SMI-target nucleic acid complex; amplifying the double-stranded SMI-target nucleic acid complex, resulting in a set of amplified SMI-target nucleic acid products; and sequencing the amplified SMI-target nucleic acid products.
In some embodiments, the method may additionally include generating an error-corrected double-stranded consensus sequence by (i) grouping the sequenced SMI-target nucleic acid products into families of paired target nucleic acid strands based on a common set of SMI sequences; and (ii) removing paired target nucleic acid strands having one or more nucleotide positions where the paired target nucleic acid strands are non-complementary (or alternatively removing individual nucleotide positions in cases where the sequence at the nucleotide position under consideration disagrees among the two strands). In further embodiments, the method confirms the presence of a true mutation by (i) identifying a mutation present in the paired target nucleic acid strands having one or more nucleotide positions that disagree; (ii) comparing the mutation present in the paired target nucleic acid strands to the error corrected double-stranded consensus sequence; and (iii) confirming the presence of a true mutation when the mutation is present on both of the target nucleic acid strands and appears in all members of a paired target nucleic acid family.
Single molecule identifier adaptors and methods for their use are provided herein. According to the embodiments described herein, a single molecule identifier (SMI) adaptor molecule is provided. Said SMI adaptor molecule is double stranded, and may include a single molecule identifier (SMI) sequence, and an SMI ligation adaptor (
The SMI adaptor molecule may form a “Y-shape” or a “hairpin shape.” In some embodiments, the SMI adaptor molecule is a “Y-shaped” adaptor, which allows both strands to be independently amplified by a PCR method prior to sequencing because both the top and bottom strands have binding sites for PCR primers FC1 and FC2 as shown in the examples below. A schematic of a Y-shaped SMI adaptor molecule is also shown in
In other embodiments, the SMI adaptor molecule is a “hairpin” shaped (or “U-shaped”) adaptor. A hairpin DNA product can be used for error correction, as this product contains both of the two DNA strands. Such an approach allows for reduction of a given sequencing error rate N to a lower rate of N*N*(1/3), as independent sequencing errors would need to occur on both strands, and the same error among all three possible base substitutions would need to occur on both strands. For example, the error rate of 1/100 in the case of Illumina sequencing [32] would be reduced to (1/100)*(1/100)*(1/3)=1/30,000.
An additional, more remarkable reduction in errors can be obtained by inclusion of a single-stranded SMI in either the hairpin adaptor or the “Y-shaped” adaptor will also function to label both of the two DNA strands. Amplification of hairpin-shaped DNA may be difficult as the polymerase must synthesize through a product containing significant regions of self-complementarity, however, amplification of hairpin-shaped structures has already been established in the technique of hairpin PCR, as described below. Amplification using hairpin PCR is further described in detail in U.S. Pat. No. 7,452,699, the subject matter of which is hereby incorporated by reference as if fully set forth herein.
According to the embodiments described herein, the SMI sequence (or “tag”) may be a double-stranded, complementary SMI sequence or a single-stranded SMI sequence. In some embodiments, the SMI adaptor molecule includes an SMI sequence (or “tag”) of nucleotides that is degenerate or semi-degenerate. In some embodiments, the degenerate or semi-degenerate SMI sequence may be a random degenerate sequence. A double-stranded SMI sequence includes a first degenerate or semi-degenerate nucleotide n-mer sequence and a second n-mer sequence that is complementary to the first degenerate or semi-degenerate nucleotide n-mer sequence, while a single-stranded SMI sequence includes a first degenerate or semi-degenerate nucleotide n-mer sequence. The first and/or second degenerate or semi-degenerate nucleotide n-mer sequences may be any suitable length to produce a sufficiently large number of unique tags to label a set of sheared DNA fragments from a segment of DNA. Each n-mer sequence may be between approximately 3 to 20 nucleotides in length. Therefore, each n-mer sequence may be approximately 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20 nucleotides in length. In one embodiment, the SMI sequence is a random degenerate nucleotide n-mer sequence which is 12 nucleotides in length. A 12 nucleotide SMI n-mer sequence that is ligated to each end of a target nucleic acid molecule, as described in the Example below, results in generation of up to 424 (i.e., 2.8×1014) distinct tag sequences.
In some embodiments, the SMI tag nucleotide sequence may be completely random and degenerate, wherein each sequence position may be any nucleotide. (i.e., each position, represented by “X,” is not limited, and may be an adenine (A), cytosine (C), guanine (G), thymine (T), or uracil (U)) or any other natural or non-natural DNA or RNA nucleotide or nucleotide-like substance or analog with base-pairing properties (e.g., xanthosine, inosine, hypoxanthine, xanthine, 7-methylguanine, 7-methylguanosine, 5,6-dihydrouracil, 5-methylcytosine, dihydouridine, isocytosine, isoguanine, deoxynucleosides, nucleosides, peptide nucleic acids, locked nucleic acids, glycol nucleic acids and threose nucleic acids). The term “nucleotide” as described herein, refers to any and all nucleotide or any suitable natural or non-natural DNA or RNA nucleotide or nucleotide-like substance or analog with base pairing properties as described above. In other embodiments, the sequences need not contain all possible bases at each position. The degenerate or semi-degenerate n-mer sequences may be generated by a polymerase-mediated method described in the Example below, or may be generated by preparing and annealing a library of individual oligonucleotides of known sequence. Alternatively, any degenerate or semi-degenerate n-mer sequences may be a randomly or non-randomly fragmented double stranded DNA molecule from any alternative source that differs from the target DNA source. In some embodiments, the alternative source is a genome or plasmid derived from bacteria, an organism other than that of the target DNA, or a combination of such alternative organisms or sources. The random or non-random fragmented DNA may be introduced into SMI adaptors to serve as variable tags. This may be accomplished through enzymatic ligation or any other method known in the art.
In some embodiments, the SMI adaptor molecules are ligated to both ends of a target nucleic acid molecule, and then this complex is used according to the methods described below. In certain embodiments, it is not necessary to include n-mers on both adapter ends, however, it is more convenient because it means that one does not have to use two different types of adaptors and then select for ligated fragments that have one of each type rather than two of one type. The ability to determine which strand is which is still possible in the situation wherein only one of the two adaptors has a double-stranded SMI sequence.
In some embodiments, the SMI adaptor molecule may optionally include a double-stranded fixed reference sequence downstream of the n-mer sequences to help make ligation more uniform and help computationally filter out errors due to ligation problems with improperly synthesized adaptors. Each strand of the double-stranded fixed reference sequence may be 4 or 5 nucleotides in length sequence, however, the fixed reference sequence may be any suitable length including, but not limited to 3, 4, 5 or 6 nucleotides in length.
The SMI ligation adaptor may be any suitable ligation adaptor that is complementary to a ligation adaptor added to a double-stranded target nucleic acid sequence including, but not limited to a T-overhang, an A-overhang, a CG overhang, a blunt end, or any other ligatable sequence. In some embodiments, the SMI ligation adaptor may be made using a method for A-tailing or T-tailing with polymerase extension; creating an overhang with a different enzyme; using a restriction enzyme to create a single or multiple nucleotide overhang, or any other method known in the art.
According to the embodiments described herein, the SMI adaptor molecule may include at least two PCR primer or “flow cell” binding sites: a forward PCR primer binding site (or a “flow cell 1” (FC1) binding site); and a reverse PCR primer binding site (or a “flow cell 2” (FC2) binding site). The SMI adaptor molecule may also include at least two sequencing primer binding sites, each corresponding to a sequencing read. Alternatively, the sequencing primer binding sites may be added in a separate step by inclusion of the necessary sequences as tails to the PCR primers, or by ligation of the needed sequences. Therefore, if a double-stranded target nucleic acid molecule has an SMI adaptor molecule ligated to each end, each sequenced strand will have two reads—a forward and a reverse read.
Adaptor 1 (shown below) is a Y-shaped SMI adaptor as described above (the SMI sequence is shown as X's in the top strand (a 4-mer), with the complementary bottom strand sequence shown as Y's):
Adaptor 2 (shown below) is a hairpin (or “U-shaped”) linker:
Following ligation of both adaptors to a double-stranded target nucleic acid, the following is structure is obtained:
When melted, the product will be of the following form (where “linker” is the sequence of adaptor 2):
This product is then PCR amplified. The reads will yield:
The sequences of the two duplex strands seen in the two sequence reads may then be compared, and sequence information and mutations will be scored only if the sequence at a given position matches in both of the reads.
This approach does not strictly require the use of an SMI tag, as the sheared ends can be used as identifiers to differentiate unique individual molecules from PCR duplicates. Thus the same concept would apply if one used any standard sequencing adaptor as “Adaptor 1” and the U-shaped linker as “Adaptor 2.” However described below, there are a limited number of shear points flanking any given genomic position and thus the power to sequence deeply is increased via inclusion of the SMI tag. A hybrid method using a combination of sheared ends and a shorter n-mer tag (such as 1 or 2 or 3 or 4 or more degenerate or semi-degenerate bases) in the adaptor may also serve as unique molecular identifiers. Another design may include use of any sequencing adaptor (such as one lacking an n-mer tag) in conjunction with an n-mer tag that is incorporated into the U-shaped linker molecule. Such a design would be of the following form (where X and Y represent complementary degenerate or semi-degenerate nucleotides):
Synthesis of such a design may be obtained in a number of ways, for example synthesizing a set of hairpin oligonucleotides in which each individual oligonucleotide encodes a complementary n-mer sequence, or alternatively by using a DNA polymerase to carry out extension from the following product (where X's represent degenerate nucleotides):
Inclusion of the SMI tag is also extremely useful for identifying correct ligation products, as the assay uses two distinct adaptors. This will yield multiple possible ligation products:
Product I.
Adaptor 1---------DNA---------Adaptor 2, which yields the desired product:
Product II.
Adaptor 1---------DNA---------Adaptor 1. This will result in the DNA being amplified as two separate strands, i.e. as occurs in the DCS approach described elsewhere in this document (the second copy of Adaptor 1 is shown below with the SMI as AAA-BBB to emphasize that every DCS adaptor has a distinct SMI sequence)
Product III.
Adaptor 2---------DNA---------Adaptor 2. This will result in a non-amplifiable circular product shown below:
Product III is non-amplifiable, given the absence of primer binding sites and thus will not be present in the final DNA sequences. Thus only Product II needs to be avoided. The formation of Product II can be minimized in the ligation step by using an excess of Adaptor 2 (relative to Adaptor 1). Then primarily Products I and III will be obtained, with minimal formation of Product II. Additionally, a variety of biochemical means of enriching for products containing adaptor 2 are possible such as using affinity probes that are complementary to the hairpin loop sequence itself. Product I results in the same SMI sequence in both the Read 1 and Read 2 sequence reads. In the example depicted above, Product I sequences can thus be identified by virtue of having matching SMIs of the form XXXX in Read 1 and XXXX in Read 2.
By contrast, in the case of Product II, the SMI sequences on either end of the sequenced molecule will arise from distinct DCS adaptors having different SMI sequences. In the example shown above, Product II sequences yield SMIs of the form XXXX (Read 1)-BBB (Read 2) upon sequencing of the top strand, and BBBB (Read 1)-XXXX (Read 2) upon sequencing of the bottom strand. Thus Product II sequences can be easily identified and computationally removed from the final sequence data.
Data resulting from Product II is useful, because Product II corresponds to the product analyzed under the approach detailed in the Example below. Product I contains a self--complementary hairpin sequence that can impair polymerase extension during amplification, however, this type of amplification has already been enabled in the technique of “Hairpin PCR” [50] which involves linking of the two strands followed by amplification with gene-specific primers. Amplification conditions that are compatible with amplification of hairpin DNA are thus already established. Moreover, ligation and amplification with circularizing “linkers” (i.e. hairpin linkers affixed to both ends of a fragment) has been demonstrated as a step in the Pacific Biosciences sample preparation workflow [49]. As the sequence of the linker itself does not matter in the workflow, the published linker sequences from either of these references would be adequate for use in the assay.
In some aspects of some embodiments, deliberate ligation of “U-shaped” adaptors or hairpin linkers containing 1) a double-stranded n-mer (or other form of degenerate or semi-degenerate double-stranded tag as enumerated above) plus 2) primer binding sites to both ends of a captured fragment may be desireable. Producing closed circles of captured material may help facilitate removal of non-captured DNA by exonuclease digestion given that circularized DNA will be protected from digestion by such enzymes. Additionally, closed circles may be pre-amplified using rolling circle amplification or serve as the substrate for continuous loop sequencing [49]. Recognition sites for restriction endonuclease digestion could be engineered into these adaptors to render closed loops open once again if more convenient for subsequent steps.
In another embodiment, flow cell sequences or PCR binding sites, again denoted as FC1 and FC2, may be included in both the PCR primers and the hairpin linker adaptor, as well as a ligatable sequence on the end of the hairpin linker (denoted as L below). The hairpin linker adaptor may additionally include one or more cleavable sequences, denoted as R in the example below (the R may be any appropriate restriction enzyme target sequence, or any other cleavable sequence). Such a hairpin linker design is shown below:
The target DNA with ligation site denoted as L is as follows:
Following ligation of the linker, the product may be amplified with PCR primers as follows:
The resultant product will be of the form:
After amplification of the product, the cleavage sites R may be cleaved to result in the following sequencable products:
and
These products may then be sequenced directly. This design has the advantage of allowing for targeted sequencing of a specific region of the genome, and furthermore avoids the need to sequence a hairpin product, as sequencing of a hairpin will be less efficient due to the self-complementarity present within the hairpin molecule.
In one embodiment, a single-stranded SMI sequence is incorporated into the single-stranded portion of the hairpin loop (regions of sequence complementarity are denoted as “=”). The SMI sequence is shown as four nucleotides in length in the following examples, but in practice an Nmer of any length, including approximately 3 to 20 nucleotides, will suffice.
Ligation of the hairpin linker and a Y-shaped sequencing adaptor (with PCR primer binding sites labeled as FC1 and FC2) yields the following product:
Melting and PCR amplification of this product yields the following DNA product:
FC1-------DNA-----NNNN------hairpin sequence------DNA′-------FC2
Following PCR duplication of the product and formation of consensus reads based upon the shared SMI sequence among all the PCR duplicates, the sequences of the two strands (denoted DNA and DNA′) can then be compared to form a duplex consensus sequence.
In another embodiment, a single-stranded SMI is incorporated into a modified “Y-shaped” sequencing adaptor in which PCR primer binding sites are located at the sites labeled FC1 and FC2 (regions of sequence complementarity are depicted as “=”)
It will be apparent to one skilled in the art that a single-stranded SMI sequence tag can be located in any of several positions within either the sequencing adaptor or the hairpin linker. The single-stranded SMI sequence tag can be synthesized as a random oligonucleotide sequence, or can be sequenced as a set of fixed sequences by synthesis on an array, or by any other suitable method known in the art.
SMI adaptors molecules containing a double-stranded, complementary, degenerate or semi-degenerate SMI tag can be made by any of a number of methods, including copying of a single-stranded SMI sequence by a DNA polymerase as described above or synthesis and annealing of two oligonucleotides containing complementary SMI sequences. An additional method involves synthesizing a set of linear oligonucleotides which will self-anneal into the appropriate form. Inclusion of a cleavable linker in each oligonucleotide will then allow for conversion of a “hairpin shaped” SMI adaptor molecule into a “Y-shaped” SMI adaptor molecule. For example, an oligonucleotide may be prepared of the following form:
In this schematic, X and Y represent complementary nucleotides, and U indicates a cleavable linker, such as uracil (which can be cleaved by combined treatment with uracil DNA glycosylase and apurinic endonuclease), although any other cleavable linker will suffice. The oligonucleotide may be designed with appropriate regions of self-complementarity to anneal into the following form:
The linker (e.g. uracil) may then be cleaved, yielding a DCS adaptor:
A double-stranded SMI hairpin linker can be constructed by an analogous method but without the need for a cleavable linker. For example, a set of nucleotides of known sequence where X and Y represent the complementary SMI sequences can be synthesized on an array, or by any other suitable method known in the art:
This oligonucleotide can then self-anneal to form a hairpin linker with complementary SMI sequences.
Any of the oligonucleotides described above can also include any ligatable sequence as overhangs on either the 5′ or 3′ end, or can be used for blunt end ligation.
DCS SMI adaptor molecules contain ligatable ends to allow attachment of the adaptor to a target DNA molecule. In some embodiments, the ligatable end may be complementary to a DNA overhang on the target DNA, for example, one generated by digestion of target DNA with a restriction endonuclease. Selective ligation of the adaptor to the targeted DNA containing the matching Single-stranded overhanging DNA sequence will then allow for partial purification of the targeted DNA. A non-limiting example of this embodiment is shown above in paragraphs [0048]-[0053]. In some embodiments, the DCS SMI adaptor molecule, or a hairpin linker SMI adaptor molecule, may additionally contain modifications such as biotin to facilitate affinity purification of target DNA that has ligated to the adaptor.
In another embodiments, specific PCR primers can selectively amplify specific regions of genome when the adaptor that is ligated to the other end of the molecule is a hairpin (or “U-shape”). Alternatively, this method may be used with or without the need for this cleavable hairpin sequence.
Another embodiment involves fragmentation of DNA at defined regions, for example by treatment of DNA with a site-specific restriction endonuclease or a mixture of such endonucleases, followed by annealing of a hairpin oligonucleotide linker, and amplification of the hairpin complex with PCR primers sufficient for amplification of the desired DNA sequence. Annealing of the hairpin linker to only one of the two ends of the DNA duplex could be accomplished by using different restriction enzymes to cut on either end of the target duplex, and then having the hairpin linker ligation adaptor being ligatable to only one of the two resultant ligatable ends.
The example shown below indicates forward and reverse PCR primers (labeled 1 and 2) in conjunction with a hairpin linker to allow linked amplification of both complementary strands of duplex DNA. Such amplification, in conjunction with a single-stranded or double-stranded SMI sequence, would allow for targeted amplification and high accuracy deep sequencing of a specific sequence of interest. In the schematic shown below, a single-stranded SMI sequence is incorporated into PCR primer FC1. It would be apparent to one skilled in the art that the SMI sequence could also be incorporated in primer FC2, or in the hairpin linker.
Amplified Product:
FC1NNNN DNA----hairpin sequence----DNA′FC2
This product can then be subjected to consensus sequencing analysis. The SMI sequence allows one to group together products of PCR amplification arising from a single molecule of duplex DNA. The sequences of the two DNA strands can then be compared for error correction.
The SMI adaptor molecules described herein have several uses. In some embodiments, the SMI adaptor molecules described herein may be used in methods to obtain the sequence or other sequence-related information of a double-stranded target nucleic acid molecule. According to the embodiments described herein, the term “double-stranded target nucleic acid molecule” includes a double-stranded DNA molecule or a double-stranded RNA molecule. Thus, the SMI adaptor molecules and methods of use described herein are applicable to genotyping and other applications related to sequencing of DNA molecules, but are also applicable to RNA sequencing applications such as for sequencing of double-stranded RNA viruses. Methods for sequencing RNA may include any of the embodiments described herein with respect to DNA sequencing, and vice-versa. For example, any double stranded target nucleic acid molecule may be ligated to an SMI adaptor molecule which includes a double-stranded RNA or DNA n-mer tag and an RNA or DNA ligation adapter as described above. Methods exist for directly sequencing RNA [51]; alternatively, the ligated product may be reverse transcribed into DNA, and then sequenced as a double-stranded target DNA molecule.
In one embodiment, the double-stranded target nucleic acid molecule may be a sheared double-stranded DNA or RNA fragment. The sheared target DNA or RNA molecule may be end repaired and a double-stranded target nucleic acid sequence ligation adaptor may be added to each end of the sheared target DNA or RNA molecule. The double-stranded target nucleic acid sequence ligation adaptor may be any suitable ligation adaptor that is complementary to the SMI ligation adaptor described above including, but not limited to a T-overhang, an A-overhang, a CG overhang, blunt end or any other ligatable sequence. In some embodiments, the double-stranded target nucleic acid sequence ligation adaptor may be made using a method for A-tailing or T-tailing with polymerase extension; adding an overhang with a different enzyme; using a restriction enzyme to create a ligatable overhang; or any other method known in the art.
Methods to obtain the sequence or other sequence-related information of a double-stranded target nucleic acid molecule may include a step of ligating the double-stranded target nucleic acid molecule to at least one SMI adaptor molecule, such as those described above, to form a double-stranded target nucleic acid complex. In one embodiment, each end of the double-stranded target nucleic acid molecule is ligated to an SMI adaptor molecule. The double-stranded target nucleic acid complex is then amplified by a method known in the art (e.g., a PCR or non-PCR method known in the art), resulting in a set of uniquely labeled, amplified SMI-target nucleic acid products. These products are then sequenced using any suitable method known in the art including, but not limited to, the Illumina sequencing platform, ABI SOliD sequencing platform, Pacific Biosciences sequencing platform, 454 Life Sciences sequencing platform, Ion Torrent sequencing platform, Helicos sequencing platform, and nanopore sequencing technology.
In certain embodiments, a method of generating an error corrected double-stranded consensus sequence is provided. Such a method, also referred to as duplex consensus sequencing (DCS), allows for a quantitative detection of sites of DNA damage. DCS analysis facilitates the detection of DNA damage signatures, in that single stranded DNA mutations that are not present in the complementary strand can be inferred to be artifactual mutations arising from damaged nucleotides. Not only can one correct for these erroneous mutations, but the ability to indirectly infer that damage is present on the DNA could be a useful biomarker (e.g. for cancer risk, cancer metabolic state, mutator phenotype related to defective damage repair, carcinogen exposure, chronic inflammation exposure, individual-specific aging, neurodegenerative diseases etc). The ability to use different polymerases during the first round(s) of PCR to mis-incorporate at damage sites could potentially add even more information. Besides polymerases, other DNA modifying/repair enzymes could be used prior to amplification to convert damage of one sort that doesn't give a specific mutagenic signature into another sort that does with whatever polymerase is used. Alternatively, DNA modifying/repair enzymes could be used to remove damaged bases, and one could sequence both strands of DNA both with and without the enzymatic treatment. Mutations in single-stranded DNA that are seen to be removed by the enzymatic treatment can thus be inferred to be arising due to DNA damage. This could be useful on human nuclear or mtDNA but also might also be useful with model organisms (mice, yeast, bacteria etc), treated with different new damaging agents, facilitating a screen for DNA damaging compounds that would be analogous to the widely used Ames test [52].
The method of generating an error corrected double-stranded consensus sequence may include a first stage termed “single strand consensus sequencing” (SSCS) followed by a second stage of duplex consensus sequencing (DCS). Therefore, the method includes steps of tagging individual duplex DNA molecules with an SMI adaptor molecule, such as those described above; generating a set of PCR duplicates of the tagged DNA molecules by performing a suitable PCR method; creating a single strand consensus sequence from all of the PCR duplicates which arose from an individual molecule of single-stranded DNA. Each DNA duplex should result in two single strand consensus sequences. The work through these three steps conclude the first stage and is termed SSCS.
The method of generating an error corrected double-stranded consensus sequence further comprises the second stance that is termed DCS. The DCS stage includes steps of comparing the sequence of the two single strand consensus sequences arising from a single duplex DNA molecule, and further reducing sequencing or PCR errors by considering only sites at which the sequences of both single-stranded DNA molecules are in agreement. The method that includes the first stage and the second stage termed Duplex Consensus Sequencing (DCS).
The step of tagging of both strands of individual duplex DNA may be accomplished by ligation of degenerate or semi-degenerate complementary DNA sequences; as the complementary nature of the two strands of such a tag sequence allows the two molecules to be grouped together for error correction. Alternatively, as described above, the two duplex DNA strands may be linked by ligation of a U-shaped SMI adaptor molecule, and the two DNA strands can thus both be tagged with a single-stranded SMI tag.
In the method described above, a set of sequenced SMI-DNA products generated in the methods described above may be grouped into families of paired target nucleic acid strands based on a common set of SMI sequences. Then, the paired target nucleic acid strands can be filtered to remove nucleotide positions where the sequences seen on both of the paired partner DNA strands are not complementary. This error corrected double-stranded consensus sequence may be used in a method for confirming the presence of a true mutation (as opposed to a PCR error or other artifactual mutation) in a target nucleic acid sequence. According to certain embodiments, such a method may include identifying one or more mutations present in the paired target nucleic acid strands that have one or more nucleotide positions that disagree between the two strands, then comparing the mutation present in the paired target nucleic acid strands to the error corrected double-stranded consensus sequence. The presence of a true mutation is confirmed when the mutation is present on both of the target nucleic acid strands and also appear in all members of a pared target nucleic acid family.
The accuracy of current approaches to next-generation sequencing is limited due to their dependence on interrogating single-stranded DNA. This dependence makes potential sources of error such as PCR amplification errors and DNA damage fundamentally limiting. However, the complementary strands of a double-stranded DNA molecule (or “DNA duplex”) contain redundant sequencing information (i.e., one molecule reciprocally encoding the sequence information of its partner) which can be utilized to eliminate such artifacts. Limitations related to sequencing single-stranded DNA (e.g., sequencing errors) may therefore be overcome using the methods described herein. This is accomplished by individually tagging and sequencing each of the two strands of a double-stranded (or duplex) target nucleic acid molecule and comparing the individual tagged amplicons derived from one half of a double-stranded complex with those of the other half of the same molecule. Duplex Consensus Sequencing (DCS), significantly lowers the error rate of sequencing. In some embodiments, the DCS method may be used in methods for high sensitivity detection of rare mutant and variant DNA as described further below.
As described above, one approach that has previously been reported for DNA sequencing involves incorporation of a random tag sequence into a PCR primer [36]. This approach results in an improvement in accuracy relative to standard Illumina sequencing, but is fundamentally limited in that it is based upon amplification and sequencing of single-stranded DNA and thus cannot overcome limitations in sensitivity owing to single-stranded DNA damage events. In the methods described herein, PCR duplicates are generated from a single strand of DNA, and the sequences of the duplicates are compared. Mutations are scored only when they are present in multiple replicates of a single starting molecule. The DCS approach overcomes the limitation of previous approaches by considering both DNA strands.
DNA damage should not be a limiting factor in DCS, because miscoding damage events at a single base-pair position occur essentially exclusively on only one of the two DNA strands. For DNA damage to result in an artifactual mutation in DCS, damage would need to be present at the same nucleotide position on both strands. Even if complementary nucleotides in a duplex were both damaged, the damage would need to result in complementary sequencing errors to result in mis-scoring of a mutation. Likewise, spontaneous PCR errors would need to result in complementary mutations at the same position on both strands; with a first-round mutation frequency of Taq polymerase of approximately 10−5 and three possible incorrect bases that could be mis-inserted, the probability of two complementary PCR errors occurring would be 10−5×10−5×1/3=3.3×10−11
According to some embodiments, the sequencing method may be performed using the Illumina or similar platforms including those enumerated above without the use of SMI adaptor molecules, but instead by using the random shear points of DNA as identifiers. For a given DNA sequence seen in sequencing read 1 with a specific set of shear points, the partner strand will be seen as a matching sequence in read two with identical shear points. In practice, this approach is limited by the limited number of possible shear points that overlap any given DNA position. However, according to some embodiments, shear points of a target nucleic acid molecule may be used as unique identifiers to identify double-stranded (or duplex) pairs, resulting in an apparent error frequency at least as low as that seen with traditional sequencing methods, but with a significantly lower loss of sequence capacity. In other embodiments, DCS based on shear points alone may have a role for confirmation that specific mutations of interest are true mutations which were indeed present in the starting sample (i.e. present in both DNA strands), as opposed to being PCR or sequencing artifacts. Overall, however, DCS is most generally applicable when randomized, complementary double-stranded SMI sequences are used. A 24 nucleotide double-stranded SMI sequence was used in the Example described below, which may yield up to 424=2.8×1014 distinct double-stranded SMI sequences. Combining information regarding the shear points of DNA with the SMI tag sequence would allow a shorter SMI to be used, thus minimizing loss of sequencing capacity due to sequencing of the SMI itself.
In certain embodiments, the SMI adaptor molecules may also be used in methods of single-molecule counting for accurate determination of DNA or RNA copy number [38]. Again, since the SMI tags are present in the adaptors, there are no altered steps required in library preparation, which is in contrast to other methods for using random tags for single-molecule counting. Single-molecule counting has a large number of applications including, but not limited to, accurate detection of altered genomic copy number (e.g., for sensitive diagnosis of genetic conditions such as trisomy 21 [47]), for accurate identification of altered mRNA copy number in transcriptional sequencing and chromatin immunoprecipitation experiments, quantification of circulating microRNAs, quantification of viral load of DNA or RNA viruses, quantification of microorganism abundance, quantification of circulating neoplastic cells, counting of DNA-labeled molecules of any variety including tagged antibodies or aptamers, and quantification of relative abundances of different individual's genomes in forensic applications.
In another embodiment, the SMI adaptor molecules may be used in methods for unambiguous identification of PCR duplicates. In order to restrict sequencing analysis to uniquely sequenced DNA fragments, many sequencing studies include a step to filter out PCR duplicates by using the shear points at the ends of DNA molecules to identify distinct molecules. When multiple molecules exhibit identical shear points, all but one of the molecules are discarded from analysis under the assumption that the molecules represent multiple PCR copies of the same starting molecule. However sequence reads with identical shear points can also reflect distinct molecules because there are a limited number of possible shear points at any given genomic location, and with increasing sequencing depth, recurrent shear points are increasingly likely to be seen [48]. Because the use of SMI tags (or “double-stranded SMI sequences”) allows every molecule to be uniquely labeled prior to PCR duplication, true PCR duplicates may be unambiguously identified by virtue of having a common (i.e., the same or identical) SMI sequence. This approach would thereby minimize the loss of data by overcoming the intrinsic limitations of using shear points to identify PCR duplicates.
Importantly, once SMI-containing adaptors are synthesized by a straightforward series of enzymatic steps or are produced through synthesis of a set of oligonucleotides containing complementary tag sequences, they may be substituted for standard sequencing adaptors. Thus, use of DCS does not require any significant deviations from the normal workflow of sample preparation for Illumina DNA sequencing. Moreover, the DCS approach can be generalized to nearly any sequencing platform because a double-stranded SMI tag can be incorporated into other existing adaptors, or for sequencing approaches that do not require adaptors, a double-stranded SMI tag can be ligated onto duplex DNA sample prior to sequencing. The compatibility of DCS with existing sequencing workflows, the potential for greatly reducing the error rate of DNA sequencing, and the multitude of applications for the double-stranded SMI sequences validate DCS as a technique that may play a general role in next generation DNA sequencing.
The following examples are intended to illustrate various embodiments of the invention. As such, the specific embodiments discussed are not to be construed as limitations on the scope of the invention. It will be apparent to one skilled in the art that various equivalents, changes, and modifications may be made without departing from the scope of invention, and it is understood that such equivalent embodiments are to be included herein. Further, all references cited in the disclosure are hereby incorporated by reference in their entirety, as if fully set forth herein.
Materials and Methods
Materials.
Oligonucleotides were from IDT and were ordered as PAGE purified. Klenow exo- was from NEB. T4 ligase was from Enzymatics.
DNA Isolation.
Genomic DNA was isolated from normal human colonic mucosa by sodium iodide extraction (Wako Chemicals USA).
Adaptor Synthesis.
The adaptors were synthesized from two oligos, designated as:
The two adaptor strands were annealed by combining equimolar amounts of each oligo to a final concentration of 50 micromolar and heating to 95° C. for 5 minutes. The oligo mix was allowed to cool to room temperature for over 1 hour. The annealed primer-template complex was extended in a reaction consisting of 40 micromolar primer-template, 25 units Klenow exo-DNA polymerase (New England Biolabs), 250 micromolar each dNTP, 50 mM NaCl, 10 mM Tris-HCl pH 7.9, 10 mM MgCl2, and 1 mM dithiothreitol (DTT) for 1 hour at 37° C. The product was isolated by ethanol precipitation. Due to the partial A-tailing property of Klenow exo-, this protocol results in a mixture of blunt-ended adapters and adapters with a single-nucleotide A hverhang. A single-nucleotide A overhang was added to residual blunt fragments by incubating the adapters with 25 units Klenow exo-, 1 mM dATP, 50 mM NaCl, 10 mM Tris-HCl pH 7.9, 10 mM MgCl2, and 1 mM dithiothreitol (DTT) for 1 hour at 37° C. The product was again ethanol precipitated and resuspended to a final concentration of 50 micromolar.
Sequencing Library Preparation.
3 micrograms of DNA was diluted into 130 microliters of TE buffer (10 mM tris-HCl, pH 8.0, 0.1 M EDTA) and was sheared on the Covaris AFA system with duty cycle 10%, intensity 5, cycles/burst 200, time 20 seconds×6, temperature 4° C. DNA was purified with 2 volumes of Agencourt AMPure XP beads per the manufacturer's protocol. After end-repair with the NEB end-repair kit per the manufacturer's protocol, DNA fragments larger than the optimal range of ˜200-500 bp were removed by adding 0.7 volumes of AMPure XP beads and transferring the supernatant to a separate tube (fragments larger than 500 bp bind to the beads and are discarded). An additional 0.65 volumes of AMPure XP beads were added (this step allows fragments of approximately 200 bp or greater to bind to the beads). The beads were washed and DNA eluted. DNA was then T-tailed in a reaction containing 5 units Klenow exo-, 1 mM dTTP, 50 mM NaCl, 10 mM Tris-HCl pH 7.9, 10 mM MgCl2, 1 mM. The reaction proceeded for 1 hour at 37 C. DNA was purified with 1.2 volumes of AMPure XP beads. The custom adaptors were ligated by combining 750 ng of T-tailed DNA with 250 pmol adaptors in a reaction containing 3000 units T4 DNA ligase, 50 mM Tris-HCl pH 7.6, 10 mM MgCl2, 5 mM DTT, 1 mM ATP. The reaction was incubated 25 C for 15 minutes, and purified with 1.2 volumes of AMPure XP beads.
Pre-Capture Amplification.
375 ng adaptor-ligated DNA was PCR amplified with primers AATGATACGGCGACCACCGAG (SEQ ID NO:3) and GTGACTGGAGTTCAGACGTGTGC (SEQ ID NO:4) using the Kappa high-fidelity PCR kit for 8 cycles with an annealing temperature of 60 C. The product was purified with 1.2 volumes of AMPure XP beads.
DNA Capture.
Target capture was performed with the Agilent SureSelect system per the manufacturer's recommendations, except that capture volumes were performed at one-half of the standard volume. The capture set targeted an arbitrary 758 kb region of the genome consisting of both coding and noncoding sequences. Capture baits were 120 nt in length, and were prepared with the Agilent eArray tool with 3× tiling.
Post-Capture Amplification.
Captured DNA was amplified with PCR primers AATGATACGGCGACCACCGAG (SEQ ID NO:3) and CAAGC AGAAGACGGCATACGAGATXXXXXXGTGACTGGAGTTCAGACGTGTGC (SEQ ID NO:5) where XXXXXX indicates the position of a fixed multiplexing barcode sequence). 2 0 fmol of DNA was used per lane for sequencing on an Illumina HiSeq 2000.
Data Processing.
Reads with intact SMI adaptors include a 12 nucleotide random sequence, followed by a 5 nucleotide fixed sequence. These reads were identified by filtering out reads that lack the expected fixed sequence at positions 13-17. The SMI sequence from both the forward and reverse sequencing reads (i.e., the first and second degenerate n-mer sequences) was computationally added to the read header, and the fixed sequence removed. The first 4 nucleotides located following the adaptor sequence were also removed due to the propensity for ligation and end-repair errors to result in an elevated error rate near the end of the DNA fragments. Reads having common (i.e., identical) SMI sequences were grouped together, and were collapsed to generate a consensus read. Sequencing positions were discounted if the consensus group covering that position consisted of fewer than 3 members, or if fewer than 90% of the sequences at that position in the consensus group had the identical sequence. Reads were aligned to the human genome with the Burrows-Wheeler Aligner (BWA). The consensus sequences were then paired with their strand-mate by grouping each 24 nucleotide tag of form AB in read 1 with its corresponding tag of form BA in read 2. Resultant sequence positions were considered only when information from both DNA strands was in perfect agreement. An overview of the data processing workflow is as follows:
Code for carrying out the workflow may be pre-existing or may involve programming within the skill of those in the art. In some embodiments, however, the Python code, which is illustrated in
Overview
To overcome limitations in the sensitivity of variant detection by single-stranded next-generation DNA sequencing, an alternative approach to library preparation and analysis was designed, which is known herein as Duplex Consensus Sequencing (DCS) (
Following tagging with a double-stranded SMI and PCR amplification, a family of molecules is obtained that arose from a single DNA molecule; members of the same PCR “family” are then grouped together by virtue of having a common (i.e., the same) SMI tag sequence. The sequences of uniquely tagged PCR duplicates are subsequently compared in order to create a PCR consensus sequence. Only DNA positions that yield the same DNA sequence in a specified proportion of the PCR duplicates in a family, such as 90% of the duplicates in one embodiment, are used to create the PCR consensus sequence. This step filters out random errors introduced during sequencing or PCR to yield the PCR consensus sequences, each of which derives from an individual molecule of single-stranded DNA. This set of PCR consensus sequences are called single strand consensus sequences (SSCSs).
Next, PCR consensus sequences arising from two complementary strands of duplex DNA can be identified by virtue of the complementary SMIs (
Results
In order to label or tag each of the strands of duplex DNA with unique complementary tags, adaptors which contain the standard sequences required for the Illumina HiSeq system were synthesized, but with addition of a double-stranded, complementary SMI sequence (or “tag”) of 12 random nucleotides (or a random “degenerate sequence”) per strand. Target DNA molecules having a random SMI sequence n-mer that is 12 nucleotides in length on each end will therefore have a unique 24 nucleotide SMI sequence. The adaptors were prepared (
DNA for sequencing was sheared and end-repaired by standard methods, with size-selection for fragments in the range of ˜200-500 bp by size-selective binding to Ampure XP beads. Standard Illumina library preparation protocols involve ligating A-tailed DNA to T-tailed adaptors. However, because A-tailed adaptors were used, the DNA was T-tailed by incubating the end-repaired DNA with Klenow exo- DNA polymerase and 1 mM dTTP. The adaptor-ligated library was PCR amplified and subjected to SureSelect capture, with targeting of an arbitrary 758 kb portion of the genome (DNA coordinates available upon request). The efficiency of adaptor ligation, PCR amplification, DNA capture, and sequencing were comparable to those seen with standard library preparation methods (data not shown). Although Agilent Sure Select probes are used in this example, any suitable method of DNA selection may be used to capture particular target double-stranded DNA sequences. For example, selection and capture may be accomplished by any selection by hybridization method (e.g., Agilent SureSelect, Primer Extension Capture, exploitation of biotinylated PCR amplicons as bait, Agilent HaloPlex) wherein probes that target the desired double-stranded DNA sequence may be recovered by an in-array capture (using probes immobilized on glass slides) or by affinity using magnetic beads in an in-solution capture. In addition, mitochondrial and some other forms of DNA may be isolated by size selection. Alternatively, in some embodiments, no enrichment is performed.
This protocol was used to sequence DNA isolated from normal colonic mucosa. Mutations were initially scored without consideration of the SMI sequences. PCR duplicates were filtered out with samtools rmdup, a standard tool which uses the shear points of DNA molecules to identify PCR duplicates, as molecules arising from duplicated DNA will have shared shear points. In order to focus specifically on non-clonal mutations, only those positions in the genome with at least 20× coverage and at which fewer than 5% of reads differed from the hg19 reference sequence were considered. This approach resulted in 70.9 million nucleotides of sequence data and 56,890 mutations, indicating an overall mutation frequency of 8.03×10−4, in accord with the error rate of Illumina next-generation sequencing of ˜0.1-1% [32].
Next, the SMI tags were used to group together PCR duplicates that arose from individual single-stranded DNA molecules and to create a consensus sequence from the family of duplicates. At least 3 PCR duplicates were required, with at least 90% agreement in sequence among all duplicates, to consider a site for mutations. Scoring the mutation frequency as above, again considering only sites with a minimum of 20× coverage and with <5% of reads differing from reference, resulted in 145 million nucleotides of sequence with 6,508 mutations and an overall mutation frequency of 4.47×10−5, consistent with prior reports [36]. Notably, far more nucleotides of DNA sequence were obtained in this approach (145 million) than in the standard Illumina sequencing approach (70 million) detailed above which is dependent on use of the shear points of single-ended reads to identify PCR duplicates. The improved sequence coverage arose from use of the SMI to identify PCR duplicates, because identifying PCR duplicates by consideration of uniquely sheared DNA ends is fundamentally limited by the small number of possible shear points that overlap a given position of the genome and the propensity for specific genomic regions to be more readily undergo shearing. Thus filtering PCR duplicates by using shear points resulted in discarding a large portion of the reads.
Finally, the complementary nature of the double-stranded SMI sequences was used to identify pairs of consensus groups that arose from complementary DNA strands. Sequence reads were considered only when the read data from each of the two strands is in perfect agreement. In a pilot experiment, after grouping of PCR duplicates as above, 29,409 SMI partner pairs were found, indicative that fewer than 1 of tags had their corresponding partner tag present in the library. The low recovery of tag pairs was most likely due to inadequate amplification of the starting DNA library. Among these tag-pairs, 24,772 duplex consensus strands were identified with an average strand length of 82 nucleotides, resulting in 2 million nucleotides of DNA consensus sequence. The sequences of the paired duplex strands disagreed at 3,585 of the nucleotide positions, indicative of single-stranded errors (i.e. PCR or sequencing errors); these sites of disagreement were removed, leaving only bases at which the sequence of both duplex strands were in perfect agreement. Next, as above, analysis of mutation frequencies was restricted to sites with at least 10× coverage and at which fewer than 10% of reads disagreed from the hg19 reference sequence. Because the 2 million nucleotides of read data were spread across a 758 kb target, our average depth was only ˜3×. Thus only 14,464 nucleotides of DNA sequence corresponded to sites with at least 10× depth. Among these sites, zero mutations were seen. To increase the number of tag pairs considered, analysis described above was repeated, but PCR duplicates were grouped with a minimum of only 1 duplicate per site. This resulted in 28,359 nucleotides of DNA sequence with at least 10× depth. Again, no mutations were detected.
Current experiments are being performed on vastly smaller target DNA molecules (ranging from ˜300 bp to ˜20 kb in size). Use of smaller DNA targets will allow for much greater sequencing depth, and far more accurate assessment of the background mutation rate of the assay. In addition, the protocol has been modified to incorporate a greater number of PCR cycles initiated off a smaller number of genome equivalents, which will increase the fraction of tags for which both of the partner tag strands have been sufficiently amplified to be represented in the final sequence data. Indeed, among the 3.6 million SMIs present in our initial library which underwent PCR duplication, 1.5 million of the SMIs were present only once, indicating insufficient amplification of the DNA due in part to the low number of PCR cycles used.
Materials and Methods
In addition to those described in Example 1 above, the following materials and methods were also used.
DNA Isolation.
Mitochondrial DNA was isolated as previously described (4).
Data Processing.
The entire human genome sequence (hg19) was used as reference for the mitochondrial DNA experiment, and reads that mapped to chromosomal DNA were removed. Reads sharing identical tag sequences were then grouped together and collapsed to consensus reads. Sequencing positions were discounted if the consensus group covering that position consisted of fewer than three members or if fewer than 90% of the sequences at that position in the consensus group had the identical sequence. A minimum group size of three was selected because next-generation sequencing systems have an average base calling error rate of ˜1/100. Requiring the same base to be identified in three distinct reads decreases the frequency of single-strand consensus sequence (SSCS) errors arising from base-call errors to (1/100)3=1×10−6, which is below the frequency of spontaneous PCR errors that fundamentally limit the sensitivity of SSCSs. The requirement for 90% of sequences to agree to score a position is a highly conservative cutoff. For example, with a group size of eight, a single disagreeing read will lead to 87.5% agreement and the position will not be scored. If all groups in an experiment are of size nine or less, this cutoff will thus require perfect agreement at any given position to score the position. Further development of our protocol may allow for less stringent parameters to be used to maximize the number of SSCS and duplex consensus sequence (DCS) reads that can be obtained from a given experiment.
Results
Having established the methodology for Duplex Sequencing with M13mp2 DNA, which is a substrate for which the mutation frequency and spectrum are fairly well established, it was desired to apply the approach to a human DNA sample. Thus, mitochondrial DNA was isolated from human brain tissue and sequenced the DNA after ligation of Duplex Sequencing adapters. A standard sequencing approach with quality filtering for a Phred score of 30 resulted in a mutation frequency of 2.7×10-3, and SSCS analysis yielded a mutation frequency of 1.5×10-4. In contrast, DCS analysis revealed a much lower overall mutation frequency of 3.5×10-5 (
When the distribution of mutations throughout the mitochondrial genome is considered, the quality filtered reads (analyzed without consideration of the tags) have many artifactual errors, such that identification of mutational hotspots is difficult or impossible (
SSCS analysis produced a strong mutational bias, with a 130-fold excess of G→T relative to C→A mutations (
Materials and Methods
In addition to those described in the Examples above, the following materials and methods were also used to demonstrate the capability of DCS analysis to remove sequencing errors
Sequencing Library Preparation.
Genomic DNA was isolated from a derivative of Saccharomyces cerevisiae strain SC288 by standard methods. The DNA was randomly sheared by the Covaris AFA system, followed by end-repair, A-tailing, and ligation of Illumina TruSeq DNA sequencing adaptors, all by standard library preparation methods. The resultant sequence data consisted of an average 32.5 fold depth of the 12 megabase S. cerevisiae genome.
Data Analysis.
The first 10 nucleotides of each sequencing read pair, corresponding to the randomly sheared DNA ends, were combined, such that the first 10 nucleotides of read 1, referred to as A, was combined with the first 10 nucleotides of read 2, referred to as B, to yield an SMI tag of form AB. Reads were grouped according to SMI sequence, and nucleotide reads were considered only if they agreed among at least 90% of family members sharing a given tag sequence. For DCS analysis, a tag of form AB1 is partnered with the corresponding tag of form BA2, and nucleotide positions are considered only when the sequence is in agreement among read pairs with matching tags AB1 and BA2.
Results
In order to demonstrate the capability of DCS analysis to remove sequencing errors, a sequencing library was prepared under standard conditions with commercially available sequencing adaptors, and the randomly sheared DNA ends were used as SMI's. First, reads were grouped by SMI with a minimum family size of 1 member. Considering only sites with a minimum of 20× coverage and with <5% of reads differing from reference, this analysis resulted in 644.8 million nucleotides of sequence data and 2,381,428 mutations, yielding an overall mutation frequency of 3.69×103.
The data was then subjected to DCS analysis with the SMI tags, searching for tags of form AB1 that have partner tags of form BA2, and considered only positions at which the sequence from the two strands was in perfect agreement. 3.1% of the tags had a matching partner present in the library, resulting in 2.9 million nucleotides of sequence data. The sequences of the duplex strands were not complementary at 40,874 nucleotide positions; these disagreeing positions, representing likely sequencing or PCR errors, were removed from analysis. Again considering positions with at least 20× coverage and <5% of reads differing from reference, 3.0 million nucleotides of sequence data and 157 mutations were obtained, with an overall mutation frequency of 5.33×10−5, indicative of removal of >98% of mutations seen in raw analysis and thereby demonstrating the capability of DCS to lower the error rate of DNA sequencing.
To compare this result to the method of Kinde et al. [36], reads were grouped into families by SMI tag as before but filtered for families with a minimum of 3 members. This resulted in 1.4 million nucleotides of sequence data and 61 mutations, with an overall mutation frequency of 4.25×10−5. Thus, the method of Kinde et al., with a minimum family size of 3, resulted in less than half as much resultant sequence data after filtering than was obtained by DCS with a minimum family size of 1. Thus, DCS lowered the error rate of sequencing to a comparable degree to a method considered state-of-the-art, but with less loss of sequencing capacity.
Discussion
It was demonstrated that DCS analysis, using sheared DNA ends as unique molecular identifiers, results in a lowering of the apparent error rate of DNA sequencing. As this proof-of-concept experiment was performed on a library that was not optimized to maximize recovery of both strands, there were not sufficient strand-pairs recovered to perform DCS analysis with a minimum family size of >1 member. Requiring family sizes >1 is expected to further reduce sequencing errors. Moreover, this analysis was limited in that it did not include ligation of degenerate SMI tag sequences; owing to the limited number of shear points flanking any given nucleotide position, use of shear points as SMIs limits the number of unique molecules that can be sequenced in a single experiment. The use of shear points as SMIs in conjunction with an exogenously ligated SMI tag sequence would allow for increased depth of sequencing at any given nucleotide position.
In addition to those described in Examples 1 and 2 above, the following materials and methods were also used.
Materials and Methods
Construction of M13mp2 Variants.
M13mp2 gapped DNA encoding the LacZ a fragment was extended by human DNA polymerase 6 [2] and the resultant products were transformed into Escherichia coli and subjected to blue-white color screening as previously described [3]. Mutant plaques were sequenced to determine the location of the mutation resulting in the color phenotype. A series of mutants, each differing from wild type by a single nucleotide change, were then mixed together with wild-type M13mp2 DNA to result in a single final mixture with distinct mutants represented at ratios of 1/10 (G6267A), 1/100 (T6299C), 1/1,000 (G6343A), and 1/10,000 (A6293T).
Oxidative Damage of M13mp2 DNA.
Induction of DNA damage was performed by minor modifications to a published protocol [5]: 300 ng of M13mp2 double-stranded DNA was incubated in 10 mM sodium phosphate buffer, pH 7.0, in the presence of 10 μM iron sulfate and 10 μM freshly diluted hydrogen peroxide. Incubation proceeded for 30 min at 37° C. in open 1.5-mL plastic microcentrifuge tubes.
DNA Isolation.
M13mp2 DNA was isolated from E. coli strain MC1061 by Qiagen Miniprep. To allow for greater sequencing depth at a defined region of the M13mp2 genome, an 840-bp fragment was enriched by complete digestion with the restriction enzymes Bsu36I and NaeI (New England Biolabs), followed by isolation of the fragment on an agarose gel by the Recochip system (Takara Bio).
Duplex Consensus Sequencing of M13 DNA Removes Artifactual Sequencing Errors.
The spontaneous mutation rate of M13mp2 DNA has been well established by a number of exquisitely sensitive genetic assays to be 3.0E-6 [53], that is, an average of one spontaneous base substitution error for every 330,000 nucleotides. Thus this substrate is well suited as a control for determining the background error frequency of DNA sequencing. M13mp2 DNA was sheared and ligated to adaptors containing double-stranded complementary SMI sequences by standard protocols, and was subjected to deep sequencing on an Illumina HiSeq 2000 followed by Consensus Sequencing analysis (
Analysis of the data by standard methods (i.e., without consideration of the double stranded SMI sequences) resulted in an error frequency of 3.8E-03, more than one thousand fold higher than the true mutation frequency of M13mp2 DNA. This indicates that >99.9% of the apparent mutations identified by standard sequencing are in fact artifactual errors.
The data were then analyzed by Single Strand Consensus Sequencing (SSCS), using the unique SMI tag affixed to each molecule to group PCR products together in order to create a consensus of all PCR products that came from an individual molecule of single-stranded DNA. This resulted in a mutation frequency of 6.4E-OS, suggesting that −98% of sequencing errors are corrected by SSCS.
Next, the data were subjected to Duplex Consensus Sequencing (DCS), which further corrects errors by using the complementary SMI tags to compare the DNA sequence arising from both of the two strands of a single molecule of duplex DNA. This approach resulted in a mutation frequency of 2.SE-06, in nearly perfect agreement with the true mutation frequency of M13mp2 DNA of 3.0E-06. The number of nucleotides of DNA sequence obtained by a standard sequencing approach, and after SSCS and DCS analysis, may be found in Table 1 below.
For an artifactual error to be scored by DCS, complementary artifactual errors must occur on both strands of a molecule of duplex DNA. Thus the background (artifactual) error frequency of DCS may be calculated as: (probability of error on one strand)*(probability of error on other strand)*(probability that both errors are complementary).
As the background error frequency of SSCS in this experiment was −6E-S, the background error frequency of DCS can be calculated as 6E-S*6E-S*(1/3)=1.2E-9. This represents a greater than 3 million fold improvement over the error rate of 3.SE-03 that was obtained by a standard sequencing approach.
Consensus Sequencing Reveals Likely Sites of DNA Damage
M13mp2 DNA was sequenced as detailed above, with DCS adaptors containing double-stranded complementary SMIs. The spectrum of mutations obtained with SSCS was determined. Data was filtered to consist of forward-mapping reads from Read 1, i.e. sequencing of the reference strand, and reverse-mapping reads from Read 1, i.e. sequencing of the anti-reference strand. True mutations would result in an equal balance between mutations on the reference strand and their complementary mutation on the anti-reference strand.
However, SSCS analysis revealed a large number of single-stranded G→T mutations on reads mapping in the forward orientation to the reference genome, with a much smaller number of C→A mutations mapping to the reverse orientation. The spectrum of mutations identified by both SSCS and DCS analysis were examined relative to literature reference values [53] for the M13mp2 substrate (
Base-specific mutagenic DNA damage is a likely explanation of these imbalances. Excess G→T mutations are consistent with the oxidative product 8-oxo-guanine (8-oxo-G) causing first round PCR errors and artifactual G→T mutations. DNA polymerases, including those commonly used in PCR, have a strong tendency to insert adenine opposite 8-oxo-G [45, 54], and misinsertion of A opposite 8-oxo-G would result in erroneous scoring of a G→T mutation. Likewise, the excess C→T mutations are consistent with spontaneous deamination of cytosine to uracil [47], a particularly common DNA damage event that results in insertion during PCR of adenine opposite uracil and erroneous scoring of a C→T mutation.
To determine whether the excess G→T mutations seen in SSCSs might reflect oxidative DNA damage at guanine nucleotides, before sequencing library preparation M13mp2 DNA was incubated with the free radical generator hydrogen peroxide in the presence of iron, a protocol that induces DNA damage [55]. This treatment resulted in a substantial further increase in G→T mutations by SSCS analysis (
Furthermore, relative to the literature reference values, DCS analysis results in a lower frequency of G→T/C→A and C→T/G→A mutations (
Consensus Sequencing Accurately Recovers Spiked-in Control Mutations.
A series of M13mp2 variants were constructed which contain known single base substitutions. These variants were then mixed together at known ratios, and the mixture was prepared for sequencing with DCS adaptors with double-stranded complementary SMIs and was sequenced on an Illumina HiSeq 2000. The data was then analyzed by consensus sequencing (
In contrast, when the data is analyzed by Single Strand Consensus Sequencing (SSCS) with −20,000 fold depth, accurate recovery of mutant sequence is seen down to one mutant molecule per 10,000 wild type molecules. Duplex Consensus Sequencing (DCS), which was not performed on this sample, would allow for detection of even rarer mutations.
The references, patents and published patent applications listed below, and all references cited in the specification above are hereby incorporated by reference in their entirety, as if fully set forth herein.
This application claims priority to U.S. Provisional Patent Application No. 61/613,413, filed Mar. 20, 2012; U.S. Provisional Patent Application No. 61/625,623, filed Apr. 17, 2012; and U.S. Provisional Patent Application No. 61/625,319, filed Apr. 17, 2012; the subject matter of all of which are hereby incorporated by reference as if fully set forth herein.
The present invention was made with government support under Grant Nos. RO1 CA115802 and RO1 CA102029 awarded by the National Institutes of Health. The Government has certain rights in the invention.
| Number | Date | Country | |
|---|---|---|---|
| 61613413 | Mar 2012 | US | |
| 61625623 | Apr 2012 | US | |
| 61625319 | Apr 2012 | US |
| Number | Date | Country | |
|---|---|---|---|
| Parent | 14386800 | Sep 2014 | US |
| Child | 15660785 | US |
