A microhap panel for kinship analysis through massively parallel sequencing technology

Pairwise kinship inference and pedigree reconstruction using 91 microhaplotypes

Kinship inference has been a major issue in forensic genetics, and it remains to be solved when there is no prior hypothesis and the relationships between multiple individuals are unknown. In this study, we genotyped 91 microhaplotypes from 46 pedigree samples using massive parallel sequencing and inferred their relatedness by calculating the likelihood ratio (LR). Based on simulated and real data, different treatments were applied in the presence and absence of relatedness assumptions. The pedigree of multiple individuals was reconstructed by calculating pedigree likelihoods based on real pedigree samples. The results showed that the 91 MHs could discriminate pairs of second-degree relatives from unrelated individuals. And more highly polymorphic loci were needed to discriminate the pairs of second-degree or more distant relative from other degrees of relationship, but correct classification could be obtained by expanding the suspected relationship searched to other relationships with lower LR values. Multiple individuals with unknown relationships can be successfully reconstructed if they are closely related. Our study provides a solution for kinship inference when there are no prior assumptions, and explores the possibility of pedigree reconstruction when the relationships of multiple individuals are unknown.

Application of a newly constructed NGS panel with 45 X-linked microhaplotypes demonstrates the unique value of X-MH for kinship testing and mixture analysis

X-linked microhaplotypes (X-MHs) have the potential to be a valuable supplementary tool in complex kinship identification or the resolution of DNA mixtures, because they bring together the distinctive genetic pattern of X chromosomal markers and the benefits of microhaplotypes (MHs). In this study, we used the 1000 Genome database to screen and select 63 X-MHs; 18 MHs were filtered out though a batch sequencing assessment of the DNA samples collected from 112 unrelated Chinese Han individuals. The resulting 45-plex panel performed well in comprehensive assessments including repeatability, sensitivity, species specificity, resistance to PCR inhibitors or degradation, mutation rate, and accuracy in detecting DNA mixture samples. The minimum amount of DNA template that can be tested with this panel is 0.5 ng. Additionally, the alleles of the minor contributor can be accurately detected when the mixture rate is larger than 1:9 in female-male mixture or 1:19 in male-male mixture. Then, we calculated population parameters on each MH based on the allele frequency data obtained from the sequence results of the aforementioned 112 unrelated samples. Combining these parameters on each MH, it can be calculated that TDPm, TDPf, CPET, CPEDFM, CPEDFF and CNCEP3 of the 45-plex system were 1-8.99×10-13, 1-1.62×10-19, 0.9999999995, 0.9999981, 0.9955, 0.9999971 and 0.99940, respectively, indicating that the panel is capable in personal identification and parentage testing. To reveal the unique advantage of X-MHs in the analyses of complex kinship and male DNA mixture, further assessments were made. For complex kinship identification, 22 types of individual pairs with different second-degree kinship were simulated and different types of likelihood ratios (LR) were calculated for each. The results revealed that the panel can achieve accuracy of approximately 70 %∼80 % when dividing each of the three types of second-degree kinships into three or four groups. Theoretically, such sub-division cannot be done by using independent autosomal markers. For male DNA mixture analysis without suspects, the maximum likelihood ratio strategy was derived and employed in the estimation of the number of male contributors (NOMC). Simulations were conducted to verify the efficacy of the 45-plex panel in the field and to compare it with autosomal markers by assuming the 45 MHs as autosomal ones. The results showed that X-MHs can achieve higher accuracy in the estimation of NOMC than autosomal ones when the mixed males were unrelated. The results highlighted the unique value of X-linked MHs in complex kinship and male mixture analyses.

Construction and evaluation of a novel set of 90 microhaplotypes for forensic applications using NGS technology

Microhaplotypes (MHs), small sets of linked single nucleotide polymorphisms (SNPs), are becoming a valuable tool for paternity testing, personal identification and other different forensic purposes due to their advantages of both short tandem repeats (STRs) and SNPs. However, only a small part of MHs with small segments have been developed and reported so far. And the current population genetic data of MHs are still insufficient. MHs with small segments possess unique advantages in mixture deconvolution, degradation material identification, noninvasive prenatal paternity testing and even medical tumor diagnostic applications. In the present study, a set of 90 autosomal MHs whose PCR amplicon lengths are from 90-150 bp, of which 58 MHs are less than or equal to 100 bp are selected, and assembled into an amplification multiplex system optimized for Ion S5™ System for forensic application. Genetic diversity study of 90 MHs in the populations from different intercontinental regions shows that the polymorphism information content (PIC) values of 83 MHs are greater than 0.4 in populations from East Asia (EAS), and the average PIC value of 90 MHs is greater than 0.5. A total of EAS populations shows the highest cumulative match probability (CMP) and cumulative probability of exclusion (CPE) values in five intercontinental populations. The CMP and CPE values of 90 MHs in EAS are 1.1688 × 10-54 and 0.999999999998954. The informativeness for assignment (In) values of the 90 MHs are calculated based on data from five intercontinental populations, and the In values of 20 MHs have greater than 0.1, indicating that the 20 MHs are high effectiveness in distinguishing different intercontinental populations, which can be used as candidate ancestry informative markers. Further, we have studied the polymorphisms of the 90 MHs based on 224 unrelated individuals of Henan Han population, China, and obtained the frequency data of the 90 MHs. In the Henan Han population, the effective number of alleles (Ae) of the 90 MHs ranges from 1.7649 (MH45) to 3.9792 (MH50), and the Ae values of 10 MHs reach to 3.0; the Ae values of 80 MHs are greater than 2, and the average Ae value for these MHs is 2.422. The average expected heterozygosity, observed heterozygosity, PIC, matching probability, discrimination power and probability of exclusion values of 90 MHs in the Henan Han population are 0.5788, 0.5851, 0.5039, 0.2608, 0.7392 and 0.2806, respectively. The CMP value of 90 MHs in the study population is less than 10-54, and their CPE value reaches 0.999999999999999923. Moreover, the results of the depth of coverage, allele coverage ratio and noise level indicate that the 90 MH amplification system has well sequencing performance, and the sequencing results are reliable. The results indicate the 90 MHs show higher polymorphisms in the study population. The present panel can be well used in paternity testing and individual identification in the study population and even the populations from EAS.

An MPS-Based 50plex Microhaplotype Assay for Forensic DNA Analysis

Microhaplotypes (MHs) are widely accepted as powerful markers in forensic studies. They have the advantage of both short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs), with no stutter and amplification bias, short fragments and amplicons, low mutation and recombination rates, and high polymorphisms. In this study, we constructed a panel of 50 MHs that are distributed on 21 chromosomes and analyzed them using the Multiseq multiple polymerase chain reaction (multi-PCR) targeted capture sequencing protocol based on the massively parallel sequencing (MPS) platform. The sizes of markers and amplicons ranged between 11–81 bp and 123–198 bp, respectively. The sensitivity was 0.25 ng, and the calling results were consistent with Sanger sequencing and the Integrative Genomics Viewer (IGV). It showed measurable polymorphism among sequenced 137 Southwest Chinese Han individuals. No significant deviations in the Hardy–Weinberg equilibrium (HWE) and linkage disequilibrium (LD) were found at all MHs after Bonferroni correction. Furthermore, the specificity was 1:40 for simulated two-person mixtures, and the detection rates of highly degraded single samples and mixtures were 100% and 93–100%, respectively. Moreover, animal DNA testing was incomplete and low depth. Overall, our MPS-based 50-plex MH panel is a powerful forensic tool that provides a strong supplement and enhancement for some existing panels.

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.

Evaluation of the MHSeqTyper47 kit for forensically challenging DNA samples

Microhaplotypes have been highly regarded for forensic mixture DNA deconvolution because they do not experience interference from stutters in the same way as short tandem repeat markers, and they tend to be more polymorphic than single nucleotide polymorphism markers. However, forensic microhaplotype kits have not been reported. The MHSeqTyper47 kit genotypes 47 microhaplotype loci. In this study, MiSeq FGx sequencing metrics for MHSeqTyper47 were presented, and the genotyping accuracy of this kit was examined. The sensitivity of MHSeqTyper47 reached 62.5 pg, and full genotyping results were obtained from degraded DNA samples with degradation indexes ≤ 3.00. Full genotypes were obtained in the presence of 100 ng/μL tannin, 50 μM heme, 25 ng/μL humic acid, and 1.25 μg/μL indigo dye. In DNA mixture studies, a minimum of 31 loci of the minor contributor were correctly genotyped at 1:99 or 99:1 mixing ratios, with the cumulative random matching probability of these loci reaching 4.54 × 10^-25. Mixing ratios could be reliably predicted from two-donor DNA mixtures based on the loci with four called alleles. Taken together, these data showed that the MHSeqTyper47 kit was effective for forensically challenging DNA analysis.

Pairwise kinship testing with microhaplotypes: Can advancements be made in kinship inference with these markers?

Kinship testing based on genetic relatedness is one of the major tasks in forensic genetics. Although short tandem repeats (STRs) are the "gold standard" biomarkers for relationship testing, microhaplotypes (MHs) have also emerged as viable options for kinship elucidation. In this work, the kinship testing efficiency of 54 highly polymorphic MHs was studied in two extended families consisting of parent-offspring, full siblings, grandparent-grandchildren, uncle/aunt-nephew/nieces, and first cousins. In addition, ten-thousand pairs of different degrees of relationships were simulated using various datasets including 54 MHs, 27 STRs plus 94 single nucleotide polymorphisms (SNPs) that were included in the ForenSeq DNA Signature Prep Kit (ForenSeq), 54 MHs plus loci in ForenSeq, and different subsets of 417-published MHs. The panels' system effectiveness in the kinship analysis were accessed by likelihood ratio distributions. The results showed that 54 MHs could be used in first-degree relationship testing with high reliability. The effectiveness of 54 MHs was slightly lower than ForenSeq but only by a narrow margin. Both 54 MHs and ForenSeq were not sufficient for distant relationship testing, and approximately 200 microhaplotypes with an average expected heterozygosity (He) = 0.79 were enough to determine second-degree relationships, but a panel of 417 MHs with an average He = 0.72 was not sufficient to first cousins testing according to the simulation analysis. In conclusion, 54 MHs could be used to serve as supplement markers for kinship testing; and well-established STR markers plus well-performing microhaplotype markers may become collective tools in forensic applications, though an enlarged pool of forensic markers is needed for distant relationship testing.

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

The flanking region variants of nonbinary SNPs and phenotype-informative SNPs (piSNPs) have been observed, which may greatly improve the discriminative ability after constituting microhaplotype. In this study, 30 microhaplotype loci based on the nonbinary SNPs and piSNPs (shown to be related to phenotypes such as hair and eye color) were selected. Genotyping were conducted on 100 unrelated northern Han Chinese, and the 26 populations from the 1000 Genome Project were also included for comparison of populations differentiation. The simulated study was conducted for evaluating the efficiency of kinship testing. These 30 microhaplotype loci we selected had good polymorphism, with a mean effective number of alleles (Ae) of 3.46. The average Ae increase was 1.27 compared with the target SNPs. The populations from the five regions worldwide could also be distinguished using these loci. The results of kinship testing showed that these microhaplotype loci had the similar ability as 15 STR loci of AmpFlSTR^R Identifiler^R PCR Amplification Kit to identify the biological parent and a stronger ability to exclude the nonbiological parents. So, these 30 microhaplotype loci may be multifunctional for forensic application, including the ability of personal identification and kinship testing equivalent to 15 STR loci, and the power of ancestry inference for distinguishing the main intercontinental population. Moreover, our selected phenotypic microhaplotype loci may theoretically have phenotype prediction capabilities. But the phenotype prediction efficiency of these phenotypic microhaplotype loci may be worse than that of piSNPs and the detailed prediction accuracy of different populations needs to be further studied.

Identification and sequencing of 59 highly polymorphic microhaplotypes for analysis of DNA mixtures

Mixture detection remains one of the major challenges within a forensic science context. In recent years, microhaplotypes were proposed to have great potential in mixture detection, although many of them are not as polymorphic as widely used short tandem repeat (STR) markers. In this study, 59 new highly polymorphic microhaplotypes were identified and sequenced with the NextSeq 500 Sequencer. Based on the whole 1000 Genomes Project dataset, the average effective number of alleles (Ae) of the 59 microhaplotypes was 5.44, and the Ae values of 36 of these microhaplotypes were > 5.00. Their genetic variations in 187 Han Chinese individuals were evaluated. The average allele coverage ratio (ACR) of heterozygotes across all loci was 0.96 ± 0.05. The number of observed alleles varied from 4 to 23, with an average of 8.8 alleles per microhaplotype locus. The average observed heterozygosity (Ho) of 59 loci was 0.77 ± 0.05, and the Ho values of 15 of these loci were > 0.80. All loci showed high polymorphisms with a discrimination power (DP) ranging from 0.80 to 0.97, and the average DP was 0.92 ± 0.03. The analysis of simulated mixtures demonstrated that the microhaplotypes reported here were highly polymorphic and performed well in forensic DNA mixture analysis. This study not only demonstrated the applicability of microhaplotypes in mixture analysis but also provided new choices for highly polymorphic microhaplotypes because after adding the markers identified here, the number of microhaplotypes with Ae values of > 4.00 will increase from ~ 50 to ~ 110.

Evaluation of a Microhaplotype-Based Noninvasive Prenatal Test in Twin Gestations: Determination of Paternity, Zygosity, and Fetal Fraction

As a novel type of genetic marker, the microhaplotype has shown promising potential in forensic research. In the present study, we analyzed maternal plasma cell-free DNA (cfDNA) samples from twin pregnancies to validate microhaplotype-based noninvasive prenatal testing (NIPT) for paternity, zygosity, and fetal fraction (FF). Paternity was determined with the combined use of the relMix package, zygosity was evaluated by examining the presence of informative loci with two fetal genome complements, and FF was assessed through fetal allele ratios. Paternity was determined in 19 twin cases, among which 13 cases were considered dizygotic (DZ) twins based on the presence of 3~10 informative loci and the remaining 6 cases were considered monozygotic (MZ) twins because no informative locus was observed. With the fetal genomic genotypes as a reference, the accuracy of paternity and zygosity determination were confirmed by standard short tandem repeat (STR) analysis. Moreover, the lower FF, higher FF, and combined FF in each DZ plasma sample were closely related to the estimated value. This present preliminary study proposes that microhaplotype-based NIPT is applicable for paternity, zygosity, and FF determination in twin pregnancies, which are expected to be advantageous for both forensic and clinical settings.

A case of heteropaternal superfecundation identified by microhap sequencing of maternal plasma cell-free DNA: A case of HS identified by microhap sequencing

Heteropaternal superfecundation (HS) refers to the fertilization of two or more oocytes by spermatozoa from different male partners during the polyovulatory period. The present study reported a newly discovered case of HS in the 10th week of gestation, in a case of disputed paternity involving a pair of female twins and two alleged fathers (AF₁ and AF₂), based on a custom-designed microhap sequencing assay and R package relMix for data interpretation. The results suggested that the twins had different biological fathers, e.g., HS, and indicated the paternity of AF₁ in relation to one of the twins while excluding AF₂ with regard to both twins. Standard short tandem repeat (STR) analysis was employed to confirm the paternity of the heteropaternal twins. The reported case indicates that HS may occur in paternity cases with dizygotic twins, and microhap, as a novel type of highly polymorphic marker proved to be suitable for mixture deconvolution, should be able to resolve this question effectively and noninvasively at the early stage of pregnancy.

Noninvasive prenatal paternity testing by target sequencing microhaps

Microhaplotypes (i.e.,microhaps or MHs) are emerging multi-allelic markers with at least two single nucleotide polymorphisms (SNPs) within ∼ 200 bp that have alleles of the same length and do not generate stutter products. Based on massively parallel sequencing (MPS) technology, microhaps have proven applicability in forensics for different application purposes. Here we evaluate the feasibility of non-invasive prenatal paternity testing (NIPPT) with a panel of polymorphic microhap markers, using cell-free DNA (cfDNA) in the maternal circulation. A custom MPS-based assay targeting 60 microhaps was developed in our previous study. Herein, we applied the developed assay to cfDNA samples in 15 NIPPT cases in the first trimester of pregnancy (6∼13 weeks). The R package relMix was employed for data interpretation, with a regression dropout estimating model. As a result, the targeted sequencing wherein target enrichment is by hybridization capture can be effectively employed for microhap sequencing with cfDNA samples. With the combined use of relMix, the paternity of the biological fathers in 15 cases was correctly determined, with the combined paternity index (CPI) value > 10¹². Moreover, the specificity of this approach was validated by the successful paternity exclusion of 3 close relatives (father, full sibling and uncle) of the biological father in one case, and further by the significant separation in CPI distribution between the biological father and 112 unrelated males in each cases. Our results indicate that this MPS-based microhap sequencing strategy could be utilized in NIPPT. This method may contribute to developments in NIPPT and to the resolution of issues related to DNA mixtures of close relatives for specific purposes.