- Short communication
- Open Access
A novel recursive algorithm for the calculation of the detailed identity coefficients
- Luis Alberto García-Cortés^{1}Email author
https://doi.org/10.1186/s12711-015-0108-6
© García-Cortés; licensee BioMed Central. 2015
- Received: 8 March 2014
- Accepted: 12 February 2015
- Published: 29 April 2015
Abstract
Background
A recursive algorithm to calculate the fifteen detailed coefficients of identity is introduced. Previous recursive procedures based on the generalized coefficients of kinship provided the detailed coefficients of identity under the assumption that the two individuals were not an ancestor of each other.
Findings
By using gametic relationships to include three, four or two pairs of gametes, we can obtain these coefficients for any pair of individuals. We have developed a novel linear transformation that allows for the calculation of pairwise detailed identity coefficients for any pedigree given the gametic relationships. We illustrate the procedure using the well-known pedigree of Julio and Mencha, which contains 20 Jicaque Indians of Honduras, to calculate their detailed coefficients.
Conclusions
The proposed algorithm can be used to calculate the detailed identity coefficients of two or more individuals with any pedigree relationship.
Keywords
- Recursive Formula
- Recursive Function
- Recursive Algorithm
- Detailed Coefficient
- Kinship Coefficient
Background
The 15 detailed states of identity were first described by Harris [1] and Gillois [2]. Throughout this paper, they will be referred to identity coefficients as described by Cockerham [3].
To circumvent the absence of recurrence rules to obtain identity coefficients, Karigl [4], following the rules for generalized kinship coefficients [1,3], obtained identity coefficients using a triangular linear transformation. This transformation provided the 9 condensed identity coefficients for any pair of individuals, but the 15 detailed identity coefficients can only be calculated using Karigl’s method under the assumption that neither of the two individuals is an ancestor of the other. To obtain generalized kinship coefficients, Lange and Sinsheimer [5] described an alternative way, which can calculate the detailed coefficients of identity without this limitation. Unfortunately, the implementation of the latter method is not a triangular linear transformation.
An alternative to using multiple kinship coefficients is the use of multiple gametic relationships. These relationships, called chromosome pedigrees by Donnelly [6], have been succesfully used to account for dominance in linear models [7-9]. Here, we use these multiple gametic relationships to develop a new method to calculate the 15 detailed identity coefficients.
Generalized gametic relationships
where a and b denote two gametes in the pedigree. Both g and h are the direct ancestral gametes of a, that is, the gametes of the father or the mother if a is a paternal or maternal gamete, respectively. Although Equations (1) and (2) are closely related, they have different interpretations.
- 1.
In any n-way relationship, merge groups sharing a given gamete, that is, ψ _{ a b,a c }=ψ _{ abac }
- 2.
In any n-way relationship, discard repeated gametes, i.e. ψ _{ aabc }=ψ _{ abc } or ψ _{ a a,b c }=ψ _{ a,b c }.
- 3.
Given that the probability of a gamete to be IBD to itself is 1, discard groups of identity including a single gamete, i.e. ψ _{ a,b c }=ψ _{ bc }.
- 4.
Identities with a single gamete are 1 and identities with two or more founder gametes at the same group are 0.
- 5.
Calculate \(\psi _{a\theta }=\frac {1}{2}\left (\psi _{g\theta }+\psi _{h\theta }\right)\), where θ stands for any identity pattern and a for a gamete of the youngest individual. For instance \(\psi _{ab,cd}=\frac {1}{2}\left (\psi _{gb,cd}+\psi _{hb,cd}\right)\) or \(\psi _{abc}=\frac {1}{2}\left (\psi _{gbc}+\psi _{hbc}\right)\).
As long as these rules are correct and regardless of the number of gametes involved, they can be used to calculate identities involving more than two individuals.
The detailed identity coefficients
The detailed identity coefficients δ _{ i } are the probabilities of each state S _{ i } for a given pair of individuals.
The coefficient matrix in Equation (4) is triangular and is equivalent to formula (8) in [4]. Nevertheless, the right hand side here includes the generalized gametic relationships instead of their kinship counterparts.
An example: the detailed identity coefficients of the Jicaque Indians
Detailed coefficients of identity for four pairs involving Julio, Mencha and their progeny
Julio-Mencha | Progeny-progeny | Julio-progeny | Mencha-progeny | |
---|---|---|---|---|
δ _{1} | 0.01025 | 0.06580 | 0.03467 | 0.03467 |
δ _{2} | 0.02393 | 0.04333 | 0.08252 | 0.00000 |
δ _{3} | 0.02490 | 0.04333 | 0.00000 | 0.08252 |
δ _{4} | 0.02393 | 0.04333 | 0.06665 | 0.06665 |
δ _{5} | 0.02490 | 0.04333 | 0.08228 | 0.08228 |
δ _{6} | 0.00708 | 0.00729 | 0.00000 | 0.00000 |
δ _{7} | 0.05103 | 0.02383 | 0.00000 | 0.00000 |
δ _{8} | 0.05103 | 0.02383 | 0.00000 | 0.00000 |
δ _{9} | 0.05127 | 0.24713 | 0.08228 | 0.06665 |
δ _{10} | 0.10937 | 0.15900 | 0.29248 | 0.00000 |
δ _{11} | 0.16992 | 0.15900 | 0.00000 | 0.29248 |
δ _{12} | 0.00708 | 0.00729 | 0.06665 | 0.08228 |
δ _{13} | 0.05103 | 0.02383 | 0.00000 | 0.29248 |
δ _{14} | 0.05103 | 0.02383 | 0.29248 | 0.00000 |
δ _{15} | 0.34326 | 0.08582 | 0.00000 | 0.00000 |
Program source codes written in Fortran 90 and examples are available at https://github.com/agarcor21/IdentityCoefficients.
Discussion
The identity by descent of two gametes [6-9] is conceptually simpler than the kinship between individuals because it avoids any random sampling of genes. Although the number of gametes is larger than the number of individuals, its simple definition results in more intuitive recursive formulae and easier calculations.
Here, we have provided a novel method to derive the detailed identity coefficients [1,2]. Karigl derived a simple transformation to calculate these coefficients from the multiple kinship coefficients, which provided the 15 identity coefficients when the individuals were not ancestor-descendant related. Lange and Sinsheimer [5] first provided exact equations for these 15 coefficients without any assumption. Depending on the purpose of the study, researchers have used 7 or 9 condensed identity coefficients. For instance, in order to apply the identity coefficients to dominance models, 7 condensed coefficients can be used [8,9]. To our knowledge, the algorithm proposed in [9] has never been used to calculate the 15 detailed coefficients, but the formulas presented in their paper could be used to obtain them in a direct way.
We have proposed using multiple gametic relationships to calculate the 15 detailed identity coefficients. The linear transformation in this paper is similar to equation 8 of Karigl’s approach [4]. Reordering the rows in Equation (4) to follow Karigl’s pattern [4], both formulae only differ on the right hand side. Note that the meaning of both formulae is rather different.
The procedure presented here has been successfully implemented on small pedigrees. In large data sets, the number of calls to the recursive function will depend on the number of generations and the structure of the pedigree. Implementing the procedure in such scenarios is beyond our goal, but path counting or graph theory-based methods [14,15] have been developed to improve the computing efficiency of this calculation.
Declarations
Acknowledgements
This work was funded by grant CGL2012-39861-C02-02 from Ministerio de Economía y Competitividad. We are most grateful to Ángeles de Cara, Luis Gómez-Raya and Andrés Legarra for comments on the manuscript. We are also grateful to Rohan Fernando for providing the source code of one of the referred algorithms.
Authors’ Affiliations
References
- Harris DL. Genotypic covariances between inbred relatives. Genetics. 1964; 50:1319–48.PubMed CentralPubMedGoogle Scholar
- Gillois M. La relation d’identité en génétique. Ann Inst Henri Poincaré. 1964; B2:1–94.Google Scholar
- Cockerham CC. Higher order probability functions of identity of alleles by descent. Genetics. 1971; 69:235–46.PubMed CentralPubMedGoogle Scholar
- Karigl G. A recursive algorithm for the calculation of identity coefficients. Ann Hum Genet. 1981; 45:299–305.View ArticlePubMedGoogle Scholar
- Lange K, Sinsheimer JS. Calculation of genetic identity coefficients. Ann Hum Genet. 1992; 56:339–46.View ArticlePubMedGoogle Scholar
- Donnelly KP. The probability that related individuals share some section of genome identical by descent. Theor Popul Biol. 1983; 23:34–63.View ArticlePubMedGoogle Scholar
- Smith SP, Allaire FR. Efficient selection rules to increase non-linear merit: application to selection. Genet Sel Evol. 1985; 17:387–406.View ArticlePubMed CentralPubMedGoogle Scholar
- de Boer IJM, Hoeschele I. Genetic evaluation methods for populations with dominance and inbreeding. Theor Appl Genet. 1993; 86:245–258.View ArticlePubMedGoogle Scholar
- Lo LL, Fernando RL, Cantet RJC, Grossman M. Theory for modelling means and covariances in a two-breed population with dominance inheritance. Theor Appl Genet. 1995; 90:49–62.View ArticlePubMedGoogle Scholar
- Emik L, Terrill CE. Systematic procedures for calculating inbreeding coefficients. J Hered. 1949; 40:51–55.PubMedGoogle Scholar
- Jacquard A. The genetic structure of populations. New York: Springer; 1973.Google Scholar
- Chapman AM, Jacquard AM. Un isolat d’Amérique Centrale: Les indiens Jicaques du Honduras. I.N.E.P., editor. Génétique et Populations. Presses Universitaires de France: 1971.Google Scholar
- Nadot R, Vaysseix G. Apparentement et identité, algorithme de calcul des coefficients d’identité. Biometrics. 1973; 29:347–59.View ArticleGoogle Scholar
- Cheng E, Elliott B, Ozsoyoglu ZM. Efficient computation of kinship and identity coefficients on large pedigrees. J Bioinform Comput Biol. 2009; 7:429–53.View ArticlePubMedGoogle Scholar
- Abney M. A graphical algorithm for fast computation of identity coefficients and generalized kinship coefficients. Bioinformatics. 2009; 25:1561–3.View ArticlePubMed CentralPubMedGoogle Scholar
Copyright
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.