Methods and systems for detecting genetic variants
Disclosed herein in are methods and systems for determining genetic variants (e.g., copy number variation) in a polynucleotide sample. A method for determining copy number variations includes tagging double-stranded polynucleotides with duplex tags, sequencing polynucleotides from the sample and estimating total number of polynucleotides mapping to selected genetic loci. The estimate of total number of polynucleotides can involve estimating the number of double-stranded polynucleotides in the original sample for which no sequence reads are generated. This number can be generated using the number of polynucleotides for which reads for both complementary strands are detected and reads for which only one of the two complementary strands is detected.
1. A method for treating a subject having lung cancer, comprising:
(a) determining that a sample from the subject comprises a genetic variant associated with a T790M EGFR mutation, wherein the sample comprises double-stranded deoxyribonucleic acid (DNA) molecules from the subject, by:
performing or having performed a diagnostic assay on the sample to determine the presence of the genetic variant associated with the EGFR mutation, wherein the diagnostic assay comprises:
i) attaching adapters comprising molecular barcodes to a plurality of the double-strand DNA molecules from the sample to produce tagged parent polynucleotides;
ii) amplifying a plurality of the tagged parent polynucleotides to produce amplified progeny polynucleotides;
iii) sequencing a plurality of the amplified progeny polynucleotides to produce a set of sequencing reads;
iv) mapping a plurality of the sequencing reads from the set of sequencing reads to a reference sequence to generate mapped sequencing reads;
v) reducing or tracking redundancy in a plurality of the mapped sequencing reads using at least sequence information from the molecular barcodes to generate a plurality of consensus sequences representative of original DNA molecules from among the tagged parent polynucleotides, wherein the plurality of consensus sequences is generated from (i) paired reads corresponding to mapped sequencing reads generated from a first tagged strand and a second tagged complementary strand derived from a DNA molecule from among the tagged parent polynucleotides, and (ii) unpaired reads corresponding to sequence reads generated from a first tagged strand having no second tagged complementary strand derived from a DNA molecule from among the tagged parent polynucleotides;
vi) detecting the genetic variant associated with the T790M EGFR mutation from among the plurality of consensus sequences; and
(b) administering a therapy responsive to the T790M EGFR mutation to the subject to treat the lung cancer based on detecting the presence of the genetic variant associated with the T790M EGFR mutation.
2. The method of claim 1 , wherein the lung cancer is non-small cell lung cancer.
3. The method of claim 1 , wherein the plurality of the double-stranded DNA molecules in i) comprises cell-free deoxyribonucleic acid (cfDNA) molecules.
4. The method of claim 2 , wherein the plurality of the double-stranded DNA molecules in i) comprises cell-free deoxyribonucleic acid (cfDNA) molecules.
5. The method of claim 3 , wherein the plurality of double-stranded cfDNA molecules comprises between 1 nanogram (ng) to 100 ng of double-stranded cfDNA molecules.
6. The method of claim 5 , wherein the attaching comprises using more than a 20× excess of adapters as compared to the double-stranded cfDNA molecules in the plurality of double-stranded cfDNA to attach at least 20% of the double-stranded cfDNA molecules with adapters.
7. The method of claim 6 , wherein the attaching comprises using more than a 50× excess of adapters as compared to the double-stranded cfDNA molecules in the plurality of double-stranded cfDNA.
8. The method of claim 6 , wherein the adapters are Y-shaped.
9. The method of claim 8 , wherein the plurality of the double-stranded cfDNA molecules is tagged with n different combinations of molecular barcodes, wherein n is at least 2 and no more than 10,000*z, wherein z is a mean of an expected number of duplicate molecules in the plurality of DNA molecules that map to identical start and stop positions on the reference sequence.
10. The method of claim 9 , wherein the plurality of the double-stranded cfDNA molecules is tagged with n different combinations of molecular barcodes, wherein n is at least 2 and no more than 1,000*z.
11. The method of claim 10 , wherein the plurality of the double-stranded cfDNA molecules is tagged with n different combinations of molecular barcodes, wherein n is at least 2 and no more than 100*z.
12. The method of claim 9 , wherein z is between 2 and 8.
13. The method of claim 6 , wherein the reducing or tracking redundancy in a plurality of the mapped sequencing reads is performed using sequence information from the start base position and the stop base position of the mapped sequencing reads in addition to the sequence information from the molecular barcodes.
14. The method of claim 1 , comprising selectively enriching the amplified progeny polynucleotides for a plurality of genomic regions of interest.
15. The method of claim 14 , wherein the plurality of genomic regions of interest comprises sequences of an EGFR gene.
16. The method of claim 1 , further comprising determining a first quantitative measure of paired reads and a second quantitative measure of unpaired reads.
17. The method of claim 13 , further comprising determining a total number of tagged parent polynucleotides that map to a locus of the reference sequence in the plurality of tagged parent polynucleotides based on the first quantitative measure and the second quantitative measure.
18. The method of claim 17 , further comprising determining a third quantitative measure of unseen sequence reads for which neither strand of a double-stranded DNA molecule is detected, wherein the third quantitative measure is calculated based on the first quantitative measure and the second quantitative measure.
19. The method of claim 18 , further comprising determining a total number of tagged parent polynucleotides in the plurality of tagged parent polynucleotides based on the first quantitative measure, the second quantitative measure, and the third quantitative measure.
20. The method of claim 1 , wherein the therapy comprises ruxolitinib.
21. A method for treating a subject having lung cancer, comprising:
(a) determining that a sample from the subject comprises a genetic variant associated with a T790M EGFR mutation, wherein the subject has previously undergone treatment with a therapy comprising erlotinib, wherein the sample comprises double-stranded deoxyribonucleic acid (DNA) molecules from the subject, by:
performing or having performed a diagnostic assay on the sample to determine the presence of the genetic variant associated with the EGFR mutation, wherein the diagnostic assay comprises:
i) attaching adapters comprising molecular barcodes to a plurality of the double-strand DNA molecules from the sample to produce tagged parent polynucleotides;
ii) amplifying a plurality of the tagged parent polynucleotides to produce amplified progeny polynucleotides;
iii) sequencing a plurality of the amplified progeny polynucleotides to produce a set of sequencing reads;
iv) mapping a plurality of the sequencing reads from the set of sequencing reads to a reference sequence to generate mapped sequencing reads;
v) reducing or tracking redundancy in a plurality of the mapped sequencing reads using at least sequence information from the molecular barcodes to generate a plurality of consensus sequences representative of original DNA molecules from among the tagged parent polynucleotides, wherein the plurality of consensus sequences is generated from (i) paired reads corresponding to mapped sequencing reads generated from a first tagged strand and a second tagged complementary strand derived from a DNA molecule from among the tagged parent polynucleotides, and (ii) unpaired reads corresponding to sequence reads generated from a first tagged strand having no second tagged complementary strand derived from a DNA molecule from among the tagged parent polynucleotides;
vi) detecting the genetic variant associated with the T790M EGFR mutation from among the plurality of consensus sequences; and
(b) administering a therapy other than erlotinib that is responsive to the T790M EGFR mutation to the subject to treat the lung cancer based on detecting the presence of the genetic variant associated with the T790M EGFR mutation, wherein the therapy other than erlotinib that is responsive to the T790M EGFR mutation is administered when the genetic variant associated with the T790M EGFR mutation rises to a predetermined clinical level that indicates resistance to erlotinib.
22. The method of claim 21 , wherein the lung cancer is non-small cell lung cancer.
23. The method of claim 22 , wherein the plurality of the double-stranded DNA molecules in i) comprises cell-free deoxyribonucleic acid (cfDNA) molecules.
24. The method of claim 23 , wherein the plurality of double-stranded cfDNA molecules comprises between 1 nanogram (ng) to 100 ng of double-stranded cfDNA molecules.
25. The method of claim 24 , wherein the attaching comprises using more than a 20× excess of adapters as compared to the double-stranded cfDNA molecules in the plurality of double-stranded cfDNA to attach at least 20% of the double-stranded cfDNA molecules with adapters.
26. The method of claim 25 , wherein the attaching comprises using more than a 50× excess of adapters as compared to the double-stranded cfDNA molecules in the plurality of double-stranded cfDNA.
27. The method of claim 21 , wherein the plurality of the double-stranded DNA molecules is tagged with n different combinations of molecular barcodes, wherein n is at least 2 and no more than 10,000*z, wherein z is a mean of an expected number of duplicate molecules in the plurality of DNA molecules that map to identical start and stop positions on the reference sequence.