Determination of nucleic acid sequence concentrations
The present invention relates to a method of determining the concentration of a detected sequence in non-fragmented nucleic acids by applying a correction coefficient to the concentration of said detected sequence measured in fragmented nucleic acids.
1 . A method of determining the concentration of a detected sequence in non-fragmented nucleic acids comprising the following steps:
i. determining a length distribution (LD) of nucleic acids in a sample comprising fragmented nucleic acids, wherein the fragmented nucleic acids are derived from said non-fragmented nucleic acids;
ii. measuring the concentration of said detected sequence in said sample comprising fragmented nucleic acids with a measuring method, wherein measuring the concentration of the detected sequence comprises amplifying a sequence to be amplified comprising the detected sequence, and
iii. correcting the measured concentration of said detected sequence in the sample comprising fragmented nucleic acids with a correction coefficient to obtain the concentration of said detected sequence in non-fragmented nucleic acids, wherein the correction coefficient is based on the length distribution (LD) and at least one parameter of the measuring method, wherein the at least one parameter of the measuring method comprises length of the sequence to be amplified (La),
wherein the correction coefficient is determined by using the formula:
1
1
l
[
(
∑
i
=
L
∞
if
(
i
)
)
+
(
1
-
L
)
∑
i
=
L
∞
if
(
i
)
]
;
wherein: L is the length of the sequence to be amplified (La), L is (L a −n) with n being an integer from 1 to 15, or L is (r·La) with r being at least 0.75 and less than 1 with the proviso that r·La is always an integer value;
wherein f(i) Where, wherein L is the length of the sequence to be amplified (La);
wherein f(i) is the probability that a fragment in the sample has a length of i base pairs; and wherein
l
ˆ
=
∑
i
=
L
∞
if
(
i
)
is the average length of the nucleic acid fragments.
2 . The method of claim 1 , wherein the correcting is applied if at least 5% of the nucleic acids in the sample have a length shorter than the length of the sequence to be amplified.
3 . The method of claim 1 , wherein the length of the sequence to be amplified is longer than 40 bp and shorter than 200 bp.
4 . The method of claim 1 , wherein the length distribution (LD) of nucleic acid fragments in the sample is comprised in a range of from 25 bp to 350 bp.
5 . The method of claim 1 , wherein the method does not comprise a step of fragmenting the non-fragmented nucleic acids to generate the sample comprising fragmented nucleic acids.
6 . The method of claim 1 , wherein the method does not comprise assembly of the fragmented nucleic acids into a contiguous sequence.
7 . The method of claim 1 , wherein the method does not comprise obtaining or predicting the genetic coordinates of the fragmented nucleic acids.
8 . The method of claim 1 , wherein the method does not comprise obtaining the sequence of the fragmented nucleic acids.
9 . The method of claim 1 , wherein said measuring method is an isothermal quantitative nucleic acid amplification method.
10 . The method of claim 1 , wherein the measuring method is a non-isothermal quantitative nucleic acid amplification method.
11 . The method of claim 1 , wherein the method further comprises:
i. measuring the concentration of a second detected sequence in said sample comprising fragmented nucleic acids with a measuring method, wherein measuring the concentration of the second detected sequence comprises amplifying a second sequence to be amplified comprising the second detected sequence; and
ii. correcting the measured concentration of said second detected sequence in the fragmented nucleic acids with a correction coefficient to obtain the concentration of said second detected sequence in non-fragmented nucleic acids, wherein the correction coefficient is based on the length distribution (LD) and at least one parameter of the measuring method;
wherein the length of the sequence to be amplified is different from the length of the second sequence to be amplified.
12 . The method of claim 11 , wherein the measuring method comprises a multiplex amplification step.
13 . The method of claim 1 , wherein the sample comprising fragmented nucleic acids is a naturally fragmented sample.
14 . The method of claim 1 , wherein the sample comprising fragmented nucleic acids is an artificially fragmented sample.
15 . The method of claim 1 , wherein the correction coefficient is determined by using the formula:
1
1
l
^
[
(
∑
i
=
L
∞
if
(
i
)
)
+
(
1
-
L
)
∑
i
=
L
∞
f
(
i
)
]
;
wherein: L is the length of the sequence to be amplified (L a ), L is (L a −n) with n being an integer from 1 to 15, or L is (r·L a ) with r being at least 0.75 and less than 1 with the proviso that r·La is always an integer value;
wherein f(i) is the probability that a fragment in the sample has a length of i base pairs;
wherein
l
ˆ
=
∑
i
=
1
∞
if
(
i
)
is the average length of the nucleic acid fragments; and
wherein the correcting comprises multiplying the concentration measured in the sample comprising fragmented nucleic acids with the correction coefficient.
16 . The method of claim 15 , wherein L is the length of the sequence to be amplified (L a ).
17 . A system configured to determine the concentration of a detected sequence in non-fragmented nucleic acids comprising
i. an amplification module configured for measuring the concentration of said detected sequence in a sample comprising fragmented nucleic acids, said fragmented nucleic acids being derived from said non-fragmented nucleic acids;
ii. a module configured to compute a correction coefficient depending on a length distribution (LD) of nucleic acids in the sample and on at least one parameter of said measuring, wherein the at least one parameter includes the length of a sequence to be amplified (La) comprising the deleted sequence; and
iii. a module configured to compute the concentration of said detected sequence in non-fragmented nucleic acids with said correction coefficient,
wherein the correction coefficient is determined by using the formula:
1
1
l
[
(
∑
i
=
L
∞
if
(
i
)
)
+
(
1
-
L
)
∑
i
=
L
∞
if
(
i
)
]
;
wherein: L is the length of the sequence to be amplified (La), L is (L a −n) with n being an integer from 1 to 15, or L is (r·La) with r being at least 0.75 and less than 1 with the proviso that r·La is always an integer value;
wherein f(i) Where, wherein L is the length of the sequence to be amplified (La);
wherein f(i) is the probability that a fragment in the sample has a length of i base pairs; and wherein
l
ˆ
=
∑
i
=
L
∞
if
(
i
)
is the average length of the nucleic acid fragments.
18 . The system configured to determine the concentration of a detected sequence in non-fragmented nucleic acids according to claim 17 , amplification wherein the module configured for measuring the concentration of said detected sequence in said fragmented nucleic acids is an isothermal quantitative nucleic acid amplification module, selected from loop mediated isothermal amplification module and quantitative nucleic acid sequence-based amplification module.
19 . The system configured to determine the concentration of a detected sequence in non-fragmented nucleic acids according to claim 17 , amplification wherein the module configured for measuring the concentration of said detected sequence in said fragmented nucleic acids is a non-isothermal quantitative nucleic acid amplification module, selected from quantitative Polymerase Chain Reaction module, real time Polymerase Chain Reaction module, digital Polymerase Chain Reaction module, multiplex Polymerase Chain Reaction module and multiplex digital Polymerase Chain Reaction module.