# Patent

## Title
Method and apparatus for mobile disaster victim identification

## Abstract

Mobile telecommunications or personal computing apparatus may be specially programmed for predicting whether an unknown biological specimen of an individual to be identified originates from or is related to a member of a particular family. The apparatus may comprise an input device including a barcode reader and a virtual touch screen keyboard whereby DNA profile or mass spectrometry test data may be entered among other data for an individual to be identified, for example, a missing person, victim, crime perpetrator or other unknown individual in hypothetical relationship to another individual or to a family pedigree of typed family members. The typed family pedigree in relation to the individual to be identified in hypothetical relation may be displayed on an output display of the apparatus including at least two different profile typing indicators for each displayed member of the pedigree.

## Background
This application is a continuation-in-part of U.S. application Ser. No. 13/633,940 filed Oct. 3, 2012, (now U.S. Pat. No. 8,775,097 issued Jul. 8, 2014), which is a continuation of U.S. application Ser. No. 12/684,539 filed Jan. 8, 2010, (now U.S. Pat. No. 8,301,392 issued Oct. 30, 2012), which claims the benefit of U.S. Provisional Application Ser. No. 61/193,927 filed Jan. 9, 2009 and is also a continuation-in part of U.S. application Ser. No. 13/621,311 filed Sep. 17, 2012, (now U.S. Pat. No. 8,788,215 issued Jul. 22, 2014), which is a continuation of U.S. application Ser. No. 11/467,834 filed Aug. 28, 2006, (now U.S. Pat. No. 8,271,201 issued Sep. 18, 2012), which claims the benefit of U.S. Provisional Application Ser. No. 60/836,941 filed Aug. 11, 2006, and is also a continuation-in-part of U.S. application Ser. No. 13/403,505, filed Feb. 23, 2012, (now U.S. Pat. No. 9,235,733 issued Jan. 12, 2016), and is also a continuation-in-part of U.S. application Ser. No. 12/823,284, filed Jun. 25, 2010, (now U.S. Pat. No. 8,392,418 issued Mar. 5, 2013), which claims the benefit of U.S. Provisional Application 61/220,248, filed Jun. 25, 2009, which applications are incorporated herein by reference as to their entire contents.


TECHNICAL FIELD
The technical field of embodiments of the present invention relate to on site collection of specimens, DNA profile data and related forensic evidence for an individual to be identified who may be a victim (or a crime perpetrator) via mobile telecommunications apparatus and to computing a match between the victim profile data and a collection of profile data for individuals that may belong to a family or to computing a ranking between the victim profile data and a plurality of families of the individual to be identified. The technical field also relates to a related method of searching for, ranking, and identifying individuals using DNA profiles and DNA database search techniques including automated decision support for identifying family members of the missing individuals for DNA typing and associating an unknown biological specimen of an individual to be identified with those of a family or a criminal.
BACKGROUND OF THE TECHNICAL ART
Identification of a biological sample using DNA profiles is an important task in forensic science. For example, a disaster caused by the terrorist attacks of Sep. 11, 2001 placed huge demands on forensic scientists to identify human remains from the collapsed World Trade Center buildings. Airplane crashes and natural disasters such as weather or earthquake related disasters are further examples of a need for forensic identification. Also, crime and accident scenes may require forensic identification of victims or perpetrators of a crime, for example, a suicide bomber terrorist or the terrorist's victims. In light of these demands, forensic scientists need more efficient and more accurate search methods to assist in identifying biological specimens by using DNA typing technologies to obtain DNA profile data.
Ideally, a forensic scientist obtains a DNA profile from a sample obtained from a personal effect of a missing person such as a toothbrush, razor, or comb, and searches for a match in a database containing DNA profiles from unknown biological specimens of a missing person or victim's remains. In theory, this approach can identify the missing person, but in practice, this approach breaks down when it encounters samples with partial profiles or when the reference origin of the personal effect cannot be obtained or verified. It is common to obtain incomplete DNA profiles from, for example, disaster areas due to harsh environmental conditions where the DNA integrity has diminished. This has forced forensic scientists to lower the match stringencies within database search engines, yielding potentially numerous false positives. In addition, incorrectly labeled personal effects can lead to inaccurate identifications. Biological specimens or forensic evidence of or associated with individuals may be subjected to additional analysis, for example, fingerprint analysis, iris scanning, biometric data collection, spectral data collection, microbody assemblage data collection and image collection in addition to DNA data collection which may provide individual identification or predict an individual's identity or predict properties of the individual or objects associated with the individual according to the principles set forth in U.S. Pat. No. 9,235,733 issued Jan. 12, 2016 (biometrics) and U.S. Pat. No. 8,392,418 issued Mar. 5, 2015 (predicting object properties).
When direct searching fails, identification using kinship analysis is often necessary. Kinship analysis comprises possibly narrowing the scope of a search by using any available DNA or non-DNA information to exclude unrelated specimens and then calculating genetic relatedness to at least one biological relative of a missing person or crime perpetrator. For example, the technology used for kinship analysis after the World Trade Center disaster of Sep. 11, 2001, relied on pair-wise comparison of a test DNA profile from an unknown biological specimen to a target DNA profile from a known biological relative, taking into account various familial relationships such as parent-child, sibling or half-sibling, and calculating the value of a function that indicates the likelihood or probability that the relationship is true (e.g., Cash et al., genecodesforensics.com/news/CashHoyleSutton.pdf, 2003). A likelihood ratio is commonly used, which indicates the likelihood that the given DNA profiles of the two samples would be obtained if they are related, relative to the likelihood or probability that these DNA profiles would be present if the individuals were unrelated. A measure of genetic similarity can also be used to indicate the likelihood that a relationship is true, such as the number of shared or matching alleles or a F.B.I. Combined DNA Indexing System (CODIS) low-stringency match. Such a measure can, for example, account for shared DNA alleles, loss of genetic information through degradation of the DNA, or the possibility of mutation of an allele. For any of these functions, the specimens are then independently sorted according to the function's value. When a likelihood function, such as probability, likelihood, or likelihood ratio is used, the specimens are sorted according to the calculated likelihood function value that assesses the likelihood that the DNA profile from an unknown biological specimen is related to the DNA profile from a biological relative. Unfortunately, this approach is cumbersome and imprecise for large cases, such as the World Trade Center disaster, because each search is for a specimen which is related to a single family member. A pair-wise comparison to the DNA profile of a single known relative can produce a large collection of candidate profiles. Multiple pair-wise comparisons can be performed and combined, but their combination does not always exclude genetically infeasible combinations of victim and family pedigree relationships. Human analysts must then sort, correlate, and analyze the matches, possibly manually with available DNA or meta data, which is a very labor intensive and time consuming process.
Software tools exist which allow the correlation of DNA match results from a single type of DNA profile, such as short tandem repeat (STR) profile, single nucleotide polymorphism (SNP) profile, mitochondrial DNA (mtDNA) profile, Y-STR DNA profile, and a mass spectrometry profile among others. Technologies can use all available DNA profile information involving a missing individual or an unknown biological specimen and his/her relatives to further enhance the ability to make an accurate identification.
U.S. Pat. No. 8,271,201 issued Sep. 18, 2012, of present inventor Dr. J. Douglas Birdwell et al., and corresponding to priority application U.S. Ser. No. 11/467,834 filed Aug. 28, 2006, (now U.S. Pat. No. 8,271,201 issued Sep. 18, 2012), describes a method for associating an unknown biological specimen with a family, for example, using a modified Elston Stewart algorithm. For example, the described method involves comparing a victim's DNA profile (an unknown specimen having been gathered, for example, from a disaster site), using one or more technologies to all available DNA profiles of selected family members and a pedigree that describes the relationship among the family members. U.S. Pat. No. 8,301,392 issued Oct. 30, 2012 of present inventor Dr. J. Douglas Birdwell et al., and corresponding to priority application U.S. Ser. No. 12/684,539 filed Jan. 8, 2010, (now U.S. Pat. No. 8,775,097 issued Jul. 8, 2015) describes an automated decision support algorithm for selecting family members of a missing person to type of available family members. U. S. Published Patent Application 2012/0148115 published Jun. 14, 2012 corresponding to priority application U.S. Ser. No. 13/403,505 filed Feb. 23, 2012 of present inventor Dr. J. Douglas Birdwell et al., (now U.S. Pat. No. 9,235,733 issued Jan. 12, 2016), describes mobile biometrics information collection and identification via a mobile telecommunications device providing iris scanning, fingerprint identification and other biometrics data collection features such that an identification of an individual and retrieval of that individual's data from a database may be obtained.
At a disaster, accident or crime site, a forensic data specialist may uncover and mark specimens taken from one or more individuals, bag the results together with completing a form identifying the sample(s) collected, and run or submit to be run, for example, individual profile tests to obtain selected DNA profiles. Various DNA technologies can be used including, but not limited to including, for example, STR (short tandem repeat), Y-STR (STR from the Y chromosome), SNP (single nucleotide polymorphisms), genetic sequence data and mass spectra. An assigned identifier, which may be encoded as a bar code that can be affixed to the sample collection bags or paperwork, may be used to correctly associate samples or other evidence with other samples or data. For example, a disaster or crime scene may have a number of unknown persons or “missing” persons for which samples may be typed and data collected. Presently, the collected data are physically returned or transmitted to a central location where software may be run to resolve a suspected DNA mixture, for example, per U.S. Pat. Nos. 7,162,372; 7,672,789; 7,860,661 or 8,140,271. In particular, the University of Tennessee has offered mixture deconvolution to police and fire and rescue organizations. This software may be utilized at a fixed location or adapted for use on portable intelligent communication devices without requiring centralized or “cloud” resources. Peak fitting algorithms analyze two dimensional graphic data representing DNA profile peaks and perform allele peak fitting and attribute extraction, for example, per U.S. Pat. No. 8,645,073 issued Feb. 4, 2014. An automated expert system also exists for performing automated expert analysis, for example, per U.S. Pat. Nos. 7,640,223; 7,624,087; 7,664,719; 7,840,519 or 7,945,526. U.S. Pat. No. 8,271,201 issued Sep. 18, 2012, as alluded to above, describes associating an unknown biological specimen with a family pedigree. U.S. Pat. No. 8,301,392 issued Oct. 30, 2012 describes an automated decision support tool for selecting family members of a family pedigree when, for example, one or more parents are unavailable for typing.
There remains a need in the art for mobile telecommunications apparatus and a method for on-site collection of DNA profile data and other forensic data such that a match may be performed locally or remotely to candidate individuals, for example, of a family pedigree or a criminal database or the like so that relatively instantaneous identification may be performed locally or transmitted to the forensic site. There is thus a need to be able to perform automated decision support for selecting available family members to type, to identify disaster victims on the scene of a disaster or crime or to collect and upload DNA profiles of family members available for typing, to fit DNA allele peaks, to perform missing persons kinship analysis, to perform DNA mixture analysis and to perform automated expert analysis locally or remotely such that results may be received at a forensic site or other location of interest using remote and local software and databases as necessary.
SUMMARY OF THE EMBODIMENTS
Several embodiments are discussed herein which provide methods of associating an unknown biological specimen with a particular family by use of intelligent telecommunications apparatus. By “unknown biological specimen” as used herein is intended a specimen of an individual to be identified, for example, by DNA profiling by known electrophoresis or mass spectrometry or other such processes and then comparing the results with similar results of processes conducted on specimens of known individuals related to the individual to be identified in a blood familial relationship. Automated decision support for identifying a family pedigree for an individual to be identified, such as a missing person, crime, disaster or accident victim or crime perpetrator, may aid forensic scientists to select enough and the correct family individuals related to such an individual to be identified for DNA typing. Once available family members are identified that are related to an individual to be identified such as a missing person, crime victim or perpetrator, there is provided mobile apparatus for identifying and displaying all such individuals including family members and the individual to be identified and typed, to compare with databases of known profiles and obtain matches to verify whether an individual is indeed the individual to be identified and so related to the family pedigree. The relative importance of different combinations of relatives to an individual to be identified and of different combinations of DNA typing technologies is used in automated decision support to replace guess work for selecting relatives to type and DNA typing technology processes to be used for typing a given selected relative. Automated decision support for selecting family relatives to a missing person maximizes the relative likelihood of obtaining a high probability of obtaining a correct match between an individual to be identified such as a missing person and a family that may be related to the individual to be identified. Mobile apparatus may be used to upload the selected family member profiles and those of individuals to be identified including missing persons and perpetrators or store and operate on the data locally at a site of interest. The output of the process may be a pedigree likelihood ratio (PLR) based on the totality of selected multiple relatives' DNA profiles or a mixture deconvolution or the like such that an unknown individual can become known, and may use one or more of the several DNA technologies mentioned above.
In one embodiment, there is provided a computer-implemented method of selecting and typing a subset of available family members for DNA typing according to a selected DNA typing technology to an individual to be identified to identify or exclude a typed unknown biological specimen for implementation on computer apparatus. The computer apparatus may be a mobile intelligent client telecommunications terminal including a processor, an input device coupled to the processor such as a bar code reader (one or two dimensional or both) and/or a camera and a virtual or physical keyboard such as a touch screen of a display, an output device coupled to the processor (such as a display), a memory for storing profile data obtained from said typed unknown biological specimen those of possible family members, the memory being coupled to the processor and a telecommunications interface for connection to a remote server or cloud service. (A bar code or a fingerprint may be read by a telephone camera which also may capture images of specimens). An example of a modern telephone is a Motorola brand Atrix 4G (fourth generation) mobile telephone among others. The computer-implemented method comprises: storing relationships in the memory between an individual to be identified and the available family members in a pedigree via the input device; either using the relative discriminating power of the pedigree with at least two combinations of the available family members to select a combination of available family members for DNA typing or specifying a pedigree consisting of a combination of available family members and their relationships; using at least one selected DNA typing technology of, for example, the typing technologies described above, for typing each available family member of the combination to obtain DNA profile data and storing the DNA profile data for the selected combination in the memory; using the pedigree and the stored DNA profile data of the selected combination to calculate a likelihood function value between the stored profile data obtained from the typed unknown biological specimen and the stored profile data for the selected combination of available family members via the processor (for example, one provided locally at a device, as a server on a network or via a network cloud service provider); and outputting a decision whether the typed unknown biological specimen originates from the individual to be identified and is related to the pedigree or to exclude the typed unknown biological specimen as unrelated to the pedigree. Introductory or subsequent processes may include but not be limited to a pedigree collection process (identifying available family members of a missing person's family pedigree) and ranking of the specimens against the at least one collected pedigree (collected profile data of selected available pedigree members) using the associated DNA profile data and a likelihood function.
In another embodiment, there is provided a computer-implemented method of identifying an unknown biological specimen as likely related to a family pedigree comprising at least first and second available family member relatives. The method is implemented on computer apparatus such as mobile telecommunications apparatus comprising a processor having memory, an input device such as a bar code reader and a touch screen of a display and a virtual or physical keyboard and an output device such as a display. The computer-implemented method comprises: determining a family pedigree for the unknown biological specimen of available family member relatives for DNA typing; using a rule base implemented as an application running on the computer apparatus to select a subset of available family members for DNA typing according to a selected DNA typing technology for each selected available family member; storing the family pedigree of the subset of available family members in the computer memory; obtaining and storing DNA typing data of the unknown biological specimen in the computer memory; obtaining and storing DNA typing data of the selected subset of available family members in the computer memory; using the family pedigree and said stored DNA typing data of the selected subset of available family members and the DNA typing data of the unknown biological specimen to calculate a likelihood function value between the stored DNA typing data obtained from the typed unknown biological specimen and the stored typing data for the selected subset of available family members via the processor; and outputting a decision whether the typed unknown biological specimen is likely to be or can be related to the pedigree or to exclude the typed unknown biological specimen as unrelated to the pedigree. For example, a likelihood ratio function (LR) may be preferably performed external to a field device, for example, at a server or via a networked cloud comprising a plurality of servers. Given sufficient computational capability in a field processor device, for example, a fully equipped personal computer, the likelihood ratio function may also be computed using a processor device in the field without having to resort to, for example, a cloud server.
In another embodiment, there is provided a computer-implemented method of supporting automated decision making for selecting relatives in a family related to a person for DNA typing to identify as related or exclude as unrelated a typed unknown biological specimen to the family. The method is implemented on computer apparatus or intelligent telecommunications apparatus comprising a processor having memory. The computer implemented method comprises: storing family pedigree information for a family composed of the relatives related to an individual to be identified from the unknown biological specimen in memory; storing DNA typing data of the typed unknown biological specimen in the memory; determining a degree of degradation for the unknown biological specimen, the unknown biological specimen being taken from the individual to be identified; if both parents of the person are available, then, select both parents for DNA typing and store DNA typing data for both parents according to a selected DNA typing technology in the memory; when one of the parents is unavailable and a child of the person is available, then, selecting the child for DNA typing along with the spouse of the individual; if one or both parents or a child is not available for DNA typing, selecting remaining available family relatives for typing according to a rule base; based on the degree of degradation of the remains and the rule base, obtaining and storing DNA typing data for selected family members according to the following hierarchy listed in order of most to least contribution of information obtainable about the individual to be identified if the following family member or family members are available for typing: a spouse and two children; child, parent and spouse; two children; one full sibling and one child; one child and a spouse; two full siblings; one parent or one child; and one full sibling; storing the DNA profile data for the determined available family members in the memory; and using the family pedigree and the stored DNA typing data of a selected subset of available family members and the DNA typing data of the unknown biological specimen of the missing person to calculate a likelihood function value between the stored DNA typing data obtained from the typed unknown biological specimen and the stored typing data for the selected subset of available family members via the processor; and outputting a decision whether the typed unknown biological specimen originates from the individual to be identified and the pedigree or to exclude the typed unknown biological specimen as unrelated to the pedigree.
A computer-implemented method of collecting DNA profile data and typing a subset of available family members for DNA profiling to an individual to be identified having an unknown biological specimen to identify or exclude a typed unknown biological specimen for implementation on telecommunications computer apparatus may be implemented on an intelligent telecommunications device such as a Blackberry, Motorola or HTC Android, Palm, an Apple iPad or other intelligent device for field forensic DNA data gathering, analysis and identification of the individual to be identified such as a missing person or perpetrator of a crime or other related use. The apparatus may comprise a processor, an input device coupled to the processor comprising a touch screen of a display and a bar rode reader, an output device coupled to the processor such as a display, a telecommunications interface coupled to the processor and a memory for storing profile data obtained from a typed unknown biological specimen coupled to the processor. The computer implemented method may comprise: a) storing relationships in said memory between said missing person and said available family members in a pedigree and storing DNA profile data for the missing person and said available family members via said input device via said bar code reader and touch screen; b) comparing DNA profile data for the individual to be identified to said available family members; and c) outputting a decision whether said typed unknown biological specimen originates from the individual and the family pedigree or to exclude the typed unknown biological specimen as unrelated to the pedigree. Automated decision support may be utilized to identify a family pedigree, the pedigree data and victim or perpetrator or mixture data input to the method and a missing persons/kinship and the missing person or perpetrator identified or least squares deconvolution resolution of a DNA mixture may result in the field. Moreover, following the principles of U.S. Published Patent Applications 2010/0332474 (predicting object properties) and 2012/0148115 (biometrics), an input device may comprise a fingerprint scanner, an iris scanner, a spectral data collector, a micro-body assemblage data collector and an image data collector for further assistance in identifying one of an individual, an object or predicting properties thereof.
Other embodiments comprise computer-readable media which store computer-executable instructions for performing any of the disclosed methods and computer apparatus for implementing the method.



BRIEF DESCRIPTION OF THE DRAWINGS
The present invention is illustrated by way of example, and not by way of limitation, in the figures of the accompanying drawings.
FIG. 1 is a first pedigree example showing family relative relationships to a missing relation.
FIG. 2 is a graph of the relative discriminatory power of a pedigree from a left-most logarithmic graph representing a low power level of three siblings in comparison to, for example, a right-most logarithmic graph showing a higher power level for one spouse and three children.
FIG. 3 is a second pedigree example showing blood relative relationships to a missing relation where both parents are deceased and no samples from them are available.
FIG. 4 is a flow diagram that illustrates an overview of one embodiment of a method of identifying an unknown biological specimen as likely related to a family comprising at least a first and second family member.
FIG. 5 is a flow diagram that illustrates an overview of one embodiment of a method of identifying an unknown biological specimen as likely related to a family pedigree comprising at least a first and second family member.
FIG. 6 is a flow diagram that illustrates an overview of one embodiment of a method of identifying at least two biological specimens as likely related to a family comprising at least a first and second family member by creating a directed graph.
FIG. 7 illustrates a block diagram of a hardware environment that may be used according to an illustrative embodiment of the invention.
FIG. 8A provides an exemplary introductory screen for a mobile disaster victim identification application for selecting an application to run on an intelligent telecommunications device such as automated decision support for selecting a family pedigree (ADS, not shown), matching an unknown biological specimen to a pedigree (MP/KN, not shown), least square deconvolution (LSD) mixture solution (LSD, not shown) and mobile disaster victim identification shown; FIG. 8B-8C provides an introductory menu for a selected mobile disaster victim identification application including specimen capture, pedigree entry, clear data, about and edit and also a mobile device trying to connect and connected to an external server (for example, to a remote server by wireless means from a disaster sire).
FIG. 9A-B provides an exemplary screen for selecting among Physical Profile, GPS Coordinates, DNA Profile and Back for a new specimen or to edit a specimen; FIG. 9C-9D provides an exemplary screen for entering physical profile shape data by scrolling through a list similar to one associated with a standard Interpol or other victim identification form silhouette sketch form part; FIG. 9E provides an exemplary screen for displaying GPS data acquired by an intelligent telecommunications device with GPS capability such that latitude and longitude may be automatically obtained, for example, for entry into a Place of Disaster field of a standard Interpol or other victim identification form; FIG. 9F shows, for example, a Place of Disaster field of an Interpol victim identification form; FIG. 9G shows a silhouette sketch form part of an Interpol victim identification form.
FIG. 10A provides exemplary laboratory or field kit analysis profile data collected for a father family member in the form of a single dimensional bar code uniquely identifying the father and two dimensional bar codes for each of STR, YSTR and mitochondrial DNA profile data, together with human-readable versions of the data represented by the bar codes, processing date and time, and analyst and (optionally) reviewer identifiers; FIG. 10B provides exemplary laboratory or field kit analysis profile data collected for a mother family member in the form of a single dimensional bar code uniquely identifying the mother and two dimensional bar codes for each of STR and mitochondrial DNA profile data, together with human-readable versions of the data represented by the bar codes, processing date and time, and analyst and (optionally) reviewer identifiers; FIG. 10C provides exemplary laboratory or field kit analysis profile data collected for a son family member in the form of a single dimensional bar code uniquely identifying the son and two dimensional bar codes for each of STR, YSTR and mitochondrial DNA profile data, together with human-readable versions of the data represented by the bar codes, processing date and time, and analyst and (optionally) reviewer identifiers.
FIG. 11A provides an exemplary screen for entering DNA profile data for a victim or a family member as per FIG. 10 and a checklist for data entered for that family member whereby each profile may be individually viewed; FIG. 11B provides a view of a one dimensional bar code being scanned for family member data; FIG. 11C provides a view of a two-dimensional bar code being scanned representing a DNA profile; FIG. 11D provides a view of an STR profile for a father family member entered as per FIG. 10A whereby a user may receive feedback that a two dimensional STR bar code has been correctly scanned and the corresponding data stored by comparison with the human-readable versions of the data shown in FIGS. 10A-10C; FIG. 11E provides corresponding data for a mitochondrial DNA profile; and FIG. 11F provides corresponding data for a YSTR profile.
FIG. 12A provides an exemplary screen for either entering a missing person (Enter Missing Person), disaster victim or other unknown biological specimen data or for browsing people (Browse People) already entered such as a mother, a father, a son, a daughter, a grandfather on a mother's side and so on for association with a pedigree and to comprise a family pedigree which may be an output of an automated decision support system and Back to take one back to a previous screen; FIG. 12B provides a screen for adding a family member such as a Father, Mother or Child and related data to that family member.
FIG. 13 provides an exemplary screen where, as shown, a father, a mother or a child of a member of a family pedigree may be edited (for example, Edit Father) or added (Add Child) and/or associated with a pedigree, and a user may edit the father, mother, child and the like, add personal information (Personal Information), collect a specimen (Collect Specimen), browse a specimen already entered (Browse Specimen) or view a pedigree chart (View Pedigree Chart) such as FIG. 1, 3 or 15A or Back to take one to a previous screen.
FIG. 14A provides an exemplary screen whereby may be entered for tentatively identifying an individual or a missing person by name gender and the like including a virtual touch screen keypad (or utilizing an actual keyboard in an alternative device) where, for example, ABC selects an alphabetic keyboard; FIG. 14B provides an exemplary screen for saving the data of FIG. 14A once input; FIG. 14C provides an exemplary screen for selecting a specimen; FIG. 14D provides a one dimensional identifier and two dimensional bar codes for each of STR, YSTR and mitochondrial DNA profile data for the missing person, together with human-readable versions of the data represented by the bar codes, processing data.
FIG. 15A provides a completed pedigree chart for mother, Emma, father, John, missing person John, Jr., son, John, III and daughter, Jane whereby collected profile data may be indicated by color; FIG. 15B shows that one may open details for a selected person; FIG. 15C provides a corresponding Select a Person screen related to FIG. 15B for selecting one of the family members of the pedigree of FIG. 15A including a further person Sarah Jane Smith, not shown in the pedigree chart of FIG. 15A; and FIG. 15C provides a corresponding Select a Person screen for an open details screen of FIG. 15B. Now the various embodiments and methods of operation will be described with reference to the drawings.



DETAILED DESCRIPTION OF THE EMBODIMENTS
Automated Decision Support for Determining which Blood Relations to Type and which DNA Typing Technology to Use
Automated decision support (ADS) may assist forensic scientists in kinship analysis to determine which blood relations or relatives to type from a set of available family members and relatives to a missing person. Determining who is missing is a first step, for example, in cases of a mass disaster, a common gravesite, war, terrorist attack, crime or natural or other phenomenon among other causes of the finding of missing human and associated remains or trying to locate a perpetrator of a crime. More than one individual may be missing in a family, (for example, if a family perishes in an airplane disaster). A next step is developing a family pedigree of relations with available DNA for developing DNA profiles. ADS may eliminate guesswork in determining a reasonable minimum number of relatives of a family to DNA type from a large pool of relatives available where DNA typing comprises obtaining at least one test DNA profile of a relative.
A mobile device which may be as small as an intelligent telephone or as large as a personal computer may provide local means for storing and/or uploading via a communications interface the selected family members, crime victim, perpetrator data and the like that results from ADS. A pedigree likelihood ratio (PLR) can result from an ADS method implemented, for example, on a personal computer via running a missing persons/kinship analysis optionally combined with a peak fitting analysis or a DNA mixture deconvolution. Challenges may arise in deciding relations to type from a large pool, for example, where most important relatives or remains are unavailable or their DNA is completely degraded, for example, where one or both parents are deceased or their artifacts from their personal use have been allowed to degrade over time (i.e. DNA samples left on articles of personal use). An ADS tool then may comprise a plurality of portions of a method and be embodied in the form of a single or plurality of computer apparatus or distributed systems. Thus, part or all of the ADS tool may be portable as may mobile disaster victim identification and used on location by inputting information about available relatives or kinship and DNA and unavailable DNA and then, from the PLR analysis results, determining a reasonable list of relatives to type rather than, for example, typing every relative available. Using ADS, one may obtain a sufficiently high degree of discrimination and match to a missing person's DNA from fewer samples without typing all available kin. This can reduce cost and the time required to identify or exclude, for example, human remains as related to other family members. ADS may determine a plurality of known family individuals to type and assess after a disaster has occurred. Mobile disaster victim identification algorithms running on intelligent devices permits STR, YSTR and mitochondrial DNA laboratory results to be collected and stored and/or uploaded for further analysis as they become available, whether from mobile laboratory processes and equipment which may be located, for example, at a disaster site or nearby, or fixed laboratory resources.
If there is a mixture of DNA as a human remain, the mixture may be resolved, for example, using the least square deconvolution (LSD) methods described in U.S. Pat. Nos. 7,162,372; 7,672,789; 7,860,661 and 8,140,271, incorporated by reference as to their entire subject matter for mixture resolution. The possible contributors to the mixture, identified for example using the methods described in U.S. Pat. Nos. 7,162,372; 7,672,789; 7,860,661 and 8,140,271, can then be assessed against the DNA profiles of the individuals identified for typing by the ADS in the same manner as with single source specimens. Any known typing methods may be employed including, for example, STR, Y-STR, mtDNA, mass spectrometry, genetic sequence, microarrays and SNP, and other methods known in the art.
Following the principles of U.S. Pat. No. 9,235,733 issued Jan. 12, 2016 (biometrics), biometric data may be collected at a disaster or crime scene utilizing an iris scanner, fingerprint scanner, a camera and related input devices 514 of a computer system 500 (FIG. 5) for assisting (in addition to DNA analysis) in the identification of human remains representing an individual. For example, a crime or disaster victim or crime perpetrator may be identified from biometric data with or without DNA profile testing from associated forensic evidence. Moreover, a mobile telecommunications device or personal computer running a biometrics application as described in this published application may run a combined or separate application for mobile disaster victim identification as described in the present application. According to FIG. 3B of U.S. Pat. No. 9,235,733 issued Jan. 12, 2016, and its description thereof, MobileDVI 114 is a selectable separate application for victim identification by DNA specimen identification that may be displayed on a touch screen. The MobileDVI application may be combined with and become one with a biometrics software application.
Following the principles of U.S. Pat. No. 8,392,418 issued Mar. 5, 2013 (predicting object properties), the clothing on a victim or of a crime perpetrator or other personal effect may comprise an object on which may be performed spectral data analysis, microbody assemblage data analysis and content-based image retrieval analysis, for example, to utilize, for example, a camera to collect an image of a tattoo or to collect data about or from an associated personal effect to predict an object property. For example, a micro-body assemblage found on a human specimen or personal effect may provide an indication of a distant crime scene (a crime victim was moved from the scene of a crime) or help tie a personal effect to the victim such as identification of dust particles found at the residence of the victim. Spectral data as defined in U.S. Pat. No. 8,392,418 issued Mar. 5, 2013 (predicting object properties) may include data indexed by frequency or time, or another coordinate system. Such spectral data, whether indexed or not, may be obtained by an input device (mass spectrometer, spectral analyzer and the like) or storage medium such as those used, as is known in the art, with mobile intelligent communications devices, computers (fixed such as desktop or mobile), digital cameras, and other devices that utilize storage media. These may be associated with a victim or perpetrator, communications may be intercepted during transmission or recorded from or during transmission or reception. These communications which may provide properties or be used to predict properties include e-mail and short message (text messaging) and voice, data, photograph or video communications.
Three DNA typing technologies are in current use in the forensic community: STR, Y-STR and mitochondrial DNA (mtDNA). In addition, mass spectrometry analysis is becoming increasingly utilized as a DNA profiling tool. See for example, U.S. Pat. Nos. 5,547,835; 6,090,558 and 6,602,662 which patents are incorporated by reference in their entirety as to the methods disclosed therein. A method and apparatus for allele peak fitting and attribute extraction from DNA sample data are described by Published U.S. Pat. No. 8,645,073 issued Feb. 4, 2014, and can be applied using any technology that creates data that exhibit peaks, and that have sizing or other reference standards or can otherwise be quantified in a standard manner. In FIGS. 1 and 3 herein, these three processes or technologies are represented as follows: STR is represented by a closed circle; mtDNA is represented by an open circle; and Y-STR is represented by an X. A fourth, SNP, is in use and a fifth, miniSTR is becoming more common. Emerging technologies include time of flight mass spectrum analysis, and genetic sequence measurement using traditional genetic analyzers and DNA amplification procedures or direct measurement using, for example, micro or nano technologies. All can have different rules of inheritance which influence the choice of relatives of missing persons to type. STR and SNP offer high powers of discrimination, the latter if array technologies are used. ADS will be exemplified by analysis of first and second pedigree examples such that an algorithm for ADS results.
Referring first to FIG. 1, there is shown a first pedigree example. The missing person, victim or crime perpetrator is indicated with a shaded box where a box represents a male and a circle a female. For example only, it is assumed in FIG. 1 that the missing person 2 is male and that an unknown specimen that may have originated front missing person 2 has been typed using all three DNA typing technologies. A diagonal line represents that the missing person 2 or other relative shown in the pedigree diagram is deceased or DNA is otherwise unavailable. The horizontal line connecting sibling 8 and sibling 10 to missing person 2 represents kinship at the same level while the vertical line connecting to parents 14 and 16 indicates the next higher level of kinship and so on. With reference to FIG. 1, missing person 2 is male and married to female 4. Female 6 represents a female child of missing person male 2 and female 4. Parent 14 has a female child 12 by another spouse than spouse 16 who is deceased, indicated by the diagonal line. So child 12 is a half-sister to missing person 2. Unavailable mother 16 has a sister 18, who is missing person 2's aunt. From FIG. 1 (and FIG. 3 using the discriminatory power calculations and graph, for example, of FIG. 2), the following rules of a rule base for ADS may be determined.
A very simple rule is utilized as a first-order approach when attempting to identify a person from remains of a missing person, victim or missing crime perpetrator. The simple rule is to type (obtain DNA profiles for) both parents if they are available. Typing both parents provides information across all DNA types and an STR profile alone (for example, using the thirteen core loci of the Federal Bureau of Investigation database) provides sufficient discrimination for disasters of moderate size such as an airplane crash when the DNA is not severely degraded. There are corollaries to the general rule. If both parents of the missing person are available and have complete typed profiles, no additional information can be gained from the parents' ancestors or the missing person's siblings. Similarly, if one has complete typed profiles of a child of the missing person, there is no reason to obtain specimens from that child's offspring. Similar properties hold for other relatives, that is, no child provides more information than a parent where the missing person is on the same or higher kinship level than the parent. Referring to FIG. 1, this simple rule of being able to type both parents and so limit the forensic scientist's work is not applicable because only one parent, father 14, is available. A more complicated case where DNA for both parents are unavailable will be discussed in connection with FIG. 3. In FIG. 1, one parent 16, the mother, is not available. For example, mother 16 may also be a victim, for example, in the case of a mass grave or has died with their male child 2 in an airplane crash and is missing as well or otherwise unavailable. Decomposition, confounding influences such as sea or other water and heat or environmental temperature where the specimen is found (typically, weather related), for example, in the case of an airplane crash or fire, degrade DNA. In such situations, it is not realistic to expect complete typing of the STR core loci. If, however, one or both parents' typed profiles are unavailable, or sufficiently degraded, then it becomes necessary to obtain profiles from other relatives, and the correct choice of relatives becomes complex. Consequently, the living male parent, father 14, is selected to be typed indicated by the dots in the square, starting with STR typing, as will be discussed further herein.
In FIG. 1, there is a living female offspring, daughter 6, who may be typed. In the case where only one female offspring, out of all immediate family members, can be typed, it is necessary to look to mother 4, grandparents, aunts 18, uncles, nieces, nephews or grandchildren for further increasing the probability or discriminatory power of a match. Two, three or four individuals may need to be typed in order to obtain a DNA pedigree of sufficiently high discriminatory power for missing person identification. Especially with the availability of multiple DNA typing technologies, for example, STR, Y-STR and mtDNA, the decision of whom to type and which DNA typing technology to use is nontrivial. Even though the mathematics, though convolved, is fairly straightforward and well-known (see, for example, R. C. Elston and J. Stewart. A General Model for the Genetic Analysis of Pedigree Data. Human Heredity 21: 523-542 (1971); C. Cannings, E. A. Thompson, and H. H. Skolnick. The Recursive Derivation of Likelihoods on Complex Pedigrees. Advances in Applied Probability, Vol 8, No 4. (December 1976), pp 622-625 and E. Lander, and P Green. Construction of Multilocus Genetic Linkage Maps in Humans. PNAS, Vol. 84, No. 8 (Apr. 15, 1987), pp. 2363-2367), it is not realistic for the average forensic scientist to make a good decision under such circumstances. Consequently, an ADS algorithm and tool as described herein provides an easy-to-use platform for the organization of information about relatives to a missing person, for example, missing person 2 and developing a recommended set of relatives' DNA specimens to type. The ADS algorithm may be run on a portable personal computer. If the ADS tool or portions thereof is implemented on a single personal computer, the portable personal computer may access remote databases and communicate with remote computer hardware via a communications interface and receive input from typing performed on selected relations to missing persons performed remotely. Moreover, the ADS tool is not limited to a personal computer embodiment but may be implemented in firmware, software and hardware on larger or smaller known computer apparatus or plural computer apparatus as will be described further herein in connection with a discussion of FIG. 7.
The rule base and algorithm will now be further defined before completion of the discussion of FIG. 1. Identification of the remains of a disaster victim or missing person MP 2 using DNA evidence is best performed using a Bayesian approach to calculate a Likelihood Ratio (LR) to gauge the relative strength of two hypotheses (1) there exists a kinship of the stipulated type to known and typed relatives, and (2) there is no relationship between the MP's remains and that of the putative relatives. The relationships between the relatives and the missing person can be captured in a pedigree, to which data that describe specimens, DNA profiles, and metadata can be associated. A fundamental question that affects cost and efficacy of the identification process is the relative discriminatory power of different combinations of family members (measured, for example, by likelihood ratios (LR's), odds, or posterior probabilities). For example, the question may be asked whether the DNA profiles of one parent and one sibling is better able to establish identity than those of two off-springs. For example, the question may be asked whether the combination of a spouse and two off-springs is better than a parent and 2 siblings. Answers to these types of questions require the knowledge of the distributions of the LR of pedigrees having various combinations of family members and the missing person. Such knowledge is best embedded in an automated decision support (ADS) environment due to the complexity of the underlying mathematics. The problem, for example, is studied in FIG. 2, whose underlying mathematical calculations for different combinations of available family members and results may be embedded in the automated decision support algorithm for implementation, as suggested above, as a tool on a personal computer.
Quantitative knowledge of the various likelihood ratio distributions of pedigree families with respect to the missing person, having different combinations of members may be required for this automated tool. It is recognized, however, that such knowledge may be imprecise and is usually obtained from relatively small populations. It is permissible that the quantitative knowledge of likelihood ratio distributions, or of likelihoods, frequencies, or probabilities be estimated. Furthermore, the distributions may be specified by one or more statistics such as a mean, median, or a mean and a variance. The tool becomes a computer-implemented method to decide which minimum combination is sufficient to identity remains with high confidence. A threshold pedigree LR value or sliding value, for example, may be based on the confidence in the typing results of a missing person 2 whose DNA may be degraded. The threshold may be established from the distribution for each combination type of pedigree family members that can be used to evaluate the LR that the DNA profile of remains must satisfy to meet desired error criteria. The incremental value of an additional family member, if brought into the pedigree, in aiding the identification can also be established, a priori, before his/her sample is acquired. Identification by pedigree likelihood ratio based on STR DNA can be further strengthened by the use of mtDNA anchor Y-STR DNA. If a disaster victim's (DV's) remains are badly degraded or fragmented, then the use of mtDNA and/or Y-STR DNA, depending on the gender, becomes necessary in matching to a suitably linked family member. Even when the remains are not degraded, the use of mtDNA and/or Y-STR DNA typing, or other typing technologies, will increase the discriminatory power over that of using STR DNA alone, for example, when selected relatives to a missing person are more distant and the number of available relatives is limited.
Prior to release of the FBI's CODIS software, version 6, pair-wise or parent-offspring trio comparison tools have been used almost exclusively in mass fatality incidents. It is generally accepted that, except for parent-offspring comparisons, a significant number of fortuitous and false positive matches can occur in pair-wise comparisons based upon either allele sharing or LRs, resulting in weak leads requiring laborious forensic analyst follow up. Parents-offspring trios require that two of the three trio members be known, a rather restrictive requirement. When genetic data are available for multiple pedigree members, pair-wise comparison is ineffective and wasteful. Exclusions can be made using two typed family members that pair-wise comparisons will not detect; further, LR based on multiple family members is predictably much more discriminatory than pair-wise LR in identifying a missing person. Some commercial software exists to perform an analysis for identification of family relationship in paternity, immigration, or inheritance disputes based on one or more typed family members, but not to provide guidance on selection of relatives for genetic typing. The present automated decision support algorithm may compute the pedigree LR of an arbitrarily specified pedigree mix and profiles.
Referring again to FIG. 1, a single parent example is presented. Referring to FIG. 3, an example where both parents are unavailable is manually analyzed to demonstrate the potential benefit and the technical soundness of ADS and to continue the development of a rule base for the ADS tool. The ADS tool implemented, for example, on a personal computer will automate the manual process that is discussed below and requires no undue experimentation for one of ordinary skill in the art to implement. One innovation is the utilization of computed distributions of pedigree LRs to recommend the acquisition of the data most valuable in the identification process.
As already introduced above, FIG. 1 shows an example pedigree of the family of a missing person 2. Remains of the missing person 2 have been recovered; the question is whether these remains correspond to the presumed identity of the missing person 2. The missing person 2 is male, has a spouse 4 and child 6, a brother 10 and sister 8, a half-sister 12, a father 14, and a maternal aunt 18. The mother 16 is deceased, and no DNA samples of mother 16 are available. A sample of the remains of missing person 2 has been typed using all three technologies (STR, mtDNA, and Y-STR), but the remains are degraded so only partial profiles are available for each technology. The degree of degradation being high, it is determined that STR, Y-STR and mtDNA test data be obtained and stored, for example, in a personal computer for the missing person 2 remains. As will be further described herein, non-DNA forensic data may also be stored to improve the probability of a match if such data are known for missing person 2. The question answered by an automated decision support (ADS) algorithm is which DNA samples from relatives should be collected and what typing should be performed, selected, for example, the three well known technologies.
Turning first to STR data, both the child 6 and the father 14 share one allele at each locus with the male missing person 2. The spouse's 4 STR profile, together with the child's 6, can be used to identify, for some loci, alleles that the victim or missing person 2 must have (but which may not be present in a profile because of the degraded condition of the remains of missing person 2). Therefore, the STR profiles of the child 6, spouse 4, and father 14 should be obtained, indicated by the depicted closed circle. Failing this, knowledge of the STR profiles of the siblings 8, 10 is next most valuable. Since the sister 8 is missing a Y chromosome, the brother 10 is preferred over the sister.
FIG. 2 shows an example of the relative discriminating power of pedigrees having three family members of various relations. The combination of three siblings of a missing family member is the least powerful in this example, with the combination of a parent and two siblings being better. Then follows, a parent and a child, three children, one spouse and two children and the best combination in the example is one spouse and three children. As a second example, Table 1, below provides a list, rank ordered from most to least, of the relative contributions of information for various combinations of typed near relatives which may be stored in computer memory.







TABLE 1





Relative contribution of information (from most to least) to identification


of a missing person's remains using STR technology by various


combinations of typed near relatives.

















1. Both Parents



2. Two Children + Spouse



3. One Child + One Parent + Spouse



4. Two Children



5. One Full Sibling + One Child



6. One Child + Spouse



7. *Two Full Siblings



8. One Parent/Child



9. *One Full Sibling







*Note that a parent or a child in a pedigree acts as a filter to eliminate infeasible candidates, whereas siblings (without parent information) do not.






Continuing with the analysis of FIG. 1, a Y-STR profile, indicated as an “X,” can be obtained from either the father 14 or the brother 10. These have almost equal value, with the brother's 10 profile slightly less valuable. The reason the brother's 10 is less valuable than the father 14 is the possibility of mutations, which occur during meiosis. One meiosis event separates the father 14 and missing person 2, while two events separate siblings such as brother 10 from missing person 2. There are no other available male relatives in the example of FIG. 1.
The missing person's 2 mtDNA profile is the same, barring mutations, as those of his siblings 8, 10, his mother 16 (who is unavailable), and his maternal aunt 18. Since the mother is not available and a relevant Y-STR profile can be obtained from his brother 10, it makes sense to obtain a sample from the brother 10 and determine Y-STR and mtDNA profiles along with the STR profile. The sister's 8 DNA is relatively less valuable because she does not have a Y chromosome. Failing this, the aunt's 18 mtDNA can be typed. Mutations are again a consideration when comparing two genetically linked mtDNA or Y-STR profiles where the aunt is removed in kinship from the missing person.
A symbol for each individual is marked in FIG. 1 to indicate the genetic material that should be typed. Thus, the result of a manual analysis of whom to type and how are presented. The sister's 8 DNA could be typed if the brother 10 is unavailable and can in any case provide additional information about the MP's STR profile. Likewise, the half sister 12 and aunt 18 can provide STR data but the relevant information content is significantly less than information provided by the siblings' 8, 10 DNA.
A pedigree likelihood ratio can be computed for each DNA typing technology using the typed genetic data from the chosen family members 4, 6, 10, 14 and 18, comparing the probability that the remains of missing person 2 correspond to the missing person's location in the pedigree (hypothesis 1 above) to the probability that the remains correspond to an unrelated person (hypothesis 2). Computational methods are known in the art. These methods are adapted in ADS to calculate the pedigree LR for arbitrary pedigrees. Since the profile information obtained from each different DNA typing technology can be assumed independent of information from the other typing technologies, a joint likelihood ratio (JLR) can be computed as the product of the likelihood ratios (LR) for each typing technology. A large value indicates the likelihood that the remains correspond to the missing person 2 (hypothesis 1). A value of zero indicates that the pedigree data, combined with the typed genetic data obtained from the remains, are not consistent and that, unless there is an error or incorrect relationship in the pedigree, the remains are not related to the family's pedigree and can be excluded.
The missing person case represented by the second example pedigree, shown in FIG. 3, is more challenging because DNA is not available from either parent, father 24 or mother 26, of missing person 22. (DNA typing of the remains of missing person 22 has been taken using all three technologies as per FIG. 1.) The missing person 22 has a living sister 30, who has children 32, 34 (a boy 34 and a girl 32), a half-brother 28 on his mother's side, and a half-brother 38 and half-sister 36 on his father's side. Aunts 40, 42 and an uncle 44 are also living.
A sample from the sibling sister 30 can provide both STR and mtDNA profiles, but not Y-STR. If the sister's 30 profiles are available, her children's 32, 34 DNA profiles add no new information according to the rule defined above in discussion of FIG. 1. STR information about the parents' genomes can be obtained from the half-brothers 28, 38 (mother's and father's sides) and half sister 36 (father's side only), and from aunts 40, 42 on both sides and the mother's brother 44. DNA profile data for most of these relatives should be used because all of the relationships between the missing person 2 and living individuals are less direct and therefore provide poorer discrimination than was seen in the first example. STR typing is shown for all selected relatives 28, 30, 36, 38, 40, 42 and 44. Typing by mtDNA is shown for relatives 28, 30, 42 and 44, that is, any of the mother's available relatives. Y-STR typing is shown for only half-brother 38 because he has a Y chromosome from his father's side, missing from half-sisters and more relevant than half-brother 28 on his mother's side. As in the first example of FIG. 1, the possibility of mutations should be considered, but, in this example of FIG. 3, the number of meiosis events between the missing person 22 and his living relatives, except for the sister's children 32, 34, is the same in all cases. The symbol for each individual is marked in FIG. 3 to indicate the presence of genetic material that should be typed and how. As before, a pedigree likelihood ratio (PLR), and comparing the probability that the remains of the missing person 22 correspond to the missing person's location in the pedigree to the probability that the remains correspond to an unrelated person, can be computed for each technology using existing techniques and, if desired, combined to generate joint likelihood ratios. Again, the hypothesis tested is whether the missing person 22 whose remains have been found is more probably or likely than not related to the available family members of FIG. 3.
Methods of Predicting Whether an Unknown Specimen is Related by Kinship
Referring to FIGS. 5-7, there will now be described several embodiments of methods of predicting whether an unknown biological specimen from an individual originates from a member of a particular family. The individual can be deceased or can be living but have an unknown or questioned genetic relationship to the family. Methods of an embodiment may use as much available DNA profile information as possible but may be reasonably limited as discussed above with respect to automated decision support (ADS). The methods permit a forensic investigator to identify an unknown biological specimen by ranking test DNA profiles from unknown biological specimens according to calculated probability, likelihood, likelihood ratio, or a likelihood function values that the test DNA profiles are related to target DNA profiles from family members of the missing individual. If an unknown biological specimen has high likelihood function values when compared to more than one target family member, this can indicate an increased likelihood (or probability) of a relationship between the unknown biological specimen and the family.
Prediction methods typically are used with humans but can also be used for any type of individuals that can be grouped into a family. A “family” according to the invention is a group of individuals who are genetically related to a missing individual. A mother, father, grandfather, grandmother, full or half siblings are examples of members of a family, as are cousins, aunts, and nephews, half-brothers and half-sisters and so on. The corresponding missing individual is also a member of the family (the hypothesis to be demonstrated). Genetic relationships exist within a family according to the Mendelian laws of inheritance. When one calculates a likelihood or likelihood ratio (LR), one takes into account known or stipulated genetic relationships such as parent-child, siblings, etc. For example, half of the genetic information contained in a father's nuclear DNA (other than the Y chromosome) will be inherited by a child. Only 1⁄4 of the genetic information contained in a grandfather's nuclear DNA (other than the Y chromosome) will be inherited by his grandchild. All of a mother's mitochondrial DNA (mtDNA) will be inherited by her child. All of a father's nuclear DNA from the Y chromosome will be inherited only by his son. The inheritance of both nuclear and mitochondrial DNA may be influenced by mutations as suggested above. The present prediction methods exploit these genetic relationships to identify an unknown biological specimen as having a probable or likely genetic relationship to the typed genotypes of members of a family.
Suitable families include those of vertebrates, most particularly mammals, including primates (e.g., gorillas, chimpanzees, baboons, squirrel monkeys, humans), companion animals (e.g., cats, rabbits, dogs, horses), farm animals (e.g., cows, sheep, swine, goats, horses), wild animals (e.g., lions, tigers, elephants), and research animals (e.g., cats, dogs, guinea pigs, rabbits, sheep, goats, swine, chimpanzees, mice, rats, and baboons). However, the methods are not limited to species of vertebrates and may be used whenever Mendelian laws of inheritance are known. For example, the methods may be used for families of a plant or other animal species.
The methods' power is based in part on comparison of DNA profiles from unknown biological specimens to DNA profiles of more than one family member, which significantly increases the methods' predictive ability. A “DNA profile” according to the embodiments comprises one or more DNA markers which together uniquely characterize(s) an individual with high probability and which contain inherited genetic information. A DNA profile can be from one or more of the following DNA types, which correspond to both the source of the DNA and the laboratory process used to obtain the information: single nucleotide polymorphism (SNP), autosomal short tandem repeat (STR), mitochondrial DNA (mtDNA), mini-STR, or STR DNA from the Y chromosome (Y-STR), and mass spectrum analysis. In the future, additional DNA types will most likely be developed, and since the methods disclosed herein rely upon the Mendelian laws of inheritance and the mathematics of probability rather than a specific DNA typing technology, these methods will also apply to any newly discovered DNA types. A DNA profile typically contains one or more DNA markers (e.g., 1, 2, 3, 4, 5, 6, 7, 8, 9, 10, 11, 12, 13, 14, 15, or more), which can be of various types. For example, polymorphic STR loci which can be included in a DNA profile include D3S1358, vWA, FGA, D8S1179, D21S11, D18S51, D5S818, D13S317, D7S820, D16S539, THO1, TPOX, CSF1PO, Penta D, Penta E, and D19S433.
DNA profiles for use in predictive methods of the embodiments are obtained from biological specimens. A “test DNA profile” according to the embodiments is a DNA profile obtained from an unknown biological specimen or an individual with a questioned relationship, such as missing person 2 of FIG. 1 and missing person 22 of FIG. 3. A “target DNA profile” as defined herein is a DNA profile obtained from a biological specimen of a personal effect of a missing person or a known family member. A “biological specimen” as defined herein is a sample from which DNA can be obtained. Suitable samples include, but are not limited to, nucleated blood cells, bones or bone fragments, skin cells, hair, saliva, cells obtained from a cheek swab, and DNA remaining on a personal effect, such as a toothbrush, bedding, a razor, a glass used for drinking, a cigarette butt, or a hairbrush. An unknown biological specimen can be derived from living tissue (e.g., a biopsy) or from deceased tissue (e.g., remains), bone fragments, hair, or fingernail scrapings. As introduced above, following the principles of U.S. Pat. No. 8,392,418 issued Mar. 5, 2013 (predicting object properties), the clothing on a victim or of a crime perpetrator or other personal effect such as a toothbrush may comprise an object on which may be performed spectral data analysis, microbody assemblage data analysis and content-based image retrieval analysis, for example, to utilize, for example, a camera to collect an image of a human specimen containing a tattoo or to collect data of from an associated personal effect to predict an object property.
One can extract DNA from an unknown biological specimen by using any DNA extraction technique. Many techniques for extracting DNA are well known in the art. See, for example, Gurvitz et al. Australas Biotechnol. 1994 March-April; 4(2):88-91; Ma et al. J Forensic Sci Soc. 1994 October-December; 34(4):231-5; Laber et al. J Forensic Sci. 1992 March; 37(2):404-24. Methods of obtaining DNA profiles from the extracted DNA are well known in the art and include, but are not limited to, DNA sequencing, restriction digestion, polymerase chain reaction followed by electrophoresis, and microarray analysis.
Referring to FIG. 5, in one embodiment, an unknown biological specimen is identified as likely related to a family of at least two family members by its position in a rank-ordered “combined list.” The combined list 140 comprises two or more family member lists, for example, obtained from relatives of missing persons. A “family member list” 130 comprises a set of unknown biological specimens rank-ordered according to a calculated likelihood function value that a genetic relationship of a specific type exists between test DNA profiles obtained from the set of unknown biological specimens and a target DNA profile obtained from a biological specimen of a family member. The combined list 140 comprises the set of unknown biological specimens rank-ordered according to a calculated probability/likelihood that a genetic relationship exists between test DNA profiles and target DNA profiles from at least two family members. A highly ranked unknown biological specimen is identified as likely related to the family at 150.
One embodiment of the present invention involves searching at least one database containing test DNA profiles derived from unknown biological specimens (for example, comprising disaster victims) and comparing these profiles to target DNA profiles from known family members (surviving family members of disaster victims). Test DNA profiles may all be of one type (e.g., a set of short tandem repeats; STR DNA) or there may be multiple types of test DNA profiles for one or more samples. Test DNA profiles can be stored in a database according to profile type (DNA typing technology), or a single database may store DNA profiles having information about multiple DNA profile types. For example, STR test DNA profiles may be stored in a STR test DNA profile database; mitochondrial test DNA (mtDNA) profiles may be stored in a corresponding mitochondrial test DNA profile database and Y-STR results stored in a corresponding database and so on.
Referring again to FIG. 4, in block 110, the DNA profile comparison can be made manually but implementation on a computer, such as a personal computer, is preferred. Thus, as will be further described with reference to FIG. 7, a database preferably contains a search engine, for example, processor 504, which implements search, match, and retrieval methods for a particular test DNA profile type. Search engines preferably are both efficient and scale sub-linearly (e.g., logarithmically) with database size.
In block 120, a probability or likelihood that a genetic relationship exists can be determined by a variety of methods. The likelihood that a relationship exists between a test DNA profile from an unknown biological specimen and a target DNA profile from a family member can be evaluated mathematically. For example, an unknown profile cannot be excluded as a child of a DNA typed parent if the unknown biological sample has at least one STR allele in common with the parent at each locus. For other family members, one uses likelihood ratios (LR's discussed above) or a likelihood function for this determination, but other evaluation methods can be used, such as expert systems or heuristic rules. An exemplary expert system is described by U.S. Pat. Nos. 7,624,087; 7,640,223; 7,664,087; 7,840,519 and 7,945,526. A probability, in this sense, is defined as the mathematical chance that a given genetic profile is observed for the test profile conditioned on the genetic profiles observed for the family relatives under the hypothetical biological relationship between missing person and family. Note that the hypothetical relationship may be in the usual sense of a relationship between two members of a family, or it may be the hypothesis that a relationship does not exist or cannot be determined (hypothesis 2). Probability is expressed as a number between zero and one. Zero means that it is impossible to observe the test profile under the hypothetical relationship, and one implies that the test profile would be observed with certainty under the hypothetical relationship and the observed profile of the family member(s). A likelihood is defined as the reciprocal of a probability. A “likelihood ratio” or LR is defined as a ratio of two probabilities. The probability of the observation of a set of profiles under the assumption that the specimens are related divided by the probability of the same set of profiles under the assumption that the specimens are not related would be one example of a likelihood ratio. The relationship may be between two specimens, or between a specimen and specimens associated by members of a family pedigree, possibly including missing persons. In any of these embodiments, a probability may be replaced by another measure of the possibility or likelihood, or a likelihood function of a relationship. A “likelihood function” may be defined herein to comprise any probability, likelihood, or likelihood ratio calculation that evaluates the degree of manifestation of a genetic relationship between at least two DNA profiles, biological specimens, graphs, or possibly anything which represents an individual or a family.
“Filtering”, as defined in these embodiments, is the exclusion of specimens, or DNA profiles that are not compatible with a family pedigree, or possibly additional information, based upon either DNA or non-DNA information. A “family pedigree” comprises information about the members of a family, where the family may include biological specimens, individuals, or missing persons and the information can be DNA information or non-DNA information, and information, either explicit or implicit, that describes the relationships among the members of a family. Filtering or excluding also refers to narrowing or limiting potential search results. It is desirable to perform such filtering using efficient database search methods. The methods of organizing data into clusters, indexed storage and retrieval of multidimensional information and partitioning data records described in U.S. Pat. Nos. 6,741,983; 7,272,612; 7,454,411; 7,769,803; 7,882,106; 8,060,522 and 8,099,733 are suitable for DNA information clustering. Efficient database search methods are well-known for non-DNA information. A database may be parallelized in order to enhance performance and/or scalability, as described, for example, in U.S. Pat. No. 7,454,411.
For example, if in a family pedigree comprising a father, mother, and two children, typed STR DNA is available for the father, and a child is missing, a search of a database of unidentified human remains can be performed to return all specimens having DNA profiles that are compatible with the father's STR DNA profile. In this case, STR DNA profiles are classified as compatible when at least one allele is shared by the father and child at each STR DNA locus, except for loci of the Y chromosome, in which case all Y-STR DNA alleles are shared if the child is male. If, for example, within a family pedigree, mtDNA is available for a mother or a relative sharing the same maternal lineage, then a search could be performed to return all specimens having mtDNA profiles that match the mtDNA of the typed mother or the maternally linked relative. Additionally, non-DNA forensic information can be used, for example, to exclude samples from remains found before the last date and time the child was seen. The benefit of using a search to filter or exclude specimens is that one has the ability to use an efficient database search method to rapidly prune the size of the population of specimens that must be further analyzed and ranked using, for example, a likelihood function, which may be less efficient than a database search.
Once the probabilities or likelihoods of genetic relationships between unknown biological specimens and family members are evaluated, the unknown biological specimens are rank-ordered according to the likelihood value of being related to a family member. Referring again to FIG. 4, a family member list 130 can be created for at least one family member. The family member list created in block 130 may contain unknown biological specimens listed in order of their likely relationship to the family member. At least two single lists from different family members are then combined into one combined list 140. For example, family member lists from a mother, father, or grandfather can be combined into one combined list 140 which represents the list for the family. Lists can be combined by summing, averaging, or using other mathematical operations to merge data from at least two lists into one combined list. Such combining can be accomplished by performing mathematical operations on the likelihood function values used to order the biological specimens of each list, by performing mathematical operations on the positions or relative positions of the biological specimens of each list, or by a combination of these methods. It can be desirable to exclude an unknown biological specimen from a combined list if the specimen can be excluded from one of the family member lists. It can also be desirable to exclude unknown biological specimen from a list if the specimen has been identified. In block 150, the top few highest-ranked unknown biological specimens in the combined list contain the unknown biological specimen that is most likely the correct missing person for this family. More generally, the unknown biological specimen that corresponds to a specified missing person is likely to be highly ranked.
In another embodiment, referring to FIG. 5, when multiple genetically related family members are available to donate their DNA, the DNA should be used jointly and simultaneously, for example, according to automated decision support described above, to derive a joint likelihood ratio (JLR), or pedigree likelihood ratio (PLR). This PLR between the missing person/victim and the family pedigree is derived under the hypothesis 1 that the missing person is related in the stipulated manner to the family pedigree, against the hypothesis that he/she is unrelated to the family pedigree (hypothesis 2). For a family pedigree under consideration, the PLR for each potential missing person candidate is to be calculated, and the missing person candidates are then ranked according to the PLR. The top ranked candidates should contain the correct missing person sought after by the family corresponding to this family pedigree. Mathematically, the more members of a family pedigree that are available and typed according to automated decision support discussed above or the closer they are related to the missing person, the more precise and effective the ranking will be. PLR is based on the set of available genotype information of the multiple family members, explicitly incorporating into its calculation the known relationships among the family members and the putative relationship to the missing person candidate.
In this embodiment, referring to FIG. 5, relationships between multiple known target DNA profiles from a family which are described by a family pedigree are used to evaluate the likelihood that the unknown specimen is from the missing person related to the family (hypothesis 1). In this embodiment, in block 210, a test DNA profile from at least one unknown biological specimen is compared to the genetic information of the family pedigree which comprises target DNA profiles obtained from biological specimens of at least two family members. As explained above, an ideal family pedigree exists when DNA typing of both parents of a missing person is available. The comparison can be made manually or, preferably, the comparison is implemented in the form of an algorithm installed on a computer. In block 220, a modified Elston Stewart algorithm can be retrieved from computer memory 506, 508, 510 as shown in FIG. 7 for calculating a pedigree likelihood ratio that a genetic relationship exists between at least one unknown biological specimen and the members of the family pedigree. In block 230, a modified Elston Stewart algorithm accomplishes this mathematical evaluation by providing a pedigree likelihood ratio that a test DNA profile from an unknown biological specimen has a genetic relationship of a specific type to a family pedigree.
The original Elston Stewart algorithm (R. Elston, J. Stewart, “A general model for the genetic analysis of pedigree data.” Hum Hered. 21 (1971) 523-542), incorporated by reference herein as to any information deemed essential to an understanding of the embodiments includes a penetrance term expressing the probability of phenotypic expression of a disease if the individual bears the genotypic disease marker. It also includes a probability measure for the degree of linkage of markers residing at different loci. The Elston and Stewart algorithm was adapted to use the penetrance term instead to represent the probability of mutation at an allele that occurs during meiosis from a parent to child lineage inheritance event, if occurrence of mutation is to be considered. If mutation occurrence is not to be considered, then the penetrance term is set to 1. In addition, the term that expresses a linked relationship between markers residing at different genetic loci is not used. The original algorithm is adapted to accommodate the assumption of independence of marker alleles residing at different loci, although should an assumption of independence of marker alleles not be valid for some future DNA type adopted by the forensics community, the term that expresses a linked relationship can be re-inserted in the algorithm. An example mathematical equation for this pedigree likelihood ratio (PLR), without considering mutational occurrence, for a family pedigree comprised of a father, mother, and sibling(s), where either the father, mother, or both can be genetically un-typed is




  
    P
    ⁢
    
      
    
    ⁢
    L
    ⁢
    
      
    
    ⁢
    R
  
  =
  
    
      
        ∑
        
          
            F
            i
          
          ∈
          F
        
      
      ⁢
      
        
          ∑
          
            
              M
              j
            
            ∈
            M
          
        
        ⁢
        
          
            p
            ⁡
            
              (
              
                F
                i
              
              )
            
          
          ⁢
          
            p
            ⁡
            
              (
              
                M
                j
              
              )
            
          
          ⁢
          
            
              ∏
              
                
                  S
                  k
                
                ∈
                S
              
            
            ⁢
            
              
            
            ⁢
            
              
                p
                ⁡
                
                  (
                  
                    
                      
                        S
                        k
                      
                      ❘
                      
                        F
                        i
                      
                    
                    ,
                    
                      M
                      j
                    
                    ,
                    
                      H
                      0
                    
                  
                  )
                
              
              ⁢
              
                p
                ⁡
                
                  (
                  
                    
                      U
                      ❘
                      
                        F
                        i
                      
                    
                    ,
                    
                      M
                      j
                    
                    ,
                    
                      H
                      0
                    
                  
                  )
                
              
            
          
        
      
    
    
      
        ∑
        
          
            F
            i
          
          ∈
          F
        
      
      ⁢
      
        
          ∑
          
            
              M
              j
            
            ∈
            M
          
        
        ⁢
        
          
            p
            ⁡
            
              (
              
                F
                i
              
              )
            
          
          ⁢
          
            p
            ⁡
            
              (
              
                M
                j
              
              )
            
          
          ⁢
          
            
              ∏
              
                
                  S
                  k
                
                ∈
                S
              
            
            ⁢
            
              
            
            ⁢
            
              
                p
                ⁡
                
                  (
                  
                    
                      
                        S
                        k
                      
                      ❘
                      
                        F
                        i
                      
                    
                    ,
                    
                      M
                      j
                    
                    ,
                    
                      H
                      1
                    
                  
                  )
                
              
              ⁢
              
                p
                ⁡
                
                  (
                  
                    U
                    ❘
                    
                      H
                      1
                    
                  
                  )
                
              
            
          
        
      
    
  




where F and M are the sets of possible DNA profiles for the father and mother, respectively, S is the set of DNA profiles for typed siblings, and U is the DNA profile for the unidentified human remain. H0 is the hypothesis that the unidentified human remain is from a missing child in this family pedigree, and H1 is the hypothesis that there is no known relationship between the unidentified human remain and the family pedigree. This equation is one example PLR for this family pedigree. Equations for family pedigrees comprised of different combinations of family members can be written using the same adaptation to the Elston Stewart algorithm or similar published works, such as those by E. Lander and P. Green (Lander and Green 1987) or by C. Cannings, E. Thompson, and E. Skolnick (Cannings et al. 1976) and others (Dudbridge F. 2003; Fernandez et al. 2001).

A PLR can be evaluated in different ways. For example, a straightforward way to evaluate the above equation is to enumerate all possible unknown profiles and multiply and sum over these enumerated unknown profiles the probabilities and conditional probabilities that depend upon them. Simplifications can be used, as indicated above, when the profile of either the father or the mother, or both, is known. Many loci are independent from each other, and in this case a term can be calculated for each independent locus, and the product of these terms can then be incorporated in the PLR. Another way is to consider the number of alleles or genotypes that are inherited by a child from one or both parents and use a probability or conditional probability that this number has occurred. This is known as “identical by descent” (IBD), and in this case the terms of the PLR are grouped according to the number of alleles or genotypes IBD. A calculation using an IBD approach can consider IBD alleles between any two genetically related biological specimens or individuals; it is not restricted to parent/child relationships. Alternatively, portions of the equation for the PLR may be pre-computed. This can be especially beneficial when the PLR is to be calculated for a number of unknown human remains; for example, the product terms over all known siblings may be pre-computed. Another way is to predetermine a sequence of machine instructions that can be executed to calculate the PLR for a given family pedigree in a highly efficient manner, and then execute these instructions for each of a set of unknown human remains. Finally, if the family pedigree incorporates multiple generations, the PLR may be computed in a recursive manner, in either a bottom-up or top-down fashion. The family pedigree may contain loops, in which case the approach of C. Cannings, E. Thompson, and E. Skolnick, or a similar approach, may be employed.
Pedigree likelihood ratios, as described above, can be used to evaluate whether there is a reasonable degree of certainty that an identification determined in this manner is correct. Optionally, as shown in block 240, the unknown biological specimens can be rank-ordered according to a pedigree likelihood ratio that at least one unknown biological specimen is related to the family or family pedigree. Higher-ranked unknown biological specimens are more likely to be candidates for the correct missing family member sought by the family represented by the family pedigree.
Another embodiment of the invention, referring to FIG. 6, takes advantage of known or hypothetical relationships between biological specimens. This embodiment is particularly suitable to large-scale mass disaster or terrorist incidents where several related members of families are missing or among the victims, such as occur in airplane crash disasters and natural disasters such as volcano eruptions or tsunamis. A database of DNA profile data from both known and unknown biological specimens can have hypothetical relationships between specimens, some of which may be known to be true, for example, samples taken from various portions of an unidentified remain, and samples from known or related sources. Hypothetical relationships may be based upon, for example, additional forensic evidence linking unidentified remains. Hypothetical relationships can also be tested using standard methods such as likelihood functions, maximum likelihood identification, or maximum likelihood estimation that provide a quantitative measure of the potential validity of each hypothesis. Discovered relationships are hypothetical relationships having a high degree of correctness, and may be based upon, for example, non-DNA forensic information (see below) or mathematical analysis such as the use of likelihood functions or maximum likelihood identification and a decision process whereby a hypothesis of relatedness is accepted and recorded or discarded. One possible decision process is comparison of the value of a likelihood function to a numeric threshold and acceptance if the value is greater than the threshold.
Referring again to FIG. 6, to determine whether a relationship may exist between two biological specimens (known or unknown), block 310, a DNA profile from one biological specimen is compared to a DNA profile from another biological specimen to determine the probability or likelihood that a genetic relationship exists. Such a relationship can be called a hypothetical relationship. The comparison can be carried out by calculation of the value of a likelihood function and comparison of the value to a threshold. A second approach is to use maximum likelihood estimation or to compute the joint conditional probability of occurrence of the two DNA profiles for a set of hypothetical relationships, and accepting and recording one or more of the relationships having greatest likelihood or probability as having support from the available information. Once two biological specimens are compared, a likelihood of a hypothetical relationship between the two specimens is calculated as described above and shown in block 320. It is important to recognize that multiple hypothetical relationships between two specimens may be accepted for recording, and that each specimen may have accepted relationships with zero, one, or more other specimens. In this context, “accepted for recording” means that the evidence for the relationship, for example, based upon forensic evidence or calculated likelihood, is sufficiently strong to indicate that the possibility of this relationship should be recorded, or maintained, for example, in a database for further evaluation or analysis. Each hypothetical or known relationship between the biological specimens can be represented, or recorded, as an edge of a directed graph 330. The term “recorded” is used in these embodiments to show that a relationship has been represented by an edge of a directed graph. A directed graph, as defined in the field of computer science, contains nodes which correspond to specific data and edges which correspond to relationships among the data. An edge has a direction from one node to another, and this edge can represent a dependence of one node upon another, such as “node A is a child of node B”. The term “directed graph” is used here to indicate that such dependencies can be represented; however, other methods of representing such dependencies are known in the art, such as entity-relationship diagrams or relations between records in tables of a relationship database. Therefore, the term “directed graph” is defined to mean any representation of entities or objects, such as are needed to represent biological specimens, DNA profiles, or individuals, and relationships between these entities or objects that are not necessarily reflexive, meaning, for example, that “node A is a child of node B” does not imply that “node B is a child of node A”. Information may be associated with any node or edge to more completely describe the data or relationship. In these methods, nodes of the directed graph correspond to individual biological specimens and/or their DNA profiles, and edges of the directed graph correspond to known or hypothetical and accepted genetic relationships between the biological specimens. In block 340, a family pedigree comprising relationships among family members can also be represented as a directed graph.
In this embodiment of FIG. 6, a directed graph that represents a hypothetical or known relationship between biological specimens or individuals and/or the corresponding test DNA profiles is compared to a directed graph that represents a family pedigree defining relationships between individuals (known or missing) at block 340. The objective of the comparison is to identify portions of the directed graph that represent hypothetical or known relationships between biological specimens or individuals and/or the corresponding test DNA profiles that at least partially correspond to the structure and information specified by the directed graph that represents a family pedigree. The comparison, as shown in block 340, can involve the structure of the directed graph in addition to the DNA profile information contained in the directed graph. For example, the directed graph that represents a family pedigree may describe parent and child relationships of a family, where many families are missing due to a mass disaster, such as occurred subsequent to the Southeast Asian/Pacific tsunami of 2005. Such a comparison may yield many possible matches, corresponding to multiple missing families. If additional known information is included in the directed graph of the family pedigree and its associated data (such as DNA profiles of known individuals or non-DNA forensic evidence), the search would be more focused and may yield one to only a few possible matches. Comparison may be done manually or, according to a computer-implemented algorithm on a computer system according to FIG. 7. One approach to implementing this method is to build upon existing directed sub-graph matching or (partial) isomorphism search and retrieval methods. A search can identify isomorphic structures in the database corresponding (partially) to a target directed graph, and compare test DNA profile information to determine any candidate matches to target DNA profiles which lie within the defined neighborhoods, as shown in block 350. Multiple test DNA profiles in the database may match simultaneously within a stored directed graph structure, corresponding to multiple test DNA profile types corresponding to an unknown biological specimen and/or the test DNA profiles corresponding to multiple biological specimens associated by a sub-graph that is at least partially isomorphic to the target directed graph, as shown in block 350. A second approach to implementing this method utilizes link discovery methods to discover clusters within the directed graph that is constructed as described herein. These clusters correspond to groups of related specimens and individuals and include the family pedigree information represented within the database. Link discovery methods can also be used to identify and visually inspect clusters of related nodes corresponding to individuals and/or remains.
One way in which this embodiment can be implemented is to construct a database that contains information about specimens and individuals. Note that a “database” may be implemented as a collection of databases or database tables, or it may be any other method of organized information storage. Each specimen or individual can be represented in the database as a node or database object, as this term is known in the literature, corresponding to a node of a directed graph. DNA and non-DNA information about each specimen or individual can be associated with this node using methods that are known in the field of computer science. Known relationships between specimens or individuals can also be stored in the database and correspond to edges of a directed graph containing the nodes. A list of possible relationships between individuals or specimen can be maintained, and the list defines the set of possible hypothetical relationships between specimen or individuals. For any pair of specimens or individuals represented in the database, a likelihood function value of each possible hypothetical relationship can be calculated based upon the information stored in the database, and this likelihood function value can be tested to determine if information about a hypothetical relationship between this pair of specimens or individuals should be stored in the database. If the result of this test determines that the hypothetical relationship should be stored or recorded, the hypothetical relationship can be recorded as an edge in the database with associated information including, for example, the hypothetical relationship and the likelihood that the hypothetical relationship is true. More than one hypothetical relationship may be recorded in this manner for any pair, although in many cases no hypothetical relationship will be recorded because the available information does not support the hypothesis that a relationship exists. In this manner, a database can be constructed that contains DNA and non-DNA information about specimens and individuals and their relationships and includes hypothetical relationships that are evaluated to be sufficiently likely to be true. The terms “likely” and “likelihood” may correspond to the mathematical definitions of likelihood functions or likelihood ratios from the fields of population genetics, statistics, or probability, but this is not necessary. Other functions that quantify the strength of a hypothetical relationship between two individuals or specimens relative to other hypothetical relationships may be used, such as probabilities, log likelihood functions, or measures of genetic similarity such as the number of genetic markers in common, the relative lengths of two genetic sequences, or the edit distance between two genetic sequences.
An advantage of this method over existing methods is its ability to simultaneously hypothesize and record multiple possible relationships. A second advantage is that a specimen may be related to multiple family pedigrees through edges of the directed graph that correspond to accepted and recorded hypotheses. Thus, unlike previous methods, the specimens are not partitioned into disjoint subsets that are hypothetically related to a single family pedigree. Such a partitioning can result in classification errors, since at the partitioning step insufficient information is known to accurately perform the partitioning. An incorrect assignment of a specimen to a partition in this case would eliminate the possibility of correctly identifying the specimen's family pedigree at a later time. In contrast, a method of one embodiment allows a specimen to hypothetically belong to more than one family pedigree, and subsequent mathematical analysis, such as by using a pedigree likelihood ratio, can be used to determine the correct family pedigree.
A family pedigree containing information about individuals and relationships can be used as a target for a query or search of this database to locate portions of the database that are likely to correspond to the structure and information provided by the family pedigree. Such correspondence can be further evaluated using a pedigree likelihood ratio such as a modified Elston Stewart algorithm or other quantitative evaluation of the degree of correspondence of that portion to the family pedigree, relative to other possible portions of the database. In this manner, multiple portions of the database, which can be described as sub-graphs, can be identified as possibly corresponding to the family pedigree, and such portions can be rank ordered using, for example, a pedigree likelihood ratio. Since there may be ambiguity within each portion as to the relationships among individuals and samples, evaluations of multiple possible family pedigrees within this portion can be performed to determine the more likely family pedigrees, providing a method of identification of possible relationships among, for example, unidentified human remains, using methods that are well-known in the field of DNA forensics, similar to those used, for example, to analyze the mass grave containing the probable remains of the Romanov family.
An alternative method is to search first for DNA profiles with likely or probable genetic relationships, followed by a search for matching directed sub-graphs to the target directed graph. Preferably, one would perform these searches simultaneously. One method for implementing this is to alternately resolve searches for potentially related DNA profiles and sub-graph structures at each level of a database index structure, where each level successively reduces the set of possible matches to the target directed graph and DNA profiles within the database.
The process of discovering relationships among biological specimens can lead to the identification of at least two biological specimens with similar DNA profiles, suggesting that these specimens are derived from the same individual. Biological specimens with similar DNA profiles can be collapsed into one specimen to enhance the speed and efficiency of the search.
An advantage of the rapid advances over the past decades in available computational capabilities is that idle computational capability is either available during periods of reduced utilization, or can be made available through the addition of computational elements or partitioning of available computational resources among several tasks. A computer-implemented program running on one or shared among many processors can allocate these unused computational capabilities to background processes that do not have to perform in a time-critical environment. It is advantageous to utilize background processes which may be likewise computer-implemented, also referred to as “spiders,” to examine stored data during idle periods, determine hypothetical relationships between samples and their corresponding DNA profiles stored in a database, perform analyses to assess the potential truth of these hypothetical relationships, and record edges in the database indicating the presence of these hypothetical relationships. Such relationships can subsequently form the basis for information that can be compared to target family pedigrees, or directed graphs representing those family pedigrees, and associated DNA profiles in order to achieve identification of unknown remains, associate such remains with missing individuals, and link these items to relatives. Such spiders, or background processes, can thus provide a substantial benefit by improving the accuracy, efficiency, and speed of identification of unknown biological specimens.
Information about unidentified human remains, known individuals, missing persons, and/or family pedigrees can be collected over a substantial period of time, such as several months to a few years. In this case, the method is applied in a continuous or intermittent manner over at least a portion of this period of time. As new information is added, spiders, or background processes, can evaluate hypothetical relationships between the new information and pre-existing information and record edges in the database indicating the presence of these hypothetical relationships. Such newly recorded relationships can subsequently form the basis for information that can be compared to target family pedigrees, or directed graphs, and associated DNA profiles in order to achieve new identification of unknown remains, associate these remains with missing individuals, and link these items to relatives. Both new information and the results of identifications can enable further identification. For example, identification of a second parent or sibling in a family can be used to assist in identification of additional siblings. New information can also lead to the exclusion of some possible identifications. For example, the identification of an unidentified human remain excludes the identification of this remain as originating from a different person.
The methods described by these embodiments can be either static or dynamic. The method is static if the method is applied at one time to a collection of information that has been gathered about two or more biological specimens from unknown remains, known individuals, and/or missing persons. It is more common for the method to be applied in a dynamic manner. In this case, information is gathered over a period of time, and the method is applied in a continuous or intermittent manner over a least a portion of this period of time. New information is added over a time period that can be fairly long, such as several months to a few years, and identifications occur as enough data accumulate to support them. Both newly arrived data and the results of identifications can enable further identification. For example, identification of a second parent or sibling in a family can be used to assist in identification of additional siblings. Newly arrived data can also lead to the exclusion of some possible identifications. For example, the identification of an unidentified human remain excludes the identification of this remain as originating from a different person.
In any of the methods described herein, in addition to DNA profile information, non-DNA forensic information can be used to filter or exclude a set of unknown biological specimens to include those which fit one or more requirements specific to the missing individual or to the family searching to find a missing individual. “Non-DNA forensic information” as defined herein includes, but is not limited to, dental records, skin markings (e.g., tattoos, scars), X-ray images, time of specimen collection, place of specimen collection, fingerprints, and gender as well as personal effects and clothing or fragments thereof. As alluded to above, following the principles of U.S. Published Patent Application 2012/0148115 (biometrics), biometric data may be collected at a disaster or crime scene utilizing an iris scanner, fingerprint scanner, a camera and related input devices 514 of a computer system 500 (FIG. 5) for assisting (in addition to DNA analysis) in the identification of human remains representing an individual. Following the principles of U.S. Published Patent Application 2010/0332474 (predicting object properties), the clothing on a victim or of a crime perpetrator or other personal effect may comprise an object on which may be performed spectral data analysis, microbody assemblage data analysis and content-based image retrieval analysis, for example, to utilize, for example, a camera to collect an image of a tattoo or to collect data of from an associated personal effect to predict an object property. For example, in a group of unknown biological specimens, there may be body parts that contain tattoos. The tattoo may be segregated from a camera image of the biological specimen using a tree of shapes method, multivariate statistical analysis or another method of content-based image retrieval. On the other hand, if a family is searching for a missing individual and knows that the missing individual does not have a tattoo, for example, this piece of information could exclude any unknown biological specimens which have a tattoo regardless of how welt the DNA profiles seem to be related. As a second example, skeletal or other anatomical features may provide gender information. On the other hand, personal effects and the like should be relied on with some degree of risk that coincidence (or intention) placed an item proximate to certain remains.
Non-DNA forensic information also can be used to confirm the identification of an unknown biological specimen. Non-DNA forensic information can either be stored in a database which also contains DNA profile information, or in a separate database. For example, a cluster of unknown biological specimens identified as likely belonging to a missing individual may contain a jawbone which includes teeth. After carrying out a search method, dental records of a database such as dental x-rays taken from the missing individual can be used to verify that the identification is correct. It is preferred to exclude an unknown biological sample from the set of unknown biological samples after identification.
Implementation Mechanisms-Hardware Overview
Methods of the first embodiment and subsequent embodiments may be utilized in connection with computer readable media which may be provided for temporary or permanent storage in a computer system 500 (FIG. 7) such as a personal computer, intelligent telecommunications device or other computer or computer system comprising one or more parallel processors known in the art. FIG. 7 is a block diagram that illustrates a computer system 500 upon which at least one embodiment of the invention may be implemented. Computer system 500 includes a bus 502 or other communication mechanism for communicating information, and at least one processor 504 coupled with bus 502 for processing information. Computer system 500 also includes a main memory 506, such as a random access memory (“RAM”) or other dynamic storage device, coupled to bus 502 for storing information and instructions to be executed by processor 504. Main memory 506 also may be used for storing temporary variables or other intermediate information during execution of instructions to be executed by processor 504. Computer system 500 may further include a read only memory (“ROM”) 508 or other static storage device coupled to bus 502 for storing static information and instructions for processor 504. A storage device 510, such as a magnetic disk, optical disk, solid-state memory, or the like, may be provided and coupled to bus 502 for storing information and instructions. A storage device may be removable using a coupling mechanism such as a universal serial bus (USB) or other hardware specific to the type of storage hardware, such as a CompactFlash, SD, or microSD card reader or port. The removable storage device (not shown) may be utilized to transfer information to or from computer system 500. Any of memories 506, 508, 510 may retain program instructions according to any embodiment of ADS, for collecting missing person profile data or profile data for a family pedigree or for associating unknown specimens with family pedigrees.
Computer system 500 may optionally be coupled via bus 502 to a display 512, such as a cathode ray tube (“CRT”), liquid crystal display (“LCD”), plasma display, television, small mobile telephone display or the like, for displaying information to a computer user. Display 512 may provide a virtual keyboard for data input or a one or two dimensional bar code reader. Alternatively, information may be delivered to or collected from a computer user or another computer system or computer program using a communication interface 518 or removable storage device. An input device 514 may comprise a barcode reader, a fingerprint scanner, a camera, an iris scanner, a physical or virtual keyboard including alphanumeric and other keys or other input device such as a spectral data or microbody assemblage data collector, may be coupled to bus 502 for communicating information and command selections to processor 504. An optional type of user input device is cursor control 516, such as a mouse, trackball, stylus, or cursor direction keys, for example, as may be found on some personal data assistants (PDA's) for communicating direction information and command selections to processor 504 and for controlling cursor movement on display 512. This input device 516 typically has two degrees of freedom in two axes, a first axis (e.g., x) and a second axis (e.g., y), that allows the device to specify positions in a plane. Such an input device may be combined with a display device such as a LCD with a touch screen, commonly found on mobile telephones or other telecommunications or presentation devices such as the Apple iPad or a computer tablet using the Android operating system or other operating system having similar features. Alternatively, information and command selections may be communicated to processor 504 using a communication interface 518. Optionally, separate communication interfaces may be used to deliver information to a computer user or another computer system or computer program, and to communicate information and command selections to processor 504.
The invention is related to the use of computer system 500 for automated decision support, for collection of specimen data at a disaster site or site of a selected family member or a laboratory and for identifying an unknown biological specimen as likely related to a family comprising at least a first and second family member. According to one embodiment of the invention, identifying an unknown biological specimen as likely related to a family comprising at least a first and second family member is provided by computer system 500 in response to processor 504 executing one or more sequences of one or more instructions contained in main memory 506, for example, in accordance with ADS and the flowcharts of FIGS. 4 through 6. Such instructions may be read into main memory 506 from another computer-readable medium, such as storage device 510. Execution of the sequences of instructions contained in main memory 506 causes processor 504 to perform the process steps described herein. In alternative embodiments, hard-wired circuitry may be used in place of or in combination with software instructions to implement embodiments of the invention. For example, a field-programmable gate array (FPGA) or application-specific integrated circuit (ASIC) may be used. Such a device can, for example, implement associative memory to aid in indexing, search, and retrieval of information stored in a database. A second example is use of a FPGA or ASIC to speed up calculation of a likelihood function used to rank specimens (block 140 of FIG. 4; block 240 of FIG. 5). Thus, embodiments of the invention are not limited to any specific combination of hardware circuitry and software.
The term “computer-readable medium” as used herein refers to any medium that participates in providing instructions to processor 504 for execution or storing information in a form that can be accessed by processor 504. Such a medium may take many forms, including but not limited to, non-volatile media, volatile media, and transmission media. Non-volatile media includes, for example, optical or magnetic disks, solid state memories, and the like, such as storage device 510. Volatile media includes dynamic memory, such as main memory 506. Transmission media includes coaxial cables, copper wire and fiber optics, including the wires that comprise bus 502. Transmission media can also take the form of acoustic or light waves, such as those generated during radio wave and infrared data communications.
Common forms of computer-readable media include, for example, a floppy disk, a flexible disk, hard disk, magnetic tape, or any other magnetic medium, a CD-ROM, any other optical medium, solid-state memory, punch cards, paper tape, any other physical medium with patterns of holes, a RAM, a PROM, and EPROM, a FLASH-EPROM, any other memory chip or cartridge, a carrier wave as described hereinafter, or any other medium from which a computer can read. Various forms of computer readable media may be involved in carrying one or more sequences of one or more instructions to processor 504 for execution.
Computer system 500 may also include a communication interface 518 coupled to bus 502. Communication interface 518 provides a two-way data communication coupling to a network link 520 that is connected to a local network 522. For example, communication interface 518 may be an integrated services digital network (“ISDN”) or digital subscriber line (“DSL”) card or a modem to provide a data communication connection to a corresponding type of telephone line, fiber optic or wireless link. Preferably, communications transmitted over such a link are encrypted or otherwise protected according to known encryption schemes and/or watermarking algorithms to uniquely identify a source. As another example, communication interface 518 may be a network card (e.g., an Ethernet card) to provide a data communication connection to a compatible local area network (“LAN”) or wide area network (“WAN”), such as the Internet or a private network. Wireless links are implemented in the about-to-be-described example of running a mobile disaster victim identification algorithm on an intelligent telecommunications device using, for example, WiFi, Bluetooth, or third generation (“3G”) or fourth generation (“4G”) wireless technologies such as WiMax or LTE. In any such implementation, communication interface 518 sends and receives electrical, electromagnetic or optical signals that carry digital data streams representing various types of information. For example, a forensic investigation may require a data communication connection to a database comprising at least DNA profile data or other forensic information. A second example is use of one or more data communication connection(s) to access at least one database used to store DNA and/or non-DNA information. Portions of the computations associated with the ADS tool and methods of FIGS. 4-6 described herein may be distributed across multiple computer systems 500 which may communicate using one or more communication interfaces 518.
Network link 520 typically provides data communication through one or more networks to other data devices. For example, network link 520 may provide a connection through local network 522 to a host computer 524 or to data equipment operated by an Internet Service Provider or private network service provider (“ISP”). It is intended that such a service provider may operate in a “cloud” computing environment such that it is a web accessible service for, for example, identifying a perpetrator of a crime, a victim or missing person. An ISP in turn provides data communication services through a packet data communication network such as the worldwide network commonly referred to as the “Internet” 528 or a private network. An example of a private network is a secure data network linking law enforcement agencies and used for transmission of DNA and/or non-DNA information. Local network 522 and Internet 528 both use electrical, electromagnetic or optical signals that carry digital data streams. The signals through the various networks and the signals on network link 520 and through communication interface 518, which carry the digital data to and from computer system 500, are exemplary forms of carrier waves transporting the information.
Computer system 500 can send messages, commands and data and receive messages, commands and data, including program code, through the network(s), network link 520 and communication interface 518. In the Internet example, a server 530 might transmit a requested code for an application program (ADS or the algorithms of FIGS. 4-6) through Internet 528, host computer 524, local network 522 and communication interface 518. In accordance with the invention, one such downloaded application provides a method of identifying an unknown biological specimen as likely related to a family comprising at least a first and second family member or a mobile disaster victim identification algorithm or a least squares deconvolution mixture algorithm for resolving DNA mixtures or a peak-fitting algorithm or an automated expert system.
The received code may be executed by processor 504 as it is received, and/or stored in storage device 510, or other tangible computer-readable medium (e.g., non-volatile storage) for later execution. In this manner, computer system 500 may obtain application code and/or data in the form of an intangible computer-readable medium such as a carrier wave, modulated data signal, or other propagated signal. An embodiment further contemplates download and execution of programs or portions of programs (software) from a digital software marketplace or a secure web site maintained by either a vendor or, for example, a government agency. The depicted network could also be a secure network used by the military during combat operations, surveillance or reconnaissance, engineering activities, or training.
Computer system 500 can be configured using the methods of this invention to provide services across a network to forensic personnel having client computers or intelligent telecommunications devices capable of connection to a network or other communication interface. Such services can include assistance in identification of unidentified remains and storage of DNA and non-DNA information about known individuals and missing persons and their relationships, or a family pedigree. These services can also be provided to other software, located in either computer system 500 or a separate computer system connected by a network, network link, or communication interface to computer system 500. The services can be protected using methods of authentication and/or encryption that are known in the fields of computer science and computer security in order to ensure data are neither compromised nor disclosed and to trace all accesses to the data. The computer system 500 and other associated information storage and communication components can be protected using devices and methods that are known in the fields of computer science and computer security, such as with firewalls, physical access controls, power conditioning equipment, and backup or redundant power sources. The protection devices and methods, embodied as either hardware, software, or a combination of hardware and software, may be incorporated in computer system 500 or exist as separate components typically connected by a network or other communications hardware. The information stored by computer system 500 and computer-readable media can be further protected using backup or redundant information storage systems, such as those that are well-known in the art. Examples include tape storage systems and RAID storage arrays.
Mobile Disaster Victim Identification
Now a method and apparatus for mobile disaster victim identification will be disclosed and described with reference to FIGS. 8-15 which are display screens for inputting data and outputting data as will be further described herein to identify and to store a family pedigree output of ADS and provide a pedigree likelihood ratio for the input family pedigree. The display screens are exemplary only and may be varied depending on desired content and order in performing a particular forensic identification. In a given forensic identification, a plurality of pedigrees are possible, for example, in a crime situation where more than one victim, perpetrator of a crime, missing person or the like has their own pedigree and are otherwise unknown or may be collected in a mixture of DNA requiring resolution to identify an unknown. The mobile disaster victim identification application may be an application for downloading to and implementation on a wireless telecommunications apparatus, for example, an iPad or iPhone available from Apple Computer, a Motorola or HTC Android telephone, a personal computer, a Blackberry or other intelligent wireless telecommunications device available for use on a public network, for example, Verizon, Sprint or AT&T. A high bandwidth is useful for, for example, fingerprint collection or other high resolution analysis and so a communications interface as describe herein may be wireless 3G or 4G, by FIOS, optical fiber, satellite or other communications interface that may be connectable to intelligent apparatus for gathering data and running algorithms as possible.
FIG. 8A provides an exemplary introductory screen for a mobile disaster victim identification application 801 (shown) and other software such as MP/KN, peak fitting, expert system, least square deconvolution, biometrics support, predicting object properties support and the like as has been introduced above and will be further described below. A mobile terminal may resolve the identity of an individual who may be a victim. The forensic analyst, for example, a member of a police department or fire and rescue team may select among Mobile DVI, LSD or other application to run on an intelligent telecommunications or personal computing device or to be remote remotely from a server or networked cloud. These applications may include automated decision support for selecting a family pedigree (ADS, not shown), matching an unknown biological specimen to a pedigree (MP/KN, not shown), DNA mixture analysis (least square deconvolution or LSD, not shown), peak fitting, not shown, an automated expert system (STRESP, not shown), remote access to a system such as the FBI CODIS system (not shown), mobile biometrics (not shown), predicting object properties (not shown), expert system (not shown) and mobile disaster victim identification (Mobile DVI shown). A first step of a missing person or DNA mixture resolution process may be to set up a hypothetical relationship and identify a pedigree relationship among a missing person or other unknown, a crime perpetrator and a family pedigree for any of these. Automated decision support (ADS) may output, for example, a simple family tree such as a mother, father and child of a missing person, crime perpetrator or contributor to a mixture. It is then necessary to input STR, Y-STR, mtDNA, mass spectrometry or other known laboratory analysis data or otherwise identify available members of each family by identification number, for example, in a CODIS or other database or local entry at a disaster site. For example, portable DNA laboratory kits are anticipated or available for local disaster or accident or fire or crime victim analysis and are anticipated for mass spectrometry. If a woman, any of STR and mitochondrial DNA profiling is available along with mass spectrometry. For a man, YSTR typing among other typing is available. Once a hypothetical pedigree or other relationship is output from ADS, that relationship may be automatically or manually input to a mobile disaster victim identification application running on an intelligent telecommunications device such as a Motorola or other telephone running the Android operating system or an Apple i-phone running the Apple operating system. Via such a device, a local police database, an FBI CODIS database and associated software may be remotely accessed or locally accessed and local software run. One may upload new data to any database or download such data for running software locally. Element 802 may represent the date, Monday, August 9 and a calendar reminder that a status report is due 1:00 PM on August 12.
FIG. 8B provides an introductory menu for a selected mobile disaster victim identification application including specimen capture 810, pedigree entry 820, clear data 830, about 840 and exit 850. Connecting 860 indicates that the device is attempting connection to a server. Referring briefly to FIG. 8C, Connected 880 indicates that the device is connected to the server. If an external server is needed for entering data to a pedigree or for uploading collected data so that MP/KN, LSD, peak fitting, expert analysis or other algorithms or databases may be used, a telecommunications link is employed. The connecting overlay 860 appears in FIG. 8B showing a connection attempt to a server for the exemplary MobileDVI application and connected 880 in FIG. 8C indicates success.
FIG. 9A provides an exemplary screen for selecting among physical profile 910, GPS coordinates 920, DNA Profile 930 and Back 940 (taking a user to an immediately previous screen) adding new specimen information while FIG. 9B relates to editing. In particular, by actuating Physical Profile touch button 910, a user may enter physical profile data. FIG. 9C-D provide exemplary screens for entering physical profile shape data by scrolling through a list similar to one associated with a standard Interpol or other victim identification form silhouette sketch form part. Referring briefly to FIG. 9G, one may see a silhouette sketch form part of an Interpol victim identification form. All the elements of the form may be scrolled to using up, down, left, right commands entered, for example, by touch to a touch screen display of a mobile telecommunications device or a personal computer. Head Shape (Front): 912 of FIG. 9C corresponds to head form, front (Shape of head from side) 912 on FIG. 9E. A Rectangular Head Shape (front) has been selected. Head Shape (Side): 914 of FIG. 9B corresponds to Head form, profile (shape of head from side) 914 on FIG. 9E. A Deep Head Shape (side) has been selected. Further, Forehead Height: 916 of FIG. 9B corresponds to Forehead—Height/Width 916 on FIG. 9G. A Medium Forehead Height has been selected. Further, Depicted partial element 918 of FIG. 9B corresponds to Forehead—Height/Width 918 on FIG. 9G. Further inputs for a Nose Curve 922, Nose Angle 924 and Earlobe 928 are selectable from screen FIG. 9D. Save 928 is used to save physical profile data 910. By scrolling up and down and by touch, a complete silhouette sketch of a standard victim identification form may be completed and stored with other victim data such as location data and hypothetical descriptors such as will be described in connection with personal identification screen, FIG. 14A.
FIG. 9E provides an exemplary screen for displaying GPS data upon actuation of GPS Coordinates touch button 920 of FIG. 9A. Latitude and longitude may be automatically acquired by an intelligent telecommunications device with GPS capability such that latitude and longitude may be automatically obtained, for example, for entry into a Place of Disaster field 950 of a standard Interpol or other victim identification form of FIG. 9F where FIG. 9F shows, for example, a Place of Disaster field 950 of an Interpol victim identification form. A standard Interpol form comprises a plurality of pages for entry of data that may be entered via a mobile telecommunications device having a mobile disaster victim identification application implemented thereon. All fields of the form may be implemented for automatic or manual input including but not limited to biometric data including fingerprint data and dental data and personal effect data. Photographs and other biometric data per U.S. Published Patent Application 2012/0148115 such as, for example, dusted fingerprints, iris or cornea scans and scars, marks, and tattoos (SMT) may be acquired using, for example, a camera or fingerprint scanner input device incorporated in the Motorola telephone and supported by the Android operating system, and form data entry as illustrated in FIGS. 9A-9D, 11A, 12, 13, and 14A. Initially, specimens are collected and bagged to be analyzed by a lab (if they cannot be analyzed on site). The specimen bags may be labeled with single dimensional bar codes identifying family members, as are the samples submitted for laboratory analysis, and thus the reports that come back from analysis have the same identity bar codes so the correct associations are made between the position when evidence is collected and the DNA profiles that result from analysis.
FIG. 10A provides exemplary laboratory or field kit analysis profile data collected for a father family member in the form of a single dimensional bar code uniquely identifying the father and two dimensional bar codes for each of STR, YSTR and mitochondrial DNA profile data, together with human-readable versions of the data represented by the bar codes, processing date and time and analyst and, optionally, the analyst's reviewer's identifiers or names. The single dimensional bar code 1000 in FIG. 10A can be used in physical evidence collection kits to provide a unique identifier for all information related to a collected specimen, and acquired by the Mobile DVI or other application using the mobile device. To input data from a single dimensional bar code or a two dimensional bar code, a bar code of FIG. 10A is manually presented to the bar code reader (an example of input device 514) and input in this case, to identify the father uniquely and profile data for the father. Further data may be input using a display touch screen similar to that of FIG. 14A whereby John Murdoch (father) may be identified by name, sex, birth date and the like and the box missing person not checked.
FIG. 10B provides exemplary laboratory or field kit analysis profile data collected for a mother family member in the form of a single dimensional bar code uniquely identifying the mother and two dimensional bar codes for each of STR and mitochondrial DNA profile data, together with human-readable versions of the data represented by the bar codes, processing date and time and analyst and, optionally, the analyst's reviewer's identifiers or names. As described above for the father, personal data identifying mother, Emma Murdoch, may be entered via the display touch screen of FIG. 14A.
FIG. 10C provides exemplary laboratory or field kit analysis profile data collected for a son family member in the form of a single dimensional bar code uniquely identifying the son and two dimensional bar codes for each of STR, YSTR and mitochondrial DNA profile data, together with human-readable versions of the data represented by the bar codes, processing date and time and analyst and, optionally, the analyst's reviewer's identifiers or names. As described above for the father, personal data identifying son, John Murdoch III, or daughter, Jane Murdoch, may be entered via the display touch screen of FIG. 14A. As taught above, having mother and father of a missing person is sufficient to identify a missing person John Murdoch, Jr. in this example without having to incorporate son or daughter data. According to FIG. 15C, further individuals such as 6—Sarah Jane Smith may be included, for example, as a known individual if a mixture deconvolution analysis is needed in the event of a rape to identify 1—John Murdoch jr. (Missing) as the perpetrator if he is missing and never been typed but is a hypothetical match to a contributor to a DNA mixture extracted from Ms. Smith.
FIG. 11A provides an exemplary screen for entering DNA profile data for a victim or a family member as shown by example in FIGS. 10A-10C and a checklist for data entered for that family member whereby each profile may be individually viewed. STR Profile 1110 remains unchecked until STR profile data is entered. If mass spectrometry data is input, a further box below MtDNA Profile box 1130 may be provided. YSTR Profile 1120 is similarly provided for a male individual and checked if profile data has been input. MtDNA Profile box 1130 is provided and checked if such data has been input. By way of example, each different typing result may be identified by color (or symbol or other method) in family pedigree FIG. 15A, for example, red, green and blue. Emma only shows red and green. Data about a specimen possibly related to a missing person is entered via an example one dimensional bar code per FIG. 11B. Profile data is preferably entered via an example two dimensional 2D bar code per FIG. 11C. Data for another individual, such as Emma, can be entered in the same manner. For example, no data connection is necessary (data stored may be locally on phone device for future upload to database), but, if connection is present, the data may be stored on the server or in the cloud.
FIG. 11D provides a view of an STR profile for a father family member entered as per FIG. 10A whereby a user may receive feedback that a two dimensional bar code has been correctly scanned and the corresponding data stored by comparison with the human-readable versions of the data shown in FIGS. 10A-10C.
FIG. 11E provides corresponding data for a mitochondrial DNA profile whereby similar feedback may be provided to a user that data has been retrieved and stored in memory accurately. FIG. 11F provides corresponding data for a YSTR profile typically collected for a male such as a son, father, uncle or grandfather.
FIG. 12A provides an exemplary screen for either entering a missing person, disaster victim or other unknown biological specimen data or for entering data (name, birth date and the like) for possible kin or browsing people already entered such as a mother, a father, a son, a daughter, a grandfather on a mother's side and so on, and for building a pedigree comprising a family pedigree which may be an output of an automated decision support system. FIG. 12B shows a screen for entering persons including father 1240, mother 1250, child 1260, their personal information 1270, specimen 1280, browse specimens 1290, view a pedigree 1292 and back 1294. FIG. 15A shows such as hypothetical family pedigree against which the likelihood or probability that the DNA information obtained from human remains or a person may be computed and used to assess whether this DNA is consistent with an indicated missing family member or can be excluded using, for example, a missing persons/kinship algorithm as described above.
FIG. 13 provides an exemplary screen where a father, mother and child of a family pedigree may be edited or have yet to be entered and associated with a pedigree with a screen similar to that of FIG. 12B. A user may edit a father 1310, mother 1320, add a child 1330 and the like, add personal information 1340, collect a specimen 1350, browse a specimen already entered 1360, view a pedigree chart 1370 such as FIG. 1, 3 or 15A or go Back 1380.
FIG. 14A provides an exemplary screen whereby personal data may be entered for actually or tentatively identifying an individual or a missing person by name, gender and the like including a virtual touch screen keypad or by using an actual keyboard. In a virtual keyboard, touch key ABC may select an alphabetic keyboard. Done indicates complete entry of a scrolled field such as name or birth date. A missing person box is checked only if the entered personal data relates to a missing person; otherwise, the personal data is for a member of a family pedigree or, for example, a crime victim such as Sarah Jane Smith of FIG. 15B (a rape victim). Once the data is entered, a save screen FIG. 14B appears with all completed data and one may edit the data or actuate Save.
FIG. 14C provides an exemplary screen for selecting an unknown specimen. As can be seen from FIG. 14B the specimen selected is 138827364880, which matches missing person laboratory report identifier one dimensional bar code FIG. 14D previously entered. Consequently, FIG. 14D provides a one dimensional identifier and two dimensional bar codes for each of STR, YSTR and mitochondrial DNA profile data for the missing person who may be personally identified by data entry touch screen FIG. 14A as well as human-readable versions of the data represented by the bar codes, processing date and time and analyst and, optionally, the analyst's reviewer's identifiers or names.
FIG. 15A provides a completed pedigree chart for mother, Emma, father, John, missing person John, Jr., son, John, III and daughter, Jane whereby collected profile data may be indicated by color. John, III has not been entered. Neither John or Jane are required to obtain a match given the presence of both profiles for mother, Emma and father, John as discussed above in the discussion of automatic decision support (ADS). FIG. 15C provides a corresponding Select a Person screen for an open details screen of FIG. 15B whereby one may select one of the five Murdoch family members of the pedigree of FIG. 15A including a further person Sarah Jane Smith, not shown in the pedigree chart of FIG. 15A. As discussed above, Sarah may be the victim of a rape and her DNA profile may be known. John Murdoch jr. may be a hypothetical perpetrator who is missing and has no profile data on file but may be identified using least square deconvolution from the pedigree, for example, including father, mother and daughter.
After a pedigree and associated typed family members are entered or defined, an unknown sample that has been that has been typed and entered or defined may be evaluated against the pedigree by a missing persons/kinship (MP/KN) algorithm as described above. A least squares deconvolution algorithm may be used if a mixture of DNA from two individuals is present in a sample, and/or a peak fitting algorithm, a mass spectrometry algorithm, or an expert system may be applied to perform analysis of data and transform it into a more useful form such as the DNA profiles of two likely contributors to a mixture or an accurate DNA profile. Other algorithms may come to mind to one of ordinary skill in the art for identifying a hypothetical missing person, crime victim, crime perpetrator or contributor to a DNA mixture. Such algorithms may utilize other data instead of or in addition to DNA profile data, such as SMT, facial profile, biometric data such as iris scans or blood vessel patterns in the skin, or fingerprints. Either the hypothetical pedigree is proven or it is not and that result may be returned by the MP/KN algorithm to for display to the user.
All patents, patent applications, articles and other references cited in this disclosure and identified below are expressly incorporated herein by reference as to their entire contents. Further, while several preferred embodiments and methods have been described using specific terms, such description is for illustrative purposes only, and it is to be understood that changes and variations may be made without departing from the spirit or scope of the following claims.
REFERENCES


    F. Dudbridge. A survey of Current Software for Linkage Analysis, Human Genomics, Vol 1, (1), 2004, 63-65.
    Soledad A. Fernandez, Rohan L. Fernando, Bernt Guldbrandtsen, Liviu R. Totir, Alicia L. Carriquiry, Sampling genotypes in large pedigrees with loops. Genet. Sel. Evol. 33 (2001) 337-367
    Howard D. Cash, Jonathon W. Hoyle, and Amy J. Sutton, “Development under extreme conditions: forensic bioinformatics in the wake of the World Trade Center disaster,” genecodesforensics.com/news/CashHoyleSutton.pdf, 2003
    C. Cannings, E. A. Thompson, and H. H. Skolnick. The Recursive Derivation of Likelihoods on Complex Pedigrees. Advances in Applied Probability, Vol 8, No 4. (December 1976), pp 622-625.
    E. Lander, and P Green. Construction of Multilocus Genetic Linkage Maps in Humans. PNAS, Vol 84, No. 8 (Apr. 15, 1987), pp. 2363-2367.
    R. C. Elston and J. Stewart. A General Model for the Genetic Analysis of Pedigree Data. Human Heredity 21: 523-542 (1971).
    http://linkage.rockefeller.edu/soft/

## Claims

1. A computer-implemented method of collecting first and second DNA profile data using first and second different DNA typing technologies about an unidentified remain and first and second DNA profile data of the same first and second DNA typing technologies of each of a selected first and a selected second family member genetically related to a hypothesized identity of the unidentified remain and the first and the second family members being selected via application of a family member selection rule base considering available members of the family for DNA typing, the family member data including first and second DNA profile data obtained from a specimen associated with a selected typed family member and a genetic relationship existing among each of the first and the second selected family member and the hypothesized identity of the unidentified remain, the method for determining one of a likelihood that the unidentified remain corresponds to a missing person genetically related to each of the selected first and second typed family member and that the unidentified remain can be excluded as not consistent with the first and the second selected typed family member, the method for implementation on telecommunications computer apparatus comprising a processor, an input device coupled to the processor comprising one of a touch screen of a display, a camera, a keyboard and a bar code reader for reading one and two dimension bar codes, an output device coupled to the processor comprising a display, a telecommunications interface coupled to the processor, and a memory for storing the first and the second DNA profile data obtained from said unidentified remain and said first and said second DNA profile data obtained from each specimen of said first selected and said second selected typed family member, the memory coupled to the processor, the computer-implemented method comprising:
(a) storing the first and the second DNA profile data obtained from typing said unidentified remain using the first and the second DNA typing technologies in said memory via said bar code reader input device for reading a one dimension bar code uniquely identifying the unidentified remain and two different associated two dimension bar codes for transformation to the first and the second DNA profile data of the two different typing technologies for the unidentified remain;
(b) storing a genetic relationship in said memory among said hypothesized identity of said unidentified remain of said missing person and each of said first and said second selected family member and family member and hypothesized identity identification data input via said input device;
(c) storing said first and said second DNA profile data of each of the first and the second selected family member in said memory, each of the first and the second selected family members uniquely identified by a one dimension bar code and their associated first and second DNA profile data represented by two different two dimension bar codes being uniquely transformable into first and second DNA profile data of the first and the second DNA typing technologies, said first and said second DNA profile data obtained from the specimen associated with each of the first and the second selected typed family member and storing said first and said second DNA profile data from the unidentified remain, said first and said second family member DNA profile data for each of said first and said second selected typed family member being identified with a different two-dimensional bar code representing the same first and second typing technologies used for typing the unidentified remain;
(d) displaying a family pedigree representing genetic relationships among the first and the second family members and the hypothesized identity of the unidentified remain, the family pedigree displayed on said display, said family pedigree being responsive to input received by said input device, said input causing information about a selected family member to be displayed including the stored DNA profile data associated with the selected family member, an identity being input by the input device and at least one of a first and a second type of collected DNA profile data each indicated by a predetermined color, each type comprising one of STR, Y-STR and mitochondrial DNA, at least one DNA typing technology comprising STR;
(e) evaluating the genetic consistency of said stored genetic relationship by computing a pedigree likelihood ratio for each of said stored first and said second DNA typing technologies using said DNA profile data of each of said selected typed family member with said first and said second DNA profile data of said unidentified remain and a joint likelihood ratio computed as the product of the pedigree likelihood ratios; and
(f) using the computed joint likelihood ratio to output a decision whether said stored genetic relationship and said stored first and said second DNA profile data obtained from the specimen associated with the first and second selected typed family member and said unidentified remain are consistent and outputting a decision of whether there exists one of said genetic relationship among the first and second family members and the hypothesized identity of the unidentified remain and exclusion of said genetic relationship among the first and second family members and the unidentified remain, including computing of the pedigree likelihood ratio according to evaluating the genetic consistency using a modified Elston Stewart algorithm for determining a likelihood that the unidentified remain corresponds to a missing person genetically related to each of the typed first and second second family member or can be excluded as not consistent with the typed first and second selected available family member, the modification expressing a probability of mutation.

2. The method of claim 1 where said first and said second DNA profile data each represented by first and second two dimensional bar codes is input by the bar code reader input device first reading a one-dimensional bar code uniquely identifying the unidentified remain and inputting the hypothesized identity of the unidentified remain via the input device and then reading each of the first and second different associated two-dimensional bar codes, each representing different of the first and second DNA profile data determined for each of the first and the second type of collected DNA profile data for the unidentified remain.

3. The method of claim 2 where said first and said second DNA profile data is obtained by using said a one-dimensional bar code to identify the unidentified remain and the first and second two-dimensional bar codes to identify the first and the second DNA profile data of the first and second different DNA profile technology types for the unidentified remain and accessing a particular database containing DNA profile data via a telecommunications interface.

4. The method of claim 3 where said database containing DNA profile data does not reside on the telecommunications computer apparatus and is accessed using the telecommunications interface.

5. The method of claim 1 where said DNA profile data about a selected typed first and second family member is supplemented to add first and second DNA profile data obtained from a specimen associated with a third family member selected using the family member selection rule base from available family members.

6. The method of claim 1 where the first and second DNA profile data for the unidentified remain are stored remotely in a database of first and second DNA profile data for each of a plurality of individuals.

7. The method of claim 1 where the device is hand-held.

8. The method of claim 1 where the device is a mobile telephone.

9. The method of claim 1 where the device is a tablet or a personal computer.

10. The method of claim 1 further comprising outputting a decision whether it is possible that said typed unknown biological specimen originated from or can be excluded as originating from said missing person of the pedigree.

11. The method of claim 10 further comprising computing and outputting a likelihood ratio determined from the stored data.

12. The method of claim 11 where said computing is performed remotely.

13. The method of claim 1 wherein one said DNA typing technology comprises STR DNA typing profile data of said unidentified remain and of said first and second selected typed family members at thirteen loci and consistency is determined from evaluating threshold values for pedigree likelihood ratios versus threshold values and determining a joint pedigree likelihood value versus a threshold value for consistency with DNA profile data of the unidentified remain.

14. The method of claim 1 wherein said DNA typing technology comprises at least two of STR, Y-STR, and mitochondrial DNA typing, each DNA typing technology resulting in DNA profile data separately given by a different two dimensional bar code.

15. The method of claim 14 wherein a DNA profile of said DNA typing technology comprising at least two of STR, Y-STR, and mitochondrial DNA typing is obtained using mass spectrometry, wherein Y-STR DNA typing is usable for typing a male unidentified remain or a selected male family member.

16. The method of claim 1 wherein typed available family members are typed according to combinations of two different DNA typing technologies selected from STR, Y-STR and mitochondrial DNA typing technologies, each DNA typing technology resulting in DNA profile data separately given by a different two dimensional bar code.

17. The method of claim 1 wherein said one-dimensional bar code entry each comprises a single dimensional bar code uniquely identifying one of the unidentified remain, the first selected available family member and the second selected family member having respective associated two dimensional bar codes giving DNA profile data for two different selected DNA typing technologies.

18. The method of claim 1 wherein said output device provides a graphical representation of the family pedigree including an identity of each selected typed family member and the hypothesized identity of the unidentified remain and the at least two different DNA typing technologies for each selected family member of the family pedigree by a different predetermined color.

19. The method of claim 1 further comprising automatic acquisition of location data via a global positioning system for storage with said first and said second DNA profile data of the unidentified human remain.

20. The method of claim 1 further comprising input device entry of silhouette data for the individual to be identified.

21. The method of claim 1 further comprising displaying a two dimensional bar code and associated, transformed DNA profile data of a given DNA technology type for evaluation to assure that the two dimensional bar code conforms identically with the transformed DNA profile data.

22. The method of claim 1 wherein said typed unknown biological specimen comprises a mixture of DNA of a plurality of contributor persons, the method further comprising resolving the mixture using least square deconvolution.

23. The method of claim 1 further comprising initially determining identities of the first and the second selected family members to type from the available family members using automated decision support.

24. Computer apparatus for collecting first and second DNA profile data using first and second different DNA typing technologies about an unidentified remain and first and second DNA profile data of the same first and second DNA typing technologies of each of a selected first and a selected second family member genetically related to a hypothesized identity of the unidentified remain and the first and second family members being selected via application of a family member selection rule base considering available family members of the family for DNA typing, the family member data including first and second DNA profile data obtained from a specimen associated with a selected typed family member and a genetic relationship existing among each of the first and the second selected family member and the hypothesized identity of the unidentified remain, the computer apparatus for calculating a pedigree likelihood ratio that either the unidentified remain corresponds to a missing person genetically related to each of the first and second selected typed family member or that the unidentified remain can be excluded as not consistent with each of the selected typed first and second family member, the computer apparatus configured for implementation on telecommunications computer apparatus comprising a processor, an input device coupled to the processor comprising one of a touch screen of a display, a camera, a keyboard and a bar code reader for reading one and two dimension bar codes, an output device coupled to the processor comprising a display, a telecommunications interface coupled to the processor, and a memory for storing the first and the second DNA profile data obtained from said unidentified remain and said first and said second DNA profile data obtained from the specimen of said first selected and said second selected typed family member, the memory coupled to the processor, the computer apparatus comprising:
(a) the memory for further storing the first and the second DNA profile data obtained from typing said unidentified remain via said bar code reader input device for reading a one dimension bar code uniquely identifying the unidentified remain and two different associated two dimension bar codes for transformation to the first and the second DNA profile data of the two different typing technologies for the unidentified remains;
(b) the memory for further storing a genetic relationship between said hypothesized identity of said unidentified remain of said missing person and of each of said first and second selected family member, the identities input via said input device;
(c) the memory for further storing said first and said second DNA profile data for each of the first and the second selected family member, each of said first and said second selected family members being uniquely identified by a one dimension bar code and identification data input by the input device, their associated first and second DNA profile data represented by two different two dimension bar codes being uniquely transformed into first and second DNA profile data of the first and second DNA typing technologies, said first and second DNA profile data obtained from the specimen associated with each of the first and the second selected typed family member and said first and said second DNA profile data obtained from the specimen of the unidentified remain via the same first and second DNA typing technologies, said first and said second DNA profile data of the first and the second selected family members being identified with a different two-dimensional bar code representing the same first and second typing technologies used for typing the unidentified remain;
(d) the display for displaying a family pedigree representing genetic relationships among the first and the second selected typed family members and the hypothesized identity of the unidentified remain, the family pedigree displayed on said display, said family pedigree being displayed responsive to input identification data for the family pedigree received by said input device, said input identification data causing information about a selected family member to be displayed including at least one of a family identity a hypothesized family identity and one of a first and a second type of collected DNA profile data indicated by a different predetermined color, the type comprising one of STR, Y-STR and mitochondrial DNA and at least one DNA typing technology comprising STR; and
(e) the processor for evaluating the genetic consistency of said first and said second DNA profile data of each of said selected typed family member with said first and said second DNA profile data of said unidentified remain, by calculating a pedigree likelihood ratio for at least the STR typing technology of the first and second typing technologies comprising the STR typing technology and outputting a consistency decision whether said stored genetic relationship and said stored first and second DNA profile data obtained from the specimen associated with the first and second selected typed family member and said unidentified remain are consistent with a predetermined threshold for the decision and a joint pedigree likelihood ratio for the first and second family members and the unidentified remain by multiplying the pedigree likelihood ratios together for the consistency decision with reference to a threshold, the processor calculating a pedigree likelihood ratio according to a modified Elston Stewart algorithm for determining a likelihood that the unidentified remain corresponds to a missing person genetically related to the selected typed first and second available family member or that the unidentified remain can be excluded as not consistent with each of the first and the second selected available typed family member, the modification expressing a probability of mutation.

25. The computer apparatus of claim 24 further comprising a global positioning system for automatically acquiring location data for storage with said first and second DNA profile data of the unidentified remain.

26. The apparatus of claim 24 wherein the display is for further displaying the genetic family pedigree information including the hypothesized identity of the missing person and the first and the second selected typed family member.

27. The apparatus of claim 24 wherein said input device comprises a touch screen, the touch screen comprising a virtual keyboard for inputting data for the unidentified remain.

28. The apparatus of claim 24 wherein said communications interface provides access to a DNA profile database of an internet cloud service provider.

29. The method of claim 1 wherein said first and said second DNA profile data comprises DNA profile data of one of a selected known available individual member of a family by genetic relationship to another known available individual member of the family and having a genetic family relationship to hypothesized identity of an unidentified remain of a missing person using the same two typing technologies for the first and second selected available family members and for the unidentified remain, one set of first and second DNA profile data taken via STR typing being uniquely transformable to loci data by a two dimension bar code by the input device.