Methods and systems for genotyping genetic samples
The invention provides methods and system for making specific base calls at specific loci using a reference sequence construct, e.g., a directed acyclic graph (DAG) that represents known variants at each locus of the genome. Because the sequence reads are aligned to the DAG during alignment, the subsequent step of comparing a mutation, vis-a-vis the reference genome, to a table of known mutations can be eliminated. The disclosed methods and systems are notably efficient in dealing with structural variations within a genome or mutations that are within a structural variation.
1. A system for genotyping a genetic sample, the system comprising:
a processor; and
a tangible, non-transitory computer-readable storage medium storing processor-executable instructions that, when executed by the processor, cause the processor to perform a method comprising:
obtaining a directed acyclic graph (DAG) representing a reference sequence of a species and genetic variation of the reference sequence, wherein a first path through the DAG represents a first sequence that includes a first allele and is associated with a first string of one or more symbols, a second path through the DAG represents a second sequence that includes a second allele associated with a second string of one or more symbols, and the second sequence represents a first genetic structural variation;
aligning at least some of a plurality of sequence reads from a genetic sample to the DAG, wherein aligning the at least some of the plurality of sequence reads to the DAG comprises:
comparing a read string of one or more symbols corresponding to a first sequence read of the plurality of sequence reads to the first and second strings associated with the first and second sequences represented by the first and second paths, respectively;
determining a first overlap score corresponding to a first alignment, at a first position in the DAG, between the read string and the first and second strings, the first overlap score being determined based on, and indicative of, a degree of overlap between the read string and the first string and a degree of overlap between the read string and the second string;
determining a second overlap score corresponding to a second alignment, at a second position in the DAG, between the read string and the first and second strings, the second overlap score being determined based on, and indicative of, a degree of overlap between the read string and the first string and a degree of overlap between the read string and the second string; and
identifying one of the first and second overlap scores and therefrom aligning the first sequence read to the respective first or second position in the DAG;
determining, based on results of the aligning, whether the first sequence read aligns to the second sequence represented by the second path through the DAG; and
upon determining that the first sequence read aligns to the second sequence represented by the second path through the DAG, identifying a presence of the first genetic structural variation within the genetic sample.
2. The system of claim 1 , further comprising:
determining, based on results of the aligning, whether two or more sequence reads of the at least some of the plurality of sequence reads align to the second path; and
upon determining that the two or more sequence reads align to the second path, reporting the presence of the first genetic structural variation within the genetic sample.
3. The system of claim 1 , wherein the DAG further comprises a plurality of nodes, and wherein the first path comprises a first node of the plurality of nodes and the second path comprises a second node of the plurality of nodes, wherein the second node is not included in the first path.
4. The system of claim 1 , wherein at least a portion of the first sequence read includes the first genetic structural variation.
5. The system of claim 1 , wherein the first genetic structural variation comprises 3 consecutive nucleotides not present in a genome of the species.
6. The system of claim 1 , wherein the first genetic structural variation comprises a deletion, an insertion, or a polymorphism.
7. The system of claim 1 , wherein the method further comprises identifying a variant in the genetic sample within 100 base pairs from the first genetic structural variation.
8. The system of claim 1 , wherein a third path through the DAG represents a third sequence, and wherein a fourth path through the DAG represents a fourth sequence, the fourth sequence representing a second genetic structural variation.
9. The system of claim 1 , wherein the method further comprises determining, based on results of the aligning, a confidence value indicative of whether the first genetic structural variation is present within the genetic sample.
10. The system of claim 1 , wherein a third path through DAG represents a third sequence that includes a third allele and is associated with a third string of one or more symbols, wherein the third sequence represents a second genetic structural variation.
11. The system of claim 1 , wherein the first genetic structural variation is associated with a disease, and wherein the method further comprises identifying the disease as being associated with the genetic sample after identifying the presence of the first genetic structural variation within the genetic sample.
12. A method of genotyping a genetic sample, comprising using at least one computer hardware processor to perform:
obtaining a directed acyclic graph (DAG) representing a reference sequence and genetic variation of the reference sequence, wherein a first path through the DAG represents a first sequence that includes a first allele and is stored as a first string of one or more symbols, a second path through the DAG represents a second sequence that includes a second allele and is stored as a second string of one or more symbols, and the second sequence represents a first genetic structural variation;
aligning at least some of a plurality of sequence reads from a genetic sample to the DAG, wherein aligning the at least some of the plurality of sequence reads to the DAG comprises:
comparing a read string of one or more symbols corresponding to a first sequence read of the plurality of sequence reads to the first and second strings storing the first and second sequences represented by the first and second paths, respectively;
determining a first overlap score corresponding to a first alignment, at a first position in the DAG, between the read string and the first and second strings, the first overlap score being determined based on, and indicative of, a degree of overlap between the read string and the first string and a degree of overlap between the read string and the second string;
determining a second overlap score corresponding to a second alignment, at a second position in the DAG, between the read string and the first and second strings, the second overlap score being determined based on, and indicative of, a degree of overlap between the read string and the first string and a degree of overlap between the read string and the second string; and
identifying one of the first and second overlap scores and therefrom aligning the first sequence read to the respective first or second position in the DAG;
determining, based on results of the aligning, whether the first sequence read aligns to the second sequence represented by the second path through the DAG; and
upon determining that the first sequence read aligns to the second sequence represented by the second path through the DAG, identifying a presence of the first genetic structural variation within the genetic sample.
13. The method of claim 12 , wherein the first genetic structural variation is between 1 Kb to 3 Mb in size.
14. 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 directed acyclic graph (DAG) representing a reference sequence of a species and genetic variation of the reference sequence, wherein a first path through the DAG represents a first sequence that includes a first allele and is stored as a first string of one or more symbols, a second path through the DAG represents a second sequence that includes a second allele and is stored as a second string of one or more symbols, and the second sequence represents a first genetic structural variation;
aligning at least some of a plurality of sequence reads to the DAG, wherein aligning the at least some of the plurality of sequence reads to the DAG comprises:
comparing a read string of one or more symbols corresponding to a first sequence read of the plurality of sequence reads to the first and second strings storing the first and second sequences represented by the first and second paths, respectively;
determining a first overlap score corresponding to a first alignment, at a first position in the DAG, between the read string and the first and second strings, the first overlap score being determined based on, and indicative of, a degree of overlap between the read string and the first string and a degree of overlap between the read string and the second string;
determining a second overlap score corresponding to a second alignment, at a second position in the DAG, between the read string and the first and second strings, the second overlap score being determined based on, and indicative of, a degree of overlap between the read string and the first string and a degree of overlap between the read string and the second string; and
identifying one of the first and second overlap scores and therefrom aligning the first sequence read to the respective first or second position in the DAG;
determining, based on results of the aligning, whether the first sequence read aligns to the second sequence represented by the second path through the DAG; and
upon determining that the first sequence read aligns to the second sequence represented by the second path through the DAG, identifying a presence of the first genetic structural variation within a genetic sample.
15. The system of claim 1 , wherein the first genetic structural variation is between 1 Kb to 3 Mb in size.
16. The at least one non-transitory computer-readable storage medium of claim 14 , wherein the first genetic structural variation is between 1 Kb to 3 Mb in size.
17. The method of claim 12 , further comprising determining, based on results of the aligning, a confidence value indicative of whether the first genetic structural variation is present within the genetic sample.
18. The at least one non-transitory computer-readable storage medium of claim 14 , further comprising determining, based on results of the aligning, a confidence value indicative of whether the first genetic structural variation is present within the genetic sample.
19. The system of claim 1 , wherein the first genetic structural variation comprises a deletion or an insertion.