The hazards of genotype imputation when mapping disease susceptibility variants

BACKGROUND

The cost-free increase in statistical power of using imputation to infer missing genotypes is undoubtedly appealing, but is it hazard-free? This case study of three type-2 diabetes (T2D) loci demonstrates that it is not; it sheds light on why this is so and raises concerns as to the shortcomings of imputation at disease loci, where haplotypes differ between cases and reference panel.

RESULTS

T2D-associated variants were previously identified using targeted sequencing. We removed these significantly associated SNPs and used neighbouring SNPs to infer them by imputation. We compared imputed with observed genotypes, examined the altered pattern of T2D-SNP association, and investigated the cause of imputation errors by studying haplotype structure. Most T2D variants were incorrectly imputed with a low density of scaffold SNPs, but the majority failed to impute even at high density, despite obtaining high certainty scores. Missing and discordant imputation errors, which were observed disproportionately for the risk alleles, produced monomorphic genotype calls or false-negative associations. We show that haplotypes carrying risk alleles are considerably more common in the T2D cases than the reference panel, for all loci.

CONCLUSIONS

Imputation is not a panacea for fine mapping, nor for meta-analysing multiple GWAS based on different arrays and different populations. A total of 80% of the SNPs we have tested are not included in array platforms, explaining why these and other such associated variants may previously have been missed. Regardless of the choice of software and reference haplotypes, imputation drives genotype inference towards the reference panel, introducing errors at disease loci.

The hazards of genotype imputation in chromosomal regions under selection: A case study using the Lactase gene region

Although imputation of missing SNP results has been widely used in genetic studies, claims about the quality and usefulness of imputation have outnumbered the few studies that have questioned its limitations. But it is becoming clear that these limitations are real-for example, disease association signals can be missed in regions of LD breakdown. Here, as a case study, using the chromosomal region of the well-known lactase gene, LCT, we address the issue of imputation in the context of variants that have become frequent in a limited number of modern population groups only recently, due to selection. We study SNPs in a 500 bp region covering the enhancer of LCT, and compare imputed genotypes with directly genotyped data. We examine the haplotype pairs of all individuals with discrepant and missing genotypes. We highlight the nonrandom nature of the allelic errors and show that most incorrect imputations and missing data result from long haplotypes that are evolutionarily closely related to those carrying the derived alleles, while some relate to rare and recombinant haplotypes. We conclude that bias of incorrectly imputed and missing genotypes can decrease the accuracy of imputed results substantially.

New Insights Into Mitochondrial Dysfunction at Disease Susceptibility Loci in the Development of Type 2 Diabetes

This study investigated the potential genetic mechanisms which underlie adipose tissue mitochondrial dysfunction in Type 2 diabetes (T2D), by systematically identifying nuclear-encoded mitochondrial genes (NEMGs) among the genes regulated by T2D-associated genetic loci. The target genes of these 'disease loci' were identified by mapping genetic loci associated with both disease and gene expression levels (expression quantitative trait loci, eQTL) using high resolution genetic maps, with independent estimates co-locating to within a small genetic distance. These co-locating signals were defined as T2D-eQTL and the target genes as T2D cis-genes. In total, 763 cis-genes were associated with T2D-eQTL, of which 50 were NEMGs. Independent gene expression datasets for T2D and insulin resistant cases and controls confirmed that the cis-genes and cis-NEMGs were enriched for differential expression in cases, providing independent validation that genetic maps can identify informative functional genes. Two additional results were consistent with a potential role of T2D-eQTL in regulating the 50 identified cis-NEMGs in the context of T2D risk: (1) the 50 cis-NEMGs showed greater differential expression compared to other NEMGs and (2) other NEMGs showed a trend towards significantly decreased expression if their expression levels correlated more highly with the subset of 50 cis-NEMGs. These 50 cis-NEMGs, which are differentially expressed and associated with mapped T2D disease loci, encode proteins acting within key mitochondrial pathways, including some of current therapeutic interest such as the metabolism of branched-chain amino acids, GABA and biotin.

A novel human pain insensitivity disorder caused by a point mutation in ZFHX2

Chronic pain is a major global public health issue causing a severe impact on both the quality of life for sufferers and the wider economy. Despite the significant clinical burden, little progress has been made in terms of therapeutic development. A unique approach to identifying new human-validated analgesic drug targets is to study rare families with inherited pain insensitivity. Here we have analysed an otherwise normal family where six affected individuals display a pain insensitive phenotype that is characterized by hyposensitivity to noxious heat and painless bone fractures. This autosomal dominant disorder is found in three generations and is not associated with a peripheral neuropathy. A novel point mutation in ZFHX2, encoding a putative transcription factor expressed in small diameter sensory neurons, was identified by whole exome sequencing that segregates with the pain insensitivity. The mutation is predicted to change an evolutionarily highly conserved arginine residue 1913 to a lysine within a homeodomain. Bacterial artificial chromosome (BAC) transgenic mice bearing the orthologous murine p.R1907K mutation, as well as Zfhx2 null mutant mice, have significant deficits in pain sensitivity. Gene expression analyses in dorsal root ganglia from mutant and wild-type mice show altered expression of genes implicated in peripheral pain mechanisms. The ZFHX2 variant and downstream regulated genes associated with a human pain-insensitive phenotype are therefore potential novel targets for the development of new analgesic drugs.awx326media15680039660001.

World-wide distributions of lactase persistence alleles and the complex effects of recombination and selection

The genetic trait of lactase persistence (LP) is associated with at least five independent functional single nucleotide variants in a regulatory region about 14 kb upstream of the lactase gene [−13910*T (rs4988235), −13907*G (rs41525747), −13915*G (rs41380347), −14009*G (rs869051967) and −14010*C (rs145946881)]. These alleles have been inferred to have spread recently and present-day frequencies have been attributed to positive selection for the ability of adult humans to digest lactose without risk of symptoms of lactose intolerance. One of the inferential approaches used to estimate the level of past selection has been to determine the extent of haplotype homozygosity (EHH) of the sequence surrounding the SNP of interest. We report here new data on the frequencies of the known LP alleles in the ‘Old World’ and their haplotype lineages. We examine and confirm EHH of each of the LP alleles in relation to their distinct lineages, but also show marked EHH for one of the older haplotypes that does not carry any of the five LP alleles. The region of EHH of this (B) haplotype exactly coincides with a region of suppressed recombination that is detectable in families as well as in population data, and the results show how such suppression may have exaggerated haplotype-based measures of past selection.

High-Resolution Genetic Maps Identify Multiple Type 2 Diabetes Loci at Regulatory Hotspots in African Americans and Europeans

Interpretation of results from genome-wide association studies for T2D is challenging. Only very few loci have been replicated in African ancestry populations and the identification of the implicated functional genes remain largely undefined. We used genetic maps that capture detailed linkage disequilibrium information in European and African Americans and applied these to large T2D case-control samples in order to estimate locations for putative functional variants in both populations. Replicated T2D locations were tested for evidence of being regulatory hotspots using adipose expression. We validated a sample of our co-location intervals using next generation sequencing and functional annotation, including enhancers, transcription, and chromatin modifications. We identified 111 additional disease-susceptibility locations, 93 of which are cosmopolitan and 18 of which are European specific. We show that many previously known signals are also risk loci in African Americans. The majority of the disease locations appear to confer risk of T2D via the regulation of expression levels for a large number (266) of cis-regulated genes, the majority of which are not the nearest genes to the disease loci. Sequencing three cosmopolitan locations provided candidate functional variants that precisely co-locate with cell-specific chromatin domains and pancreatic islet enhancers. These variants have large effect sizes and are common across populations. Results show that disease-associated loci in different populations, gene expression, and cell-specific regulatory annotation can be effectively integrated by localizing these effects on high-resolution genetic maps. The cis-regulated genes provide insights into the complex molecular pathways involved and can be used as targets for sequencing and functional molecular studies.

The genomic and phenotypic diversity of Schizosaccharomyces pombe

Natural variation within species reveals aspects of genome evolution and function. The fission yeast Schizosaccharomyces pombe is an important model for eukaryotic biology, but researchers typically use one standard laboratory strain. To extend the usefulness of this model, we surveyed the genomic and phenotypic variation in 161 natural isolates. We sequenced the genomes of all strains, finding moderate genetic diversity (π = 3 × 10(-3) substitutions/site) and weak global population structure. We estimate that dispersal of S. pombe began during human antiquity (∼340 BCE), and ancestors of these strains reached the Americas at ∼1623 CE. We quantified 74 traits, finding substantial heritable phenotypic diversity. We conducted 223 genome-wide association studies, with 89 traits showing at least one association. The most significant variant for each trait explained 22% of the phenotypic variance on average, with indels having larger effects than SNPs. This analysis represents a rich resource to examine genotype-phenotype relationships in a tractable model.

ABCC5 transporter is a novel type 2 diabetes susceptibility gene in European and African American populations

Numerous functional studies have implicated PARL in relation to type 2 diabetes (T2D). We hypothesised that conflicting human association studies may be due to neighbouring causal variants being in linkage disequilibrium (LD) with PARL. We conducted a comprehensive candidate gene study of the extended LD genomic region that includes PARL and transporter ABCC5 using three data sets (two European and one African American), in relation to healthy glycaemic variation, visceral fat accumulation and T2D disease. We observed no evidence for previously reported T2D association with Val262Leu or PARL using array and fine-map genomic and expression data. By contrast, we observed strong evidence of T2D association with ABCC5 (intron 26) for European and African American samples (P = 3E-07) and with ABCC5 adipose expression in Europeans [odds ratio (OR) = 3.8, P = 2E-04]. The genomic location estimate for the ABCC5 functional variant, associated with all phenotypes and expression data (P = 1E-11), was identical for all samples (at Chr3q 185,136 kb B36), indicating that the risk variant is an expression quantitative trait locus (eQTL) with increased expression conferring risk of disease. That the association with T2D is observed in populations of disparate ancestry suggests the variant is a ubiquitous risk factor for T2D.

Refinement in localization and identification of gene regions associated with Crohn disease

The risk of Crohn disease (CD) has a large genetic component. A recent meta-analysis of 6 genome-wide association studies reported 71 chromosomal intervals but does not account for all of the known genetic contribution. Here, we refine localization of the previously reported intervals and also identify additional CD susceptibility genes using a mapping approach that localizes causal variants based on genetic maps in linkage disequilibrium units (LDU maps). Using 2 of the 6 cohorts, 66 of the 71 previously reported loci are confirmed and more precise location estimates for these intervals are given. We identify 78 additional gene regions that pass genome-wide significance, providing strong evidence for 144 genes. Additionally, 56 nominally significant signals, but with more stringent and precise colocalization, are identified. In total, we provide evidence for 200 gene regions confirming that CD is truly multifactorial and complex in nature. Many identified genes have functions that are compatible with involvement in immune/inflammatory processes and seem to have a large effect in individuals with extra ileal as well as ileal inflammation. The precise locations and the evidence that some genes reflect phenotypic subgroups will help identify functional variants and will lead to greater insight of CD etiology.

A candidate-gene association study for berry colour and anthocyanin content in Vitis vinifera L

Anthocyanin content is a trait of major interest in Vitis vinifera L. These compounds affect grape and wine quality, and have beneficial effects on human health. A candidate-gene approach was used to identify genetic variants associated with anthocyanin content in grape berries. A total of 445 polymorphisms were identified in 5 genes encoding transcription factors and 10 genes involved in either the biosynthetic pathway or transport of anthocyanins. A total of 124 SNPs were selected to examine association with a wide range of phenotypes based on RP-HPLC analysis and visual characterization. The phenotypes were total skin anthocyanin (TSA) concentration but also specific types of anthocyanins and relative abundance. The visual assessment was based on OIV (Organisation Internationale de la Vigne et du Vin) descriptors for berry and skin colour. The genes encoding the transcription factors MYB11, MYBCC and MYC(B) were significantly associated with TSA concentration. UFGT and MRP were associated with several different types of anthocyanins. Skin and pulp colour were associated with nine genes (MYB11, MYBCC, MYC(B), UFGT, MRP, DFR, LDOX, CHI and GST). Pulp colour was associated with a similar group of 11 genes (MYB11, MYBCC, MYC(B), MYC(A), UFGT, MRP, GST, DFR, LDOX, CHI and CHS(A)). Statistical interactions were observed between SNPs within the transcription factors MYB11, MYBCC and MYC(B). SNPs within LDOX interacted with MYB11 and MYC(B), while SNPs within CHI interacted with MYB11 only. Together, these findings suggest the involvement of these genes in anthocyanin content and on the regulation of anthocyanin biosynthesis. This work forms a benchmark for replication and functional studies.

Dissecting the genetics of complex inheritance: linkage disequilibrium mapping provides insight into Crohn disease

Family studies for Crohn disease (CD) report extensive linkage on chromosome 16q and pinpoint NOD2 as a possible causative locus. However, linkage is also observed in families that do not bear the most frequent NOD2 causative mutations, but no other signals on 16q have been found so far in published genome-wide association studies. Our aim is to identify this missing genetic contribution. We apply a powerful genetic mapping approach to the Wellcome Trust Case-Control Consortium and the National Institute of Diabetes and Digestive and Kidney Diseases genome-wide association data on CD. This method takes into account the underlying structure of linkage disequilibrium (LD) by using genetic distances from LD maps and provides a location for the causal agent. We find genetic heterogeneity within the NOD2 locus and also show an independent and unsuspected involvement of the neighboring gene, CYLD. We find associations with the IRF8 region and the region containing CDH1 and CDH3, as well as substantial phenotypic and genetic heterogeneity for CD itself. The genes are known to be involved in inflammation and immune dysregulation. These findings provide insight into the genetics of CD and suggest promising directions for understanding disease heterogeneity. The application of this method thus paves the way for understanding complex inheritance in general, leading to the dissection of different pathways and ultimately, personalized treatment.

Comparative analysis of genome-wide association studies signals for lipids, diabetes, and coronary heart disease: Cardiovascular Biomarker Genetics Collaboration

Aims

To evaluate the associations of emergent genome-wide-association study-derived coronary heart disease (CHD)-associated single nucleotide polymorphisms (SNPs) with established and emerging risk factors, and the association of genome-wide-association study-derived lipid-associated SNPs with other risk factors and CHD events.

Methods and results

Using two case–control studies, three cross-sectional, and seven prospective studies with up to 25 000 individuals and 5794 CHD events we evaluated associations of 34 genome-wide-association study-identified SNPs with CHD risk and 16 CHD-associated risk factors or biomarkers. The Ch9p21 SNPs rs1333049 (OR 1.17; 95% confidence limits 1.11–1.24) and rs10757274 (OR 1.17; 1.09–1.26), MIA3 rs17465637 (OR 1.10; 1.04–1.15), Ch2q36 rs2943634 (OR 1.08; 1.03–1.14), APC rs383830 (OR 1.10; 1.02, 1.18), MTHFD1L rs6922269 (OR 1.10; 1.03, 1.16), CXCL12 rs501120 (OR 1.12; 1.04, 1.20), and SMAD3 rs17228212 (OR 1.11; 1.05, 1.17) were all associated with CHD risk, but not with the CHD biomarkers and risk factors measured. Among the 20 blood lipid-related SNPs, LPL rs17411031 was associated with a lower risk of CHD (OR 0.91; 0.84–0.97), an increase in Apolipoprotein AI and HDL-cholesterol, and reduced triglycerides. SORT1 rs599839 was associated with CHD risk (OR 1.20; 1.15–1.26) as well as total- and LDL-cholesterol, and apolipoprotein B. ANGPTL3 rs12042319 was associated with CHD risk (OR 1.11; 1.03, 1.19), total- and LDL-cholesterol, triglycerides, and interleukin-6.

Conclusion

Several SNPs predicting CHD events appear to involve pathways not currently indexed by the established or emerging risk factors; others involved changes in blood lipids including triglycerides or HDL-cholesterol as well as LDL-cholesterol. The overlapping association of SNPs with multiple risk factors and biomarkers supports the existence of shared points of regulation for these phenotypes.

Rohrer's constant, K2, as a factor of determining inspiratory resistance of common adult endotracheal tubes

The aim of the study was to calculate the in vitro inspiratory resistance (R(ETT)) of adult endotracheal tubes (ETT), via the end-inspiratory occlusion method, and to apply this method in vivo in order to estimate R(ETT) value in real time. By plotting R(ETT) over inspiratory flow (V) and calculating Rohrer's coefficients of linear and nonlinear resistance, K1 and K2 respectively, we determined the resistive behaviour of each ETT. Peak and plateau pressures were recorded at both proximal and distal sites of the ETT after applying a three-second occlusion under constant flow. Distal pressure was obtained via an intraluminal catheter R(ETT) was calculated as (P(peak) - P(plateau))/(V), at both sites. R(ETT) value resulted from the difference R(proximal) - R(distal). Graph R(ETT) over (V) was plotted and Rohrer's constants were calculated by the method of least squares. For ETTs with inner diameter 9.0, 8.5, 8.0, 7.5, 7.0 and 6.5 mm, K2 was 2.42, 3.05, 4.65, 6.01, 9.17 and 12.80 cmH2O/l/s, respectively. The intraluminal catheter increased R(ETT) No.7.0 by an average of 49%. Finally, ten patients with partially obstructed ETTs were tested and K2 in vivo constants found to be higher than their corresponding in vitro values (P value 0.00012). Therefore, knowing the performing size of an ETT may help the clinicians identify ETT obstruction and deal with weaning problems.

APOE, CETP and LPL genes show strong association with lipid levels in Greek children

BACKGROUND AND AIMS

Studies have consistently demonstrated that variants in a number of candidate genes are significant determinants of lipid levels in adults. However, few studies have investigated the impact of these variants in children. Therefore, in the present investigation we examined the influence of ten common variants in the genes for lipoprotein lipase (LPL-S447X), cholesterol ester transfer protein (CETP-Taq1B) apolipoprotein (APO) E (epsilon2, epsilon3, epsilon4), APOA5 (-1131C>T and S19W), APOA4 (S347T) and APOC3 (-482C>T; 1100C>T and 3238G>C) on lipoprotein levels children from the Gene-Diet Attica Investigation on childhood obesity (GENDAI).

METHODS AND RESULTS

The ten variants selected were genotyped in 882 Greek children, mean age: 11.2+/-0.7 years (418 females and 464 males). Genotypes were assessed using TaqMan technology. Significantly higher total cholesterol (TC) (p=0.0001) and low-density lipoprotein cholesterol (LDL-C) (p<0.0001) were observed in APOE epsilon4 carriers compared to epsilon3/epsilon3 homozygotes and epsilon2 carriers. The association of APOE genotype with TC and high-density lipoprotein cholesterol (HDL-C) ratio (p=0.0008) was further modulated by body mass index. Carriers of the CETP TaqIB B2 allele had significantly higher HDL-C (p<0.0001) and significantly lower TC: HDL-C ratio (p<0.0001) compared to B1/B1 individuals. No significant associations were observed between APOA4, APOA5 and APOC3 variants and serum lipids.

CONCLUSION

This study demonstrates that these common variants are associated with lipid levels in this healthy paediatric cohort, suggesting that even in these young children there may be potential in predicting their lifelong exposure to an adverse lipid profile.

Effects of single SNPs, haplotypes, and whole-genome LD maps on accuracy of association mapping

We describe an association mapping approach that utilizes linkage disequilibrium (LD) maps in LD units (LDU). This method uses composite likelihood to combine information from all single marker tests, and applies a model with a parameter for the location of the causal polymorphism. Previous analyses of the poor drug metabolizer phenotype provided evidence of the substantial utility of LDU maps for disease gene association mapping. Using LDU locations for the 27 single nucleotide polymorphisms (SNPs) flanking the CYP2D6 gene on chromosome 22, the most common functional polymorphism within the gene was located at 15 kb from its true location. Here, we examine the performance of this mapping approach by exploiting the high-density LDU map constructed from the HapMap data. Expressing the locations of the 27 SNPs in LDU from the HapMap LDU map, analysis yielded an estimated location that is only 0.3 kb away from the CYP2D6 gene. This supports the use of the high marker density HapMap-derived LDU map for association mapping even though it is derived from a much smaller number of individuals compared to the CYP2D6 sample. We also examine the performance of 2-SNP haplotypes. Using the same modelling procedures and composite likelihood as for single SNPs, the haplotype data provided much poorer localization compared to single SNP analysis. Haplotypes generate more autocorrelation through multiple inclusions of the same SNPs, which could inflate significance in association studies. The results of the present study demonstrate the great potential of the genome HapMap LDU maps for high-resolution mapping of complex phenotypes.

Genome scanning by composite likelihood

Ambitious programs have recently been advocated or launched to create genomewide databases for meta-analysis of association between DNA markers and phenotypes of medical and/or social concern. A necessary but not sufficient condition for success in association mapping is that the data give accurate estimates of both genomic location and its standard error, which are provided for multifactorial phenotypes by composite likelihood. That class includes the Malecot model, which we here apply with an illustrative example. This preliminary analysis leads to five inferences: permutation of cases and controls provides a test of association free of autocorrelation; two hypotheses give similar estimates, but one is consistently more accurate; estimation of the false-discovery rate is extended to causal genes in a small proportion of regions; the minimal data for successful meta-analysis are inferred; and power is robust for all genomic factors except minor-allele frequency. An extension to meta-analysis is proposed. Other approaches to genome scanning and meta-analysis should, if possible, be similarly extended so that their operating characteristics can be compared.

Linkage disequilibrium maps and disease-association mapping

Over the last few years, association mapping of disease genes has developed into one of the most dynamic research areas of human genetics. It focuses on identifying functional polymorphisms that predispose to complex diseases. Population-based approaches are concerned with exploiting linkage disequilibrium (LD) between single-nucleotide polymorphism (SNPs) and disease-predisposing loci. The utility of SNPs in association mapping is now well established and the interest in this field has been escalated by the discovery of millions of SNPs across the genome. This chapter reviews an association-mapping method that utilizes metric LD maps in LD units and employs a composite likelihood approach to combine information from all single SNP tests. It applies a model that incorporates a parameter for the location of the causal polymorphism. A proof-of-principle application of this method to a small region is given and its potential properties to large-scale datasets are discussed.

Refined Association Mapping for a Quantitative Trait: Weight in the H19-IGF2-INS-TH Region

Previous analyses have provided evidence for one or more loci affecting body weight in the H19-IGF2-INS-TH region on chromosome 11p15. To identify the location of a possible causal locus or loci we applied association analysis by composite likelihood to a large cohort under the Malecot model for body weight. A random sample of 2731 men in the UK were typed for eleven single nucleotide polymorphisms (SNPs) in IGF2, two SNPs in H19, one SNP in INS and one microsatellite marker in the TH genes. Using F tests appropriate to small marker sets, the superiority of regression over correlation was confirmed. All the evidence for association came from IGF2, with P= 0.007 for height-adjusted weight and P= 0.019 for weight additionally adjusted for smoking and alcohol drinking. Although the estimated point location for the suspected causal variant was close to IGF2 ApaI, the 95% confidence and support intervals covered most of IGF2 but none of the other loci. Identification of the causal SNP or SNPs within IGF2 will require typing of more variants in this region.

A map of the human genome in linkage disequilibrium units

Two genetic maps with additive distances contribute information about recombination patterns, recombinogenic sequences, and discovery of genes affecting a particular phenotype. Recombination is measured in morgans (w) over a single generation in a linkage map but may cover thousands of generations in a linkage disequilibrium (LD) map measured in LD units (LDU). We used a subset of single nucleotide polymorphisms from the HapMap Project to create a genome-wide map in LDU. Recombination accounts for 96.8% of the LDU variance in chromosome arms and 92.4% in their deciles. However, deeper analysis shows that LDU/w, an estimate of the effective bottleneck time (t), is significantly variable among chromosome arms because (i) the linkage map is approximated from the Haldane function, then adjusted toward the Kosambi function that is more accurate but still exaggerates w for all chromosomes, especially shorter ones; (ii) the non-pseudoautosomal region of the X chromosome is subject to hemizygous selection; and (iii) at resolution less than approximately 40,000 markers per w, there are indeterminacies (holes) in the LD map reflecting intervals of very high recombination. Selection and stochastic variation in small regions must have effects, which remain to be investigated by comparisons among populations. These considerations suggest an optimal strategy to eliminate holes quickly, greatly enhance the resolution of sex-specific linkage maps, and maximize the gain in association mapping by using LD maps.

Polymorphisms in a disintegrin and metalloprotease 33 (ADAM33) predict impaired early-life lung function

RATIONALE

Asthma commonly originates in early life in association with impaired lung function, which tracks to adulthood.

OBJECTIVES

Within the context of a prospective birth cohort study, we investigated the association between single nucleotide polymorphisms (SNPs) in a disintegrin and metalloprotease 33 (ADAM33) gene and early-life lung function.

METHODS

Children were genotyped for 17 SNPs in ADAM33. Lung function at age 3 (n = 285) and 5 years (n = 470) was assessed using plethysmographic measurement of specific airway resistance (sRaw). At age 5, we also measured FEV(1). SNPs were analyzed individually using logistic regression, followed by linkage disequilibrium mapping to identify the causal locus.

MAIN RESULTS

Carriers of the rare allele of F+1 SNP had reduced lung function at age 3 years (p = 0.003). When the recessive model was considered, four SNPs (F+1, S1, ST+5, V4) showed association with sRaw at age 5 years (p < 0.04). Using linkage disequilibrium mapping, we found evidence of a significant causal location between BC+1 and F1 SNPs, at the 5' end of the gene. Four SNPs were associated with lower FEV(1) (F+1, M+1, T1, and T2; p < or = 0.04). The risk of transient early wheezing more than doubled among children homozygous for the A allele of F+1 (odds ratio, 2.39; 95% confidence intervals, 1.18-4.86; p = 0.02), but there was no association between any SNP and allergic sensitization or physician-diagnosed asthma.

CONCLUSIONS

Polymorphisms in ADAM33 predict impaired early-life lung function. The functionally relevant polymorphism is likely to be at the 5' end of the gene.

The optimal measure of linkage disequilibrium reduces error in association mapping of affection status

We have developed a simple yet powerful approach for disease gene association mapping by linkage disequilibrium (LD). This method is unique because it applies a model with evolutionary theory that incorporates a parameter for the location of the causal polymorphism. The method exploits LD maps, which assign a location in LD units (LDU) for each marker. This approach is based on single marker tests within a composite likelihood framework, which avoids the heavy Bonferroni correction through multiple testing. As a proof of principle, we tested an 890 kb region flanking the CYP2D6 gene associated with poor drug-metabolizing activity in order to refine the localization of a causal mutation. Previous LD mapping studies using single markers and haplotypes have identified a 390 kb significant region associated with the poor drug-metabolizing phenotype on chromosome 22. None of the 27 Single nucleotide polymorphisms was within the gene. Using a metric LDU map, the commonest functional polymorphism within the gene was located at 14.9 kb from its true location, surrounded within a 95% confidence interval of 172 kb. The kb map had a relative efficiency of 33% compared with the LDU map. Our findings indicate that the support interval and location error are smaller than any published results. Despite the low resolution and the strong LD in the region, our results provide evidence of the substantial utility of LDU maps for disease gene association mapping. These tests are robust to large numbers of markers and are applicable to haplotypes, diplotypes, whole-genome association or candidate region studies.

Positional cloning by linkage disequilibrium

Recently, metric linkage disequilibrium (LD) maps that assign an LD unit (LDU) location for each marker have been developed (Maniatis et al. 2002). Here we present a multiple pairwise method for positional cloning by LD within a composite likelihood framework and investigate the operating characteristics of maps in physical units (kb) and LDU for two bodies of data (Daly et al. 2001; Jeffreys et al. 2001) on which current ideas of blocks are based. False-negative indications of a disease locus (type II error) were examined by selecting one single-nucleotide polymorphism (SNP) at a time as causal and taking its allelic count (0, 1, or 2, for the three genotypes) as a pseudophenotype, Y. By use of regression and correlation, association between every pseudophenotype and the allelic count of each SNP locus (X) was based on an adaptation of the Malecot model, which includes a parameter for location of the putative gene. By expressing locations in kb or LDU, greater power for localization was observed when the LDU map was fitted. The efficiency of the kb map, relative to the LDU map, to describe LD varied from a maximum of 0.87 to a minimum of 0.36, with a mean of 0.62. False-positive indications of a disease locus (type I error) were examined by simulating an unlinked causal SNP and the allele count was used as a pseudophenotype. The type I error was in good agreement with Wald's likelihood theorem for both metrics and all models that were tested. Unlike tests that select only the most significant marker, haplotype, or haploset, these methods are robust to large numbers of markers in a candidate region. Contrary to predictions from tagging SNPs that retain haplotype diversity, the sample with smaller size but greater SNP density gave less error. The locations of causal SNPs were estimated with the same precision in blocks and steps, suggesting that block definition may be less useful than anticipated for mapping a causal SNP. These results provide a guide to efficient positional cloning by SNPs and a benchmark against which the power of positional cloning by haplotype-based alternatives may be measured.

A Metric Linkage Disequilibrium Map of a Human Chromosome

We used LDMAP (Maniatis et al. 2002) to analyse SNP data spanning chromosome 22 (Dawson et al. 2002), to obtain a whole-chromosome metric LD map. The LD map, with map distances analogous to the centiMorgan scale of linkage maps, identifies regions of high LD as plateaus ('blocks') and characterises steps which define the relationship between these regions. From this map we estimate that block regions comprise between 32% and 55% of the euchromatic portion of chromosome 22 and that increasing marker density within steps may increase block coverage. Steps are regions of low LD which correspond to areas of variable recombination intensity. The intensity of recombination is related to the height of the step and thus intense recombination hot-spots can be distinguished from more randomly distributed historical events. The LD maps are more closely related to the high-resolution linkage map (Kong et al. 2002) than average measures of rho with recombination accounting for between 34% and 52% of the variance in patterns of LD (r=0.58 - 0.71, p=0.0001). Step regions are closely correlated with a range of sequence motifs including GT/CA repeats. The LD map identifies holes in which greater marker density is required and defines the optimal SNP spacing for positional cloning, which suggests that some multiple of around 50,000 SNPs will be required to efficiently screen Caucasian genomes. Further analyses which investigate selection of informative SNPs and the effect of SNP allele frequency and marker density will refine this estimate.

Linkage disequilibrium in human populations

Whereas the human linkage map appears on limited evidence to be constant over populations, maps of linkage disequilibrium (LD) vary among populations that differ in gene history. The greatest difference is between populations of sub-Saharan origin and populations remotely derived from Africa after a major bottleneck that reduced their heterozygosity and altered their Malecot parameters, increasing the intercept M that reflects association in founders and decreasing the exponential decline epsilon. Variation among populations within this ethnic dichotomy is much smaller. These observations validate use of a cosmopolitan LD map based on a sizeable sample representing a large population reliably typed for markers at high density. Then an LD map for a region or isolate within an ethnic group may be created by fitting the sample LD to the cosmopolitan map, estimating Malecot parameters simultaneously. The cosmopolitan map scaled by epsilon recovers 95% of the information that a local map at the same density gives and therefore more than the information in a low-resolution local map. Relative to a Eurasian cosmopolitan map the scaling factors are estimated to be 0.82 for isolates of European descent, 1.53 for Yorubans, and 1.74 for African Americans. These observations are consistent with a common bottleneck (perhaps but not necessarily speciation) approximately 173,500 years ago, if the bottleneck associated with migration out of Africa was 100,000 years ago. Eurasian populations (especially isolates with numerous cases) are efficient for genome scans, and populations of recent African origin (such as African Americans) are efficient for identification of causal polymorphisms within a candidate sequence.

The impact of data structure on genetic (co)variance components of early growth in sheep, estimated using an animal model with maternal effects

Several studies have noted high negative correlations between maternal genetic and direct additive effects and their influence on additive and maternal heritability of early growth traits in sheep. Multigeneration data from the Suffolk Sire Reference Scheme (SSRS) were used to investigate the effect of data structure on estimates of direct and maternal (co)variances for lamb 8-wk weight. In all analyses the additive, maternal genetic, maternal environmental, and residual effects were fitted along with the covariance between direct and maternal additive effects. The contributions of particular genetic relationships to the estimates were studied by analyzing subsets of the SSRS data. A further eight subsets were formed having 10% or 50% of the dams with their own records and having one or two, three or four, five or six, and more than six offspring per dam. Analysis of data having only 10% of the dams with their own record and one or two offspring records yielded a high negative correlation (-0.99) between direct and maternal genetic effects. However, the seven other data sets with more records per dam or a higher proportion of dams with their own records produced values of -0.35 to -0.51. Data structure and the number of dams and granddams with records are important determinants of estimated direct and maternal effects in early growth traits.

Properties of linkage disequilibrium (LD) maps

A linkage disequilibrium map is expressed in linkage disequilibrium (LD) units (LDU) discriminating blocks of conserved LD that have additive distances and locations monotonic with physical (kb) and genetic (cM) maps. There is remarkable agreement between LDU steps and sites of meiotic recombination in the one body of data informative for crossing over, and good agreement with another method that defines blocks without assigning an LD location to each marker. The map may be constructed from haplotypes or diplotypes, and efficiency estimated from the empirical variance of LD is substantially greater for the rho metric based on evolutionary theory than for the absolute correlation r, and for the LD map compared with its physical counterpart. The empirical variance is nearly three times as great for the worst alternative (r and kb map) as for the most efficient approach (rho and LD map). According to the empirical variances, blocks are best defined by zero distance between included markers. Because block size is algorithm-dependent and highly variable, the number of markers required for positional cloning is minimized by uniform spacing on the LD map, which is estimated to have approximately equal 1 LDU per locus, but with much variation among regions. No alternative representation of linkage disequilibrium (some of which are loosely called maps) has these properties, suggesting that LD maps are optimal for positional cloning of genes determining disease susceptibility.

The first linkage disequilibrium (LD) maps: delineation of hot and cold blocks by diplotype analysis

Linkage disequilibrium (LD) provides information about positional cloning, linkage, and evolution that cannot be inferred from other evidence, even when a correct sequence and a linkage map based on more than a handful of families become available. We present theory to construct an LD map for which distances are additive and population-specific maps are expected to be approximately proportional. For this purpose, there is only a modest difference in relative efficiency of haplotypes and diplotypes: resolving the latter into 2-locus haplotypes has significant cost or error and increases information by about 50%. LD maps for a cold spot in 19p13.3 and a more typical region in 3q21 are optimized by interval estimates. For a random sample and trustworthy map the value of LD at large distance can be predicted reliably from information over a small distance and does not depend on the evolutionary variance unless the sample size approaches the population size. Values of the association probability that can be distinguished from the value at large distance are determined not by population size but by time since a critical bottleneck. In these examples, omission of markers with significant Hardy-Weinberg disequilibrium does not improve the map, and widely discrepant draft sequences have similar estimates of the genetic parameters. The LD cold spot in 19p13.3 gives an unusually high estimate of time, supporting an argument that this relationship is general. As predicted for a region with ancient haplotypes or uniformly high recombination, there is no clear evidence of LD clustering. On the contrary, the 3q21 region is resolved into alternating blocks of stable and decreasing LD, as expected from crossover clustering. Construction of a genomewide LD map requires data not yet available, which may be complemented but not replaced by a catalog of haplotypes.

Nuclear, cytoplasmic, and environmental effects on growth, fat, and muscle traits in suffolk lambs from a sire referencing scheme

Maternal effects are an important source of variation in early growth and body traits in sheep but are often excluded from genetic analyses. Maternal additive genetic, maternal environmental, and cytoplasmic effects were investigated in a large Suffolk breeding scheme using a range of models involving different combinations of these effects with the direct additive genetic effect. Weights at 8 wk of age and at scanning (mean age 146 d) and ultrasonically measured muscle and fat depth were analyzed using an animal model on 55,683 (8-wk weight) and 28,947 (scanning traits) lamb records. Simple additive models always overestimated the heritability of all traits when compared to more complex models. The successive inclusion of maternal environmental, maternal genetic, and the covariance between direct and maternal additive effects in the model significantly improved the fit for almost all models and all traits, as indicated by a likelihood ratio test. Under the full model, the heritability of both weight traits was low (0.14 and 0.20 for 8-wk and scanning weight, respectively). The maternal additive and maternal environmental effects, as a proportion of the phenotypic variance, were similar (0.10 and 0.08 for 8-wk weight and 0.07 and 0.06 for scanning weight). The two scanning traits had higher heritabilities (0.29 and 0.27 for muscle depth and fat depth, respectively) with low levels of maternal genetic and maternal environmental variance. No evidence was found of a cytoplasmic effect on any of the traits studied under the full model. Breeding schemes for early growth and body traits in sheep should account for maternal effects in their genetic evaluations in order to improve their accuracy. The exact model to use will depend on the trait and individual circumstances of the scheme.

Allelic association and disease mapping

The application of allelic association to map genes for complex traits, particularly using high-density maps of single nucleotide polymorphisms in candidate regions, is an area of very active research. Here we present some aspects of the methodology and applications to both major gene mapping, which illustrates the effectiveness of the method, and oligogenes, where methods are still in flux and for which there have been relatively few successes to date. Several important considerations emerge, including the selection of the optimal metric for measuring association and the importance of modelling the decline in association with distance given the variability in association in a candidate region. The Malecot model of association with distance is shown to have a resolution of greater than 50 kilobases but the available evidence suggests that considerably higher resolution might be achieved with dense single nucleotide polymorphism (SNP) maps.