A study of contemporary levels and temporal trends in inbreeding in the Tangier Island, Virginia, population using pedigree data and isonymy

In this study we describe inbreeding in a large pedigree from Tangier Island, Virginia, in which we compare two commonly used methods to estimate inbreeding in humans: pedigree and isonymy (identical surnames of spouses). Genealogical data on 3,512 individuals dating back to 1722 were used. Using the pedigree method, we determined an average inbreeding coefficient (F) of 0.00873 for the community as a whole, and 0.018 for inbred individuals. Analysis of temporal trends showed that inbreeding began around 1800 and peaked at 0.0109 in 1824-1849 and 1875-1899. Thereafter, inbreeding steadily declined to 0.00565 in 1975-1997. Analysis of pedigree structure complexity over time showed that close consanguinity contributes to inbreeding in the earlier cohorts, and remote consanguinity accounts for much of the inbreeding in the later cohorts. The number of common ancestors increases over time, as does the number of paths connecting inbred individuals to these common ancestors. Inbreeding estimates based on the isonymy approach yielded a 2.2-fold higher value of F (0.01945) compared to the pedigree method. Total isonymy estimates over 25-year cohorts overestimated inbreeding values from pedigree data between 1. 5-8-fold. We speculate that the overestimation is probably due to the inability of our data to satisfy the method's assumption of monophyletic origin of each surname. In conclusion, inbreeding in the Tangier Island population is consistent with the isolated nature of its population, and temporal trends reflect patterns in emigration and a breakdown in isolation over time.     

Demographic genetic approach in anthropologic research. II. A method for calculating the coefficient of consanguinity in a complex nongraphic pedigree

A simplified method of manual calculation for inbreeding coefficient (and/or relationship) from complicated pedigree tables presented in a nongraphic mode is offered. The method is based on the use of positional dual codes of ancestors in which all genealogical links of the proband are reflected.     

Quantification of inbreeding due to distant ancestors and its detection using dense single nucleotide polymorphism data

Inbreeding depression, which refers to reduced fitness among offspring of related parents, has traditionally been studied using pedigrees. In practice, pedigree information is difficult to obtain, potentially unreliable, and rarely assessed for inbreeding arising from common ancestors who lived more than a few generations ago. Recently, there has been excitement about using SNP data to estimate inbreeding (F) arising from distant common ancestors in apparently "outbred" populations. Statistical power to detect inbreeding depression using SNP data depends on the actual variation in inbreeding in a population, the accuracy of detecting that with marker data, the effect size, and the sample size. No one has yet investigated what variation in F is expected in SNP data as a function of population size, and it is unclear which estimate of F is optimal for detecting inbreeding depression. In the present study, we use theory, simulated genetic data, and real genetic data to find the optimal estimate of F, to quantify the likely variation in F in populations of various sizes, and to estimate the power to detect inbreeding depression. We find that F estimated from runs of homozygosity (Froh), which reflects shared ancestry of genetic haplotypes, retains variation in even large populations (e.g., SD=0.5% when Ne=10,000) and is likely to be the most powerful method of detecting inbreeding effects from among several alternative estimates of F. However, large samples (e.g., 12,000-65,000) will be required to detect inbreeding depression for likely effect sizes, and so studies using Froh to date have probably been underpowered.     

Inbreeding and genetic disease in Sottunga, Finland

The contribution of inbreeding to the prevalence of recessive genetic diseases in the Aland Island parish of Sottunga is investigated. Genealogical data for 3,030 individuals spanning up to 15 generations were used to estimate inbreeding. This small island community shows a low average inbreeding value of .0031 for the period 1725-1975. A cohort analysis shows that inbreeding increased from 1750 to 1900, when maximum inbreeding for those born in Sottunga reached .0057. A sharp decline in inbreeding occurred thereafter. Individuals with island-born parents made the largest contributions to inbreeding in all time periods compared to those with one or two migrant parents. These trends are consistent with changing migration patterns and isolate breakdown in Aland since 1900. An analysis of pedigree development demonstrates that remote consanguinity contributed more to inbreeding through time than close consanguinity. Both the number of common ancestors and the number of paths of relationship between spouses increased dramatically through time, the latter at a much faster rate. The contribution to average inbreeding per path, however, diminished rapidly through time. This analysis indicates that inbreeding does not account for the high incidence of autosomal recessive disorders, such as tapetoretinal disease, found in the parish.     

Analysis of genealogical structure of populations. I. A method of collecting genealogical data in numerical form

A method for collecting genealogical data with respect to an individual, a family, and members of the whole population is suggested. The essence of vertical pedigree construction consists of the same type of steps for filling in data (in the fixed order which excludes skips in the enumeration of lines of descent) about the father and the mother of the next ancestor. Each number in the received ordered list of ancestors uniquely determines a path (line of descent) to the given pedigree member. The path is explicitly described by a sequence of digits 0 and 1 (that corresponds to the sequence of fathers and mothers in the line of descent) at binary notation of this number. As a result, a pedigree is presented as a set of numbered rows that contain information, which uniquely identifies direct ancestors as individual persons. Results of joining separate pedigrees are recorded as a family list that contains lists of children for each parental pair. A pair of parents (more exactly, pointers of their families in the previous generation and numbers of pair members in their families) plays the role of the family "heading." Such a family list permits one to trace lines of descent and relationships for any population members presented in the list. It contains all genealogical information within the bounds of the study in a compact form. Here the process of collection requires considerably less time than traditional graphic representation of pedigrees. In addition, due to repeated checks of data during accumulation of material, error is minimized. Using pedigrees that have been collected, it is possible to calculate the coefficient of inbreeding manually. In connection with the wide prevalence of personal computers at present, it is also important that the data received are in fact ready to direct input to a computer for further automated data processing.     

Efficient path-based computations on pedigree graphs with compact encodings

A pedigree is a diagram of family relationships, and it is often used to determine the mode of inheritance (dominant, recessive, etc.) of genetic diseases. Along with rapidly growing knowledge of genetics and accumulation of genealogy information, pedigree data is becoming increasingly important. In large pedigree graphs, path-based methods for efficiently computing genealogical measurements, such as inbreeding and kinship coefficients of individuals, depend on efficient identification and processing of paths. In this paper, we propose a new compact path encoding scheme on large pedigrees, accompanied by an efficient algorithm for identifying paths. We demonstrate the utilization of our proposed method by applying it to the inbreeding coefficient computation. We present time and space complexity analysis, and also manifest the efficiency of our method for evaluating inbreeding coefficients as compared to previous methods by experimental results using pedigree graphs with real and synthetic data. Both theoretical and experimental results demonstrate that our method is more scalable and efficient than previous methods in terms of time and space requirements.     

Purging of inbreeding depression within the Irish Holstein-Friesian population

The objective of this study was to investigate whether inbreeding depression in milk production or fertility performance has been partially purged due to selection within the Irish Holstein-Friesian population. Classical, ancestral (i.e., the inbreeding of an individual''s ancestors according to two different formulae) and new inbreeding coefficients (i.e., part of the classical inbreeding coefficient that is not accounted for by ancestral inbreeding) were computed for all animals. The effect of each coefficient on 305-day milk, fat and protein yield as well as calving interval, age at first calving and survival to second lactation was investigated. Ancestral inbreeding accounting for all common ancestors in the pedigree had a positive effect on 305-day milk and protein yield, increasing yields by 4.85 kg and 0.12 kg, respectively. However, ancestral inbreeding accounting only for those common ancestors, which contribute to the classical inbreeding coefficient had a negative effect on all milk production traits decreasing 305-day milk, fat and protein yields by -8.85 kg, -0.53 kg and -0.33 kg, respectively. Classical, ancestral and new inbreeding generally had a detrimental effect on fertility and survival traits. From this study, it appears that Irish Holstein-Friesians have purged some of their genetic load for milk production through many years of selection based on production alone, while fertility, which has been less intensely selected for in the population demonstrates no evidence of purging.     

Inbreeding avoidance in an isolated indigenous Zapotec community in the valley of Oaxaca, southern Mexico

We analyzed inbreeding using surname isonymy in an indigenous genetic isolate. The subjects were residents of a rural Zapotec-speaking community in the valley of Oaxaca, southern Mexico. The community can be classified as a genetic isolate with an average gene flow of < or = 3% per generation. Surnames were collected for individuals in each household in pedigree form using the culturally traditional patronym-matronym naming. Estimation of inbreeding from surname isonymy is facilitated by the traditional patronym-matronym name assignment among indigenous Mexican populations. A total of 2,149 individuals had valid surname patronym-matronym pairings, including 484 deceased ancestors. Surname isonymy analysis methods were used to estimate total inbreeding and to segregate it into random and nonrandom components. The surname isonymy coefficient computed from 119 isonymous surname pairings (119/2,149) was 0.0554. The estimated inbreeding coefficient from surname isonymy was 0.0138 (0.0554/4). The random and nonrandom components of inbreeding were F(r) = 0.0221 and F(n) = -0.0091, respectively. The results suggest that consanguinity is culturally avoided. Nonrandom inbreeding decreased total inbreeding by about 41%. Total estimated inbreeding by surname isonymy was 0.0138, which is similar to inbreeding estimated from a sample of pedigrees, 0.01. Socially prescribed inbreeding avoidance substantially lowered total F through negative nonrandom inbreeding. Even in the situation of genetic isolation and small effective population size (N(e)), estimated inbreeding is lower than may have otherwise occurred if inbreeding were only random. However, among the poorest individuals, socially prescribed jural rules for inbreeding avoidance failed to operate. Thus the preponderance of inbreeding appears to occur among the poor, economically disadvantaged in the community.     

Inbreeding avoidance behaviors

Inbreeding is defined as mating between individuals related by common ancestry. Thus, the degree to which a particular mating is inbred depends on how far back in a pedigree one begins counting common ancestors. In general practice, the term inbreeding is used to describe mating between close relatives (first cousins or closer). Animal breeders have known for centuries that inbreeding causes a loss of constitutional vigor and fertility in domestic livestock. A growing literature now demonstrates that the offspring of matings between close relatives in species of undomesticated birds and mammals are less fit than outbred offspring. The deleterious consequences of inbreeding suggest the possibility that many species have evolved behaviors that lower the frequency of inbreeding.     

Inbreeding, pedigree size, and the most recent common ancestor of humanity

How many generations ago did the common ancestor of all present-day individuals live, and how does inbreeding affect this estimate? The number of ancestors within family trees determines the timing of the most recent common ancestor of humanity. However, mating is often non-random and inbreeding is ubiquitous in natural populations. Rates of pedigree growth are found for multiple types of inbreeding. This data is then combined with models of global population structure to estimate biparental coalescence times. When pedigrees for regular systems of mating are constructed, the growth rates of inbred populations contain Fibonacci n-step constants. The timing of the most recent common ancestor depends on global population structure, the mean rate of pedigree growth, mean fitness, and current population size. Inbreeding reduces the number of ancestors in a pedigree, pushing back global common ancestry times. These results are consistent with the remarkable findings of previous studies: all humanity shares common ancestry in the recent past.     

Quantifying realized inbreeding in wild and captive animal populations

Most molecular measures of inbreeding do not measure inbreeding at the scale that is most relevant for understanding inbreeding depression—namely the proportion of the genome that is identical-by-descent (IBD). The inbreeding coefficient F_Ped obtained from pedigrees is a valuable estimator of IBD, but pedigrees are not always available, and cannot capture inbreeding loops that reach back in time further than the pedigree. We here propose a molecular approach to quantify the realized proportion of the genome that is IBD (propIBD), and we apply this method to a wild and a captive population of zebra finches (Taeniopygia guttata). In each of 948 wild and 1057 captive individuals we analyzed available single-nucleotide polymorphism (SNP) data (260 SNPs) spread over four different genomic regions in each population. This allowed us to determine whether any of these four regions was completely homozygous within an individual, which indicates IBD with high confidence. In the highly nomadic wild population, we did not find a single case of IBD, implying that inbreeding must be extremely rare (propIBD=0–0.00094, 95% CI). In the captive population, a five-generation pedigree strongly underestimated the average amount of realized inbreeding (F_Ped=0.013<propIBD=0.064), as expected given that pedigree founders were already related. We suggest that this SNP-based technique is generally useful for quantifying inbreeding at the individual or population level, and we show analytically that it can capture inbreeding loops that reach back up to a few hundred generations.     

FamSeq: A Variant Calling Program for Family-Based Sequencing Data Using Graphics Processing Units

Various algorithms have been developed for variant calling using next-generation sequencing data, and various methods have been applied to reduce the associated false positive and false negative rates. Few variant calling programs, however, utilize the pedigree information when the family-based sequencing data are available. Here, we present a program, FamSeq, which reduces both false positive and false negative rates by incorporating the pedigree information from the Mendelian genetic model into variant calling. To accommodate variations in data complexity, FamSeq consists of four distinct implementations of the Mendelian genetic model: the Bayesian network algorithm, a graphics processing unit version of the Bayesian network algorithm, the Elston-Stewart algorithm and the Markov chain Monte Carlo algorithm. To make the software efficient and applicable to large families, we parallelized the Bayesian network algorithm that copes with pedigrees with inbreeding loops without losing calculation precision on an NVIDIA graphics processing unit. In order to compare the difference in the four methods, we applied FamSeq to pedigree sequencing data with family sizes that varied from 7 to 12. When there is no inbreeding loop in the pedigree, the Elston-Stewart algorithm gives analytical results in a short time. If there are inbreeding loops in the pedigree, we recommend the Bayesian network method, which provides exact answers. To improve the computing speed of the Bayesian network method, we parallelized the computation on a graphics processing unit. This allowed the Bayesian network method to process the whole genome sequencing data of a family of 12 individuals within two days, which was a 10-fold time reduction compared to the time required for this computation on a central processing unit.

This is a PLOS Computational Biology Software Article

     

Inbreeding as measured by isonymy, pedigrees, and population size in Törbel, Switzerland.

Törbel provides an interesting test case for the study of the relationship between inbreeding measured by pedigrees and inbreeding measured by isonymy. At the start of this investigation, we were aware that isonymy could introduce biases into the calculation of the inbreeding coefficient in either direction. However, it was expected that in Switzerland, inbreeding from isonymy would be an overestimate due to patrilocal residence and polyphyletic names. One way of dealing with this problem [13] was not to be concerned with the absolute value of inbreeding but only in the difference between estimates. Any bias introduced in the estimate itself disappears in such comparisons, so that a trend of inbreeding can be ascertained correctly. However, it was considered equally important to subject several populations to both a complete pedigree analysis and an isonymic analysis to determine the relationship between estimates of inbreeding. Despite the fact that several authors (Swedlund [18], for example) warned users of isonymy to exercise caution, the careless application of isonymy still persists. In the present study, estimates of inbreeding from isonymy were brought into line with other methods based on pedigree analysis and population size. However, it was possible to do this only in Törbel where pedigree depth was extensive and relatively complete. Similar corrections are possible only when the distribution of mono- and polyphyletic names is known and when migration data are reliable. If the trouble is taken to make these corrections, the same time and effort might as well be spent in pedigree analysis (when fairly complete ascertainment is possible) to achieve the same end result.     

Evaluation of genetic variation in the international Brown Swiss population

The international Brown Swiss cattle population pedigree was studied to measure genetic variations and to identify the most influential animals. Twenty-two countries provided pedigree information on 71 497 Brown Swiss bulls used for artificial insemination (AI). The total number of animals with the pedigree is 181 094. The mean inbreeding coefficient for the pedigree population was 0.77%. There was, in most cases, an increase in the mean inbreeding coefficient, with the highest value at 2.89% during the last 5-year period (2000 to 2004). The mean average relatedness for the pedigree population was 1.1%. The effective population size in 2004 was 204. There was notable variation between average generation intervals for the four parental pathways. The longest average generation interval, at 8.73 years, was observed in the sire–son pathway. The average generation interval for the whole population was 6.53 years. Most genetically influential individuals were sires. The highest contributing founder was a sire with a 3.22% contribution, and the highest contributing founder dam made a contribution of 1.75%. The effective number of founders and the effective number of ancestors were 141 and 88, respectively. The study showed that genetic variation within the pedigree population has been decreasing over recent years. Increasing the number of AI bulls with a low individual coefficient of inbreeding could help to maintain a good level of genetic variation in the Brown Swiss population.     

Population structure of the Trakehner Horse breed

The objective of this study was to examine the population structure of the Trakehner Horse breed. A total of 13 793 pedigree records were used for analysing the active breeding population and their ancestors dating back to 1950. Ancestors that were born before 1950 were called as base animals. The average generation interval was calculated as 10.2 years. The effective population size (Ne) was estimated by the increase in average year-wise inbreeding coefficient and average coancestry, respectively. Two methods were applied to estimate the effective population size: 1. Numerator-relationship-matrix (NRM), which did not consider missing ancestries. 2. Uncertain-parentage-matrix (UPM), which considered a probabilistic correction for unknown ancestors. There were no major differences between these two methods with respect to the rate of increase in inbreeding although the global levels using the UPM method were observed to be higher. Estimates for the inbreeding coefficients and the average coancestries varied little between both methods. The estimates of the effective population size per generation based on the rate of inbreeding ranged from 169 (NRM) to 150 (UPM) and 158 (NRM) to 144 (UPM) calculated by the average coancestry. From the early 1990s onwards, a strong increase in the rate of inbreeding was observed. This may be due to an increasing variance of the family size of sires and may be interpreted as a consequence of the growing use of artificial insemination. Analysing coancestries within and between the centrally managed regional breeding societies in Germany further revealed the Trakehner horse breed to be a genetically fragmented population with a main partition corresponding to formerly divided East and West Germany. The average rate of gene contributions (Thoroughbred (xx), Arab Horse breed (ox)) to the defined actual breeding population was calculated to be 22.3% xx-genes and 11.7% ox-genes.     

A FORTRAN subroutine to compute inbreeding and kinship coefficients according to the number of ancestral generations

This paper presents a FORTRAN IV subroutine to calculate inbreedingand kinship coefficients from pedigree information in a diploidpopulation without self-fertilization. The user can specifythe number of ancestral generations to be taken into account.It is thus possible to determine contributions of succeedingancestral generations to the inbreeding and kinship coefficientsunder consideration. The subroutine is based on a recursiveprocedure that generates systematically all paths connectingtwo individuals, NP and NM, whose kinship coefficient is tobe calculated (or between the father NP and the mother NM ofthe individual whose inbreeding coefficient is to be calculated).These paths obey the following conditions: (i) a given pathdoes not contain the same parent—offspring link more thanonce; (ii) the vertex of a path is an ancestor common to individualsNP and NM, with a rank lower or equal to the parameter specifiedin input. Constraints regarding the size of the corpus of genealogicaldata and the storage method are discussed, as well as the interestof this subroutine compared to the existing ones. An exampleof application is given. Received on October 20, 1988; accepted on March 21, 1989     

Quantitative Genetics Model as the Unifying Model for Defining Genomic Relationship and Inbreeding Coefficient

The traditional quantitative genetics model was used as the unifying approach to derive six existing and new definitions of genomic additive and dominance relationships. The theoretical differences of these definitions were in the assumptions of equal SNP effects (equivalent to across-SNP standardization), equal SNP variances (equivalent to within-SNP standardization), and expected or sample SNP additive and dominance variances. The six definitions of genomic additive and dominance relationships on average were consistent with the pedigree relationships, but had individual genomic specificity and large variations not observed from pedigree relationships. These large variations may allow finding least related genomes even within the same family for minimizing genomic relatedness among breeding individuals. The six definitions of genomic relationships generally had similar numerical results in genomic best linear unbiased predictions of additive effects (GBLUP) and similar genomic REML (GREML) estimates of additive heritability. Predicted SNP dominance effects and GREML estimates of dominance heritability were similar within definitions assuming equal SNP effects or within definitions assuming equal SNP variance, but had differences between these two groups of definitions. We proposed a new measure of genomic inbreeding coefficient based on parental genomic co-ancestry coefficient and genomic additive correlation as a genomic approach for predicting offspring inbreeding level. This genomic inbreeding coefficient had the highest correlation with pedigree inbreeding coefficient among the four methods evaluated for calculating genomic inbreeding coefficient in a Holstein sample and a swine sample.     

The genotypic probability structure of an individual at the head of a pedigree

A proper probabilistical proof of a generalization of Wright's classical formula relating the coefficient of inbreeding of an individual at the head of a pedigree to the genotypic probability structure of this individual at one gene locus is presented. It is shown that in general the knowledge of gene frequencies realized within the initial populations from which individuals entering the pedigree are selected at random is not sufficient to predict expected genotypic frequencies in the resulting inbred population. To treat any arbitrary situation concerning the choice of individuals and genotypes to enter the pedigree, it is necessary to determine an additional set of coefficients, which merely depend on the type of the pedigree. A basic method for computing these coefficients is outlined briefly.     

Population genetics of the population of the European north of the RSFSR. IV. The level of blood relationship in 5 villages of Pinega District, Archangel Province]

Computer program "RODAN-1" is used for inbreeding coefficient estimation. The population studied consists of two communities of 5 villages. 385 marriages were computed. The coefficient of inbreeding is 0.00145 for pedigree for rural Russian population (the Arkhangelsk region). The inverse dependence between a village size and corresponding data of inbreeding coefficient is suggested. An attempt was undertaken to estimate the genealogical information value for each pedigree and average information value for a village.