Non-identifiability of identity coefficients at biallelic loci

Inferring biological kinship in ancient datasets: comparing the response of ancient DNA-specific software packages to low coverage data

BACKGROUND

The inference of biological relations between individuals is fundamental to understanding past human societies. Caregiving, resource sharing and sexual behaviours are often mediated by biological kinship and yet the identification and interpretation of kin relationships in prehistoric human groups is difficult. In recent years, the advent of archaeogenetic techniques have offered a fresh approach, and when combined with more traditional osteological and interpretive archaeological methods, allows for improved interpretation of the burial practices, cultural behaviours, and societal stratification in ancient societies. Although archaeogenetic techniques are developing at pace, questions remain as to their accuracy, particularly when applied to the low coverage datasets that results from the sequencing of DNA derived from highly degraded ancient material.

RESULTS

The performance of six of the most commonly used kinship identifcation software methods was explored at a range of low and ultra low genome coverages. An asymmetrical response was observed across packages, with decreased genome coverage resulting in differences in both direction and degree of change of calculated kinship scores and thus pairwise relatedness estimates are dependant on both package used and genome coverage. Methods reliant upon genotype likelihoods methods (lcMLkin, NGSrelate and NGSremix) show a decreased level of prediction at coverage below 1x, although were consistent in the particular relationships identified at these coverages when compared to the pseudohaploid reliant methods tested (READ, the Kennett 2017 method and TKGWV2.0). The three pseudohaploid methods show predictive potential at coverages as low as 0.05x, although the accuracy of the relationships identified is questionable given the increase in the number of relationships identifIed at the low coverage (type I errors).

CONCLUSION

Two pseudohaploid methods (READ and Kennett 2017) show relatively consistent inference of kin relationships at low coverage (0.5x), with READ only showing a significant performance drop off at ultralow coverages (< 0.2x). More generally, our results reveal asymmetrical kinship classifications in some software packages even at high coverages, highlighting the importance of applying multiple methods to authenticate kin relationships in ancient material, along with the continuing need to develop laboratory methods that maximise data output for downstream analyses.

Estimating FST and kinship for arbitrary population structures

F_ST and kinship are key parameters often estimated in modern population genetics studies in order to quantitatively characterize structure and relatedness. Kinship matrices have also become a fundamental quantity used in genome-wide association studies and heritability estimation. The most frequently-used estimators of F_ST and kinship are method-of-moments estimators whose accuracies depend strongly on the existence of simple underlying forms of structure, such as the independent subpopulations model of non-overlapping, independently evolving subpopulations. However, modern data sets have revealed that these simple models of structure likely do not hold in many populations, including humans. In this work, we analyze the behavior of these estimators in the presence of arbitrarily-complex population structures, which results in an improved estimation framework specifically designed for arbitrary population structures. After generalizing the definition of F_ST to arbitrary population structures and establishing a framework for assessing bias and consistency of genome-wide estimators, we calculate the accuracy of existing F_ST and kinship estimators under arbitrary population structures, characterizing biases and estimation challenges unobserved under their originally-assumed models of structure. We then present our new approach, which consistently estimates kinship and F_ST when the minimum kinship value in the dataset is estimated consistently. We illustrate our results using simulated genotypes from an admixture model, constructing a one-dimensional geographic scenario that departs nontrivially from the independent subpopulations model. Our simulations reveal the potential for severe biases in estimates of existing approaches that are overcome by our new framework. This work may significantly improve future analyses that rely on accurate kinship and F_ST estimates.

Kinship coefficients and F_ST, which measure relatedness and population structure, respectively, are important quantities needed to accurately perform various analyses on genetic data, including genome-wide association studies and heritability estimation. However, existing estimators require restrictive assumptions of independence that are not met by real human and other datasets. In this work we find that existing estimators can be severely biased under reasonable scenarios, first by theoretically determining their properties, and then using an admixture simulation to illustrate our findings. In particular, we find that existing F_ST estimators are downwardly biased, and that existing kinship matrix estimators have related biases that are on average downward and of similar magnitude but vary for every pair of individuals. These insights led us to a new estimation framework for kinship and F_ST that is practically unbiased for any population structure, as demonstrated by theory and simulations. Our new approaches—available as open-source R packages—are easy to use and are more widely applicable than existing approaches, and they are likely to improve downstream analyses that require accurate kinship and F_ST estimates.

Fast and accurate relatedness estimation from high-throughput sequencing data in the presence of inbreeding

BACKGROUND

The estimation of relatedness between pairs of possibly inbred individuals from high-throughput sequencing (HTS) data has previously not been possible for samples where we cannot obtain reliable genotype calls, as in the case of low-coverage data.

RESULTS

We introduce ngsRelateV2, a major revision of ngsRelateV1, a program that originally allowed for estimation of relatedness from HTS data among non-inbred individuals only. The new revised version takes into account the possibility of individuals being inbred by estimating the 9 condensed Jacquard coefficients along with various other relatedness statistics. The program is threaded and scales linearly with the number of cores allocated to the process.

CONCLUSION

The program is available as an open source C/C++ program under the GPL license and hosted at https://github.com/ANGSD/ngsRelate. To facilitate easy analysis, the program is able to work directly on the most commonly used container formats for raw sequence (BAM/CRAM) and summary data (VCF/BCF).

Estimation of the relatedness coefficients from biallelic markers, application in plant mating designs

The problem of inferring the relatedness distribution between two individuals from biallelic marker data is considered. This problem can be cast as an estimation task in a mixture model: at each marker the latent variable is the relatedness state, and the observed variable is the genotype of the two individuals. In this model, only the prior proportions are unknown, and can be obtained via ML estimation using the EM algorithm. When the markers are biallelic and the data unphased, the identifiability of the model is known not to be guaranteed. In this article, model identifiability is investigated in the case of phased data generated from a crossing design, a classical situation in plant genetics. It is shown that identifiability can be guaranteed under some conditions on the crossing design. The adapted ML estimator is implemented in an R package called Relatedness. The performance of the ML estimator is evaluated and compared to that of the benchmark moment estimator, both on simulated and real data. Compared to its competitor, the ML estimator is shown to be more robust and to provide more realistic estimates.