A 472-SNP panel for pairwise kinship testing of second-degree relatives

Advancing human genotyping: The Infinium HTS iSelect Custom microarray panel (Rita) development study

Single Nucleotide Polymorphisms (SNPs), as the most prevalent type of variation in the human genome, play a pivotal role in influencing human traits. They are extensively utilized in diverse fields such as population genetics, forensic science, and genetic medicine. This study focuses on the 'Rita' BeadChip, a custom SNP microarray panel developed using Illumina Infinium HTS technology. Designed for high-throughput genotyping, the panel facilitates the analysis of over 4000 markers efficiently and cost-effectively. After careful clustering performed on a set of 1000 samples, an evaluation of the Rita panel was undertaken, assessing its sensitivity, repeatability, reproducibility, precision, accuracy, and resistance to contamination. The panel's performance was evaluated in various scenarios, including sex estimation and parental relationship assessment, using GenomeStudio data analysis software. Findings show that over 95 % of the custom BeadChip assay markers were successful, with better performance of transitions over other mutations, and a considerably lower success rate for Y chromosome loci. An exceptional call rate exceeding 99 % was demonstrated for control samples, even with DNA input as low as 0.781 ng. Call rates above 80 % were still obtained with DNA quantities under 0.1 ng, indicating high sensitivity and suitability for forensic applications where DNA quantity is often limited. Repeatability, reproducibility, and precision studies revealed consistency of the panel's performance across different batches and operators, with no significant deviations in call rates or genotyping results. Accuracy assessments, involving comparison with multiple available genetic databases, including the 1000 Genome Project and HapMap, denoted over 99 % concordance, establishing the Rita panel's reliability in genotyping. The contamination study revealed insights into background noise and allowed the definition of thresholds for sample quality evaluation. Multiple metrics for differentiating between negative controls and true samples were highlighted, increasing the reliability of the obtained results. The sex estimation tool in GenomeStudio proved highly effective, correctly assigning sex in all samples with autosomal loci call rates above 97 %. The parental relationship assessment of family trios highlighted the utility of GenomeStudio in identifying genotyping errors or potential Mendelian inconsistencies, promoting the application of arrays such as Rita in kinship testing. Overall, this evaluation confirms the Rita microarray as a robust, high-throughput genotyping tool, underscoring its potential in genetic research and forensic applications. With its custom content and adaptable design, it not only meets current genotyping demands but also opens avenues for further research and application expansion in the field of genetic analysis.

A novel approach of kinship determination based on the physical length of genetically shared regions of chromosomes

BACKGROUND

Determination of genetic relatedness between individuals plays a crucial role in resolving numerous civil cases involving familial relationships and in forensic investigation concerning missing persons. Short tandem repeats (STRs), known for their high degree of DNA polymorphism, have traditionally been the primary choice of DNA markers in genetic testing, but their application for kinships testing is limited to cases involving close kinship. SNPs have emerged as promising supplementary markers for kinship determination. Nevertheless, the challenging remains in discriminating between third-degree or more distant relatives, such as first cousins, using SNPs.

OBJECTIVE

To investigate a kinship analysis method for distant degree of familial relationships using high-density SNP data.

METHODS

A high-density SNP data from 337 individuals of Korean families using Affymetrix Axiom KORV1.0-96 Array was obtained for this study. SNPs were aligned by chromosomal positions, and identity-by-state (IBS) was determined, and then shared regions as consecutive SNPs with IBS of 1 or 2 were investigated. The physical lengths of these IBS segments were measured and summed them to create an Index, as a measure of kinship.

RESULTS

The kinship was determined by the physical length of shared chromosomal regions that are distinguished by each kinship. Using this method, the relationship was able be distinguished up to the fourth degree of kinship, and non-relatives were clearly distinguished from true relatives. We also found a potential for this approach to be used universally, regardless of microarray platforms for SNP genotyping and populations.

CONCLUSION

This method has a potential to determine the different degree of kinship between individuals and to distinguish non-relatives from true relatives, which can be of great help for practical applications in kinship determination.

An approach to unified formulae for likelihood ratio calculation in pairwise kinship analysis

Introduction: The likelihood ratio (LR) can be an efficient means of distinguishing various relationships in forensic fields. However, traditional list-based methods for derivation and presentation of LRs in distant or complex relationships hinder code editing and software programming. This paper proposes an approach for a unified formula for LRs, in which differences in participants' genotype combinations can be ignored for specific identification. This formula could reduce the difficulty of by-hand coding, as well as running time of large-sample-size simulation. Methods: The approach is first applied to a problem of kinship identification in which at least one of the participants is alleged to be inbred. This can be divided into two parts: i) the probability of different identical by descent (IBD) states according to the alleged kinship; and ii) the ratio of the probability that specific genotype combination can be detected assuming the alleged kinship exists between the two participants to the similar probability assuming that they are unrelated, for each state. For the probability, there are usually recognized results for common identification purposes. For the ratio, subscript letters representing IBD alleles of individual A's alleles are used to eliminate differences in genotype combinations between the two individuals and to obtain a unified formula for the ratio in each state. The unification is further simplified for identification cases in which it is alleged that both of the participants are outbred. Verification is performed to show that the results obtained with the unified and list-form formulae are equivalent. Results: A series of unified formulae are derived for different identification purposes, based on which an R package named KINSIMU has been developed and evaluated for use in large-size simulations for kinship analysis. Comparison between the package with two existing tools indicated that the unified approach presented here is more convenient and time-saving with respect to the coding process for computer applications compared with the list-based approach, despite appearing more complicated. Moreover, the method of derivation could be extended to other identification problems, such as those with different hypothesis sets or those involving multiple individuals. Conclusion: The unified approach of LR calculation can be beneficial in kinship identification field.

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.

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.

Identification of novel SNP markers for kinship analysis in the Korean population

Kinship testing using genetic markers such as short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) is crucial for forensic analysis. Although STR markers have superior discriminatory power due to their highly polymorphic properties, they have several weak points in determining extended distant or complex relationships because of high mutation rates and low success rates in degraded samples. Therefore, SNPs are regarded as promising tools in forensic science because they have low mutation rates and small amplicon sizes. Herein, we propose an SNP panel consisting of 1400 autosomal SNPs obtained from the Korean National Standard Reference Variome (KoVariome) database. To evaluate its performance, in-silico analysis was performed using whole-genome sequencing (WGS) data from 21 Korean families. Subsequently, to estimate pairwise relatedness, kinship coefficients were calculated using PLINK, and Welch's one-way ANOVA test with Games-Howell's pairwise comparison test was performed. As a result, the average kinship coefficients of first- (parent-offspring and full siblings), second- (grandparent-grandchildren and aunt/uncle-niece/nephew), and third- (first cousin and grandniece/grandnephew) degree relatives, and unrelated were 0.24, 0.11, - 0.054, and - 0.0082, respectively. Consequently, relatives (first and second degree) were distinguished from non-relatives; however, further studies are required to investigate more effective SNP markers for discriminating extended kinship. Nevertheless, the results of this study go beyond the scope of screening using the discovered 1400 SNPs in Korean families and suggest the applicability of kinship analysis in the Korean population.

Classical and Non-Classical HLA Alleles as Supplementary Markers in Indirect Kinship Parentage Testing

A civil paternity investigation involving the parents of the deceased alleged father in order to establish a family relationship is presented. On the basis of the 23 autosomal short tandem repeat (aSTR) genotyping results, conclusive proof of paternity was not achieved, as the probability of paternity (W) was calculated to 0.99988. Additional genetic data of 17 classical and non-classical human leukocyte alleles (HLA) typing by next-generation sequencing (NGS) at a high-resolution level supported the hypothesis of grandpaternity over the hypothesis of coincidental paternal obligate allele (POA) sharing (total WaSTR&HLA = 0.9999998). The present study demonstrates the utility of 17 HLA genetic markers-typing in the solution of deficiency cases of disputed parentage.

An Application of ITO Analysis in Secondary Kinship Identification

Objective

As the methods of the paternity and kinship testing have been developed, the second-degree and more distant relationships remain challenging in forensic science. Currently, the ITO method is the mainstream method to clarify the kinship between two individuals.

Methods

In this study, the ITO algorithm was used to calculate the uncle-nephew index based on 55 autosomal short tandem repeats (STRs) loci that were universally used for forensic identification. 19 STRs loci in Y chromosome were used for verification of the kinship.

Results

The cumulative uncle-nephew index between A and B was calculated to 0.993 by the analysis of the genotyping results of 21 STRs. When genotyping results of the other 34 STRs were added to the calculation algorithm, the cumulative uncle-nephew index between A and B was promoted to 227.928. Meanwhile, genotyping results of 17 Y-STRs loci showed that A and B shared the same Y-STRs haplotype that was in accord with the paternal inheritance law.

Conclusion

The biological uncle-nephew relationship between A and B are identified by applying the statistical principles and genetic technologies.

Kinship assignment with the ForenSeq™ DNA Signature Prep Kit: Sources of error in simulated and real cases

Kinship recognition between anonymous DNA samples is becoming a relevant issue in forensics, more so with the increasing number of DNA profiles in databanks. Also, NGS-based genotyping is being increasingly used in routine personal identification, to simultaneously type large numbers of markers of different kind. In the present work, we explored computationally and experimentally the performance of the ForenSeq™ DNA Signature Prep Kit in identifying the true relationship between two anonymous samples, distinguishing it from other possible relationships. We analyzed with Familias R series of 10,000 pairs with 9 different simulated relationships, corresponding to different degrees of autosomal sharing. For each pair we obtained likelihood ratios for five kinship hypotheses vs. unrelatedness, and used their ranking to identify the preferred relationship. We also typed 21 subjects from two pedigrees, representing from parent-child to 4th cousins relationships. As expected, the power for identifying the true relationship decays in the order of autosomal sharing. Parent-child and full siblings can be robustly identified against other relationships. For half-siblings the chance of reaching a significant conclusion is already small. For more distant relationships the proportion of cases correctly and significantly identified is 10% or less. Bidirectional errors in kinship attribution include the suggestion of relatedness when this does not exist (10-50%), and the suggestion of independence in pairs of individuals more than 4 generations apart (25-60%). The real cases revealed a relevant effect of genotype miscalling at some loci, which could only be partly avoided by modulating the analysis parameters. In conclusion, with the exception of first degree relatives, the kit can be useful to inform additional investigations, but does not usually provide probatory results.

Research progress on application of microhaplotype in forensic genetics

As a novel genetic marker, microhaplotype can be applied in the field of forensic genetics. Microhaplotype has the advantages of high polymorphism, low mutation rate, no stutter products and short amplification fragments. Microhaplotype can effectively detect mixture, and quantitatively analyze the contributors of mixture. DNA with severe fragmentation can be successfully genotyped by microhaplotype. It can be used as ancestry informative marker to effectively divide the global continental population according to genetic structure. Microhaplotype system can provide more information than traditional short tandem repeat and help to identify complex relationships. It can provide new ideas for tumor source identification, cell line identification and prenatal paternity testing. Here we review the applications of microhaplotype, intending to provide references for forensic practice.

BGISEQ-500RS sequencing of a 448-plex SNP panel for forensic individual identification and kinship analysis

Next generation sequencing (NGS)-based single nucleotide polymorphism (SNP) genotyping is widely used in the field of forensics. SNP genotyping data from several NGS platforms have been published, but forensic application trials of DNA nanoball sequencing platforms have been very limited. In this work, we developed a 448-plex SNP panel on the BGISEQ-500RS platform. The sequencing metrics of a total of 261 samples that were sequenced with this panel are reported in detail. The average sequencing depth was 8373 × and the average heterozygosity of the 448-plex assay was 0.85. Sensitivity analysis showed that 325 SNPs were successfully genotyped with as little as 50 pg of genomic DNA, with the mean quality score of the sequencing data above Q30. Forensic parameters were calculated based on the data of 142 unrelated Chinese Han individuals and the combined matching probability was as low as 5.21 × 10^-101. Kinship analyses based on experiments and computer simulations showed that the 448-panel was as effective as the ForenSeq™ DNA Signature Prep Kit for second-degree kinship identification, and when the two panels were merged, the related pairs were almost completely distinguished from unrelated pairs. The 448-plex SNP panel on the BGISEQ-500RS platform provides a powerful tool for forensic individual identification and kinship 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.

Forensic nanopore sequencing of STRs and SNPs using Verogen's ForenSeq DNA Signature Prep Kit and MinION

The MinION nanopore sequencing device (Oxford Nanopore Technologies, Oxford, UK) is the smallest commercially available sequencer and can be used outside of conventional laboratories. The use of the MinION for forensic applications, however, is hindered by the high error rate of nanopore sequencing. One approach to solving this problem is to identify forensic genetic markers that can consistently be typed correctly based on nanopore sequencing. In this pilot study, we explored the use of nanopore sequencing for single nucleotide polymorphism (SNP) and short tandem repeat (STR) profiling using Verogen’s (San Diego, CA, USA) ForenSeq DNA Signature Prep Kit. Thirty single-contributor samples and DNA standard material 2800 M were genotyped using the Illumina (San Diego, CA, USA) MiSeq FGx and MinION (with R9.4.1 flow cells) devices. With an optimized cutoff for allelic imbalance, all 94 identity-informative SNP loci could be genotyped reliably using the MinION device, with an overall accuracy of 99.958% (1 error among 2926 genotypes). STR typing was notably error prone, and its accuracy was locus dependent. We developed a custom-made bioinformatics workflow, and finally selected 13 autosomal STRs, 14 Y-STRs, and 4 X-STRs showing high consistency between nanopore and Illumina sequencing among the tested samples. These SNP and STR loci could be candidates for panel design for forensic analysis based on nanopore sequencing.

Parentage analysis using genome-wide high-density SNP microarray

OBJECTIVE

Parentage analysis is a technology that uses genetic methods to verify or exclude relationships between individuals. STR technology is often used in parentage analysis. We received three sets of samples from three families. Each set of samples consisted of a male individual and a female individual. Their test requirements were meant to determine whether they were a paternity relationship, a sibling relationship, or grandparent-grandchild relationship. However, only one STR locus mismatch was detected in each group. Other family members to assist in testing could not be identified; therefore, other methods were needed to assist in judgment. Using high-density SNP microarrays, we analyzed the feasibility of its application in paternity analysis.

RESULTS

A total of 180 samples were tested, including 100 unrelated samples, and 74 samples from 30 families, and six samples from three families. The data were analyzed, grouped according to the chromosome of SNP, and the mismatching rate was counted. The total mismatching rate of SNP in unrelated individuals was 8-10 times higher than that of parent-child individuals. Individuals with a total mismatch rate of more than 5.3% were defined as individuals with no kinship, and the individuals with a total mismatch rate of less than 0.6% were defined as the individuals with a parent-child relationship.

CONCLUSIONS

Through the use of high-density gene chips for analysis, we also completed an auxiliary analysis of the kinship of the three families. The gene chip is a better method for auxiliary analysis of the kinship between individuals.

Developmental validation of the MGIEasy Signature Identification Library Prep Kit, an all-in-one multiplex system for forensic applications

Analyzing genetic markers in nuclear and mitochondrial genomes is helpful in various forensic applications, such as individual identifications and kinship analyses. However, most commercial kits detect these markers separately, which is time-consuming, laborious, and more error-prone (mislabelling, contamination, ...). The MGIEasy Signature Identification Library Prep Kit (hereinafter "MGIEasy identification system"; MGI Tech, Shenzhen, China) has been designed to provide a simple, fast, and robust way to detect appropriate markers in one multiplex PCR reaction: 52 autosomal STRs, 27 X-chromosomal STRs, 48 Y-chromosomal STRs, 145 identity-informative SNPs, 53 ancestry-informative SNPs, 29 phenotype-informative SNPs, and the hypervariable regions of mitochondrial DNA (mtDNA). Here, we validated the performance of MGIEasy identification system following the guidelines of the Scientific Working Group on DNA Analysis Methods (SWGDAM), assessing species specificity, sensitivity, mixture identification, stability under non-optimal conditions (degraded samples, inhibitor contamination, and various substrates), repeatability, and concordance. Libraries prepared using MGIEasy identification system were sequenced on a MGISEQ-2000 instrument (MGI Tech). MGIEasy-derived STR, SNP, and mtDNA genotypes were highly concordant with CE-based STR genotypes (99.79%), MiSeq FGx-based SNP genotypes (99.78%), and Sanger-based mtDNA genotypes (100%), respectively. This system was strongly human-specific, resistant to four common PCR inhibitors, and reliably amplified both low quantities of DNA (as low as 0.125 ng) and degraded DNA (~ 150 nt). Most of the unique alleles from the minor contributor were detected in 1:10 male-female and male-male mixtures; some minor Y-STR alleles were even detected in 1:1000 male-female mixtures. MGIEasy also successfully directly amplified markers from blood stains on FTA cards, filter papers, and swabs. Thus, our results demonstrated that MGIEasy identification system was suitable for use in forensic analyses due to its robust and reliable performance on samples of varying quality and quantity.

A new digital approach to SNP encoding for DNA identification

Identification of individuals has become an urgent problem for mankind. In the last three decades, STR-based DNA identification has actively evolved along with traditional biometric methods. Nonetheless, single-nucleotide polymorphisms (SNPs) are now of great interest and a number of relevant SNP panels have been proposed for DNA identification. Here, a simple approach to SNP data digitization that can provide assigning a unique genetic identification number (GIN) to each person is proposed. The key points of this approach are as follows: 1) SNP data are digitized as whole 4-bit boxes in the most convenient binary format, where character "1" (YES) is assigned to revealed nucleotides, and character "0" (NO) to missing nucleotides after SNP-typing; 2) all SNPs should be considered tetra-allelic. Calculations showed that a 72-plex SNP panel is enough to provide the population with unique GINs, which can be represented in digital (binary or hexadecimal) or graphic (linear or two-dimensional) formats. Simple software for SNP data processing and GINs creation in any format was written. It is likely that the national and global GIN databases will facilitate the solution of problems related to identification of individuals or human biological materials. The proposed approach may be extended to other living organisms as well.

A Highly Polymorphic Panel Consisting of Microhaplotypes and Compound Markers with the NGS and Its Forensic Efficiency Evaluations in Chinese Two Groups

Novel genetic markers like microhaplotypes and compound markers show promising potential in forensic research. Based on previously reported single nucleotide polymorphism (SNP) and insertion/deletion (InDel) polymorphism loci, 29 genetic markers including 22 microhaplotypes and seven compound markers were identified. Genetic distributions of the 29 loci in five continental populations, Kazak and Mongolian groups in China were investigated. We found that the expected heterozygosity values of these 29 loci were >0.4 in these populations, indicating these loci were relatively high polymorphisms. Population genetic analyses of five continental populations showed that five loci displayed relatively high genetic variations among these continental populations and could be useful markers for ancestry analysis. In summary, the 29 loci displayed relatively high genetic diversities in continental populations and Chinese two groups and could be informative loci for forensic research.

The limitations of kinship determinations using STR data in ill-defined populations

The likelihood ratio (LR) method is commonly used to determine kinship in civil, criminal, or forensic cases. For the past 15 years, our research group has also applied LR to ancient STR data and obtained kinship results for collections of graves or necropolises. Although we were able to reconstruct large genealogies, some pairs of individuals showed ambiguous results. Second-degree relationships, half-sibling pairs for example, were often inconsistent with detected first-degree relationships, such as parent/child or brother/sister pairs. We therefore set about providing empirical estimations of the error rates for the LR method in living populations with STR allelic diversities comparable to that of the ancient populations we had previously studied. We collected biological samples in the field in North-Eastern Siberia and West Africa and studied more than 800 pairs of STR profiles from individuals with known relationships. Because commercial STR panels were constructed for specific regions (namely Europe and North America), their allelic makeup showed a significant deficit in diversity when compared to European populations, replicating a situation often faced in ancient DNA studies. We assessed the capacity of the LR method to confirm known relationships (effectiveness) and its capacity to detect those relationships (reliability). Concerns over the effectiveness of LR determinations are mostly an issue in forensic studies, while the reliability of the detection of kinship is an issue for the study of necropolises or other large gatherings of unidentified individuals, such as disaster victims or mass graves. We show that the application of LR to both test populations highlights specific issues (both false positives and false negatives) that prevent the confirmation of second-degree kinship or even full siblingship in small populations. Up to 29% of detected full sibling relationships were either overestimated half-sibling relationships or underestimated parent-offspring relationships. The error rate for detected half-sibling relationships was even higher, reaching 41%. Only parent-offspring pairs were reliably detected or confirmed. This implies that, in populations that are small, ill-defined, or for which the STR loci analyzed are inappropriate, an examiner might not be able to distinguish a pair of full siblings from a pair of half-siblings. Furthermore, half-sibling pairs might be overlooked altogether, an issue that is exacerbated by the common confusion, in many languages and cultures, between half-siblings and full siblings. Consequently, in the study of ancient populations, human remains of unknown origins, or poorly surveyed modern populations, we recommend a conservative approach to kinship determined by LR. Next-generation sequencing data should be used when possible, but the costs and technology involved might be prohibitive. Therefore, in potentially contentious situations or cases lacking sufficient external information, uniparental markers should be analyzed: ideally, complete mitochondrial genomes and Y-chromosome haplotypes (STR, SNP, and/or sequencing).

Development of an NGS panel containing 42 autosomal STR loci and the evaluation focusing on secondary kinship analysis

High-throughput next-generation sequencing (NGS) is a feasible technique to detect considerably more markers and simultaneously obtain length and sequence information in a single reaction. In this study, we developed an NGS panel including 42 commonly used autosomal short tandem repeats (STRs) and amelogenin on the Illumina MiSeq FGx™. Sequencing accuracy was validated by the consistency of 2800M Control DNA detected using the ForenSeq™ DNA Signature Prep Kit and Sanger sequencing. Nomenclature incompatibility was found between NGS-STR and CE-STR typing at 9 loci (D3S3045, D6S477, D7S3048, D9S925, D14S608, D17S1290, D18S535, D21S1270, GATA198B05), despite the correct sequence. The difference was caused by the two different methods of identifying motif sequence and a one-to-one correspondence can be found. We evaluated the panel by investigating consistency, sequencing sensitivity and the effectiveness of the 2nd-degree relationship identification. Herein, we present sequencing results from 58 unrelated individuals of the Hebei Han population. The total discrimination power (TDP) and cumulative probability of exclusion for trio paternity testing (CPE_trio) of the 42 NGS-STR panels reached 1-2.84 × 10^-57 and 1-9.87 × 10^-21, respectively. By family simulation and likelihood ratio (LR) calculation, this panel was shown to have effectiveness for the 2nd-degree kinship identification similar to the ForenSeq™ DNA Signature Prep Kit and certain advantages compared with it due to the relatively small number of loci. As expected, it provides new data for the development of NGS-STR typing technology.

Large-scale identification of human bone remains via SNP microarray analysis with reference SNP database

Single nucleotide polymorphisms (SNPs) are valuable markers complementary to conventional forensic short tandem repeat (STR) markers in genetic typing, with potential advantages in challenging forensic casework. With the advent of high-throughput technologies, such as microarrays and massively parallel sequencing, the use of SNP typing has now expanded to large-scale forensic applications. Herein, a forensic case is presented to demonstrate the usefulness of SNP typing in identifying large-scale human bone remains with reference database construction. A total of 402 bone remains were recovered from an island in the Jeju Province of Korea where a massive disaster occurred in 1948. The first phase of the identification process was accomplished via conventional DNA typing methods including autosomal and Y-chromosomal STR typing, and mitochondrial DNA sequencing, which resulted in the identification of 74 of 402 remains. The second phase of the identification involved the remaining 327 unidentified remains using SNP typing as a supplementary tool based on Affymetrix resequencing array. The SNP markers of 782 family members were also analyzed and a reference database was constructed for comparison. An additional 51 bone remains were identified in the second phase. SNP data obtained from the supplementary genotyping yielded additional genetic information as well as contributed to kinship testing to determine the second degrees of relationship. In addition SNPs are useful in discriminating ambiguous relationship when only STR data are available. A software program developed for SNP typing system enabled efficient kinship analysis for large-scale forensic identification. The results and the casework are described and discussed.