Self-reported walking pace, polygenic risk scores and risk of coronary artery disease in UK biobank

BACKGROUND AND AIMS

Both polygenic risk scores (PGS) and self-reported walking pace have been shown to predict cardiovascular disease; whether combining both factors produces greater risk differentiation is, however, unknown.

METHODS AND RESULTS

We estimated the 10-year absolute risk of coronary artery disease (CAD), adjusted for traditional risk factors, and the C-index across nine PGS and self-reported walking pace in UK Biobank study participants between Mar/2006-Feb/2021. In 380,693 individuals (54.8% women), over a median (5th, 95th percentile) of 11.9 (8.3, 13.4) years, 2,603 (1.2%) CAD events occurred in women and 8,259 (4.8%) in men. Both walking pace and genetic risk were strongly associated with CAD. The absolute 10-year risk of CAD was highest in slow walkers at high genetic risk (top 20% of PGS): 2.72% (95% CI: 2.30-3.13) in women; 9.60% (8.62-10.57) in men. The risk difference between slow and brisk walkers was greater at higher [1.26% (0.81-1.71) in women; 3.63% (2.58-4.67) in men] than lower [0.76% (0.59-0.93) and 2.37% (1.96-2.78), respectively] genetic risk. Brisk walkers at high genetic risk had equivalent (women) or higher (men) risk than slow walkers at moderate-to-low genetic risk (bottom 80% of PGS). When added to a model containing traditional risk factors, both factors separately improved risk discrimination; combining them resulted in the greatest discrimination: C-index of 0.801 (0.793-0.808) in women; 0.732 (0.728-0.737) in men.

CONCLUSION

Self-reported slow walkers at high genetic risk had the greatest risk of CAD, identifying a potentially important population for intervention. Both PGS and walking pace contributed to risk discrimination.

Mendelian randomization with Egger pleiotropy correction and weakly informative Bayesian priors

Background

The MR-Egger (MRE) estimator has been proposed to correct for directional pleiotropic effects of genetic instruments in an instrumental variable (IV) analysis. The power of this method is considerably lower than that of conventional estimators, limiting its applicability. Here we propose a novel Bayesian implementation of the MR-Egger estimator (BMRE) and explore the utility of applying weakly informative priors on the intercept term (the pleiotropy estimate) to increase power of the IV (slope) estimate.

Methods

This was a simulation study to compare the performance of different IV estimators. Scenarios differed in the presence of a causal effect, the presence of pleiotropy, the proportion of pleiotropic instruments and degree of 'Instrument Strength Independent of Direct Effect' (InSIDE) assumption violation. Based on empirical plasma urate data, we present an approach to elucidate a prior distribution for the amount of pleiotropy.

Results

A weakly informative prior on the intercept term increased power of the slope estimate while maintaining type 1 error rates close to the nominal value of 0.05. Under the InSIDE assumption, performance was unaffected by the presence or absence of pleiotropy. Violation of the InSIDE assumption biased all estimators, affecting the BMRE more than the MRE method.

Conclusions

Depending on the prior distribution, the BMRE estimator has more power at the cost of an increased susceptibility to InSIDE assumption violations. As such the BMRE method is a compromise between the MRE and conventional IV estimators, and may be an especially useful approach to account for observed pleiotropy.

Comparison of variance estimators for meta-analysis of instrumental variable estimates

Background: Mendelian randomization studies perform instrumental variable (IV) analysis using genetic IVs. Results of individual Mendelian randomization studies can be pooled through meta-analysis. We explored how different variance estimators influence the meta-analysed IV estimate.

Methods: Two versions of the delta method (IV before or after pooling), four bootstrap estimators, a jack-knife estimator and a heteroscedasticity-consistent (HC) variance estimator were compared using simulation. Two types of meta-analyses were compared, a two-stage meta-analysis pooling results, and a one-stage meta-analysis pooling datasets.

Results: Using a two-stage meta-analysis, coverage of the point estimate using bootstrapped estimators deviated from nominal levels at weak instrument settings and/or outcome probabilities ≤ 0.10. The jack-knife estimator was the least biased resampling method, the HC estimator often failed at outcome probabilities ≤ 0.50 and overall the delta method estimators were the least biased. In the presence of between-study heterogeneity, the delta method before meta-analysis performed best. Using a one-stage meta-analysis all methods performed equally well and better than two-stage meta-analysis of greater or equal size.

Conclusions: In the presence of between-study heterogeneity, two-stage meta-analyses should preferentially use the delta method before meta-analysis. Weak instrument bias can be reduced by performing a one-stage meta-analysis.

Association analysis of 31 common polymorphisms with type 2 diabetes and its related traits in Indian sib pairs

AIMS/HYPOTHESIS

Evaluation of the association of 31 common single nucleotide polymorphisms (SNPs) with fasting glucose, fasting insulin, HOMA-beta cell function (HOMA-β), HOMA-insulin resistance (HOMA-IR) and type 2 diabetes in the Indian population.

METHODS

We genotyped 3,089 sib pairs recruited in the Indian Migration Study from four cities in India (Lucknow, Nagpur, Hyderabad and Bangalore) for 31 SNPs in 24 genes previously associated with type 2 diabetes in European populations. We conducted within-sib-pair analysis for type 2 diabetes and its related quantitative traits.

RESULTS

The risk-allele frequencies of all the SNPs were comparable with those reported in western populations. We demonstrated significant associations of CXCR4 (rs932206), CDKAL1 (rs7756992) and TCF7L2 (rs7903146, rs12255372) with fasting glucose, with β values of 0.007 (p = 0.05), 0.01 (p = 0.01), 0.007 (p = 0.05), 0.01 (p = 0.003) and 0.08 (p = 0.01), respectively. Variants in NOTCH2 (rs10923931), TCF-2 (also known as HNF1B) (rs757210), ADAM30 (rs2641348) and CDKN2A/B (rs10811661) significantly predicted fasting insulin, with β values of -0.06 (p = 0.04), 0.05 (p = 0.05), -0.08 (p = 0.01) and -0.08 (p = 0.02), respectively. For HOMA-IR, we detected associations with TCF-2, ADAM30 and CDKN2A/B, with β values of 0.05 (p = 0.04), -0.07 (p = 0.03) and -0.08 (p = 0.02), respectively. We also found significant associations of ADAM30 (β = -0.05; p = 0.01) and CDKN2A/B (β = -0.05; p = 0.03) with HOMA-β. THADA variant (rs7578597) was associated with type 2 diabetes (OR 1.5; 95% CI 1.04, 2.22; p = 0.03).

CONCLUSIONS/INTERPRETATION

We validated the association of seven established loci with intermediate traits related to type 2 diabetes in an Indian population using a design resistant to population stratification.

Associations of FTO and MC4R Variants with Obesity Traits in Indians and the Role of Rural/Urban Environment as a Possible Effect Modifier

Few studies have investigated the association between genetic variation and obesity traits in Indian populations or the role of environmental factors as modifiers of these relationships. In the context of rapid urbanisation, resulting in significant lifestyle changes, understanding the aetiology of obesity is important. We investigated associations of FTO and MC4R variants with obesity traits in 3390 sibling pairs from four Indian cities, most of whom were discordant for current dwelling (rural or urban). The FTO variant rs9939609 predicted increased weight (0.09 Z-scores, 95% CI: 0.03, 0.15) and BMI (0.08 Z-scores, 95% CI: 0.02, 0.14). The MC4R variant rs17782313 was weakly associated with weight and hip circumference (P < .05). There was some indication that the association between FTO and weight was stronger in urban than that in rural dwellers (P for interaction = .03), but no evidence for effect modification by diet or physical activity. Further studies are needed to investigate ways in which urban environment may modify genetic risk of obesity.

Genes related to sex steroids, neural growth, and social-emotional behavior are associated with autistic traits, empathy, and Asperger syndrome

Genetic studies of autism spectrum conditions (ASC) have mostly focused on the "low functioning" severe clinical subgroup, treating it as a rare disorder. However, ASC is now thought to be relatively common ( approximately 1%), and representing one end of a quasi-normal distribution of autistic traits in the general population. Here we report a study of common genetic variation in candidate genes associated with autistic traits and Asperger syndrome (AS). We tested single nucleotide polymorphisms in 68 candidate genes in three functional groups (sex steroid synthesis/transport, neural connectivity, and social-emotional responsivity) in two experiments. These were (a) an association study of relevant behavioral traits (the Empathy Quotient (EQ), the Autism Spectrum Quotient (AQ)) in a population sample (n=349); and (b) a case-control association study on a sample of people with AS, a "high-functioning" subgroup of ASC (n=174). 27 genes showed a nominally significant association with autistic traits and/or ASC diagnosis. Of these, 19 genes showed nominally significant association with AQ/EQ. In the sex steroid group, this included ESR2 and CYP11B1. In the neural connectivity group, this included HOXA1, NTRK1, and NLGN4X. In the socio-responsivity behavior group, this included MAOB, AVPR1B, and WFS1. Fourteen genes showed nominally significant association with AS. In the sex steroid group, this included CYP17A1 and CYP19A1. In the socio-emotional behavior group, this included OXT. Six genes were nominally associated in both experiments, providing a partial replication. Eleven genes survived family wise error rate (FWER) correction using permutations across both experiments, which is greater than would be expected by chance. CYP11B1 and NTRK1 emerged as significantly associated genes in both experiments, after FWER correction (P<0.05). This is the first candidate-gene association study of AS and of autistic traits. The most promising candidate genes require independent replication and fine mapping.

Bone structural effects of variation in the TNFRSF1B gene encoding the tumor necrosis factor receptor 2

The 1p36 region of the human genome has been identified as containing a QTL for BMD in multiple studies. We analysed the TNFRSF1B gene from this region, which encodes the TNF receptor 2, in two large population-based cohorts. Our results suggest that variation in TNFRSF1B is associated with BMD.

INTRODUCTION

The TNFRSF1B gene, encoding the TNF receptor 2, is a strong positional and functional candidate gene for impaired bone structure through the role that TNF has in bone cells. The aims of this study were to evaluate the role of variations in the TNFRSF1B gene on bone structure and osteoporotic fracture risk in postmenopausal women.

METHODS

Six SNPs in TNFRSF1B were analysed in a cohort of 1,190 postmenopausal Australian women, three of which were also genotyped in an independent cohort of 811 UK postmenopausal women. Differences in phenotypic means for genotype groups were examined using one-way ANOVA and ANCOVA.

RESULTS

Significant associations were seen for IVS1+5580A>G with BMD and QUS parameters in the Australian population (P = 0.008 - 0.034) and with hip BMD parameters in the UK population (P = 0.005 - 0.029). Significant associations were also observed between IVS1+6528G>A and hip BMD parameters in the UK cohort (P = 0.0002 - 0.003). We then combined the data from the two cohorts and observed significant associations between both IVS1+5580A>G and IVS1+6528G>A and hip BMD parameters (P = 0.002 - 0.033).

CONCLUSIONS

Genetic variation in TNFRSF1B plays a role in the determination of bone structure in Caucasian postmenopausal women, possibly through effects on osteoblast and osteoclast differentiation.

Mapping the platelet profile for functional genomic studies and demonstration of the effect size of the GP6 locus

BACKGROUND

Evidence suggests the wide variation in platelet response within the population is genetically controlled. Unraveling the complex relationship between sequence variation and platelet phenotype requires accurate and reproducible measurement of platelet response.

OBJECTIVE

To develop a methodology suitable for measuring signaling pathway-specific platelet phenotype, to use this to measure platelet response in a large cohort, and to demonstrate the effect size of sequence variation in a relevant model gene.

METHODS

Three established platelet assays were evaluated: mobilization of [Ca(2+)](i), aggregometry and flow cytometry, each in response to adenosine 5'-diphosphate (ADP) or the glycoprotein (GP) VI-specific crosslinked collagen-related peptide (CRP). Flow cytometric measurement of fibrinogen binding and P-selectin expression in response to a single, intermediate dose of each agonist gave the best combination of reproducibility and inter-individual variability and was used to measure the platelet response in 506 healthy volunteers. Pathway specificity was ensured by blocking the main subsidiary signaling pathways.

RESULTS

Individuals were identified who were hypo- or hyper-responders for both pathways, or who had differential responses to the two agonists, or between outcomes. 89 individuals, retested three months later using the same methodology, showed high concordance between the two visits in all four assays (r(2) = 0.872, 0.868, 0.766 and 0.549); all subjects retaining their phenotype at recall. The effect of sequence variation at the GP6 locus accounted for approximately 35% of the variation in the CRP-XL response.

CONCLUSION

Genotyping-phenotype association studies in a well-characterized, large cohort provides a powerful strategy to measure the effect of sequence variation in genes regulating the platelet response.

Definition of novel GP6 polymorphisms and major difference in haplotype frequencies between populations by a combination of in-depth exon resequencing and genotyping with tag single nucleotide polymorphisms

BACKGROUND

Common genetic variants of cell surface receptors contribute to differences in functional responses and disease susceptibility. We have previously shown that single nucleotide polymorphisms (SNPs) in platelet glycoprotein VI (GP6) determine the extent of response to agonist. In addition, SNPs in the GP6 gene have been proposed as risk factors for coronary artery disease.

METHODS

To completely characterize genetic variation in the GP6 gene we generated a high-resolution SNP map by sequencing the promoter, exons and consensus splice sequences in 94 non-related Caucasoids. In addition, we sequenced DNA encoding the ligand-binding domains of GP6 from non-human primates to determine the level of evolutionary conservation.

RESULTS

Eighteen SNPs were identified, six of which encoded amino acid substitutions in the mature form of the protein. The single non-synonymous SNP identified in the exons encoding the ligand-binding domains, encoding for a 103Leu > Val substitution, resulted in reduced ligand binding. Two common protein isoforms were confirmed in Caucasoid with frequencies of 0.82 and 0.15. Variation at the GP6 locus was characterized further by determining SNP frequency in over 2000 individuals from different ethnic backgrounds.

CONCLUSIONS

The SNPs were polymorphic in all populations studied although significant differences in allele frequencies were observed. Twelve additional GP6 protein isoforms were identified from the genotyping results and, despite extensive variation in GP6, the sequence of the ligand-binding domains is conserved. Sequences from non-human primates confirmed this observation. These data provide valuable information for the optimal selection of genetic variants for use in future association studies.

Association of the paternally transmitted copy of common Valine allele of the Val66Met polymorphism of the brain-derived neurotrophic factor (BDNF) gene with susceptibility to ADHD

Attention deficit hyperactivity disorder (ADHD) is a common, highly heritable, neurodevelopmental disorder with onset in early childhood. Genes involved in neuronal development and growth are, thus, important etiological candidates and brain-derived neurotrophic factor (BDNF), has been hypothesized to play a role in the pathogenesis of ADHD. BDNF is a member of the neurotrophin family and is involved in the survival and differentiation of dopaminergic neurons in the developing brain (of relevance because drugs that block the dopamine transporter can be effective therapeutically). The common Val66Met functional polymorphism in the human BDNF gene (rs 6265) was genotyped in a collaborative family-based sample of 341 white UK or Irish ADHD probands and their parents. We found evidence for preferential transmission of the valine (G) allele of BDNF (odds ratio, OR=1.6, P=0.02) with a strong paternal effect (paternal transmissions: OR=3.2, P=0.0005; maternal transmissions: OR=1.00; P=1.00). Our findings support the hypothesis that BDNF is involved in the pathogenesis of ADHD. The transmission difference between parents raises the possibility that an epigenetic process may be involved.

Genotype effects and epistasis in type 1 diabetes and HLA-DQ trans dimer associations with disease

Alleles of HLA class II genes DQB1, DQA1, and DRB1 in the MHC region are major determinants of genetic predisposition to type 1 diabetes (T1D). Several alleles of each of these three loci are associated with susceptibility or protection from disease. In addition, relative risks for some DR-DQ genotypes are not simply the sum or product of the single haplotype relative risks. For example, the risk of the DRB1*03-DQB1*02/DRB1*0401-DQB1*0302 genotype is often found to be higher than for the individual DRB1*03-DQB1*02 and DRB1*0401-DQB1*0302 homozygous genotypes. It has been hypothesized that this synergy or epistasis occurs through formation of highly susceptible trans-encoded HLA-DQ(alpha 1, beta 1) heterodimers. Here, we evaluated this hypothesis by estimating the disease associations of the range of DR-DQ genotypes and their inferred dimers in a large collection of nuclear families. We determined whether the risk of haplotypes in DRB1*0401-DQB1*0302-positive genotypes relative to the DRB1*03-DQB1*02-positive genotypes is different from that of DRB1*01-DQB1*0501, which we used as a baseline reference. Several haplotypes showed a different risk compared to DRB1*01-DQB1*0501, which correlated with their ability to form certain trans-encoded DQ dimers. This result provides new evidence for the potential importance of trans-encoded HLA DQ molecules in the determination of HLA-associated risk in T1D.

Mitochondrial dysfunction in schizophrenia: evidence for compromised brain metabolism and oxidative stress

The etiology and pathophysiology of schizophrenia remain unknown. A parallel transcriptomics, proteomics and metabolomics approach was employed on human brain tissue to explore the molecular disease signatures. Almost half the altered proteins identified by proteomics were associated with mitochondrial function and oxidative stress responses. This was mirrored by transcriptional and metabolite perturbations. Cluster analysis of transcriptional alterations showed that genes related to energy metabolism and oxidative stress differentiated almost 90% of schizophrenia patients from controls, while confounding drug effects could be ruled out. We propose that oxidative stress and the ensuing cellular adaptations are linked to the schizophrenia disease process and hope that this new disease concept may advance the approach to treatment, diagnosis and disease prevention of schizophrenia and related syndromes.

Functional genetic polymorphisms in cytokines and metabolic genes as additional genetic markers for susceptibility to develop type 1 diabetes

Genetic association with type 1 diabetes (T1D) has been established for two chromosomal regions: HLA DQ/DR (IDDM1) and INS VNTR (IDDM2). To identify additional genetic markers, we tested polymorphisms in regulatory regions of several cytokine and important metabolic genes. These polymorphisms exhibit functional consequences for expression and function. Functional genetic polymorphisms of proinflammatory (T-helper-1: IL-2, IL-12 and IFN-gamma), anti-inflammatory (T-helper-2: IL-4, IL-6 and IL-10) and metabolic (IGF-I, VDR and INS) genes were determined in 206 Dutch simplex families with juvenile onset T1D and the results were analysed using the transmission disequilibrium test. Significantly increased transmission to T1D probands was observed for the loci IDDM1, IDDM2 and the vitamin D receptor. Although none of the other individual polymorphisms was associated with disease individually, the combination of T-helper-2 and metabolic/growth alleles IL-10(*)R2, IL-4(*)C, VDR(*)C and IGF-I(*)wt was found to be transmitted more frequently than expected (67%, P(c)=0.015). We conclude that additional genetic predisposition to T1D is defined by combinations of markers (eg Th2 and metabolic) rather than by a single marker. The consequences of the increased transmission of a low Th2 expressing genotypes together with a normal Th1 profile may result in a net proinflammatory cytokine expression pattern.

Linkage analysis of Genetic Analysis Workshop 12 simulated data based on affected individuals only

A variety of allele-sharing methods was used to assess linkage to disease for a whole-genome screen in each of ten replicates from the Genetic Analysis Workshop 12 simulated data. Analysis of each replicate produced many false positive results in addition to a few "true" positives. Conditioning on the strongest initial linkage and reanalyzing using a variety of conditional methods did not improve the power for detection or help discriminate between true and false positive signals.

Haplotype tagging for the identification of common disease genes

Genome-wide linkage disequilibrium (LD) mapping of common disease genes could be more powerful than linkage analysis if the appropriate density of polymorphic markers were known and if the genotyping effort and cost of producing such an LD map could be reduced. Although different metrics that measure the extent of LD have been evaluated, even the most recent studies have not placed significant emphasis on the most informative and cost-effective method of LD mapping-that based on haplotypes. We have scanned 135 kb of DNA from nine genes, genotyped 122 single-nucleotide polymorphisms (SNPs; approximately 184,000 genotypes) and determined the common haplotypes in a minimum of 384 European individuals for each gene. Here we show how knowledge of the common haplotypes and the SNPs that tag them can be used to (i) explain the often complex patterns of LD between adjacent markers, (ii) reduce genotyping significantly (in this case from 122 to 34 SNPs), (iii) scan the common variation of a gene sensitively and comprehensively and (iv) provide key fine-mapping data within regions of strong LD. Our results also indicate that, at least for the genes studied here, the current version of dbSNP would have been of limited utility for LD mapping because many common haplotypes could not be defined. A directed re-sequencing effort of the approximately 10% of the genome in or near genes in the major ethnic groups would aid the systematic evaluation of the common variant model of common disease.

The HLA-DPB1--associated component of the IDDM1 and its relationship to the major loci HLA-DQB1, -DQA1, and -DRB1

The major histocompatibility complex (MHC) HLA region on chromosome 6p21 contains the major locus of type 1 diabetes (IDDM1). Common allelic variants at the class II HLA-DRB1, -DQA1, and -DQB1 loci account for the major part of IDDM1. Previous studies suggested that other MHC loci are likely to contribute to IDDM1, but determination of their relative contributions and identities is difficult because of strong linkage disequilibrium between MHC loci. One prime candidate is the polymorphic HLA-DPB1 locus, which (with the DPA1 locus) encodes the third class II antigen-presenting molecule. However, the results obtained in previous studies appear to be contradictory. Therefore, we have analyzed 408 white European families (200 from Sardinia and 208 from the U.K.) using a combination of association tests designed to directly compare the effect of DPB1 variation on the relative predisposition of DR-DQ haplotypes, taking into account linkage disequilibrium between DPB1 and the DRB1, DQA1, and DQB1 loci. In these populations, the overall contribution of DPB1 to IDDM1 is small. The main component of the DPB1 contribution to IDDM1 in these populations appears to be the protection associated with DPB1*0402 on DR4-negative haplotypes. We suggest that the HLA-DP molecule itself contributes to IDDM1.

Typhoid fever and genetic polymorphisms at the natural resistance-associated macrophage protein 1

Control of Salmonella enterica serovar Typhimurium (S. typhimurium) infection in the mouse model of typhoid fever is critically dependent on the natural resistance-associated macrophage protein 1 (Nramp1). In this study, we examined the role of genetic polymorphisms in the human homologue, NRAMP1, in resistance to typhoid fever in southern Vietnam. Patients with blood-culture-confirmed typhoid fever and healthy control subjects were genotyped for 6 polymorphic markers within and near NRAMP1 on chromosome 2q35. Four single base-pair polymorphisms (274 C/T, 469+14 G/C, 1465-85 G/A, and D543N), a (GT)(n) repeat in the promoter region of NRAMP1 and D2S1471, and a microsatellite marker approximately 130-kb downstream of NRAMP1 were examined. The allelic and genotypic frequencies for each polymorphism were compared in case patients and control subjects. No allelic association was identified between the NRAMP1 alleles and typhoid fever susceptibility. In addition, neither homozygotes nor heterozygotes for any NRAMP1 variants were at increased risk of typhoid fever.

Genes of the class II and class III major histocompatibility complex are associated with typhoid fever in Vietnam

The influence of genes of the major histocompatibility complex (MHC) class II and class III loci on typhoid fever susceptibility was investigated. Individuals with blood culture-confirmed typhoid fever and control subjects from 2 distinct geographic locations in southern Vietnam were genotyped for HLA-DRB1 and HLA-DQB1 alleles, the gene that encodes tumor necrosis factor (TNF)-alpha (TNFA [-238] and TNFA [-308]), the gene that encodes lymphotoxin-alpha, and alleles of the TNF-alpha microsatellite. HLA-DRB1*0301/6/8, HLA-DQB1*0201-3, and TNFA*2 (-308) were associated with susceptibility to typhoid fever, whereas HLA-DRB1*04, HLA-DQB1*0401/2, and TNFA*1 (-308) were associated with disease resistance. The frequency of all possible haplotypes of the 3 individually associated loci were estimated and were found to be significantly different in typhoid case patients and control subjects (chi2=55.56, 32 df; P=.006). Haplotypes that were either protective (TNFA*1 [-308].DRB1*04) or predisposed individuals to typhoid fever (TNFA*2 [-308].DRB1*0301) were determined. This report identifies a genetic association in humans between typhoid fever and MHC class II and III genes.

Unbiased application of the transmission/disequilibrium test to multilocus haplotypes

When the transmission/disequilibrium test (TDT) is applied to multilocus haplotypes, a bias may be introduced in some families for which both parents have the same heterozygous genotype at some locus. The bias occurs because haplotypes can only be deduced from certain offspring, with the result that the transmissions of the two parental haplotypes are not independent. We obtain an unbiased TDT for individual haplotypes by calculating the correct variance for the transmission count within a family, using information from multiple siblings if they are available. An existing correction for dependence between siblings in the presence of linkage is retained. To obtain an unbiased multihaplotype TDT, we must either count transmissions from one randomly chosen parent or count all transmissions and estimate the significance level empirically. Alternatively, we may use missing-data techniques to estimate uncertain haplotypes, but these methods are not robust to population stratification. An illustration using data from the insulin-gene region in type 1 diabetes shows that the validity and power of the TDT may vary by an order of magnitude, depending on the method of analysis.

Evaluation of fine mapping strategies for a multifactorial disease locus: systematic linkage and association analysis of IDDM1 in the HLA region on chromosome 6p21

The positional cloning of multifactorial disease genes is a major challenge in human genetics. We have therefore empirically tested the utility of the available polymorphic microsatellite map to locate the already identified type 1 diabetes locus IDDM1 (sibling risk/population prevalence ratio lambda(s)= 2.7) within a 14 Mb region of chromosome 6p21 linked to disease. In a two-stage approach to fine mapping, linkage was evaluated in 385 affected sib-pair families using 13 evenly spaced polymorphic microsatellite markers. The whole 14 Mb showed strong linkage. Then, each marker was analysed for evidence of allelic association, revealing evidence of disease association at one marker located within the 95% confidence interval of 1.7 cM obtained by linkage. Analysis of an additional 12 markers flanking this marker revealed a highly specific region of 570 kb associated with disease ( P = 7.5 x 10(-35)), which included the HLA class II genes, known to be the primary determinants of IDDM1. The peak of association was as close as 85 kb centromeric of the disease-predisposing class II gene HLA-DQB1. We investigated the importance of the underlying inter-marker linkage disequilibrium, marker informativity and recombination for fine mapping and demonstrate that the majority of disease association in the region can be explained by linkage disequilibrium with the class II susceptibility genes. Recombination within the major histocompatibility complex was rare and nearly absent in the class III region. We demonstrate that fine mapping of a multifactorial disease gene is possible with high accuracy even in a region with extraordinary linkage disequilibrium across distances of several Mb. The results will be applicable to association studies of disease loci with lambda(s)values <2.7 except that much larger data sets will be required.

Conditional ETDT analysis of the human leukocyte antigen region in type 1 diabetes

Several studies have indicated that additional genes in the major histocompatibility complex (MHC) region, other than the class II genes HLA-DQB1 and -DRB1 (the IDDM1 locus), may contribute to susceptibility and resistance to type 1 diabetes. The relative magnitude of these non- DR/DQ effects is uncertain and their map location is unknown owing to the extraordinary linkage disequilibrium that extends over the 3.5 Mb of the MHC. The homozygous parent test has been proposed as a method for detection of additional risk factors conditional on HLA-DQB1 and -DRB1. However, this method is inefficient since it uses only parents homozygous for the primary disease locus, the DQB1-DRB1 haplotype. To overcome this limitation, Conditional ETDT was used in the present report to test for association conditional on the DQB1-DRB1 haplotype, thereby allowing all parents to be included in the analysis. First, we confirm in UK and Sardinian type 1 diabetic families that allelic variation at HLA-DRB1 has a very significant effect on the association of DQB1 and vice versa. The Conditional ETDT was then applied to the HLA TNF (tumour necrosis factor) region and microsatellite marker D6S273 region, both of which have been reported to contribute to IDDM1 independent of the HLA-DQB1-DRB1 genes. We found no evidence for a major role for either of these two regions in IDDM1.

Adaptation of the extended transmission/disequilibrium test to distinguish disease associations of multiple loci: the Conditional Extended Transmission/Disequilibrium Test

Linkage and association studies in complex diseases are used to identify and fine map disease loci. The process of identifying the aetiological polymorphism, the molecular variant responsible for the linkage and association of the chromosome region with disease, is complicated by the low penetrance of the disease variant, the linkage disequilibrium between physically-linked polymorphic markers flanking the disease variant, and the possibility that more than one polymorphism in the most associated region is aetiological. It is important to be able to detect additional disease determinants in a region containing a cluster of genes, such as the major histocompatibility complex (MHC) region on chromosome 6p21. Some methods have been developed for detection of additional variants, such as the Haplotype Method, Marker Association Segregation Chi-squares (MASC) Method, and the Homozygous Parent Test. Here, the Extended Transmission/Disequilibrium Test is adapted to test for association conditional on a previously associated locus. This test is referred to as the Conditional Extended TDT (CETDT). We discuss the advantages of the CETDT compared to existing methods and, using simulated data, investigate the effect of polymorphism, inheritance, and linkage disequilibrium on the CETDT.