Systems and methods for determining structural variation and phasing using variant call data

Abstract

Systems and methods for determining structural variation and phasing using variant call data obtained from nucleic acid of a biological sample are provided. Sequence reads are obtained, each comprising a portion corresponding to a subset of the test nucleic acid and a portion encoding a barcode independent of the sequencing data. Bin information is obtained. Each bin represents a different portion of the sample nucleic acid. Each bin corresponds to a set of sequence reads in a plurality of sets of sequence reads formed from the sequence reads such that each sequence read in a respective set of sequence reads corresponds to a subset of the nucleic acid represented by the bin corresponding to the respective set. Binomial tests identify bin pairs having more sequence reads with the same barcode in common than expected by chance. Probabilistic models determine structural variation likelihood from the sequence reads of these bin pairs.

Description

CROSS REFERENCE TO RELATED APPLICATIONS[0001]This application is a continuation of U.S. patent application Ser. No. 15 / 019,928, entitled “Systems and Methods for Determining Structural Variation and Phasing Using Variant Call Data,” filed Feb. 9, 2016, now U.S. Pat. No. 10,854,315, which is hereby incorporated by reference in its entirety.[0002]This application claims priority to U.S. Provisional Patent Application 62 / 238,077, entitled “Systems and Methods for Determining Structural Variation Using Probabilistic Models,” filed Oct. 6, 2015, which is hereby incorporated by reference in its entirety.[0003]This application also claims priority to U.S. Provisional Patent Application 62 / 113,693, entitled “Systems and Methods for Determining Structural Variation,” filed Feb. 9, 2015, which is hereby incorporated by reference in its entirety.[0004]This application also claims priority to U.S. Provisional Patent Application 62 / 120,247, entitled “Systems and Methods for Implementing Linked R...

AI Technical Summary

Benefits of technology

[0062]In some embodiments, a blackout list is maintained, where the blackout list comprises a plurality of blackout regions of the test nucleic acid, the method further comprising eliminating a sequence read from the plurality of sequence reads when the first portion of the sequence read overlaps a blackout region in the plurality of blackout regions.

[0077]In some alternative embodiments, p∈{0, 1, −1} in which label “0” assigns a respective variant call in Ai;p to HO, label “1” assigns the respective variant call to H1, and label “−1” provides the advantageous possibility of denoting an error condition in the assignment of the variant call to a haplotype. This alternative embodiment takes into consideration that standard variant calling algorithms that are relied upon to provide calls at positions Ai;p between HO and H1 may, in fact, call such positions incorrectly on occasion. For instance, consider the case in which there are twenty sequence reads from the same sequenced nucleic acid fragment f, each with the same barcode 132, for a position i of fragment f and further suppose that conventional variant calling algorithms call position i heterozygous with seventeen of the sequence reads called HO at position i and the remaining three called H1 at position i. Further suppose that ground truth for position i from fragment f is, in fact homozygous HO meaning that the standard variant calling algorithms should have called position i HO for all twenty sequence reads. As such, the conventional haplotype assignment has miscalled three of the sequence reads. The disclosed alternative phasing embodiment p∈{0, 1, −1} advantageously accounts for the possibility of this form of error. In the example of the twenty sequence reads above where three of the twenty sequence reads are miscalled at position i, this position i in all twenty of the sequence reads is assigned −1 (H−1), the error state, when such error at position i is being sampled (to see if it provides a better phasing solution) by the disclosed phasing algorithms. The ability to selectively sample for this error state advantageously protects the phasing algorithm from error in the input data that arises, for instance due to error in the sequencing process, weak sequencing signal, and the like.

[0147]The present disclosure further provides a non-transitory computer readable storage medium storing one or more programs configured for execution by a computer, the one or more programs comprising instructions for carrying out any of the disclosed methods. Thus, these methods, systems, and non-transitory computer readable storage medium provide improved methods for detecting a structural variant in sequencing data of a test nucleic acid obtained from a biological sample.

Problems solved by technology

This limitation makes existing algorithms not applicable to targeted sequencing data, such as whole exome sequencing (WES) data.

Such a problem has been shown to be NP-hard.

Images