Evaluation of the allele-sharing approach, known as the IBS method, in kinship analysis

VCF2PCACluster: a simple, fast and memory-efficient tool for principal component analysis of tens of millions of SNPs

Principal component analysis (PCA) is an important and widely used unsupervised learning method that determines population structure based on genetic variation. Genome sequencing of thousands of individuals usually generate tens of millions of SNPs, making it challenging for PCA analysis and interpretation. Here we present VCF2PCACluster, a simple, fast and memory-efficient tool for Kinship estimation, PCA and clustering analysis, and visualization based on VCF formatted SNPs. We implemented five Kinship estimation methods and three clustering methods for its users to choose from. Moreover, unlike other PCA tools, VCF2PCACluster possesses a clustering function based on PCA result, which enabling users to automatically and clearly know about population structure. We demonstrated the same accuracy but a higher performance of this tool in performing PCA analysis on tens of millions of SNPs compared to another popular PLINK2 software, especially in peak memory usage that is independent of the number of SNPs in VCF2PCACluster.

Development of a new screening method for faster kinship analyses in mass disasters: a proof of concept study

Kinship analysis in forensics is based on the calculation of the respective kinship indices. However, this calculation is only possible when the subject under identification has been associated with a particular population, whose allele frequency data is available for the particular set of STR markers used in the forensic practices. In the case of mass disasters, where a large number of individuals are to be identified, gathering the population frequency data and calculating the kinship indices can be an intricate process which requires a lot of time and huge resources. The new method of allele matching cut off score (AMCOS) developed in this study is based on the allele sharing approach. This approach simply refers to the number of shared alleles (1 or 2) between the two individuals; also known as identical by state (IBS) alleles which might have been inherited from a recent common ancestor in which the alleles are identical by descendent (IBD). In case of mass disasters, this method can be used to narrow down the number of pairs (dead and alive) to be matched for kinship without using the allele frequency data. The results obtained from this method could further be confirmed by LR based method, which uses the allele frequency data of the respective population of the pairs being tested for kinship. AMCOS method has been tested for its sensitivity, specificity and various other statistical parameters and has shown promising values for the same in various types of kinship analyses. This ascertains the authenticity and potential use of this method in forensic practice but only after its validation in a larger sample size. AMCOS method has been tested on siblings and grandparent-grandchildren by using autosomal and X-STR markers both, as the reference samples from the parents cannot always be available for the identification. The present study also compared the results shown by the autosomal and X-STR markers in siblings and grandparent-grandchildren identification, thereby suggesting the use of better set of markers on the basis of obtained values of various statistical parameters.

Pairwise kinship analysis of 17 pedigrees using massively parallel sequencing

With the tremendous development of massively parallel sequencing (MPS) in the last decade, it has been widely applied in basic science, clinical diagnostics, microbial genomics, as well as forensic genetics. MPS has lots of advantages that may facilitate the kinship analysis. In this study, 243 Chinese Han individuals from 17 families were involved and sequenced using the ForenSeq™ DNA Signature Prep Kit (Verogen, Inc., San Diego, USA), which provided the sequence information of 27 autosomal STRs (A-STRs), 7 X chromosomal STRs (X-STRs), 24 Y chromosomal STRs (Y-STRs) and 94 identity-informative SNPs (iSNPs). A total of 275 pairs of parent-child, 123 pairs of full siblings, 1 pair of twins, 1 pair of half siblings, 158 pairs of grandparent-grandchild, 222 pairs of uncle/aunt-nephew/niece and 121 pairs of first cousins, as well as 701 pairs of unrelated individuals were identified. Using both likelihood ratio (LR) and identical by state (IBS) methods, the kinship analysis was conducted among these relative and non-relative pairs based on the A-STRs and SNPs. As a result, the ForenSeq Signature Kit could solve the analysis of parent-child (t1 = -4, t2 = 4), full siblings (t1 = -2, t2 = 2) and most second-degree kinships (t1 = -1, t2 = 1) using the LR method. When the IBS method was applied, 123 full sibling pairs had a higher average IBS value than other kinship groups in this study. And the IBS method could play a role in the testing of parent-child and full siblings.

Combined effects of multiple linked loci on pairwise sibling tests

The advanced multiplex STR system, PowerPlex Fusion, includes four linked locus pairs. The conventional Identifiler system has one pair of linked loci. Therefore, sibling tests conducted using the advanced system might be more affected by linkage than those conducted using the conventional system. This study simulated single and combined effects of the four linked locus pairs on pairwise sibling tests. Simulated genotypes of 100,000 pairs of full siblings and nonrelatives were constructed according to allele frequencies of the Japanese population. The single linkage effect was evaluated for simulated genotype data by calculating both the likelihood ratio accounting for the linkage between two loci and the likelihood ratio ignoring the linkage. The combined effect was obtained by multiplication of the respective single effects. Furthermore, we investigated the possibility that ignoring the linkage affects subject classification by introducing a scale of the likelihood ratio into sibling tests. The single effect in the Identifiler analysis was 0.645-1.746 times if the linkage was ignored. Overestimations and underestimations were predictable from the identical-by-state status at two linked loci. The combined effect in the PowerPlex Fusion analysis was 0.217-7.390 times. Ignoring the linkage rarely caused a false conclusive or inconclusive result, even from PowerPlex Fusion analysis. Application of the advanced system improved sibling tests considerably. The additional examined loci were more beneficial than the adverse effect of the linkage derived from the four linked locus pairs.

Allelic Imbalance of mRNA Associated with α2-HS Glycoprotein (Fetuin-A) Polymorphism

Alpha 2-HS glycoprotein (AHSG), also designated as fetuin-A, exhibits polymorphism in population genetics consisting of two major alleles of AHSG^∗1 and AHSG^∗2. The serum level in the AHSG^∗1 homozygote is significantly higher than that of the AHSG^∗2 homozygote. This study examined the molecular mechanism for the cis-regulatory expression. To quantitate allele-specific mRNA in intra-assays of the heterozygote, RT-PCR method employing primers that were incorporated to the two closely located SNPs was developed. The respective magnitudes of AHSG^∗1 to AHSG^∗2 in the liver tissues and hepatic culture cells of PLC/PRF/5 were determined quantitatively as 2.5-fold and 6.2-fold. The mRNA expressional difference of two major alleles was observed, which is consistent with that in the serum level. The culture cells carried heterozygous genotypes in rs4917 and rs4918, but homozygous one in rs2248690. It was unlikely that the imbalance was derived from the SNP located in the promotor site. Furthermore, to investigate the effect of mRNA degradation, RNA synthesis in the cell culture was inhibited potently by the addition of actinomycin-D. No marked change was apparent between the two alleles. The results indicated that the cis-regulatory expressional difference is expected to occur at the level of transcription or splicing of mRNA.