Patent Application
20220145380
Traditional genetic sequencing methodologies, such as whole genome (WGS) and whole exome (WES), have focused on the important contribution of germline mutations that are present in all cells throughout the human body. However, recent studies have shown numerous examples of mutations occurring after fertilization (i.e. postzygotic mutations), which are only present in a fraction of the cells. Postzygotic mutations, or somatic mutations, have been heavily studied in cancers where clinical diagnostic testing for somatic mutations in tumor and blood samples are becoming a standard practice due to improved detection sensitivities when most cells in the sample carry a given mutation.
Beyond technical errors, an important consideration for skewed alternate allelic fraction (AAFs), false negatives, and false positives are allelic imbalances caused by inherent differences in the genome content around a mutation. These issues, such as additional mutations, repeat content, methylation, or copy number changes, can have dramatic impacts on AAFs, resulting in the commonly recognized issue of allelic dropout. To avoid allelic dropout, many methods avoid placing primers in areas with known genetic variation in the general population. However, these methods remain susceptible to allelic skewing from ultra-rare or private alleles and other loci specific causes of allelic imbalance. Cost-effective methods are needed for the detection and characterization of rare alleles and other genetic variants.
As described below, the present disclosure features methods for detecting and quantifying genetic variants in a sample.
In one aspect of the present disclosure, a method is provided for determining alternate allele frequency, the method involves performing two or more parallel amplification reactions on a single sample, thereby generating overlapping amplicons, where each amplification reaction includes a unique pair of forward and reverse primers, where the forward or reverse primer includes an index sequence, and where the forward and reverse primers include different adapter sequences. The method also involves sequencing the overlapping amplicons to produce sequence reads, segregating the sequencing reads into bins by index sequence, and detecting the presence or absence of one or more genetic variants within sequencing reads within a bin, where the frequency of detection of the variant determines the alternate allele frequency.
Another aspect provides a method for determining alternate allele frequency, the method involves a) performing three amplification reactions on a single sample, thereby generating three overlapping amplicons, where each amplification reaction includes a unique pair of forward and reverse primers, where each primer includes a nucleic acid sequence complementary to a portion of a target nucleic acid sequence, where the forward or reverse primer includes an index sequence, where the forward and reverse primers include different adapter sequences at or near the 5′ terminus of the primer and upstream of the sequence complementary to the target, and where at least one adapter sequence is complementary to a nucleic acid sequence used in sequencing; b) sequencing the overlapping amplicons to produce sequence reads; c) segregating the sequencing reads into bins by index sequence; and d) detecting the presence or absence of one or more genetic variants within sequencing reads within a bin, where the frequency of detection of the variant determines the alternate allele frequency.
Another aspect of the present invention provides a method for method for determining alternate allele frequency, the method involving a) performing three amplification reactions on a single sample, thereby generating three overlapping amplicons, where each amplification reaction includes a unique pair of forward and reverse primers, where the forward or reverse primer comprises an index sequence and/or a unique molecular identifier (UMI); and each primer includes i. a nucleotide sequence complementary to a portion of a target nucleic acid sequence; ii. an adapter at or near its 5′ terminus, where the adapter is upstream of the sequence complementary to the target and wherein the forward and reverse primers include different adapter sequences, and where at least one adapter sequence is complementary to a nucleic acid sequence used in sequencing; b) sequencing the overlapping amplicons to produce sequence reads; c) segregating the sequencing reads into bins by index sequence; d) detecting the UMI and removing duplicate reads from the bin, where the detecting can be simultaneous with step c or subsequent to step c; and e) detecting the presence or absence of one or more genetic variants within sequencing reads within a bin, where the frequency of detection of the variant determines the alternate allele frequency.
In some embodiments, the methods disclosed herein further involve pooling the amplicons prior to sequencing. In some embodiments of the methods disclosed herein, sequencing the amplicons involves contacting the amplicons with a nucleic acid complementary to the adapter sequence. In some embodiments, the amplicons include a nucleotide having a label, and in some embodiments, the label is biotin. In some embodiments, the methods disclosed herein also involve contacting the label with a capture agent that specifically binds the label. In some embodiments, the methods also involve enzymatically digesting the primers. In some embodiments of the present disclosure, the methods also involve amplifying the amplicons, thereby generating enriched populations of amplicons. In some embodiments, the genetic variation to be detected is known or unknown. In some embodiments, the genetic variant has an alternate allele fraction of at least 0.1%. In some embodiments, the genetic variant has an alternate allele fraction of at least 0.025%. In some embodiments, the genetic variant is a mosaic variant. In some embodiments, detection of the genetic variant identifies the presence of a disease or a predisposition to a disease in a subject from whom the sample was derived. In some embodiments, the disease is cancer. In some embodiments, the sample includes circulating tumor cells or cell free DNA. In some embodiments, the genetic variant originated from a somatic event or a germline event. In some embodiments, the alternate allele frequency is compared to the allele frequency of a reference sample to determine if the subject's disease is progressing, regressing, or in remission. In some embodiments, the methods further involve averaging the alternate allele frequencies determined for each bin. In some embodiments, the methods further involve determining the error rate of the nucleic acid sequences flanking the alternate allele.
Methods defined by the present disclosure were performed in connection with the examples provided below. Other features and advantages of the disclosure will be apparent from the detailed description and from the claims.
Unless defined otherwise, all technical and scientific terms used herein have the meaning commonly understood by a person skilled in the art to which this disclosure relates. The following references provide one of skill with a general definition of many of the terms used in this disclosure: Singleton et al., Dictionary of Microbiology and Molecular Biology (2nd ed. 1994); The Cambridge Dictionary of Science and Technology (Walker ed., 1988); The Glossary of Genetics, 5th Ed., R. Rieger et al. (eds.), Springer Verlag (1991); and Hale & Marham, The Harper Collins Dictionary of Biology (1991). As used herein, the following terms have the meanings ascribed to them below, unless specified otherwise.
As used herein, “adapter” refers to a nucleic acid sequence in an amplification primer that is complementary to the sequence of a nucleic acid molecule used to prime downstream sequencing reactions.
The term “allelic dropout” refers to the loss of one allele during amplification, resulting in apparent homozygosity. Nucleotide variation, cytosine methylation, and nucleic acid structure in the primer binding site of only one allele can cause allelic dropout when primer binding to the to the primer binding site is inhibited or reduced. For example, G-quadruplexes (secondary structures formed from stacks of G-quartets) present in the primer binding sites of an allele can prevent efficient priming of the template nucleic acid and lead to allelic dropout.
By “alternative allele” is meant an allele other than a reference allele. An alternative allele will have genetic variation that is not present in the reference allele. In some embodiments, a reference allele is a wildtype allele. A reference allele may differ between different populations, races, or ethnicities. Genetic variation present in an alternative allele can be nucleotide variation (i.e., a transition or a transversion), an insertion, or a deletion. An alternative allele may have a silent variant or mutation, a missense variant or mutation, or a nonsense variant or mutation.
By “alternative allele fraction” is meant the frequency of an allele, other than a reference allele, in a population of cells in an individual. The alternative allele fraction is often less than that of the reference allele fraction, especially when the reference allele is a wildtype allele.
By “amplicon” is meant the product of an amplification reaction.
By “amplification bias” is meant a tendency for a nucleic acid amplification reaction to yield a particular amplicon. Amplification bias is often associated with inefficient primer binding. For example, if a primer's nucleic acid sequence is less complementary to the sequence of a template nucleic acid, the primer will be less likely to bind to the template than a primer having a more complementary sequence. Variants present in the primer binding site of a template nucleic acid may result in conformational or structural changes to the nucleic acid molecule that inhibit primer binding. Other variants or modifications (e.g., methylated nucleic acid residues) present in the primer binding site or elsewhere in the nucleic acid molecule can also cause to amplification bias. Amplification bias may result in underrepresentation of an allele or allelic dropout.
By “analog” is meant a molecule that is not identical, but has analogous functional or structural features to a naturally occurring molecule. For example, a polynucleotide analog retains the biological activity of a corresponding naturally-occurring polynucleotide while having certain modifications that enhance the analog's function relative to a naturally occurring polynucleotide. Such modifications could increase the polynucleotide's affinity for DNA, half-life, and/or nuclease resistance, an analog may include an unnatural nucleotide or amino acid.
By “bin” is meant a collection of sequencing reads that are substantially identical. In some instances, a bin comprises sequences reads that have the same index sequence or UMI sequence.
The phrase “biological sample” as used herein refers to a sample taken from a biological source and includes, but is not limited to, blood, serum, plasma, sputum, lavage fluid, cerebrospinal fluid, urine, semen, sweat, tears, tissue biopsy, and saliva. As used herein, the terms “blood,” “plasma,” and “serum” expressly encompass fractions or processed portions thereof.
In this disclosure, “comprises,” “comprising,” “containing” and “having” and the like can have the meaning ascribed to them in U.S. Patent law and can mean “includes,” “including,” and the like; “consisting essentially of” or “consists essentially” likewise has the meaning ascribed in U.S. Patent law and the term is open-ended, allowing for the presence of more than that which is recited so long as basic or novel characteristics of that which is recited is not changed by the presence of more than that which is recited, but excludes prior art embodiments.
By “demultiplex” is meant a process in which sequence reads generated from different amplicons are segregated into groups based on at least one characteristic unique to each group. For example, the index sequence of a primer can be used to segregate the sequence reads.
The term “denaturing,” as contemplated herein, refers to removing impediments to primer binding from a nucleic acid. For example, denaturing includes removing conformational or structural properties of a nucleic acid or separating a nucleic acid duplex into single strands. Denaturing is facilitated by exposing the duplex to at least one denaturing condition or agent. Denaturing conditions are well known in the art. In one embodiment, a nucleic acid duplex is denatured by exposing it to a temperature that is above the melting temperature (Tm) of the duplex. In certain embodiments, a nucleic acid may be denatured by exposing it to a temperature of at least 90° C. for a sufficient amount of time to denature the nucleic acid molecule. In some embodiments, a denaturing agent may include a chemical additive that facilitates denaturation, for example, sodium hydroxide or urea.
“Detect” refers to discovering or identifying the presence, absence, or amount of an analyte (e.g., genetic variation) to be detected.
By “detectable label” is meant a composition that when linked to a molecule of interest renders the latter detectable, via spectroscopic, photochemical, biochemical, immunochemical, or chemical means. For example, useful labels include radioactive isotopes, magnetic beads, metallic beads, colloidal particles, fluorescent dyes, electron-dense reagents, enzymes (for example, as commonly used in an ELISA), biotin, digoxigenin, or haptens.
“DMSO” refers to dimethyl sulfoxide, which has the following structure:
The term “enrich,” as used herein, refers to the process of further amplifying nucleic acid amplicons. In some embodiments, enrichment of nucleic acid amplicon allows for more efficient detection and quantifying of genetic variants having very low alternative allele frequency relative to detecting and quantifying genetic variants with very low alternative allele frequency in non-enriched nucleic acid amplicons.
By “GC buffer” is meant a reagent designed to optimize the ionic environment of an amplification reaction of a nucleic acid molecule having an enriched guanine/cytosine sequence.
“Germline allele” means an allele specific to germ cells or progenitors thereof.
“Hybridization” means hydrogen bonding, which may be Watson-Crick, Hoogsteen or reversed Hoogsteen hydrogen bonding, between complementary nucleobases. For example, adenine and thymine are complementary nucleobases that pair through the formation of hydrogen bonds.
By “index sequence” or “barcode” is meant a portion of a nucleic acid molecule that allows grouping or demultiplexing of sequencing reads. For example, an index sequence enables the segregation of sequence reads into bins, wherein each bin comprises sequence reads of amplicons generated from the primer pair having the index sequence. In some embodiments, each primer pair used in the presently disclosed methods has a unique index sequence.
As used herein, “interrogate” refers to obtaining nucleotide sequence information for a nucleic acid molecule.
The terms “isolated,” “purified,” or “biologically pure” refer to material that is free to varying degrees from components which normally accompany it as found in its native state. “Isolate” denotes a degree of separation from original source or surroundings. “Purify” denotes a degree of separation that is higher than isolation. A “purified” or “biologically pure” nucleic acid is sufficiently free of other materials such that any impurities do not materially affect the biological properties of the nucleic acid or cause other adverse consequences. That is, a nucleic acid of this disclosure is purified if it is substantially free of cellular material, viral material, or culture medium. Purity and homogeneity are typically determined using analytical chemistry techniques, for example, polyacrylamide gel electrophoresis or high-performance liquid chromatography. The term “purified” can denote that a nucleic acid gives rise to essentially one band in an electrophoretic gel.
By “isolated polynucleotide” is meant a nucleic acid (e.g., a DNA) that is free of the genes which, in the naturally-occurring genome of the organism from which the nucleic acid molecule of the disclosure is derived, flank the gene. The term therefore includes, for example, a recombinant DNA that is incorporated into a vector; into an autonomously replicating plasmid or virus; or into the genomic DNA of a prokaryote or eukaryote; or that exists as a separate molecule (for example, a cDNA or a genomic or cDNA fragment produced by PCR or restriction endonuclease digestion) independent of other sequences. In addition, the term includes an RNA molecule that is transcribed from a DNA molecule, as well as a recombinant DNA that is part of a hybrid gene encoding additional polypeptide sequence.
“Isothermal” refers to a process incubated at about a constant temperature. For example, some isothermal amplification reactions are carried out at about 65° C. An isothermal temperature may depart from an intended temperature by not more than about 10% or 5° C., whichever is greater. An isothermal reaction may include an initial incubation at a higher temperature (“a hot start”). A hot start may comprise incubating the amplification reaction at a temperature sufficient to denature a region of interest on a nucleic acid molecule or to active a reagent (i.e., a polymerase).
By “marker” is meant any protein or polynucleotide associated with a disease or disorder.
As used herein, “mosaic” refers to two or more cells or populations of cells with different genotypes within an individual subject. For example, “somatic mosaicism” refers to two or more genotypically distinct somatic cells or populations of somatic cells in an individual. “Germline mosaicism” occurs when two or more genotypically distinct germ cells or populations of germ cells are present in an individual. Germline mosaicism generally arises after a mutation gives rise to a genotypically distinct gamete.
The term “Next Generation Sequencing (NGS)” refers to massive parallel sequencing of clonally amplified molecules or single nucleic acid molecules. “Massive parallel sequencing” refers to simultaneously performing more than 1000 separate, parallel sequencing reactions. Non-limiting examples of NGS include sequencing-by-synthesis using reversible dye terminators, sequencing-by-ligation, and electronic detection sequencing methods. Electronic detection sequencing methods include those used in the Ion Torrent sequencing strategy (ThermoFisher Scientific) or MiSeq platform (Illumina), wherein changes in pH are detected when a nucleotide is incorporated into a nucleic acid strand resulting in release of a hydrogen ion.
The terms “nucleic acid” and “nucleic acid molecule,” are used interchangeably herein and refer to a compound comprising a nucleobase and an acidic moiety, e.g., a nucleoside, a nucleotide, or a polymer of nucleotides. Typically, polymeric nucleic acids, e.g., nucleic acid molecules comprising three or more nucleotides are linear molecules, in which adjacent nucleotides are linked to each other via a phosphodiester linkage. In some embodiments, “nucleic acid” refers to individual nucleic acid residues (e.g. nucleotides and/or nucleosides). In some embodiments, “nucleic acid” refers to an oligonucleotide chain comprising three or more individual nucleotide residues. As used herein, the terms “oligonucleotide” and “polynucleotide” can be used interchangeably to refer to a polymer of nucleotides (e.g., a string of at least three nucleotides). In some embodiments, “nucleic acid” encompasses RNA as well as single and/or double-stranded DNA. Nucleic acids may be naturally occurring, for example, in the context of a genome, a transcript, an mRNA, tRNA, rRNA, siRNA, snRNA, a plasmid, cosmid, chromosome, chromatid, or other naturally occurring nucleic acid molecule. On the other hand, a nucleic acid molecule may be a non-naturally occurring molecule, e.g., a recombinant DNA or RNA, an artificial chromosome, an engineered genome, or fragment thereof, or a synthetic DNA, RNA, DNA/RNA hybrid, or including non-naturally occurring nucleotides or nucleosides. Furthermore, the terms “nucleic acid,” “DNA,” “RNA,” and/or similar terms include nucleic acid analogs, e.g., analogs having other than a phosphodiester backbone. Nucleic acids can be purified from natural sources, produced using recombinant expression systems and optionally purified, chemically synthesized, etc. Where appropriate, e.g., in the case of chemically synthesized molecules, nucleic acids can comprise nucleoside analogs such as analogs having chemically modified bases or sugars, and backbone modifications. A nucleic acid sequence is presented in the 5′ to 3′ direction unless otherwise indicated. In some embodiments, a nucleic acid is or comprises natural nucleosides (e.g. adenosine, thymidine, guanosine, cytidine, uridine, deoxyadenosine, deoxythymidine, deoxyguanosine, and deoxycytidine); nucleoside analogs (e.g., 2-aminoadenosine, 2-thiothymidine, inosine, pyrrolo-pyrimidine, 3-methyl adenosine, 5-methylcytidine, 2-aminoadenosine, C5-bromouridine, C5-fluorouridine, C5-iodouridine, C5-propynyl-uridine, C5-propynyl-cytidine, C5-methylcytidine, 2-aminoadenosine, 7-deazaadenosine, 7-deazaguanosine, 8-oxoadenosine, 8-oxoguanosine, 0(6)-methylguanine, and 2-thiocytidine); chemically modified bases; biologically modified bases (e.g., methylated bases); intercalated bases; modified sugars (e.g., 2′-fluororibose, ribose, 2′-deoxyribose, arabinose, and hexose); and/or modified phosphate groups (e.g., phosphorothioates and 5′-N-phosphoramidite linkages).
Nucleic acid molecules assayed using the methods described herein need not be 100% identical with an endogenous nucleic acid sequence, but will typically exhibit substantial identity. Polynucleotides having “substantial identity” to an endogenous sequence are typically capable of hybridizing with at least one strand of a double-stranded nucleic acid molecule. Nucleic acid molecules useful in the methods of the disclosure include any nucleic acid molecule that encodes a polypeptide of the disclosure or a fragment thereof. Such nucleic acid molecules need not be 100% identical with an endogenous nucleic acid sequence, but will typically exhibit substantial identity. Polynucleotides having “substantial identity” to an endogenous sequence are typically capable of hybridizing with at least one strand of a double-stranded nucleic acid molecule. By “hybridize” is meant pair to form a double-stranded molecule between complementary polynucleotide sequences (e.g., a gene described herein), or portions thereof, under various conditions of stringency. (See, e.g., Wahl, G. M. and S. L. Berger (1987) Methods Enzymol. 152:399; Kimmel, A. R. (1987) Methods Enzymol. 152:507).
For example, stringent salt concentration will ordinarily be less than about 750 mM NaCl and 75 mM trisodium citrate, less than about 500 mM NaCl and 50 mM trisodium citrate, or about 250 mM NaCl and 25 mM trisodium citrate. Low stringency hybridization can be obtained in the absence of organic solvent, e.g., formamide, while high stringency hybridization can be obtained in the presence of at least about 35% formamide, and in some embodiments, at least about 50% formamide. Stringent temperature conditions will ordinarily include temperatures of at least about 30° C. at least about 37° C., or at least about 42° C. Varying additional parameters, such as hybridization time, the concentration of detergent, e.g., sodium dodecyl sulfate (SDS), and the inclusion or exclusion of carrier DNA, are well known to those skilled in the art. Various levels of stringency are accomplished by combining these various conditions as needed. In one embodiment, hybridization will occur at 30° C. in 750 mM NaCl, 75 mM trisodium citrate, and 1% SDS. In another embodiment, hybridization will occur at 37° C. in 500 mM NaCl, 50 mM trisodium citrate, 1% SDS, 35% formamide, and 100 μg/ml denatured salmon sperm DNA (ssDNA). In yet another embodiment, hybridization will occur at 42° C. in 250 mM NaCl, 25 mM trisodium citrate, 1% SDS, 50% formamide, and 200 μg/ml ssDNA. Useful variations on these conditions will be readily apparent to those skilled in the art.
For most applications, washing steps that follow hybridization will also vary in stringency. Wash stringency conditions can be defined by salt concentration and by temperature. As above, wash stringency can be increased by decreasing salt concentration or by increasing temperature. For example, stringent salt concentration for the wash steps will comprise less than about 30 mM NaCl and 3 mM trisodium citrate or less than about 15 mM NaCl and 1.5 mM trisodium citrate. Stringent temperature conditions for the wash steps will ordinarily include a temperature of at least about 25° C., at least about 42° C., or at least about 68° C. In some embodiments, wash steps will occur at 25° C. in 30 mM NaCl, 3 mM trisodium citrate, and 0.1% SDS. In other embodiments, wash steps will occur at 42° C. in 15 mM NaCl, 1.5 mM trisodium citrate, and 0.1% SDS. In other embodiments, wash steps will occur at 68° C. in 15 mM NaCl, 1.5 mM trisodium citrate, and 0.1% SDS. Additional variations on these conditions will be readily apparent to those skilled in the art. Hybridization techniques are well known to those skilled in the art and are described, for example, in Benton and Davis (Science 196:180, 1977); Grunstein and Hogness (Proc. Natl. Acad. Sci., USA 72:3961, 1975); Ausubel et al. (Current Protocols in Molecular Biology, Wiley Interscience, New York, 2001); Berger and Kimmel (Guide to Molecular Cloning Techniques, 1987, Academic Press, New York); and Sambrook et al., Molecular Cloning: A Laboratory Manual, Cold Spring Harbor Laboratory Press, New York.
As used herein, “obtaining” as in “obtaining an agent” includes synthesizing, purchasing, or otherwise acquiring the agent.
By “overlapping amplicons” is meant two or more amplicons that comprise a shared nucleic acid sequence but have at least one different terminal sequence.
“Polymerase” refers to an enzyme capable of catalyzing nucleic acid synthesis. A polymerase can be a DNA polymerase or an RNA polymerase. A polymerase can be characterized by its error rate, or the rate at which the polymerase inserts an incorrect nucleotide into the nucleic acid molecule it is synthesizing. In some embodiments, a polymerase can be a high-fidelity polymerase, which has a much lower error rate than a reference polymerase. A non-limiting example of a reference polymerase is Taq polymerase.
“Pooling,” as used herein, means combining multiple amplification reactions or groups of reactions. Pooling is synonymous with multiplexing.
By “portion” is meant a segment of an intact nucleic acid molecule. This portion contains, in some embodiments, at least 10%, 20%, 30%, 40%, 50%, 60%, 70%, 80%, or 90% of the entire length of the reference nucleic acid molecule. A portion may contain 10, 20, 30, 40, 50, 60, 70, 80, 90, or 100, 200, 300, 400, 500, 600, 700, 800, 900, or 1000 nucleotides.
The term “read,” “sequence read,” or “sequencing read” refers to sequencing data from a region of a nucleic acid molecule obtained from a single nucleic acid molecule. A read represents a short sequence of contiguous bases in the nucleic acid molecule and may be depicted, for example, as a chromatogram or as a linear string of letters that represent the nitrogenous bases of the nucleotide sequence, wherein A=adenine; G=guanine; C=cytosine; T=thymine; U=uracil; R=purine (A or G); Y=pyrimidine (C or T); N=any nucleotide; W=A or T; S=G or C; K=G or T; B=Not A; H=Not G; D=Not C; and V=Not T.
“Reduces” or “increases” refers to a negative or positive alteration, respectively, of at least 10%, 25%, 50%, 75%, or 100%.
By “reference” is meant a standard or control condition.
A “reference sequence” is a defined sequence used for sequence comparison. A reference sequence may be a subset of or the entirety of a specified sequence; for example, a segment of a full-length gene sequence, or the complete gene sequence. For nucleic acids, the length of the reference nucleic acid sequence will generally be at least about 50 nucleotides, at least about 60 nucleotides, at least about 75 nucleotides, about 100 nucleotides, or even about 300, 400, or 500 nucleotides or any integer thereabout or therebetween. In some embodiments, the length of the reference nucleic acid sequence will be less than 50 nucleotides. In some embodiments, the reference nucleic acid sequence will be more than 500 nucleotides.
The term “sequence variant,” as used herein, refers to an alteration in a sequence relative to a reference sequence. In one embodiment, a nucleotide sequence variant comprises one or more alterations relative to a reference nucleotide sequence. In some embodiments, the reference sequence is a consensus sequence. Optimally aligned sequencing reads obtained from multiple individuals of the same species or a population thereof, or multiple sequencing reads for the same individual, may be used to produce a consensus sequence. As contemplated herein, a “consensus sequence” refers to a nucleotide sequence that comprises the base most in common among all the sequencing reads at each nucleotide in the sequence.
In some embodiments, a sequence variant represents a variation relative to corresponding sequences in the same sample. In some embodiments, the sequence variant occurs with a low frequency (i.e., at least <1%) in the population (also referred to as a “rare variant”). For example, the sequence variant may occur with a frequency of about or less than about 5%, 4%, 3%, 2%, 1.5%, 1%, 0.75%, 0.5%, 0.25%, 0.1%, 0.075%, 0.05%, 0.04%, 0.03%, 0.02%, 0.01%, 0.005%, 0.001%, or lower. In some embodiments, the sequence variant occurs with a frequency above about 0.1%. In some embodiments, the sequence variant occurs at a frequency of above about 0.0025%.
By “somatic allele” is meant an allele specific to a non-germline cell (i.e., somatic cell).
By “somatic event” is meant the acquisition of a genetic variant by a somatic cell.
By “subject” is meant a mammal, including a human or a non-human mammal, such as a bovine, equine, canine, ovine, feline, or rodent (e.g., mouse, rat).
By “substantially identical” is meant a polypeptide or nucleic acid molecule exhibiting at least 50% identity to a reference amino acid sequence (for example, any one of the amino acid sequences described herein) or nucleic acid sequence (for example, any one of the nucleic acid sequences described herein). In some embodiments, such a sequence is at least 60%, 80% or 85%, 90%, 95% or even 99% identical at the amino acid level or nucleic acid to the sequence used for comparison.
Sequence identity is typically measured using sequence analysis software (for example, Sequence Analysis Software Package of the Genetics Computer Group, University of Wisconsin Biotechnology Center, 1710 University Avenue, Madison, Wis. 53705, BLAST, BESTFIT, GAP, or PILEUP/PRETTYBOX programs). Such software matches identical or similar sequences by assigning degrees of homology to various substitutions, deletions, and/or other modifications. In an exemplary approach to determining the degree of identity, a BLAST program may be used, with a probability score between e−3 and e−100 indicating a closely related sequence.
The term “tissue” refers to a group or layer of similarly specialized cells, which together perform certain special functions. The term “tissue-specific” refers to a source or defining characteristic of cells from a specific tissue.
By “unique molecular identifier (UMI)” is meant a distinct nucleic acid sequence that individualizes each primer used in an amplification reaction. For example, 500 primers having identical complementary nucleic acid sequences will have 500 different UMIs. UMIs facilitate the detection and removal of redundant sequencing reads.
Unless specifically stated or obvious from context, as used herein, the term “or” is understood to be inclusive. Unless specifically stated or obvious from context, as used herein, the terms “a,” “an,” and “the” are understood to be singular or plural.
Unless specifically stated or obvious from context, as used herein, the term “about” is understood as within a range of normal tolerance in the art, for example within 2 standard deviations of the mean. About can be understood as within 10%, 9%, 8%, 7%, 6%, 5%, 4%, 3%, 2%, 1%, 0.5%, 0.1%, 0.05%, or 0.01% of the stated value. Unless otherwise clear from context, all numerical values provided herein are modified by the term about.
Ranges provided herein are understood to be shorthand for all the values within the range. For example, a range of 1 to 50 is understood to include any number, combination of numbers, or sub-range from the group consisting 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24, 25, 26, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45, 46, 47, 48, 49, or 50.
The recitation of a listing of chemical groups in any definition of a variable herein includes definitions of that variable as any single group or combination of listed groups. The recitation of an embodiment for a variable or aspect herein includes that embodiment as any single embodiment or in combination with any other embodiments or portions thereof.
Any compositions or methods provided herein can be combined with one or more of any of the other compositions and methods provided herein.
The present disclosure features methods for detecting and quantifying genetic variants in a sample.
The invention is based, at least in part, on the discovery of triple primer PCR sequencing (“TriPP-seq”), which provides a highly sensitive, low-cost approach for detecting and validating mutation on a highly scalable system. Mosaic mutations in somatic or germline cells contribute to a wide range of human disorders. As such, their identification and accurate allelic fraction quantification from tissue-derived and cell-free DNA are essential for clinical diagnoses and early detection of cancers. However, rapid, low-cost detection and validation of ultra-low alternate allelic fraction (AAF) mutations has traditionally required expensive and low throughput methods that have limited widespread testing. Recent methods, (e.g., ddPCR) have shown great promise for detection and validating known mutations at very low AAFs, but remain low-throughput due to allele-specific optimization.
Accordingly, the present disclosure features methods for detecting low frequency genetic variation. The present disclosure's novel approach is based on generating deep coverage of overlapping amplicons of a target nucleic acid sequence. Because the primers used in the reactions are designed to allow discernment and segregation of the overlapping amplicons, the sequencing data can be segregated into groups, and analysis of the sequencing data can be performed in parallel. The methods provide not only deep coverage of the target nucleic acid, but also a cost-effective means of characterizing and validating sequencing results.
Recently, the important roles of somatic mutations beyond cancer are becoming more appreciated with discoveries of somatic mutations across a wide range of neurodevelopmental, overgrowth, and hematological disorders. Even more, the presence of somatic mutations in healthy cells and individuals are associated with normal development and aging and are, therefore, a powerful tool for understanding how cells divide and form complex organs like the human brain. Finally, with the detection of cell-free DNA (e.g., fetal and tumor), it is becoming possible for early detection of disease, tracking of disease recurrence in cancers, and even non-invasive prenatal genetic testing where mutations of the placenta are detected in the pregnant mother's blood sample. The rapid advancements in sequencing technologies and interest in genetic mutation present at low alternate allelic fraction (i.e., ratio of DNA fragments carrying the mutation to those with the wild-type allele in a given samples; AAF) poses some major challenges for both the clinical and research communities related to the sensitivity to detect mutations, false positives, and the precision of the assessed AAFs. These challenges are often confounded by the inability to directly assess tissues with the highest AAFs, as is the case with brain tissue, or by limited or degraded DNA samples, as is typical for cell free DNA.
While germline mutations are relatively easy to detect with small amounts of DNA with variable qualities using WES, WGS, targeted gene panels, and traditional Sanger sequencing due to the equal fractions of mutant to wild-type alleles (50% AAF) in a given DNA sample, the AAF of a somatic mutation will depend on the given tissue, cell type, and the stage in development at which the mutation arose. Traditional WGS and WES sequencing in both the research and clinical diagnostic settings are optimized to identify germline events, but often lack the sequencing depth to robustly detect low-AAF variants. However, many recently improvements allow for robust detection of mutations present at greater than 0.1% AAF. These tools often employ strategies such as molecular barcoding, increased read depth, and reduced use of PCR to mitigate sequencing-induced errors while improving sensitivity. Despite these measures, the identification of somatic alleles, particularly those at very low AAFs has an elevated false positive rate compared to germline mutations. Therefore, while essential, the validation of large numbers of somatic alleles is often challenging due to many factors like assay costs, throughput, and sensitivity limitations.
The methodology utilized to accurately detect or validate somatic mutations have rapidly advanced in the last few years. The challenge for validating or measuring low AAFs is multifaceted, spanning sequencing platforms, inherent error rates of polymerases, and locus specific challenges. Each of these result in additional errors and skewing of AAFs, which can mask or alter the detected AAF in each assay. The utilization of PCR to amplify the genomic loci without inducing additional mutations and maintain the original AAFs has been improved using improved polymerases with proofreading capabilities and, in some cases, unique molecular barcodes for each DNA fragment. Additionally, errors can occur during sequencing on both the Illumina and Ion Torrent platforms. For example, in one study, the Ion Torrent had an error rate ˜0.05% for SNVs but ˜1.5% for indels while the on the Illumina MiSeq had 0.1% errors for SNVs and 0.7% for indels.
The original methods used employed either pyrosequencing or bacterial cloning followed by sanger sequencing of hundreds or thousands of individual bacterial colonies to measure a single mutation. These methods, while accurate and robust, were often cost-prohibitive, less scalable to large numbers of mutations, and were less sensitive for mutations below 5% AAF. These methods were recently succeeded by the advancement of digital droplet PCR, ddPCR, where an allele-specific PCR conditions are designed to allow for the measurement of mutation positive and negative DNA fragments in thousands of droplets. This method is routinely considered a gold standard for validation of somatic alleles in both research and clinical settings, but each allele requires the development of a custom assay, validation and optimization prior to use. The ddPCR assay can accurately detect AAFs below 0.5%, but its sensitivity relies on the quantity and concentration of input DNA and the number of positive droplets formed in each reaction. Despite its great success, the use of ddPCR is somewhat limited as it remains limited by scalability, the potential for allelic dropout, and the ability to design allele-specific primers, which is more challenging in repetitive regions and for small indels.
The growing consensus that somatic mutations might underly a wide range of clinical phenotypes ranging from cancer risk to severe neurodevelopmental and overgrowth conditions suggests that a robust method for both detection and validation of alleles and their mosaic fraction in the body is essential. Here, an improved strategy that aims to mitigate the previously stated limitations for assessing somatic mutations is presented. This strategy, which can be referred to as triple-primer PCR, relies on the power of designing and running at least 3 unique, nonoverlapping amplicons over a suspected mutation. Through independently analyzing each amplicon, the impact of allelic dropout, amplification bias, sequencing and PCR induced artifacts, and general optimization challenges, are markedly reduced while achieving the highest sensitivity to accurately detect ultra-low allelic fractions below 0.1% regardless of tissue origin. As described, below, this triple-primer PCR sequencing method allows for additional improvements to future improve accuracy through incorporations of molecular barcoding and improved purification processes.
Nucleic acid amplification according to the presently disclosed methods requires at least two pairs of primers and in some embodiments, at least three pairs of primers. Each pair of primers comprises a forward and a reverse primer, and each primer comprises a complementary nucleic acid sequence that is at least 85% complementary to a nucleic acid sequence (i.e., the primer binding site) on a template nucleic acid molecule. The primers of each pair define the termini of an amplicon that is generated by an amplification reaction, and the region of the amplicon between the termini comprises the target nucleic acid sequence. The combined length of the primers and the target sequence is referred to as the amplicon length. Amplicon length is typically between about 150 and about 500 nucleotides. In some embodiments, the length of the amplicon is about 150, 200, 250, 300, 350, 400, 450, 500, or any integer in-between, nucleotides. In some embodiments, the length of the amplicon is less than 150 nucleotides. In some embodiments, the length of the amplicon is greater than 500 nucleotides. Each primer has a unique nucleic acid sequence that can bind to a complementary primer binding site on the template nucleic acid.
Amplicons generated by amplification reactions using one of the primer pairs will be distinguishable from other amplicons generated by amplification reactions that use different primer pairs due to the length and sequence of the amplicon (
A primer binding site in a template nucleic acid sequence may harbor a variant that impairs primer biding, which results in decreased amplification of the template harboring the variant and a loss of sequencing coverage of the allele. The resulting loss of coverage of a particular variant is allelic dropout. Referring to
In some embodiments, the complementary nucleic acid sequence of a primer is about 15, 16, 17, 18, 19, 20, 25, 30, 35, or even 40 nucleotides long. In some embodiments, the complementary nucleic acid sequence of a primer is between about 85% and about 100% complementary to a nucleic acid sequence in the template nucleic acid molecule. In some embodiments, the complementary nucleic acid sequence of the primer is between about 85%, 86%, 87%, 88%, 89%, 90%, 91%, 92%, 93%, 94%, 95%, 96%, 97%, 98%, 99%, and 100% complementary to a nucleic acid sequence in the template nucleic acid molecule. In some embodiments, wherein the complementary nucleic acid sequence of the primer is less than 100% complementary with a primer binding site in the template nucleic acid molecule, the mismatch nucleotide or nucleotides in the primer reside at least three bases from the 3′ terminus of the primer. This allows for efficient binding at the terminus of the primer to the template molecule, which facilitates polymerase binding to the primer:template hybrid and extending the primer.
In some embodiments, a primer is comprised of DNA or RNA nucleotides. In some embodiments, a primer comprises at least one modified base. A modified base includes, but is not limited to, those nucleotide analogs described herein or a labeled nucleotide. In some embodiments, a primer may have a modified backbone comprising at least one phosphorothioate linkage. In some embodiments, the primer comprises a label, such as, but not limited to, a fluorescent label, a radiolabel, a nanoparticle label, and/or a biotin label.
In some embodiments, each primer will have an adapter upstream from the complementary nucleic acid sequence. The adapter has a nucleic acid sequence that is complementary to a sequence of a nucleic acid molecule used in a downstream sequencing reaction. For example, the adapters used in some embodiments are designed to be compatible with Next Generation Sequencing including, but not limited to, Ion Torrent and MiSeq platforms. In some embodiments, the length of the adapter is between 8 and 20 nucleotides. In some embodiments, the length of the adapter is 8, 9, 10, 11, 12, 13, 14, 15, 16, 17, 18, 19, or 20 nucleotides. The adapter's sequence is designed to reduce or eliminate nonspecific binding of the adapter to the template nucleic acid molecule. In some embodiments, the adapter is designed to have a sequence that is not substantially complementary to any nucleic acid sequence present in the template nucleic acid molecule. In some embodiments, the adapter is designed to diverge from perfect complementarity with the template by 2, 3, or 4 or more nucleotides.
At least one primer in each pair also has an index sequence, or barcode (
In some embodiments, at least one primer in each pair comprises a unique molecular identifier (UMI) (
There are approximately 1,000 possible sequences for a 5-nucleotide UMI, approximately 65,000 possible sequences for an 8-nucleotide UMI, approximately 1×106 possibilities for a 10-nucleotide UMI, and approximately 1×1012 possibilities for a 20-nucleotide UMI. Even if some UMIs are not suitable for the reasons given above, large UMI libraries can be produced for use in the presently disclosed methods. Use of nucleotide analogs increases the number of possible sequences for a UMI.
Table 1 characterizes five primer pairs used in the disclosed methods. In this table, “Chr. No.” means chromosome number; “Ref” refers to the reference nucleotide; and “Alt” refers to the alternate nucleotide. Each of the primer pairs is designed to amplify a region containing a single nucleotide variant (the “allele start” and “allele end” are the same locus number). Three of the primer pairs on Table 1 (X:153579431-153579431/T/C-F1; X:153579431-153579431/T/C-F2; and X:153579431-153579431/T/C-F3) are used to interrogate a single nucleotide variant in the Filamin A (FLNA) gene on the X chromosome. The remaining two primer pairs (X:153579431-153579431/T/C-F1 and X:153579431-153579431/T/C-F2) are used to interrogate a single nucleotide variant in the SR-Related CTD Associated Factor 11 (SCAF-11) gene on chromosome 12. The amplicons generated in amplification reactions comprising the primer pairs disclosed in Table 1 will be about 220 to 260 nucleotides in length.
Samples comprising template nucleic acid molecules to be assayed using the methods disclosed herein can be obtained from a variety of sources including, but not limited to, tissue biopsies, blood draws, buccal swabs, hair, sweat, skin, semen, and mucus. In some embodiments, the sample comprises cells from a subject, for example, circulating tumor cells, blood cells, skin cells, and the like. In some embodiments, the sample comprises cell free nucleic acid, such as, but not limited to, cell free tumor nucleic acid and cell free fetal nucleic acid. In some embodiments, the template nucleic acid molecule is isolated or purified before amplification. Methods of isolating and purifying nucleic acids are well known in the art. Template nucleic acid molecules comprise at least one target nucleic acid sequence. The target sequence is flanked by primer binding sites. In some embodiments, the template is a DNA molecule. In some embodiments, the template is an RNA molecule. In some embodiments, the template may be double-stranded, while in other embodiments, the template is single-stranded.
In some embodiments, the target nucleic acid is a portion of a gene such as, but not limited to, ABCC8, ABLIM3, ACBD3, ACIN1, ACSL5, ACTA2, ACVR1, ACVR1B, ACVR1C, ACVR2B, ADAMTSL3, ADORA2A, AEBP2, AES, AFAP1, AGAP1, AKR7A2, AKT1, ALK, AMHR2, AMPD3, ANGPTL6, ANO7, APC, APOL2, AQP4-AS1, ARHGEF3, ARID1A, ARIDSA, ARIH1, ARNT, ATM, ATP5A1, ATP9B, ATXN7L1, AX747372, BAG1, BAIAP2L1, BECN2, BMP4, BMP8A, BMP8B, BMPR1A, BMPR1B, C12orf60, C17orf89, C1ORF210, C6ORF10, C6orf211, C9orf40, CACNA1A, CACNA1H, CACNA2D4, CAMK1D, CAMKMT, CARM1, CAST, CBS, CCBE1, CDC40, CDH23, CDH4, CDKN2B, CHRNA4, CLASP1, CLCA1, CLDN2, CLIC3, CNN3, CNTN1, COL11A2, COL3A1, COL3A2, COL4A1, COL4A5, COL4A6, COL5A1, COL5A2, COL6A2, COL6A3, COX7A2L, CRADD, CREBBP, CRY2, CSGALNACT2, CTBP2, CYP2S1, DAG1, DCAF8, DCAF8,DCAF8, DLAT, DLGS, DLGAP4-AS1, DNAH3, DOCK4, DOCK8, DOPEY1, DPYSLS, DYNC1H1, DYNC1I2, DYRK2, E2F4, E2F6, ECI2, EEF1DP3, EHD4, EIF2B5, EIF4G3, ELAC2, ELK3, EMD, EMX20S, EPPK1, EPT1, ERBB4, ERCCS, ETS2, ETV4, FAM107B, FAM13B, FAM175A, FAM83E, FAV, FBN1, FBN2, FBN3, FBXO28, FGFR2, FHL2, FIRRE, FLNA, FLT3, FOXA3, FOXG1-AS1, FST, GABRG1, GALM, GAPDH, GDF6, GDF7, GLI2, GLI3, GLRXS, GLT8D2, GOLPH3, GPD2, GPR68, GPRASP1, H2AFX, HDAC4, HHAT, HIST1H2AH, HIST2H2AB, HK1, HMCN1, HMSD, HNF4A, HNRNPU, HOXD3, HPS3, HS3ST3A1, IDH1, IFNG, IKBKAP, IMP3, INHBA, INPP4B, INPP5A, IQCK, JAG1, JWT213-1, JWT213-2, JWT213-3, JWT213-4, JWT213-5, JWT213-6, JWT213-7, JWT213-8, JWT213-9, JWT307_1, JWT307_2, JWT307_3, JWT307_4, JWT307_5, JWT307_6, JWT307_7, JWT310-1, JWT310-2, JWT310-3, JWT310-4, JWT310-5, JWT310-6, JWT310-7, JWT311-1, JWT311-2, JWT311-3, JWT311-4, JWT311-5, JWT311-6, JWT311-7, JWT312-1, JWT312-2, JWT312-3, JWT312-4, JWT312-5, JWT312-6, JWT312-7, JWT312-8, JWT312-9, JWT313-1, JWT313-2, JWT313-3, JWT313-4, JWT313-5, JWT313-6, JWT313-7, JWT313-8, JWT313-9, JWT364_1, JWT364_2, JWT364_3, JWT364_4, JWT364_5, JWT364_6, JWT364_7, KANSL1, KCNQ1, KDM3A, KDR, KIRREL3, KLF13, KLHL14, KMTD2, L3MBTL1, LACTB2, LAMA2, LAMA3, LEFTY1, LINGO4, LMAN2L, LRRC4C, LSAMP, LTBP1, LTBP2, LTBP3, LZTS2, MAD1L1, MAD2L1, MAEA, MAGI2, MAML2, MAP3K7, MAPK1, MAPK3, MAPK8IP2, MARK3, MAT2A, MATR3, MBNL2, MCL1, MCU, MECP2, MED12, MED29, MEF2A, MEGF6, MESD, METTL17, MIER2, MIR181A1HG, MKL1, MKL2, MLH1, MOB2, MPRIP, MRPL32, MRS2, MTCH1, MTOR, MUC16, MUC3A, MYC, MYH11, MYH11,NDE1, MYH11; MYH11, MYLK, MYLK-AS1, MYOCD, NA, NDFIP2, NDUFC1, NEK9, NF1, NFKB1, NGEF, NME4, NME4,DECR2, NOL9, NOTCH1, NOTCH3, NPLOC4, NRG4, NRM, NRTN, NTM, NUCB1, NUDT16, NUDT16L1, OAS3, OR4K3, OSTC, PAG1, PCDH15, PDCD6, PDE4DIP, PDSSA, PHC1, PHF12, PHKG1, PIK3R1, PLEKHG6, PLXDC2, PMM2, POLG2, POLR3B, PPARGC1A, PPHLN1, PPP1R14A, PPP1R15B, PRAF2, PRDM16, PRKG1, PRPH2, PRTG, PTGDR, PTPN12, PTPN14, PTPRC, PTPRS, PUS7, RABL6, RALGAPA1, RAPGEF4, RBM10, REPS2, RHBDF2, RIN2, RNF175, RNU1-35P, RNU1-35P, RP11-149P24.1, ROCK1, ROCK2, RPRD2, RSF1, RUSC1, SAFB2, SASH1, SCAF11, SCARF1, SEPT11, SH3GLB2, SHPK, SHPK, SHPK, SHROOM3, SIKE1, SIPA1L1, SIRPA, SK213, SK215, SLAIN1, SLC1A4, SLC25A48, SLC2A10, SLC4A1AP, SLMO2, SLTM, SLX4, SMAD3, SMAD4, SMAD5, SMAD6, SMAD7, SMARCA4, SMLR1, SMTNL1, SMURF1, SNK307, SNK310, SNK311, SNK312, SNK313, SNK364, SNK380, SNK382, SNK383, SNK384, SNK385, SNK386, SOX21-AS1, SOX9, SPOCK2, SPRED1, SPSB2, SRGN, SRP68, SRRM2-AS1, ST6GAL1, STK16, STRN3, SUCLA2, SUCO, SWI5, SYNE2, TAB1, TBC1D13, TBCE, TCERG1, TCF4, TERT, TFB2M, TFDP1, TGFB1, TGFB3, TGFBR1, TGFBR2, THBS1, TMEFF2, TMEM132C, TMEM2, TMEM268, TNPO1, TPCN2, TPM3, TPRX1, TRAM1, TRAPPC9, TRPM1, TSC2, TSHZ2, TTN, TUBG1, TUBGCP3, TULP4, UBAP2, UBE2I, UBE2W, UHRF1, UNC45A, UNG, UROC1, USP24, USP34, USP8, VANGL1, VIPR2, VPS13D, WDR35, WDR45B, WDR77, WDSUB1, WHSC1, YARS2, YIPF3, ZFHX4, ZFYVE16, ZFYVE9, ZMIZ1, ZNF223, ZNF292, ZNF3, ZNF362, ZNF451, ZNF517, ZNF593, ZNF630, ZNRF3, or ZSCAN5A.
The subject from whom the template nucleic acid molecule sample is obtained can be any organism. In some embodiments, the subject is a vertebrate. In some embodiments, the subject is a mammal such as a human, mouse, rat, dog, cat, horse, cow, sheep, or other domesticated mammal. In some embodiments, the mammal is a human. In some embodiments, the subject from whom the sample is obtained has or is suspected of having a disease or condition associated at least in part with a genetic variant or variants.
The methods provided herein use a nucleic acid polymerase to amplify a target nucleic acid sequence. Because some polymerases have high error rates (incorporating the wrong nucleotide at a position in a synthesized nucleic acid), selection of a suitable polymerase is an important concern. Sequence errors introduced by a polymerase confound authentic sequence data, making discernment of low frequency variants unreliable or expensive due to the amount of coverage necessary to overcome the polymerase's error rate. High-fidelity polymerases, are particularly well-suited for use in the presently disclosed methods, and can be used to synthesize copies of a target nucleic acid sequence that potentially harbors a low-frequency variant. Such high-fidelity polymerases introduce fewer nucleotide sequence errors than non-high-fidelity polymerases. Thus, in some embodiments, the nucleic acid amplification reactions comprise a high-fidelity nucleic acid polymerase. For example, in some embodiments, nucleic acid reactions comprise a Phusion high-fidelity DNA polymerase (New England Biolabs (NEB)). This polymerase has a reported error rate of 4.4×10−7 errors per base in Phusion HF buffer and 9.5×10−7 errors per base in GC buffer. Thermus aquaticus (Taq) polymerase has a 50-fold higher error rate than the error rate of the Phusion high-fidelity polymerase. Other polymerases may be used to amplify nucleic acids according to the presently disclosed methods, but an increase in polymerase error rates may decrease the reliability of the method. Table 2 provides a summary of the differences between the high-fidelity Phusion DNA polymerase and the Pyrococcus furiosus and the Taq DNA polymerases (HF=high-fidelity; “GC Buffer” refers to a buffer suited for reactions amplifying a target rich in G and/or C). To overcome such errors generated by non-high-fidelity polymerases, additional coverage of the interrogated nucleic acid may be necessary, resulting in increased costs.
Pyrococcus furiosus DNA polymerase
The methods disclosed herein are suitable for detecting low frequency variants. The methods described herein involve detecting the presence or absence of low frequency genetic variation in a nucleic acid molecule by amplifying the nucleic acid sequence of interest using multiple pairs of primers. Each pair of primers comprises a forward primer and a reverse primer, each having a unique binding sequence complementary to a target polynucleotide, wherein the intervening sequences between each pair of primers (i.e., the amplified nucleic acid sequence) at least partially overlap. The resulting overlapping amplicons are sequenced using a Next Generation Sequencing platform, which provides the deep coverage necessary to validate low frequency variants. The sequencing reads are aligned, and determinations regarding the presence or absence of genetic variation are made. The sequencing data can be used for further characterization of any detected genetic variation (i.e., alternative allele fraction).
In some embodiments, the low frequency variant is a known variant, and the methods disclosed herein may be used to confirm the variant's presence and/or characteristics (i.e., its alternate allele frequency). In some embodiments, the low frequency variant originated during a germline event, while in other embodiments, the low frequency variant to be interrogated originated during a somatic event. In some embodiments, the low frequency variant is a silent variant, a missense variant, or a nonsense variant. In some embodiments, the low frequency variant alters a splice site or is an insertion or deletion.
In some embodiments, nucleic acid amplification reactions comprise a template nucleic acid molecule having a target nucleic acid sequence, at least three primer pairs suitable for interrogating the target nucleic acid, nucleotides, and a polymerase. Due to the use of at least three primer pairs in the amplification, the overall method described herein can be referred to a triple-primer PCR sequencing. In some embodiments of the present disclosure, the reaction further comprises a buffer that provides a suitable ionic environment for the polymerase to synthesize a nucleic acid molecule. In some embodiments, the reaction comprises a buffer having essential cofactors (e.g., magnesium) necessary for polymerase function. In some embodiments, the cofactors necessary for proper polymerase function are added to the reaction independently of the buffer.
In some embodiments, the amplification reaction comprises labeled nucleotides, wherein the labeled nucleotides facilitate efficient capture of any amplicon that comprises one or more labeled nucleotides. Referring to
In some embodiments, separate nucleic acid amplification reactions are prepared for each pair of primers. For example, amplifying a target nucleic acid sequence may comprise at least three reactions according to the methods described herein, wherein each reaction comprises one of three different pairs of primers. The primers, as discussed supra, are used in amplification reactions that generate overlapping amplicons (i.e., semi-redundant interrogation of the target nucleic acid sequence), thereby reducing the probability of impaired detection of variants or skewed downstream determination of alternate allele fractions due to amplification bias. In some embodiments, a single amplification reaction will comprise all pairs of primers. Combining the different primers into a single amplification reaction will generate a greater number of distinct amplicons.
In some embodiments, the amplification reactions are polymerase chain reactions (PCR). PCR reactions undergo multiple thermocycles, wherein each thermocycle comprises a denaturing step, an annealing step, and an extension step. During the denaturation step, the reaction is incubated at or above 90° C., which is a sufficient temperature, in some embodiments, to cause a double-stranded DNA molecule to denature into single DNA strands or to cause the nucleic acid molecule to undergo a conformational change that is more conducive for an amplification reaction.
The annealing step comprises complementary binding of the primers to the template nucleic acid and occurs at a lower temperature than that used in the denaturing step. In some embodiments, each primer will be designed to anneal to a complementary nucleic acid sequence at a temperature of between about 50° C. and about 65° C. In some embodiments, the annealing temperature is about 50° C., 51° C., 52° C., 53° C., 54° C., 55° C., 56° C., 57° C., 58° C., 59° C., 60° C., 61° C., 62° C., 63° C., 64° C., or 65° C. about 50° C., 51° C., 52° C., 53° C., 54° C., 55° C., 56° C., 57° C., 58° C., 59° C., 60° C., 61° C., 62° C., 63° C., 64° C., or 65° C. In some embodiments, the temperature at which the primers anneal to the nucleic acid template can be modified by adjusting conditions (e.g., salt concentration) in the sample or in the amplification reaction. One skilled in the art will understand how changing sample or reaction conditions can affect the temperature at which a primer binds to template nucleic acid.
In the extension step of a PCR cycle, the primers annealed to the template nucleic acid's primer binding sites are extended by a polymerase to produce a nucleic acid molecule that is complementary to a portion of the template nucleic acid molecule. A proper extension temperature is at or about the optimal temperature for the polymerase to synthesize a nucleic acid molecule. In some embodiments, the extension temperature is between about 65° C. and 75° C. In some embodiments, the extension temperature is about 65° C., 66° C., 67° C., 68° C., 69° C., 70° C., 71° C., 72° C., 73° C., 74° C., or 75° C. In some embodiments, the extension temperature may be 5, 10, 15, 20, or 25% higher or lower than the optimal temperature of the polymerase. Those skilled in the art will understand how to adjust the temperatures, or other reaction conditions, necessary for successful PCR amplification of a nucleic acid sequence.
In some embodiments, the template nucleic acid is amplified isothermally. For example, helicase dependent amplification is an isothermal amplification method that utilizes a helicase, rather than high temperatures, to separate the strands of a duplex nucleic acid. By not requiring a denaturation step, the isothermal reaction can be incubated at or about the optimal temperature of the polymerase. However, in some embodiments, the isothermal amplification reaction comprises an initial heat denaturation step. Exponential amplification is achieved by incubating the reaction at an isothermal temperature, which obviates the need for thermocycling equipment. Other isothermal amplification techniques are known in the art, and one skilled in the art would understand how to optimize these techniques to comport with the methods described herein.
Referring to
In some embodiments, the amplification reaction products are purified or isolated before pooling. Methods for isolating and purifying nucleic acids are well known in the art, and there are many commercially available kits for purifying or isolating amplicons. In some embodiments, purifying or isolating amplicons occurs after pooling. In some embodiments, enriched amplicons resulting from biotin:streptavidin capture and reamplification, can be purified using streptavidin to bind and separate all biotin labeled amplicons.
In some embodiments, the amplicons are assessed prior to being sequenced. Assessing the amplicons can include, for example, gel electrophoresis, real time detection, or spectrophotometric determination of amplicon concentration. For example, amplicons may be assessed using a TapeStation (Agilent) or Bioanalyzer 2100 (Agilent). These analyses allow an investigator to determine if the amplification reaction generated sufficient amounts of high quality amplicons for subsequent sequencing.
Sequencing of the overlapping amplicons provides multiple independent interrogations of a variant nucleotide or nucleic acid sequence compared to using a single pair of primers. Traditional Sanger sequencing platforms can be used to sequence the overlapping amplicons, but this approach is inefficient for detecting rare variants. Conversely, Next Generation Sequencing (NGS) platforms can generally accommodate thousands of sequencing reactions run in parallel, thereby providing deeper coverage than is possible with Sanger sequencing. For example, referring to
The amplicons to be sequenced are, by design, generally less than 300 nucleotides in length, and there are several NGS platforms that can cost-effectively generate sequencing data at the desired coverage level. For example, ThermoFisher's Ion Torrent and Illumina's MiSeq can each generate maximum read lengths of approximately 250 nucleotides. Other NGS approaches are available for shorter or longer read lengths. For example, Illumina's HiSeq platform has a maximum read length of about 150 nucleotides, while the Roche 454 platform can generate at least 400 nucleotide reads. One skilled in the art will be to determine which platform can be used to generate the desired sequencing data, and will optimize the adapters on each primer to comport with that platform.
In some embodiments, the sequencing data is assessed for quality before alignment, and those reads not possessing the required quality characteristics are removed from the data set. Typically, quality control of sequencing reactions comprises establishing a signal-to-noise threshold, and reads that do not meet the threshold are discarded. Such quality control lessens the probability of erroneous base calls in a read that would decrease reliability of the assay.
Sequencing data generated using the disclosed methods can be processed to accurately determine alternate allele frequencies. Referring to
The data in each bin is aligned to provide maximal sequence identity between the individual reads. For example, if a read has a single nucleotide deletion, the alignment will incorporate the deletion into the read's aligned sequence so that the nucleotide sequences on either side of the deletion align with other reads that do not have the deletion. Referring to
Primer binding sites are also identified (
In some embodiments, all but one read having the same unique molecular identifiers will be removed from the data set, which indicates multiple amplification reactions that used the exact same primer. These duplicated amplifications reactions are not considered independent interrogations of the nucleic acid. Retention of such redundant data could impact alternate allele fraction determination. In some embodiments, accurate determination or validation of alternate allele frequencies of about 0.025% comprise removing redundant reads from the data. In some embodiments, wherein the alternate allele fraction is known to be 0.1% or greater, removal of redundant reads may not be necessary due to the deep coverage available in Next Generation Sequencing platforms. Once the alignment is set in each bin, the alternate allele frequencies for variants in each bin are determined.
The methods provided can distinguish between germline and somatic events resulting in genetic variation. Referring to
A somatic event occurring in a single subject will likely have a much lower allele frequency than an inherited allele, and a subject having a genetic variant derived from a somatic event is said to be mosaic for the variant. As shown in Table 3, the alternate allele frequencies (AAF) observed in three different amplicon samples are about 1%, well below the frequency expected in an individual for an inherited allele, which suggests the variant is a somatic mosaic variant. For example, for the sequencing reads of amplicons generated using the Primer 1 set of primers, 416 reads out of 37,779 total reads contained the alternate allele (
Two methods are currently used to detect and quantify rare variants, droplet digital PCR (ddPCR) and Sanger sequencing of TOPO (Topoisomerase-based) cloned nucleic acids. Referring to Table 4, the estimated cost of the method described herein (“mosaic validation method”) is about 90% less expensive than ddPCR and 85× less expensive than the Sanger sequencing/TOPO cloning method. Furthermore, the Sanger sequencing/TOPO cloning method is much less sensitive as its lowest level of reliable detection is an alternate allele fraction of 0.5%. While the purported resolution of ddPCR is an alternate allele fraction of 0.1%, it is not reliable for alternate allele fractions of 0.02% that are within the reliable range of the presently disclosed methods.
Additionally, high-throughput Next Generation Sequencing platforms used in the presently disclosed methods can run massive parallel reactions. Conversely, both Sanger Sequencing/TOPO cloning and ddPCR have relatively limited throughput, thereby increasing cost and time requirements. ddPCR, while having higher throughput than the Sanger sequencing/TOPO cloning method, does not enjoy the throughput of the presently described methods. Additionally, ddPCR primers are labeled with a relatively expensive fluorophore.
The methods described herein can be used for the detection and/or monitoring of a disease. The detection and characterization of disease-associated variants, including somatic mosaic variants, can provide information relevant for diagnosing a disease, determining the progression or regression of disease, and treating disease. For example, when a cancer cell arises after a somatic event, or when circulating tumor cells are present in a subject, the methods described herein can be used to detect of these cells.
A subject having a disease may undergo periodic testing to determine if the number of a diseased cells is increasing, decreasing, or static. For example, a subject that has cancer may determine the alternative allele frequency of a cancer marker present in samples after the cancer is detected or after treatment has begun. Changes in the alternative allele frequency of the cancer marker would indicate a change in the number of cells carrying the marker (e.g., cancer cells) present in the sample. If the alternative allele frequency is greater than that observed in a previous sample, the subject's cancer is likely progressing or not responding effectively to treatment. If the alternative allele frequency remains static relative to an earlier sample, the disease may be responding treatment sufficiently to stop disease progression, but perhaps not to a level sufficient for disease regression or remission. If the alternative allele frequency decreases relative to an earlier sample, the subject's disease may be regressing, and the absence of such cells (i.e., AAF=0) may signify remission.
In another embodiment, kits and compositions are provided that advantageously allow for the detection and/or quantification of the presence of low frequency genetic variation in a subject sample (e.g., blood or serum). In one embodiment, the kit includes a composition comprising reagents for performing an amplification reaction, including multiple pairs of forward and reverse primers as described herein. In some embodiments, the reagents include nucleotides, labeled nucleotides, a buffer, a cofactor, and/or a polymerase. In some embodiments, the kit comprises a sterile container that contains the amplification reaction reagents; such containers can be boxes, ampoules, bottles, vials, tubes, bags, pouches, blister-packs, or other suitable container forms known in the art. Such containers can be made of plastic, glass, laminated paper, metal foil, or other materials suitable for holding amplification reagents.
In one embodiment, the kit comprises high-quality (PAGE-purified) RNA or DNA-based primers, premixed at proper concentrations. In some embodiments, the kit comprises reagents for biotin labeling for higher sensitivity assays. In some embodiments, the kit comprises a preselected polymerase (e.g., Phusion U if using RNA primers, or another option) with high fidelity (100× improved error rates compared to a reference polymerase (Taq polymerase). In some embodiments, the kit comprises duplicate primers with differing barcodes for testing case/control samples side-by-side. In some embodiments, the kit comprises preselected primers to avoid other mutation sites, non-overlapping binding sites, and the like. In some embodiments, the kit comprises control DNA (e.g., for negative controls). In some embodiments, the kit comprises ddPCR probes for performing ddPCR and sequencing from the same reaction—(i.e., to obtain copy/expression values and genotype correlation).
In another embodiment, the kit includes a composition comprising reagents for performing a sequencing reaction, including nucleic acid molecules that can specifically bind to an adapter as described above. The reagents, in some embodiments, include nucleotides, labeled nucleotides, a buffer, a cofactor, ion spheres comprising the nucleic acid molecule to be sequenced, and/or enzymes for catalyzing the sequencing reaction. In some embodiments, the kit comprises a sterile container that contains the sequencing reaction reagents; such containers are described above.
In some embodiments, the kit comprises compositions for amplification and sequencing as described above. Kits may also include instructions for performing the reactions.
The practice of the present disclosure teaches, unless otherwise indicated, conventional techniques of molecular biology (including recombinant techniques), microbiology, cell biology, biochemistry and immunology, which are well within the purview of the skilled artisan. Such techniques are explained fully in the literature, such as, “Molecular Cloning: A Laboratory Manual”, second edition (Sambrook, 1989); “Oligonucleotide Synthesis” (Gait, 1984); “Animal Cell Culture” (Freshney, 1987); “Methods in Enzymology” “Handbook of Experimental Immunology” (Weir, 1996); “Gene Transfer Vectors for Mammalian Cells” (Miller and Calos, 1987); “Current Protocols in Molecular Biology” (Ausubel, 1987); “PCR: The Polymerase Chain Reaction”, (Mullis, 1994); “Current Protocols in Immunology” (Coligan, 1991). These techniques are applicable to the production of the polynucleotides and polypeptides of the disclosure, and, as such, may be considered in making and practicing the compositions and methods disclosed herein. Particularly useful techniques for particular embodiments will be discussed in the sections that follow.
The following examples are put forth to provide those of ordinary skill in the art with a complete disclosure and description of how to perform the amplification, sequencing, and quantifying methods presently disclosed, and are not intended to limit the scope of what the inventors regard as their invention.
To identify low frequency genetic variation in a target nucleic acid sequence with an alternate allele fraction (AAF) of 0.1% or greater, three pairs of primers were designed to yield overlapping amplicons. Each pair of primers comprised a forward and a reverse primer, with each primer having a nucleotide sequence complementary to a portion of the target nucleic acid sequence. Each primer had an adapter at or near its 5′ terminus and upstream from its complementary nucleic acid sequence. The adapter's nucleic acid sequence was complementary to a nucleic acid sequence used in a Next Generation Sequencing (NGS) platform, such as Ion Torrent or Illumina's MiSeq. Additionally, the reverse primer for each pair of primers further comprised an index sequence upstream from the primer's complementary nucleic acid sequence that was unique to the pair.
Three distinct amplification reactions were prepared, each comprising one of the three pairs of primers. The reactions comprised 1.0 μM primers, 1× final concentration of 5× Phusion High-fidelity Buffer (NEB), 200 μM dNTPs, 1.0 units of Phusion High-fidelity Polymerase (NEB), and about 25 to 50 ng of template DNA. The reactions were subjected to an initial denaturation step of 30 seconds at 98° C. followed by 20 cycles of 98° C. (denaturing the template DNA) for 10 seconds, 62° C. (annealing the primers to the template nucleic acid) for 20 seconds, and 72° C. (to extend the DNA product) for 30 seconds. After cycling, the reactions were subjected to an additional 10 minutes at 72° C. as a final extension step.
5 μl of each PCR product were then pooled and purified using a ThermoFisher MagJet purification kit (any kit that removes products <100 base pairs in length can be used). The purified reaction products were resuspended in 20 μl of water, mixed, and incubated for two minutes. The reactions were then placed on a magnet for two minutes, and the eluted DNA was removed. About 1 μl was run on a TapeStation or a Bioanalyzer 2100 to confirm quality.
Aliquots of the amplicons generated from a single round of amplification were analyzed on a Bioanalyzer 2100. This amplification strategy yielded detectable amplicons at the expected time point (i.e., between 50 and 60 seconds for the control (
The purified PCR reaction products were sequenced using the Ion Torrent system (ThermoFisher Scientific) to generate sequencing reads that comprise the nucleic acid sequence of the target nucleic acid. The sequencing reads were demultiplexed, or segregated, into different bins depending on the detected index sequence. Table 5 provides a summary of the observed alternate allele fractions detected using this method.
To identify low frequency genetic variation in a target nucleic acid sequence with an alternate allele fraction of 0.025% or greater, three pairs of primers were designed to yield overlapping amplicons. Each pair of primers comprised a forward and a reverse primer, with each primer having a nucleotide sequence complementary to a portion of the target nucleic acid sequence. Each primer had an adapter at or near its 5′ terminus and upstream from its complementary nucleic acid sequence. The adapter's nucleic acid sequence was complementary to a nucleic acid sequence used in an NGS platform, such as Ion Torrent or Illumina's MiSeq. Each individual reverse primer further comprised an index sequence upstream from the primer's complementary nucleic acid sequence. Additionally, each individual forward or reverse primer in each pair of primers further comprised a unique molecular identifier (UMI). No two primers had the same UMI.
Three distinct amplification reactions were prepared, each comprising one of the three pairs of primers. The reactions comprised 1.0 μM primers, 1× final concentration of 5× Phusion High-fidelity Buffer (NEB), 200 μM deoxynucleotide triphosphates (dNTPs), 0.1 μl of 0.4 mM Biotin-14-dCTP, 1.0 units of Phusion High-fidelity Polymerase (NEB), and about 25 to 50 ng of template DNA. The reactions were subjected to an initial denaturation step of 30 seconds at 98° C. followed by 8 cycles of 98° C. (denaturing the template DNA) 10 seconds, 62° C. (annealing the primers to the template nucleic acid) for 20 seconds, and 72° C. (to extend the DNA product) for 30 seconds. After cycling, the reactions were subjected to an additional 10 minutes at 72° C. as a final extension step. The reaction products, or amplicons, were purified by washing 5 μl of MyOne C1 streptavidin beads two times with 1× Binding-Washing (B&W) buffer and then resuspending the beads in 25 μl of 2×B&W buffer. 25 μl of the MyOne C1 streptavidin beads was then added to 25 μl of the PCR amplicon and incubated at room temperature for 15 minutes with mixing. The mixture was exposed to a magnet, which isolates the beads with the amplicons bound thereto. The supernatant was removed, and 500 μl 1× B&W buffer was added to the beads, mixed, and exposed to the magnet. Again, the supernatant is removed, and the wash was repeated. The beads were finally resuspended in 28 μl water. Some reaction products were purified using an exonuclease 1/shrimp alkaline phosphatase (ExoSap) enzymatic purification protocol, wherein 8 μl of the commercially available ExoSap-It reagent (ThermoFisher) was added to the 20 μl amplification reaction and incubated at 37° C. for 15 minutes followed by 80° C. for 15 minutes.
While the amplicons were attached to the streptavidin beads, an additional amplification was performed to enhance the copy number of the bound amplicons. Briefly, the additional amplification reactions comprised 1.0 μM primers, 1× final concentration of 5× Phusion High-fidelity Buffer (NEB), 200 μM deoxynucleotide triphosphates (dNTPs), 0.1 μl of 0.4 mM Biotin-14-dCTP, 1.0 units of Phusion High-fidelity Polymerase (NEB), and about 25 to 50 ng of template DNA. The reactions were subjected to an initial denaturation step of 30 seconds at 98° C. followed by 20 cycles of 98° C. (denaturing the template DNA) 10 seconds, 62° C. (annealing the primers to the template nucleic acid) for 20 seconds, and 72° C. (to extend the DNA product) for 30 seconds. After cycling, the reactions were subjected to an additional 10 minutes at 72° C. as a final extension step, and 5 μl of the PCR reactions were pooled. A ThermoFisher MagJet purification kit that removes products <100 base pairs in length was used to purify the amplicons. Specifically, the amplicons in the pooled reactions were bound to streptavidin beads, and the supernatant was removed. The beads were then resuspended in 200 of water, mixed, and incubated for two minutes. The mixture was then exposed to a magnet for two minutes, and the eluted DNA was captured.
Referring to
After determining the concentration of the eluted DNA, it was diluted to 100 pM, and the purified PCR reaction products were sequenced using the Ion Torrent system (ThermoFisher Scientific).
The sensitivity and reproducibility of the methods described herein were assessed through serial dilutions of known germline mutations and known somatic mutations across a spectrum of alternative allele fractions. A comparison of alternative allele fractions with other known detections strategies including whole genome sequencing, whole exome sequencing, targeted sequencing, Sanger sequencing with Topo-cloning, and ddPCR was performed. First, triplicate primers (i.e., 3 unique pairs of primers) were designed as described in the methods for known germline mutations occurring in both the autosomal and X-chromosomal regions, including both heterozygous and hemizygous alleles. Twelve serial dilutions were sequenced on the Ion Torrent S5 with 400 base pair reads using six unique barcodes per primer. All reads were processed using custom analytical scripts (described in methods), allowing the for comparison of assessed and expected allelic fractions.
Referring to
Given that input DNA is often limited but is also known as an important factor for sensitivity for somatic alleles, decreased inputs of DNA were tested to determine if they could achieve a similar level of precision under the same dilution curve. Indeed, while decreased input DNA does impact the sensitivity, alternative allele fractions down to 0.05% remain detectable, though at a slightly elevated standard deviation among the triplicate primes for the lowest alternative allele fractions of 0.05%, indicating that when validating alleles below 0.1% alternative allele fractions, increased input DNA could improve precision. Furthermore, the impact of total sequencing depth on the accuracy was assessed to identify the minimum depth needed for accurate determination of alternative allele fractions. Using random sampling of the initial raw unmapped data, a strong correlation of read depths above threshold level can be made, and sequencing beyond this threshold will provide minimal benefits on the precision of the alternative allele fraction assessment.
Frozen postmortem human brain specimens from 61 autism spectrum disorder cases and 15 neurotypical controls were obtained for analysis. DNA was extracted from dorsolateral prefrontal cortex where available (or generic cortex in a minority of cases) using lysis buffer from the QIAamp DNA Mini kit (Qiagen) followed by phenol chloroform extraction and isopropanol cleanup. Samples UMB4334, UMB4899, UMB4999, UMB5027, UMB5115, UMB5176, UMB5297, UMB5302, UMB1638, UMB4671, and UMB797 were processed using TruSeq Nano DNA library preparation (Illumina) followed by Illumina HiSeq X Ten sequencing to a minimum 200× depth. All remaining samples were processed using TruSeq DNA PCR-Free library preparation (Illumina) followed by minimum 30× sequencing of seven separate libraries on the Illumina HiSeq X Ten, for a total minimum coverage of 210× per sample. An average of 251× depth was achieved across all samples, using 150 base pair paired-end reads. Two samples, UMB5771 and UMB5939, had parental saliva-derived DNA available, and DNA from both parents for these two cases was obtained and sequenced to about 50× depth. Parental DNA was not available for any other samples. Additionally, DNA was extracted from Brodmann Area 17 (occipital lobe) for cases UMB4638 and UMB4643 and sequenced at Macrogen to a minimum 210× depth following PCR-free library preparation. Bulk heart and liver sequencing data, as well as single-cell sequencing data from three individuals (UMB1465, UMB4643, and UMB4638) were used in this study.
All paired-end FASTQ files were aligned using BWA-MEM version 0.7.8 to the GRCh37 human reference genome including the hs37d5 decoy sequence from the Broad Institute, following GATK best practices (software.broadinstitute.org/gatk/best-practices/). Mutect2-PoN was used to generate two pairs of panel-of-normals (PoN) by using 60 autism spectrum disorder samples or 15 control samples to remove sequencing artifacts and germline variants from the other group. Rare variants were further selected by filtering out any variant with a maximum population minor allele frequency >0.001 in any of Kaviar, 1000 Genomes, EVS6500 (evs.gs.washington.edu/EVS/), ExACnonpsych, or gnomAD (gnomad.broadinstitute.org/). Repetitive region variants were removed using RepeatMasker (www.repeatmasker.org/), and variants within segmental duplication regions or shared between multiple individuals were also removed. Low-quality calls tagged “t_lod_fstar,” “str_contraction,” and “triallelic_site” were removed. For analysis of damaging heterozygous variants, variants were identified in the 78 risk genes previously used.
For somatic mutation detection, a minimum alternate (or variant) allele fraction (AAF or VAF) of 0.03 was required unless a variant was phasable by Mutect2, which allowed for rescue of variants down to an alternate allele fraction of 0.02. Low-quality calls tagged “triallelic_site” were removed. A minimum alternate read depth of four reads was required. Only private events among the population were analyzed. An upper alternate allele fraction threshold of 0.40 was set and heterozygous germline variants were removed. Variants within repetitive regions were also removed, leaving 14,984 candidate somatic mutations. MosaicForecast was then used to perform read-backed phasing and identify high-confidence mosaics from the candidate call set. Briefly, features likely to be correlated with mosaic detection specificity were selected: mapping quality, base quality, clustering of mutations, read depth, number of mismatches per read, read1/read2 bias, strand bias, base position, read position, trinucleotide context, sequencing cycle, library preparation method, and genotype likelihood. Based on these features a random forest model was trained using phased variants. Further training was conducted using parental whole genome sequencing data from two cases UMB5771 and UMB5939 as well as single cell whole genome sequencing data from three control brains, UMB1465, UMB4643, and UMB4638 for which inherited germline mutations or variants present in multiple single cells at a low alternate allele fraction (averaging alternate allele fraction <0.30, likely representing sequencing or alignment artifact), supplied a training set of false positives. Predicted mosaics were further filtered by removing genomic regions enriched for low-alternate allele fraction variants and by removing variants with unusually high sequencing depth that also occurred in regions marked as copy number variants (CNVs) by Meerkat. Following all training and filtration, 1143 putative mosaic variants were identified. One autism spectrum disorder sample, MSSM007, was eliminated from the study due to very high noise suggestive of contamination or sequencing artifact.
Pathogenicity prediction scores were calculated for functional mosaic and germline variants using SIFT, PolyPhen-2, MutationTaster, and CADD. To be considered damaging, a variant had to be predicted as damaging or probably damaging (or CADD phred score >20) by at least three out of four prediction tools. Mutations in genes were checked for overlap with the Simons Foundation Autism Research Initiative (SFARI) database of autism spectrum disorder—relevant genes (gene.sfari.org/), and with the Online Inheritance in Man (OMIM) database of genes with relevance to any human disease (www.omim.org/).
Targeted validation was attempted on 243 of 1143 possible mosaic variants. PCR primers were designed for each variant and synthesized with Ion Torrent adapters P and A, with barcodes added for unique identification. PCR amplification was performed using Phusion HotStart II DNA Polymerase (Thermo) as described by the manufacturer, with 20-25 cycles of amplification. Reactions were pooled and purified with AMPure XP technology (Agencourt), then sequenced on the Ion Torrent Personal Genome Machine using the Ion 530 chip with 400 base pair reads, reaching an average coverage of 118,000 reads per variant amongst reactions that yielded mappable reads. Following demultiplexing and trimming, reads were mapped using BWAMEM (a Burrows-Wheeler aligner algorithm) and locally realigned using GATK. BAM files were then imported into a CLC Genomics workbench (Qiagen) and mosaic variants were identified using the following filters: minimum frequency 0.05%, minimum depth 10,000× per reaction, minimum count 50, required significance 0.1%, central and neighborhood base quality of >15, and 3-nucleotide homopolymer filtration. Variants were then classified as validated true mosaics (198 variants), homozygous reference with variant not present (21 variants), germline heterozygous (1 variant), PCR reactions failed to amplify (19 variants), or undetermined (4 variants). The “undetermined” designation was used for variants for which the originally sequenced DNA was not available, so validation was conducted on a separate DNA extraction that could have slightly different clonal architecture. It was also used to classify two variants in which sequencing noise precluded validation interpretation. Validation success rates were calculated as the number of true mosaics divided by the sum of true mosaics, homozygous reference, and germline heterozygous. Weighted averaging across PCR and PCR-free variant validation was used to determine a comprehensive validation rate of 93%. Five variants from UMB5771 and UMB5939 were also re-sequenced in parent DNA, which confirmed a mosaic state in the offspring and homozygous reference in parents.
A deleterious missense C to A change in the autism spectrum disorder risk gene CACNA1A was called in 5.2% of sequencing reads in case UMB1174 (
Ion Torrent amplicon resequencing for 34 germline heterozygous mutations revealed that alternate allele frequencies were slightly over-dispersed compared to a binomial distribution (
The triple-primer PCR sequencing method substantially increases the throughput and sensitivity for the detection and validation of somatic mutations (
While numerous studies have sought to define the error rates for the Ion Torrent platform due to the potential increased rate of insertion and deletion errors, particularly at homopolymers, the exact error rate appears to vary from sample to sample. Even more, while the rate of indel errors is likely elevated in the Ion Torrent platform over Illumina technology, the rates of SNV errors appear to be similar. It is likely that many estimates of errors are compounded by the combined effects of polymerase induced errors, mapping issues, and sequencing artifacts, all of which are known to reduce the sensitivity of detecting somatic mutations present in low fractions of a sample. Therefore, triple-primer PCR sequencing was developed to assess and partially mitigate these errors, while leveraging the rates to provide statistical confidence about a given mutation.
Prior studies have demonstrated the method of validating low AAF alleles using ultra-deep amplicon sequencing. However, technical issues including allelic dropout, artifacts (e.g., PCR- and sequencing platform-induced) and PCR duplicates can reduce the accuracy detected AAFs and possible result in both false negative calls as well as skewed AAFs. Triple-primer PCR sequencing overcomes these limitations through the use of multiple unique primers that are specifically designed to prevent sharing binding sites while avoiding known mutations (i.e., individual specific and general population) but are within 250 nucleotides (nts) of the target mutation. Once designed, unique primer-specific barcodes are appended to the reverse primers, along with Ion Torrent adapters. Optionally, Illumina adapters and/or 10 nt molecular barcodes can be appended to the primers to improve sensitivity or usage on the Illumina platform. Customized primers amplify targets including the mutation or region of interest using reduced cycling and minimal amounts of DNA, and amplification products are sequenced on either the Ion Torrent S5 or Illumina MiSeq platform for ultra-deep coverage. This optimized process allows for independent analyses of each primer pair, determination of error rates bases on amplicon-specific error rates (i.e., level of PCR and sequencing induced artifacts across the amplicon), identification of allelic imbalances from additional mutations affecting primer binding or chromatin structure, and the assessment of the variation in AAF among primers. Together, these steps provide a robust and low-cost strategy for extremely precise estimation of AAFs which is broadly applicable to studies of somatic and germline mutations.
As the utility of the presently described invention relies on overcoming the previously described limitations of somatic mutation detection, triplicate unique primer sets were first designed around 5 known germline mutations (Tables 6A-6C) previously identified in bulk genomic DNA for testing the error rates of the method. The reduced PCR cycling conditions with a high-fidelity polymerase (4.4×10−7; Phusion HS, ThermoFisher) is estimated to result in an error rate of 8.8×10−6 at any given nucleotide position (ThermoFisher PCR Fidelity Calculator). Given that error rates vary amongst amplicons due to the specific nucleotide content of each amplicon, an internal control was designed for assigning the significance of each identified mutation. Using these primers, background error rates from PCR and sequencing, the sensitivity to detect extremely low AAFs, accuracy of the ascertained AAF measurement, and required DNA input and sequencing depths were assessed.
First, reads and nucleotides were stringently filtered for nucleotide and mapping qualities (q>20 and Q>20), resulting in the removal of an average of 10% of bases at any given nucleotide position. Relaxing these parameters (e.g., q10, Q10) did not decrease the fraction of excluded sites or assessed AAF, supporting that most nucleotide positions are of high quality. Next, the rate of artifacts in the region of the amplicon surrounding the mutation of interest was assessed by the AAF of all alternate alleles at each position under the assumption that all non-reference high-quality alleles present at sites not known to have a mutation represent errors. Across all amplicons, a low average background mutation frequency (0.018% AAF+/−0.0067%) was found for nucleotides located in the flanking 50 nt on either site of a mutation. Consistent with prior studies, some amplicons exhibited positional variability in error rates due to mapping errors around indels, including artifacts arising during sequencing.
To further reduce the rate of indel-associated errors, a computational modeling approach that detects and corrects sequencing platform errors was incorporated. Specifically, Pollux, a recent error modeling algorithm that screens for and corrects an estimated >95% of all indel associated errors, was used. The correction of indel-associated errors resulted in nearly a 5-fold reduction in nucleotide error frequency (0.0034%+/−0.0009%), allowing for mutations at extremely low AAFs to be distinguished from background sequencing and PCR-induced artifacts.
The AAF of somatic mutations can vary dramatically across tissues, where they can be nearly undetectable in tissues such as blood, but higher frequency in tissues like the brain. Given that most genetic testing is performed on blood or cell free DNA samples with anticipated low AAFs, the ability of the presently described methods to accurately detect AAFs at extremely low levels, which are often difficult or impossible to accurately assess by other methods.
The sensitivity of triple-primer PCR sequencing was assessed through serial dilution of a genomic control DNA sample containing the same 5 known germline mutations described above (Tables 6A-6C) with a control DNA lacking these mutations, thereby generating AAFs ranging from 50% down to 0.01%. The dilutions were amplified with primers for each mutation and sequenced on the Ion Torrent S5 with sequencing reads of 400 bp in length. All reads were processed using custom analytical scripts (described in methods), allowing the comparison of assessed and expected allelic fractions.
The presently described method accurately measures AAFs as low as 0.01% when using a 50 ng of genomic DNA, although for significant detection above the amplicon-specific error rates, AAFs were typically required to be above 0.05% (
The measured AAFs (average across triple primer sets) were linearly correlated with the expected AAFs down to 0.01% (R2>0.999), though as expected, individual AAFs do vary amongst individual primers (R2>0.98). Therefore, while individual primer sets are prone to biases in AAFs, the utility of multiple primer provides a robust and accurate measurement.
DNA is often limited, particularly in clinical contexts, but is also known as an important factor for sensitivity for somatic alleles due to the presence of fewer DNA fragments containing the targeted allele. Therefore, the sensitivity of using 50 ng was compared to using a reduced concentration of 25 ng (˜3800 cells) (PMID: 30813969). With 3800 cells, the accurate detection of the lowest dilution of 0.01% AAF is unlikely as it would likely only be represented by a single fragment. Surprisingly, AAFs down to 0.05% remained detectable with 25 ng DNA (
Furthermore, the impact of total sequencing depth on the accuracy was assessed to identify the minimum depth needed for accurate determination of AAFs. Sequencing data for each amplicon were randomly sampled to create artificial datasets containing a wide range of depths ranging from 10,000 to 150,000× coverage. Increasing read depths above 10,000× did not have a substantial impact on the background error rates within the amplicons. Even more, a minimum depth of 10,000× was able to accurately measure AAFs down to 0.1% with no improvement with elevated coverage. However, accurate measurement of AAFs below 0.1% required depths of 25,000× to ensure significance over the background errors. Overall, a strong correlation was found of AAFs measured across a wide range of read depths, indicating that detection of AAFs of 0.01% is possible at depths greater than above 25,000×.
The assessment of error rates and the potential for false positive allele calls was extended by performing similar sequencing on DNA samples lacking mutations. As expected, these alleles were not detectable, with only the typical background error rate being detected, which is often not the same allele as the mutation, supporting the specificity of this method.
As some tissues are more difficult to work with, the ability was assessed of the method to accurately detect known mosaic alleles that were previously identified in blood and brain tissue by a range of methods including WGS, WES, and targeted Illumina sequencing. Even more, given the importance of validating indels and the elevated indels error rates on Ion Torrent data, >50 somatic indels were tested using the method of the present invention with a direct comparison of the sites between the DNA sample containing the mutation and a control sample. It was demonstrated that AAFs of SNVs (R=0.93, (
The known increased error rates for indel in Ion Torrent data and the inability to utilize PCR duplicate information may limit the ability to quantitate some ultra-rare alleles (<0.05% AAF) and indels. Even more, the Pollux software is known to overcorrect for indels and has difficulty distinguishing rare indels from artifacts. Despite these limitations, it was assessed how the method performs on a wide range of indels occurring at AAFs from 1% to 30% and 1 to 21 base pairs in length, including 40 insertions and 60 deletions previously identified using 200× whole genome sequencing. Even more importantly, these mutations were not identified in control DNA, where at these sites very low error rates for indels (0.010%±0.05%) were found, supporting that even the single base indels are not being introduced by PCR or the Ion Torrent. These data indicate a sensitivity to accurately quantitate AAFs of indels down to 0.05% in many instances. Despite that many of these mutations were detected using only a few reads in the WGS data, a strong correlation was found between the predicted AAFs in the WGS and the measured values by the method described in this example (
To further improve the sensitivity for low AAFs, a modified version of the protocol was performed (
The incorporation of biotin into the PCR product did not impact the overall measured AAFs, but slightly reduced the error rate (0.0023%±0.0011% AAF), possibly due to the ability to perform better purification and the use of a common primer for the majority of the amplifications. These indicate that a 2-step UMI approach for the method is valuable in situations requiring reduced error rates for ultra-low AAFs or where PCR duplicates may be of particular concern.
The increased sensitivity of the the presently described approach can be further applied for the detection of novel ultra-low AAFs variants with Illumina-based sequencing. Overlapping primers were developed so that all regions of the PRNP gene was covered by at least 3 independent amplicons, each containing Illumina sequencing adapters and UMIs. Using the 2-step PCR approach, sequencing libraries were prepared for a dilution series of a known mutation (5%, 0.5%, and 0.05% AAFs) and additional samples were screened for novel alleles. While any given amplicon can have some errors, as outlined above and previously documented in amplicon-based sequencing studies, it was contemplated whether the method could reduce such effects to identify high-confidence mutations. By requiring consistent AAFs across multiple unique primer sets, the AAFs of mutations were accurately measured down to at least 0.05% (
The following materials and methods were used in carrying out this example.
At least three unique sets of primers were designed for each mutation by extracting the flanking sequence around each mutation so that the mutation is located at different positions within each of the three sequences. Next, common alleles are masked, along with the targeted mutation and flanking 5bps on each site using the bedtools maskfasta tool. The masked multi-fasta file containing all sequences for targeted alleles are input into BatchPrimer webtool to design primers for each sequence. Primers are designed to an average TM of 60° C., with a minimum of 59° C. and maximum of 62° C. The amplicon length is dependent on the specific mutation and DNA sources. For example, difficult to map regions may have longer products while degraded DNA samples may require shorter amplicons. In general, to ensure that all primers are likely unique and of similar amplicon length, amplicons have a target length of 225-300 bp in length. The primer sequences are checked by BLAT and in-silico PCR to ensure both their unique amplificon in the genome and that the primer binding sites do not overlap between any set of primers. The final set of primers are then uniquely barcoded using 10 nt barcodes and if desired, an additional 10 nt UMI is added. Finally, Ion Torrent specific adapter sequences are appended to the forward and reverse primers, allowing for their direct sequencing.
For the standard, single step PCR sequencing method described above, PCR was performed using 20 cycles on a 25 μl reaction mix containing either 25 or 50 ng of input DNA sample, Phusion Hot-Start polymerase, dNTPs, HC-Buffer, and the primers. For initial testing, 30 cycles of enrichment were used to ensure only a single amplicon is produced. The high-sensitivity method modifies this process by reduction of the PCR cycling to 5 and the incorporation of 0.1 μL of 0.4 mM biotin-14-dCTP into the reaction mix. Biotinylated PCR amplicons are captured by adding 5 μl of washed Strepatvidin Myone beads resuspended in 25 μl of 2× binding and washing buffer. The mixture is incubated at room temperature with gentle mixing for 15 minutes and placed on a 96-well magnetic plate. The liquid was removed and the beads were washed one time with 1× binding and washing buffer. Then beads are then resuspended in 25 μl PCR reaction mixture containing custom primers which preserve the original UMI sequences, Phusion Hot-Start polymerase, dNTPs, and HC-Buffer. The biotin labeled product was amplified with an additional 20 cycles of enrichment before the beads were removed. Enriched products were pools at equal volumes and purified using the MagJet purification kit.
Purified library pools are analyzed for enrichment efficiency and the complete removal of primers through by either the Agilent Bioanalyzer Hi-sensitivity chip or the TapeStation. The concentration was determined using PicoGreen. Pools were diluted to a final concentration of 100 pM prior to sequencing on the 430 chip for the Ion Torrent S5.
Raw unmapped bam files were obtained for each run and were processed using our custom analyses pipeline. First, all BAMs are converted to a fastq fiel using bedtools bamtofastq tool. Then, quality and adapter trimming was performed using cutadapt tool. Next, samples lacking UMIs, are demultiplexed using fastx_barcode_splitter, resulting in separate fastq files for each primer set. The barcode sequences are removed from the sequences using cutadapt. If the allele being tested in an SNV, indel correction is performed using Pollux. Finally, all samples are aligned to the reference genome using BWA-mem.
Variants are then called across the length of each amplicon though the use of samtools mPileup with the settings: q=20, Q=20. The resulting vcfs are parsed into a file containing the flanking 50 nt positions on each side of the variant and a separate file for the allele of interest. The average allele frequency across the flanking regions are then compared to the average AAF of the mutation across the 3 unique primers.
From the foregoing description, it will be apparent that variations and modifications may be made to the invention described herein to adopt it to various usages and conditions. Such embodiments are also within the scope of the following claims.
The recitation of a listing of elements in any definition of a variable herein includes definitions of that variable as any single element or combination (or subcombination) of listed elements. The recitation of an embodiment herein includes that embodiment as any single embodiment or in combination with any other embodiments or portions thereof.
All patents and publications mentioned in this specification are herein incorporated by reference to the same extent as if each independent patent and publication was specifically and individually indicated to be incorporated by reference.
This application claims the benefit of the following U.S. Provisional Application No. 62/799,671, filed Jan. 31, 2019, the entire contents of which are incorporated herein by reference.
This invention was made with government support under Grant Nos. R01NS032457 and U01MH106883 awarded by the National Institutes of Health. The government has certain rights in the invention.
| Filing Document | Filing Date | Country | Kind |
|---|---|---|---|
| PCT/US2019/063447 | 11/26/2019 | WO | 00 |
| Number | Date | Country | |
|---|---|---|---|
| 62799671 | Jan 2019 | US |
