Dan Wen1†Hao Xing1†
Ying Liu2Jienan Li1Weifeng Qu1Wei He1Chudong Wang1Ruyi Xu1Yi Liu1
Hongtao Jia1
Lagabaiyila Zha1*- 1Department of Forensic Medicine, School of Basic Medical Sciences, Central South University, Changsha, China
- 2Xiangya Stomatological Collage, Central South University, Changsha, China
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.
1 Introduction
Kinship testing is a major area of forensic research, and often includes paternity testing and sibling testing. Paternity testing refers to the identification of suspected relationships between parent and child, and sibling testing refers to the identification of suspected relationships between brothers and (or) sisters. Siblings are called full siblings if they are of the same father and mother, and half-siblings if they are of the half-blood. Paternity testing and sibling testing are mainly applied to disputed kinship testing, kinship testing in immigration and property inheritance, and kinship testing in major disasters and accidents (Wenk, 2004). Paternity testing can be easily solved by finding genetic exclusions. But siblings share only part of their genetic material, sibling testing is full of risk and uncertainty. Moreover, if the sample is degraded, it will make kinship testing more special and complicated.
At present, the most common method for paternity testing and sibling testing is based on the short tandem repeats (STR) typing technology, which has the advantages of high sensitivity, strong identification ability, high result accuracy, and comprehensive database (Butler, 2006). However, mismatch loci are frequently observed in relatives due to high mutation rates, and degraded samples cannot be conducive to analysis due to long amplicons (Lai and Sun, 2003). Some forensic genetics experts recommend the use of single nucleotide polymorphism (SNP) genetic markers to supplement kinship testing. SNP has obvious advantages, such as low mutation rate, large number, and short length. However, due to limited genetic information, a large number of SNPs are needed to find true genetic exclusion (Pakstis et al., 2007). Combining the advantages of STR and SNP, Kidd’s research group proposed a new genetic marker-microhaplotype, which has special advantages including lower mutation rate than STR, more polymorphic than SNP, short length, and no stutter peak. These microhaplotype loci can also be used for ancestry inference, personal identification, kinship testing, and mixture sample analysis (Kidd et al., 2014). So, short and highly polymorphic microhaplotype loci may be promising genetic markers for paternity testing and sibling testing in degraded samples.
Some studies have reported the use of microhaplotype loci in kinship testing (Zhu et al., 2019a; Zhu et al., 2019b; de la Puente et al., 2020; Kureshi et al., 2020; Qu et al., 2020; Sun et al., 2020; Staadig and Tillmar, 2021; Wen et al., 2021; Wu et al., 2021; Bai et al., 2022). In 2019, Zhu published two kinship testing studies, but a limited number of microhaplotype loci were reported, which may only be suitable for paternity testing (Zhu et al., 2019a; Zhu et al., 2019b). Then, to improve the effectiveness of kinship testing, Sun (Sun et al., 2020), Kureshi (Kureshi et al., 2020), Wen (Wen et al., 2021), and Wu (Wu et al., 2021) previously reported greatly highly polymorphic microhaplotype loci for kinship testing, which obtained good kinship detection ability. But most of the microhaplotype loci in these four studies were larger than 100 bp in length, for example, the mean length of 216 bp for Sun, 123 bp for Kureshi, 215 bp for Wen, and 164 bp for Wu, of which some microhaplotype loci may not be useful for kinship testing in severely degraded samples. Meanwhile, Staadig (Staadig and Tillmar, 2021), de la Puente (de la Puente et al., 2020), Qu (Qu et al., 2020), and Bai (Bai et al., 2022) published many short microhaplotype loci smaller than 100 bp in length, but the polymorphisms of some loci were limited. To obtain sufficient effectiveness of kinship testing, a large number of loci need to be detected, which may lead to multiplex amplification difficulties, linkage disequilibrium, as well as time-consuming and labor-intensive. Therefore, it is necessary to develop a novel panel containing an appropriate number of short and highly polymorphic microhaplotype loci for paternity testing and sibling testing in degraded samples.
In our previous studies, we screened many multi-allelic SNPs (Zha et al., 2012), some of which could form microhaplotype loci with closely linked SNPs nearby. According to this phenomenon, Sun (Sun et al., 2020), Kureshi (Kureshi et al., 2020), and Wen (Wen et al., 2021) from our laboratory reported some greatly highly polymorphic microhaplotype loci for kinship testing, and Li (Zhao et al., 2022) from our laboratory reported some short microhaplotype loci for personal identification in forensic challenging samples. So, the highly polymorphic SNP-based microhaplotype loci may have high polymorphism and short length, which may be potential genetic markers for paternity testing and sibling testing in degraded samples. In addition, next-generation sequencing (NGS) is widely accepted by the forensic community. Illumina sequencing devices have high throughput and appropriate microhaplotype reading lengths, and NGS can directly determine the phase between SNP alleles (Bruijns et al., 2018). Therefore, NGS is considered to be the optimal strategy for microhaplotype genotyping, making short and highly polymorphic microhaplotype loci suitable for kinship analysis in degraded samples. In conclusion, this study intended to select an appropriate number of highly polymorphic SNP-based microhaplotype loci, detect them by the NGS technology, and evaluate their effectiveness in paternity testing and sibling testing, which were also compared to 15 STR loci.
2 Methods and materials
2.1 Sample collection and DNA extraction
A total of 96 whole blood samples were collected from unrelated Shandong Han Chinese. The collected samples were extracted using the universal Genomic DNA kit (CWBIO, China). The extracted DNA was quantified using NanoDrop™ one (Thermo Scientific, America). The 96 blood samples were named Sample1 to 96, and Sample8 was extracted twice to create a duplicate sample (Sample8-duplicate). The extracted DNA of Sample11 was incubated with 100°C hot water for 40 and 60 min, resulting in two degraded samples (Sample11-40, and Sample11-60). Two degraded samples were also genotyped using the AmpFLSTR® Identifiler® Plus PCR amplification kit (Applied Biosystems, America) and the AGCU Expressmarker 16CS PCR amplification kit (AGCU ScienTech Incorporation, China). Both kits contain the same 15 autosomal STR loci, but the AGCU Expressmarker 16CS PCR amplification kit has smaller amplicons, which is suitable for the detection of degraded samples. In addition, 2,504 individuals from 26 different populations were included in this study based on the data of the 1000 Genomes Project (Sudmant et al., 2015) (Supplementary Table S1). Written informed consent from each participant was obtained, and ethical approval was received from the Ethics Committee of Central South University (2018-S194).
2.2 Candidate loci
The candidate microhaplotype loci were screened based on the data of CHB of the 1000 Genomes Project (Sudmant et al., 2015) according to the following criteria: 1) Each microhaplotype locus contained two or more SNPs; 2) The allelic frequencies of SNPs within the same microhaplotype locus were different; 3) The length of each microhaplotype locus was smaller than 100 bp; 4) The Ae of each microhaplotype locus was larger than 3.00; 5) The heterozygosity of each microhaplotype locus was greater than 0.65; 6) The genetic distance between adjacent microhaplotype loci within the same chromosome was larger than 5 Mb. All candidate microhaplotype loci were named according to the criteria proposed by Kidd (Kidd, 2016). The details of candidate microhaplotype loci are shown in Table 1. There was a total of 36 microhaplotype loci, of which 22 loci were from Li’s study (Zhao et al., 2022), four loci from Kureshi’s study (Kureshi et al., 2020), and two loci from Wen’s study (Wen et al., 2021). Moreover, to meet the selected criteria for this study, some reported loci had SNPs deleted or SNPs added to form novel SNP combinations, which were added lower-case letters (a, b, c, … ) to the names for distinguishing them from the original combination. The other eight microhaplotype loci (loci mh01zha018, mh02zha025, mh07zha018, mh07zha027, mh10zha010, mh12zha012, mh13zha008 and mh14zha010) were firstly reported in this study.
TABLE 1. The details of the selected microhaplotype loci.
2.3 MiSeq sequencing
The multiplex amplified PCR primers for the selected loci were designed by Thermo Fisher Scientific Life Technologies. The extracted 99 DNA samples, including 96 unrelated samples and one duplicate sample and two degraded samples, were subjected to two rounds of PCR amplification (multiplex amplified PCR and index PCR) to complete the library construction, and then the constructed library was sequenced by the Illumina MiSeq sequencing platform. Reads containing adapter contamination and low-quality reads were removed from the raw data using bcl2fastq software and BBMap (version 37.75)’s BBDuk software. The clean reads were compared to the human genome (GRCh37.p13) using BWA software, and the sequencing results were analyzed using Freebayes software (Li and Durbin, 2009; Garrison and Marth, 2012). The Integrative Genomics Viewer (IGV) software was used to view the sequencing results (Thorvaldsdóttir et al., 2013). The values of GQ > = 20 and GQ > = 30 for each sample were greater than 0.99 and 0.90, respectively. Because a slight imbalance was observed between these loci, the sequencing reads were filtered by two different researchers to call the genotype for each locus.
2.4 Statistical analysis
The Log 10 values of total reads for each sample were analyzed by the histogram, and the Log 10 values of mean reads for each locus were also analyzed by the histogram. The reproducibility was analyzed by comparing the sequencing results between Sample8 and Sample8-duplicate, and the ability to detect degraded samples was evaluated by comparing the sequencing results between Sample11 and Sample11–40 and Sample11–60. The Hardy–Weinberg equilibrium was analyzed based on the exact test using a Markov chain (Guo and Thompson, 1992), and linkage equilibrium in genotypic data was analyzed based on the permutation test using the EM algorithm (Slatkin and Excoffier, 1996), which were all performed using the Arlequin version 3.5 software (Excoffier and Lischer, 2010). The forensic parameters, including allelic frequency, power of discrimination (PD), probability of exclusion (PE), and observed heterozygosity (Ho) were calculated using the modified Powerstats version 1.2 software (Zhao et al., 2003). Ae was also calculated according to the formula reported by Kidd (Kidd and Speed, 2015).
The 100,000 parent-child vs. 100,000 unrelated individual pairs, and 100,000 full-sibling vs. 100,000 unrelated individual pairs, and 100,000 half-sibling vs. 100,000 unrelated individual pairs, and 100,000 full-sibling vs. 100,000 half-sibling pairs were simulated using Families 3 software based on data of 15 STR loci, 29 microhaplotype loci and 79 microhaplotype loci, respectively (Kling et al., 2014). The allelic frequencies of 15 STR loci were from Luo’s study (Luo et al., 2020), the allelic frequencies of 29 microhaplotype loci were from our studied population, and the allelic frequencies of 79 microhaplotype loci were from CHB. The mutation rate of 10–3 and extended stepwise mutation model was applied for STR. The mutation rate of 10–8 and equal probability mutation model was applied for microhaplotype. The likelihood ratio (LR) values of the above four kinds of relationships were recorded as paternity index (PI), full-sibling index (FSI), half-sibling index (HSI) and full/half-sibling index (FHSI), separately. LR involves two alternative hypotheses (Hp and Hd), where Hp represents a true relationship (parent-child, full-sibling, or half-sibling) and Hd represents unrelated individuals. But for FHSI, Hp represents full-sibling and Hd represents half-sibling. The distributions of Log 10 of PI, FSI, HSI, and FHSI were analyzed, and the uncovered rates (UCR) were also calculated for these four kinds of relationship pairs. The UCR was calculated as the following formula: The number of simulated true relationship Hp (Hd) pairs larger (smaller) than the maximum (minimum) LR of simulated true relationship Hd (Hp) pairs/Total simulated true relationship Hp (Hd) pairs. The system power based on the data of 15 STR loci, 29 microhaplotype loci and 79 microhaplotype loci for the above four kinds of relationships simulated pairs at different threshold values (t1, t2) was also calculated, including sensitivity, specificity, positive predictive value (PPV), negative predictive value (NPV), error rate and effectiveness. When the Log 10 LR was larger than t1, the relationship Hp was supported, but when the Log 10 LR was smaller than t2, the relationship Hd was supported. The relationship was uncertain when the Log10 LR was between t1 and t2. The sensitivity was calculated by the formula: Number of relatives correctly judging as relatives/Number of relatives; the specificity was calculated by the formula: Number of non-relatives correctly judging as non-relatives/Number of non-relatives; the PPV was calculated by the formula: Number of relatives correctly judging as relatives/Number of judging as relatives; the NPV was calculated by the formula: Number of non-relatives correctly judging as non-relatives/Number of judging as non-relatives; the error rate was calculated by the formula: (Number of relatives judging as non-relatives + Number of non-relatives judging as relatives)/(Total relatives + Total non-relatives); the effectiveness was calculated by the formula: (Number of relatives correctly judging as relatives + Number of non-relatives correctly judging as non-relatives)/(Total relatives + Total non-relatives).
3 Results
3.1 The general information
The 36 microhaplotype loci were successfully sequenced in 96 unrelated samples (Sample1–96), one duplicate sample (Sample8-duplicate) and two degraded samples (Sample 11–40 and Sample11–60). These microhaplotype loci were located on 16 different chromosomes. One microhaplotype locus included five SNPs, one microhaplotype locus included four SNPs, 14 microhaplotype loci included three SNPs, and the other 20 microhaplotype loci included two SNPs. The length of these microhaplotype loci ranged from 7 to 98 bp, and the mean length was 45.19 bp. The example sequencing raw data of locus mh10zha010 for Sample8 is shown in Supplementary Figure S1, according to which it was genotyped as AC/GA. The Log 10 values of total reads for each sample are shown in Supplementary Figure S2, which were larger than 4.00. The Log 10 values for mean reads for each locus are presented in Supplementary Figure S3, and except for locus mh04zha020, the Log 10 values of other loci were greater than 2.00. The detailed genotyping profiles of 99 samples are listed in Supplementary Table S2. The genotyping profile of Sample8 was consistent with the duplicate sample (Sample8-duplicate), which indicated good reproducibility of our panel. The genotyping profiles of Sample11–40 and Sample11–60 were identical to that of Sample11, which suggested these microhaplotype loci had a good ability to detect degraded samples. However, when these two degraded samples were examined using the AmpFLSTR® Identifiler® Plus PCR amplification kit, the allele dropout was observed, and the dropout number gradually increased with increasing incubation time (Supplementary Figure S4). Even after using the AGCU Expressmarker 16CS PCR amplification kit with smaller amplicons, the allele dropout was also observed in Sample11-60 (Supplementary Figure S5). So, our microhaplotype panel may be more suitable for the detection of degraded samples than universal STR genetic markers.
3.2 The forensic parameters analysis
3.2.1 The forensic parameters based on the data of our studied population
For total of 36 microhaplotype loci based on the data of our studied population, after the Bonferroni correction (p < 0.05/36 = 0.0014), seven microhaplotype loci (mh02zha025, mh02zha033, mh07zha018, mh10zha020, mh11zha010, mh12zAha012, and mh12zha014) showed significant deviations from Hardy-Weinberg equilibrium but the other 29 microhaplotype loci did not (Supplementary Table S3). The seven microhaplotype loci with significant deviations may be affected by genotyping errors (Hosking et al., 2004; Attia et al., 2010), but for the other 29 microhaplotype loci, the signals for disequilibrium may also be undetected due to the conservativeness of Bonferroni correction (Ye et al., 2020; Graffelman and Weir, 2022). The above 29 microhaplotype loci also did not observe the significant linkage disequilibrium after the Bonferroni correction (p < 0.05/406 = 0.0001), which is presented in Supplementary Table S4. So, only 29 microhaplotype loci were included in the subsequent analysis.
The forensic parameters of 29 microhaplotype loci based on the data of our studied population are listed in Table 2. A total of 140 alleles were observed, and the locus mh04zha032a had the largest number of 13 alleles. The smallest PD value was obtained in the locus mh05zha004a (0.83), and the largest PD value was obtained in the locus mh04zha032a (0.91). The PE values had the range of 0.27 (mh14zha010) to 0.58 (mh06zha025). The combined power of discrimination (CPD) for these 29 microhaplotype loci was 1-2.96E-26, while the combined probability of exclusion (CPE) was 1-5.45E-09. The Ho values ranged from 0.58 (mh14zha010) to 0.79 (mh06zha025), and the mean Ho was 0.73. The mean Ae was 3.61, and the Ae values of 29 microhaplotype loci were all larger than 3.00. These results indicated that our microhaplotype panel had a good potential for personal identification and kinship testing.
TABLE 2. Allelic frequencies and forensic parameters of the 29 microhaplotype loci based on the dataset of 96 unrelated Shandong Han Chinese from our study.
3.2.2 The forensic parameters based on the data of the 1000 Genomes Project
The heatmap of Ae distribution of 29 microhaplotype loci in 26 populations based on the data of the 1000 Genomes Project is shown in Figure 1. The loci mh04zha032a and mh18zha010a were highly polymorphic in all 26 populations with Ae larger than 3.00. The populations on the same continent had similar Ae distributions, for example, ACB, ASW, ESN, GWD, LWK, MSL, and YRI in AFR; CLM, MXL, PEL and PUR in AMR; CDX, CHB, CHS, JPT, and KHV in EAS; CEU, FIN, GBR, IBS, and TSI in EUR; BEB, GIH, ITU, PJL, and STU in SAS. But the populations CLM and PUR of AMR were more similar to the Ae distribution of EUR. The polymorphism of EAS was higher than that of the other four continents, and the polymorphism of AFR was the worst. The CPD values of 29 microhaplotype loci in 26 populations ranged from 1-1.60E-19 (YRI) to 1-4.89E-27 (CHS), and the CPE values ranged from 1-2.62E-05 (YRI) to 1-2.28E-09 (CHS). These results suggested our microhaplotype panel was more polymorphic in EAS and can discriminated between different populations.
FIGURE 1. The heatmap of Ae distribution of 29 microhaplotype loci in 26 populations based on the data of the 1000 Genomes Project.
3.3 The pairwise kinship testing
3.3.1 The pairwise kinship testing based on 15 STR loci
The distributions of Log 10 of PI, FSI, HSI and FHSI based on the data of 15 STR loci are presented in Figure 2. For relationship parent-child vs. unrelated individuals, a slight overlap was observed, and some degree of overlap was obtained in relationship full-sibling vs. unrelated individuals. There was a significant overlap in relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The UCR values for these four kinds of relationship pairs based on the data of 15 STR loci are shown in Table 3. The UCR for relationship parent-child vs. unrelated individuals was 86.70% of true relationship Hp and 98.51% of true relationship Hd. The UCR was smaller than 60% for relationship full-sibling vs. unrelated individuals, and smaller than 10% for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The system power based on the data of 15 STR loci for relationship parent-child vs. unrelated individuals at different threshold values is listed in Figure 5. When the threshold was set as 4/−4 for relationship parent-child vs. unrelated individuals, the sensitivity, specificity, error rate and effectiveness were 74.50%, 98.99%, 0.00%, and 86.75%, separately. The sensitivity, specificity, error rate and effectiveness were 83.51%, 79.95%, 0.07%, and 81.73% at the threshold of 2/−2 for relationship full-sibling vs. unrelated individuals (Figure 6). But for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling (Figure 7 and Figure 8), the effectiveness values were about 50% and the error rates reached 1% even at the threshold of 1/−1. These 15 STR loci were only suitable for paternity testing, but the effectiveness of paternity testing at the threshold of 4/−4 was also smaller than 90%.
FIGURE 2. The distributions of Log 10 of PI, FSI, HSI, and FHSI based on the data of 15 STR loci. (A) PI; (B) FSI; (C) HSI; (D) FHSI.
TABLE 3. The UCR values of four kinds of relationship pairs using the 15 STR loci, 29 microhaplotype loci and 79 microhaplotype loci.
3.3.2 The pairwise kinship testing based on 29 microhaplotype loci
The distributions of Log 10 of PI, FSI, HSI and FHSI based on the data of 29 microhaplotype loci are presented in Figure 3. After using the 29 microhaplotype loci, the mean Log 10 LR values for four kinds of relationship pairs with true relationship Hp were larger than that of 15 STR loci, and mean Log 10 LR values for four kinds of relationship pairs with true relationship Hd were smaller than that of 15 STR loci, especially for unrelative pairs in relationship parent-child vs. unrelated individuals due to the lower mutation rate. For relationship parent-child vs. unrelated individuals, no overlap was observed, and a slight overlap was obtained in relationship full-sibling vs. unrelated individuals, which were greatly smaller than that of 15 STR loci. There was also a significant overlap in relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The UCR values for these four kinds of relationship pairs based on the data of 29 microhaplotype loci are shown in Table 3. The UCR for relationship parent-child vs. unrelated individuals was larger than 99%. The UCR was about 60% for relationship full-sibling vs. unrelated individuals. The UCR values were also smaller than 10% for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The system power based on the data of 29 microhaplotype loci for relationship parent-child vs. unrelated individuals at different threshold values is listed in Figure 5. When the threshold was set as 4/−4 for relationship parent-child vs. unrelated individuals, the sensitivity, specificity, error rate and effectiveness were 96.79%, 99.99%, 0.00% and 98.39%, separately. The sensitivity, specificity, error rate and effectiveness were 93.30%, 92.72%, 0.03%, and 93.01% at the threshold of 2/−2 for relationship full-sibling vs. unrelated individuals (Figure 6). But for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling (Figure 7 and Figure 8), the effectiveness values were 61.33% and 71.98% even at the threshold of 1/−1. These results suggested the system power of 29 microhaplotype loci was greater than that of 15 STR loci, and our microhaplotype panel had a good ability in paternity testing and full sibling testing. But for the identification of more distant kinship relationship, our system still needs to be supplemented by other loci.
FIGURE 3. The distributions of Log 10 of PI, FSI, HSI and FHSI based on the data of 29 microhaplotype loci. (A) PI; (B) FSI; (C) HSI; (D) FHSI.
3.3.3 The pairwise kinship testing based on 79 microhaplotype loci
To further improve the ability of half sibling testing and to distinguish full sibling from half sibling in degraded samples, our 29 microhaplotype loci were combined with other 50 short and highly polymorphic microhaplotype loci (Supplementary Table S5). A total of 79 short and highly polymorphic microhaplotype loci were included in the simulated kinship testing, of which 29 microhaplotype loci were reported by our study, and 9 microhaplotype loci were reported by Staading’s study (Staadig and Tillmar, 2021), and 41 microhaplotype loci were reported by the Puente’s study (de la Puente et al., 2020). For simulated kinship analysis, the allelic frequencies of 79 microhaplotype loci were from CHB, and the linkage equilibrium was assumed. The distributions of Log 10 of PI, FSI, HSI and FHSI based on the data of 79 microhaplotype loci are presented in Figure 4. For relationship parent-child vs. unrelated individuals and full-sibling vs. unrelated individuals, no overlap was observed, and a slight overlap was obtained in relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The UCR values for these four kinds of relationship pairs based on the data of 79 microhaplotype loci are shown in Table 3. The UCR values for relationship parent-child vs. unrelated individuals and full-sibling vs. unrelated individuals reached 100%. The UCR values were also about 50% for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. The system power based on the data of 79 microhaplotype loci for relationship parent-child vs. unrelated individuals and relationship full-sibling vs. unrelated individuals at different threshold values is listed in Figure 5 and Figure 6. When the threshold was set as 4/−4 for relationship parent-child vs. unrelated individuals and full-sibling vs. unrelated individuals, the sensitivity, specificity, PPV, NPV and effectiveness values were larger than 99%, and the error rate values were 0.00%. The sensitivity, specificity, error rate and effectiveness were 82.40%, 82.44%, 0.07%, and 82.42% at the threshold of 2/-2 for relationship half-sibling vs. unrelated individuals (Figure 7). For relationship full-sibling vs. half-sibling (Figure 8), the effectiveness was about 90.89% and the error rate reached 0.05% at the threshold of 2/−2. These 79 microhaplotype loci can completely distinguish the parent-child from unrelated individuals and full-sibling from unrelated individuals, and had a strong ability to identify half-sibling and distinguish full-sibling from half-sibling. The combined short and highly polymorphic microhaplotype panel may be very useful for the complex kinship analysis in degraded samples.
FIGURE 4. The distributions of Log 10 of PI, FSI, HSI and FHSI based on the data of 79 microhaplotype loci. (A) PI; (B) FSI; (C) HSI; (D) FHSI.
FIGURE 5. The system power based on the data of 15 STR loci, 29 microhaplotype loci and 79 microhaplotype loci for relationship parent-child vs. unrelated individuals at different threshold values. (A) Threshold t1 = 1/t2 = −1; (B) Threshold t1 = 2/t2 = −2; (C) Threshold t1 = 3/t2 = −3; (D) Threshold t1 = 4/t2 = −4.
FIGURE 6. The system power based on the data of 15 STR loci, 29 microhaplotype loci and 79 microhaplotype loci for relationship full-sibling vs. unrelated individuals at different threshold values. (A) Threshold t1 = 1/t2 = −1; (B) Threshold t1 = 2/t2 = −2; (C) Threshold t1 = 3/t2 = −3; (D) Threshold t1 = 4/t2 = −4.
FIGURE 7. The system power based on the data of 15 STR loci, 29 microhaplotype loci and 79 microhaplotype loci for relationship half-sibling vs. unrelated individuals at different threshold values. (A) Threshold t1 = 1/t2 = −1; (B) Threshold t1 = 2/t2 = −2; (C) Threshold t1 = 3/t2 = −3; (D) Threshold t1 = 4/t2 = −4.
FIGURE 8. The system power based on the data of 15 STR loci, 29 microhaplotype loci and 79 microhaplotype loci for relationship full-sibling vs. half-sibling at different threshold values. (A) Threshold t1 = 1/t2 = −1; (B) Threshold t1 = 2/t2 = −2; (C) Threshold t1 = 3/t2 = −3; (D) Threshold t1 = 4/t2 = −4.
4 Discussion
In this study, we developed a short and highly polymorphic microhaplotype panel containing 36 highly polymorphic SNP-based microhaplotype loci with the length 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, and no allele dropout was observed in degraded samples incubated with 100°C hot water for 40 and 60 min. The developed microhaplotype panel may be suitable for personal identification and kinship testing in degraded samples. 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 50 short and highly polymorphic microhaplotype loci reported by Staading’s study and Puente’s study, the effectiveness values were 82.42% and 90.89% at the threshold of 2/−2 for relationship half-sibling vs. unrelated individuals and full-sibling vs. half-sibling. Our developed short and highly polymorphic microhaplotype panel 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.
Although the 15 autosomal STR loci included in the AmpFlSTRTM IdentifilerTM Plus PCR Amplification Kit and the AGCU Expressmarker 16CS PCR amplification kit were still the main loci used in paternity testing and personal identification according to Luo’s study (Luo et al., 2020), Hill’s study also reported a better-powered combination of 29 autosomal STR loci with a mean Ho of 0.81 (Hill et al., 2013). The combination of 29 microhaplotype loci in our developed panel had better performance than the combination of 15 autosomal STR loci, but the microhaplotype loci in our panel should be combined with other short and highly polymorphic microhaplotype loci to achieve the performance of other kits containing a large number of STR loci. Our developed panel may be very useful for first-degree relationship testing in degraded samples, and when combined with other 50 short and highly polymorphic microhaplotype loci, may be helpful for second-degree relationship testing. But the third-degree relationship testing, such as first cousin testing, can also be observed in forensic cases. The effectiveness of relationship first-cousin vs. unrelated individuals was 11.48% at the threshold of 2/−2 after simulated kinship analysis using the Families 3 software based on the data of 79 microhaplotype loci. To simplify the simulation, the mutation rates of 79 microhaplotype loci were set to 0. The combination of 79 microhaplotype loci, including 29 microhaplotype loci in our developed panel and the other 50 short and highly polymorphic microhaplotype loci, had limited performance in third-degree relationship testing. To address these complex and distant kinship relationship analyses, more microhaplotype loci are needed. After analyzing the Ae distribution of 29 microhaplotype loci in 26 populations based on the data of the 1000 Genomes Project, it was found that five continents had different polymorphisms, and the polymorphism of EAS was higher than that of the other four continents. Therefore, the developed panel is more suitable for paternity testing and personal identification in EAS, while the construction of population-specific microhaplotype panels may also be useful for other populations. Moreover, the detection of degraded samples collected from real cases can provide deeper insight into the applicability of our panel, so we will use our panel in further studies to detect samples exposed to various degradable conditions and degraded samples collected from real cases.
Data availability statement
The original contributions presented in the study are publicly available. This data can be found here: https://www.ncbi.nlm.nih.gov/ PRJNA858268
Ethics statement
The studies involving human participants were reviewed and approved by the ethnics approval code: 2018-S194 and granted by ethics committee of Central South University. The patients/participants provided their written informed consent to participate in this study.
Author contributions
DW and HX performed the experiments and wrote the manuscript, YL, JL, WH, and WQ contributed to data interpretation and revised the whole manuscript, CW, RX, YL, and HJ helped with data acquisition and manuscript modification, LZ designed this research and modified the manuscript. All authors gave final approval for publication.
Funding
This work was supported by the National Natural Science Foundation of China (NSFC) [Grant Number 81871533] and the Natural Science Foundation of Hunan Province [Grant Number 2020JJ4779].
Acknowledgments
We thank the volunteers who contributed samples for this study.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2022.983811/full#supplementary-material
References
Attia, J., Thakkinstian, A., McElduff, P., Milne, E., Dawson, S., Scott, R. J., et al. (2010). Detecting genotyping error using measures of degree of Hardy-Weinberg disequilibrium. Stat. Appl. Genet. Mol. Biol. 9, 5. doi:10.2202/1544-6115.1463
Bai, Z., Zhang, N., Liu, J., Ding, H., Zhang, Y., Wang, T., et al. (2022). Identification of missing persons through kinship analysis by microhaplotype sequencing of single-source DNA and two-person DNA mixtures. Forensic Sci. Int. Genet. 58, 102689. doi:10.1016/j.fsigen.2022.102689
Bruijns, B., Tiggelaar, R., and Gardeniers, H. (2018). Massively parallel sequencing techniques for forensics: A review. Electrophoresis 39, 2642–2654. doi:10.1002/elps.201800082
Butler, J. M. (2006). Genetics and genomics of core short tandem repeat loci used in human identity testing. J. Forensic Sci. 51, 253–265. doi:10.1111/j.1556-4029.2006.00046.x
de la Puente, M., Phillips, C., Xavier, C., Amigo, J., Carracedo, A., Parson, W., et al. (2020). Building a custom large-scale panel of novel microhaplotypes for forensic identification using MiSeq and Ion S5 massively parallel sequencing systems. Forensic Sci. Int. Genet. 45, 102213. doi:10.1016/j.fsigen.2019.102213
Excoffier, L., and Lischer, H. E. (2010). Arlequin suite ver 3.5: A new series of programs to perform population genetics analyses under linux and windows. Mol. Ecol. Resour. 10, 564–567. doi:10.1111/j.1755-0998.2010.02847.x
Garrison, E., and Marth, G., Haplotype-based variant detection from short-read sequencing. arXiv:1207.3907 (2012).
Graffelman, J., and Weir, B. S. (2022). The transitivity of the Hardy-Weinberg law. Forensic Sci. Int. Genet. 58, 102680. doi:10.1016/j.fsigen.2022.102680
Guo, W. S., and Thompson, E. A. (1992). Performing the exact test of Hardy-Weinberg proportion for multiple alleles. Biometrics 48, 361–372. doi:10.2307/2532296
Hill, C. R., Duewer, D. L., Kline, M. C., Coble, M. D., and Butler, J. M. U. S. (2013). U.S. population data for 29 autosomal STR loci. Forensic Sci. Int. Genet. 7, e82–e83. doi:10.1016/j.fsigen.2012.12.004
Hosking, L., Lumsden, S., Lewis, K., Yeo, A., McCarthy, L., Bansal, A., et al. (2004). Detection of genotyping errors by Hardy-Weinberg equilibrium testing. Eur. J. Hum. Genet. 12, 395–399. doi:10.1038/sj.ejhg.5201164
Kidd, K. K., and Speed, W. C. (2015). Criteria for selecting microhaplotypes: mixture detection and deconvolution. Investig. Genet. 6, 1. doi:10.1186/s13323-014-0018-3
Kidd, K. K., Pakstis, A. J., Speed, W. C., Lagacé, R., Chang, J., Wootton, S., et al. (2014). Current sequencing technology makes microhaplotypes a powerful new type of genetic marker for forensics. Forensic Sci. Int. Genet. 12, 215–224. doi:10.1016/j.fsigen.2014.06.014
Kidd, K. K. (2016). Proposed nomenclature for microhaplotypes. Hum. Genomics 10, 16. doi:10.1186/s40246-016-0078-y
Kling, D., Tillmar, A. O., and Egeland, T. (2014). Familias 3 - extensions and new functionality. Forensic Sci. Int. Genet. 13, 121–127. doi:10.1016/j.fsigen.2014.07.004
Kureshi, A., Li, J., Wen, D., Sun, S., Yang, Z., and Zha, L. (2020). Construction and forensic application of 20 highly polymorphic microhaplotypes. R. Soc. Open Sci. 7, 191937. doi:10.1098/rsos.191937
Lai, Y., and Sun, F. (2003). The relationship between microsatellite slippage mutation rate and the number of repeat units. Mol. Biol. Evol. 20, 2123–2131. doi:10.1093/molbev/msg228
Li, H., and Durbin, R. (2009). Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754–1760. doi:10.1093/bioinformatics/btp324
Luo, L., Gao, H., Yao, L., Liu, H., Zhang, H., Wu, J., et al. (2020). Updated population genetic data of 15 autosomal STR loci in a Shandong Han population from East China and genetic relationships among 26 Chinese populations. Ann. Hum. Biol. 47, 472–477. doi:10.1080/03014460.2020.1749928
Pakstis, A. J., Speed, W. C., Kidd, J. R., and Kidd, K. K. (2007). Candidate SNPs for a universal individual identification panel. Hum. Genet. 121, 305–317. doi:10.1007/s00439-007-0342-2
Qu, N., Lin, S., Gao, Y., Liang, H., Zhao, H., and Ou, X. (2020). A microhap panel for kinship analysis through massively parallel sequencing technology. Electrophoresis 41, 246–253. doi:10.1002/elps.201900337
Slatkin, M., and Excoffier, L. (1996). Testing for linkage disequilibrium in genotypic data using the Expectation-Maximization algorithm. Hered. (Edinb) 76, 377–383. doi:10.1038/hdy.1996.55
Staadig, A., and Tillmar, A. (2021). Evaluation of microhaplotypes in forensic kinship analysis from a Swedish population perspective. Int. J. Leg. Med. 135, 1151–1160. doi:10.1007/s00414-021-02509-y
Sudmant, P. H., Rausch, T., Gardner, E. J., Handsaker, R. E., Abyzov, A., Huddleston, J., et al. (2015). An integrated map of structural variation in 2, 504 human genomes. Nature 526, 75–81. doi:10.1038/nature15394
Sun, S., Liu, Y., Li, J., Yang, Z., Wen, D., Liang, W., et al. (2020). Development and application of a nonbinary SNP-based microhaplotype panel for paternity testing involving close relatives. Forensic Sci. Int. Genet. 46, 102255. doi:10.1016/j.fsigen.2020.102255
Thorvaldsdóttir, H., Robinson, J. T., and Mesirov, J. P. (2013). Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief. Bioinform. 14, 178–192. doi:10.1093/bib/bbs017
Wen, D., Sun, S., Liu, Y., Li, J., Yang, Z., Kureshi, A., et al. (2021). Considering the flanking region variants of nonbinary SNP and phenotype-informative SNP to constitute 30 microhaplotype loci for increasing the discriminative ability of forensic applications. Electrophoresis 42, 1115–1126. doi:10.1002/elps.202000341
Wenk, R. E. (2004). Testing for parentage and kinship. Curr. Opin. Hematol. 11, 357–361. doi:10.1097/01.moh.0000137914.80855.8a
Wu, R., Chen, H., Li, R., Zang, Y., Shen, X., Hao, B., et al. (2021). Pairwise kinship testing with microhaplotypes: Can advancements be made in kinship inference with these markers? Forensic Sci. Int. 325, 110875. doi:10.1016/j.forsciint.2021.110875
Ye, Z., Wang, Z., and Hou, Y. (2020). Does Bonferroni correction "rescue" the deviation from Hardy-Weinberg equilibrium? Forensic Sci. Int. Genet. 46, 102254. doi:10.1016/j.fsigen.2020.102254
Zha, L., Yun, L., Chen, P., Luo, H., Yan, J., and Hou, Y. (2012). Exploring of tri-allelic SNPs using pyrosequencing and the SNaPshot methods for forensic application. Electrophoresis 33, 841–848. doi:10.1002/elps.201100508
Zhao, F., Wu, X. Y., Cai, G. Q., and Xu, C. C. (2003). The application of modified-powerstates software in forensic biostatistics. Chin. J. Forensic Med. 8, 297–298. doi:10.13618/j.issn.1001-5728.2003.05.017
Zhao, X., Fan, Y., Zeye, M. M. J., He, W., Wen, D., Wang, C., et al. (2022). A novel set of short microhaplotypes based on non-binary SNPs for forensic challenging samples. Int. J. Leg. Med. 136, 43–53. doi:10.1007/s00414-021-02719-4
Zhu, J., Chen, P., Qu, S., Wang, Y., Jian, H., Cao, S., et al. (2019a). Evaluation of the microhaplotype markers in kinship analysis. Electrophoresis 40, 1091–1095. doi:10.1002/elps.201800351
Keywords: forensic, microhaplotype, degraded samples, paternity testing, sibling testing
Citation: Wen D, Xing H, Liu Y, Li J, Qu W, He W, Wang C, Xu R, Liu Y, Jia H and Zha L (2022) The application of short and highly polymorphic microhaplotype loci in paternity testing and sibling testing of temperature-dependent degraded samples. Front. Genet. 13:983811. doi: 10.3389/fgene.2022.983811
Received: 04 July 2022; Accepted: 07 September 2022;
Published: 26 September 2022.
Edited by:
Guanglin He, Sichuan University, ChinaReviewed by:
Jan Graffelman, Universitat Politecnica de Catalunya, SpainAbhishek Singh, National Forensic Sciences University, India
Copyright © 2022 Wen, Xing, Liu, Li, Qu, He, Wang, Xu, Liu, Jia and Zha. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Lagabaiyila Zha, 40409716@qq.com
†These authors contributed equally to this work and should be considered co-first authors