Systems and methods for analyzing sequence data
The invention provides methods for comparing one set of genetic sequences to another without discarding any information within either set. A set of genetic sequences is represented using a directed acyclic graph (DAG) avoiding any unwarranted reduction to a linear data structure. The invention provides a way to align one sequence DAG to another to produce an alignment that can itself be stored as a DAG. DAG-to-DAG alignment is a natural choice wherever a set of genomic information consisting of more than one string needs to be compared to any non-linear reference. For example, a subpopulation DAG could be compared to a population DAG in order to compare the genetic features of that subpopulation to those of the population.
1. A method, comprising:
using at least one computer hardware processor to perform:
obtaining a plurality of sequence reads previously obtained from one or more genetic samples;
obtaining, from at least one non-transitory computer-readable storage medium, a directed acyclic graph (DAG) data structure, wherein the DAG data structure comprises a plurality of nodes connected by a plurality of edges, the plurality of nodes including a first node included in a first path of a plurality of paths in the DAG data structure, the first node representing a first nucleotide sequence stored as a first string of one or more symbols;
aligning at least some of the plurality of sequence reads to the DAG data structure;
identifying, based on results of the aligning, one or more paths of the plurality of paths, the identifying comprising:
determining support values for nodes in the first path based on the results of the aligning, the support values including a support value for the first node that is indicative of a number of sequence reads aligned to the first string;
determining a first support value for the first path based on the support values for the nodes in the first path; and
including the first path in the one or more paths when the first support value exceeds a first threshold; and
outputting at least one of the one or more identified paths.
2. The method of claim 1 , wherein the first support value represents a minimum of the support values.
3. The method of claim 1 , wherein the at least one of the one or more identified paths represents at least one genetic haplotype included in the one or more genetic samples.
4. The method of claim 1 , further comprising:
determining whether a number of the one or more identified paths exceeds a second threshold; and
outputting each of the one or more identified paths when the number of the one or more identified paths does not exceed the second threshold.
5. The method of claim 4 , further comprising:
outputting a subset of the one or more identified paths when the number of the one more identified paths exceeds the second threshold.
6. The method of claim 4 , wherein the second threshold is indicative of a number of human genomes represented in the one or more genetic samples.
7. The method of claim 5 , wherein outputting the subset of the one or more identified paths comprises outputting a number of paths equivalent to the second threshold.
8. The method of claim 5 , wherein the subset of the one or more identified paths includes paths associated with support values that are greater than support values associated with paths that are not included in the subset of the one or more identified paths.
9. The method of claim 1 , wherein the first threshold is indicative of a ratio of a number paths included in the plurality of paths and a number of the one or more genetic samples.
10. A system, comprising:
at least one computer hardware processor; and
at least one non-transitory computer-readable storage medium storing processor executable instructions that, when executed by the at least one computer hardware processor, cause the at least one computer hardware processor to perform:
obtaining a plurality of sequence reads previously obtained from one or more genetic samples;
obtaining, from the at least one non-transitory computer-readable storage medium, a directed acyclic graph (DAG) data structure, wherein the DAG data structure comprises a plurality of nodes connected by a plurality of edges, the plurality of nodes including a first node included in a first path of a plurality of paths in the DAG data structure, the first node representing a first nucleotide sequence stored as a first string of one or more symbols;
aligning at least some of the plurality of sequence reads to the DAG data structure;
identifying, based on results of the aligning, one or more paths of the plurality of paths, the identifying comprising:
determining support values for nodes in the first path based on the results of the aligning, the support values including a support value for the first node that is indicative of a number of sequence reads aligned to the first string;
determining a first support value for the first path based on the support values for the nodes in the first path; and
including the first path in the one or more paths when the first support value exceeds a first threshold; and
outputting at least one of the one or more identified paths.
11. The system of claim 10 , wherein the first support value represents a minimum of the support values.
12. The system of claim 10 , further comprising:
determining whether a number of the one or more identified paths exceeds a second threshold; and
outputting each of the one or more identified paths when the number of the one or more identified paths does not exceed the second threshold.
13. The system of claim 12 , further comprising:
outputting a subset of the one or more identified paths when the number of the one more identified paths exceeds the second threshold.
14. The system of claim 12 , wherein the second threshold is indicative of a number of human genomes represented in the one or more genetic samples.
15. The system of claim 13 , wherein outputting the subset of the one or more identified paths comprises outputting a number of paths equivalent to the second threshold.
16. At least one non-transitory computer-readable storage medium storing processor executable instructions that, when executed by at least one computer hardware processor, cause the at least one computer hardware processor to perform:
obtaining a plurality of sequence reads previously obtained from one or more genetic samples;
obtaining, from the at least one non-transitory computer-readable storage medium, a directed acyclic graph (DAG) data structure, wherein the DAG data structure comprises a plurality of nodes connected by a plurality of edges, the plurality of nodes including a first node included in a first path of a plurality of paths in the DAG data structure, the first node representing a first nucleotide sequence stored as a first string of one or more symbols;
aligning at least some of the plurality of sequence reads to the DAG data structure;
identifying, based on results of the aligning, one or more paths of the plurality of paths, the identifying comprising:
determining support values for nodes in the first path based on the results of the aligning, the support values including a support value for the first node that is indicative of a number of sequence reads aligned to the first string;
determining a first support value for the first path based on the support values for the nodes in the first path; and
including the first path in the one or more paths when the first support value exceeds a first threshold; and
outputting at least one of the one or more identified paths.
17. The at least one non-transitory computer-readable storage medium of claim 16 , wherein the first support value represents a minimum of the support values.
18. The at least one non-transitory computer-readable storage medium of claim 16 , further comprising:
determining whether a number of the one or more identified paths exceeds a second threshold; and
outputting each of the one or more identified paths when the number of the one or more identified paths does not exceed the second threshold.
19. The at least one non-transitory computer-readable storage medium of claim 18 , further comprising:
outputting a subset of the one or more identified paths when the number of the one more identified paths exceeds the second threshold.
20. The at least one non-transitory computer-readable storage medium of claim 18 , wherein the second threshold is indicative of a number of human genomes represented in the one or more genetic samples.