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.

Group-based pharmacogenetic prediction: is it feasible and do current NHS England ethnic classifications provide appropriate data?

Inter-individual variation of drug metabolising enzymes (DMEs) leads to variable efficacy of many drugs and even adverse drug responses. Consequently, it would be desirable to test variants of many DMEs before drug treatment. Inter-ethnic differences in frequency mean that the choice of SNPs to test may vary across population groups. Here we examine the utility of testing representative groups as a way of assessing what variants might be tested. We show that publicly available population information is potentially useful for determining loci for pre-treatment genetic testing, and for determining the most prevalent risk haplotypes in defined groups. However, we also show that the NHS England classifications have limitations for grouping for these purposes, in particular for people of African descent. We conclude: (1) genotyping of hospital patients and people from the hospital catchment area confers no advantage over using samples from appropriate existing ethnic group collections or publicly available data, (2) given the current NHS England Black African grouping, a decision as to whether to test, would have to apply to all patients of recent Black African ancestry to cover reported risk alleles and (3) the current scarcity of available genome and drug effect data from Africans is a problem for both testing and treatment decisions.

Adaptation to milking agropastoralism in Chilean goat herders and nutritional benefit of lactase persistence.

The genetic trait of lactase persistence (LP) evolved as an adaptation to milking pastoralism in the Old World and is a well-known example of positive natural selection in humans. However, the specific mechanisms conferring this selective advantage are unknown. To understand the relationship between milk drinking, LP, growth, reproduction, and survival, communities of the Coquimbo Region in Chile, with recent adoption of milking agropastoralism, were used as a model population. DNA samples and data on stature, reproduction, and diet were collected from 451 participants. Lactose tolerance tests were done on 41 of them. The European -13,910*T (rs4988235) was the only LP causative variant found, showing strong association (99.6%) with LP phenotype. Models of associations of inferred LP status and milk consumption, with fertility, mortality, height, and weight were adjusted with measures of ancestry and relatedness to control for population structure. Although we found no statistically significant effect of LP on fertility, a significant effect (P = 0.002) was observed of LP on body mass index (BMI) in males and of BMI on fertility (P = 0.003). These results fail to support a causal relationship between LP and fertility yet suggest the idea of a nutritional advantage of LP. Furthermore, the proportion of European ancestry around the genetic region of -13,910*T is significantly higher (P = 0.008) than the proportion of European ancestry genome-wide, providing evidence of recent positive selection since European-Amerindian admixture. This signature was absent in nonpastoralist Latin American populations, supporting the hypothesis of specific adaptation to milking agropastoralism in the Coquimbo communities.

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.

In Vitro Functional Analyses of Infrequent Nucleotide Variants in the Lactase Enhancer Reveal Different Molecular Routes to Increased Lactase Promoter Activity and Lactase Persistence.

The genetic trait that allows intestinal lactase to persist into adulthood in some 35% of humans worldwide operates at the level of transcription, the effect being caused by cis-acting nucleotide changes upstream of the lactase gene (LCT). A single nucleotide substitution, -13910 C>T, the first causal variant to be identified, accounts for lactase persistence over most of Europe. Located in a region shown to have enhancer function in vitro, it causes increased activity of the LCT promoter in Caco-2 cells, and altered transcription factor binding. Three other variants in close proximity, -13907 C>G, -13915 T>C and -14010 G>C, were later shown to behave in a similar manner. Here, we study four further candidate functional variants. Two, -14009 T>G and -14011 C>T, adjacent to the well-studied -14010 G>C variant, also have a clear effect on promoter activity upregulation as assessed by transfection assays, but notably are involved in different molecular interactions. The results for the two other variants (-14028 T>C, -13779 G>C) were suggestive of function, -14028*C showing a clear change in transcription factor binding, but no obvious effect in transfections, while -13779*G showed greater effect in transfections but less on transcription factor binding. Each of the four variants arose on independent haplotypic backgrounds with different geographic distribution.

Diversity of lactase persistence in African milk drinkers.

The genetic trait of lactase persistence is attributable to allelic variants in an enhancer region upstream of the lactase gene, LCT. To date, five different functional alleles, -13910*T, -13907*G, -13915*G, -14009*G and -14010*C, have been identified. The co-occurrence of several of these alleles in Ethiopian lactose digesters leads to a pattern of sequence diversity characteristic of a 'soft selective sweep'. Here we hypothesise that throughout Africa, where multiple functional alleles co-exist, the enhancer diversity will be greater in groups who are traditional milk drinkers than in non-milk drinkers, as the result of this sort of parallel selection. Samples from 23 distinct groups from 10 different countries were examined. Each group was classified 'Yes 'or 'No' for milk-drinking, and ethnicity, language spoken and geographic location were recorded. Predicted lactase persistence frequency and enhancer diversity were, as hypothesised, higher in the milk drinkers than the non-milk-drinkers, but this was almost entirely accounted for by the Afro-Asiatic language speaking peoples of east Africa. The other groups, including the 'Nilo-Saharan language speaking' milk-drinkers, show lower frequencies of LP and lower diversity, and there was a north-east to south-west decline in overall diversity. Amongst the Afro-Asiatic (Cushitic) language speaking Oromo, however, the geographic cline was not evident and the southern pastoralist Borana showed much higher LP frequency and enhancer diversity than the other groups. Together these results reflect the effects of parallel selection, the stochastic processes of the occurrence and spread of the mutations, and time depth of milk drinking tradition.

Genetic variation underlying common hereditary hyperbilirubinaemia (Gilbert's syndrome) and respiratory health in the 1946 British birth cohort.

BACKGROUND & AIMS: Bilirubin has potent antioxidant properties in vitro and raised serum levels have been associated with lower rates of respiratory disease. The enzyme uridine diphosphate glucuronosyltransferase polypeptide 1A1 (UGT1A1) is solely responsible for clearing bilirubin from the blood and homozygosity for seven thymine-adenine (TA) repeats in the TATA box regulatory element of the UGT1A1 gene underlies a mild hereditary unconjugated hyperbilirubinaemia (Gilbert's syndrome). Our aim was to investigate whether this genetic variation is associated with differences in respiratory health. METHODS: The relationship between the promoter genotype underlying Gilbert's syndrome (UGT1A1 rs8175347 [TA]7/7) and respiratory outcomes assessed at ages 43, 53, and 60-64 were examined in 2190 members of the 1946 British birth cohort. RESULTS: The (TA)7/7 genotype, present in 9% of the cohort, was associated with higher forced expiratory volume (FEV1) and forced vital capacity (FVC). The relationship was strongest for heavy smokers (⩾20 cigarettes per day) at age 53 with mean FEV1 409 ml higher (191 to 627; p<0.001) and mean FVC 530 ml higher (95% CI 262-798; p<0.001) for UGT1A1 (TA)7/7 Gilbert's syndrome participants than for all others, indicating a protection from the pulmonary consequences of heavy smoking. The odds of respiratory disease (chronic obstructive pulmonary disease, self-reported asthma, or prescription of respiratory drugs) were half in those with Gilbert's syndrome genotype (odds ratio 0.49 [95% CI 0.39-0.74]; p<0.001) compared to those without this genotype. CONCLUSIONS: Genetically raised unconjugated serum bilirubin is associated with higher adult respiratory function and protection from respiratory disease.

Diversity of lactase persistence alleles in Ethiopia: signature of a soft selective sweep.

The persistent expression of lactase into adulthood in humans is a recent genetic adaptation that allows the consumption of milk from other mammals after weaning. In Europe, a single allele (-13910(∗)T, rs4988235) in an upstream region that acts as an enhancer to the expression of the lactase gene LCT is responsible for lactase persistence and appears to have been under strong directional selection in the last 5,000 years, evidenced by the widespread occurrence of this allele on an extended haplotype. In Africa and the Middle East, the situation is more complicated and at least three other alleles (-13907(∗)G, rs41525747; -13915(∗)G, rs41380347; -14010(∗)C, rs145946881) in the same LCT enhancer region can cause continued lactase expression. Here we examine the LCT enhancer sequence in a large lactose-tolerance-tested Ethiopian cohort of more than 350 individuals. We show that a further SNP, -14009T>G (ss 820486563), is significantly associated with lactose-digester status, and in vitro functional tests confirm that the -14009(∗)G allele also increases expression of an LCT promoter construct. The derived alleles in the LCT enhancer region are spread through several ethnic groups, and we report a greater genetic diversity in lactose digesters than in nondigesters. By examining flanking markers to control for the effects of mutation and demography, we further describe, from empirical evidence, the signature of a soft selective sweep.

MUC5AC and inflammatory mediators associated with respiratory outcomes in the British 1946 birth cohort.

BACKGROUND AND OBJECTIVE: Dysregulation of respiratory mucins, MUC5AC in particular, has been implicated in respiratory disease and MUC5AC expression is up-regulated in response to environmental challenges and inflammatory mediators. The aim of this study was to examine the effect of genetic variation on susceptibility to common respiratory conditions. METHODS: The association of MUC5AC and the closely linked genes MUC2 and MUC5B with respiratory outcomes was tested in the MRC National Survey of Health and Development, a longitudinal birth cohort of men and women born in 1946. Also examined were the functional variants of the genes encoding inflammatory mediators, IL13, IL1B, IL1RN, TNFA and ERBB1, for which there is a likely influence on MUC5AC expression and were explored potential gene-gene interactions with these inflammatory mediators. RESULTS: Statistically significant associations between the 3'ter MUC5AC simple nucleotide polymorphism (SNP) rs1132440 and various non-independent respiratory outcomes (bronchitis, wheeze, asthma, hay fever) were reported while the adjacent loci show slight (but largely non-statistically significant) differences, presumably reflective of linkage disequilibrium (allelic association) across the region. A novel association between bronchitis and a non-synonymous functional ERBB1 SNP, rs2227983 (aka epidermal growth factor receptor:R497K, R521K) is also reported and evidence presented of interaction between MUC5AC and ERBB1 and between MUC5AC and IL1RN with respect to bronchitis. The ERBB1 result suggests a clear mechanism for a biological interaction in which the allelic variants of epidermal growth factor receptor differentially affect mucin expression. CONCLUSIONS: The MUC5AC association and the interactions with inflammatory mediators suggest that genetically determined differences in MUC5AC expression alter susceptibility to respiratory disease.

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.

Mucin methods: genes encoding mucins and their genetic variation with a focus on gel-forming mucins.

Mucin genes encode the polypeptide backbone of the mucin glycoproteins which are expressed on all epithelial surfaces and are major constituents of the mucus layer. Mucins are, thus, expressed at the interface between the external and the internal environment of the organism, and represent the first line of defence of our body. These genes often have an extensive region of repetitive exonic sequence which codes for the heavily glycosylated domain, whose roles include bacterial interactions and gel hydration. This region shows, in several of the genes, considerable inter-individual variation in repeat number and sequence. Because of their site of expression and their high variability in this important domain, mucin genes are good candidates for conferring differences in genetic susceptibility to multifactorial epithelial and inflammatory disease. However, progress in characterizing the genes has been considerably slower than the rest of the genome because of their size and the GC-rich content of the large, repetitive variable region. Some of the issues relating to the study of these genes are discussed in this chapter. In addition, methods and approaches that have been used successfully are described.

Herders of Indian and European cattle share their predominant allele for lactase persistence.

Milk consumption and lactose digestion after weaning are exclusively human traits made possible by the continued production of the enzyme lactase in adulthood. Multiple independent mutations in a 100-bp region--part of an enhancer--approximately 14-kb upstream of the LCT gene are associated with this trait in Europeans and pastoralists from Saudi Arabia and Africa. However, a single mutation of purported western Eurasian origin accounts for much of observed lactase persistence outside Africa. Given the high levels of present-day milk consumption in India, together with archaeological and genetic evidence for the independent domestication of cattle in the Indus valley roughly 7,000 years ago, we sought to determine whether lactase persistence has evolved independently in the subcontinent. Here, we present the results of the first comprehensive survey of the LCT enhancer region in south Asia. Having genotyped 2,284 DNA samples from across the Indian subcontinent, we find that the previously described west Eurasian -13910 C>T mutation accounts for nearly all the genetic variation we observed in the 400- to 700-bp LCT regulatory region that we sequenced. Geography is a significant predictor of -13910*T allele frequency, and consistent with other genomic loci, its distribution in India follows a general northwest to southeast declining pattern, although frequencies among certain neighboring populations vary substantially. We confirm that the mutation is identical by descent to the European allele and is associated with the same>1 Mb extended haplotype in both populations.

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.

Mucin variable number tandem repeat polymorphisms and severity of cystic fibrosis lung disease: significant association with MUC5AC.

Variability in cystic fibrosis (CF) lung disease is partially due to non-CFTR genetic modifiers. Mucin genes are very polymorphic, and mucins play a key role in the pathogenesis of CF lung disease; therefore, mucin genes are strong candidates as genetic modifiers. DNA from CF patients recruited for extremes of lung phenotype was analyzed by Southern blot or PCR to define variable number tandem repeat (VNTR) length polymorphisms for MUC1, MUC2, MUC5AC, and MUC7. VNTR length polymorphisms were tested for association with lung disease severity and for linkage disequilibrium (LD) with flanking single nucleotide polymorphisms (SNPs). No strong associations were found for MUC1, MUC2, or MUC7. A significant association was found between the overall distribution of MUC5AC VNTR length and CF lung disease severity (p = 0.025; n = 468 patients); plus, there was robust association of the specific 6.4 kb HinfI VNTR fragment with severity of lung disease (p = 6.2×10(-4) after Bonferroni correction). There was strong LD between MUC5AC VNTR length modes and flanking SNPs. The severity-associated 6.4 kb VNTR allele of MUC5AC was confirmed to be genetically distinct from the 6.3 kb allele, as it showed significantly stronger association with nearby SNPs. These data provide detailed respiratory mucin gene VNTR allele distributions in CF patients. Our data also show a novel link between the MUC5AC 6.4 kb VNTR allele and severity of CF lung disease. The LD pattern with surrounding SNPs suggests that the 6.4 kb allele contains, or is linked to, important functional genetic variation.

The -14010*C variant associated with lactase persistence is located between an Oct-1 and HNF1α binding site and increases lactase promoter activity.

In most people worldwide intestinal lactase expression declines in childhood. In many others, particularly in Europeans, lactase expression persists into adult life. The lactase persistence phenotype is in Europe associated with the -13910*T single nucleotide variant located 13,910 bp upstream the lactase gene in an enhancer region that affects lactase promoter activity. This variant falls in an Oct-1 binding site and shows greater Oct-1 binding than the ancestral variant and increases enhancer activity. Several other variants have been identified very close to the -13910 position, which are associated with lactase persistence in the Middle East and Africa. One of them, the -14010*C, is associated with lactase persistence in Africa. Here we show by deletion analysis that the -14010 position is located in a 144 bp region that reduces the enhancer activity. In transfections the -14010*C allele shows a stronger enhancer effect than the ancestral -4010*G allele. Binding sites for Oct-1 and HNF1α surrounding the -14010 position were identified by gel shift assays, which indicated that -14010*C has greater binding affinity to Oct-1 than -14010*G.

Evolution of lactase persistence: an example of human niche construction.

Niche construction is the process by which organisms construct important components of their local environment in ways that introduce novel selection pressures. Lactase persistence is one of the clearest examples of niche construction in humans. Lactase is the enzyme responsible for the digestion of the milk sugar lactose and its production decreases after the weaning phase in most mammals, including most humans. Some humans, however, continue to produce lactase throughout adulthood, a trait known as lactase persistence. In European populations, a single mutation (-13910*T) explains the distribution of the phenotype, whereas several mutations are associated with it in Africa and the Middle East. Current estimates for the age of lactase persistence-associated alleles bracket those for the origins of animal domestication and the culturally transmitted practice of dairying. We report new data on the distribution of -13910*T and summarize genetic studies on the diversity of lactase persistence worldwide. We review relevant archaeological data and describe three simulation studies that have shed light on the evolution of this trait in Europe. These studies illustrate how genetic and archaeological information can be integrated to bring new insights to the origins and spread of lactase persistence. Finally, we discuss possible improvements to these models.

Serum bilirubin and risk of respiratory disease and death.

CONTEXT: Serum total bilirubin levels in healthy patients reflect genetic and environmental factors that could influence the risk of developing respiratory disease. OBJECTIVE: To examine the relationship between bilirubin levels and respiratory disease. DESIGN, SETTING, AND PARTICIPANTS: Cohort study among 504,206 adults from a UK primary care research database (the Health Improvement Network) with serum bilirubin levels recorded but no evidence of hepatobiliary or hemolytic disease. Data were recorded between January 1988 and December 2008. MAIN OUTCOME MEASURES: Incidence of chronic obstructive pulmonary disease (COPD), lung cancer, and all-cause mortality. RESULTS: Median bilirubin levels were 0.64 mg/dL (interquartile range, 0.47-0.88 mg/dL) in men and 0.53 mg/dL (interquartile range, 0.41-0.70 mg/dL) in women. There were 1341 cases of lung cancer, 5863 cases of COPD, and 23,103 deaths, with incidence rates of 2.5, 11.9, and 42.5 per 10,000 person-years, respectively. The incidence of lung cancer per 10,000 person-years in men was 5.0 (95% confidence interval [CI], 4.2-6.0) in the first decile category of bilirubin compared with 3.0 (95% CI, 2.3-3.8) in the fifth decile. The corresponding incidences for COPD in men were 19.5 (95% CI,17.7-21.4) and 14.4 (95% CI, 12.7-16.2). The mortality rates per 10,000 person-years in men were 51.3 (95% CI, 48.5-54.2) in the first decile category compared with 38.1 (95% CI, 35.5-40.8) in the fifth decile. The associations were similar for women. After adjusting for other important health indicators, regression estimates for incidence rate of lung cancer per 0.1-mg/dL increase in bilirubin level were an 8% decrease (95% CI, 5%-11%) for men and an 11% decrease (95% CI, 7%-14%) for women. The regression estimate for COPD in men per 0.1-mg/dL increase in bilirubin level was a 6% decrease (95% CI, 5%-7%) and for mortality in men was a 3% decrease (95% CI, 2%-3%) after accounting for other health indicators. The results for COPD and mortality in women were very similar. CONCLUSION: Among patients with normal-range bilirubin levels in primary care practices, relatively higher levels of bilirubin were associated with a lower risk of respiratory disease and all-cause mortality.

Prevalence of clinically relevant UGT1A alleles and haplotypes in African populations.

Variation of a short (TA)(n) repeat sequence (rs8175347) covering the TATA box of UGT1A1 (UDP-glucuronosyltransferase1A1) is associated with hyperbilirubinaemia (Gilbert's syndrome) and adverse drug reactions, and is used for dosage advice for irinotecan. Several reports indicate that the low-activity (risk) alleles ((TA)(7) and (TA)(8) )) are very frequent in Africans but the patterns of association with other variants in the UGT1A gene complex that may modulate these responses are not well known. rs8175347 and two other clinically relevant UGT1A variants (rs11692021 and rs10929302) were assayed in 2616 people from Europe and Africa. Low-activity (TA)(n) alleles frequencies were highest in equatorial Africa, (TA)(7,) being the most common in Cameroon, Ghana, southern Sudan, and in Ethiopian Anuak. Haplotypic diversity was also greatest in equatorial Africa, but in Ethiopia was very variable across ethnic groups. Resequencing of the promoter of a sample subset revealed no novel variations, but rs34547608 and rs887829 were typed and shown to be tightly associated with (TA)(n) . Our results illustrate the need for investigation of the effect of UGT1A variants other than (TA)(n) on the risk of irinotecan toxicity, as well as hyperbilirubinaemia due to hemolytic anaemia or human immunodeficiency virus protease inhibitors, so that appropriate pharmacogenetic advice can be given.