A computer system for determination of paternity or maternity in a trio of a subject, comprising a processor operably connected to a memory storing a plurality of instructions, wherein the processor, upon processing the instructions, performs the following steps:
- (i) aligning at least two sequence reads from low-pass genome sequencing of genomic DNA of biological samples to a human genome reference according to an aligned chromosome and one or more genomic coordinates yielding a respective aligned sequence read;
- (ii) identifying a multiplicity of single-nucleotide variants (SNVs) in each respective aligned sequence read, wherein an SNV at each site has a mutant base type different from a base type at the corresponding site from the human genome reference;
- (iii) identifying a number of homozygous SNVs and a number of diploid heterozygous SNVs from the multiplicity of SNVs identified in step (ii), wherein
- a homozygous SNV is identified as an SNV where a percentage of sequence reads supporting the mutant base type different from a base type at a corresponding site from the human genome reference is 100%, and
- a diploid heterozygous SNV is identified as an SNV where a percentage of sequence reads supporting a mutant base type different from a base type at a corresponding site from the human genome reference is at least about 25% and no larger than 75%;
- and
- (iv) determining an inconsistent rate of base-type inheritance from the number of homozygous SNVs and the number of diploid heterozygous SNVs identified in step (iii) for an analytical model comprising a trio-based analysis comprising a proband and two parents, a mother and a father, wherein:
- for the trio-based analysis, (a) a number of loci that both the mother and the father are in homozygous alignment but with different genotypes in the i-th chromosome is denoted as Ai;
- (b) among the members denoted in (a), the number of SNVs that are homozygous in the proband but with different genotypes from the presumed father is denoted as pi for a paternity test and the number of SNVs that are homozygous in the proband but with different genotypes from the presumed mother is denoted as mi:
- for a maternity test;
- (c) the inconsistent rate of paternal inheritance αi in chromosome i is calculated based on the formula (1), and the inconsistent rate of maternal inheritance βi in chromosome i is calculated based on the formula (2);
- α
-
- i
- =
- pi
- Ai
- (
- 1
- )
- β
-
- i
- =
- mi
- Ai
- (
- 2
- )
- (d) the paternity is determined by using an average rate λpat across all autosomal chromosomes based on formula (3), and the maternity is determined by using an average rate λmat across all autosomal chromosomes based on formula (4);
- λ
- _
-
- pat
- =
- ∑
- i
- =
- 1
- n
- =
- 22
-
- α
-
- i
- 22
- (
- 3
- )
- λ
- _
-
- mat
- =
- ∑
- i
- =
- 1
- n
- =
- 22
-
- β
-
- i
- 22
- (
- 4
- )
- for providing the determination of paternity or maternity in a trio.