mixIndependR: a R package for statistical independence testing of loci in database of multi-locus genotypes

Background Multi-locus genotype data are widely used in population genetics and disease studies. In evaluating the utility of multi-locus data, the independence of markers is commonly considered in many genomic assessments. Generally, pairwise non-random associations are tested by linkage disequilibrium; however, the dependence of one panel might be triplet, quartet, or other. Therefore, a compatible and user-friendly software is necessary for testing and assessing the global linkage disequilibrium among mixed genetic data. Results This study describes a software package for testing the mutual independence of mixed genetic datasets. Mutual independence is defined as no non-random associations among all subsets of the tested panel. The new R package “ mixIndependR ” calculates basic genetic parameters like allele frequency, genotype frequency, heterozygosity, Hardy–Weinberg equilibrium, and linkage disequilibrium (LD) by mutual independence from population data, regardless of the type of markers, such as simple nucleotide polymorphisms, short tandem repeats, insertions and deletions, and any other genetic markers. A novel method of assessing the dependence of mixed genetic panels is developed in this study and functionally analyzed in the software package. By comparing the observed distribution of two common summary statistics (the number of heterozygous loci [K] and the number of share alleles [X]) with their expected distributions under the assumption of mutual independence, the overall independence is tested. Conclusion The package “ mixIndependR ” is compatible to all categories of genetic markers and detects the overall non-random associations. Compared to pairwise disequilibrium, the approach described herein tends to have higher power, especially when number of markers is large. With this package, more multi-functional or stronger genetic panels can be developed, like mixed panels with different kinds of markers. In population genetics, the package “ mixIndependR ” makes it possible to discover more about admixture of populations, natural selection, genetic drift, and population demographics, as a more powerful method of detecting LD. Moreover, this new approach can optimize variants selection in disease studies and contribute to panel combination for treatments in multimorbidity. Application of this approach in real data is expected in the future, and this might bring a leap in the field of genetic technology. Availability The R package mixIndependR, is available on the Comprehensive R Archive Network (CRAN) at: https://cran.r-project.org/web/packages/mixIndependR/index.html .

[1]  Manuel A. R. Ferreira,et al.  PLINK: a tool set for whole-genome association and population-based linkage analyses. , 2007, American journal of human genetics.

[2]  Richard R. Hudson,et al.  Generating samples under a Wright-Fisher neutral model of genetic variation , 2002, Bioinform..

[3]  D N Stivers,et al.  The utility of short tandem repeat loci beyond human identification: Implications for development of new DNA typing systems , 1999, Electrophoresis.

[4]  M. Gymrek,et al.  A reference haplotype panel for genome-wide imputation of short tandem repeats , 2018, Nature Communications.

[5]  Hadley Wickham,et al.  ggplot2 - Elegant Graphics for Data Analysis (2nd Edition) , 2017 .

[6]  R. Chakraborty The distribution of the number of heterozygous Loci in an individual in natural populations. , 1981, Genetics.

[7]  Bridget F B Algee-Hewitt,et al.  Linkage disequilibrium matches forensic genetic records to disjoint genomic marker sets , 2017, Proceedings of the National Academy of Sciences.

[8]  E. Thompson,et al.  Performing the exact test of Hardy-Weinberg proportion for multiple alleles. , 1992, Biometrics.

[9]  S. Cichon,et al.  Can long-range microsatellite data be used to predict short-range linkage disequilibrium? , 2002, Human molecular genetics.

[10]  Peter M Vallone,et al.  STRs vs. SNPs: thoughts on the future of forensic DNA testing , 2007, Forensic science, medicine, and pathology.

[11]  Mutation rates at 23 different short tandem repeat loci , 2003 .

[12]  M. Monteros,et al.  Development of an Alfalfa SNP Array and Its Use to Evaluate Patterns of Population Structure and Linkage Disequilibrium , 2014, PloS one.

[13]  David Levine,et al.  A high-performance computing toolset for relatedness and principal component analysis of SNP data , 2012, Bioinform..

[14]  G. Trotoux,et al.  High-density SNP-based genetic map development and linkage disequilibrium assessment in Brassica napus L , 2013, BMC Genomics.

[15]  H. Gregorius The probability of losing an allele when diploid genotypes are sampled. , 1980, Biometrics.

[16]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[17]  Bruce S Weir,et al.  Matching and partially-matching DNA profiles. , 2004, Journal of forensic sciences.

[18]  Chao Xiao,et al.  A novel multiplex assay of SNP-STR markers for forensic purpose , 2018, PloS one.

[19]  Peter M Visscher,et al.  Recent human effective population size estimated from linkage disequilibrium. , 2007, Genome research.

[20]  Jeffrey Ross-Ibarra,et al.  Genetic Data Analysis II. Methods for Discrete Population Genentic Data , 2002 .