SNPs for a universal individual identification panel

Malignancy following solid organ transplantation: Current techniques for determination of donor versus recipient origin

Among the post-transplantation complications that patients may encounter, the transmission of a donor-derived malignant neoplasm is uncommon but potentially life threatening. The determination of donor versus recipient origin is essential particularly in the setting of multiple transplant recipients from the donor. Advances in molecular biology now allow accurate discrimination utilizing routine tissue samples in a timely and cost-effective manner. The techniques are routinely performed in hospital molecular biology laboratories and are also available in commercial labs. The current methodologies are discussed and future possibilities are presented for clinicians caring for solid organ recipients.

Identification of individuals from low template blood samples using whole transcriptome shotgun sequencing

Biological trace samples consisting of very few cells pose a challenge to conventional forensic genetic DNA analysis. RNA may be an alternative to DNA when handling low template samples. Whereas each cell only contains two copies of an autosomal DNA segment, the transcriptome retains much of the genomic variation replicated in abundant RNA fragments. In this study, we describe the development of a prototype RNA-based SNP selection set for forensic human identification from low template samples (50 pg gDNA). Whole blood from a subset of the Danish population (41 individuals) and blood stains subjected to degradation at room temperature for up to two weeks were analysed by whole transcriptome shotgun sequencing. Concordance was determined by DNA genotyping with the Infinium Omni5-4 SNP chip. In the 100 protein-coding genes with the most reads, 5214 bi-allelic SNPs with gnomAD minor allele frequencies > 0.1 in the African/African American, East Asian, and (non-Finnish) European populations were identified. Of these, 24 SNPs in 21 genes passed screening in whole blood and degraded blood stains, with a resulting mean match probability of 4.5 ∙ 10-9. Additionally, ancestry informative SNPs and SNPs in genes useful for body fluid identification were identified in the transcriptome. Consequently, shotgun sequencing of RNA from low template samples may be used for a vast host of forensic genetics purposes, including simultaneous human and body fluid identification, leading to direct donor identification in the identified body fluid.

Single nucleotide polymorphism-based biomarker in primary hypertension

Primary hypertension is a multiplex and multifactorial disease influenced by various strong components including genetics. Extensive research such as Genome-wide association studies and candidate gene studies have revealed various single nucleotide polymorphisms (SNPs) related to hypertension, providing insights into the genetic basis of the condition. This review summarizes the current status of SNP research in primary hypertension, including examples of hypertension-related SNPs, their location, function, and frequency in different populations. The potential clinical implications of SNP research for primary hypertension management are also discussed, including disease risk prediction, personalized medicine, mechanistic understanding, and lifestyle modifications. Furthermore, this review highlights emerging technologies and methodologies that have the potential to revolutionize the vast understanding of the basis of genetics in primary hypertension. Gene editing holds the potential to target and correct any kind of genetic mutations that contribute to the development of hypertension or modify genes involved in blood pressure regulation to prevent or treat the condition. Advances in computational biology and machine learning enable researchers to analyze large datasets and identify complex genetic interactions contributing to hypertension risk. In conclusion, SNP research in primary hypertension is rapidly evolving with emerging technologies and methodologies that have the potential to transform the knowledge about genetic basis related to the condition. These advances hold promise for personalized prevention and treatment strategies tailored to an individual's genetic profile ultimately improving patient outcomes and reducing healthcare costs.

SWGTS-a platform for stream-based host DNA depletion

MOTIVATION

Microbial sequencing data from clinical samples is often contaminated with human sequences, which have to be removed prior to sharing. Existing methods for human read removal, however, are applicable only after the target dataset has been retrieved in its entirety, putting the recipient at least temporarily in control of a potentially identifiable genetic dataset with potential implications under regulatory frameworks such as the GDPR. In some instances, the ability to carry out stream-based host depletion as part of the data transfer process may be preferable.

RESULTS

We present SWGTS, a client-server application for the transfer and stream-based host depletion of sequencing reads. SWGTS enforces a robust upper bound on the maximum amount of human genetic data from any one client held in memory at any point in time by storing all incoming sequencing data in a limited-size, client-specific intermediate processing buffer, and by throttling the rate of incoming data if it exceeds the speed of host depletion carried out on the SWGTS server in the background. SWGTS exposes a HTTP-REST interface, is implemented using docker-compose, Redis and traefik, and requires less than 8 Gb of RAM for deployment. We demonstrate high filtering accuracy of SWGTS; incoming data transfer rates of up to 1.65 megabases per second in a conservative configuration; and mitigation of re-identification risks by the ability to limit the number of SNPs present on a popular population-scale genotyping array covered by reads in the SWGTS buffer to a low user-defined number, such as 10 or 100.

AVAILABILITY AND IMPLEMENTATION

SWGTS is available on GitHub: https://github.com/AlBi-HHU/swgts (https://doi.org/10.5281/zenodo.10891052). The repository also contains a jupyter notebook that can be used to reproduce all the benchmarks used in this article. All datasets used for benchmarking are publicly available.

Assessing Privacy Vulnerabilities in Genetic Data Sets: Scoping Review

BACKGROUND

Genetic data are widely considered inherently identifiable. However, genetic data sets come in many shapes and sizes, and the feasibility of privacy attacks depends on their specific content. Assessing the reidentification risk of genetic data is complex, yet there is a lack of guidelines or recommendations that support data processors in performing such an evaluation.

OBJECTIVE

This study aims to gain a comprehensive understanding of the privacy vulnerabilities of genetic data and create a summary that can guide data processors in assessing the privacy risk of genetic data sets.

METHODS

We conducted a 2-step search, in which we first identified 21 reviews published between 2017 and 2023 on the topic of genomic privacy and then analyzed all references cited in the reviews (n=1645) to identify 42 unique original research studies that demonstrate a privacy attack on genetic data. We then evaluated the type and components of genetic data exploited for these attacks as well as the effort and resources needed for their implementation and their probability of success.

RESULTS

From our literature review, we derived 9 nonmutually exclusive features of genetic data that are both inherent to any genetic data set and informative about privacy risk: biological modality, experimental assay, data format or level of processing, germline versus somatic variation content, content of single nucleotide polymorphisms, short tandem repeats, aggregated sample measures, structural variants, and rare single nucleotide variants.

CONCLUSIONS

On the basis of our literature review, the evaluation of these 9 features covers the great majority of privacy-critical aspects of genetic data and thus provides a foundation and guidance for assessing genetic data risk.

Integrative lncRNA, circRNA, and mRNA analysis reveals expression profiles of six forensic body fluids/tissue

RNAs have attracted much attention in forensic body fluid/tissue identification (BFID) due to their tissue-specific expression characteristics. Among RNAs, long RNAs (e.g., mRNA) have a higher probability of containing more polymorphic sites that can be used to assign the specific donor of the body fluid/tissue. However, few studies have characterized their overall profiles in forensic science. In this study, we sequenced the transcriptomes of 30 samples from venous blood, menstrual blood, semen, saliva, vaginal secretion, and skin tissue, obtaining a comprehensive picture of mRNA, lncRNA, and circRNA profiles. A total of 90,305 mRNAs, 102,906 lncRNAs (including 19,549 novel lncRNAs), and 40,204 circRNAs were detected. RNA type distribution, length distribution, and expression distribution were presented according to their annotation and expression level, and many novel body fluid/tissue-specific RNA markers were identified. Furthermore, the cognate relations among the three RNAs were analyzed according to gene annotations. Finally, SNPs and InDels from RNA transcripts were genotyped, and 21,611 multi-SNP and 4,471 multi-InDel transcriptomic microhaplotypes (tMHs) were identified. These results provide a comprehensive understanding of transcriptome profiles, which could provide new avenues for tracing the origin of the body fluid/tissue and identifying an individual.

Engraftment and Measurable Residual Disease Monitoring after Hematopoietic Stem Cell Transplantation: Comparison of Two Chimerism Test Strategies, Next-Generation Sequencing versus a Combination of Short-Tandem Repeats and Quantitative PCR

Chimerism testing supports the study of engraftment and measurable residual disease (MRD) in patients after allogeneic hematopoietic stem cell transplant. In chimerism MRD, relapse can be predicted by increasing mixed chimerism (IMC), recipient increase ≥0.1% in peripheral blood, and proliferating recipient cells as a surrogate of tumor activity. Conventionally, the combination of short-tandem repeat (STR) and quantitative PCR (qPCR) was needed to ensure assay sensitivity and accuracy in all chimerism status. We evaluated the use of next-generation sequencing (NGS) as an alternate technique. The median numbers of informative markers in unrelated/related cases were 124/82 (NGS; from 202 single-nucleotide polymorphism), 5/3 (qPCR), and 17/10 (STR). Assay sensitivity was 0.22% (NGS), 0.1% (qPCR), and 1% (STR). NGS batch (4 to 48 samples) required 19.60 to 24.80 hours and 1.52 to 2.42 hours of hands-on time (comparable to STR/qPCR). NGS assay cost/sample was $91 to $151, similar to qPCR ($99) but higher than STR ($27). Using 56 serial DNAs from six post-transplant patients monitored by the qPCR/STR, the correlation with NGS was strong for percentage recipient (y = 1.102x + 0.010; R2 = 0.968) and percentage recipient change (y = 0.892x + 0.041; R2 = 0.945). NGS identified all 17 IMC events detected by qPCR (100% sensitivity). The NGS chimerism provides sufficient sensitivity, accuracy, and economical/logistical feasibility in supporting engraftment and MRD monitoring.

Genetic sex validation for sample tracking in next-generation sequencing clinical testing

OBJECTIVE

Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups.

RESULTS

Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors (49.09%), samples from transgender participants (3.64%) and stem cell or bone marrow transplant patients (7.27%) along with undetermined sample mix-ups (40%) for which sample swaps occurred prior to arrival at genome centers, however the exact cause of the events at the sampling sites resulting in the mix-ups were not able to be determined.

Population genetic analyses of Eastern Chinese Han nationality using ForenSeq™ DNA Signature Prep Kit

Currently, the most commonly used method for human identification and kinship analysis in forensic genetics is the detection of length polymorphism in short tandem repeats (STRs) using polymerase chain reaction (PCR) and capillary electrophoresis (CE). However, numerous studies have shown that considerable sequence variations exist in the repeat and flanking regions of the STR loci, which cannot be identified by CE detection. Comparatively, massively parallel sequencing (MPS) technology can capture these sequence differences, thereby enhancing the identification capability of certain STRs. In this study, we used the ForenSeq™ DNA Signature Prep Kit to sequence 58 STRs and 94 individual identification SNPs (iiSNPs) in a sample of 220 unrelated individuals from the Eastern Chinese Han population. Our aim is to obtain MPS-based STR and SNP data, providing further evidence for the study of population genetics and forensic applications. The results showed that the MPS method, utilizing sequence information, identified a total of 486 alleles on autosomal STRs (A-STRs), 97 alleles on X-chromosome STRs (X-STRs), and 218 alleles on Y-chromosome STRs (Y-STRs). Compared with length polymorphism, we observed an increase of 260 alleles (157, 31, and 72 alleles on A-STRs, X-STRs, and Y-STRs, respectively) across 36 STRs. The most substantial increments were observed in DYF387S1 and DYS389II, with increases of 287.5% and 250%, respectively. The most increment in the number of alleles was found at DYF387S1 and DYS389II (287.5% and 250%, respectively). The length-based (LB) and sequence-based (SB) combined random match probability (RMP) of 27 A-STRs were 6.05E-31 and 1.53E-34, respectively. Furthermore, other forensic parameters such as total discrimination power (TDP), cumulative probability of exclusion of trios (CPEtrio), and duos (CPEduo) were significantly improved when using the SB data, and informative data were obtained for the 94 iiSNPs. Collectively, these findings highlight the advantages of MPS technology in forensic genetics, and the Eastern Chinese Han genetic data generated in this study could be used as a valuable reference for future research in this field.

Solving the Arizona search problem by imputation

An "Arizona search" is an evaluation of the numbers of pairs of profiles in a forensic-genetic database that possess partial or complete genotypic matches; such a search assists in establishing the extent to which a set of loci provides unique identifications. In forensic genetics, however, the potential for performing Arizona searches is constrained by the limited availability of actual forensic profiles for research purposes. Here, we use genotype imputation to circumvent this problem. From a database of genomes, we impute genotypes of forensic short-tandem-repeat (STR) loci from neighboring single-nucleotide polymorphisms (SNPs), searching for partial STR matches using the imputed profiles. We compare the distributions of the numbers of partial matches in imputed and actual profiles, finding close agreement. Despite limited potential for performing Arizona searches with actual forensic STR profiles, the questions that such searches seek to answer can be posed with imputation-based Arizona searches in increasingly large SNP databases.

Feasibility for non-invasive prenatal fetal blood group and platelet genotyping by massively parallel sequencing: A single test system for multiple atypical red cell, platelet and quality control markers

Non-invasive prenatal tests (NIPT) to predict fetal red cell or platelet antigen status for alloimmunised women are provided for select antigens. This study reports on massively parallel sequencing (MPS) using a red cell and platelet probe panel targeting multiple nucleotide variants, plus individual identification single nucleotide polymorphisms (IISNPs). Maternal blood samples were provided from 33 alloimmunised cases, including seven with two red cell antibodies. Cell-free and genomic DNA was sequenced using targeted MPS and bioinformatically analysed using low-frequency variant detection. The resulting maternal genomic DNA allele frequency was subtracted from the cell-free DNA counterpart. Outcomes were matched against validated phenotyping/genotyping methods, where available. A 2.5% subtractive allele frequency threshold was set after comparing MPS predictions for K, RhC/c, RhE/e and Fya /Fyb against expected outcomes. This threshold was used for subsequent predictions, including HPA-15a, Jka /Jkb , Kpa /Kpb and Lua . MPS outcomes were 97.2% concordant with validated methods; one RhC case was discordantly negative and lacked IISNPs. IISNPs were informative for 30/33 cases as controls. NIPT MPS is feasible for fetal blood group genotyping and covers multiple blood groups and control targets in a single test. Noting caution for the Rh system, this has the potential to provide a personalised service for alloimmunised women.

Genetic landscape of forensic DNA phenotyping markers among Mediterranean populations

Forensic DNA Phenotyping can reveal the appearance of an unknown individual by predicting the External Visible Characteristics (EVC) from DNA obtained at the crime scene. Our aim is to characterize the genetic landscape of Human identification markers responsible for EVC among Mediterranean populations compared to other worldwide groups. We conducted an exhaustive search for genes involved in EVC variation. Then, variants located on these genes were extracted from public genotypic data of Mediterranean, American, African and East Asiatic populations. The genetic landscape of these Human identification markers, their allelic distribution and admixture analyses, were determined using plink, R and ADMIXTURE softwares. Our results showed that the Mediterranean populations appear close to the Mexican populations and distinguished from sub Saharan African populations living in the USA and from East Asiatic populations. We highlighted a total of 103454 common variants shared between the studied populations and among them, 25 common variants associated with EVC. Interestingly, genotype frequencies results showed that the rs17646946, rs13016869, rs977588, rs1805008 and rs2240751 variants located respectively in the TCHH, PRKCE, OCA2, MC1R and MFSD12 genes are significantly different between the Mediterranean and Asiatic populations. The genotype frequencies of the variants rs977589 and rs7179994 located in the OCA2 gene, and of rs12913832 and rs2240751 located respectively in HERC2 and MFSD12 genes are significantly different between the Mediterranean and American populations. Our work generates a large number of EVC variants that could be a valuable resource for future studies in the forensic field.

Insertion/deletion polymorphism for genetic background and forensic performance exploration of the Sui group from Guizhou

Insertion/deletion polymorphisms (InDels) as ideal genetic markers for forensic genetics are appreciated by scholars both nationally and internationally because they integrated the favorable features of single nucleotide polymorphisms (SNPs) and short tandem repeats (STRs). Nevertheless, with the limited identification efficiency of InDels, the multiplex amplification systems of InDels might just be applied as the supplementary methods in paternity testing with respect to commonly used STRs. In the current research, we successfully genotyped 105 unrelated individuals from the Guizhou Sui population based on a six-color fluorescence multiplex panel that could simultaneously detect 64 genetic markers (59 autosomal InDels, two autosomal miniSTRs and three Y chromosomal genetic markers). In addition, frequency distributions and forensic statistical parameters of these loci in the Sui group were assessed using the STRAF software. Phylogenetic relationships among the Sui group and other reference populations were dissected by two methods (principal component analysis and phylogenetic trees) based on 59 InDels. The combined discrimination power and probability of exclusion values of 61 autosomal genetic markers in the Sui group were nearly equal to 1-1.90063 × 10-27 and 0.999998272, respectively. Furthermore, we observed that the Sui group from Guizhou had closer genetic affinities with East Asian populations with respect to other continental populations. In summary, we stated that the multiplex amplification system might be utilized as a prospective independent tool for human individual identification and parentage testing in the Sui group residing in Guizhou.

Forensic genetic analysis of single-nucleotide polymorphisms and microhaplotypes in Koreans through next-generation sequencing using precision ID identity panel

BACKGROUND

Forensic DNA analysis has seen remarkable advancements with the advent of Next Generation Sequencing (NGS). In particular, NGS analysis of single nucleotide polymorphisms (SNPs) offers significant advantages in the analysis of challenging samples compared to conventional STR analysis.

OBJECTIVE

This study aimed to investigate the SNPs of the Precision ID Identity Panel, a commercially available NGS panel for personal identification, by generating genetic profiles of 298 Koreans and comparing them with other global populations.

METHODS

A total of 124 SNPs, including 90 autosomal and 34 Y-SNPs, were analyzed using the Precision ID Identity Panel, and forensic parameters, microhaplotypes, and population differences were investigated.

RESULTS

The NGS data were successfully obtained from 298 Koreans. The analysis of forensic parameters exhibited a low combined match probability of 1.532 × 10- 34, which is comparable to that obtained from commonly used STR analysis. Additionally, the microhaplotype analysis revealed that the use of 16 microhaplotypes provided higher discriminatory power compared to single target SNPs. Furthermore, the adoption of microhaplotype data resulted in an increase of over 20% in expected heterozygosity at five loci. Inter-population analysis showed a close genetic relationship between Koreans and individuals from China and Myanmar in East and Southeast Asia, which are geographically adjacent to Korea.

CONCLUSIONS

The results of this study show that the Precision ID Identity panel can be a useful alternative where traditional STR typing is not feasible. Also, the data from our study will be useful as a reference for Koreans in forensic investigations and the prosecution of criminal justice.

Genetic Sex Validation for Sample Tracking in Clinical Testing

Objective

Data from DNA genotyping via a 96-SNP panel in a study of 25,015 clinical samples were utilized for quality control and tracking of sample identity in a clinical sequencing network. The study aimed to demonstrate the value of both the precise SNP tracking and the utility of the panel for predicting the sex-by-genotype of the participants, to identify possible sample mix-ups.

Results

Precise SNP tracking showed no sample swap errors within the clinical testing laboratories. In contrast, when comparing predicted sex-by-genotype to the provided sex on the test requisition, we identified 110 inconsistencies from 25,015 clinical samples (0.44%), that had occurred during sample collection or accessioning. The genetic sex predictions were confirmed using additional SNP sites in the sequencing data or high-density genotyping arrays. It was determined that discrepancies resulted from clerical errors, samples from transgender participants and stem cell or bone marrow transplant patients along with undetermined sample mix-ups.

Genomic screening methodology not requiring barcoding: Single nucleotide polymorphism-based, mixed-cell screening (SMICS)

Although high-throughput, cancer cell-line screening is a time-honored, important tool for anti-cancer drug development, this process involves the testing of each, individual drug in each, individual cell-line. Despite the availability of robotic liquid handling systems, this process remains a time-consuming and costly investment. The Broad Institute developed a new method called Profiling Relative Inhibition Simultaneously in Mixtures (PRISM) to screen a mixture of barcoded, tumor cell-lines. Although this methodology significantly improved the efficiency of screening large numbers of cell-lines, the barcoding process itself was tedious that requires gene transfection and subsequent selection of stable cell-lines. In this study, we developed a new, genomic approach for screening multiple cancer cell-lines using endogenous "tags" that did not require prior barcoding: single nucleotide polymorphism-based, mixed-cell screening (SMICS). The code for SMICS is available at https://github.com/MarkeyBBSRF/SMICS.

US Population Data for 94 Identity-Informative SNP Loci

The US National Institute of Standards and Technology (NIST) analyzed a set of 1036 samples representing four major US population groups (African American, Asian American, Caucasian, and Hispanic) with 94 single nucleotide polymorphisms (SNPs) used for individual identification (iiSNPs). The compact size of iiSNP amplicons compared to short tandem repeat (STR) markers increases the likelihood of successful amplification with degraded DNA samples. Allele frequencies and relevant forensic statistics were calculated for each population group as well as the aggregate population sample. Examination of sequence data in the regions flanking the targeted SNPs identified additional variants, which can be combined with the target SNPs to form microhaplotypes (multiple phased SNPs within a short-read sequence). Comparison of iiSNP performance with and without flanking SNP variation identified four amplicons containing microhaplotypes with observed heterozygosity increases of greater than 15% over the targeted SNP alone. For this set of 1036 samples, comparison of average match probabilities from iiSNPs with the 20 CODIS core STR markers yielded an estimate of 1.7 × 10^-38 for iiSNPs (assuming independence between all 94 SNPs), which was four orders of magnitude lower (more discriminating) than STRs where internal sequence variation was considered, and 10 orders of magnitude lower than STRs using established capillary electrophoresis length-based genotypes.

Genetic investigation of 124 SNPs in a Myanmar population using the Precision ID Identity Panel and the Illumina MiSeq

Single nucleotide polymorphisms (SNPs) have become popular in forensic genetics as an alternative to short tandem repeats (STRs). The Precision ID Identity Panel (Thermo Fisher Scientific), consisting of 90 autosomal SNPs and 34 Y-chromosomal SNPs, enabled human identification studies on global populations through next-generation sequencing (NGS). However, most previous studies on the panel have used the Ion Torrent platform, and there are few reports on the Southeast Asian population. Here, a total of 96 unrelated males from Myanmar (Yangon) were analyzed with the Precision ID Identity Panel on a MiSeq (Illumina) using an in-house TruSeq compatible universal adapter and a custom variant caller, Visual SNP. The sequencing performance evaluated by locus balance and heterozygote balance was comparable to that of the Ion Torrent platform. For 90 autosomal SNPs, the combined match probability (CMP) was 6.994 × 10^-34, lower than that of 22 PowerPlex Fusion autosomal STRs (3.130 × 10^-26). For 34 Y-SNPs, 14 Y-haplogroups (mostly O2 and O1b) were observed. We found 51 cryptic variations (42 haplotypes) around target SNPs, of which haplotypes corresponding to 33 autosomal SNPs decreased CMP. Interpopulation analysis revealed that the Myanmar population is genetically closer to the East and Southeast Asian populations. In conclusion, the Precision ID Identity Panel can be successfully analyzed on the Illumina MiSeq and provides high discrimination power for human identification in the Myanmar population. This study broadened the accessibility of the NGS-based SNP panel by expanding the available NGS platforms and adopting a robust NGS data analysis tool.

Genetic Ethnic Differences in Human 2'-5'-Oligoadenylate Synthetase and Disease Associations: A Systematic Review

Recently, several studies have highlighted a skewed prevalence of infectious diseases within the African continent. Furthermore, a growing number of studies have demonstrated unique genetic variants found within the African genome are one of the contributing factors to the disease severity of infectious diseases within Africa. Understanding the host genetic mechanisms that offer protection against infectious diseases provides an opportunity to develop unique therapeutic interventions. Over the past two decades, several studies have linked the 2'-5'-oligoadenylate synthetase (OAS) family with a range of infectious diseases. More recently, the OAS-1 gene has also been associated with disease severity caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), which led to a global pandemic. The OAS family serves as an antiviral factor through the interaction with Ribonuclease-Latent (RNase-L). This review explores the genetic variants observed within the OAS genes and the associations with various viral infections and how previously reported ethnic-specific polymorphisms drive clinical significance. This review provides an overview of OAS genetic association studies with a particular focus on viral diseases affecting individuals of African descent.

Privacy-aware estimation of relatedness in admixed populations

BACKGROUND

Estimation of genetic relatedness, or kinship, is used occasionally for recreational purposes and in forensic applications. While numerous methods were developed to estimate kinship, they suffer from high computational requirements and often make an untenable assumption of homogeneous population ancestry of the samples. Moreover, genetic privacy is generally overlooked in the usage of kinship estimation methods. There can be ethical concerns about finding unknown familial relationships in third-party databases. Similar ethical concerns may arise while estimating and reporting sensitive population-level statistics such as inbreeding coefficients for the concerns around marginalization and stigmatization.

RESULTS

Here, we present SIGFRIED, which makes use of existing reference panels with a projection-based approach that simplifies kinship estimation in the admixed populations. We use simulated and real datasets to demonstrate the accuracy and efficiency of kinship estimation. We present a secure federated kinship estimation framework and implement a secure kinship estimator using homomorphic encryption-based primitives for computing relatedness between samples in two different sites while genotype data are kept confidential. Source code and documentation for our methods can be found at https://doi.org/10.5281/zenodo.7053352.

CONCLUSIONS

Analysis of relatedness is fundamentally important for identifying relatives, in association studies, and for estimation of population-level estimates of inbreeding. As the awareness of individual and group genomic privacy is growing, privacy-preserving methods for the estimation of relatedness are needed. Presented methods alleviate the ethical and privacy concerns in the analysis of relatedness in admixed, historically isolated and underrepresented populations.

SHORT ABSTRACT

Genetic relatedness is a central quantity used for finding relatives in databases, correcting biases in genome wide association studies and for estimating population-level statistics. Methods for estimating genetic relatedness have high computational requirements, and occasionally do not consider individuals from admixed ancestries. Furthermore, the ethical concerns around using genetic data and calculating relatedness are not considered. We present a projection-based approach that can efficiently and accurately estimate kinship. We implement our method using encryption-based techniques that provide provable security guarantees to protect genetic data while kinship statistics are computed among multiple sites.

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.

Forensic and genetic landscape explorations of Chinese Kyrgyz group based on autosomal SNPs, Y-chromosomal SNPs and STRs

To assess the effect of population genetic polymorphism on forensic research, we investigated the genetic polymorphisms of Chinese Kyrgyz group (n = 98) and evaluated forensic application values in Chinese Kyrgyz group and other 26 reference populations at 90 autosomal SNPs, and then combined with 34 SNPs and 37 STRs on Y chromosome to reveal the genetic background of Kyrgyz group in autosomal and Y-chromosomal inheritances, respectively. The 90 autosomal SNPs and 34 Y-chromosomal SNPs were sequenced base on next generation sequencing technology, and 37 Y-chromosomal STRs were analyzed by capillary electrophoresis platform. The results showed that cumulative power of discrimination and cumulative power of exclusion of 90 autosomal SNPs in the panel met the application need of forensic genetics in Kyrgyz group. The forensic effectivenesses of the panel were high in all 27 populations, although there were genetic differences among these populations. The forensic effectiveness of the panel was relatively higher in the European populations, but relatively lower in the African populations. The population genetic results indicated that the Kyrgyz group had the relatively closer genetic relationships with the reference East Asian populations at autosomal SNPs, and there were gene exchanges between the Kyrgyz group and East Asian, European populations based on the analytical results of autosomal SNPs, Y-chromosomal SNPs and STRs.

P2GT: Fine-Grained Genomic Data Access Control With Privacy-Preserving Testing in Cloud Computing

With the rapid development of bioinformatics and the availability of genetic sequencing technologies, genomic data has been used to facilitate personalized medicine. Cloud computing, features as low cost, rich storage and rapid processing can precisely respond to the challenges brought by the emergence of massive genomic data. Considering the security of cloud platform and the privacy of genomic data, we first introduce P2GT which utilizes key-policy attribute-based encryption to realize genomic data access control with unbounded attributes, and employs equality test algorithm to achieve personalized medicine test by matching digitized single nucleotide polymorphisms (SNPs) directly on the users' ciphertext without encrypting multiple times. We then propose an enhanced scheme P2GT+, which adopts identity-based encryption with equality test supporting flexible joint authorization to realize privacy-preserving paternity test, genetic compatibility test and disease susceptibility test over the encrypted SNPs with P2GT. We prove the security of proposed schemes and conduct extensive experiments with the 1,000 Genomes dataset. The results show that P2GT and P2GT+ are practical and scalable enough to meet the privacy-preserving and authorized genetic testing requirements in cloud computing.

A New Computational Deconvolution Algorithm for the Analysis of Forensic DNA Mixtures with SNP Markers

Single nucleotide polymorphisms (SNPs) support robust analysis on degraded DNA samples. However, the development of a systematic method to interpret the profiles derived from the mixtures is less studied, and it remains a challenge due to the bi-allelic nature of SNP markers. To improve the discriminating power of SNPs, this study explored bioinformatic strategies to analyze mixtures. Then, computer-generated mixtures were produced using real-world massively parallel sequencing (MPS) data from the single samples processed with the Precision ID Identity Panel. Moreover, the values of the frequency of major allele reads (F_MAR) were calculated and applied as key parameters to deconvolve the two-person mixtures and estimate mixture ratios. Four custom R language scripts (three for autosomes and one for Y chromosome) were designed with the K-means clustering method as a core algorithm. Finally, the method was validated with real-world mixtures. The results indicated that the deconvolution accuracy for evenly balanced mixtures was 100% or close to 100%, which was the same as the deconvolution accuracy of inferring the genotypes of the major contributor of unevenly balanced mixtures. Meanwhile, the accuracy of inferring the genotypes of the minor contributor decreased as its proportion in the mixture decreased. Moreover, the estimated mixture ratio was almost equal to the actual ratio between 1:1 and 1:6. The method proposed in this study provides a new paradigm for mixture interpretation, especially for inferring contributor profiles of evenly balanced mixtures and the major contributor profile of unevenly balanced mixtures.

Systematic Evaluation of a Novel 6-dye Direct and Multiplex PCR-CE-Based InDel Typing System for Forensic Purposes

Insertion/deletion (InDel) polymorphisms, combined desirable characteristics of both short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs), are considerable potential in the fields of forensic practices and population genetics. However, most commercial InDel kits designed based on non-Asians limited extensive forensic applications in East Asian (EAS) populations. Recently, a novel 6-dye direct and multiplex PCR-CE-based typing system was designed on the basis of genome-wide EAS population data, which could amplify 60 molecular genetic markers, consisting of 57 autosomal InDels (A-InDels), 2 Y-chromosomal InDels (Y-InDels), and Amelogenin in a single PCR reaction and detect by capillary electrophoresis, simultaneously. In the present study, the DNA profiles of 279 unrelated individuals from the Hainan Li group were generated by the novel typing system. In addition, we collected two A-InDel sets to evaluate the forensic performances of the novel system in the 1,000 Genomes Project (1KG) populations and Hainan Li group. For the Universal A-InDel set (UAIS, containing 44 A-InDels) the cumulative power of discrimination (CPD) ranged from 1-1.03 × 10-14 to 1-1.27 × 10-18, and the cumulative power of exclusion (CPE) varied from 0.993634 to 0.999908 in the 1KG populations. For the East Asia-based A-InDel set (EAIS, containing 57 A-InDels) the CPD spanned from 1-1.32 × 10-23 to 1-9.42 × 10-24, and the CPE ranged from 0.999965 to 0.999997. In the Hainan Li group, the average heterozygote (He) was 0.4666 (0.2366-0.5448), and the polymorphism information content (PIC) spanned from 0.2116 to 0.3750 (mean PIC: 0.3563 ± 0.0291). In total, the CPD and CPE of 57 A-InDels were 1-1.32 × 10-23 and 0.999965, respectively. Consequently, the novel 6-dye direct and multiplex PCR-CE-based typing system could be considered as the reliable and robust tool for human identification and intercontinental population differentiation, and supplied additional information for kinship analysis in the 1KG populations and Hainan Li group.

Performance characteristics of chimerism testing by next generation sequencing

Chimerism testing provides informative clinical data regarding the status of a biological sample mixture. For years, this testing was achieved by measuring the peaks of informative short tandem repeat (STR) loci using capillary electrophoresis (CE). With the advent of next generation sequencing (NGS) technology, the quantification of the percentage of donor/recipient mixtures is more easily done using sequence reads in large batches of samples run on a single flow cell. In this study, we present data on using a FORENSIC NGS chimerism platform to accurately measure the percentage of donor/recipient mixtures. We were able to detect chimerism to a limit threshold of 1% using both STR and single nucleotide polymorphism (SNP) informative loci. Importantly, a significant correlation was observed between NGS and CE chimerism methods when compared at donor detection ranges from 1% to 10%. Furthermore, 100% accuracy was achieved through proficiency testing over six surveys. Its usefulness was expanded beyond this to help identify suitable donors for solid organ transplant patients using ancestry SNP profiles. In summary, the NGS method provides a sensitive and reliable alternative to traditional CE for chimerism testing of clinical samples.

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.

Donor-derived Cell-free DNA in Solid-organ Transplant Diagnostics: Indications, Limitations, and Future Directions

The last few years have seen an explosion in clinical research focusing on the use of donor-derived cell-free DNA (dd-cfDNA) in solid-organ transplants (SOT). Although most of the literature published so far focuses on kidney transplants, there are several recent as well as ongoing research studies on heart, lung, pancreas, and liver transplants. Though initially studied as a noninvasive means of identifying subclinical or acute rejection in SOT, it is rapidly becoming clear that instead of being a specific marker for allograft rejection, dd-cfDNA is more appropriately described as a marker of severe injury, although the most common cause of this injury is allograft rejection. Multiple studies in kidney transplants have shown that although sensitivity for the diagnosis of antibody-mediated rejection is excellent, it is less so for T-cell-mediated rejection. It is possible that combining dd-cfDNA with other novel urine- or blood-based biomarkers may increase the sensitivity for the diagnosis of rejection. Irrespective of the cause, though, elevated dd-cfDNA seems to portend adverse allograft prognosis and formation of de novo donor-specific antibody. Although current data do not lend themselves to a clear conclusion, ongoing studies may reveal the utility of serial surveillance for the management of SOT as following levels of dd-cfDNA over time may provide windows of opportunity to intervene early and before irreversible allograft injury. Finally, cost-effectiveness studies will be needed to guide the ideal incorporation of dd-cfDNA into routine clinical practice.

Human sample authentication in biomedical research: comparison of two platforms

Samples used in biomedical research are often collected over years, in some cases from subjects that may have died and thus cannot be retrieved in any way. The value of these samples is priceless. Sample misidentification or mix-up are unfortunately common problems in biomedical research and can eventually result in the publication of incorrect data. Here we have compared the Fluidigm SNPtrace and the Agena iPLEX Sample ID panels for the authentication of human genomic DNA samples. We have tested 14 pure samples and simulated their cross-contamination at different percentages (2%, 5%, 10%, 25% and 50%). For both panels, we report call rate, allele intensity/probability score, performance in distinguishing pure samples and contaminated samples at different percentages, and sex typing. We show that both panels are reliable and efficient methods for sample authentication and we highlight their advantages and disadvantages. We believe that the data provided here is useful for sample authentication especially in biorepositories and core facility settings.

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.

Selection and evaluation of bi-allelic autosomal SNP markers for paternity testing in Koreans

Due to the advantages of single-nucleotide polymorphisms (SNPs) in forensic science, many forensic SNP panels have been developed. However, the existing SNP panels have a problem that they do not reflect allele frequencies in Koreans or the number of markers is not sufficient to perform paternity testing. Here, we filtered candidate SNPs from the Ansan-Ansung cohort data and selected 200 SNPs with high allele frequencies. To reduce the risk of false inclusion and false exclusion, we calculated likelihood ratios of alleged father-child pairs from simulated families when the alleged father is the true father, the close relative of the true father, and the random man. As a result, we estimated that 160 SNPs were needed to perform paternity testing. Furthermore, we performed validation using Twin-Family cohort data. When 160 selected SNPs were used to calculate the likelihood ratio, paternity and non-paternity were accurately distinguished. Our set of 160 SNPs could be useful for paternity testing in Koreans.

Parallel sequencing of 87 STR and 294 SNP markers using the prototype of the SifaMPS panel on the MiSeq FGx™ system

Massively parallel sequencing (MPS), or next generation sequencing (NGS), is a promising methodology for the detection of short tandem repeats (STRs) and single nucleotide polymorphisms (SNPs) in forensic genetics. Here, the prototype SifaMPS Panel is designed to simultaneously target 87 STRs and 294 SNPs with forensic interest in a single multiplex in conjunction with the TruSeq™ Custom Amplicon workflow and MiSeq FGx™ System. Two in-house python scripts are adopted for the fastq-to-genotype interpretation of MPS data concerning STR and SNP, respectively. In the present study, by sequencing 50 Chinese Hans and many other DNA samples involved in validation studies, system parameters including the depth of coverage (DoC), heterozygote balance (Hb) and sequence coverage ratios (SCRs), as well as different forensic parameters of STRs and SNPs in a population study, were calculated to evaluate the overall performance of this new panel and its practicality in forensic application. In general, except for two STRs (DYS505 and DYS449) and one SNP (rs4288409) that performed poorly, the other 85 STRs and 293 SNPs in our panel had good performance that could strengthen efficiency for human identification and paternity testing. In addition, discordant STR genotype results between those generated from capillary electrophoresis (CE) and from the MPS platform were clearly illustrated, and these results could be a useful reference for applying these particular non-CODIS STRs in forensic practice.

An alternate workflow for preparing Precision ID Ancestry and Identity Panel libraries for Illumina sequencing

Single nucleotide polymorphisms (SNPs) are well-established for forensic applications. Although they are not compatible with existing criminal databases, they offer some advantages over short tandem repeat (STR) markers including smaller amplicons, no stutter artifacts, and biogeographic ancestry and phenotype predictions. The Precision ID NGS System, a commercial workflow by Thermo Fisher Scientific, offers a streamlined solution for genotyping forensically relevant SNPs using next-generation sequencing. The Precision ID Ancestry and Identity Panels combined target 289 SNPs, and their sensitivity, reproducibility, and accuracy have been evaluated by the forensic community. The aim of this study was to develop an alternative workflow to genotype these SNP panels using Illumina chemistry. Commercial genomic DNAs (gDNAs) (n, 3) were amplified using three uracil-tolerant polymerase master mixes. Resulting amplicons were prepared into libraries using the KAPA Hyper Prep Kit (KAPA Biosystems) and sequenced via Illumina's MiniSeq. Reads were analyzed using a published analysis pipeline to compile final genotypes with read depth information. Phusion U Multiplex PCR Master Mix (Thermo Fisher Scientific) statistically outperformed the other master mixes tested (P <0.0001), with respect to the number of SNPs genotyped. To ensure a workflow using Phusion U would be compatible across diverse samples, we optimized PCR cycle number using the same commercial gDNAs (n, 3), reference buccal swabs (n, 3), and environmental (household dust) samples (n, 6). Using the developed workflow, 93.9% of all SNPs were successfully genotyped across sample types. Implementation of the developed workflow should be straightforward for forensic laboratories and suitable for processing reference and casework samples.

Donor-derived Cell-free DNA: Advancing a Novel Assay to New Heights in Renal Transplantation

Despite advances in transplant immunosuppression, long-term renal allograft outcomes remain suboptimal because of the occurrence of rejection, recurrent disease, and interstitial fibrosis with tubular atrophy. This is largely due to limitations in our understanding of allogeneic processes coupled with inadequate surveillance strategies. The concept of donor-derived cell-free DNA as a signal of allograft stress has therefore rapidly been adopted as a noninvasive monitoring tool. Refining it for effective clinical use, however, remains an ongoing effort. Furthermore, its potential to unravel new insights in alloimmunity through novel molecular techniques is yet to be realized. This review herein summarizes current knowledge and active endeavors to optimize cell-free DNA-based diagnostic techniques for clinical use in kidney transplantation. In addition, the integration of DNA methylation and microRNA may unveil new epigenetic signatures of allograft health and is also explored in this report. Directing research initiatives toward these aspirations will not only improve diagnostic precision but may foster new paradigms in transplant immunobiology.

Evaluation of OpenArray™ as a Genotyping Method for Forensic DNA Phenotyping and Human Identification

A custom plate of OpenArray™ technology was evaluated to test 60 single-nucleotide polymorphisms (SNPs) validated for the prediction of eye color, hair color, and skin pigmentation, and for personal identification. The SNPs were selected from already validated subsets (Hirisplex-s, Precision ID Identity SNP Panel, and ForenSeq DNA Signature Prep Kit). The concordance rate and call rate for every SNP were calculated by analyzing 314 sequenced DNA samples. The sensitivity of the assay was assessed by preparing a dilution series of 10.0, 5.0, 1.0, and 0.5 ng. The OpenArray™ platform obtained an average call rate of 96.9% and a concordance rate near 99.8%. Sensitivity testing performed on serial dilutions demonstrated that a sample with 0.5 ng of total input DNA can be correctly typed. The profiles of the 19 SNPs selected for human identification reached a random match probability (RMP) of, on average, 10⁻⁸. An analysis of 21 examples of biological evidence from 8 individuals, that generated single short tandem repeat profiles during the routine workflow, demonstrated the applicability of this technology in real cases. Seventeen samples were correctly typed, revealing a call rate higher than 90%. Accordingly, the phenotype prediction revealed the same accuracy described in the corresponding validation data. Despite the reduced discrimination power of this system compared to STR based kits, the OpenArray™ System can be used to exclude suspects and prioritize samples for downstream analyses, providing well-established information about the prediction of eye color, hair color, and skin pigmentation. More studies will be needed for further validation of this technology and to consider the opportunity to implement this custom array with more SNPs to obtain a lower RMP and to include markers for studies of ancestry and lineage.

Development and performance of a next generation sequencing (NGS) assay for monitoring of mixed chimerism

The aim of this study was to evaluate the performance of a novel NGS-based assay to monitor mixed chimerism (MC) and compare its technical capacity to established techniques for chimerism analysis. Artificial and clinical samples with increasing amounts of patient DNA were compared using real-time PCR detection of indels and SNP, fragment analysis of short-tandem repeats (STR) and NGS analysis of indels. Real-time PCR displayed excellent sensitivity (>0,01%) but poor accuracy (>20 CV% at MC > 20%), while fragment analysis exhibited good accuracy (<5 CV% at MC > 20%) with limited sensitivity (>2,5%). In contrast, NGS chimerism demonstrated a sensitivity (>0,1%) equal to real-time PCR and an accuracy equal or better than STR analysis throughout an extensive range of mixed chimerism (0,1 - 100%). To evaluate performance of the separate techniques for chimerism determination, 75 retrospective patient monitoring samples (3-7 weeks post-HSCT) with low (<5%), intermediate (5-20%) or high mixed chimerism (>20%) were analyzed. The between run precision for the NGS assay varied from 0,72% (>20% MC) to 7,38% (MC < 5%). In conclusion, NGS displayed a combination of high sensitivity with good accuracy in both artificial and clinical chimerism samples.

The first GHEP-ISFG collaborative exercise on forensic applications of massively parallel sequencing

One of the main goals of the Spanish and Portuguese-Speaking Working Group of the International Society for Forensic Genetics (GHEP-ISFG) is to promote and contribute to the development and dissemination of scientific knowledge in the field of forensic genetics. The GHEP-ISFG supports several Working Commissions which develop different scientific activities. One of them, the Working Commission on "Massively Parallel Sequencing (MPS): Forensic Applications", organized its first collaborative exercise on forensic applications of MPS technology in 2019. The aim of this exercise was to assess the concordance between the MPS results and those obtained with conventional technologies (capillary electrophoresis and Sanger sequencing), as well as to compare the results obtained within the different MPS platforms and/or the different kits/panels and analysis software packages (commercial and open-access) available on the market. The seven participating laboratories analyzed some samples of the annual GHEP-ISFG proficiency test (EIADN No. 27 (2019)), using Ion Torrent™ or MiSeq FGx® platforms. Six of them sent autosomal STR sequence data, five laboratories performed MPS analysis of individual identification SNPs, four laboratories reported MPS data of Y-chromosomal STRs, and X-chromosomal STRs, three laboratories performed MPS analysis of ancestry informative SNPs and phenotype informative SNPs, two labs performed MPS analysis of the mitochondrial DNA control region, and only one lab produced MPS data of lineage informative SNPs. Autosomal STR sequencing results were highly concordant to the consensus obtained by capillary electrophoresis in the EIADN No. 27 (2019) exercise. Furthermore, in general, a high level of concordance was observed between the results of the participating laboratories, regardless of the platform used. The main discordances were due to errors during the analysis process or from sequence data obtained with low depth of coverage. In this paper we highlight some issues that still arise, such as standardization of the nomenclature for STRs analyzed by sequencing with MPS, the universal uptake of a nomenclature framework by the analysis software, and well established validation and accreditation of the new MPS platforms for use in routine forensic case-work.

A novel forensic panel of 186-plex SNPs and 123-plex STR loci based on massively parallel sequencing

The MiSeq® FGX Forensic system and the HID-Ion AmpliSeq Panel were previously developed for massively parallel sequencing (MPS) for forensic casework. Among the three major sequencing platforms, BGISEQ-500^TM, which is based on multiple PCRs, is still lacking in forensics. Here, a novel forensic panel was constructed to detect 186 single-nucleotide polymorphisms (SNPs) and 123 short tandem repeats (STRs) with MPS technology on the BGISEQ-500™ platform. First, the library preparation, sequencing process, and data analysis were performed, focusing on the average depth of coverage and heterozygote balance. We calculated the allelic frequencies and forensic parameters of STR and SNP loci in 73 unrelated Chinese Han individuals. In addition, performance was evaluated with accuracy, uniformity, sensitivity, PCR inhibitor, repeatability and reproducibility, mixtures, degraded samples, case-type samples, and pedigree analyses. The results showed that 100% accurate and concordant genotypes can be obtained, and the loci with an abundance in the interquartile range accounted for 92.90% of the total, suggesting reliable uniformity in this panel. We obtained a locus detection rate that was higher than 98.78% from 78 pg of input DNA, and the optimal amount was 1.25-10 ng. The maximum concentrations of hematin and humic acid were 200 and 100 μM, respectively (the ratios of detected loci were 96.52% and 92.41%), in this panel. As a mixture, compared with those of SNPs, minor-contributor alleles of STRs could be detected at higher levels. For the degraded sample, the ratio of detected loci was 98.41%, and most profiles from case-type samples were not significantly different in abundance in our studies. As a whole, this panel showed high-performance, reliable, robust, repeatable, and reproducible results, which are sufficient for paternity testing, individual identification, and use for potentially degraded samples in forensic science.

Validation of novel forensic DNA markers using multiplex microhaplotype sequencing

Microhaplotypes (MH) are comprised of multiple single nucleotide polymorphisms (SNPs) that are located within 300 bases of genomic sequence. Improved tools are needed to facilitate broader application of microhaplotypes in a diverse range of populations and forensic settings. We designed an assay for multiplex sequencing of 90 microhaplotypes (mMHseq) that include 46 MH loci with high Effective Number of Alleles (A_e) from previous studies [1], and 44 high A_e MH loci containing between four to fourteen SNPs that were identified from the 1000 Genomes (1KG) Project. The unique design of mMHseq integrates a novel method for multiplex amplification from small DNA amounts, and multiplex sequencing of 48 samples in a single MiSeq run to detect all relevant MH variation. Assay performance was evaluated in a cohort of 156 individuals from seven different world populations from Africa, Asia, and Europe. Three of those populations from East Africa (Chagga, Sandawe, and Zaramo) and one from Eastern Europe (Adygei) had sufficient individuals sequenced by the assay to be included in statistical analyses with the 26 1KG populations. For those 30 populations the mean global average A_e was 5.08 (range: 2.7-11.54) and mean informativeness for biogeographic variation (I_n) was 0.30 (range: 0.08-0.70). Eighty-five novel SNPs were detected in 58 of the 90 microhaplotypes. Open-source, web-based software was developed to visualize haplotype phase data for each microhaplotype and individual. Our approach for multiplex microhaplotype sequencing can be customized and expanded as novel loci are being discovered.

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.

Next-generation sequencing methods to detect donor-derived cell-free DNA after transplantation

Following the initial technical challenge of implanting an organ, maintaining the organ against a vast array of pathologies for years to come, remains a colossal challenge for all clinicians working in transplantation. Drug toxicity, opportunistic infection, primary disease recurrence, and the constant battle against organ rejection are all differentials that are considered when graft dysfunction is observed, promoting a lifetime of laborious surveillance. Cell free DNA (cfDNA) since its discovery in 1948 has made an impactful change in transplantation. A growing body of evidence in transplantation (109 manuscripts from 55 studies) shows the promise of this tool as an early and accurate detection of allograft injury rejection as well the benefit to rule out injury as part of screening and routine monitoring. With next generation sequencing rapidly becoming the standard of care in quantifying DNA, understanding this science in the context of transplantation is critical to ensure studies, outcomes and care is improved.

Analysis of complex DNA mixtures using massively parallel sequencing of SNPs with low minor allele frequencies

DNA mixtures from 3 or more contributors have proven difficult to analyze using the current state-of-the-art method of short-tandem repeat (STR) amplification followed by capillary electrophoresis (CE). Here we analyze samples from both laboratory-defined mixtures and complex multi-contributor touch samples using a single nucleotide polymorphism (SNP) panel comprised of 2311 low-minor-allele-frequency loci, combined with massively parallel sequencing (MPS). This approach demonstrates that as many as 10 people can be identified in touch samples using a threshold of -Log P(RMNE) of 6, and a detection rate of 18-94 % across 10 different materials using a threshold of -Log P(RMNE) of 2. Thirty-two false positives were observed in 100 touch samples.

Potential of Next-Generation Sequencing in Noninvasive Fetal Molecular Blood Group Genotyping

Hemolytic disease of the fetus and newborn and fetal and neonatal alloimmune thrombocytopenia are caused by maternal antibodies against fetal alloantigens on red blood cells or platelets that are inherited from the father. After transplacental transport to the fetal circulation, antibodies of the IgG class may cause severe fetal anemia or bleeding complications. The indication for noninvasive fetal blood group genotyping is given if a clinically relevant antibody is detected in a pregnant woman and if the father is heterozygous (or unknown) for the implicated blood group allele. This mini-review will focus on the advantages and current limitations of next-generation sequencing (NGS) for noninvasive diagnosis of fetal blood groups which is, in contrast to fetal aneuploidy screening, proposed only by some research groups. Targeted massively parallel sequencing of short DNA fragments from maternal cell-free plasma samples enables counting of fetal alleles for many single nucleotide polymorphisms in parallel. This information can be utilized for estimation of the fetal fraction of cell-free DNA (cfDNA) as well as detection of the paternal blood group allele in question. Adherence to a cut-off of ≥4% fetal fraction for reporting conclusive results is recommended to avoid false-negative results due to low fetal fraction. For screening purposes of fetal RHD in RhD-negative pregnant women, real-time PCR methods are very well established. However, for diagnostic purposes, the targeted amplicon-based NGS approach has the inherent capability to estimate the fetal fraction of cfDNA. In the future, improving the accuracy of NGS by consensus sequencing of single cfDNA molecules may enable reliable fetal blood group genotyping already in the first trimester of pregnancy.

A new strategy to confirm the identity of tumour tissues using single-nucleotide polymorphisms and next-generation sequencing

With growing cancer morbidity, forensics cases in which archived tumour tissues can be used as biological samples are increasing, and an effective method to identify the body source of tumour tissues is needed. Single nucleotide polymorphisms (SNPs) may be a promising biomarker to identify the source of tumour tissues because of their low mutation rate and small amplicon size. Next-generation sequencing techniques offers the ability to detect hundreds of SNPs in a single run. The Precision ID Identity Panel (Thermo Fisher Scientific, Waltham, MA, USA) detects 90 autosomal SNPs for individual identification and 34 lineage-informative SNPs on Y chromosome using the Ion PGM system (Thermo Fisher Scientific). In this study, we evaluated performance of the panel for individual identification of tumour tissues. One hundred and fifty pairs of tumour tissues and corresponding normal tissues were analysed. Loss of heterozygosity was detected only in tumour tissues. The identity-by-state (IBS) scoring system was adopted to identify the body source of tumour tissues. The IBS score, as well as the number of loci with 2 alleles (A₂), 1 allele (A₁) and 0 alleles (A₀) shared, were analysed within each tumour-normal pair, unrelated individual pairs, parent-offspring pairs and full-sibling pairs. According to the probability distribution, threshold of A₂ in the range of 69 to 89 could achieve accuracy > 99% in identifying the source of tumour tissues. Thus, we developed a new strategy (process and criteria) to identify the source of tumour tissues that could be used in practice.

Estimation of the number of contributors to mixed samples of DNA by mitochondrial DNA analyses using massively parallel sequencing

We evaluated whether the number of contributors to mixed DNA samples can be estimated by analyzing the D-loop of mitochondrial DNA using massively parallel sequencing. The A- (positions 16,209-16,400) and B- (positions 30-284) amplicons in hypervariable regions 1 and 2, respectively, were sequenced using MiSeq with 2 × 251 cycles. Sequence extraction and trimming were performed using CLC Genomics Workbench 11 and the number of observed haplotypes was counted for each amplicon type using Microsoft Excel. The haplotype ratios were calculated by dividing the number of counted reads of the corresponding haplotype by the total number of sequence reads. Haplotypes that were over the threshold (5%) were defined as positive haplotypes. The number of larger positive haplotypes in either of the two amplicon types was defined as the number of contributors. Samples were collected from seven individuals. Seventeen mixed samples were prepared by mixing DNA from two to five contributors at various ratios. The number of contributors was correctly estimated from almost all of the mixed samples containing equal amounts of DNA from two to five people. In mixed samples of two or three people, the minor components were detected down to a ratio of 20:1 or 8:2:1. However, heteroplasmy, base deletions, and sharing of the same haplotypes caused incorrect estimations of the number of contributors. Although this method still has room for improvement, it may be useful for estimating the number of contributors in a mixed sample, as it does not rely on forensic mathematics.