CN115206424B - Method, system, equipment and storage medium for identifying full sibling relationship - Google Patents
Method, system, equipment and storage medium for identifying full sibling relationship Download PDFInfo
- Publication number
- CN115206424B CN115206424B CN202211021134.2A CN202211021134A CN115206424B CN 115206424 B CN115206424 B CN 115206424B CN 202211021134 A CN202211021134 A CN 202211021134A CN 115206424 B CN115206424 B CN 115206424B
- Authority
- CN
- China
- Prior art keywords
- allele
- identified person
- homomorphic
- identified
- person
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 15
- 108700028369 Alleles Proteins 0.000 claims abstract description 89
- 108090000623 proteins and genes Proteins 0.000 claims description 7
- 238000004590 computer program Methods 0.000 claims description 6
- 238000004891 communication Methods 0.000 claims description 2
- 238000004364 calculation method Methods 0.000 abstract description 11
- 230000006870 function Effects 0.000 description 9
- 230000002068 genetic effect Effects 0.000 description 3
- 239000000463 material Substances 0.000 description 3
- 238000012408 PCR amplification Methods 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000001962 electrophoresis Methods 0.000 description 2
- 238000007400 DNA extraction Methods 0.000 description 1
- 238000007792 addition Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 239000008280 blood Substances 0.000 description 1
- 210000004369 blood Anatomy 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000005251 capillar electrophoresis Methods 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000004513 sizing Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16B—BIOINFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR GENETIC OR PROTEIN-RELATED DATA PROCESSING IN COMPUTATIONAL MOLECULAR BIOLOGY
- G16B20/00—ICT specially adapted for functional genomics or proteomics, e.g. genotype-phenotype associations
- G16B20/20—Allele or variant detection, e.g. single nucleotide polymorphism [SNP] detection
Landscapes
- Bioinformatics & Cheminformatics (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Genetics & Genomics (AREA)
- Biotechnology (AREA)
- Biophysics (AREA)
- Chemical & Material Sciences (AREA)
- Molecular Biology (AREA)
- Proteomics, Peptides & Aminoacids (AREA)
- Bioinformatics & Computational Biology (AREA)
- Analytical Chemistry (AREA)
- Evolutionary Biology (AREA)
- General Health & Medical Sciences (AREA)
- Medical Informatics (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Theoretical Computer Science (AREA)
- Measuring Or Testing Involving Enzymes Or Micro-Organisms (AREA)
Abstract
The invention discloses a method for identifying a holomorphic relationship, which comprises the following steps: obtaining genotype PQ of the identified person A and allele frequencies P and Q corresponding to the alleles P and Q respectively, obtaining genotype FG of the identified person B, and calculating the homomorphic relation index (FSI) of the identified person A and the identified person B by adopting a formula FSI = [0.25mpq +0.5 (aq + bnp)/y +0.25c ]/(mpq). The method simplifies a calculation formula of the holomorphic relation index (FSI) under the condition of 10 different types of genotype combinations in the technical specification of biological holomorphic identification (SF/T0117-2021) into 1 function model, greatly reduces the judgment workload and is convenient for batch processing.
Description
Technical Field
The invention relates to the technical field of forensic material evidence identification, in particular to a method, a system, equipment and a storage medium for identifying a holomorphic relationship.
Background
Full sibling, abbreviated FS, refers to multiple offspring individuals having the same biological father and mother (i.e., brother, sister, etc. of the same father and mother).
Identification of the holomorphic relationship, identification of full sizing relationship, refers to the process of determining whether the holomorphic relationship exists between two disputed individuals by detecting human genetic markers and analyzing according to genetic rules.
The forensic material evidence identification process mainly comprises the following steps: 1. case acceptance and sample collection (blood, fine spots or hair); DNA extraction (laboratory, obtaining the amount of DNA template required for subsequent detection); PCR amplification (PCR amplification instrument, using specific STR locus primers of STR kit combination to amplify or copy DNA template to reach the amount required by electrophoresis); 4. capillary electrophoresis (genetic analyzer); 5. typing the electrophoresis results (software); 6. calculating a full sibling relationship index; 7. forming an appraisal opinion; 8. the authentication script is written.
At present, the technical specification of the identification of the holomorphic relationship of the forensic physical evidence is "biological holomorphic identification technical specification" (SF/T0117-2021).
The problems existing in the prior art mainly comprise: in judicial identification practice, one case generally has at least 15 to 20 STR loci, and more 40 to 50 loci, and tens of combinations are present, and tens of calculations are required, and if tens of cases exist in one batch, hundreds of calculations are required. As the number of cases increases, the computational effort increases dramatically with the increase in STR loci detected. If the calculation is purely manual, the calculation is easy to make mistakes and is also very labor-consuming.
Disclosure of Invention
The embodiment of the invention provides a method for identifying a holomorphic relationship, which comprises the following steps:
obtaining the allele frequencies P and Q corresponding to the genotype PQ and the alleles P and Q of the identified human A respectively;
acquiring genotype FG of the identified human B;
calculating the holomorphic relation index FSI of the identified person A and the identified person B by adopting the formula (1):
FSI＝[0.25mpq+0.5(aq+bnp)/y+0.25c]/(mpq) (1)
calculating the value of each parameter of the formula (1) according to the following judgment conditions:
if the allele P of the identified person a is identical to the allele Q, i.e. P = Q, then m =1,n =0, otherwise m =2,n =1;
if allele F of identified person B is identical to allele G, i.e. F = G, then y =1, otherwise y =2;
a =1 if allele P of identified person a matches allele F or allele G of identified person B, i.e. P = F or P = G, otherwise a =0;
b =1 if allele Q of identified person a matches allele F or allele G of identified person B, i.e. Q = F or Q = G, otherwise B =0;
c =1 if allele P of identified person a matches allele F of identified person B, while allele Q of identified person a matches allele G of identified person B, i.e. P = F and Q = G, otherwise c =0;
and (3) judging whether the identified person A and the identified person B are in the homomorphic relation or not according to the homomorphic relation index FSI of the identified person A and the identified person B calculated according to the formula (1) and the biological homomorphic identification technical specification.
The invention also provides a system for identifying the full sibling relationship, which comprises the following steps:
the gene acquisition module is used for acquiring the allele frequencies P and Q corresponding to the genotype PQ and the alleles P and Q of the identified person A respectively; acquiring genotype FG of the identified human B;
the homomorphic relation index calculating module is used for calculating the homomorphic relation index FSI of the identified person A and the identified person B by adopting a formula (1):
FSI＝[0.25mpq+0.5(aq+bnp)/y+0.25c]/(mpq) (1)
calculating the value of each parameter of the formula (1) according to the following judgment conditions:
if the allele P of the identified person a is identical to the allele Q, i.e. P = Q, then m =1,n =0, otherwise m =2,n =1;
if allele F of identified person B is identical to allele G, i.e. F = G, then y =1, otherwise y =2;
a =1 if allele P of identified person a matches allele F or allele G of identified person B, i.e. P = F or P = G, otherwise a =0;
b =1 if allele Q of identified person a matches allele F or allele G of identified person B, i.e. Q = F or Q = G, otherwise B =0;
c =1 if allele P of identified person a matches allele F of identified person B, while allele Q of identified person a matches allele G of identified person B, i.e. P = F and Q = G, otherwise c =0;
and the homomorphic relation judging module is used for judging whether the homomorphic relation between the identified person A and the identified person B is established according to the homomorphic relation index FSI of the identified person A and the identified person B calculated by the formula (1) and the biological homomorphic identification technical specification.
The invention also provides a device for identifying and identifying the full sibling relationship, which adopts a system for identifying the full sibling relationship and also comprises an input device, an output device, a memory, a processor and a communication part.
The present invention also provides a storage medium for identifying a homomorphic relationship, on which a computer program is stored, the computer program being a homomorphic relationship identification program, the computer program, when executed by a processor, implementing the steps of a method for identifying a homomorphic relationship.
Compared with the prior art, the embodiment of the invention provides a method, a system, equipment and a storage medium for identifying the full sibling relationship, which have the following beneficial effects:
1. the calculation formula under the condition of 10 different genotype combinations in the whole sibling relation index calculation part in the technical specification of biological whole sibling identification (SF/T0117-2021) is simplified into 1 function model, and the judgment and calculation workload is greatly reduced.
2. The parameters of the function can be preset in the computer, and after the allele is input, the parameter index is automatically selected. The function can be realized by using common office software (such as EXCEL or WPS), and is convenient and easy to implement.
3. The function model developed by the patent can be used as a core algorithm, can be used for subsequently developing professional software by software developers, and can also be used for subsequently and autonomously developing mass processing software by forensic material evidence identification users through common office software (such as EXCEL or WPS).
Drawings
Fig. 1 is a flowchart of a method for identifying a homomorphic relationship according to an embodiment of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be obtained by a person skilled in the art without making any creative effort based on the embodiments in the present invention, belong to the protection scope of the present invention.
Referring to fig. 1, the present patent is optimized with respect to a functional model (computational method). For the DNA typing result obtained in the previous detection link in the forensic physical evidence identification process, through the setting of related parameters, the calculation formula (see the above content) under the condition of 10 different genotype combinations in the holomorphic relationship index calculation part in the technical specification of biological holomorphic identification (SF/T0117-2021) is simplified into 1 function model for the holomorphic relationship index calculation (FSI value). The method comprises the following specific steps:
description of symbols of the function model:
PQ, the genotype of the identified human A (the allele frequencies P and Q may be the same or different);
FG, the genotype of identified human B (corresponding to allele frequencies f and g, however, this gene frequency was not used directly in this functional model);
the expression of genotypes (PQ, FG) may be the same or different for the 2 alleles constituting the genotype of the same individual at a given locus, the numerical values (representing the number of repeats of the STR allele repeat unit). When the values are different, the smaller value is arranged in front. For example, an individual having 8 and 10 alleles at a locus is expressed as 8/10; when the allele combinations were 8 and 8, the expression was 8/8.
m, n, the genotype PQ of identified person a, if allele P is identical to Q (i.e., P = Q), then m =1,n =0, otherwise m =2,n =1;
y, genotype FG of identified human B, if allele F is identical to G (i.e. F = G), then y =1, otherwise y =2;
a, whether the allele P of the identified person a matches the allele F or the allele G of the identified person B, if there is one allele match (i.e. P = F or P = G), then a =1, otherwise a =0;
b, whether the allele Q of the identified person A matches the allele F or the allele G of the identified person B, if one allele is matched (i.e. Q = F or Q = G), then B =1, otherwise B =0;
c, genotype PQ of identified person A and genotype FG of identified person B, c =1 if allele P of identified person A matches allele F of identified person B (i.e. P = F), or c =0 otherwise.
The function model is reduced to 1:
FSI＝[0.25mpq+0.5(aq+bnp)/y+0.25c]/(mpq)
example 1: it is assumed that at a locus, the genotype of identified person A is 8/8 (homozygote, which can also be expressed as 8) and the genotype of identified person B is 8/10. At this locus, allele 8 had a gene frequency of 0.175 and allele 10 had a gene frequency of 0.221. Then, m =1, n =0, y =2, a =1, b =1, c =0, and the function is (0.25 × 1 × 0.175+0.5 (1 × 0.175+1 × 0.175)/2 +0.25 × 0)/(1 × 0.175) =1.6788.
Example 2: it is assumed that at a locus, the genotype of the identified person A is 8/10 (heterozygote) and the genotype of the identified person B is 8/10 (heterozygote). At this locus, allele 8 had a gene frequency of 0.175 and allele 10 had a gene frequency of 0.221. Then, m =2, n =1, y =2, a =1, b =1, c =1, and the function is (0.25 × 2 × 0.175 × 0.221+0.5 (1 × 0.175+1 × 0.221)/2 +0.25 × 1)/(2 × 0.175) =4.7620.
Although the embodiments of the present invention have been disclosed in the foregoing for illustrative purposes, those skilled in the art will appreciate that various modifications, additions and substitutions are possible, without departing from the scope and spirit of the invention as disclosed in the accompanying drawings.
Claims (4)
1. A method for identifying a homoeocellular relationship, comprising:
obtaining the allele frequencies P and Q corresponding to the genotype PQ and the alleles P and Q of the identified human A respectively;
acquiring genotype FG of the identified human B;
calculating the holomorphic relation index FSI of the identified person A and the identified person B by adopting the formula (1):
FSI＝[0.25mpq+0.5(aq+bnp)/y+0.25c]/(mpq) (1)
calculating the value of each parameter of the formula (1) according to the following judgment conditions:
if the allele P of the identified person a is identical to the allele Q, i.e. P = Q, then m =1,n =0, otherwise m =2,n =1;
if allele F of identified person B is identical to allele G, i.e. F = G, then y =1, otherwise y =2;
a =1 if allele P of identified person a matches allele F or allele G of identified person B, i.e. P = F or P = G, otherwise a =0;
b =1 if allele Q of identified person a matches allele F or allele G of identified person B, i.e. Q = F or Q = G, otherwise B =0;
c =1 if allele P of identified person a matches allele F of identified person B, while allele Q of identified person a matches allele G of identified person B, i.e. P = F and Q = G, otherwise c =0;
and (3) judging whether the identified person A and the identified person B are in the homomorphic relation or not according to the homomorphic relation index FSI of the identified person A and the identified person B calculated according to the formula (1) and the biological homomorphic identification technical specification.
2. A system for identifying a homonymy relationship, comprising:
the gene acquisition module is used for acquiring the genotype PQ of the identified person A and the allele frequencies P and Q corresponding to the alleles P and Q respectively; obtaining genotype FG of the identified person B;
the homomorphic relation index calculating module is used for calculating the homomorphic relation index FSI of the identified person A and the identified person B by adopting a formula (1):
FSI＝[0.25mpq+0.5(aq+bnp)/y+0.25c]/(mpq) (1)
calculating the value of each parameter of the formula (1) according to the following judgment conditions:
if the allele P of the identified person a is identical to the allele Q, i.e. P = Q, then m =1,n =0, otherwise m =2,n =1;
if allele F of identified person B is identical to allele G, i.e. F = G, then y =1, otherwise y =2;
a =1 if allele P of identified person a matches allele F or allele G of identified person B, i.e. P = F or P = G, otherwise a =0;
b =1 if allele Q of identified person a matches allele F or allele G of identified person B, i.e. Q = F or Q = G, otherwise B =0;
c =1 if allele P of identified person a matches allele F of identified person B, while allele Q of identified person a matches allele G of identified person B, i.e. P = F and Q = G, otherwise c =0;
and the homomorphic relation judging module is used for judging whether the homomorphic relation between the identified person A and the identified person B is established according to the homomorphic relation index FSI of the identified person A and the identified person B calculated by the formula (1) and the biological homomorphic identification technical specification.
3. A homonymy relationship authentication device using the system according to claim 2, further comprising an input device, an output device, a memory, a processor, and a communication section.
4. A storage medium for use in a homomorphic relationship authentication, having a computer program stored thereon, wherein the computer program is a homomorphic relationship authentication program, and wherein the computer program when executed by a processor implements the steps of a method for homomorphic relationship authentication as recited in claim 1.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211021134.2A CN115206424B (en) | 2022-08-24 | 2022-08-24 | Method, system, equipment and storage medium for identifying full sibling relationship |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211021134.2A CN115206424B (en) | 2022-08-24 | 2022-08-24 | Method, system, equipment and storage medium for identifying full sibling relationship |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115206424A CN115206424A (en) | 2022-10-18 |
| CN115206424B true CN115206424B (en) | 2023-04-07 |
Family
ID=83572334
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211021134.2A Expired - Fee Related CN115206424B (en) | 2022-08-24 | 2022-08-24 | Method, system, equipment and storage medium for identifying full sibling relationship |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN115206424B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105861668A (en) * | 2016-04-21 | 2016-08-17 | 昆明医科大学 | Forensic physical evidence paternity testing and individual identification parameter calculating method |
| CN108491691A (en) * | 2018-03-23 | 2018-09-04 | 河北医科大学 | Relationship iden- tification method and terminal device |
| CN109273046A (en) * | 2018-10-19 | 2019-01-25 | 上海晶准生物医药有限公司 | A kind of biology full sibs identification method based on probability statistics model |
| WO2022047413A2 (en) * | 2020-08-31 | 2022-03-03 | Etalon Diagnostics | Systems and methods for producing or identifying non-human animals with a predetermined phenotype or genotype |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR102446162B1 (en) * | 2020-07-21 | 2022-09-22 | 대한민국 | Multiplex PCR Kit for identifying human genotype profile using new combination of mini STRs and method for identifying human genotype profile using the same |
-
2022
- 2022-08-24 CN CN202211021134.2A patent/CN115206424B/en not_active Expired - Fee Related
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105861668A (en) * | 2016-04-21 | 2016-08-17 | 昆明医科大学 | Forensic physical evidence paternity testing and individual identification parameter calculating method |
| CN108491691A (en) * | 2018-03-23 | 2018-09-04 | 河北医科大学 | Relationship iden- tification method and terminal device |
| CN109273046A (en) * | 2018-10-19 | 2019-01-25 | 上海晶准生物医药有限公司 | A kind of biology full sibs identification method based on probability statistics model |
| WO2022047413A2 (en) * | 2020-08-31 | 2022-03-03 | Etalon Diagnostics | Systems and methods for producing or identifying non-human animals with a predetermined phenotype or genotype |
Non-Patent Citations (1)
| Title |
|---|
| 李梅等.复杂亲缘关系鉴定的研究进展.《法医学杂志》.2020,第691-698页. * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN115206424A (en) | 2022-10-18 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Cooper et al. | Quantitative estimates of sequence divergence for comparative analyses of mammalian genomes | |
| US20220101944A1 (en) | Methods for detecting copy-number variations in next-generation sequencing | |
| Tatsumoto et al. | Direct estimation of de novo mutation rates in a chimpanzee parent-offspring trio by ultra-deep whole genome sequencing | |
| CN112885408B (en) | Method and device for detecting SNP marker loci based on low-depth sequencing | |
| US20180018422A1 (en) | Systems and methods for nucleic acid-based identification | |
| CN111292802A (en) | Method, electronic device, and computer storage medium for detecting sudden change | |
| Stone | Joint genotyping on the fly: identifying variation among a sequenced panel of inbred lines | |
| CN110491446B (en) | Method and system for rapid mass SNP/Indel primer design | |
| CN107075565B (en) | Individual single nucleotide polymorphism site typing method and device | |
| CN118098348B (en) | Method and device for detecting genotype of hybrid parent, electronic equipment and medium | |
| KR101770962B1 (en) | A method and apparatus of providing information on a genomic sequence based personal marker | |
| CN115206424B (en) | Method, system, equipment and storage medium for identifying full sibling relationship | |
| US20220020449A1 (en) | Vector-based haplotype identification | |
| Rice et al. | Distinguishing multiple-merger from Kingman coalescence using two-site frequency spectra | |
| CN115762628A (en) | Detection method and detection device for gene progressive infiltration among biological populations | |
| Heyduk et al. | Targeted DNA region re-sequencing | |
| CN113981070B (en) | Method, device, equipment and storage medium for detecting embryo chromosome microdeletion | |
| CN115198024A (en) | Grandfather-grandson relationship identification method, system, equipment and medium for mother participation | |
| Martin et al. | Distribution of parental genome blocks in recombinant inbred lines | |
| CN116508105A (en) | Genomic marker interpolation based on haplotype blocks | |
| CN115206425B (en) | Triplet paternity testing method, system, equipment and storage medium | |
| CN115273976B (en) | Method, system, equipment and storage medium for identifying semi-sibling relation | |
| CN117953969B (en) | Mitochondrial disease prediction method and mitochondrial disease prediction system | |
| US11001880B2 (en) | Development of SNP islands and application of SNP islands in genomic analysis | |
| CN110570908A (en) | Sequencing sequence polymorphic identification method and device, storage medium and electronic equipment |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant | ||
| CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20230407 |
|
| CF01 | Termination of patent right due to non-payment of annual fee |