Incorporating genotypes of relatives into a test of linkage disequilibrium.

Genetic data from autosomal loci in diploids generally consist of genotype data for which no phase information is available, making it difficult to implement a test of linkage disequilibrium. In this paper, we describe a test of linkage disequilibrium based on an empirical null distribution of the likelihood of a sample. Information on the genotypes of related individuals is explicitly used to help reconstruct the gametic phase of the independent individuals. Simulation studies show that the present approach improves on estimates of linkage disequilibrium gathered from samples of completely independent individuals but only if some offspring are sampled together with their parents. The failure to incorporate some parents sharply decreases the sensitivity and accuracy of the test. Simulations also show that for multiallelic data (more than two alleles) our testing procedure is not as powerful as an exact test based on known haplotype frequencies, owing to the interaction between departure from Hardy-Weinberg equilibrium and linkage disequilibrium.     

AWclust: point-and-click software for non-parametric population structure analysis

Background  

Population structure analysis is important to genetic association studies and evolutionary investigations. Parametric approaches, e.g. STRUCTURE and L-POP, usually assume Hardy-Weinberg equilibrium (HWE) and linkage equilibrium among loci in sample population individuals. However, the assumptions may not hold and allele frequency estimation may not be accurate in some data sets. The improved version of STRUCTURE (version 2.1) can incorporate linkage information among loci but is still sensitive to high background linkage disequilibrium. Nowadays, large-scale single nucleotide polymorphisms (SNPs) are becoming popular in genetic studies. Therefore, it is imperative to have software that makes full use of these genetic data to generate inference even when model assumptions do not hold or allele frequency estimation suffers from high variation.     

Genetic Structure and the Search for Genotype-Phenotype Relationships: An Example from Disequilibrium in the Apo B Gene Region

We analyzed allelic associations (disequilibria) for four restriction fragment length polymorphisms (RFLPs) in the region of the 43-kb Apo B gene in a sample of 233 unrelated individuals from Montreal, Canada, sampled for health. This total sample (T) included 160 individuals of known French Canadian (FC) ancestry. We present a rigorous application of current methodology to these samples, including estimation of type II error probabilities and correlations between markers for estimates of disequilibria. We then consider the utility of these estimates of allelic disequilibria for the interpretation of genotype-phenotype relations. Significant deviations from Hardy-Weinberg equilibrium were not predicted by proximity to other markers in disequilibrium. We found significant quadri-allelic disequilibrium for two marker pairs despite absence of significant deviations from Hardy-Weinberg equilibrium for either marker or tri-allelic disequilibrium, respectively. Altogether these results underscore the complexity of the genotypic structure of the data. A combination of nonevolutionary factors, including sampling for health, small sample size and data exclusion due to methodological constraints of not successfully typing all members of the sample for every RFLP, is a likely explanation for this complexity. These types of factors are common to many RFLP studies. Patterns of composite di-allelic disequilibrium indicated that some RFLP allele pairs may have a longer shared evolutionary history than others and that disequilibrium is not predicted by distance between RFLPs. Type II error probabilities were generally much higher than those for type I errors. Correlations between marker pairs for disequilibria were generally not high. We show from a review of 14 published studies of association between the XbaI RFLP and variation in a total of 15 lipid traits that deviations from Hardy-Weinberg equilibrium can cause substantial differences in the estimation of variability associated with phenotypic differences among marker genotypes relative to Hardy-Weinberg conditions.     

Genome Scan for Predisposing Loci for Distal Interphalangeal Joint Osteoarthritis: Evidence for a Locus on 2q

The genetic contribution to common forms of osteoarthritis (OA) is well established but poorly understood. We performed a genome scan, using 302 markers for loci predisposing to distal interphalangeal joint (DIP) OA. To minimize genetic heterogeneity in our study sample, we identified siblings with a severe, radiologically defined phenotype from the nationwide registers of Finland. In the initial genome scan, linkage analysis in 27 sibships gave a pairwise LOD score (Z) >1.00 with nine of the screening markers. In the second stage, additional markers and family members were genotyped in these chromosomal regions. On 2q12-q13, IL1R1 resulted in Z=2.34 at recombination fraction (theta) 0, allowing a dominant mode of inheritance. Association analysis of markers D2S2264, IL1R1, D2S373, and D2S1789 jointly provided some evidence for a shared haplotype among the affected individuals (P value of.012). Also, multipoint nonparametric linkage analysis yielded a P value of.0001 near the locus IL1R1 and P=.0007 approximately 20 cM telomeric near marker D2S1399, which, in two-point analysis, gave Z=1.48 (straight theta=. 02). This chromosomal region on 2q harbors the interleukin 1 gene cluster and, thus, represents a good candidate region for inflammatory and autoimmune disorders. Three additional chromosomal regions-4q26-q27, 7p15-p21, and Xcen-also provided some evidence for linkage, and further analyses would be justified to clarify their potential involvement in the genetic predisposition to DIP OA.     

A Two-Stage Test for Distinguishing Random,Pseudo-Random and Nonrandom Mating Populations

There is no such an implication that a population in Hardy-Weinberg equilibrium must have undergone random mating. Therefore, it is unequivocal that the usual tests for “Hardy-Weinberg equilibrium” are indeed tests for “random union of gametes” but not for “random mating”. In this paper, utilizing population characteristics expressed in equilibrium state (equilibrium or disequilibrium) and mating behavior (random or nonrandom), a two-stage testing procedure for distinguishing random, pseudo-random and nonrandom mating populations is proposed. At the first stage, a population is tested for Hardy-Weinberg equilibrium. If insignificant result (i.e., in equilibrium) is obtained, then to a second stage the population is further tested for mating behavior. Random mating-pairs data are needed here for analysis instead of random individuals for usual Hardy-Weinberg equilibrium tests. Since distinguishing the three types of mating populations depends on the combined results of two stages, the probability of correct determination of the two-stage tests is discussed by simulation studies.     

Blood groups in the Chueta community (Majorcan Jews)

A sample of 443 Chuetas (descendents of Majorcan Jews) was typed for the ABO, Rh, Lewis, Duffy, MNSs, Kell and P blood groups. Significant deviations from the Hardy-Weinberg equilibrium were observed in the MNSs and Duffy systems with a deficiency of heterozygotes. The gene frequencies were compared with those of the Balearic non-Jewish populations and significant differences were found. Genetic distances and cluster analysis demonstrated that the Chuetas resemble more their neighboring non-Jewish populations than other Jewish populations. Discriminant analyses showed that the Chuetas and other Jews considered in this study do not resemble each other but their host peoples with respect to these markers.     

The interaction between TPH2 promoter haplotypes and clinical-demographic risk factors in suicide victims with major psychoses

Tryptophan hydroxylase isoform 2 (TPH2) is a rate-limiting enzyme in the biosynthesis of serotonin (5-HT) and is predominantly localized in the brain. Previous studies have suggested that there is an association between serotonergic dysfunction in the brain and suicidality. This study was designed to examine whether the -473T > A and -8396G > C polymorphisms of the TPH2 gene may be associated with completed suicide in subjects with major psychoses from the Stanley Foundation Brain Bank sample. TPH2 genotypes were determined in 69 subjects with a diagnosis of schizophrenia or bipolar disorder, among which 22 died by suicide. Genomic DNA was amplified by polymerase chain reaction and typed by automated methods. Both markers were found to be in Hardy-Weinberg equilibrium and in strong linkage disequilibrium. No association with history of suicide was found for either polymorphism. Haplotype analysis with EHAP showed no association between completed suicide and haplotype distribution (chi2 = 1.877; 3 df; P = 0.598). Nor was there any association between suicide and these genetic markers even when clinical-demographic factors were considered as covariates in the haplotype analysis. These findings suggest that these 5' marker haplotypes in the TPH2 gene do not influence suicidal behaviour.     

Genetic evidence for a distinct subtype of schizophrenia characterized by pervasive cognitive deficit

A novel phenotyping strategy in schizophrenia, targeting different neurocognitive domains, neurobehavioral features, and selected personality traits, has allowed us to identify a homogeneous familial subtype of the disease, characterized by pervasive neurocognitive deficit. Our genome scan data indicate that this subtype, which accounts for up to 50% of our sample, has a distinct genetic basis and explains linkage to chromosome 6p24 reported previously. If representative of other populations, the ratio of schizophrenia subtypes observed in our families could have a profound impact on sample heterogeneity and on the power of genetic studies to detect linkage and association. Our proposed abbreviated battery of tests should facilitate phenotype characterization for future genetic analyses and allow a focus on a crisply defined schizophrenia subtype, thus promoting a more informed search for susceptibility genes.     

Linkage disequilibrium testing when linkage phase is unknown

Linkage disequilibrium, the nonrandom association of alleles from different loci, can provide valuable information on the structure of haplotypes in the human genome and is often the basis for evaluating the association of genomic variation with human traits among unrelated subjects. But, linkage phase of genetic markers measured on unrelated subjects is typically unknown, and so measurement of linkage disequilibrium, and testing whether it differs significantly from the null value of zero, requires statistical methods that can account for the ambiguity of unobserved haplotypes. A common method to test whether linkage disequilibrium differs significantly from zero is the likelihood-ratio statistic, which assumes Hardy-Weinberg equilibrium of the marker phenotype proportions. We show, by simulations, that this approach can be grossly biased, with either extremely conservative or liberal type I error rates. In contrast, we use simulations to show that a composite statistic, proposed by Weir and Cockerham, maintains the correct type I error rates, and, when comparisons are appropriate, has similar power as the likelihood-ratio statistic. We extend the composite statistic to allow for more than two alleles per locus, providing a global composite statistic, which is a strong competitor to the usual likelihood-ratio statistic.     

Quantitative-Trait Loci Influencing Body-Mass Index Reside on Chromosomes 7 and 13: The National Heart, Lung, and Blood Institute Family Heart Study

Obesity is a risk factor for many chronic diseases, including glucose intolerance, lipid disorders, hypertension, and coronary heart disease. Even though the body-mass index (BMI) is a heterogeneous phenotype reflecting the amount of fat, lean mass, and body build, several studies have provided evidence of one or two major loci contributing to the variation in this complex trait. We sought to identify loci with potential influence on BMI in the data obtained from National Heart, Lung, and Blood Institute Family Heart Study. Two complementary samples were studied: (a) 1,184 subjects in 317 sibships, with 243 markers typed by the Utah Molecular Genetics Laboratory (UMGL) and (b) 3,027 subjects distributed among 401 three-generation families, with 404 markers typed by the Mammalian Genotyping Service (MGS). A genome scan using a variance-components-based linkage approach was performed for each sample, as well as for the combined sample, in which the markers from each analysis were placed on a common genetic map. There was strong evidence for linkage on chromosome 7q32.3 in each sample: the maximum multipoint LOD scores were 4.7 (P<10-5) at marker GATA43C11 and 3.2 (P=.00007) at marker D7S1804, for the MGS and UMGL samples, respectively. The linkage result is replicated by the consistent evidence from these two complementary subsets. Furthermore, the evidence for linkage was maintained in the combined sample, with a LOD score of 4.9 (P<10-5) for both markers, which map to the same location. This signal is very near the published location for the leptin gene, which is the most prominent candidate gene in this region. For the combined-sample analysis, evidence of linkage was also found on chromosome 13q14, with D13S257 (LOD score 3.2, P=.00006), and other, weaker signals (LOD scores 1.5-1.9) were found on chromosomes 1, 2, 3, 5, 6, 14, and 15.     

Metapopulation genetic structure in the water vole, Arvicola terrestris, in NE Scotland

The genetic structure of nine colonies of the water vole, Arvicola terrestris , in one area of NE Scotland was studied. Non-destructive samples from 478 individuals (mostly immature animals) were typed for 12 microsatellites. Cases of Hardy-Weinberg disequilibrium were frequent within all colonies. The five colonies in the inland part of the study area also showed frequent cases of linkage disequilibrium. All colonies showed high levels of genetic diversity (unbiased H = 0.52-0.74). All five colonies sampled in successive years showed significant annual changes in genetic composition. All colonies showed genetic differentiation from each other, whether measured by average 6, pairwise 9 or pairwise Nei's genetic distance. The spatial pattern of genetic differentiation was consistent with either a stepping-stone model over the whole study area or an island model within the coastal and inland parts and an intervening barrier to gene flow. The study suggested that the genetic structure of colonies of A. terrestris often departs from the equilibrium states assumed by traditional mefhods for the study of gene flow, and that a parentage-based approach would be fruitful.     

Linkage maps of microsatellite DNA markers for the Pacific oyster Crassostrea gigas

We constructed male and female consensus linkage maps for the Pacific oyster Crassostrea gigas, using a total of 102 microsatellite DNA markers typed in 11-day-old larvae from three families. We identified 11 and 12 linkage groups in the male and female consensus maps, respectively. Alignment of these separate maps, however, suggests 10 linkage groups, which agrees with the haploid chromosome number. The male linkage map comprises 88 loci and spans 616.1 cM, while the female map comprises 86 loci and spans 770.5 cM. The male and the female maps share 74 loci; 2 markers remain unlinked. The estimated coverages for the consensus linkage maps are 79% for the male and 70-75% for the female, on the basis of two estimates of genome length. Ninety-five percent of the genome is expected to lie within 16 and 21 cM of markers on the male and female maps, respectively, while 95% of simulated minimum distances to the male and female maps are within 10.1 and 13.6 cM, respectively. Females have significantly more recombination than males, across 118 pairs of linked markers in common to the parents of the three families. Significant differences in recombination and orders of markers are also evident among same-sex parents of different families as well as sibling parents of opposite sex. These observations suggest that polymorphism for chromosomal rearrangements may exist in natural populations, which could have profound implications for interpreting the evolutionary genetics of the oyster. These are the first linkage maps for a bivalve mollusc that use microsatellite DNA markers, which should enable them to be transferred to other families and to be useful for further genetic analyses such as QTL mapping.     

Assignment of canine MSS1 microsatellite markers to chromosomes by linkage data.

Recent advances in mapping the canine genome have led to an increase in the number of linkage studies aimed at dissecting the genetic causes of many hereditary diseases that affect the domestic dog. The first step in developing molecular tools for a whole genome scan was the characterization of a set of microsatellite markers, termed minimal screening set 1 (MSS1), that provided an estimated coverage of 10 cM. A limiting factor in use of the MSS1 is not all of the 172 MSS1 markers have been localized to specific chromosomes. Seventy-five of the markers were positioned on a total of 15 chromosomes with the original publication of the MSS1. The localization based on linkage data of 14 additional MSS1 markers to chromosomes using CRIMAP v. 2.4 to build a linkage map of 113 MSS1 markers that were polymorphic in a kindred of Dalmatians is reported here.     

A Genome Scan among Nigerians Linking Resting Energy Expenditure to Chromosome 16

Energy requirements at rest account for 50% to 75% of total energy expenditure. Interindividual variation in resting energy expenditure (REE) has been studied for potential links to obesity and hypertension. REE is a modestly heritable trait, and yet virtually nothing is known about the genetic factors that might influence the familial patterns. The objectives of this study were to identify the genomic regions showing genetic linkage to REE variation in a Nigerian population. For linkage analysis across the genome, three hundred seventy‐seven microsatellite markers were typed on DNA from 995 individuals in 153 families. A genome scan was performed using a multipoint variance component method. Heritability of REE was 0.30 after adjustment for body size. The strongest linkage signal was detected on chromosome 16 (16q22.3) with a likelihood of odds of 2.96 (p = 0.08). Linkage evidence (likelihood of odds > 1) was detected on another three chromosomal regions, namely 2q12.1, 8q21.2, and 15p11.2.     

The multiple-locus symmetric fertility model

Selection due to variation in the fecundity among matings of genotypes with respect to many loci each with two alleles is studied. The fitness of a mating depends only on the genotypic distinction between homozygote and heterozygote at each locus in the two individuals, and differences among loci are allowed. This symmetric fertility model is therefore a generalization of the multiple-locus symmetric viability model. The phenomena seen in the two-locus symmetric fertility model generalize—e.g., the possibility of joint stability of equilibria with linkage equilibrium and with linkage disequilibrium, and the existence of different types of totally polymorphic equilibria with the gametic proportions in linkage equilibrium. The central equilibrium with genotypic frequencies in Hardy-Weinberg proportions and gametic frequencies in Robbins proportions exists for all symmetric fertility models. For some symmetric fertility regimes additional equilibria exist with gametic frequencies in linkage equilibrium and with genotypic frequencies in Hardy-Weinberg proportions at all except one locus. These equilibria may exist in the dioecious symmetric viability model, and then they will be locally stable. For free recombination the stable equilibria show linkage equilibrium, but several of these with different numbers of polymorphic loci may be stable simultaneously.     

Replication of autism linkage: fine-mapping peak at 17q21

Autism is a heritable but genetically complex disorder characterized by deficits in language and in reciprocal social interactions, combined with repetitive and stereotypic behaviors. As with many genetically complex disorders, numerous genome scans reveal inconsistent results. A genome scan of 345 families from the Autism Genetic Resource Exchange (AGRE) (AGRE_1), gave the strongest evidence of linkage at 17q11-17q21 in families with no affected females. Here, we report a full-genome scan of an independent sample of 91 AGRE families with 109 affected sibling pairs (AGRE_2) that also shows the strongest evidence of linkage to 17q11-17q21 in families with no affected females. Taken together, these samples provide a replication of linkage to this chromosome region that is, to our knowledge, the first such replication in autism. Fine mapping at 2-centimorgan (cM) intervals in the combined sample of families with no affected females reveals a linkage peak at 66.85 cM, which places this locus at 17q21.     

Isolation and characterization of 11 tetranucleotide microsatellite loci in the Egyptian mongoose (Herpestes ichneumon)

We report the isolation of 11 polymorphic tetranucleotide microsatellite loci in the Egyptian mongoose (Herpestes ichneumon). In a sample of 27 individuals, we observed between 4 and 7 alleles per locus and their observed and expected heterozygosities ranged from 0.37 to 0.85 and from 0.44 to 0.79, respectively. All genotypic frequencies conformed to Hardy-Weinberg equilibrium expectations and there were no instances of linkage disequilibrium detected between pairs of loci.     

The effect of close relatives on unsupervised Bayesian clustering algorithms in population genetic structure analysis

The inference of population genetic structures is essential in many research areas in population genetics, conservation biology and evolutionary biology. Recently, unsupervised Bayesian clustering algorithms have been developed to detect a hidden population structure from genotypic data, assuming among others that individuals taken from the population are unrelated. Under this assumption, markers in a sample taken from a subpopulation can be considered to be in Hardy-Weinberg and linkage equilibrium. However, close relatives might be sampled from the same subpopulation, and consequently, might cause Hardy-Weinberg and linkage disequilibrium and thus bias a population genetic structure analysis. In this study, we used simulated and real data to investigate the impact of close relatives in a sample on Bayesian population structure analysis. We also showed that, when close relatives were identified by a pedigree reconstruction approach and removed, the accuracy of a population genetic structure analysis can be greatly improved. The results indicate that unsupervised Bayesian clustering algorithms cannot be used blindly to detect genetic structure in a sample with closely related individuals. Rather, when closely related individuals are suspected to be frequent in a sample, these individuals should be first identified and removed before conducting a population structure analysis.     

Asymptotic properties of affected-sib-pair linkage analysis.

The likelihood-ratio method for affected-sib-pair analysis, introduced by Risch, is a powerful method for detecting linkage when the marker is not perfectly polymorphic, as is often the case. The power of this method can be improved by restricting maximization to the set of possible haplotype-sharing probabilities--denoted the "possible triangle" method. The asymptotic distributions of the resulting distributions are derived, enabling test criteria to be found for any required test size (i.e., the probability of falsely detecting linkage when none exists) and enabling p values to be assigned to results. The criteria were found to be approximately constant when the PIC of the marker varies, making them applicable to any marker. The asymptotic power approximations were used to investigate the relative performance of pairs with typed parents, relative to those without, by comparing the sample sizes necessary for a given power. Under certain circumstances, typing the parents proved to be inefficient, even when PIC was low.     

Linkage disequilibrium in synthetic varieties of perennial ryegrass

Synthetic varieties obtained after three to four panmictic generations are variable, not structured and so can be used for association studies. The pattern of linkage disequilibrium (LD) decay determines whether a genome scan or a candidate gene approach can be used for an association study between genotype and phenotype. Our goal was to evaluate the effect of the number of parents used to build the synthetic varieties on the pattern of LD decay. LD was investigated in the gibberelic acid insensitive gene (GAI) region in three synthetic varieties of perennial ryegrass (Lolium perenne L.) chosen for their contrasted number of parents in the initial polycrosses. Results were compared with those obtained from a core collection. STS and SSR markers were used to evaluate variation, structuration and LD in each variety. As expected, the varieties variability increased with the number of parents almost up to the core collection variability. No structuration was observed in the varieties. Significant LDs were observed up to 1.6 Mb in a variety originated from six related parents and not above 174 kb in a variety originated from 336 parents. These results suggest that a candidate gene approach can be used when varieties have a large number of parents and a genome scan approach can be envisaged in specific regions when varieties have a low number of parents. Nevertheless, we strongly recommend to estimate the pattern of LD decay in the population and in the genomic region studied before performing an association study.