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.