Pairwise Kinship Research Articles

Spatial genetic structure in the Laperrine's olive (Olea europaea subsp. laperrinei), a long-living tree from the central Saharan mountains

The Laperrine's olive (Olea europaea subsp. laperrinei) is an emblematic species of the Sahelo-Saharan Mountains. Populations of this tree are locally threatened by extinction due to climatic vicissitudes and human activities, particularly in Niger and Algeria. In order to study the spatial genetic structure and the dynamics of O. e. laperrinei populations, we sampled trees in four isolated mountain ranges (Tassili n'Ajjer and Hoggar (Algeria), Tamgak and Bagzane (Niger)). A total of 237 genets were identified using nuclear microsatellites. Phylogenetic reconstruction based on plastid DNA data supported a maternal origin of O. e. laperrinei populations in South Algeria, where a higher allelic richness was observed. Based on nuclear microsatellite data, two levels of structure were revealed: first, individuals from Niger and Algeria were separated in two distinct groups; second, four less differentiated clusters corresponded to the four studied mountain ranges. These results give support to the fact that desert barriers have greatly limited long distance gene flow. Within populations, pairwise kinship coefficients were significantly correlated to geographical distance for Niger populations but not for Algerian mountains. Historical factors and habitat heterogeneity may explain the differences observed. We conclude that the Hoggar acts as an important genetic reservoir that has to be taken into account in future conservation programmes. Moreover, very isolated endangered populations (for example, Bagzane) displaying evident genetic particularities have to be urgently considered for their endemism.

Variation in demographic and fine‐scale genetic structure with population‐history stage of Hemerocallis taeanensis (Liliaceae) across the landscape

AbstractFine‐scale genetic structure (FSGS) in plant populations is expected to be influenced by variation in demographic processes across space and over time. I chose Hemerocallis taeanensis (Liliaceae), a perennial herb with a rapid population turnover, to quantify how demographic structure and FSGS change with a population's history (i.e., density). Nonaccumulative O‐ring statistic and spatial autocorrelation analysis (kinship coefficient, Fij) were used to quantify spatial patterns of individuals and FSGS in four populations belonging to two population stages (expansion and maturation) in west‐central Korea. The O‐ring function revealed that significant aggregation of individuals occurs at short spatial scales during the earlier stage of population expansion, which reflects restricted seed dispersal around maternal individuals. However, this pattern disappears as the population density increases during population maturation, probably due to a high population density. Significant evidence of FSGS was found in two populations at the stage of population expansion (Sp, a statistic which describes the rate of decrease of pairwise kinship with distance, was 0.018 and 0.029). The results show that most seeds fall around maternal plants when initially established colonists proliferate at suitable microhabitats. In contrast to this, much lower Sp values (−0.003 and 0.004) were estimated for two populations at the stage of population maturation, which may result from the overlapping of seed shadows due to high adult density. All of these results demonstrate considerable variation in within‐population demographic and genetic structures of H. taeanensis with respect to population temporal stage across the landscape.

An Arabidopsis Example of Association Mapping in Structured Samples

A potentially serious disadvantage of association mapping is the fact that marker-trait associations may arise from confounding population structure as well as from linkage to causative polymorphisms. Using genome-wide marker data, we have previously demonstrated that the problem can be severe in a global sample of 95 Arabidopsis thaliana accessions, and that established methods for controlling for population structure are generally insufficient. Here, we use the same sample together with a number of flowering-related phenotypes and data-perturbation simulations to evaluate a wider range of methods for controlling for population structure. We find that, in terms of reducing the false-positive rate while maintaining statistical power, a recently introduced mixed-model approach that takes genome-wide differences in relatedness into account via estimated pairwise kinship coefficients generally performs best. By combining the association results with results from linkage mapping in F2 crosses, we identify one previously known true positive and several promising new associations, but also demonstrate the existence of both false positives and false negatives. Our results illustrate the potential of genome-wide association scans as a tool for dissecting the genetics of natural variation, while at the same time highlighting the pitfalls. The importance of study design is clear; our study is severely under-powered both in terms of sample size and marker density. Our results also provide a striking demonstration of confounding by population structure. While statistical methods can be used to ameliorate this problem, they cannot always be effective and are certainly not a substitute for independent evidence, such as that obtained via crosses or transgenic experiments. Ultimately, association mapping is a powerful tool for identifying a list of candidates that is short enough to permit further genetic study.

Association genetics in Pinus taeda L. I. Wood property traits.

Genetic association is a powerful method for dissecting complex adaptive traits due to (i) fine-scale mapping resulting from historical recombination, (ii) wide coverage of phenotypic and genotypic variation within a single experiment, and (iii) the simultaneous discovery of loci and alleles. In this article, genetic association among single nucleotide polymorphisms (58 SNPs) from 20 wood- and drought-related candidate genes and an array of wood property traits with evolutionary and commercial importance, namely, earlywood and latewood specific gravity, percentage of latewood, earlywood microfibril angle, and wood chemistry (lignin and cellulose content), was tested using mixed linear models (MLMs) that account for relatedness among individuals by using a pairwise kinship matrix. Population structure, a common systematic bias in association studies, was assessed using 22 nuclear microsatellites. Different phenotype:genotype associations were found, some of them confirming previous evidence from collocation of QTL and genes in linkage maps (for example, 4cl and percentage of latewood) and two that involve nonsynonymous polymorphisms (cad SNP M28 with earlywood specific gravity and 4cl SNP M7 with percentage of latewood). The strongest genetic association found in this study was between allelic variation in alpha-tubulin, a gene involved in the formation of cortical microtubules, and earlywood microfibril angle. Intragenic LD decays rapidly in conifers; thus SNPs showing genetic association are likely to be located in close proximity to the causative polymorphisms. This first multigene association genetic study in forest trees has shown the feasibility of candidate gene strategies for dissecting complex adaptive traits, provided that genes belonging to key pathways and appropriate statistical tools are used. This approach is of particular utility in species such as conifers, where genomewide strategies are limited by their large genomes.

Extensive spatial genetic structure revealed by AFLP but not SSR molecular markers in the wind‐pollinated tree, Fagus sylvatica

Studies of fine-scale spatial genetic structure (SGS) in wind-pollinated trees have shown that SGS is generally weak and extends over relatively short distances (less than 30-40 m) from individual trees. However, recent simulations have shown that detection of SGS is heavily dependent on both the choice of molecular markers and the strategy used to sample the studied population. Published studies may not always have used sufficient markers and/or individuals for the accurate estimation of SGS. To assess the extent of SGS within a population of the wind-pollinated tree Fagus sylvatica, we genotyped 200 trees at six microsatellite or simple sequence repeat (SSR) loci and 250 amplified fragment length polymorphisms (AFLP) and conducted spatial analyses of pairwise kinship coefficients. We re-sampled our data set over individuals and over loci to determine the effect of reducing the sample size and number of loci used for SGS estimation. We found that SGS estimated from AFLP markers extended nearly four times further than has been estimated before using other molecular markers in this species, indicating a persistent effect of restricted gene flow at small spatial scales. However, our SSR-based estimate was in agreement with other published studies. Spatial genetic structure in F. sylvatica and similar wind-pollinated trees may therefore be substantially larger than has been estimated previously. Although 100-150 AFLP loci and 150-200 individuals appear sufficient for adequately estimating SGS in our analysis, 150-200 individuals and six SSR loci may still be too few to provide a good estimation of SGS in this species.

Analysis of genetic diversity and the determination of relationships among western Mediterranean horse breeds using microsatellite markers

The distribution of genetic diversity and the genetic relationships among western Mediterranean horse breeds were investigated using microsatellite markers. The examined sample included seven Spanish and three Italian local horse breeds and populations, plus a Spanish Thoroughbred outgroup. The total number of animals examined was 682 (on average 62 animals per breed; range 20-122). The microsatellite marker set analysed provided 128 alleles (10.7 alleles per locus). Within-breed genetic diversity was always high (>0.70), with breeds contributing about 8% of the total genetic variability. The mean molecular coancestry of the entire population examined was 0.205, Losino being the breed that contributed most. In addition to Nei's standard and Reynolds' genetic distances, pair-wise kinship distance and molecular coancestry were estimated. Remarkably similar breed rankings were obtained with all methods. Clustering analysis provided an accurate representation of the current genetic relationships among the breeds. Determining coancestry is useful for analysing genetic diversity distribution between and within breeds and provides a good framework for jointly analysing molecular markers and pedigree information. An integrated analysis was undertaken to obtain information on the population dynamics in western Mediterranean native horse populations, and to better determine conservation priorities.

Fine‐scale genetic structure and gene dispersal inferences in 10 Neotropical tree species

The extent of gene dispersal is a fundamental factor of the population and evolutionary dynamics of tropical tree species, but directly monitoring seed and pollen movement is a difficult task. However, indirect estimates of historical gene dispersal can be obtained from the fine-scale spatial genetic structure of populations at drift-dispersal equilibrium. Using an approach that is based on the slope of the regression of pairwise kinship coefficients on spatial distance and estimates of the effective population density, we compare indirect gene dispersal estimates of sympatric populations of 10 tropical tree species. We re-analysed 26 data sets consisting of mapped allozyme, SSR (simple sequence repeat), RAPD (random amplified polymorphic DNA) or AFLP (amplified fragment length polymorphism) genotypes from two rainforest sites in French Guiana. Gene dispersal estimates were obtained for at least one marker in each species, although the estimation procedure failed under insufficient marker polymorphism, limited sample size, or inappropriate sampling area. Estimates generally suffered low precision and were affected by assumptions regarding the effective population density. Averaging estimates over data sets, the extent of gene dispersal ranged from 150 m to 1200 m according to species. Smaller gene dispersal estimates were obtained in species with heavy diaspores, which are presumably not well dispersed, and in populations with high local adult density. We suggest that limited seed dispersal could indirectly limit effective pollen dispersal by creating higher local tree densities, thereby increasing the positive correlation between pollen and seed dispersal distances. We discuss the potential and limitations of our indirect estimation procedure and suggest guidelines for future studies.

Genetic diversity and population structure of Valdensinia heterodoxa, a potential biocontrol agent for salal in coastal British Columbia

Valdensinia heterodoxa is an ascomycete fungus currently being considered as a potential biocontrol agent against salal (Gaultheria shallon). To assess its effectiveness and risks as a biocontrol agent, we investigated its population structure in three geographically separate populations on Vancouver Island and in coastal mainland, British Columbia. Valdensinia heterodoxa was cultured from infected salal leaves, and single-spore cultures were obtained prior to DNA isolation. Amplified fragment length polymorphisms (AFLPs) were used to generate individual DNA fingerprints for each V. heterodoxa isolate. We found low genetic diversity, with just 30 of the 214 loci being polymorphic. There were many shared haplotypes within each population, and analysis of pairwise kinship coefficients showed that as spatial distance increased, genetic similarity decreased. Analysis of molecular variance among populations revealed significant genetic differentiation, with an F ST of 0.182, suggesting limited gene flow. The findings suggest that with low diversity and high population differentiation, the effectiveness of V. heterodoxa as a biocontrol agent to manage salal in forested areas may be limited to localized salal populations or to use in combination with other control methods.

Genetic structure and mating system of northernArbutus menziesii(Ericaceae) populations

Arbutus ( Arbutus menziesii Pursh. (Ericaceae)) is the only broadleaved evergreen tree native to Canada. It occurs in three red-listed (endangered) plant communities in British Columbia (BC), threatened by urban encroachment, fire suppression, grazing, and exotic invasive species. Its growth is sensitive to environmental changes: more severe summer drought caused by climate change could further threaten this species. Amplified fragment length polymorphisms (AFLPs) were assayed in 10 populations in BC and 1 in Washington to obtain baseline population genetic and mating system data. We found that genetic diversity within populations was low (mean H = 0.094) for a woody perennial. Genetic variation among populations (FST= 0.15) was comparable to woody perennials on average (0.19). Pairwise kinship coefficients were significantly associated with distance (p < 0.01). The multilocus outcrossing rate for one BC island population was high (tm= 0.97), but inbreeding due to consanguineous matings was also quite high (tm– ts= 0.10). These data can be used to guide conservation strategies and future research priorities for arbutus.

Within‐population spatial genetic structure, neighbourhood size and clonal subrange in the seagrass Cymodocea nodosa

Abstract The extent of clonality within populations strongly influences their spatial genetic structure (SGS), yet this is hardly ever thoroughly analysed. We employed spatial autocorrelation analysis to study effects of sexual and clonal reproduction on dispersal of the dioecious seagrass Cymodocea nodosa. Analyses were performed both at genet level (i.e. excluding clonal repeats) and at ramet level. Clonal structure was characterized by the clonal subrange, a spatial measure of the linear limits where clonality still affects SGS. We show that the clonal subrange is equivalent to the distance where the probability of clonal identity approaches zero. This combined approach was applied to two meadows with different levels of disturbance, Cadiz (stable) and Alfacs (disturbed). Genotypic richness, the proportion of the sample representing distinct genotypes, was moderate (0.38 Cadiz, 0.46 Alfacs) mostly due to dominance of a few clones. Expected heterozygosities were comparable to those found in other clonal plants. SGS analyses at the genet level revealed extremely restricted gene dispersal in Cadiz (Sp = 0.052, a statistic reflecting the decrease of pairwise kinship with distance), the strongest SGS found for seagrass species, comparable only to values for selfing herbaceous land plants. At Cadiz the clonal subrange extended across shorter distances (20-25 m) than in Alfacs (30-35 m). Comparisons of sexual and vegetative components of gene dispersal suggest that, as a dispersal vector within meadows, clonal spread is at least as important as sexual reproduction. The restricted dispersal and SGS pattern in both meadows indicates that the species follows a repeated seedling recruitment strategy.

Fine-Scale Genetic Structure and Gene Dispersal in Centaurea corymbosa (Asteraceae). II. Correlated Paternity Within and Among Sibships

The fine-scale pattern of correlated paternity was characterized within a population of the narrow-endemic model plant species, Centaurea corymbosa, using microsatellites and natural progeny arrays. We used classical approaches to assess correlated mating within sibships and developed a new method based on pairwise kinship coefficients to assess correlated paternity within and among sibships in a spatio-temporal perspective. We also performed numerical simulations to assess the relative significance of different mechanisms promoting correlated paternity and to compare the statistical properties of different estimators of correlated paternity. Our new approach proved very informative to assess which factors contributed most to correlated paternity and presented good statistical properties. Within progeny arrays, we found that about one-fifth of offspring pairs were full-sibs. This level of correlated mating did not result from correlated pollen dispersal events (i.e., pollen codispersion) but rather from limited mate availability, the latter being due to limited pollen dispersal distances, the heterogeneity of pollen production among plants, phenological heterogeneity and, according to simulations, the self-incompatibility system. We point out the close connection between correlated paternity and the "TwoGener" approach recently developed to infer pollen dispersal and discuss the conditions to be met when applying the latter.

New insights from fine‐scale spatial genetic structure analyses in plant populations

Many empirical studies have assessed fine-scale spatial genetic structure (SGS), i.e. the nonrandom spatial distribution of genotypes, within plant populations using genetic markers and spatial autocorrelation techniques. These studies mostly provided qualitative descriptions of SGS, rendering quantitative comparisons among studies difficult. The theory of isolation by distance can predict the pattern of SGS under limited gene dispersal, suggesting new approaches, based on the relationship between pairwise relatedness coefficients and the spatial distance between individuals, to quantify SGS and infer gene dispersal parameters. Here we review the theory underlying such methods and discuss issues about their application to plant populations, such as the choice of the relatedness statistics, the sampling scheme to adopt, the procedure to test SGS, and the interpretation of spatial autocorrelograms. We propose to quantify SGS by an 'Sp' statistic primarily dependent upon the rate of decrease of pairwise kinship coefficients between individuals with the logarithm of the distance in two dimensions. Under certain conditions, this statistic estimates the reciprocal of the neighbourhood size. Reanalysing data from, mostly, published studies, the Sp statistic was assessed for 47 plant species. It was found to be significantly related to the mating system (higher in selfing species) and to the life form (higher in herbs than trees), as well as to the population density (higher under low density). We discuss the necessity for comparing SGS with direct estimates of gene dispersal distances, and show how the approach presented can be extended to assess (i) the level of biparental inbreeding, and (ii) the kurtosis of the gene dispersal distribution.

Fine scale genetic structure in a wild soybean (Glycine soja) population and the implications for conservation.

• The fine-scale spatial distribution of a wild soybean (Glycine soja Sieb. et Zucc) population was measured using kinship coefficients and spatial autocorrelation coefficients to reveal fine-scale spatial patterns of genetic structure for guiding collection in ex situ conservation. • Individual families (100) of a single population were sampled in Jiangwan (Shanghai, China) and their specific locations mapped. Individual plants were genotyped for 331 ISSR loci from 15 selective ISSR primers; 45% were demonstrated to be polymorphic. • Pairwise kinship coefficients between individuals were calculated, and then regressed based on the spatial distance between individuals. The results indicated a close genetic relationship among individual families in the patch with a diameter of c. 30m. Spatial autocorrelation analysis suggested that there were significant clusters of genes within the population; c. 81.4% of the loci were found to be positively correlated in the first two distance classes (0-10m). • The results suggest, for conservation genetics, sampling within a wild soybean population should be conducted at 10 m intervals across the entire population to optimize the genetic diversity in collections taken from that population.

The inheritance of hand osteoarthritis in Iceland.

To assess the contribution of genetics to hand osteoarthritis (HOA) and its subsets in the Icelandic population. A list of 2,919 HOA patients, constituting 1% of the Icelandic population, was compiled through nationwide sources. This patient list was cross-referenced with a comprehensive Icelandic genealogy database, enabling the use of algorithms to assess familiality of HOA. Two methods were used: the average pairwise kinship coefficient (KC) of the patients, and the relative risk (RR) of HOA in relatives of patients. In each case, the results were compared with 1,000 control sets of similar composition with regard to number, age, and sex, generated from the genealogy database. The KC for patients was significantly higher than for the control sets and was proportional to the degree of both interphalangeal (IP) and thumb base (first carpometacarpal [CMC] joint) involvement. The RR of HOA in sisters of women in the study was 2.0 (P < 0.001), while the RR in spouses was not significantly different from that in controls. The RR increased with the severity of the disease. Thus, sisters of women with severe IP HOA had an RR of 5.0 and sisters of those with severe first CMC involvement had an RR of 6.9. The increased risk also extended beyond the nuclear family, with significantly increased risk in cousins. Patients seeking medical services for HOA are more related to each other than matched controls, supporting the role of a genetic component in the disease. The genetic influence in both IP and first CMC HOA appears to be similar and increases with increasing severity of the disease.

QUANTIFYING GENE FLOW FROM SPATIAL GENETIC STRUCTURE DATA IN A METAPOPULATION OF CHAMAECRISTA FASCICULATA (LEGUMINOSAE)

An extensive allozyme survey was conducted within a natural “meta” population of the native North American annual legume, Chamaecrista fasciculata (Leguminosae) to quantify genetic structure at different spatial scales. Gene flow was then estimated by a recently developed indirect method based on a continuous population model, using pairwise kinship coefficients between individuals. The indirect estimates of gene flow, quantified in terms of neighborhood size, with an average value on the order of 150 individuals, were concordant among different spatial scales (subpopulation, population, metapopulation). This gene-flow value lies within the range of direct estimates previously documented from observations of pollen and seed dispersal for the same metapopulation. Monte Carlo simulations using the direct measures of gene flow as parameters further demonstrated that the observed spatial pattern of allozyme variation was congruent with a model of isolation by distance. Combining previously published estimates of pollen dispersal distances with kinship coefficients from this study, we quantified biparental inbreeding relative to either a single subpopulation or the whole metapopulation. At the level of a neighborhood, little biparental inbreeding was observed and most departure from Hardy-Weinberg genotypic proportions was explained by self-fertilization, whereas both selfing and biparental inbreeding contributed to nonrandom mating at the metapopulation level. Gene flow was also estimated from indirect methods based on a discontinuous population structure model. We discuss these results with respect to the effect of a patchy population structure on estimation of gene flow.

Prevalence and inheritance of hip osteoarthritis in Iceland

The purpose of this study was to: (I) assess the prevalence of hip osteoarthritis (OA) in Iceland and compare it with that in Southern Scandinavia, (II) determine the incidence of total hip replacement (THR) for primary OA in Iceland, (III) compare two different methods for defining radiographic hip OA, (IV) assess in a population-wide study in Iceland the genetic contribution t hip OA leading to THR, and (V) perform a genome-wide scan of a large Icelandic family to identify a chromosomal susceptibility locus for hip OA leading to THR. many Icelandic patients with hip OA have been well aware that this disease "goes in the family". by examining a large proportion of all Icelandic colon radiographs taken 1990-1996 the prevalence of radiographic hip OA in Iceland was found to be at least five-fold higher compared to Swedish and Danish studies that have used the same methods. A comparison of two methods for estimating hip OA from colon radiographs showed that a simple quantitative method of measuring joint space was more reliable than a qualitative method. The age-standardized incidence of THR for primary hip OA in Iceland between 1982 and 1996 was estimated and found to be about 50 percent higher than for Sweden. The higher Icelandic prevalence of hip OA may explain most of this difference. To investigate the contribution of heritability to hip OA leading to THR, information from two population-wide database in Iceland was combined: A national registry of THR between 1972 and 1996, and a genealogy database of all Icelandic genealogy records for the last 11 centuries made available by deCode Genetics. The genetic contribution to THR for OA was assessed by (a) identifying familial clusters of THR for OA, (b) applying the minimum founder test (MFT) to estimate the minimum number of ancestors to account for all patients with THR for OA, compared to the average number of founders for control lists, (c) calculating an average pairwise kinship coefficient (KC) for the patient and control lists, (d) estimating the relative risk (RR) for relatives of patients with THR for OA. A large number of familial clusters of patients with THR for OA were identified. MFT showed that OA patients descended from fewer founders than the control groups. The average pairwise KC among patients with OA was greater than in the population. RR for siblings of THR for OA patients was 3.05 (2.52, 3.10). Icelandic patients with THR for OA are thus significantly more related to each other than are matched controls. These findings support a significant genetic contribution to a common form of OA and encourages the search for genes conferring an increased susceptibility to OA. New techniques now make it possible to search the whole human genome for chromosomal susceptibility loci associating with OA. A genome wide scan was done to identify susceptibility loci for hip OA leading to THR, using DNA from a large Icelandic family with a very high prevalence of primary hip OA. A genome locus with a lod score of 2.58 was identified on chromosome 16p. A similar locus has been reported on from England. This is the first instance where what may be the same susceptibility locus for OA is independently described in two different populations with hip OA. We have identified other families with hip OA which link to the studied family and are continuing an expanded genome-wide scan. Continued studies of the kind outlined here will clarify the complex genetic background of OA and identify genetic variation associated with the disease. In addition to improving our understanding of the pathogenesis of OA and identifying new molecular targets for treatment, this will allow a better insight into the interactions between genetic background and environmental factors that initiate and drive OA.

Prevalence and inheritance of hip osteoarthritis in Iceland

The purpose of this study was to: (I) assess the prevalence of hip osteoarthritis (OA) in Iceland and compare it with that in Southern Scandinavia, (II) determine the incidence of total hip replacement (THR) for primary OA in Iceland, (III) compare two different methods for defining radiographic hip OA, (IV) assess in a population-wide study in Iceland the genetic contribution t hip OA leading to THR, and (V) perform a genome-wide scan of a large Icelandic family to identify a chromosomal susceptibility locus for hip OA leading to THR. many Icelandic patients with hip OA have been well aware that this disease "goes in the family". by examining a large proportion of all Icelandic colon radiographs taken 1990-1996 the prevalence of radiographic hip OA in Iceland was found to be at least five-fold higher compared to Swedish and Danish studies that have used the same methods. A comparison of two methods for estimating hip OA from colon radiographs showed that a simple quantitative method of measuring joint space was more reliable than a qualitative method. The age-standardized incidence of THR for primary hip OA in Iceland between 1982 and 1996 was estimated and found to be about 50 percent higher than for Sweden. The higher Icelandic prevalence of hip OA may explain most of this difference. To investigate the contribution of heritability to hip OA leading to THR, information from two population-wide database in Iceland was combined: A national registry of THR between 1972 and 1996, and a genealogy database of all Icelandic genealogy records for the last 11 centuries made available by deCode Genetics. The genetic contribution to THR for OA was assessed by (a) identifying familial clusters of THR for OA, (b) applying the minimum founder test (MFT) to estimate the minimum number of ancestors to account for all patients with THR for OA, compared to the average number of founders for control lists, (c) calculating an average pairwise kinship coefficient (KC) for the patient and control lists, (d) estimating the relative risk (RR) for relatives of patients with THR for OA. A large number of familial clusters of patients with THR for OA were identified. MFT showed that OA patients descended from fewer founders than the control groups. The average pairwise KC among patients with OA was greater than in the population. RR for siblings of THR for OA patients was 3.05 (2.52, 3.10). Icelandic patients with THR for OA are thus significantly more related to each other than are matched controls. These findings support a significant genetic contribution to a common form of OA and encourages the search for genes conferring an increased susceptibility to OA. New techniques now make it possible to search the whole human genome for chromosomal susceptibility loci associating with OA. A genome wide scan was done to identify susceptibility loci for hip OA leading to THR, using DNA from a large Icelandic family with a very high prevalence of primary hip OA. A genome locus with a lod score of 2.58 was identified on chromosome 16p. A similar locus has been reported on from England. This is the first instance where what may be the same susceptibility locus for OA is independently described in two different populations with hip OA. We have identified other families with hip OA which link to the studied family and are continuing an expanded genome-wide scan. Continued studies of the kind outlined here will clarify the complex genetic background of OA and identify genetic variation associated with the disease. In addition to improving our understanding of the pathogenesis of OA and identifying new molecular targets for treatment, this will allow a better insight into the interactions between genetic background and environmental factors that initiate and drive OA.

The inheritance of hip osteoarthritis in Iceland.

To assess, in a population-wide study in Iceland, the genetic contribution to hip osteoarthritis (OA) leading to total hip replacement (THR). Information from 2 population-based databases in Iceland was combined: a national registry of all THRs performed between 1972 and 1996, and a genealogy database of all available Icelandic genealogy records for the last 11 centuries. A genetic contribution to THR for OA was assessed by 1) identifying familial clusters of OA patients with THR, 2) applying the minimum founder test (MFT) to estimate the minimum number of ancestors ("founders") that would account for the genealogy of all 2,713 patients with THR for OA, compared with the average number of founders for control lists, 3) calculating an average pairwise kinship coefficient (KC) for the patient list and control lists, and 4) estimating the relative risk (RR) for THR among relatives of OA patients who have undergone the procedure. One thousand matched control lists, each the same size as the patient list, were created using the genealogy database. A large number of familial clusters of patients with THR for OA were identified. The MFT showed that OA patients descended from fewer founders than did subjects in the control groups (P < 0.001). The average pairwise KC among patients with OA was greater than in the control population (P < 0.001). The RR for THR among siblings of OA patients was 3.05 (95% confidence interval 2.52-3.10). This population-based study shows that Icelandic patients with hip replacement for OA are significantly more related to each other than are matched controls drawn from the Icelandic population. These findings support a significant genetic contribution to a common form of OA and encourage the search for genes conferring an increased susceptibility to OA.

A Linear Time Algorithm for Calculation of Multiple Pairwise Kinship Coefficients and the Genetic Index of Familiality

We present an algorithm that calculates the mean kinship coefficient between all pairs of individuals in a given set relative to a known pedigree. The time required to make the calculations is proportional to the size of the set. The constant of proportionality is determined by the number of generations that the pedigree spans. A simple alteration to the algorithm allows efficient calculation of all pairwise kinship coefficients within a set of individuals.