Methods and systems for detecting sequence variants
The invention provides methods for identifying rare variants near a structural variation in a genetic sequence, for example, in a nucleic acid sample taken from a subject. The invention additionally includes methods for aligning reads (e.g., nucleic acid reads) to a reference sequence construct accounting for the structural variation, methods for building a reference sequence construct accounting for the structural variation or the structural variation and the rare variant, and systems that use the alignment methods to identify rare variants. The method is scalable, and can be used to align millions of reads to a construct thousands of bases long, or longer.
1. At least one non-transitory computer-readable storage medium storing processor-executable instructions that, when executed by a computer hardware processor, cause the computer hardware processor to perform a method for updating a genomic reference graph representing a plurality of genomic sequences, the plurality of genomic sequences including a first genomic variant, the method comprising:
accessing, from the at least one non-transitory computer-readable storage medium, the genomic reference graph, wherein the genomic reference graph comprises nodes and edges connecting the nodes, the nodes including a first node representing the first genomic variant;
obtaining nucleic acid sequence reads;
aligning the nucleic acid sequence reads to the genomic reference graph by, for a nucleic acid sequence read of the nucleic acid sequence reads:
determining overlap scores for the nucleic acid sequence read against the genomic reference graph, and
aligning the nucleic acid sequence read to the genomic reference graph based on the overlap scores;
identifying, based on the aligned nucleic acid sequence reads, a second genomic variant, wherein a number of base pairs separating locations of the first genomic variant and the second genomic variant in the genomic reference graph is below a threshold number of base pairs; and
updating the genomic reference graph to include a second node representing the second genomic variant.
2. The at least one non-transitory computer-readable storage medium of claim 1 , wherein the second genomic variant comprises fewer base pairs than the first genomic variant.
3. The at least one non-transitory computer-readable storage medium of claim 1 , wherein the first genomic variant is at least 100 base pairs long.
4. The at least one non-transitory computer-readable storage medium of claim 1 , wherein at least one nucleic acid sequence read of the nucleic acid sequence reads comprises the second genomic variant.
5. The at least one non-transitory computer-readable storage medium of claim 4 , wherein the second genomic variant occurs with a frequency of less than 5% in the nucleic acid sequence reads.
6. The at least one non-transitory computer-readable storage medium of claim 1 , wherein the threshold number of base pairs is 100 base pairs.
7. The at least one non-transitory computer-readable storage medium of claim 1 , wherein the genomic reference graph includes a third node representing:
a genomic sequence lacking the first genomic variant, at a same location as the first node in the updated genomic reference graph; or
a genomic sequence lacking the second genomic variant, at a same location as the second node in the updated genomic reference graph.
8. The at least one non-transitory computer-readable storage medium of claim 1 , wherein:
the first genomic variant comprises a deletion, duplication, copy-number variation, insertion, inversion, translocation, or any combination thereof; and
the second genomic variant comprises a single nucleotide polymorphism.
9. The at least one non-transitory computer-readable storage medium of claim 1 , wherein a path through the genomic reference graph represents a genome or a chromosome of an organism.
10. The at least one non-transitory computer-readable storage medium of claim 1 , wherein the method further comprises:
identifying, based on the aligned nucleic acid sequence reads, at least one additional genomic variant; and
updating the genomic reference graph to include at least one additional node representing the at least one additional genomic variant.
11. A system comprising:
a computer hardware processor; and
at least one non-transitory computer-readable storage medium storing processor-executable instructions that, when executed by the computer hardware processor, cause the computer hardware processor to perform a method for updating a genomic reference graph representing a plurality of genomic sequences, the plurality of genomic sequences including a first genomic variant, the method comprising:
accessing, from the at least one non-transitory computer-readable storage medium, the genomic reference graph, wherein the genomic reference graph comprises nodes and edges connecting the nodes, the nodes including a first node representing the first genomic variant;
obtaining nucleic acid sequence reads;
aligning the nucleic acid sequence reads to the genomic reference graph by, for a nucleic acid sequence read of the nucleic acid sequence reads:
determining overlap scores for the nucleic acid sequence read against the genomic reference graph, and
aligning the nucleic acid sequence read to the genomic reference graph based on the overlap scores;
identifying, based on the aligned nucleic acid sequence reads, a second genomic variant, wherein a number of base pairs separating locations of the first genomic variant and the second genomic variant in the genomic reference graph is below a threshold number of base pairs; and
updating the genomic reference graph to include a second node representing the second genomic variant.
12. The system of claim 11 , wherein the second genomic variant comprises fewer base pairs than the first genomic variant.
13. The system of claim 11 , wherein the threshold number of base pairs is 100 base pairs.
14. The system of claim 11 , wherein the genomic reference graph includes a third node representing:
a genomic sequence lacking the first genomic variant, at a same location as the first node in the updated genomic reference graph; or
a genomic sequence lacking the second genomic variant, at a same location as the second node in the updated genomic reference graph.
15. The system of claim 11 , wherein the method further comprises:
identifying, based on the aligned nucleic acid sequence reads, at least one additional genomic variant; and
updating the genomic reference graph to include at least one additional node representing the at least one additional genomic variant.
16. A method for updating a genomic reference graph representing a plurality of genomic sequences, the plurality of genomic sequences including a first genomic variant, the method comprising:
using at least one processor to perform:
accessing, from at least one non-transitory computer-readable storage medium, the genomic reference graph, wherein the genomic reference graph comprises nodes and edges connecting the nodes, the nodes including a first node representing the first genomic variant;
obtaining nucleic acid sequence reads;
aligning the nucleic acid sequence reads to the genomic reference graph at least in part by, for a nucleic acid sequence read of the nucleic acid sequence reads:
determining overlap scores for the nucleic acid sequence read against the genomic reference graph, and
aligning the nucleic acid sequence read to the genomic reference graph based on the overlap scores;
identifying, based on the aligned nucleic acid sequence reads, a second genomic variant, wherein a number of base pairs separating locations of the first genomic variant and the second genomic variant in the genomic reference graph is below a threshold number of base pairs; and
updating the genomic reference graph to include a second node representing the second genomic variant.
17. The method of claim 16 , wherein the second genomic variant comprises fewer base pairs than the first genomic variant.
18. The method of claim 16 , wherein the threshold number of base pairs is 100 base pairs.
19. The method of claim 16 , wherein the genomic reference graph includes a third node representing:
a genomic sequence lacking the first genomic variant, at a same location as the first node in the updated genomic reference graph; or
a genomic sequence lacking the second genomic variant, at a same location as the second node in the updated genomic reference graph.
20. The method of claim 16 , wherein the method further comprises:
identifying, based on the aligned nucleic acid sequence reads, at least one additional genomic variant; and
updating the genomic reference graph to include at least one additional node representing the at least one additional genomic variant.