Methods and systems for evaluating microsatellite instability status
Methods for evaluating microsatellite instability (MSI) analyze nucleic acid sequence reads corresponding to a plurality of marker regions for MSI. The marker regions may include long homopolymers and/or short tandem repeats (STRs). For a target homopolymer, a histogram of homopolymer signal values is calculated based on flow space signal measurements for the homopolymer region in the sequence reads. A score per marker based on features of the histogram of homopolymer signal values is determined for each marker region corresponding to the target homopolymers. For a target STR, the method includes calculating a histogram of repeat lengths for sequence reads corresponding to the marker region of the target STR. A score per STR marker is calculated based on features of the histogram of repeat lengths. A plurality of per marker scores may be combined to form a total MSI score for the sample.
1 . A method for determining control values for use in detecting microsatellite instability (MSI), comprising:
receiving a plurality of nucleic acid sequence reads corresponding to a plurality of marker regions for MSI, wherein the plurality marker regions include a plurality of target homopolymer regions, wherein a portion of the plurality of nucleic acid sequence reads correspond to the plurality of target homopolymer regions, the plurality of nucleic acid sequence reads resulting from a current sequencing run;
selecting a set of the homopolymer marker regions corresponding to the sequence reads having coverage values greater than or equal to a minimum coverage value;
calculating homopolymer signal values for each homopolymer in the selected set of the homopolymer marker regions;
calculating a first set of means of the homopolymer signal values for each homopolymer corresponding to the selected set;
calculating a second set of means of homopolymer signal values for each homopolymer corresponding to a second set of homopolymer signal values for same target homopolymer regions in a previous sequencing run;
determining a transformation function based on the first set of means and the second set of means; and
applying the transformation function to estimate a third set of means of homopolymer signal values for each homopolymer corresponding to a third set of homopolymer signal values, wherein third set of homopolymer signal values corresponds to sequence reads having coverage values less than the minimum coverage value, wherein the first set of means and the third set of means provide the control values for use in detecting MSI.
2 . The method of claim 1 , wherein the minimum coverage value is 50 homopolymer sequence reads.
3 . The method of claim 1 , wherein the step of calculating homopolymer signal values further comprises calculating a sum of M signal measurements corresponding to M nucleotide flows of a sequence of nucleotide flows having a same nucleotide type as a target homopolymer.
4 . The method of claim 1 , wherein homopolymer lengths in the selected set of the homopolymer marker regions are between 8 and 14.
5 . The method of claim 1 , wherein for the step of calculating a first set of means, the homopolymer signal values correspond to monomorphic homopolymer regions in a test sample for the current sequencing run.
6 . The method of claim 1 , wherein for the step of calculating a first set of means, the homopolymer signal values correspond to a control sample in the current sequencing run.
7 . The method of claim 1 , further comprising calculating a fourth set of means of homopolymer signal values for each homopolymer corresponding to a fourth set of homopolymer signal values for same target homopolymer regions in a test sample.
8 . The method of claim 7 , further comprising, for a given marker, calculating a difference between a corresponding mean of the fourth set of means and a corresponding control value.
9 . The method of claim 8 , further comprising, for the given marker, dividing the difference by a standard deviation of corresponding homopolymer signal values of the same target homopolymer region in the previous sequencing run.
10 . The method of claim 9 , wherein the second set of means of homopolymer signal values from the previous sequencing run are stored in a configuration file.
11 . A computer-readable media comprising machine-readable instructions that, when loaded in a machine-readable memory and executed by a processor, are configured to cause a system to perform a method for determining control values for use in detecting microsatellite instability (MSI), said method comprising:
receiving a plurality of nucleic acid sequence reads corresponding to a plurality of marker regions for MSI, wherein the plurality marker regions include a plurality of target homopolymer regions, wherein a portion of the plurality of nucleic acid sequence reads correspond to the plurality of target homopolymer regions, the plurality of nucleic acid sequence reads resulting from a current sequencing run;
selecting a set of the homopolymer marker regions corresponding to the sequence reads having coverage values greater than or equal to a minimum coverage value;
calculating homopolymer signal values for each homopolymer in the selected set of the homopolymer marker regions;
calculating a first set of means of the homopolymer signal values for each homopolymer corresponding to the selected set;
calculating a second set of means of homopolymer signal values for each homopolymer corresponding to a second set of homopolymer signal values for same target homopolymer regions in a previous sequencing run;
determining a transformation function based on the first set of means and the second set of means; and
applying the transformation function to estimate a third set of means of homopolymer signal values for each homopolymer corresponding to a third set of homopolymer signal values, wherein third set of homopolymer signal values corresponds to sequence reads having coverage values less than the minimum coverage value, wherein the first set of means and the third set of means provide the control values for use in detecting MSI.
12 . The computer-readable media of claim 11 , wherein the minimum coverage value is 50 homopolymer sequence reads.
13 . The computer-readable media of claim 11 , wherein the step of calculating homopolymer signal values further comprises calculating a sum of M signal measurements corresponding to M nucleotide flows of a sequence of nucleotide flows having a same nucleotide type as a target homopolymer.
14 . The computer-readable media of claim 11 , wherein homopolymer lengths in the selected set of the homopolymer marker regions are between 8 and 14.
15 . The computer-readable media of claim 11 , wherein for the step of calculating a first set of means, the homopolymer signal values correspond to monomorphic homopolymer regions in a test sample for the current sequencing run.
16 . The computer-readable media of claim 11 , wherein for the step of calculating a first set of means, the homopolymer signal values correspond to a control sample in the current sequencing run.
17 . The computer-readable media of claim 11 , further comprising calculating a fourth set of means of homopolymer signal values for each homopolymer corresponding to a fourth set of homopolymer signal values for same target homopolymer regions in a test sample.
18 . The computer-readable media of claim 17 , the method further comprising, for a given marker, calculating a difference between a corresponding mean of the fourth set of means and a corresponding control value.
19 . The computer-readable media of claim 18 , the method further comprising, for the given marker, dividing the difference by a standard deviation of corresponding homopolymer signal values of the same target homopolymer region in the previous sequencing run.
20 . The computer-readable media of claim 19 , wherein the second set of means of homopolymer signal values from the previous sequencing run are stored in a configuration file.