- Research article
- Open Access
A SNP panel for identification of DNA and RNA specimens
BMC Genomics volume 19, Article number: 90 (2018)
SNP panels that uniquely identify an individual are useful for genetic and forensic research. Previously recommended SNP panels are based on DNA profiles and mostly contain intragenic SNPs. With the increasing interest in RNA expression profiles, we aimed for establishing a SNP panel for both DNA and RNA-based genotyping.
To determine a small set of SNPs with maximally discriminative power, genotype calls were obtained from DNA and blood-derived RNA sequencing data belonging to healthy, geographically dispersed, Dutch individuals. SNPs were selected based on different criteria like genotype call rate, minor allele frequency, Hardy–Weinberg equilibrium and linkage disequilibrium. A panel of 50 SNPs was sufficient to identify an individual uniquely: the probability of identity was 6.9 × 10− 20 when assuming no family relations and 1.2 × 10− 10 when accounting for the presence of full sibs. The ability of the SNP panel to uniquely identify individuals on DNA and RNA level was validated in an independent population dataset. The panel is applicable to individuals from European descent, with slightly lower power in non-Europeans. Whereas most of the genes containing the 50 SNPs are expressed in various tissues, our SNP panel needs optimization for other tissues than blood.
This first DNA/RNA SNP panel will be useful to identify sample mix-ups in biomedical research and for assigning DNA and RNA stains in crime scenes to unique individuals.
Over the past years, DNA profiles have found increasing use in the identification of human samples: they are ideal for sample tracking in biomedical studies and forensic investigations. In recent years, joint analysis of DNA and RNA has proven to be valuable: 1) Forensic investigations where RNA profiles may complement DNA profiles, and may be used to establish the tissue origin of the specimen , wound age, post-mortem interval, and the age of stains [2,3,4,5]; 2) Population research in which the genetic component of gene expression is studied.
Single nucleotide polymorphisms (SNPs) and other genetic markers like mitochondrial haplotypes, Y chromosomal markers and short tandem repeats (STRs) are all used for individual identification. Mitochondrial DNA (mtDNA) is found in both females and males but it is inherited only through the mother, which makes it impossible to differentiate between mothers and offspring [1, 6]. Y chromosome, a male-specific part is widely used in genetics studies and forensic data analysis particularly in cases of sexual assault, although identification of individuals using Y chromosome DNA analysis is limited to non-related subjects [6,7,8]. STRs, regions with short repeat units (usually 2–6 base pairs in length), are highly informative because of the large number of alleles present even in genetically rather homogeneous populations [9,10,11]. Despite their high discriminatory power, they have some limitations such as required large amplicon sizes, high mutation rates, and the presence of artefacts, which can negatively influence the downstream analysis [12, 13]. To overcome these limitations, SNPs have been more recently introduced for individual identification .
SNPs are defined here as single nucleotide substitutions that occur in more than 1 % of the general population. SNP assays can be used for multiple types of genetics studies. A recent review by Kayser and de Knijff provides an overview of recent advances in the use of SNPs for forensic investigation . Many studies have discussed the advantages of SNPs compared to STRs, including low mutation rates, fast genotyping, high abundance in the genome, and straightforward detection using high-throughput technologies [16,17,18,19,20]. Kidd et al.  described a set of SNPs with high heterozygosity and low frequency variation in different populations, both helpful characteristics for an individual identification panel. In the last decade, several research groups have found valuable alternative individual identification SNP (IISNP) panels for different populations in the world [14, 22,23,24,25,26]. Pakstis et al.  selected 86 unlinked candidate IISNPs for 44 major populations across the world. Also, recently an IISNP panel for global forensic casework was established (Illumina, 2015). Moreover, research groups have tried to develop new SNP markers to improve their discrimination power using high throughput data sets [13, 27, 28]. However, these panels mostly consist of intragenic SNPs or SNPs in genes that are not expressed in blood or other tissues relevant for forensic identification. Therefore, there is no suitable panel that contains a minimum of informative SNPs that can be used on both DNA and RNA specimens.
To address the limitations of current IISNP panels, we present a small and powerful IISNP panel. When turned into a multiplex assay, this panel can be exploited for unequivocal identification of individuals in both DNA and RNA specimens in forensic investigations. Moreover the panel can be used to identify sample mix-ups in human gene expression studies, which have been demonstrated to severely affect the power of such studies [PMID: 21653519].
Sample collection and sequencing
Several Dutch biobanks contributed to sample collection of Dutch ancestry (with parents born in the Netherlands) within the Biobanking and Biomolecular Research Infrastructure-Netherlands (BBMRI-NL), established as a national node of the European BBMRI infrastructure in the Netherlands . Our DNA datasets are derived from Genome of the Netherlands (GoNL), a whole-genome-sequencing effort within BBMRI-NL consisting of 250 representative parent-offspring families widely dispersed across the Netherlands (231 trios and 19 quartets, of which 11 had monozygotic twins and 8 had dizygotic twins), which aims to characterize DNA variation in the Dutch population. DNA-based genotype calls were derived from DNA isolated from blood .
The Biobank-based Integrative Omics Studies Consortium (BIOS Consortium) is also part of BBMRI-NL. Its aim is to create a large-scale data infrastructure and to bring together BBMRI researchers focusing on integrative omics studies in Dutch biobanks. The BIOS Consortium applies a functional genomics approach that integrates genome-wide genetic data with data on the epigenome and transcriptome. The RNA data were derived from individuals from seven Dutch biobanks participating in the BIOS Consortium (LL, LifeLines Cohort Study; LLS, Leiden Longevity Study; RS, Rotterdam Study; CODAM, Cohort on Diabetes and Atherosclerosis Maastricht; NTR, the Dutch Twin Registry; PAN, Prospective ALS study Netherlands). Globin RNAs were removed from whole blood and the polyA fraction of the remaining RNA was subjected to RNA sequencing using HiSeq2000 sequencers and were analyzed as described by . RNA-Seq data are available in the European Genome-Phenome Archive (EGA) under accession: EGAD00001001623. Additionally, for each sample microarray-derived SNP data (Immunochip on all samples and at least one other GWAS array per sample) were generated by the biobanks [30,31,32,33]. Sequencing and the primary analysis of the data was performed within the BIOS and GoNL working groups. Variant calling was done using Samtools (v1.1) and Varscan (v2.3.7) and genotype calling then was performed at the SNP sites. We have carried out further analysis based on RNA-based genotypes calling on 2115 samples (after removing related samples) (Additional file 1: Table S1).
Filtering and identifying SNPs
Our strategy and criteria are shown in Fig. 1. Briefly, the allele frequencies for each SNP were calculated by genotype counting within population assuming each marker is a two-allele, co-dominant system. Genotype frequencies and maximum/minimum allele frequencies were then calculated based on allele frequencies. Reference alleles were extracted from Ensembl GRCh37 release 84 and alternative alleles were acquired for each SNP. Alternative allele frequencies (AAF) were measured and, minor allele frequency (MAF) was calculated based on the frequency of least common allele for each SNP in population. Then we have eliminated: 1) SNPs with genotype call rate less than 90%, 2) SNPs with more than three genotypes or more than two alleles, 3) SNPs with MAF less than 0.2 in both DNA-Seq and RNA-Seq data. To further filtration, the Hardy–Weinberg equilibrium (HWE), and linkage disequilibrium (LD) tests were used. Also, we have ignored SNPs located on human leukocyte antigen (HLA) loci where the SNP calls are prone to artefacts, particularly in NGS-derived datasets (Fig. 1). Finally, consistency was determined by analyzing the AAF in the RNA and DNA data from the same set of 2115 individuals (Additional file 1: Table S1), and 50 independent SNPs were selected with high MAF, high heterozygosity and no LD between them in our population.
The selected SNPs were compared with previously published panels such as SNPforID 52-plex, 75 Chinese SNPs, 30 Korean SNPs and 92 IISNPs [20, 23, 34, 35]. Furthermore, the AAF of the 50 SNPs were compared to 1000 Genomes, Exome Aggregation Consortium (ExAC), NHLBI GO Exome Sequencing Project (Go.ESP) and the 1000 Genomes phase-3 populations. Moreover, the predicted effect of our variants on protein sequence was extracted from Ensembl. R scripts implementing all these filtering steps.
Observed heterozygosity and expected heterozygosity were estimated based on genotype counts and allele frequencies, respectively. So, deviation from HWE was determined by comparing the expected and observed number of individuals with each possible genotype using Fisher exact test. The p-values of the HW tests were corrected for multiple testing using the method developed by Benjamin and Hochberg , implementing a false discovery rate (FDR) of 1%.
To evaluate statistical independence of the SNPs, r2 was calculated for all unique pairwise combinations of the common SNPs between DNA-Seq and RNA-Seq data. The LD values were used to determine whether there was any evidence for significant linkage among the markers. In addition, heterozygosity, fixation index (also called the inbreeding coefficient, is defined as (He – Ho) / He (where He is expected heterozygosity and Ho is observed heterozygosity). It may range from − 1 to + 1. Values close to zero are expected under random mating, while substantial positive values indicate inbreeding or undetected null alleles. Negative values indicate excess of heterozygosity, due to negative assortative mating, or selection for heterozygotes.), polymorphism, Hardy–Weinberg equilibrium and probability of identity (provides an estimate of the average probability that two unrelated individuals, drawn from the same randomly mating population, will by chance have the same multilocus genotype. It is also called Population Match Probability), were checked using the Excel add-in GenAlEx [37, 38].
The performance of the 50 SNP panel was examined by two strategies. First, the ability of the 50 SNP panel to identify sample mix-ups was checked in the second batch of samples from the BIOS project with 1357 independent samples, not used for the panel identification, using the same genotype calling methods as described in the first phase of the project. In the second strategy, the SNP panel was compared with a much larger set of 2622 SNPs. The 2622 SNP panel was selected as follows: 1) All ExAC v0.3 biallelic SNPs were selected ; 2) overlapped with exonic regions from Ensembl v75 and 1000 Genomes phase 1 high confidence SNPs, and 3) filtered on MAF > 0.4. Simultaneously, 1) SNPs that were called using Unified Genotyper in the LL subset of GoNL RNASeq samples were selected, and 2) filtered on a Unified Genotyper quality score of > 100,000. The final set of 2622 SNPs was constructed by overlapping the ExAC based and RNA-Seq based SNP lists. Both 2622 SNP and 50 SNP panels were used to evaluate the ability for sample and mismatch identification in samples from the Diabetes Care System (DCS) cohort , where there were 562 RNA-Seq as well as 3428 GWAS samples imputed using the HapMap reference panel. Genotype calling at all genomics coordinates was performed for the two SNP panels in both RNA-Seq and imputed GWAS samples. Briefly, our RNA-Seq pipeline (http://biopet-docs.readthedocs.io/en/latest/pipelines/gentrap/) was used to obtain genotype calls for 562 RNA-Seq samples. Further, for imputed GWAS samples, genotyping using the Human Core Exome chip was performed according the manufacturers protocol (Illumina Inc. San Diego, Ca, USA). Then, quality control was demonstrated using following settings: a cut-off of 99% for genotyping call rate, gentrain and cluster score < 0.6 and 0.4, respectively, and p-value cut-off for Hardy-Weinberg equilibrium was set at 10− 4. Consequently, imputation was done using SHAPEIT (v2.r644) and IMPUTE (v2.3.2). These two files were merged into one multi-sample VCF file that contained genotype information of total 3990 samples. All genomics coordinates specified by those two SNP lists were included in this VCF file. To compare these two SNP panels, we first calculated the pair-wise allele concordance scores for all 3990 by 3990 sample pairs by examining only the genomics coordinates specified in the 50 SNP panel. The allele concordance score (ranging from 0 to 1) for each sample pair was defined as the ratio between the number of identical alleles and the total number of alleles (100 alleles in total) through all 50 SNPs coordinates. In addition, the identification of the best matching GWAS samples for each RNA-Seq sample was examined with a minimal allele concordance score of 0.8. Multiple best GWAS sample hits for one RNA-Seq sample was possible as there have been repeated GWAS measurements performed on the same person. Then we performed the same steps using the long 2622 SNP panel to identify the best matching GWAS samples for each RNA-Seq sample.
Data and SNP identification
In previous studies, most of the recommended individual identification SNP panels were generated based on DNA profiles and they mostly contain intragenic SNPs. Therefore, there is no efficient panel with informative SNPs which can be used for both DNA and RNA. To overcome these limitations, in this study both DNA and RNA-based genotype calls were used to find a small set of SNPs that can be used for identification of individuals. DNA-based genotype calls (19,562,004 SNP positions) and RNA-based genotype calls were made on 2115 individuals from four Dutch biobanks. Reliable RNA-based genotype calls were obtained for 507,975 SNP positions across four cohorts (Additional file 1: Table S1).
To find the most informative SNPs, we applied a number of filtering steps: details on the different filters applied and number of SNPs remaining can be found in the Method section, Fig. 1 and Table 1. After selection of SNPs with reliable genotype calls and high discriminative power (i.e. high MAF) and SNPs in Hardy–Weinberg equilibrium, 100 SNPs remained (Table 1). A final step to select the smallest informative set of SNPs for individual identification is absence of linkage disequilibrium between the SNP positions. To evaluate statistical independence of the SNPs, r2 was calculated for all unique pairwise combinations of 100 SNPs (Fig. 2a). SNP positions with r2 less than 0.01 were selected (Fig. 2b). We removed the least heterozygous SNP of any two SNPs, with LD and close genetic distance based on criteria of Kidd et al.  and Hwa et al. . In addition, SNPs located in the HLA region were removed. Consequently, 50 SNPs were selected as a final panel to identify Dutch individuals (Table 1; Additional file 2: Table S2). The r2 values for the final list of 50 SNPs are very close to zero (median: 2.4 × 10− 4; average: 5.5 × 10− 4).
Variant consequences for each position were extracted from Ensembl. As expected, most SNPs are located in exonic regions, because we selected SNPs which can be consistently detected in both RNA and DNA data. Twenty-four SNPs are located in the 3′-UTR, 22 SNPs in the coding region (14 synonymous, 8 missense), 3 SNPs in the 5′-UTR and 1 SNP in an intron (Additional file 2: Table S2).
Characteristics and quality control of the 50 SNP set
The final selection of SNPs showed highly concordant AAF, for DNA and RNA-derived genotypes suggestive of the absence of bias in the genotype calls of these SNPs (Fig. 3a, r2 = 0.98). The correlation of MAF between RNA and DNA genotype calls was much higher for 50 selected SNPs (r2 = 0.81) than all confidently called SNP genotypes in both datasets (r2 = 0.6) (Additional file 3: Figure S1). The average MAF in the 50 selected SNPs was 0.35. The final selection of SNPs also had AAF > 0.15 in 1000 Genomes, but the AAF was generally higher in the Dutch population (Fig. 3b).
Observed heterozygosity and expected heterozygosity were measured based on genotype counts and allele frequencies, respectively. There was a high positive correlation (r2 = 0.96) with nearly equal frequency between expected heterozygosity and observed heterozygosity of 50 selected SNPs, suggesting that there was no large bias in detection of these SNPs. (Additional file 4: Figure S2 and Additional file 5: Figure S3).
Population genetic parameters were calculated to further characterize the SNPs in the panel (details in the Methods section) using the Excel add-in GenAlEx [37, 38]. The fixation indices were slightly negative for most SNPs (average: − 0.019), indicating some negative assortative mating with proper heterozygosity (Additional file 6: Figure S4 and Additional file 7: Figure S5).
The probability of identity was analyzed for the panel of 50 SNPs in 2115 Dutch samples. The PI provides an estimate of the average probability of two independent samples having identical genotype calls. It is used to determine the minimum number of SNPs which are needed for identity calling. PI can be calculated under the assumption that all individuals are unrelated or under that the assumption that individuals may be related (PI-sibs). Figure 4 shows that at least 17 SNPs are required to achieve uniqueness in 2115 Dutch samples (PI is 3.3 × 10− 7 for unrelated individuals and PI-sibs is 4.4 × 10− 4). For our final list of 50 SNPs, the PI was 6.9 × 10− 20 and PI-sibs were 1.2 × 10− 10 (Fig. 4). This makes the marker set appropriate for tagging and tracking samples in large biomedical, association, and epidemiological studies.
To investigate the utility of our 50 SNP panel in other than the Dutch population, AAF of 50 SNPs were compared with the AAF in 1000 Genomes, ExAC and Go.ESP (Fig. 5a; Additional file 4: Figure S2). Figure 5a shows the AAF in Dutch SNPs are mostly consistent with those from ExAC. They are overall slightly higher in the Dutch population than in the populations catalogued in these three other databases, but consistently high for these databases which contain mostly individuals from European descent (Fig. 5a).
Even in non-European populations, such as the South Asian and American populations, most of the 50 SNPs had MAF > 0.2 (AAF between 0.2 and 0.8) (Fig. 5b; Additional file 2: Table S2). However, for the African population, our SNP panel was less effective, because 18 out of the 50 SNPs had MAF < 0.2 (Fig. 5b; Additional file 2: Table S2).
The selected SNPs were compared with previous published panels such as SNPforID 52-plex, 75 Chinese SNPs, 30 Korean SNPs and 92 IISNPs [20, 23, 34, 35]. There were no common SNPs between these panels and our 50 SNP panel.
Validation of the 50 SNP panel
The ability of the SNP to uniquely identify individuals in both RNA and DNA level was evaluated by studying genotype concordance in an independent, paired set of 1357 DNA and blood-derived RNA samples. The distribution of concordant genotype calls for matched and unmatched DNA and RNA samples was clearly distinct and non-overlapping, with all matching samples having concordant genotype calls for at least 38 out of 50 of the SNPs (Fig. 6). Forty out of 50 SNPs demonstrated more than 90% concordant DNA and RNA genotype calls in this validation set of 1357 samples, whereas the minimum concordance observed was 65% for SNP rs2230267.
The ability of the 50 SNP panel to identify DNA and RNA sample concordance was compared with a bigger panel of 2622 SNPs in an independent study. The list of best GWAS hits for all 538 RNA-Seq samples with the RNA-Seq to GWAS mapping list provided by the study coordination center showed that the shorter 50 SNP panel detected more matching samples, while still reporting potential sample mix-ups (Table 2). Specifically, the number of samples where the genotype concordance test was indecisive was larger for the 2622 SNP panel.
We have established a pipeline to select a Dutch-specific SNP panel based on both DNA-Seq and RNA-Seq data. This panel consists of 50 SNPs with high heterozygosity, high MAF, low LD and robust detection in blood DNA and RNA (Table 1; Additional file 2: Table S2).
During the past years, various SNP panels have been published for individual identification [20,21,22,23,24,25, 41,42,43,44,45,46,47]. The SNPforID consortium developed 52 SNPs for individual identification . Also, Kim et al.  developed a SNP-based individual assignment system containing 30 SNP loci for Korean individuals. Lou et al.  reported the performance of a 44 SNPs individual identification assay for Chinese. In addition, studies indicated considerable potential of high throughput platforms for SNP detection which could increase unprecedented discriminative power for human identification [49,50,51]. These panels are all based on DNA profiles and mostly contain intragenic SNPs which disqualify them for RNA-based genotype calling. Our set of SNPs is 98% exonic and can uniquely identify individuals in DNA and RNA profiles. Although use of coding SNPs in forensic DNA analysis may be restricted due to specific legislation in certain countries, this should not apply to highly polymorphic SNPs without any associations with appearance phenotypes.
There is no overlap between Dutch selected IISNPs and the established IISNP panels [13, 20, 23, 35]. One of our selected SNPs (rs1866141) is located in an intronic region of the highly expressed GNLY gene (MIM # 188855). Observing intronic SNPs among RNA-based calls analysis is common because of the presence of pre-mRNA, incomplete splicing or intron retention.
The 50 SNP panel showed better identification performance when it was compared with another panel containing 2622 SNPs (Table 2). Although the 50 SNP panel had superior performance in the Dutch and other populations of European descent, it is less optimal for individuals from African populations (Fig. 5b). There was no information for two SNPs (rs1950943 and rs1042558) in non-European populations.
RNA profiles can complement DNA profiles in research projects and forensic applications. In forensics, despite the presumed low stability of RNA, RNA profiles could not only identify individuals but also provide information about the type of tissues found at the crime scene, wound age determination, determination of the post-mortem interval and the functional status of cells as well as organs [2,3,4,5, 41, 52,53,54,55,56,57,58,59]. Unlike previous studies [20, 21, 25], 98% of our final SNPs located in exonic regions. For this reason, a SNP profiling assay for this 50 SNP panel can be an efficient method for individual identification in RNA and DNA stains from crime scenes. Analysis of RNA from crime scenes has demonstrated differences in RNA degradation rates. Although the SNPs from the 50 SNP panel are in high expressed genes and are robustly detected in biobank samples, their robust detection in forensic specimens containing partially to severely degraded RNA still needs to be demonstrated.
To address whether the SNP panel could also be used on RNA samples from tissues other than blood, the expression level of genes in which the SNPs are located was surveyed using GTEx portal . While 25 genes were expressed ubiquitously and 17 genes were expressed in multiple other tissues, the expression of 8 genes (MNDA (MIM # 159553), SELL (MIM # 153240), CSF3R (MIM # 138971), IFITM2 (MIM # 605578), FPR1 (MIM # 136537), CXCR2 (MIM # 146928), GNLY (MIM # 188855), FCN1 (MIM # 601252)) was rather specific for whole blood, as their expression levels were near to zero in other tissues. Further, IFI30 (MIM # 604664), contained the lowest gene expression in different tissues. As 17 SNPs are sufficient to identify an individual uniquely in the Dutch population it is assumed that our SNP panel could be effective for identification of individuals in RNA from different tissues than blood. However, when testing our panel on another SNP-chip genotype and RNA-Seq dataset, consisting of 36 samples from brain tissue, it appeared that many of the 50 SNPs had insufficient coverage for reliable RNA-based genotype calling in brain. This indicates that our 50 SNP panel needs optimization to be used for other tissues than blood and/or high sequencing depth.
Our SNP panel compared favorably to previously published panels in terms of discrimination power (probability of identity (PI)), even for closely related individuals. The FBI (USA) has selected 13 STR loci to serve as a panel for forensic investigations (CODIS, Combined DNA Index System). With this set of loci, the probability of a match between the profiles of two unrelated persons in a randomly mating population of Caucasian Americans PI is 2.97 × 10− 15 . Also, the 52-plex SNP assay, which is now more routinely used in Europe, has a mean PI of 5.0 × 10− 21 for the European population . The PI of our SNP panel was 6.9 × 10− 20 and 1.2 × 10− 10 in unrelated and related individuals, respectively, similar to the 52-plex SNP assay, far more discriminative than the 13 CODIS markers  and similar to the recently published study of 20 CODIS markers in Caucasian Americans [61,62,63] .
We developed a first SNP panel based on both DNA and RNA data for the Dutch population. This panel contains 50 informative SNPs with high heterozygosity, low PI and close MAF and AAF frequencies in DNA and RNA. It will be useful for efficient sample identification/tagging in large biomedical, association, and epidemiologic studies, and for developing forensic profiling and kinship assays. Our panel will be useful for other European populations and can be considered in conjunction with other panels to develop a global IISNP panel with more markers.
Alternative Allele Frequencies
Biobanking and Biomolecular Research Infrastructure-Netherlands
- BIOS Consortium:
The Biobank-based Integrative Omics Studies Consortium
Cohort on Diabetes and Atherosclerosis Maastricht
Combined DNA Index System
Diabetes Care System
European Genome-Phenome Archive
Exome Aggregation Consortium
False Discovery Rate
NHLBI GO Exome Sequencing Project
Genome of the Netherlands
Human Leukocyte Antigen
Individual Identification SNP
Life Lines Cohort Study
Leiden Longevity Study
Minor Allele Frequency
the Dutch Twin Registry
Prospective ALS study Netherlands
Probability of Identity
Single Nucleotide Polymorphisms
Short Tandem Repeats
Phillips ML. Crime scene genetics: transforming forensic science through molecular technologies. Bioscience. 2008;58(6):484–9.
Bauer M. RNA in forensic science. Forensic Sci Int Genet. 2007;1(1):69–74.
Sijen T. Molecular approaches for forensic cell type identification: on mRNA, miRNA, DNA methylation and microbial markers. Forensic Sci Int Genet. 2015;18:21–32.
den Berge M, Sijen T. Extended specificity studies of mRNA assays used to infer human organ tissues and body fluids. Electrophoresis. 2017;38:3155–60.
Zhao H, Wang C, Yao L, Lin Q, Xu X, Hu L, et al. Identification of aged bloodstains through mRNA profiling: experiments results on selected markers of 30- and 50-year-old samples. Forensic Sci Int. 2017;272(Supplement C):e1–6.
Ambers AD, Churchill JD, King JL, Stoljarova M, Gill-King H, Assidi M, et al. Erratum to: more comprehensive forensic genetic marker analyses for accurate human remains identification using massively parallel DNA sequencing. BMC Genomics. 2017;18(1):312.
Kayser M. Forensic use of Y-chromosome DNA: a general overview. Hum Genet. 2017;136(5):621–35.
Jobling MA, Tyler-Smith C. Human Y-chromosome variation in the genome-sequencing era. Nat Rev Genet. 2017;18(8):485–97.
Brinkmann B, Rand S, Bajanowski T. Forensic identification of urine samples. Int J Legal Med. 1992;105(1):59–61.
Budowle B, Moretti TR, Baumstark AL, Defenbaugh DA, Keys KM. Population data on the thirteen CODIS core short tandem repeat loci in African Americans, US Caucasians, Hispanics, Bahamians, Jamaicans, and Trinidadians. J Forensic Sci. 1999;44(6):1277–86.
Urquhart A, Kimpton C, Downes T, Gill P. Variation in short tandem repeat sequences—a survey of twelve microsatellite loci for use as forensic identification markers. Int J Legal Med. 1994;107(1):13–20.
Butler JM, Shen Y, McCord BR. The development of reduced size STR amplicons as tools for analysis of degraded DNA. J Forensic Sci. 2003;48(5):1054–64.
Kim SM, Yoo SY, Nam SH, Lee JM, Chung KW. Identification of Korean-specific SNP markers from whole-exome sequencing data. Int J Legal Med. 2016;130(3):669–77.
Dixon LA, Murray CM, Archer EJ, Dobbins AE, Koumi P, Gill P. Validation of a 21-locus autosomal SNP multiplex for forensic identification purposes. Forensic Sci Int. 2005;154(1):62–77.
Kayser M, de Knijff P. Improving human forensics through advances in genetics, genomics and molecular biology. Nat Rev Genet. 2011;12(3):179–92.
Amorim A, Pereira L. Pros and cons in the use of SNPs in forensic kinship investigation: a comparative analysis with STRs. Forensic Sci Int. 2005;150(1):17–21.
Borsting C, Sanchez JJ, Hansen HE, Hansen AJ, Bruun HQ, Morling N. Performance of the SNPforID 52 SNP-plex assay in paternity testing. Forensic Sci Int Genet. 2008;2(4):292–300.
Kidd KK, Kidd JR, Speed WC, Fang R, Furtado MR, Hyland FC, et al. Expanding data and resources for forensic use of SNPs in individual identification. Forensic Sci Int Genet. 2012;6(5):646–52.
Krjutskov K, Viltrop T, Palta P, Metspalu E, Tamm E, Suvi S, et al. Evaluation of the 124-plex SNP typing microarray for forensic testing. Forensic Sci Int Genet. 2009;4(1):43–8.
Pakstis AJ, Speed WC, Fang R, Hyland FC, Furtado MR, Kidd JR, et al. SNPs for a universal individual identification panel. Hum Genet. 2010;127(3):315–24.
Kidd KK, Pakstis AJ, Speed WC, Grigorenko EL, Kajuna SL, Karoma NJ, et al. Developing a SNP panel for forensic identification of individuals. Forensic Sci Int. 2006;164(1):20–32.
Hou G, Jiang X, Yang Y, Jia F, Li Q, Zhao J, et al. A 21-locus autosomal SNP multiplex and its application in forensic science. J Forensic Sci. 2014;59(1):5–14.
Hwa H-L, Wu LSH, Lin C-Y, Huang T-Y, Yin H-I, Tseng L-H, et al. Genotyping of 75 SNPs using arrays for individual identification in five population groups. Int J Legal Med. 2016;130(1):81–9.
Ibarra A, Freire-Aradas A, Martínez M, Fondevila M, Burgos G, Camacho M, et al. Comparison of the genetic background of different Colombian populations using the SNPforID 52plex identification panel. Int J Legal Med. 2014;128(1):19–25.
Pakstis AJ, Speed WC, Kidd JR, Kidd KK. Candidate SNPs for a universal individual identification panel. Hum Genet. 2007;121(3–4):305–17.
Wei Y-L, Qin C-J, Liu H-B, Jia J, Hu L, Li C-X. Validation of 58 autosomal individual identification SNPs in three Chinese populations. Croat Med J. 2014;55(1):10–3.
Pakstis AJ, Haigh E, Cherni L, ElGaaied AB, Barton A, Evsanaa B, et al. 52 additional reference population samples for the 55 AISNP panel. Forensic Sci Int Genet. 2015;19:269–71.
Seo SB, King JL, Warshauer DH, Davis CP, Ge J, Budowle B. Single nucleotide polymorphism typing with massively parallel sequencing for human identification. Int J Legal Med. 2013;127(6):1079–86.
Genome of the Netherlands C. Whole-genome sequence variation, population structure and demographic history of the Dutch population. Nat Genet. 2014;46(8):818–25.
Zhernakova DV, Deelen P, Vermaat M, van Iterson M, van Galen M, Arindrarto W, et al. Identification of context-dependent expression quantitative trait loci in whole blood. Nat Genet. 2017;49(1):139–45.
Scholtens S, Smidt N, Swertz MA, Bakker SJ, Dotinga A, Vonk JM, et al. Cohort profile: LifeLines, a three-generation cohort study and biobank. Int J Epidemiol. 2015;44(4):1172–80.
Tigchelaar EF, Zhernakova A, Dekens JA, Hermes G, Baranska A, Mujagic Z, et al. Cohort profile: LifeLines DEEP, a prospective, general population cohort study in the northern Netherlands: study design and baseline characteristics. BMJ Open. 2015;5(8):e006772.
van Dam RM, Boer JM, Feskens EJ, Seidell JC. Parental history of diabetes modifies the association between abdominal adiposity and hyperglycemia. Diabetes Care. 2001;24(8):1454–9.
Kim J-J, Han B-G, Lee H-I, Yoo H-W, Lee J-K. Development of SNP-based human identification system. Int J Legal Med. 2010;124(2):125–31.
Sanchez JJ, Phillips C, Borsting C, Balogh K, Bogus M, Fondevila M, et al. A multiplex assay with 52 single nucleotide polymorphisms for human identification. Electrophoresis. 2006;27(9):1713–24.
Benjamini Y, Hochberg Y. Controlling the false discovery rate: a practical and powerful approach to multiple testing. J R Stat Soc Ser B Methodol. 1995;57(1):289–300.
Peakall R, Smouse PE. GenAlEx 6.5: genetic analysis in excel. Population genetic software for teaching and research--an update. Bioinformatics. 2012;28(19):2537–9.
Peakall ROD, Smouse PE. genalex 6: genetic analysis in excel. Population genetic software for teaching and research. Mol Ecol Notes. 2006;6(1):288–95.
Karczewski KJ, Weisburd B, Thomas B, Solomonson M, Ruderfer DM, Kavanagh D, et al. The ExAC browser: displaying reference data information from over 60 000 exomes. Nucleic Acids Res. 2017;45(D1):D840–5.
van der Heijden AA, Rauh SP, Dekker JM, Beulens JW, Elders P, t Hart LM, et al. The Hoorn Diabetes Care System (DCS) cohort. A prospective cohort of persons with type 2 diabetes treated in primary care in the Netherlands. BMJ Open. 2017;7(5):e015599.
Amigo J, Salas A, Phillips C, Carracedo A. SPSmart: adapting population based SNP genotype databases for fast and comprehensive web access. BMC Bioinformatics. 2008;9:428.
Lee HY, Park MJ, Yoo JE, Chung U, Han GR, Shin KJ. Selection of twenty-four highly informative SNP markers for human identification and paternity analysis in Koreans. Forensic Sci Int. 2005;148(2–3):107–12.
Zeng Z, Wang L, Feng Q, Zhang L, Lee L, Wang L, et al. Evaluation of 96 SNPs in 14 populations for worldwide individual identification. J Forensic Sci. 2012;57(4):1031–5.
Zeng Z, Yan H, Wang L, Yuan E, Yang W, Liao Z, et al. Genome-wide screen for individual identification SNPs (IISNPs) and the confirmation of six of them in Han Chinese with pyrosequencing technology. J Forensic Sci. 2010;55(4):901–7.
Sarkar A, Nandineni MR. Development of a SNP-based panel for human identification for Indian populations. Forensic Sci Int Genet. 2017;27:58–66.
Zhang S, Bian Y, Chen A, Zheng H, Gao Y, Hou Y, et al. Developmental validation of a custom panel including 273 SNPs for forensic application using Ion Torrent PGM. Forensic Sci Int Genet. 2017;27:50–7.
Giardina E, Pietrangeli I, Martone C, Asili P, Predazzi I, Marsala P, et al. In silico and in vitro comparative analysis to select, validate and test SNPs for human identification. BMC Genomics. 2007;8(1):457.
Lou C, Cong B, Li S, Fu L, Zhang X, Feng T, et al. A SNaPshot assay for genotyping 44 individual identification single nucleotide polymorphisms. Electrophoresis. 2011;32(3–4):368–78.
Li L, Wang Y, Yang S, Xia M, Yang Y, Wang JP, et al. Genome-wide screening for highly discriminative SNPs for personal identification and their assessment in world populations. Forensic Sci Int Genet. 2017;28:118–27.
Zhang S, Bian Y, Chen A, Zheng H, Gao Y, Hou Y, et al. Massively parallel sequencing of 231 autosomal SNPs with a custom panel: a SNP typing assay developed for human identification with Ion Torrent PGM. Forensic Sci Res. 2017;2(1):26–33.
Zhang S, Bian Y, Zhang Z, Zheng H, Wang Z, Zha L, et al. Parallel analysis of 124 universal SNPs for human identification by targeted semiconductor sequencing. Sci Rep. 2015;5:18683.
Xu Y, Xie J, Cao Y, Zhou H, Ping Y, Chen L, et al. Development of highly sensitive and specific mRNA multiplex system (XCYR1) for forensic human body fluids and tissues identification. PLoS One. 2014;9(7):e100123.
Danaher P, White RL, Hanson EK, Ballantyne J. Facile semi-automated forensic body fluid identification by multiplex solution hybridization of NanoString® barcode probes to specific mRNA targets. Forensic Sci Int Genet. 2015;14(Supplement C):18–30.
Haas C, Hanson E, Anjos MJ, Ballantyne KN, Banemann R, Bhoelai B, et al. RNA/DNA co-analysis from human menstrual blood and vaginal secretion stains: results of a fourth and fifth collaborative EDNAP exercise. Forensic Sci Int Genet. 2014;8(1):203–12.
van den Berge M, Sijen T. A male and female RNA marker to infer sex in forensic analysis. Forensic Sci Int Genet. 2017;26:70–6.
van den Berge M, Wiskerke D, Gerretsen R, Tabak J, Sijen T. DNA and RNA profiling of excavated human remains with varying postmortem intervals. Int J Legal Med. 2016;130(6):1471–80.
Lech K, Ackermann K, Revell VL, Lao O, Skene DJ, Kayser M. Dissecting daily and circadian expression rhythms of clock-controlled genes in human blood. J Biol Rhythm. 2016;31(1):68–81.
Lech K, Liu F, Ackermann K, Revell VL, Lao O, Skene DJ, et al. Evaluation of mRNA markers for estimating blood deposition time: towards alibi testing from human forensic stains with rhythmic biomarkers. Forensic Sci Int Genet. 2016;21:119–25.
Zubakov D, Liu F, Kokmeijer I, Choi Y, van Meurs JBJ, van Ijcken WFJ, et al. Human age estimation from blood using mRNA, DNA methylation, DNA rearrangement, and telomere length. Forensic Sci Int Genet. 2016;24:33–43.
Consortium GT. The genotype-tissue expression (GTEx) project. Nat Genet. 2013;45(6):580–5.
Butler JM. Hill CR, Coble MD. Variability of new STR loci and kits in U.S. population groups. Profiles in DNA. 2012. https://worldwide.promega.com/resources/profiles-in-dna/2012/variability-of-new-str-loci-and-kits-in-us-population-groups/. Accessed 12 Dec 2017.
Hares DR. Expanding the CODIS Core Loci in the United States. Forensic Sci Int Genet. 2012;6:52–4.
Hares DR. Selection and implementation of expanded CODIS core loci in the United States. Forensic Sci Int Genet. 2015;17:33–4.
BIOS consortium: Jansen R, van Meurs JB, Heijmans BT, Boomsma DI, van Dongen J, Hottenga JJ, Slagboom PE, Suchiman HED, van Iterson M, van Zwet EW, ‘t Hoen PA, Pool R, van Greevenbroek MM, Stehouwer CD, van der Kallen CJ, Schalkwijk CG, Wijmenga C, Zhernakova A, Tigchelaar EF, Beekman M, Deelen J, van Heemst D, Veldink JH, van den Berg LH, van Duijn CM, Hofman BA, Uitterlinden AG, Jhamai PM, Verbiest M, Verkerk M, van der Breggen R, van Rooij J, Lakenberg N, Mei H, Bot J, Zhernakova DV, Van’t Hof P, Deelen P, Nooren I, Moed M, Vermaat M, Bonder MJ, van Dijk F, van Galen M, Arindrarto W, Kielbasa SM, Swertz MA, Isaacs A, Franke L.
This work was performed within the framework of the Biobank-Based Integrative Omics Studies (BIOS) Consortium and the GoNL Project which are funded by BBMRI-NL, a research infrastructure financed by the Netherlands Organization for Scientific Research (NWO project 184.021.007).
Availability of data and materials
All steps are presented in the main paper and additional files. Also, we have deposited data on the following repositories: Accession ID: EGAD00001001623, Databank URL: https://ega-archive.org/, Accession ID: EGAD00001000743, Databank URL: https://ega-archive.org/.
Ethics approval and consent to participate
The ethical approval for this study lies with the individual participating cohorts (CODAM, LL, LLS and RS). A broad consent for participation in research, including research on genotypes, was obtained from all participants. Given the privacy-sensitive nature of the DNA and RNA data, the data has been deposited at European Genome-Phenome Archive (EGA) and is under controlled access. Requests for the data can be filed in the EGA system and will be handled by the BIOS data access committee. The committee will provide access to researchers for studies with a solid scientific background.
Consent for publication
A consent for publication of the results from research projects, including research on genotypes, lies with the individual participating cohorts (CODAM, LL, LLS and RS).
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Four biobanks (RNA-Seq) with number of different SNP positions and number of samples. (DOC 33 kb)
Characteristics of 50 selected SNPs. (CSV 15 kb)
Distribution of MAF calculated from DNA and RNA data. (DOC 58 kb)
Plot of the expected heterozygosity (x-axis) and observed heterozygosity (y-axis) for the 50 SNPs in the panel. Pearson correlation is 0.98. (DOC 54 kb)
Individual heterozygosity across loci for each sample determined by GenAlEx software. (DOC 46 kb)
The fixation index values of 50 selected SNPs. (DOC 54 kb)
The average of population genetic parameters for 50 selected SNPs. (DOC 46 kb)