Prospects for whole genome linkage disequilibrium mapping in domestic dog breeds

Linkage disequilibrium (LD) mapping is commonly used as a fine mapping tool in human genome mapping and has been used with some success for initial disease gene isolation in certain isolated inbred human populations. An understanding of the population history of domestic dog breeds suggests that LID mapping could be routinely utilized in this species for initial genome-wide scans. Such an approach offers significant advantages over traditional linkage analysis. Here, we demonstrate, using canine copper toxicosis in the Bedlington terrier as the model, that LID mapping could be reasonably expected to be a useful strategy in low-resolution, genome-wide scans in pure-bred dogs. Significant LID was demonstrated over distances up to 33.3 cM. It is very unlikely, for a number of reasons discussed, that this result could be extrapolated to the rest of the genome. It is, however, consistent with the expectation given the population structure of canine breeds and, in this breed at least, with the hypothesis that it may be possible to utilize LID in a genome-wide scan. In this study, LD mapping confirmed the location of the copper toxicosis in Bedlington terrier gene (CT-BT) and was able to do so in a population that was refractory to traditional linkage analysis.

New concepts for distinguishing the hidden patterns of linkage disequilibrium which underlie association between genotypes and complex phenotypes

A disappointing feature of conventional methods for detecting association between DNA variation and a phenotype of interest is that they tell us little about the hidden pattern of linkage disequilibrium (LD) with the functional variant that is actually responsible for the association. This limitation applies to case-control studies and also to the transmission/disequilibrium test (TDT) and other family-based association methods. Here we present a fresh perspective on genetic association based on two novel concepts called 'LD squares' and 'equi-risk alleles'. These describe and characterize the different patterns of gametic LD which underlie genetic association. These concepts lead to a general principle - the Equi-Risk Allele Segregation Principle - which captures the way in which underlying LD patterns affect the transmission patterns of genetic variants associated with a phenotype. This provides a basis for distinguishing the hidden LD patterns and might help to locate the functional variants responsible for the association.

Genomewide significant linkage to migrainous Hheadache on chromosome 5q21

Familial typical migraine is a common, complex disorder that shows strong familial aggregation. Using latent-class analysis (LCA), we identified subgroups of people with migraine/severe headache in a community sample of 12,245 Australian twins (60% female), drawn from two cohorts of individuals aged 23-90 years who completed an interview based on International Headache Society criteria. We report results from genomewide linkage analyses involving 756 twin families containing a total of 790 independent sib pairs ( 130 affected concordant, 324 discordant, and 336 unaffected concordant for LCA-derived migraine). Quantitative-trait linkage analysis produced evidence of significant linkage on chromosome 5q21 and suggestive linkage on chromosomes 8, 10, and 13. In addition, we replicated previously reported typical-migraine susceptibility loci on chromosomes 6p12.2-p21.1 and 1q21-q23, the latter being within 3 cM of the rare autosomal dominant familial hemiplegic migraine gene (ATP1A2), a finding which potentially implicates ATP1A2 in familial typical migraine for the first time. Linkage analyses of individual migraine symptoms for our six most interesting chromosomes provide tantalizing hints of the phenotypic and genetic complexity of migraine. Specifically, the chromosome 1 locus is most associated with phonophobia; the chromosome 5 peak is predominantly associated with pulsating headache; the chromosome 6 locus is associated with activity-prohibiting headache and photophobia; the chromosome 8 locus is associated with nausea/vomiting and moderate/severe headache; the chromosome 10 peak is most associated with phonophobia and photophobia; and the chromosome 13 peak is completely due to association with photophobia. These results will prove to be invaluable in the design and analysis of future linkage and linkage disequilibrium studies of migraine.

Association of a glucocorticoid receptor gene marker with human essential hypertension.

Recently, a bi-allelic polymorphism in the glucocorticoid receptor gene (GRL) has been shown to be associated with individuals at high risk of developing hypertension and accumulation of abdominal visceral fat, a known risk factor for cardiovascular disease. The evaluate the role of GRL in essential hypertension and obesity, case-control studies were conducted using 88 hypertensive, 123 normotensive, 150 lean and 94 obese subjects. Genotypes for a highly polymorphic microsatellite marker (D5S207) located within 200 kb of the glucocorticoid receptor gene, were determined by PCR. Allele frequencies between hypertensive and normotensive groups were significantly (P = 0.0005) different whereas no significant differences were observed between lean and obese populations. In conclusion, the results suggest that the glucocorticoid receptor gene or perhaps another gene located in close proximity and in linkage disequilibrium with D5S207, is involved in hypertension development

A high-resolution genetic map of the familial Mediterranean fever candidate region allows identification of haplotype-sharing among ethnic groups

Familial Mediterranean fever (FMF) is a recessive disorder of inflammation caused by mutations in a gene (designated MEFV) on chromosome 16p13.3, We have recently constructed a 1-Mb cosmid contig that includes the FMF critical region. Here we show genotype data for 12 markers from our physical map, including 5 newly identified microsatellites, in FMF families. Intrafamilial recombinations placed MEFV in the similar to 285 kb between D16S468/D16S3070 and D16S3376. We observed significant linkage disequilibrium in the North African Jewish population, and historical recombinants in the founder haplotype placed MEFV between D16S3082 and D16S3373 (similar to 200 kb). In smaller panels of Iraqi Jewish, Arab, and Armenian families, there were significant allelic associations only for D16S3370 and D16S2617 among the Armenians. A sizable minority of Iraqi Jewish and Armenian carrier chromosomes appeared to be derived from the North African Jewish ancestral haplotype. We observed a unique FMF haplotype common to Iraqi Jews, Arabs, and Armenians and two other haplotypes restricted to either the Iraqi Jewish or the Armenian population. These data support the view that a few major mutations account for a large percentage of the cases of FMF and suggest that same of these mutations arose before the affected Middle Eastern populations diverged from one another. (C) 1997 Academic Press.

HLA class II antigens positively and negatively associated with hepatosplenic schistosomiasis in a Chinese population

To identify possible associations between host genetic factors and the onset of liver fibrosis following Schistosoma japonicum infection, the major histocompatibility class II alleles of 84 individuals living on an island (Jishan) endemic for schistosomiasis japonica in the Poyang Lake Region of Southern China were determined. Forty patients exhibiting advanced schistosomiasis, characterised by extensive liver fibrosis, and 44 age and sex-matched control subjects were assessed for the class II haplotypes HLA-DRBI and HLA-DQB1. Two HLA-DRB1 alleles, HLA-DRB1*0901 (P = 0.012) and *1302 (P = 0.039), and two HLA-DQB1 alleles, HLA-DQB1*0303 (P = 0.012) and *0609 (P = 0.037), were found to be significantly associated with susceptibility to fibrosis. These associated DRB1 and DQB1 alleles are in very strong linkage disequilibrium, with DRB1*0901-DQB1*0303 and DRB1*1302-DQB1*0609 found as: common haplotypes in this population. In contrast, the alleles HLA-DRB1*1501 (P = 0.025) and HLA-DQB 1*0601 (P = 0.022) were found to be associated with resistance to hepatosplenic disease. Moreover, the alleles DQB1*0303 and DRB1*0901 did not increase susceptibility in the presence of DQB1*0601, indicating that DQB1*0601 is dominant over DQB1*0303 and DRB1*0901. The study has thus identified both positive and negative associations between HLA class II alleles and the risk of individuals developing moderate to severe liver fibrosis following schistosome infection. (C) 2001 Australian Society for Parasitology Inc. Published by Elsevier Science Ltd. All rights reserved.

The human melanocortin-1 receptor locus: analysis of transcription unit, locus polymorphism and haplotype evolution

The complete sequence of the MCIR locus has been assembled, the coding region of the gene is intronless and placed within a 12 kb region flanked by the NULP1 and TUBB4 genes. The immediate promoter region has an E-box site with homology to the M-box consensus known to bind the microphthalmia transcription factor (MITF), however, promoter deletion analysis and transactivation studies have failed to show activation through this element by MITF. Polymorphism within the coding region, immediate 5' promoter region and a variable number tandem repeat (VNTR) minisatellite within the locus have been examined in a collection of Caucasian families and African individuals. Haplotype analysis shows linkage disequilibrium between the VNTR and MCIR coding region red hair variant alleles which can be used to estimate the age of these missense changes. Assuming a mean VNTR mutation rate of 1% and a star phylogeny, we estimate the Arg151Cys variant arose 7500 years before the present day, suggesting these variants may have arisen in the Caucasian population more recently than previously thought. (C) 2001 Published by Elsevier Science B.V.

The progesterone receptor exon 4 Va1660Leu G/T polymorphism and risk of breast cancer in Australian women

An Alu insertion polymorphism of the progesterone receptor (PR) was reported recently to be associated with a reduced risk of breast cancer, with risks of 0.8- and 0.3-fold associated with the heterozygote and homozygote genotypes, respectively. This intronic variant is considered to be in linkage disequilibrium with an exon 4 hinge region G to T Val660Leu polymorphism. We investigated whether the exon 4 PR polymorphism was associated with breast cancer in Australian women, using a population-based study of 1452 cases and 793 controls, half of whom were

Repeated blood pressure measurements in a sample of Swedish twins: heritabilities and associations with polymorphisms in the renin-angiotensin-aldosterone system Iliadou, Anastasia, Lichtenstein, Paul, Morgenstern, Ralf, Forsberg, Lena, Svensson, Richard, de Faire, Ulf, Martin, Nicholas, Pedersen, Nancy

Background Twin and family studies have shown that genetic effects explain a relatively high amount of the phenotypic variation in blood pressure. However, many studies have not been able to replicate findings of association between specific polymorphisms and diastolic and systolic blood pressure. Methods In a structural equation-modelling framework the authors investigated longitudinal changes in repeated measures of blood pressures in a sample of 298 like-sexed twin pairs from the population-based Swedish Twin Registry. Also examined was the association between blood pressure and polymorphisms in the angiotensin-I converting enzyme and the angiotensin 11 receptor type 1 with the 'Fulker' test Both linkage and association were tested simultaneously revealing whether the polymorphism is a Quantitative Trait Locus (QTL) or in linkage disequilibrium with the QTL. Results Genetic influences explained up to 46% of the phenotypic variance in diastolic and 63% of the phenotypic variance in systolic blood pressure. Genetic influences were stable over time and contributed up to 78% of the phenotypic correlation in both diastolic and systolic blood pressure. Non-shared environmental effects were characterised by time specific influences and little transmission from one time point to the next. There was no significant linkage and association between the polymorphisms and blood pressure. Conclusions There is a considerable genetic stability in both diastolic and systolic blood pressure for a 6-year period of time in adult life. Non-shared environmental influences have a small long-term effect Although associations with the polymorphisms could not be replicated, results should be interpreted with caution due to power considerations. (C) 2002 Lippincott Williams Wilkins.

The impact of genotyping error on haplotype reconstruction and frequency estimation

The choice of genotyping families vs unrelated individuals is a critical factor in any large-scale linkage disequilibrium (LD) study. The use of unrelated individuals for such studies is promising, but in contrast to family designs, unrelated samples do not facilitate detection of genotyping errors, which have been shown to be of great importance for LD and linkage studies and may be even more important in genotyping collaborations across laboratories. Here we employ some of the most commonly-used analysis methods to examine the relative accuracy of haplotype estimation using families vs unrelateds in the presence of genotyping error. The results suggest that even slight amounts of genotyping error can significantly decrease haplotype frequency and reconstruction accuracy, that the ability to detect such errors in large families is essential when the number/complexity of haplotypes is high (low LD/common alleles). In contrast, in situations of low haplotype complexity (high LD and/or many rare alleles) unrelated individuals offer such a high degree of accuracy that there is little reason for less efficient family designs. Moreover, parent-child trios, which comprise the most popular family design and the most efficient in terms of the number of founder chromosomes per genotype but which contain little information for error detection, offer little or no gain over unrelated samples in nearly all cases, and thus do not seem a useful sampling compromise between unrelated individuals and large families. The implications of these results are discussed in the context of large-scale LD mapping projects such as the proposed genome-wide haplotype map.

Evolution of additive and nonadditive genetic variance in development time along a cline in Drosophila serrata

Latitudinal clines provide natural systems that may allow the effect of natural selection on the genetic variance to be determined. Ten clinal populations of Drosophila serrata collected from the eastern coast of Australia were used to examine clinal patterns in the trait mean and genetic variance of the life-history trait egg-to-adult development time. Development time significantly lengthened from tropical areas to temperate areas. The additive genetic variance for development time in each population was not associated with latitude but was associated with the population mean development time. Additive genetic variance tended to be larger in populations with more extreme development times and appeared to be consistent with allele frequency change. In contrast, the nonadditive genetic variance was not associated with the population mean but was associated with latitude. Levels of nonadditive genetic variance were greatest in the region of the cline where the gradient in the change in mean was greatest, consistent with Barton's (1999) conjecture that the generation of linkage disequilibrium may become an important component of the genetic variance in systems with a spatially varying optimum.

Dizygotic twinning is not associated with methylenetetrahydrofolate reductase haplotypes

Background: Folate metabolism is critical to embryonic development, influencing neural tube defects (NTD) and recurrent early pregnancy loss. Polymorphisms in 5,10-methylenetetrahydrofolate reductase (MTHFR) have been associated with dizygotic (DZ) twinning through pregnancy loss. Methods: The C677T and A1298C polymorphisms in MTHFR were genotyped in 258 Australasian families (1016 individuals) and 118 Dutch families (462 individuals) of mothers of DZ twins and a population sample of 462 adolescent twin families (1861 individuals). Haplotypes were constructed from the alleles, and transmission of the MTHFR haplotypes to mothers of DZ twins and from parents to twins in the adolescent twin families analysed. Results: The C677T and A1298C were common in all three populations (frequencies > 0.29). There was strong linkage disequilibrium (D'=1) between the variants, showing that specific combinations of alleles (haplotypes) were transmitted together. Three haplotypes accounted for nearly all the variation. There was no evidence of any association between MTHFR genotype and twinning in mothers of twins, or of the loss of specific MTHFR genotypes during twin pregnancies. Conclusions: It is concluded that variation in twinning frequency is not associated with MTHFR genotype.

Polymorphisms in the 5 '-untranslated region of the human serotonin receptor 1B (HTR1B) gene affect gene expression

We present evidence of complex balancing regulation of HTR1B transcription by common polymorphisms in its promoter. Computational analysis of the HTR1B gene predicted that a 50 segment, spanning common DNA sequence variations, T-261G, A-161T, and -182INS/DEL-181, contained a putative functional promoter. Using a secreted alkaline phosphatase (SEAP) reporter gene system, we found that the haplotype -261G_-182INS-181_A-161 enhanced transcriptional activity 2.3-fold compared with the haplotype T-261_-182INS-181_A-161. Conversely, -161T reversed this, and the net effect when -261G and -161T were in the same haplotype (-261G_-182INS-181_-161T) was equivalent to the major haplotype (T-261_-182INS-181_A-161). Electrophoretic mobility shift experiments showed that -261G and -161T modify the binding of transcription factors (TFs): -261G generates a new AP2 binding site, while alleles A-161 and -161T exhibit different binding characteristics to AP1. T-261G and A-161T were found to be in linkage disequilibrium (LD) with G861C in a European ancestry population. Interestingly, G861C has been reported to be associated with several psychiatric disorders. Our results indicate that HTR1B is the target of substantial transcriptional genetic regulation by common haplotypes, which are in LD with the HTR1B single-nucleotide polymorphism (SNP) most commonly used in association studies.

Genetic assignment methods for the direct, real-time estimation of migration rate: a simulation-based exploration of accuracy and power

Genetic assignment methods use genotype likelihoods to draw inference about where individuals were or were not born, potentially allowing direct, real-time estimates of dispersal. We used simulated data sets to test the power and accuracy of Monte Carlo resampling methods in generating statistical thresholds for identifying F-0 immigrants in populations with ongoing gene flow, and hence for providing direct, real-time estimates of migration rates. The identification of accurate critical values required that resampling methods preserved the linkage disequilibrium deriving from recent generations of immigrants and reflected the sampling variance present in the data set being analysed. A novel Monte Carlo resampling method taking into account these aspects was proposed and its efficiency was evaluated. Power and error were relatively insensitive to the frequency assumed for missing alleles. Power to identify F-0 immigrants was improved by using large sample size (up to about 50 individuals) and by sampling all populations from which migrants may have originated. A combination of plotting genotype likelihoods and calculating mean genotype likelihood ratios (D-LR) appeared to be an effective way to predict whether F-0 immigrants could be identified for a particular pair of populations using a given set of markers.