Gene-Gene Interactions among SPRYs for Nonsyndromic Cleft Lip/Palate.

Nonsyndromic cleft lip with or without cleft palate (NSCL/P) is a common birth defect with a complex genetic architecture. Gene-gene interactions have been increasingly regarded as contributing to the etiology of NSCL/P. A recent genome-wide association study revealed that a novel single-nucleotide polymorphism at SPRY1 in 4q28.1 showed a significant association with NSCL/P. In the current study, we explored the role of 3 SPRY genes in the etiology of NSCL/P by detecting gene-gene interactions: SPRY1, SPRY2, and SPRY4-with SPRY3 excluded due to its special location on the X chromosome. We selected markers in 3 SPRY genes to test for gene-gene interactions using 1,908 case-parent trios recruited from an international consortium established for a genome-wide association study of nonsyndromic oral clefts. As the trios came from populations with different ancestries, subgroup analyses were conducted among Europeans and Asians. Cordell's method based on conditional logistic regression models was applied to test for potential gene-gene interactions via the statistical package TRIO in R software. Gene-gene interaction analyses yielded 10 pairs of SNPs in Europeans and 6 pairs in Asians that achieved significance after Bonferroni correction. The significant interactions were confirmed in the 10,000-permutation tests (empirical P = 0.003 for the most significant interaction). The study identified gene-gene interactions among SPRY genes among 1,908 NSCL/P trios, which revealed the importance of potential gene-gene interactions for understanding the genetic architecture of NSCL/P. The evidence of gene-gene interactions in this study also provided clues for future biological studies to further investigate the mechanism of how SPRY genes participate in the development of NSCL/P.

Evidence of interaction between genes in the folate/homocysteine metabolic pathway in controlling risk of non-syndromic oral cleft.

OBJECTIVE: Little consistent evidence is available for the association between the risk of non-syndromic cleft lip with or without cleft palate (NSCL/P) and any of the individual genes in the folate/homocysteine metabolic pathway. We investigated the genes in the folate pathway to further clarify its potential influence on the risk of NSCL/P considering gene-gene (G×G) interaction. SUBJECTS AND METHODS: We selected markers in 18 genes from the pathway and applied Cordell's method to test for G×G interaction using 1,908 NSCL/P case-parent trios ascertained in an international consortium where a genomewide association study (GWAS) of oral clefts was conducted. RESULTS: We found intriguing signals among Asian and European ancestry groups for G×G interaction between markers in betaine-homocysteine methyltransferase gene (BHMT/BHMT2) and dimethylglycine dehydrogenase gene (DMGDH) attaining genomewide significance. In the pooled data, the top significant interaction was found between rs13158309 (BHMT) and rs10514154 (DMGDH, p = 1.45 × 10-12 ). CONCLUSIONS: Our study illustrated the importance of taking into account potential G×G interaction for genetic association analysis in NSCL/P, and this study suggested both BHMT/BHMT2 and DMGDH should be considered as candidate genes for NSCL/P in future studies.

GCKR and PPP1R3B identified as genome-wide significant loci for plasma lactate: the Atherosclerosis Risk in Communities (ARIC) study.

AIM: To investigate the genetic influence of circulating lactate level, a marker of oxidative capacity associated with diabetes. METHODS: We conducted a genome-wide association study of log-transformed plasma lactate levels in 6901 European-American participants in the Atherosclerosis Risk in Communities study. For regions that achieved genome-wide significance in European-American participants, we conducted candidate region analysis in African-American subjects and tested for interaction between metformin use and the index single nucleotide polymorphisms for plasma lactate in European-American subjects. RESULTS: The genome-wide association study in European-American subjects identified two genome-wide significant loci, GCKR (rs1260326, T allele ß=0.08; P=1.8×10(-47) ) and PPP1R3B/LOC157273 (rs9987289, A allele ß=0.06; P=1.6×10(-9) ). The index single nucleotide polymorphisms in these two loci explain 3.3% of the variance in log-transformed plasma lactate levels among the European-American subjects. In the African-American subjects, based on a region-significant threshold, the index single nucleotide polymorphism at GCKR was associated with plasma lactate but that at PPP1R3B/LOC157273 was not. Metformin use appeared to strengthen the association between the index single nucleotide polymorphism at PPP1R3B/LOC157273 and plasma lactate in European-American subjects (P for interaction=0.01). CONCLUSIONS: We identified GCKR and PPP1R3B/LOC157273 as two genome-wide significant loci of plasma lactate. Both loci are associated with other diabetes-related phenotypes. These findings increase our understanding of the genetic control of lactate metabolism.

A genome-wide analysis of the response to inhaled ß2-agonists in chronic obstructive pulmonary disease.

Short-acting ß2-agonist bronchodilators are the most common medications used in treating chronic obstructive pulmonary disease (COPD). Genetic variants determining bronchodilator responsiveness (BDR) in COPD have not been identified. We performed a genome-wide association study (GWAS) of BDR in 5789 current or former smokers with COPD in one African-American and four white populations. BDR was defined as the quantitative spirometric response to inhaled ß2-agonists. We combined results in a meta-analysis. In the meta-analysis, single-nucleotide polymorphisms (SNPs) in the genes KCNK1 (P=2.02 × 10(-7)) and KCNJ2 (P=1.79 × 10(-7)) were the top associations with BDR. Among African Americans, SNPs in CDH13 were significantly associated with BDR (P=5.1 × 10(-9)). A nominal association with CDH13 was identified in a gene-based analysis in all subjects. We identified suggestive association with BDR among COPD subjects for variants near two potassium channel genes (KCNK1 and KCNJ2). SNPs in CDH13 were significantly associated with BDR in African Americans.The Pharmacogenomics Journal advance online publication, 27 October 2015; doi:10.1038/tpj.2015.65.

IRF6 mutation screening in non-syndromic orofacial clefting: analysis of 1521 families.

Van der Woude syndrome (VWS) is an autosomal dominant malformation syndrome characterized by orofacial clefting (OFC) and lower lip pits. The clinical presentation of VWS is variable and can present as an isolated OFC, making it difficult to distinguish VWS cases from individuals with non-syndromic OFCs. About 70% of causal VWS mutations occur in IRF6, a gene that is also associated with non-syndromic OFCs. Screening for IRF6 mutations in apparently non-syndromic cases has been performed in several modestly sized cohorts with mixed results. In this study, we screened 1521 trios with presumed non-syndromic OFCs to determine the frequency of causal IRF6 mutations. We identified seven likely causal IRF6 mutations, although a posteriori review identified two misdiagnosed VWS families based on the presence of lip pits. We found no evidence for association between rare IRF6 polymorphisms and non-syndromic OFCs. We combined our results with other similar studies (totaling 2472 families) and conclude that causal IRF6 mutations are found in 0.24-0.44% of apparently non-syndromic OFC families. We suggest that clinical mutation screening for IRF6 be considered for certain family patterns such as families with mixed types of OFCs and/or autosomal dominant transmission.

Confirming genes influencing risk to cleft lip with/without cleft palate in a case-parent trio study.

A collection of 1,108 case-parent trios ascertained through an isolated, nonsyndromic cleft lip with or without cleft palate (CL/P) was used to replicate the findings from a genome-wide association study (GWAS) conducted by Beaty et al. (Nat Genet 42:525-529, 2010), where four different genes/regions were identified as influencing risk to CL/P. Tagging SNPs for 33 different genes were genotyped (1,269 SNPs). All four of the genes originally identified as showing genome-wide significance (IRF6, ABCA4 and MAF, plus the 8q24 region) were confirmed in this independent sample of trios (who were primarily of European and Southeast Asian ancestry). In addition, eight genes classified as 'second tier' hits in the original study (PAX7, THADA, COL8A1/FILIP1L, DCAF4L2, GADD45G, NTN1, RBFOX3 and FOXE1) showed evidence of linkage and association in this replication sample. Meta-analysis between the original GWAS trios and these replication trios showed PAX7, COL8A1/FILIP1L and NTN1 achieved genome-wide significance. Tests for gene-environment interaction between these 33 genes and maternal smoking found evidence for interaction with two additional genes: GRID2 and ELAVL2 among European mothers (who had a higher rate of smoking than Asian mothers). Formal tests for gene-gene interaction (epistasis) failed to show evidence of statistical interaction in any simple fashion. This study confirms that many different genes influence risk to CL/P.

Genetic association between human chitinases and lung function in COPD.

Two primary chitinases have been identified in humans--acid mammalian chitinase (AMCase) and chitotriosidase (CHIT1). Mammalian chitinases have been observed to affect the host's immune response. The aim of this study was to test for association between genetic variation in the chitinases and phenotypes related to chronic obstructive pulmonary disease (COPD). Polymorphisms in the chitinase genes were selected based on previous associations with respiratory diseases. Polymorphisms that were associated with lung function level or rate of decline in the Lung Health Study (LHS) cohort were analyzed for association with COPD affection status in four other COPD case-control populations. Chitinase activity and protein levels were also related to genotypes. In the caucasian LHS population, the baseline forced expiratory volume in one second (FEV(1)) was significantly different between the AA and GG genotypic groups of the AMCase rs3818822 polymorphism. Subjects with the GG genotype had higher AMCase protein and chitinase activity compared with AA homozygotes. For CHIT1 rs2494303, a significant association was observed between rate of decline in FEV(1) and the different genotypes. In the African American LHS population, CHIT1 rs2494303 and AMCase G339T genotypes were associated with rate of decline in FEV(1). Although a significant effect of chitinase gene alleles was found on lung function level and decline in the LHS, we were unable to replicate the associations with COPD affection status in the other COPD study groups.

Polymorphisms in IL10 are associated with total Immunoglobulin E levels and Schistosoma mansoni infection intensity in a Brazilian population.

Interleukin (IL)-10 is a regulatory cytokine of the helper T cell type 2 (TH2) pathway, which underlies both the host defense to helminthic infection and atopic diseases, including asthma. Although IL10 promoter polymorphisms are associated with increased atopy risk, IL10 variation has not been thoroughly explored in schistosomiasis-endemic populations. Three atopy-related IL10 promoter polymorphisms (rs1800896, rs1800871 and rs1800872), complemented by six tagging single-nucleotide polymorphisms (SNPs), were genotyped in 812 individuals in 318 nuclear families from a schistosomiasis-endemic area in Brazil. Associations between markers and total serum Immunoglobulin E (tIgE) levels, indicating non-specific activation of the TH2 pathway, and Schistosoma mansoni fecal egg counts, indicating burden of infection reflecting effectiveness of schistosomiasis host immunity, were performed using family-based transmission disequilibrium tests for quantitative traits (QTDTs). Alleles A, T and A at the three promoter SNPs rs1800896, rs1800871 and rs1800872 were associated with high tIgE levels in the same direction as in atopy populations (P=0.0008, 0.026 and 0.045), but not with egg counts. IL10 promoter polymorphisms appear to influence non-specific tIgE levels, but not schistosomiasis-specific immunity. The tagging SNP rs3024495 was associated with high S. mansoni egg counts (P=0.005), suggesting a novel locus in IL10 may influence clinically relevant burden of infection.

CD14, a key candidate gene associated with a specific immune response to cockroach.

BACKGROUND: Sensitization to cockroach allergen is one of the strongest predictors of asthma morbidity, especially among African Americans. OBJECTIVE: Our aims were to determine the genomic basis of cockroach sensitization and the specific response to cockroach antigen. METHODS: We investigated the Th1/Th2 cytokine profile of co-cultured plasmacytoid dendritic cells (pDCs) and CD4+ T cells and the 'transcript signature' of the immune response to cockroach antigen using high-throughput expression profiling of co-cultured cells. RESULTS: We observed significantly elevated levels of IL-13, IL-10, and TNF-alpha, but undetectable levels of IL-12p70 and IFN-alpha, when cultures were exposed to crude cockroach antigen. A significant difference was observed for IL-13 between cockroach-allergic and non-allergic individuals (P=0.039). Microarray analyses demonstrated a greater response at 48 h compared with 4 h, with 50 genes being uniquely expressed in cockroach antigen-treated cells, including CD14, S100A8, CCL8, and IFI44L. The increased CD14 expression was further observed in purified pDCs, human monocytic THP-1 cells, and the supernatant of co-cultured pDCs and CD4+ T cells on exposure to cockroach extract. Furthermore, the most differential expression of CD14 between cockroach allergy and non-cockroach allergy was only observed among individuals with the CC 'high-risk' genotype of the CD14-260C/T. Ingenuity Pathways Analysis analyses suggested the IFN signalling as the most significant canonical pathway. CONCLUSION: Our results suggest that these differentially expressed genes, particularly CD14, and genes in the IFN signalling pathway may be important candidates for further investigation of their role in the immune response to cockroach allergen.

TSLP polymorphisms are associated with asthma in a sex-specific fashion.

BACKGROUND: Single nucleotide polymorphisms (SNPs) in thymic stromal lymphopoietin (TSLP) have been associated with IgE (in girls) and asthma (in general). We sought to determine whether TSLP SNPs are associated with asthma in a sex-specific fashion. METHODS: We conducted regular and sex-stratified analyses of association between SNPs in TSLP and asthma in families of children with asthma in Costa Rica. Significant findings were replicated in whites and African-American participants in the Childhood Asthma Management Program, in African-Americans in the Genomic Research on Asthma in the African Diaspora study, in whites and Hispanics in the Children's Health Study, and in whites in the Framingham Heart Study (FHS). MAIN RESULTS: Two SNPs in TSLP (rs1837253 and rs2289276) were significantly associated with a reduced risk of asthma in combined analyses of all cohorts (P values of 2 × 10(-5) and 1 × 10(-5) , respectively). In a sex-stratified analysis, the T allele of rs1837253 was significantly associated with a reduced risk of asthma in males only (P = 3 × 10(-6) ). Alternately, the T allele of rs2289276 was significantly associated with a reduced risk of asthma in females only (P = 2 × 10(-4) ). Findings for rs2289276 were consistent in all cohorts except the FHS. CONCLUSIONS: TSLP variants are associated with asthma in a sex-specific fashion.

Association between IRF6 SNPs and oral clefts in West China.

Analyses of previous data have confirmed the contribution of the IRF6 gene to susceptibility to non- syndromic oral clefts (NSOC) in some populations. We tested for associations between the rs2013162, rs2235375, and rs2235371 polymorphisms in IRF6 and the risk of NSOC, using both case-parent trio and case-control designs on samples from western China. Our study group consisted of 332 persons with NSOC, their parents (289 mothers and 243 fathers for 206 complete trios for these three SNPs), and 174 control individuals. We found strong evidence of over- and under-transmission of the C allele (the Val allele) at rs2235371, and the C allele at rs2235375 in cleft case-parent trios (P = 0.013 and P = 0.000, respectively). There were significant differences in the frequency distributions of both genotypes and alleles when cases were compared with control infants at rs2235371 and rs2235375. Five specific haplotypes showed significant over- and under-transmission. These results further support a role for IRF6 variants in western Chinese populations.

Leptin receptor polymorphisms and lung function decline in COPD.

Only a fraction of all smokers develop chronic obstructive pulmonary disease (COPD), suggesting a large role for genetic susceptibility. The leptin receptor (LEPR) is present in human lung tissue and may play a role in COPD pathogenesis. The present study examined the association between genetic variants in the LEPR gene and lung function decline in COPD. In total, 429 European Americans were randomly selected from the National Heart Lung and Blood Institute Lung Health Study. 36 single nucleotide polymorphisms (SNPs) in LEPR were genotyped using the Illumina GoldenGate platform (Broad Institute, Cambridge, MA, USA). Mean annual decline in forced expiratory volume in 1 s % predicted over the 5-yr period was calculated using linear regression. Linear regression models were also used to adjust for potential confounders. In addition, in vivo expression of the receptor gene was assessed with immunohistochemistry on lungs from smoke-exposed inbred mice. We identified significant associations (p<0.05) between lung function decline and 21 SNPs. Haplotype analyses confirmed several of these associations seen with individual markers. Immunohistochemistry results in inbred mice strains support a potential role of LEPR in COPD pathogenesis. We identified genetic variants in the LEPR gene significantly associated with lung function decline in a population of smokers with COPD. Our results support a role for LEPR as a novel candidate gene for COPD.

Association of a nonsynonymous single-nucleotide polymorphism of matrix metalloproteinase 9 with giant cell arteritis.

OBJECTIVE: Giant cell arteritis (GCA) is the most common type of primary vasculitis. Matrix metalloproteinase 9 (MMP-9) is present in arterial lesions of GCA and may be involved in its pathogenesis. We investigated whether certain genotypes of 4 single-nucleotide polymorphisms (SNPs) of MMP-9 are overrepresented in patients with histologically confirmed GCA. METHODS: Four SNPs of MMP-9, rs3918242 in the promoter region and 3 nonsynonymous coding SNPs (rs3918252, rs17576, and rs2250889) were genotyped by polymerase chain reaction-restriction fragment length polymorphism analysis in 58 white patients for whom there was a clinical suspicion of GCA. Thirty of these patients had histologically confirmed GCA (group 1), and 28 patients had negative results of a temporal artery biopsy for GCA (group 2). Estimates of the genotype distributions of each of these SNPs in a white population were determined using publicly available genotype data for a panel of 23 individuals (group 3). RESULTS: Although 1 SNP was monomorphic in all 3 groups, we observed statistically significant differences in the genotype distributions for rs2250889 between group 1 and group 2 (P = 0.005) and between group 1 and group 3 (P = 0.009), but not between groups 2 and 3 (P = 0.965). CONCLUSION: These data derived from a sample of patients with GCA suggest that the G allele of MMP-9 polymorphism rs2250889 is overrepresented in patients with histologically confirmed GCA. Clearly, larger sample sizes will be necessary to confirm this suggestive association and better characterize a possible linkage disequilibrium structure among polymorphisms.

Genes at human chromosome 5q31.1 regulate delayed-type hypersensitivity responses associated with Leishmania chagasi infection.

Visceral leishmaniasis (VL) caused by Leishmania chagasi is endemic to northeast Brazil. A positive delayed-type hypersensitivity skin test response (DTH+) is a marker for acquired resistance to disease, clusters in families and may be genetically controlled. Twenty-three single nucleotide polymorphisms (SNPs) were genotyped in the cytokine 5q23.3-q31.1 region IRF1-IL5-IL13-IL4-IL9-LECT2-TGFBI in 102 families (323 DTH+; 190 DTH-; 123 VL individuals) from a VL endemic region in northeast Brazil. Data from 20 SNPs were analyzed for association with DTH+/- status and VL using family-based, stepwise conditional logistic regression analysis. Independent associations were observed between the DTH+ phenotype and markers in separate linkage disequilibrium blocks in LECT2 (OR 2.25; P=0.005; 95% CI=1.28-3.97) and TGFBI (OR 1.94; P=0.003; 95% CI=1.24-3.03). VL child/parent trios gave no evidence of association, but the DTH- phenotype was associated with SNP rs2070874 at IL4 (OR 3.14; P=0.006; 95% CI=1.38-7.14), and SNP rs30740 between LECT2 and TGFBI (OR 3.00; P=0.042; 95% CI=1.04-8.65). These results indicate several genes in the immune response gene cluster at 5q23.3-q31.1 influence outcomes of L. chagasi infection in this region of Brazil.

Invited commentary: two studies of genetic control of birth weight where large data sets were available.

Physical size at birth (primarily birth weight) is a key predictor of infant mortality and morbidity and may serve as a predictor of risk to chronic disease in adulthood. By use of birth records, it is sometimes possible to accumulate large, even massive, data sets that can permit analyses to separate genetic and environmental factors controlling variation in the complex phenotype, birth weight. Two studies of birth records, one from Norway and one from Lebanon, demonstrate how hospital- or registry-based data sets can be used to address fundamental questions about genetic control. The Medical Birth Registry of Norway has birth certificate data dating back to 1967 and allowed nuclear families to be reconstructed by linking children and their parents. Path models were used to estimate heritability of head circumference, along with birth weight and length for over 100,000 reconstructed families. A hospital-based study of birth records in Beirut, Lebanon, collected information on over 10,000 births, including sufficient numbers of marriages between first and second cousins to estimate inbreeding effects. Both of these studies confirm that birth weight is not simply due to the direct effects of the baby's genes but is a complex phenotype reflecting the effects of maternal genes and environments.

Analysis of candidate genes on chromosome 2 in oral cleft case-parent trios from three populations.

Isolated oral clefts, including cleft lip with/without cleft palate (CL/P) and cleft palate (CP), have a complex and heterogeneous etiology. Case-parent trios from three populations were used to study genes spanning chromosome 2, where single nucleotide polymorphic (SNP) markers were analyzed individually and as haplotypes. Case-parent trios from three populations (74 from Maryland, 64 from Singapore and 95 from Taiwan) were genotyped for 962 SNPs in 104 genes on chromosome 2, including two well-recognized candidate genes: TGFA and SATB2. Individual SNPs and haplotypes (in sliding windows of 2-5 SNPs) were used to test for linkage and disequilibrium separately in CL/P and CP trios. A novel candidate gene (ZNF533) showed consistent evidence of linkage and disequilibrium in all three populations for both CL/P and CP. SNPs in key regions of ZNF533 showed considerable variability in estimated genotypic odds ratios and their significance, suggesting allelic heterogeneity. Haplotype frequencies for regions of ZNF533 were estimated and used to partition genetic variance into among-and within-population components. Wright's fixation index, a measure of genetic diversity, showed little difference between Singapore and Taiwan compared with Maryland. The tensin-1 gene (TNS1) also showed evidence of linkage and disequilibrium among both CL/P and CP trios in all three populations, albeit at a lower level of significance. Additional genes (VAX2, GLI2, ZHFX1B on 2p; WNT6-WNT10A and COL4A3-COL4A4 on 2q) showed consistent evidence of linkage and disequilibrium only among CL/P trios in all three populations, and TGFA showed significant evidence in two of three populations.

High throughput SNP and expression analyses of candidate genes for non-syndromic oral clefts.

BACKGROUND: Recent work suggests that multiple genes and several environmental risk factors influence risk for non-syndromic oral clefts, one of the most common birth defects in humans. Advances in high-throughput genotyping technology now make it possible to test multiple markers in many candidate genes simultaneously. METHODS: We present findings from family based association tests of single nucleotide polymorphism (SNP) markers in 64 candidate genes genotyped using the BeadArray approach in 58 case-parent trios from Maryland (USA) to illustrate how multiple markers in multiple genes can be analysed. To assess whether these genes were expressed in human craniofacial structures relevant to palate and lip development, we also analysed data from the Craniofacial and Oral Gene Expression Network (COGENE) consortium, and searched public databases for expression profiles of these genes. RESULTS: Thirteen candidate genes showed significant evidence of linkage in the presence of disequilibrium, and ten of these were found to be expressed in relevant embryonic tissues: SP100, MLPH, HDAC4, LEF1, C6orf105, CD44, ALX4, ZNF202, CRHR1, and MAPT. Three other genes showing statistical evidence (ADH1C, SCN3B, and IMP5) were not expressed in the embryonic tissues examined here. CONCLUSIONS: This approach demonstrates how statistical evidence on large numbers of SNP markers typed in case-parent trios can be combined with expression data to identify candidate genes for complex disorders. Many of the genes reported here have not been previously studied as candidates for oral clefts and warrant further investigation.

Variants in the gene encoding C3 are associated with asthma and related phenotypes among African Caribbean families.

Proinflammatory and immunoregulatory products from C3 play a major role in phagocytosis, respiratory burst, and airways inflammation. C3 is critical in adaptive immunity; studies in mice deficient in C3 demonstrate that features of asthma are significantly attenuated in the absence of C3. To test the hypothesis that the C3 gene on chromosome 19p13.3-p13.2 contains variants associated with asthma and related phenotypes, we genotyped 25 single nucleotide polymorphism (SNP) markers distributed at intervals of approximately 1.9 kb within the C3 gene in 852 African Caribbean subjects from 125 nuclear and extended pedigrees. We used the multiallelic test in the family-based association test program to examine sliding windows comprised of 2-6 SNPs. A five-SNP window between markers rs10402876 and rs366510 provided strongest evidence for linkage in the presence of linkage disequilibrium for asthma, high log[total IgE], and high log[IL-13]/[log[IFN-gamma] in terms of global P-values (P = 0.00027, 0.00013, and 0.003, respectively). A three-SNP haplotype GGC for the first three of these markers showed best overall significance for the three phenotypes (P = 0.003, 0.007, 0.018, respectively) considering haplotype-specific tests. Taken together, these results implicate the C3 gene as a priority candidate controlling risk for asthma and allergic disease in this population of African descent.

Haplotype diversity in 11 candidate genes across four populations.

Analysis of haplotypes based on multiple single-nucleotide polymorphisms (SNP) is becoming common for both candidate gene and fine-mapping studies. Before embarking on studies of haplotypes from genetically distinct populations, however, it is important to consider variation both in linkage disequilibrium (LD) and in haplotype frequencies within and across populations, as both vary. Such diversity will influence the choice of "tagging" SNPs for candidate gene or whole-genome association studies because some markers will not be polymorphic in all samples and some haplotypes will be poorly represented or completely absent. Here we analyze 11 genes, originally chosen as candidate genes for oral clefts, where multiple markers were genotyped on individuals from four populations. Estimated haplotype frequencies, measures of pairwise LD, and genetic diversity were computed for 135 European-Americans, 57 Chinese-Singaporeans, 45 Malay-Singaporeans, and 46 Indian-Singaporeans. Patterns of pairwise LD were compared across these four populations and haplotype frequencies were used to assess genetic variation. Although these populations are fairly similar in allele frequencies and overall patterns of LD, both haplotype frequencies and genetic diversity varied significantly across populations. Such haplotype diversity has implications for designing studies of association involving samples from genetically distinct populations.