Forensic Y-SNP analysis beyond SNaPshot: High-resolution Y-chromosomal haplogrouping from low quality and quantity DNA using Ion AmpliSeq and targeted massively parallel sequencing

UYSD: a novel data repository accessible via public website for worldwide population frequencies of Y-SNP haplogroups

For decades, there has been scientific interest in the variation and geographic distribution of paternal lineages associated with the human Y chromosome. However, the relevant data have been dispersed across numerous publications, making it difficult to consolidate. Additionally, understanding the relationships between different variants, and the tools used to analyze them, have evolved over time, further complicating efforts to harmonize this information. The Universal Y-SNP Database (UYSD) marks a substantial advancement by providing a comprehensive and accessible platform for Y-SNP and haplogroup data from populations around the world. UYSD harmonizes diverse datasets into a unified repository, facilitating the exploration of global Y-chromosomal variation. The platform handles data generated with both high- and low-throughput technology and is compatible with the automated analysis software tool, Yleaf v3. Key functionalities include the ability to: i) visualize haplogroup distributions on an interactive world map, ii) estimate haplogroup frequencies in geographic regions with sparse data through interpolation, and iii) display detailed phylogenetic trees of Y-chromosomal haplogroups. Currently, UYSD encompasses data from over 6,600 males across 27 populations. This dataset largely aligns with known global Y-haplogroup patterns, but also reveals unexplored finer-scale geographic variations. While the present dataset is largely European-centered, UYSD is designed for ongoing expansion by the scientific community, aiming to include more global data and higher-resolution population sequencing data. The platform thus offers valuable insights into human genetic diversity and migration patterns, serving several fields of research such as: human population genetics, genetic anthropology, ancient DNA analysis and forensic genetics.

Forensic DNA phenotyping: a review on SNP panels, genotyping techniques, and prediction models

In the past few years, forensic DNA phenotyping has attracted a strong interest in the forensic research. Among the increasing publications, many have focused on testing the available panels to infer biogeographical ancestry on less represented populations and understanding the genetic mechanisms underlying externally visible characteristics. However, there are currently no publications that gather all the existing panels limited to forensic DNA phenotyping and discuss the main technical limitations of the technique. In this review, we performed a bibliographic search in Scopus database of phenotyping-related literature, which resulted in a total of 48, 43, and 15 panels for biogeographical ancestry, externally visible characteristics, and both traits inference, respectively. Here we provide a list of commercial and non-commercial panels and the limitations regarding the lack of harmonization in terms of terminology (i.e., categorization and measurement of traits) and reporting, the lack of genetic knowledge and environment influence to select markers and develop panels, and the debate surrounding the selection of genotyping technologies and prediction models and algorithms. In conclusion, this review aims to be an updated guide and to present an overview of the current related literature.

YHSeqY3000 panel captures all founding lineages in the Chinese paternal genomic diversity database

BACKGROUND

The advancements in second-/third-generation sequencing technologies, alongside computational innovations, have significantly enhanced our understanding of the genomic structure of Y-chromosomes and their unique phylogenetic characteristics. These researches, despite the challenges posed by the lack of population-scale genomic databases, have the potential to revolutionize our approach to high-resolution, population-specific Y-chromosome panels and databases for anthropological and forensic applications.

OBJECTIVES

This study aimed to develop the highest-resolution Y-targeted sequencing panel, utilizing time-stamped, core phylogenetic informative mutations identified from high-coverage sequences in the YanHuang cohort. This panel is intended to provide a new tool for forensic complex pedigree search and paternal biogeographical ancestry inference, as well as explore the general patterns of the fine-scale paternal evolutionary history of ethnolinguistically diverse Chinese populations.

RESULTS

The sequencing performance of the East Asian-specific Y-chromosomal panel, including 2999-core SNP variants, was found to be robust and reliable. The YHSeqY3000 panel was designed to capture the genetic diversity of Chinese paternal lineages from 3500 years ago, identifying 408 terminal lineages in 2097 individuals across 41 genetically and geographically distinct populations. We identified a fine-scale paternal substructure that was correlating with ancient population migrations and expansions. New evidence was provided for extensive gene flow events between minority ethnic groups and Han Chinese people, based on the integrative Chinese Paternal Genomic Diversity Database.

CONCLUSIONS

This work successfully integrated Y-chromosome-related basic genomic science with forensic and anthropological translational applications, emphasizing the necessity of comprehensively characterizing Y-chromosome genomic diversity from genomically under-representative populations. This is particularly important in the second phase of our population-specific medical or anthropological genomic cohorts, where dense sampling strategies are employed.

Comparison of commercial targeted amplicon sequencing assays for human remains identification casework

Targeted amplicon sequencing (TAS) facilitates the genotyping of forensically informative single nucleotide polymorphisms (SNPs) using massively parallel sequencing (MPS). For human remains identification, where any extracted DNA is likely to be degraded, TAS may succeed when short tandem repeat (STR) profiling using capillary electrophoresis fails. Further, as well as yielding identity information, SNPs can provide information about ancestry, phenotype, kinship and paternal lineage (Y chromosome haplotypes). Two TAS platforms were compared in this study: Ion AmpliSeq™ panels coupled with Ion Torrent sequencing on an Ion GeneStudio™ S5 Plus System, manufactured by Thermo Fisher Scientific, and the ForenSeq® Kintelligence Kit coupled with Illumina sequencing on the MiSeq FGx® Sequencing System, manufactured by QIAGEN. Four Ion AmpliSeq™ panels (Precision ID Identity, Precision ID Ancestry, DNA Phenotyping and HID Y-SNP) share 177 SNPs with the ForenSeq® Kintelligence Kit and all five were used to profile the DNA extracted from the petrous part of the temporal bone from six skeletonised cadavers. Of the 6 × 177 = 1,062 SNP genotype comparisons, 1,055 (99%) were concordant between the Ion AmpliSeq™ panels and Kintelligence Kit. Of the seven (< 1%) non-concordant SNPs, only three of them (0.3%) would have resulted in erroneous genotypes being reported as a result of allele dropout by either assay, using our optimised relative variant frequency windows for allele calling. We conclude that both the Ion AmpliSeq™ panels and the ForenSeq® Kintelligence Kit were suitable for TAS applied to the human remains in this study.

Ancestry analysis using a self-developed 56 AIM-InDel loci and machine learning methods

Insertion/deletion (InDel) polymorphisms can be used as one of the ancestry-informative markers in ancestry analysis. In this study, a self-developed panel consisting of 56 ancestry-informative InDels was used to investigate the genetic structures and genetic relationships between Chinese Inner Mongolia Manchu group and 26 reference populations. The Inner Mongolia Manchu group was closely related in genetic background to East Asian populations, especially the Han Chinese in Beijing. Moreover, populations from northern and southern East Asia displayed obvious variations in ancestral components, suggesting the potential value of this panel in distinguishing the populations from northern and southern East Asia. Subsequently, four machine learning models were performed based on the 56 AIM-InDel loci to evaluate the performance of this panel in ancestry prediction. The random forest model presented better performance in ancestry prediction, with 91.87% and 99.73% accuracy for the five and three continental populations, respectively. The individuals of the Inner Mongolia Manchu group were assigned to the East Asian populations by the random forest model, and they exhibited closer genetic affinities with northern East Asian populations. Furthermore, the random forest model distinguished 87.18% of the Inner Mongolia Manchus from the East Asian populations, suggesting that the random forest model based on the 56 ancestry-informative InDels could be a potential tool for ancestry analysis.

Application of Targeted Y-Chromosomal Capture Enrichment to Increase the Resolution of Native American Haplogroup Q

Y-chromosomal haplogroups and the Y-SNPs defining them are relevant for the exploration of male lineages, inference of paternal ancestry, and reconstruction of migration pathways, to name a few. Currently, over 300,000 Y-SNPs have been reported, defining 20 main haplogroups. However, ascertainment bias in the investigations has led to some haplogroups being overlooked, which hinders a representative depiction of certain populations and their migration events. For migration pattern analyses of the first settlers of the Americas, the Native American main founding lineage Q-M3 needs to be further investigated to allow clear genetic differentiation of individuals of different ethnogeographic origins. To increase the resolution within this haplogroup, a total of 7.45 Mb of the Y chromosome of 59 admixed South Americans of haplogroup Q was targeted for sequencing using hybridization capture enrichment. Data were combined with 218 publicly available sequences of Central and South Americans of haplogroup Q. After rigorous data processing, variants not meeting the quality criteria were excluded and 4128 reliable Y-SNPs were reported. A total of 2224 Y-SNPs had previously unknown positions in the phylogenetic tree, and 1291 of these are novel. The phylogenetic relationships between the Y-SNPs were established using the software SNPtotree in order to report a redesigned phylogenetic tree containing 300 branches, defined by 3400 Y-SNPs. The new tree introduces 117 previously undescribed branches and is the most comprehensive phylogenetic tree of the Native American haplogroup Q lineages to date. The 214 sequences were assigned to 135 different low- to high-resolution branches, while in the previous phylogenetic tree, only 195 sequences could be sorted into 14 low-resolution branches with the same quality criteria. The improved genetic differentiation of subhaplogroup Q-M3 has a great potential to resolve migration patterns of Native Americans.

Development and validation of YARN: A novel SE-400 MPS kit for East Asian paternal lineage analysis

Y-chromosomal short tandem repeat polymorphisms (Y-STRs) and Y-chromosomal single nucleotide polymorphisms (Y-SNPs) are valuable genetic markers used in paternal lineage identification and population genetics. Currently, there is a lack of an effective panel that integrates Y-STRs and Y-SNPs for studying paternal lineages, particularly in East Asian populations. Hence, we developed a novel Y-chromosomal targeted panel called YARN (Y-chromosome Ancestry and Region Network) based on multiplex PCR and a single-end 400 massive parallel sequencing (MPS) strategy, consisting of 44 patrilineage Y-STRs and 260 evolutionary Y-SNPs. A total of 386 reactions were validated for the effectiveness and applicability of YARN according to SWGDAM validation guidelines, including sensitivity (with a minimum input gDNA of 0.125 ng), mixture identification (ranging from 1:1-1:10), PCR inhibitor testing (using substances such as 50 μM hematin, 100 μM hemoglobin, 100 μM humic acid, and 2.5 mM indigo dye), species specificity (successfully distinguishing humans from other animals), repeatability study (achieved 100% accuracy), and concordance study (with 99.91% accuracy for 1121 Y-STR alleles). Furthermore, we conducted a pilot study using YARN in a cohort of 484 Han Chinese males from Huaiji County, Zhaoqing City, Guangdong, China (GDZQHJ cohort). In this cohort, we identified 52 different Y-haplogroups and 73 different surnames. We found weak to moderate correlations between the Y-haplogroups, Chinese surnames, and geographical locations of the GDZQHJ cohort (with λ values ranging from 0.050 to 0.340). However, when we combined two different categories into a new independent variable, we observed stronger correlations (with λ values ranging from 0.617 to 0.754). Overall, the YARN panel, which combines Y-STR and Y-SNP genetic markers, meets forensic DNA quality assurance guidelines and holds potential for East Asian geographical origin inference and paternal lineage analysis.

Simultaneous sequencing of 102 Y-STRs on Ion Torrent ™ GeneStudio ™ S5 System

The Precision ID NGS System from Thermo Fisher Scientific is a mainstream next-generation sequencing (NGS) platform used in forensic laboratories to detect almost all commonly used forensic markers, except for Y-chromosomal short tandem repeats (Y-STRs). This study aimed to: 1) develop a Y-STR panel compatible with the automatic workflow of the NGS system using Ion AmpliSeq Technology, 2) evaluate the panel performance following the SWGDAM guidelines, and 3) explore the possibility of using a combination workflow to detect autosomal STRs and Y-STRs (AY-STR NGS workflow). The GrandFiler Y-STR Panel was successfully designed using the 'separating' and 'merging' strategies, including 102 Y-STRs and Amelogenin with an average amplicon length of 133 bp. It is a mega Y-STR multiplex system in which up to 16 samples can be sequenced simultaneously on an Ion 530 ™ Chip. Developmental validation studies of the performance of the NGS platform, species specificity, reproducibility, concordance, sensitivity, degraded samples, case-type samples, and mixtures were conducted to unequivocally determine whether the GrandFiler Y-STR Panel is suitable for real scenarios. The newly developed Y-STR panel showed compelling run metrics and NGS performance, including 92.47% bases with ≥ Q20, 91.80% effective reads, 2106 × depth of coverage (DoC), and 97.09% inter-locus balance. Additionally, it showed high specificity for human males and 99.40% methodological and bioinformatical concordance, generated complete profiles at ≥ 0.1 ng input DNA, and recovered more genetic information from severely degraded and diverse case samples. Although the outcome when used on mixtures was not as expected, more genetic information was obtained compared to that from capillary electrophoresis (CE) methods. The AY-STR NGS workflow was established by combining the GrandFiler Y-STR Panel with the Precision ID GlobalFiler ™ NGS STR Panel v2 at a 2:1 concentration ratio. The combination workflow on NGS performance, reproducibility, concordance, and sensitivity was as stable as the single Y-STR NGS workflow, providing more options for forensic scientists when dealing with different case scenarios. Overall, the GrandFiler Y-STR Panel was confirmed as the first to effectively detect a large number of Y-STR markers on the Precision ID NGS System, which is compatible with 51 Y-STRs in commercial CE kits and 51 Y-STRs in commercial NGS kits and the STRBase. The panel is as robust, reliable, and sensitive as current CE/NGS kits, and is suitable for solving real cases, especially for severely degraded samples (degradation index > 10).

Evolutionary profiles and complex admixture landscape in East Asia: New insights from modern and ancient Y chromosome variation perspectives

Human Y-chromosomes are characterized by nonrecombination and uniparental inheritance, carrying traces of human history evolution and admixture. Large-scale population-specific genomic sources based on advanced sequencing technologies have revolutionized our understanding of human Y chromosome diversity and its anthropological and forensic applications. Here, we reviewed and meta-analyzed the Y chromosome genetic diversity of modern and ancient people from China and summarized the patterns of founding lineages of spatiotemporally different populations associated with their origin, expansion, and admixture. We emphasized the strong association between our identified founding lineages and language-related human dispersal events correlated with the Sino-Tibetan, Altaic, and southern Chinese multiple-language families related to the Hmong-Mien, Tai-Kadai, Austronesian, and Austro-Asiatic languages. We subsequently summarize the recent advances in translational applications in forensic and anthropological science, including paternal biogeographical ancestry inference (PBGAI), surname investigation, and paternal history reconstruction. Whole-Y sequencing or high-resolution panels with high coverage of terminal Y chromosome lineages are essential for capturing the genomic diversity of ethnolinguistically diverse East Asians. Generally, we emphasized the importance of including more ethnolinguistically diverse, underrepresented modern and spatiotemporally different ancient East Asians in human genetic research for a comprehensive understanding of the paternal genetic landscape of East Asians with a detailed time series and for the reconstruction of a reference database in the PBGAI, even including new technology innovations of Telomere-to-Telomere (T2T) for new genetic variation discovery.

Global and Local Ancestry and its Importance: A Review

The fastest way to significantly change the composition of a population is through admixture, an evolutionary mechanism. In animal breeding history, genetic admixture has provided both short-term and long-term advantages by utilizing the phenomenon of complementarity and heterosis in several traits and genetic diversity, respectively. The traditional method of admixture analysis by pedigree records has now been replaced greatly by genome-wide marker data that enables more precise estimations. Among these markers, SNPs have been the popular choice since they are cost-effective, not so laborious, and automation of genotyping is easy. Certain markers can suggest the possibility of a population's origin from a sample of DNA where the source individual is unknown or unwilling to disclose their lineage, which are called Ancestry-Informative Markers (AIMs). Revealing admixture level at the locus-specific level is termed as local ancestry and can be exploited to identify signs of recent selective response and can account for genetic drift. Considering the importance of genetic admixture and local ancestry, in this mini-review, both concepts are illustrated, encompassing basics, their estimation/identification methods, tools/software used and their applications.

Unraveling the paternal genetic structure and forensic traits of the Hui population in Liaoning Province, China using Y-chromosome analysis

The Hui people are the second-largest ethnic minority in China, and they are distributed throughout the country. A previous study explored the paternal genetic structure of the Hui population in nine different regions of China, but it overlooked the Liaoning province. In this study, we examined the paternal genetic makeup and forensic traits of the Hui population in Liaoning province by analyzing 157 Y-chromosome single nucleotide polymorphisms (Y-SNPs) and 26 short tandem repeats (Y-STRs). We successfully genotyped 282 unrelated male individuals from the Hui population of Liaoning province using the SNaPshot® single base extension assay and Goldeneye™ Y26 system kit (PEOPLESPOT R&D, Beijing, China). The results revealed high haplotypic diversity (0.9998) and identified 46 terminal haplogroups for the Hui population. Additional analyses, such as heat maps, principal component analysis (PCA), genetic distance (FST), Multidimensional scaling (MDS) analysis, and median-joining network (MJ) analysis, showed that the Hui population could be classified into three groups: Northwest Hui populations (NWH), including Liaoning, Xinjiang, Qinghai, Gansu, Ningxia, Shaanxi, and Henan; Hui populations from Sichuan and Shandong (SSH); and Yunnan Hui populations (YNH). Pairwise genetic distance (Rst) comparisons with other Chinese populations revealed that the Hui population displayed genetic affinity with the Han population. The comprehensive understanding of the Hui population in Liaoning province, explored by Y-SNPs and Y-STRs, can be utilized to interpret their genetic structure and enhance the accuracy of forensic databases.

Developmental validation of a high-resolution panel genotyping 639 Y-chromosome SNP and InDel markers and its evolutionary features in Chinese populations

Uniparental-inherited haploid genetic marker of Y-chromosome single nucleotide polymorphisms (Y-SNP) have the power to provide a deep understanding of the human evolutionary past, forensic pedigree, and bio-geographical ancestry information. Several international cross-continental or regional Y-panels instead of Y-whole sequencing have recently been developed to promote Y-tools in forensic practice. However, panels based on next-generation sequencing (NGS) explicitly developed for Chinese populations are insufficient to represent the Chinese Y-chromosome genetic diversity and complex population structures, especially for Chinese-predominant haplogroup O. We developed and validated a 639-plex panel including 633 Y-SNPs and 6 Y-Insertion/deletions, which covered 573 Y haplogroups on the Y-DNA haplogroup tree. In this panel, subgroups from haplogroup O accounted for 64.4% of total inferable haplogroups. We reported the sequencing metrics of 354 libraries sequenced with this panel, with the average sequencing depth among 226 individuals being 3,741×. We illuminated the high level of concordance, accuracy, reproducibility, and specificity of the 639-plex panel and found that 610 loci were genotyped with as little as 0.03 ng of genomic DNA in the sensitivity test. 94.05% of the 639 loci were detectable in male-female mixed DNA samples with a mix ratio of 1:500. Nearly all of the loci were genotyped correctly when no more than 25 ng/μL tannic acid, 20 ng/μL humic acid, or 37.5 μM hematin was added to the amplification mixture. More than 80% of genotypes were obtained from degraded DNA samples with a degradation index of 11.76. Individuals from the same pedigree shared identical genotypes in 11 male pedigrees. Finally, we presented the complex evolutionary history of 183 northern Chinese Hans and six other Chinese populations, and found multiple founding lineages that contributed to the northern Han Chinese gene pool. The 639-plex panel proved an efficient tool for Chinese paternal studies and forensic applications.

Forensic biogeographical ancestry inference: recent insights and current trends

BACKGROUND

As a powerful complement to the paradigmatic DNA profiling strategy, biogeographical ancestry inference (BGAI) plays a significant part in human forensic investigation especially when a database hit or eyewitness testimony are not available. It indicates one's biogeographical profile based on known population-specific genetic variations, and thus is crucial for guiding authority investigations to find unknown individuals. Forensic biogeographical ancestry testing exploits much of the recent advances in the understanding of human genomic variation and improving of molecular biology.

OBJECTIVE

In this review, recent development of prospective ancestry informative markers (AIMs) and the statistical approaches of inferring biogeographic ancestry from AIMs are elucidated and discussed.

METHODS

We highlight the research progress of three potential AIMs (i.e., single nucleotide polymorphisms, microhaplotypes, and Y or mtDNA uniparental markers) and discuss the prospects and challenges of two methods that are commonly used in BGAI.

CONCLUSION

While BGAI for forensic purposes has been thriving in recent years, important challenges, such as ethics and responsibilities, data completeness, and ununified standards for evaluation, remain for the use of biogeographical ancestry information in human forensic investigations. To address these issues and fully realize the value of BGAI in forensic investigation, efforts should be made not only by labs/institutions around the world independently, but also by inter-lab/institution collaborations.

Benchmarking of human Y-chromosomal haplogroup classifiers with whole-genome and whole-exome sequence data

In anthropological, medical, and forensic studies, the nonrecombinant region of the human Y chromosome (NRY) enables accurate reconstruction of pedigree relationships and retrieval of ancestral information. Using high-throughput sequencing (HTS) data, we present a benchmarking analysis of command-line tools for NRY haplogroup classification. The evaluation was performed using paired Illumina data from whole-genome sequencing (WGS) and whole-exome sequencing (WES) experiments from 50 unrelated donors. Additionally, as a validation, we also used paired WGS/WES datasets of 54 individuals from the 1000 Genomes Project. Finally, we evaluated the tools on data from third-generation HTS obtained from a subset of donors and one reference sample. Our results show that WES, despite typically offering less genealogical resolution than WGS, is an effective method for determining the NRY haplogroup. Y-LineageTracker and Yleaf showed the highest accuracy for WGS data, classifying precisely 98% and 96% of the samples, respectively. Yleaf outperforms all benchmarked tools in the WES data, classifying approximately 90% of the samples. Yleaf, Y-LineageTracker, and pathPhynder can correctly classify most samples (88%) sequenced with third-generation HTS. As a result, Yleaf provides the best performance for applications that use WGS and WES. Overall, our study offers researchers with a guide that allows them to select the most appropriate tool to analyze the NRY region using both second- and third-generation HTS data.

NuMY-A qPCR Assay Simultaneously Targeting Human Autosomal, Y-Chromosomal, and Mitochondrial DNA

The accurate quantification of DNA in forensic samples is of utmost importance. These samples are often present in limited amounts; therefore, it is indicated to use the appropriate analysis route with the optimum DNA amount (when possible). Also, DNA quantification can inform about the degradation stage and therefore support the decision on which downstream genotyping method to use. Consequently, DNA quantification aids in getting the best possible results from a forensic sample, considering both its DNA quantity and quality limitations. Here, we introduce NuMY, a new quantitative real-time PCR (qPCR) method for the parallel quantification of human nuclear (n) and mitochondrial (mt) DNA, assessing the male portion in mixtures of both sexes and testing for possible PCR inhibition. NuMY is based on previous work and follows the MIQE guidelines whenever applicable. Although quantification of nuclear (n)DNA by simultaneously analyzing autosomal and male-specific targets is available in commercial qPCR kits, tools that include the quantification of mtDNA are sparse. The quantification of mtDNA has proven relevant for samples with low nDNA content when conventional DNA fingerprinting techniques cannot be followed. Furthermore, the development and use of new massively parallel sequencing assays that combine multiple marker types, i.e., autosomal, Y-chromosomal, and mtDNA, can be optimized when precisely knowing the amount of each DNA component present in the input sample. For high-quality DNA extracts, NuMY provided nDNA results comparable to those of another quantification technique and has also proven to be a reliable tool for challenging, forensically relevant samples such as mixtures, inhibited, and naturally degraded samples.

Evaluation of Library Preparation Workflows and Applications to Different Sample Types Using the PowerSeq® 46GY System with Massively Parallel Sequencing

This project evaluated the prototype PowerSeq^® 46GY System using donor DNA and casework-type samples. The goal of this study was to determine whether modifications to the manufacturer's protocol could increase read coverage and improve sample results. Buccal and casework-type libraries were prepared using the TruSeq^® DNA PCR-Free HT kit or the KAPA HyperPrep kit. Both kits were evaluated unmodified, and by substituting AMPure^® XP beads for the beads of the most optimal kit. Two qPCR kits, the PowerSeq^® Quant MS System and KAPA Library Quantification Kit, were also evaluated along with a KAPA size-adjustment workbook, which was compared as a third quantification method. Libraries were sequenced using the MiSeq^® FGx and data were analyzed with STRait Razor. Results suggested that all three quantification methods overestimated library concentration, but the PowerSeq kit was most accurate. Samples prepared with the TruSeq library kit provided the highest coverage and the fewest instances of dropout and below-threshold alleles compared with the KAPA kit. Additionally, all bone and hair samples demonstrated full profile completeness, with bone samples yielding a higher average coverage than hair samples. Overall, our study demonstrated that the 46GY manufacturer's protocol produced the best quality results compared to alternative library preparation options.

Recent advances in Forensic DNA Phenotyping of appearance, ancestry and age

Forensic DNA Phenotyping (FDP) comprises the prediction of a person's externally visible characteristics regarding appearance, biogeographic ancestry and age from DNA of crime scene samples, to provide investigative leads to help find unknown perpetrators that cannot be identified with forensic STR-profiling. In recent years, FDP has advanced considerably in all of its three components, which we summarize in this review article. Appearance prediction from DNA has broadened beyond eye, hair and skin color to additionally comprise other traits such as eyebrow color, freckles, hair structure, hair loss in men, and tall stature. Biogeographic ancestry inference from DNA has progressed from continental ancestry to sub-continental ancestry detection and the resolving of co-ancestry patterns in genetically admixed individuals. Age estimation from DNA has widened beyond blood to more somatic tissues such as saliva and bones as well as new markers and tools for semen. Technological progress has allowed forensically suitable DNA technology with largely increased multiplex capacity for the simultaneous analysis of hundreds of DNA predictors with targeted massively parallel sequencing (MPS). Forensically validated MPS-based FDP tools for predicting from crime scene DNA i) several appearance traits, ii) multi-regional ancestry, iii) several appearance traits together with multi-regional ancestry, and iv) age from different tissue types, are already available. Despite recent advances that will likely increase the impact of FDP in criminal casework in the near future, moving reliable appearance, ancestry and age prediction from crime scene DNA to the level of detail and accuracy police investigators may desire, requires further intensified scientific research together with technical developments and forensic validations as well as the necessary funding.

Multiple founding paternal lineages inferred from the newly-developed 639-plex Y-SNP panel suggested the complex admixture and migration history of Chinese people

BACKGROUND

Non-recombining regions of the Y-chromosome recorded the evolutionary traces of male human populations and are inherited haplotype-dependently and male-specifically. Recent whole Y-chromosome sequencing studies have identified previously unrecognized population divergence, expansion and admixture processes, which promotes a better understanding and application of the observed patterns of Y-chromosome genetic diversity.

RESULTS

Here, we developed one highest-resolution Y-chromosome single nucleotide polymorphism (Y-SNP) panel targeted for uniparental genealogy reconstruction and paternal biogeographical ancestry inference, which included 639 phylogenetically informative SNPs. We genotyped these loci in 1033 Chinese male individuals from 33 ethnolinguistically diverse populations and identified 256 terminal Y-chromosomal lineages with frequency ranging from 0.0010 (singleton) to 0.0687. We identified six dominant common founding lineages associated with different ethnolinguistic backgrounds, which included O2a2b1a1a1a1a1a1a1-M6539, O2a1b1a1a1a1a1a1-F17, O2a2b1a1a1a1a1b1a1b-MF15397, O2a2b2a1b1-A16609, O1b1a1a1a1b2a1a1-F2517, and O2a2b1a1a1a1a1a1-F155. The AMOVA and nucleotide diversity estimates revealed considerable differences and high genetic diversity among ethnolinguistically different populations. We constructed one representative phylogenetic tree among 33 studied populations based on the haplogroup frequency spectrum and sequence variations. Clustering patterns in principal component analysis and multidimensional scaling results showed a genetic differentiation between Tai-Kadai-speaking Li, Mongolic-speaking Mongolian, and other Sinitic-speaking Han Chinese populations. Phylogenetic topology inferred from the BEAST and Network relationships reconstructed from the popART further showed the founding lineages from culturally/linguistically diverse populations, such as C2a/C2b was dominant in Mongolian people and O1a/O1b was dominant in island Li people. We also identified many lineages shared by more than two ethnolinguistically different populations with a high proportion, suggesting their extensive admixture and migration history.

CONCLUSIONS

Our findings indicated that our developed high-resolution Y-SNP panel included major dominant Y-lineages of Chinese populations from different ethnic groups and geographical regions, which can be used as the primary and powerful tool for forensic practice. We should emphasize the necessity and importance of whole sequencing of more ethnolinguistically different populations, which can help identify more unrecognized population-specific variations for the promotion of Y-chromosome-based forensic applications.

Development and evaluations of the ancestry informative markers of the VISAGE Enhanced Tool for Appearance and Ancestry

The VISAGE Enhanced Tool for Appearance and Ancestry (ET) has been designed to combine markers for the prediction of bio-geographical ancestry plus a range of externally visible characteristics into a single massively parallel sequencing (MPS) assay. We describe the development of the ancestry panel markers used in ET, and the enhanced analyses they provide compared to previous MPS-based forensic ancestry assays. As well as established autosomal single nucleotide polymorphisms (SNPs) that differentiate sub-Saharan African, European, East Asian, South Asian, Native American, and Oceanian populations, ET includes autosomal SNPs able to efficiently differentiate populations from Middle East regions. The ability of the ET autosomal ancestry SNPs to distinguish Middle East populations from other continentally defined population groups is such that characteristic patterns for this region can be discerned in genetic cluster analysis using STRUCTURE. Joint cluster membership estimates showing individual co-ancestry that signals North African or East African origins were detected, or cluster patterns were seen that indicate origins from central and Eastern regions of the Middle East. In addition to an augmented panel of autosomal SNPs, ET includes panels of 85 Y-SNPs, 16 X-SNPs and 21 autosomal Microhaplotypes. The Y- and X-SNPs provide a distinct method for obtaining extra detail about co-ancestry patterns identified in males with admixed backgrounds. This study used the 1000 Genomes admixed African and admixed American sample sets to fully explore these enhancements to the analysis of individual co-ancestry. Samples from urban and rural Brazil with contrasting distributions of African, European, and Native American co-ancestry were also studied to gauge the efficiency of combining Y- and X-SNP data for this purpose. The small panel of Microhaplotypes incorporated in ET were selected because they showed the highest levels of haplotype diversity amongst the seven population groups we sought to differentiate. Microhaplotype data was not formally combined with single-site SNP genotypes to analyse ancestry. However, the haplotype sequence reads obtained with ET from these loci creates an effective system for de-convoluting two-contributor mixed DNA. We made simple mixture experiments to demonstrate that when the contributors have different ancestries and the mixture ratios are imbalanced (i.e., not 1:1 mixtures) the ET Microhaplotype panel is an informative system to infer ancestry when this differs between the contributors.

Comparative evaluation of the MAPlex, Precision ID Ancestry Panel, and VISAGE Basic Tool for biogeographical ancestry inference

Biogeographical ancestry (BGA) inference from ancestry-informative markers (AIMs) has strong potential to support forensic investigations. Over the past two decades, several forensic panels composed of AIMs have been developed to predict ancestry at a continental scale. These panels typically comprise fewer than 200 AIMs and have been designed and tested with a limited set of populations. How well these panels recover patterns of genetic diversity relative to larger sets of markers, and how accurately they infer ancestry of individuals and populations not included in their design remains poorly understood. The lack of comparative studies addressing these aspects makes the selection of appropriate panels for forensic laboratories difficult. In this study, the model-based genetic clustering tool STRUCTURE was used to compare three popular forensic BGA panels: MAPlex, Precision ID Ancestry Panel (PIDAP), and VISAGE Basic Tool (VISAGE BT) relative to a genome-wide reference set of 10k SNPs. The genotypes for all these markers were obtained for a comprehensive set of 3957 individuals from 228 worldwide human populations. Our results indicate that at the broad continental scale (K=6) typically examined in forensic studies, all forensic panels produced similar genetic structure patterns compared to the reference set (G'≈90%) and had high classification performance across all regions (average AUC-PR > 97%). However, at K= 7 and K= 8, the forensic panels displayed some region-specific clustering deviations from the reference set, particularly in Europe and the region of East and South-East Asia, which may be attributed to differences in the design of the respective panels. Overall, the panel with the most consistent performance in all regions was VISAGE BT with an average weighted AUC̅_W score of 96.26% across the three scales of geographical resolution investigated.

Developmental validation of a 381 Y-chromosome SNP panel for haplogroup analysis in the Chinese populations

Y-chromosome single nucleotide polymorphism (Y-SNP) shows great variation in geographical distribution and population heterogeneity and can be used to map population genetics around the world. Massive parallel sequencing (MPS) methodology enables high-resolution Y-SNP haplogrouping for a certain male and is widely used in forensic genetics and evolutionary studies. In this present study, we used MPS to develop a customized 381 Y-SNP panel (SifaMPS 381 Y-SNP panel) to investigate the basic structure and subbranches of the haplogroup tree of the Chinese populations. The SifaMPS 381 Y-SNP panel covers all the Y-SNPs from our previously designed 183 Y-SNP panel and additional SNPs under the predominant haplogroups in the Chinese populations based on certain criteria. We also evaluated the sequencing matrix, concordance, sensitivity, repeatability of this panel and the ability to analyze mixed and case-type samples based on the Illumina MiSeq System. The results demonstrated that the novel MPS Y-SNP panel possessed good sequencing performance and generated accurate Y-SNP genotyping results. Although the recommended DNA input was greater than 1.25 ng, we observed that a lower DNA amount could still be used to analyze haplogroups correctly. In addition, this panel could handle mixed samples and common case-type samples and had higher resolution among Chinese Han males than previously reported. In conclusion, the SifaMPS 381 Y-SNP panel showed an overall good performance and offers a better choice for Y-SNP haplogrouping of the Chinese population, thereby facilitating paternal lineage classification, familial searching and other forensic applications.

Multiplex Snapshot minisequencing for the detection of common PAH gene mutations in Iranian patients with Phenylketonuria

Background

Phenylketonuria is a common inborn defect of amino acid metabolism in the world. This failure is caused by an autosomal recessive insufficiency of the hepatic enzyme hyperphenylalaninemia (PAH), which catalyzes the irreversible hydroxylation of phenylalanine to tyrosine. More than 1,040 different disease-causing mutations have already been identified in the PAH gene. The most prominent complication of Phenylketonuria, if not diagnosed and treated, is severe mental retardation. Hence, early diagnosis and initiation of nutritional therapy are the most significant measures in preventing this mental disorder. Given these data, we developed a simple and rapid molecular test to detect the most frequent PAH mutations.

Methods

Multiplex assay was developed based on the SNaPshot minisequencing approach to simultaneously perform genotyping of the 10 mutations at the PAH gene. We optimized detection of these mutations in one multiplex PCR, followed by 10 single-nucleotide extension reactions. DNA sequencing assay was also used to verify genotyping results obtained by SNaPshot minisequencing.

Result

All 10 genotypes were determined based on the position and the fluorescent color of the peaks in a single electropherogram. Sequencing results of these frequent mutations showed that by using this method, a 100% detection rate could be achieved in the Iranian population.

Conclusion

SNaPshot minisequencing can be useful as a secondary test in neonatal screening for HPA in neonates with a positive screening test, and it is also suitable for carrier screening. The assay can be easily applied for accurate and time- and cost-efficient genotyping of the selected SNPs in various population.

Recent advances in forensic biology and forensic DNA typing: INTERPOL review 2019-2022

This review paper covers the forensic-relevant literature in biological sciences from 2019 to 2022 as a part of the 20th INTERPOL International Forensic Science Managers Symposium. Topics reviewed include rapid DNA testing, using law enforcement DNA databases plus investigative genetic genealogy DNA databases along with privacy/ethical issues, forensic biology and body fluid identification, DNA extraction and typing methods, mixture interpretation involving probabilistic genotyping software (PGS), DNA transfer and activity-level evaluations, next-generation sequencing (NGS), DNA phenotyping, lineage markers (Y-chromosome, mitochondrial DNA, X-chromosome), new markers and approaches (microhaplotypes, proteomics, and microbial DNA), kinship analysis and human identification with disaster victim identification (DVI), and non-human DNA testing including wildlife forensics. Available books and review articles are summarized as well as 70 guidance documents to assist in quality control that were published in the past three years by various groups within the United States and around the world.

Applications of massively parallel sequencing in forensic genetics

Massively parallel sequencing, also referred to as next-generation sequencing, has positively changed DNA analysis, allowing further advances in genetics. Its capability of dealing with low quantity/damaged samples makes it an interesting instrument for forensics. The main advantage of MPS is the possibility of analyzing simultaneously thousands of genetic markers, generating high-resolution data. Its detailed sequence information allowed the discovery of variations in core forensic short tandem repeat loci, as well as the identification of previous unknown polymorphisms. Furthermore, different types of markers can be sequenced in a single run, enabling the emergence of DIP-STRs, SNP-STR haplotypes, and microhaplotypes, which can be very useful in mixture deconvolution cases. In addition, the multiplex analysis of different single nucleotide polymorphisms can provide valuable information about identity, biogeographic ancestry, paternity, or phenotype. DNA methylation patterns, mitochondrial DNA, mRNA, and microRNA profiling can also be analyzed for different purposes, such as age inference, maternal lineage analysis, body-fluid identification, and monozygotic twin discrimination. MPS technology also empowers the study of metagenomics, which analyzes genetic material from a microbial community to obtain information about individual identification, post-mortem interval estimation, geolocation inference, and substrate analysis. This review aims to discuss the main applications of MPS in forensic genetics.

HIV-1 Drug Resistance Assay Using Ion Torrent Next Generation Sequencing and On-Instrument End-to-End Analysis Software

HIV-1 antiretroviral therapy management requires sequencing the protease, reverse transcriptase, and integrase portions of the HIV-1 pol gene. Most resistance testing is performed with Sanger sequencing, which has limited ability to detect minor variants. Next generation sequencing (NGS) platforms enable variant detection at frequencies as low as 1% allowing for earlier detection of resistance and modification of therapy. Implementation of NGS assays in the clinical laboratory is hindered by complicated assay design, cumbersome wet bench procedures, and the complexity of data analysis and bioinformatics. We developed a complete NGS protocol and companion analysis and reporting pipeline using AmpliSeq multiplex PCR, Ion Torrent S5 XL sequencing, and Stanford's HIVdb resistance algorithm. Implemented as a Torrent Suite software plugin, the pipeline runs automatically after sequencing. An optimum variant frequency threshold of 10% was determined by comparing Sanger sequences of archived samples from ViroSeq testing, resulting in a sensitivity of 98.2% and specificity of 99.0%. The majority (91%) of drug resistance mutations were detected by both Sanger and NGS, with 1.7% only by Sanger and 7.3% only by NGS. Variant calls were highly reproducible and there was no cross-reactivity to VZV, HBV, CMV, EBV, and HCV. The limit of detection was 500 copies/mL. The NGS assay performance was comparable to ViroSeq Sanger sequencing and has several advantages, including a publicly available end-to-end analysis and reporting plugin. The assay provides a straightforward path for implementation of NGS for HIV drug resistance testing in the laboratory setting without additional investment in bioinformatics infrastructure and resources.

On the Forensic Use of Y-Chromosome Polymorphisms

Nowadays, the use of Y-chromosome polymorphisms forms an essential part of many forensic DNA investigations. However, this was not always the case. Only since 1992 have we seen that some forensic scientists started to have an interest in this chromosome. In this review, I will sketch a brief history focusing on the forensic use of Y-chromosome polymorphisms. Before describing the various applications of short-tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) on the Y-chromosome, I will discuss a few often ignored aspects influencing proper use and interpretation of Y-chromosome information: (i) genotyping Y-SNPs and Y-STRs, (ii) Y-STR haplotypes shared identical by state (IBS) or identical by descent (IBD), and (iii) Y-haplotype database frequencies.

A multi-dimensional evaluation of the 'NIST 1032' sample set across four forensic Y-STR multiplexes

This manuscript reports Y-chromosomal short tandem repeat (Y-STR) haplotypes for 1032 male U.S. population samples across 30 Y-STR loci characterized by three capillary electrophoresis (CE) length-based kits (PowerPlex Y23 System, Yfiler Plus PCR Amplification Kit, and Investigator Argus Y-28 QS Kit) and one sequence-based kit (ForenSeq DNA Signature Prep Kit): DYF387S1, DYS19, DYS385 a/b, DYS389I, DYS389II, DYS390, DYS391, DYS392, DYS393, DYS437, DYS438, DYS439, DYS448, DYS449, DYS456, DYS458, DYS460, DYS481, DYS505, DYS518, DYS522, DYS533, DYS549, DYS570, DYS576, DYS612, DYS627, DYS635, DYS643, and Y-GATA-H4. The length-based Y-STR haplotypes include six loci that are not reported in the sequence-based kit (DYS393, DYS449, DYS456, DYS458, DYS518, and DYS627), whereas three loci included in the sequence-based kit are not present in length-based kits (DYS505, DYS522, and DYS612). For the latter, a custom multiplex was used to generate CE length-based data, allowing 1032 samples to be evaluated for concordance across the 30 Y-STR loci included in these four commercial Y-STR typing kits. Discordances between typing methods were analyzed further to assess underlying causes such as primer binding site mutations and flanking region insertions/deletions. Allele-level frequency and statistical information is provided for sequenced loci, excluding the multi-copy loci DYF387S1 and DYS385 a/b, for which locus-specific haplotype-level frequencies are provided instead. The resulting data reveals the degree of information gained through sequencing: 88% of sequenced Y-STR loci contain additional sequence-based alleles compared to length-based data, with the DYS389II locus containing the most additional alleles (51) observed by sequencing. Despite these allelic increases, only minimal improvement was observed in haplotype resolution by sequence, with all four commercial kits providing a similar ability to differentiate length-based haplotypes in this sample set. Finally, a subset of 369 male samples were compared to their corresponding additionally sequenced father samples, revealing the sequence basis for the 50 length-based changes observed, and no additional sequence-based mutations. GenBank accession numbers are reported for each unique sequence, and associated records are available in the STRSeq Y-Chromosomal STR Loci National Center for Biotechnology Information (NCBI) BioProject, accession PRJNA380347. Haplotype data is updated in the Y-STR Haplotype Reference Database (YHRD) for the 'NIST 1032' data set to now achieve the level of maximal haplotype of YHRD. All supplementary files including revisions to previously published Y-STR data are available in the NIST Public Data Repository: U.S. population data for human identification markers, DOI 10.18434/t4/1500024.

Improving the regional Y-STR haplotype resolution utilizing haplogroup-determining Y-SNPs and the application of machine learning in Y-SNP haplogroup prediction in a forensic Y-STR database: A pilot study on male Chinese Yunnan Zhaoyang Han population

Improving the resolution of the current widely used Y-chromosomal short tandem repeat (Y-STR) dataset is of great importance for forensic investigators, and the current approach is limited, except for the addition of more Y-STR loci. In this research, a regional Y-DNA database was investigated to improve the Y-STR haplotype resolution utilizing a Y-SNP Pedigree Tagging System that includes 24 Y-chromosomal single nucleotide polymorphism (Y-SNP) loci. This pilot study was conducted in the Chinese Yunnan Zhaoyang Han population, and 3473 unrelated male individuals were enrolled. Based on data on the male haplogroups under different panels, the matched or near-matching (NM) Y-STR haplotype pairs from different haplogroups indicated the critical roles of haplogroups in improving the regional Y-STR haplotype resolution. A classic median-joining network analysis was performed using Y-STR or Y-STR/Y-SNP data to reconstruct population substructures, which revealed the ability of Y-SNPs to correct misclassifications from Y-STRs. Additionally, population substructures were reconstructed using multiple unsupervised or supervised dimensionality reduction methods, which indicated the potential of Y-STR haplotypes in predicting Y-SNP haplogroups. Haplogroup prediction models were built based on nine publicly accessible machine-learning (ML) approaches. The results showed that the best prediction accuracy score could reach 99.71% for major haplogroups and 98.54% for detailed haplogroups. Potential influences on prediction accuracy were assessed by adjusting the Y-STR locus numbers, selecting Y-STR loci with various mutabilities, and performing data processing. ML-based predictors generally presented a better prediction accuracy than two available predictors (Nevgen and EA-YPredictor). Three tree models were developed based on the Yfiler Plus panel with unprocessed input data, which showed their strong generalization ability in classifying various Chinese Han subgroups (validation dataset). In conclusion, this study revealed the significance and application prospects of Y-SNP haplogroups in improving regional Y-STR databases. Y-SNP haplogroups can be used to discriminate NM Y-STR haplotype pairs, and it is important for forensic Y-STR databases to develop haplogroup prediction tools to improve the accuracy of biogeographic ancestry inferences.

The FORCE Panel: An All-in-One SNP Marker Set for Confirming Investigative Genetic Genealogy Leads and for General Forensic Applications

The FORensic Capture Enrichment (FORCE) panel is an all-in-one SNP panel for forensic applications. This panel of 5422 markers encompasses common, forensically relevant SNPs (identity, ancestry, phenotype, X- and Y-chromosomal SNPs), a novel set of 3931 autosomal SNPs for extended kinship analysis, and no clinically relevant/disease markers. The FORCE panel was developed as a custom hybridization capture assay utilizing ~20,000 baits to target the selected SNPs. Five non-probative, previously identified World War II (WWII) cases were used to assess the kinship panel. Each case included one bone sample and associated family reference DNA samples. Additionally, seven reference quality samples, two 200-year-old bone samples, and four control DNAs were processed for kit performance and concordance assessments. SNP recovery after capture resulted in a mean of ~99% SNPs exceeding 10X coverage for reference and control samples, and 44.4% SNPs for bone samples. The WWII case results showed that the FORCE panel could predict first to fifth degree relationships with strong statistical support (likelihood ratios over 10,000 and posterior probabilities over 99.99%). To conclude, SNPs will be important for further advances in forensic DNA analysis. The FORCE panel shows promising results and demonstrates the utility of a 5000 SNP panel for forensic applications.

An Introductory Overview of Open-Source and Commercial Software Options for the Analysis of Forensic Sequencing Data

The top challenges of adopting new methods to forensic DNA analysis in routine laboratories are often the capital investment and the expertise required to implement and validate such methods locally. In the case of next-generation sequencing, in the last decade, several specifically forensic commercial options became available, offering reliable and validated solutions. Despite this, the readily available expertise to analyze, interpret and understand such data is still perceived to be lagging behind. This review gives an introductory overview for the forensic scientists who are at the beginning of their journey with implementing next-generation sequencing locally and because most in the field do not have a bioinformatics background may find it difficult to navigate the new terms and analysis options available. The currently available open-source and commercial software for forensic sequencing data analysis are summarized here to provide an accessible starting point for those fairly new to the forensic application of massively parallel sequencing.

CSYseq: The first Y-chromosome sequencing tool typing a large number of Y-SNPs and Y-STRs to unravel worldwide human population genetics

Male-specific Y-chromosome (chrY) polymorphisms are interesting components of the DNA for population genetics. While single nucleotide polymorphisms (Y-SNPs) indicate distant evolutionary ancestry, short tandem repeats (Y-STRs) are able to identify close familial kinships. Detailed chrY analysis provides thus both biogeographical background information as paternal lineage identification. The rapid advancement of high-throughput massive parallel sequencing (MPS) technology in the past decade has revolutionized genetic research. Using MPS, single-base information of both Y-SNPs as Y-STRs can be analyzed in a single assay typing multiple samples at once. In this study, we present the first extensive chrY-specific targeted resequencing panel, the 'CSYseq', which simultaneously identifies slow mutating Y-SNPs as evolution markers and rapid mutating Y-STRs as patrilineage markers. The panel was validated by paired-end sequencing of 130 males, distributed over 65 deep-rooted pedigrees covering 1,279 generations. The CSYseq successfully targets 15,611 Y-SNPs including 9,014 phylogenetic informative Y-SNPs to identify 1,443 human evolutionary Y-subhaplogroup lineages worldwide. In addition, the CSYseq properly targets 202 Y-STRs, including 81 slow, 68 moderate, 27 fast and 26 rapid mutating Y-STRs to individualize close paternal relatives. The targeted chrY markers cover a high average number of reads (Y-SNP = 717, Y-STR = 150), easy interpretation, powerful discrimination capacity and chrY specificity. The CSYseq is interesting for research on different time scales: to identify evolutionary ancestry, to find distant family and to discriminate closely related males. Therefore, this panel serves as a unique tool valuable for a wide range of genetic-genealogical applications in interdisciplinary research within evolutionary, population, molecular, medical and forensic genetics.

The analysis of ancestry with small-scale forensic panels of genetic markers

In the last 10 years, forensic genetic analysis has been extended beyond identification tests that link a suspect to crime scene evidence using standard DNA profiling, to new supplementary tests that can provide information to investigators about a suspect in the absence of a database hit or eyewitness testimony. These tests now encompass the prediction of physical appearance, ancestry and age. In this review, we give a comprehensive overview of the full range of DNA-based ancestry inference tests designed to work with forensic contact traces, when the level of DNA is often very low or highly degraded. We outline recent developments in the design of ancestry-informative marker sets, forensic assays that use capillary electrophoresis or massively parallel sequencing, and the statistical analysis frameworks that examine the test profile and compares it to reference population variation. Three casework ancestry analysis examples are described which were successfully accomplished in the authors' laboratory, where the ancestry information obtained was critical to the outcome of the DNA analyses made.

The impact of correlations between pigmentation phenotypes and underlying genotypes on genetic prediction of pigmentation traits

Predicting appearance phenotypes from genotypes is relevant for various areas of human genetic research and applications such as genetic epidemiology, human history, anthropology, and particularly in forensics. Many appearance phenotypes, and thus their underlying genotypes, are highly correlated, with pigmentation traits serving as primary examples. However, all available genetic prediction models, including those for pigmentation traits currently used in forensic DNA phenotyping, ignore phenotype correlations. Here, we investigated the impact of appearance phenotype correlations on genetic appearance prediction in the exemplary case of three pigmentation traits. We used data for categorical eye, hair and skin colour as well as 41 DNA markers utilized in the recently established HIrisPlex-S system from 762 individuals with complete phenotype and genotype information. Based on these data, we performed genetic prediction modelling of eye, hair and skin colour via three different strategies, namely the established approach of predicting phenotypes solely based on genotypes while not considering phenotype correlations, and two novel approaches that considered phenotype correlations, either incorporating truly observed correlated phenotypes or DNA-predicted correlated phenotypes in addition to the DNA predictors. We found that using truly observed correlated pigmentation phenotypes as additional predictors increased the DNA-based prediction accuracies for almost all eye, hair and skin colour categories, with the largest increase for intermediate eye colour, brown hair colour, dark to black skin colour, and particularly for dark skin colour. Outcomes of dedicated computer simulations suggest that this prediction accuracy increase is due to the additional genetic information that is implicitly provided by the truly observed correlated pigmentation phenotypes used, yet not covered by the DNA predictors applied. In contrast, considering DNA-predicted correlated pigmentation phenotypes as additional predictors did not improve the performance of the genetic prediction of eye, hair and skin colour, which was in line with the results from our computer simulations. Hence, in practical applications of DNA-based appearance prediction where no phenotype knowledge is available, such as in forensic DNA phenotyping, it is not advised to use DNA-predicted correlated phenotypes as predictors in addition to the DNA predictors. In the very least, this is not recommended for the pigmentation traits and the established pigmentation DNA predictors tested here.

Development and validation of a custom panel including 183 Y-SNPs for Chinese Y-chromosomal haplogroups dissection using a MALDI-TOF MS system

Y-chromosome SNP haplogroups exhibit geographic structuring in many populations around the world. Therefore, Y-chromosome haplogroups have been widely used to infer paternal biogeographical ancestry and high-resolution paternal lineage classification. In the present study, we designed a customized panel containing 183 Y-SNPs based on previous studies and evaluated the genotyping performance and repeatability, concordance, sensitivity, and ability of analysing case-type samples using a MALDI-TOF MS platform. The average call rate for duplicate typing of any one SNP in the panel was 97.0%. In the concordance and accuracy study, the results of haplogroup designation obtained from MALDI-TOF MS platform were fully consistent with those obtained from the next-generation sequencing (NGS) platform. The optimal amount of template DNA in the PCR seemed to be 10 ng. However, if less DNA (≥156.25 pg) was available, it was still possible to obtain meaningful haplogroup information. For the application part, this panel could be applied for the detection of blood, semen, and buccal swabs samples. Particularly, blood stain on FTA card samples could be dissected by direct PCR amplification on the MALDI-TOF MS platform. Besides, 371 unrelated male individuals from four Chinese ethnic groups (Han, Hui, Mongolian, and Kazak) were detected using this panel. Total 78 terminal haplogroups were found and the haplogroup diversity was 0.933576. The results demonstrate that this panel could be an accurate, fast, and cost-effective method for database construction where the amount of sample material is less of a concern and when the cost of the assay is taken into consideration.

Predicting haplogroups using a versatile machine learning program (PredYMaLe) on a new mutationally balanced 32 Y-STR multiplex (CombYplex): Unlocking the full potential of the human STR mutation rate spectrum to estimate forensic parameters

We developed a new mutationally well-balanced 32 Y-STR multiplex (CombYplex) together with a machine learning (ML) program PredYMaLe to assess the impact of STR mutability on haplogourp prediction, while respecting forensic community criteria (high DC/HD). We designed CombYplex around two sub-panels M1 and M2 characterized by average and high-mutation STR panels. Using these two sub-panels, we tested how our program PredYmale reacts to mutability when considering basal branches and, moving down, terminal branches. We tested first the discrimination capacity of CombYplex on 996 human samples using various forensic and statistical parameters and showed that its resolution is sufficient to separate haplogroup classes. In parallel, PredYMaLe was designed and used to test whether a ML approach can predict haplogroup classes from Y-STR profiles. Applied to our kit, SVM and Random Forest classifiers perform very well (average 97 %), better than Neural Network (average 91 %) and Bayesian methods (< 90 %). We observe heterogeneity in haplogroup assignation accuracy among classes, with most haplogroups having high prediction scores (99-100 %) and two (E1b1b and G) having lower scores (67 %). The small sample sizes of these classes explain the high tendency to misclassify the Y-profiles of these haplogroups; results were measurably improved as soon as more training data were added. We provide evidence that our ML approach is a robust method to accurately predict haplogroups when it is combined with a sufficient number of markers, well-balanced mutation rate Y-STR panels, and large ML training sets. Further research on confounding factors (such as CNV-STR or gene conversion) and ideal STR panels in regard to the branches analysed can be developed to help classifiers further optimize prediction scores.

The construction and application of a new 17-plex Y-STR system using universal fluorescent PCR

Y-chromosomal short tandem repeat (Y-STR) polymorphisms are useful in forensic identification, population genetics, and human structures. However, the current Y-STR systems are limited in discriminating distant relatives in a family with a low discrimination power. Increasing the capacity of detecting Y chromosomal polymorphisms will drastically narrow down the matching number of genealogy populations or pedigrees. In this study, we developed a system containing 17 Y-STRs that are complementary to the current commercially available Y-STR kits. This system was constructed by multiplex PCR with expected sizes of 126-400 bp labeled by different fluorescence molecules (DYS715, DYS709, DYS716, DYS713, and DYS607 labeled by FAM; DYS718, DYS723, DYS708, and DYS714 labeled by JOE; DYS712, DYS717, DYS721, and DYS605 labeled by TAMRA; and DYS719, DYS726, DYS598, and DYS722 labeled by ROX). The system was extensively tested for sensitivity, male specificity, species specificity, mixture, population genetics, and mutation rates following the Scientific Working Group on DNA Analysis Methods (SWGDAM) guidelines. The genetic data were obtained from eight populations with a total of 1260 individuals. Our results showed that all the 17 Y-STRs are human- and male-specific and include only one copy of the Y-chromosome. The 17 Y-STR system detects 143 alleles and has a high discrimination power (0.996031746). Mutation rates were different among the 17 Y-STRs, ranging from 0.30 to 3.03%. In conclusion, our study provides a robust, sensitive, and cost-effective genotyping method for human identification, which will be beneficial for narrowing the search scope when applied to genealogy searching with the Y-STR DNA databank.

The MASTiFF panel-a versatile multiple-allele SNP test for forensics

Forensic identification tests often need recourse to markers that can successfully type highly degraded DNA, and binary single nucleotide polymorphisms (SNPs) have become the variants of choice for such analyses because of their short amplified fragment lengths. The two main drawbacks of SNPs are their reduced power of discrimination per marker compared with mainstream forensic STRs and an inability to robustly detect mixed DNA-particularly using capillary electrophoresis genotyping systems such as SNaPshot™, where the dye signals are much more imbalanced than those of STR profiles. This study compiled a compact set of multiple-allele SNPs consisting of loci that had three or four nucleotide variants at the same site in order to address the lack of mixture detection capability with binary SNP tests, as well as improving levels of polymorphism per SNP by transitioning to a maximum of six or ten genotypes per locus. We report the development and optimisation of a SNaPshot-based forensic test comprising 27 tri-allelic and 2 tetra-allelic SNPs, which we named MASTiFF: a multiple-allele SNP test for forensics. Assessments of the MASTiFF panel's levels of discrimination power in the five main population groups indicate random match probabilities ranging from 10^-15 down to 10^-20-improving the levels possible from an equivalent number of binary SNPs. The SNaPshot test was able to detect simple mixtures successfully with more than two alleles observed in 30% of SNPs. From allele frequency data, it is estimated that more than two alleles will be present in at least one MASTiFF SNP in 99.8% of two-person mixtures, making this panel an ideal supplementary test when SNPs are chosen for the analysis of degraded forensic DNA.

Update on the predictability of tall stature from DNA markers in Europeans

Predicting adult height from DNA has important implications in forensic DNA phenotyping. In 2014, we introduced a prediction model consisting of 180 height-associated SNPs based on data from 10,361 Northwestern Europeans enriched with tall individuals (770 > 1.88 standard deviation), which yielded a mid-ranged accuracy (AUC = 0.75 for binary prediction of tall stature and R² = 0.12 for quantitative prediction of adult height). Here, we provide an update on DNA-based height predictability considering an enlarged list of subsequently-published height-associated SNPs using data from the same set of 10,361 Europeans. A prediction model based on the full set of 689 SNPs showed an improved accuracy relative to previous models for both tall stature (AUC = 0.79) and quantitative height (R² = 0.21). A feature selection analysis revealed a subset of 412 most informative SNPs while the corresponding prediction model retained most of the accuracy (AUC = 0.76 and R² = 0.19) achieved with the full model. Over all, our study empirically exemplifies that the accuracy for predicting human appearance phenotypes with very complex underlying genetic architectures, such as adult height, can be improved by increasing the number of phenotype-associated DNA variants. Our work also demonstrates that a careful sub-selection allows for a considerable reduction of the number of DNA predictors that achieve similar prediction accuracy as provided by the full set. This is forensically relevant due to restrictions in the number of SNPs simultaneously analyzable with forensically suitable DNA technologies in the current days of targeted massively parallel sequencing in forensic genetics.