Identification of missing persons through kinship analysis by microhaplotype sequencing of single-source DNA and two-person DNA mixtures

Improved individual identification in DNA mixtures of unrelated or related contributors through massively parallel sequencing

DNA mixtures are a common sample type in forensic genetics, and we typically assume that contributors to the mixture are unrelated when calculating the likelihood ratio (LR). However, scenarios involving mixtures with related contributors, such as in family murder or incest cases, can also be encountered. Compared to the mixtures with unrelated contributors, the kinship within the mixture would bring additional challenges for the inference of the number of contributors (NOC) and the construction of probabilistic genotyping models. To evaluate the influence of potential kinship on the individual identification of the person of interest (POI), we conducted simulations of two-person (2 P) and three-person (3 P) DNA mixtures containing unrelated or related contributors (parent-child, full-sibling, and uncle-nephew) at different mixing ratios (for 2 P: 1:1, 4:1, 9:1, and 19:1; for 3 P: 1:1:1, 2:1:1, 5:4:1, and 10:5:1), and performed massively parallel sequencing (MPS) using MGIEasy Signature Identification Library Prep Kit on MGI platform. In addition, in silico simulations of mixtures with unrelated and related contributors were also performed. In this study, we evaluated 1): the MPS performance; 2) the influence of multiple genetic markers on determining the presence of related contributors and inferring the NOC within the mixture; 3) the probability distribution of MAC (maximum allele count) and TAC (total allele count) based on in silico mixture profiles; 4) trends in LR values with and without considering kinship in mixtures with related and unrelated contributors; 5) trends in LR values with length- and sequence-based STR genotypes. Results indicated that multiple numbers and types of genetic markers positively influenced kinship and NOC inference in a mixture. The LR values of POI were strongly dependent on the mixing ratio. Non- and correct-kinship hypotheses essentially did not affect the individual identification of the major POI; the correct kinship hypothesis yielded more conservative LR values; the incorrect kinship hypothesis did not necessarily lead to the failure of POI individual identification. However, it is noteworthy that these considerations could lead to uncertain outcomes in the identification of minor contributors. Compared to length-based STR genotyping, using sequence-based STR genotype increases the individual identification power of the POI, concurrently improving the accuracy of mixing ratio inference using EuroForMix. In conclusion, the MGIEasy Signature Identification Library Prep kit demonstrated robust individual identification power, which is a viable MPS panel for forensic DNA mixture interpretations, whether involving unrelated or related contributors.

Application of SNPs with low minor allele frequencies in missing person identification (MPI) through kinship analysis of DNA mixtures

The need to identify a missing person (MP) through kinship analysis of DNA samples found at a crime scene has become increasingly prevalent. DNA samples from MPs can be severely degraded, contain little DNA and mixed with other contributors, which often makes it difficult to apply conventional methods in practice. This study developed a massively parallel sequencing-based panel that contains 1661 single-nucleotide polymorphisms (SNPs) with low minor allele frequencies (MAFs) (averaged at 0.0613) in the Chinese Han population, and the strategy for relationship inference from DNA mixtures comprising different numbers of contributors (NOCs) and of varying allele dropout probabilities. Based on the simulated dataset and genotyping results of 42 artificial DNA mixtures (NOC = 2-4), it was observed that the present SNP panel was sufficient for balanced mixtures when referenced to the closest relatives (parents/offspring and full siblings). When the mixture profiles suffered from dropout, incorrect assignments were markedly associated with relatedness, NOC and the dropout level. We, therefore, indicate that SNPs with low MAFs could be reliably interpreted for MP identification through the kinship analysis of complex DNA mixtures. Further studies should be extended to more possible scenarios to test the feasibility of this present approach.

Massively parallel sequencing of 74 microhaplotypes and forensic characteristics in three Chinese Sino-Tibetan populations

Microhaplotype (MH), as an emerging type of forensic genetic marker in recent years, has the potential to support multiple forensic applications, especially for mixture deconvolution and biogeographic ancestry inference. Herein, we investigated the genotype data of 74 MHs included in a novel MH panel, the Ion AmpliSeq MH-74 Plex Microhaplotype Research Panel, in three Chinese Sino-Tibetan populations (Han, Tibetan, and Yi) using the Ion Torrent semiconductor sequencing. The sequencing performance, allele frequencies, effective number of alleles (Ae), informativeness (In), and forensic parameters were subsequently estimated and calculated. In addition, principal component analysis (PCA) and structure analysis were performed to explore the population relationships among the three populations and the ancestry component distribution. Overall, this novel MH panel is robust and reliable, and has an excellent sequencing performance. The Ae values ranged from 1.0126 to 7.0855 across all samples, and 75.68 % of MHs had Ae values >2.0000. Allele frequencies at some loci varied considerably among the three studied populations, and the mean In value was 0.0195. Moreover, the genetic affinity between Tibetans and Yis was closer than that between Tibetans and Hans. The aforementioned results suggest that the Ion AmpliSeq MH-74 Plex Microhaplotype Research Panel is highly polymorphic in three investigated populations and could be used as an effective tool for human forensics. Although these 74 MHs have demonstrated the competency in continental population stratification, a higher resolution for distinguishing intracontinental subpopulations and a more comprehensive database with sufficient reference population data still remain to be accomplished.

Applications of 1993 single nucleotide polymorphism loci in forensic pairwise kinship identifications and inferences

Kinship testing plays critical roles in criminal investigations, missing person searches, civil disputes, as well as identifying disaster victims. The existing commonly used short tandem repeat (STR) loci have limited effectiveness in the identification of second-degree and more distant kinships. In this study, a total of 1993 SNP loci of 119 Chinese Han individuals from eight families were sequenced on the MGISEQ-2000RS platform. The system powers of this panel for kinship identifications were evaluated based on both the likelihood ratio (LR) and identical by state (IBS) methods. The results indicated that this panel could be used as an effective tool to kinship analyses including paternity testing, full sibling testing, second-degree kinships, and first cousin kinship analyses. Both the LR and IBS methods could be applied in distinguishing first-degree and second-degree pairs from unrelated individuals. Based on the 1993 SNP loci, LR>1000 and LR<0.001 are recommended as the thresholds of identifying first-cousin kinships from unrelated individuals, and the system power of such thresholds was 0.9470. Besides, kinship coefficients for different kinship pairs were estimated and then were used to predict the kinships for pairwise individuals. This panel performs an effective kinship inference power for the predictions of first-degree, second-degree kinships and unrelated individual pairs, while presenting low sensitivity in the prediction of first-cousin kinships.

Exploring the forensic effectiveness and population genetic differentiation by self-constructed 41 multi-InDel panel in Yunnan Zhuang group

Yunnan is one of the main residences of the Zhuang group which is one of the 55 ethnic minorities in China. At present, there are relatively few researches on population genetics and forensic science of the Yunnan Zhuang group. Therefore, this study used a self-constructed panel containing 41 multi-InDel markers to analyze the genetic polymorphisms of 173 individuals from Yunnan Zhuang group. The results indicated that these 41 multi-InDels in Yunnan Zhuang group were highly polymorphic markers expect for three markers. The cumulative match probability and combined exclusion probability values of the 40 multi-InDels (MI38 marker was excluded) were 8.0671E-26 and 0.9999995959, respectively. In addition, population genetic analyses were performed on genotyping data of 41 multi-InDel markers among the Yunnan Zhuang and 26 reference populations, revealing that the Yunnan Zhuang group was genetically close to the five populations in East Asia. According to the STRUCTURE analysis, the Yunnan Zhuang group presented similar ancestral compositions to the five populations from East Asia, and when the K value was three, the five intercontinental populations showed their different genetic structures. In conclusion, the 41 multi-InDel markers could be used as an effective tool for individual identification and paternity testing of the Zhuang group in Yunnan province, as well as for their ancestry information inference studies.

Development and evaluation of a novel panel containing 188 microhaplotypes for 2nd-degree kinship testing in the Hebei Han population

Distant kinship identification is one of the critical problems in forensic genetics. As a new type of genetic marker defined and discussed in the last decade, the microhaplotype (MH) has drawn much attention in such identification owing to its specific advantages to traditional short tandem repeat (STR) or single nucleotide polymorphism (SNP) markers. In this study, MH markers were screened step by step from the 1000 Genomes Project database, and a novel multiplex panel containing 188 MHs (in which 181 are reported the first time, while 1 was reported in a previous study and the other 6 have partial overlaps with known markers) was constructed for application in 2nd- and 3rd-degree kinship identification. Along with the construction, a novel MH nomenclature was proposed, in which the SNP position information they contained was taken into account to eliminate the possibility that the same locus was named differently interlaboratory. After a series of evaluations, the panel was shown to have good sequencing accuracy, high sensitivity, species specificity, and resistance to anti-PCR inhibitors or degradation. Population data of the 188 MHs were calculated based on the genetic information of 221 unrelated Hebei Han individuals, and the effective number of alleles (Ae) ranged from 2.0925 to 8.2634 (with an average of 2.9267). For the whole system, the cumulative matching probability (CMP), the cumulative power of exclusion in paternity testing of duos (CPEduo) and that of trios (CPEtrio) reached 2.8422 × 10^-137, 1-1.3109 × 10^-21, and 1-2.8975 × 10^-39, respectively, indicating that this panel was satisfactory for individual identification and paternity testing. Then, the efficiency of the 188 MHs in 2nd- and 3rd-degree kinship testing was studied based on 30 extended families consisting of 179 2nd-degree and 121 3rd-degree relatives, as well as simulations of 0.5 million pairs of those two kinships. The results showed that clear opinions would be given in 83.36% of 2nd-degree identifications with a false rate less than 10^-5, when the confirming and excluding thresholds of cumulative likelihood ratio (CLR) were set as 10⁴ and 10^-4, respectively. This panel is still not sufficient to solve the problem of 3rd-degree kinship identification alone, and approximately 300 or 870 MH loci would be needed in 2nd- or 3rd-degree kinship identification, respectively, to achieve a system efficiency not less than 0.99 with such a threshold set; such necessary numbers would be used only as a reference in further research.

The application of short and highly polymorphic microhaplotype loci in paternity testing and sibling testing of temperature-dependent degraded samples

Paternity testing and sibling testing become more complex and difficult when samples degrade. But the commonly used genetic markers (STR and SNP) cannot completely solve this problem due to some disadvantages. The novel genetic marker microhaplotype proposed by Kidd’s research group combines the advantages of STR and SNP and is expected to become a promising genetic marker for kinship testing in degraded samples. Therefore, in this study, we intended to select an appropriate number of highly polymorphic SNP-based microhaplotype loci, detect them by the next-generation sequencing technology, analyze their ability to detect degraded samples, calculate their forensic parameters based on the collected 96 unrelated individuals, and evaluate their effectiveness in paternity testing and sibling testing by simulating kinship relationship pairs, which were also compared to 15 STR loci. Finally, a short and highly polymorphic microhaplotype panel was developed, containing 36 highly polymorphic SNP-based microhaplotype loci with lengths smaller than 100 bp and Ae greater than 3.00, of which 29 microhaplotype loci could not reject the Hardy-Weinberg equilibrium and linkage equilibrium after the Bonferroni correction. The CPD and CPE of these 29 microhaplotype loci were 1-2.96E-26 and 1-5.45E-09, respectively. No allele dropout was observed in degraded samples incubated with 100°C hot water for 40min and 60min. According to the simulated kinship analysis, the effectiveness at the threshold of 4/−4 reached 98.39% for relationship parent-child vs. unrelated individuals, and the effectiveness at the threshold of 2/−2 for relationship full-sibling vs. unrelated individuals was 93.01%, which was greater than that of 15 STR loci (86.75% for relationship parent-child vs. unrelated individuals and 81.73% for relationship full-sibling vs. unrelated individuals). After combining our 29 microhaplotype loci with other 50 short and highly polymorphic microhaplotype loci, the effectiveness values at the threshold of 2/−2 were 82.42% and 90.89% for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The short and highly polymorphic microhaplotype panel we developed may be very useful for paternity testing and full sibling testing in degraded samples, and in combination with short and highly polymorphic microhaplotype loci reported by other researchers, may be helpful to analyze more distant kinship relationships.

State of the Art for Microhaplotypes

In recent years, the number of publications on microhaplotypes has averaged more than a dozen papers annually. Many have contributed to a significant increase in the number of highly polymorphic microhaplotype loci. This increase allows microhaplotypes to be very informative in four main areas of forensic uses of DNA: individualization, ancestry inference, kinship analysis, and mixture deconvolution. The random match Probability (RMP) can be as small as 10−100 for a large panel of microhaplotypes. It is possible to measure the heterozygosity of an MH as the effective number of alleles (Ae). Ae > 7.5 exists for African populations and >4.5 exists for Native American populations for a smaller panel of two dozen selected microhaplotypes. Using STRUCTURE, at least 10 different ancestral clusters can be defined by microhaplotypes. The Ae for a locus is also identical to the Paternity Index (PI), the measure of how informative a locus will be in parentage testing. High Ae loci can also be useful in missing persons cases. Finally, high Ae microhaplotypes allow the near certainty of seeing multiple additional alleles in a mixture of two or more individuals in a DNA sample. In summary, a panel of higher Ae microhaplotypes can outperform the standard CODIS markers.