Loss of function mutations in HARS cause a spectrum of inherited peripheral neuropathies.

Inherited peripheral neuropathies are a genetically heterogeneous group of disorders characterized by distal muscle weakness and sensory loss. Mutations in genes encoding aminoacyl-tRNA synthetases have been implicated in peripheral neuropathies, suggesting that these tRNA charging enzymes are uniquely important for the peripheral nerve. Recently, a mutation in histidyl-tRNA synthetase (HARS) was identified in a single patient with a late-onset, sensory-predominant peripheral neuropathy; however, the genetic evidence was lacking, making the significance of the finding unclear. Here, we present clinical, genetic, and functional data that implicate HARS mutations in inherited peripheral neuropathies. The associated phenotypic spectrum is broad and encompasses axonal and demyelinating motor and sensory neuropathies, including four young patients presenting with pure motor axonal neuropathy. Genome-wide linkage studies in combination with whole-exome and conventional sequencing revealed four distinct and previously unreported heterozygous HARS mutations segregating with autosomal dominant peripheral neuropathy in four unrelated families (p.Thr132Ile, p.Pro134His, p.Asp175Glu and p.Asp364Tyr). All mutations cause a loss of function in yeast complementation assays, and p.Asp364Tyr is dominantly neurotoxic in a Caenorhabditis elegans model. This study demonstrates the role of HARS mutations in peripheral neuropathy and expands the genetic and clinical spectrum of aminoacyl-tRNA synthetase-related human disease.

A mild form of SLC29A3 disorder: a frameshift deletion leads to the paradoxical translation of an otherwise noncoding mRNA splice variant

We investigated two siblings with granulomatous histiocytosis prominent in the nasal area, mimicking rhinoscleroma and Rosai-Dorfman syndrome. Genome-wide linkage analysis and whole-exome sequencing identified a homozygous frameshift deletion in SLC29A3, which encodes human equilibrative nucleoside transporter-3 (hENT3). Germline mutations in SLC29A3 have been reported in rare patients with a wide range of overlapping clinical features and inherited disorders including H syndrome, pigmented hypertrichosis with insulin-dependent diabetes, and Faisalabad histiocytosis. With the exception of insulin-dependent diabetes and mild finger and toe contractures in one sibling, the two patients with nasal granulomatous histiocytosis studied here displayed none of the many SLC29A3-associated phenotypes. This mild clinical phenotype probably results from a remarkable genetic mechanism. The SLC29A3 frameshift deletion prevents the expression of the normally coding transcripts. It instead leads to the translation, expression, and function of an otherwise noncoding, out-of-frame mRNA splice variant lacking exon 3 that is eliminated by nonsense-mediated mRNA decay (NMD) in healthy individuals. The mutated isoform differs from the wild-type hENT3 by the modification of 20 residues in exon 2 and the removal of another 28 amino acids in exon 3, which include the second transmembrane domain. As a result, this new isoform displays some functional activity. This mechanism probably accounts for the narrow and mild clinical phenotype of the patients. This study highlights the"rescue" role played by a normally noncoding mRNA splice variant of SLC29A3, uncovering a new mechanism by which frameshift mutations can be hypomorphic.

Little evidence for association between the TGFBR1*6A variant and colorectal cancer: a family-based association study on non-syndromic family members from Australia and Spain.

Genome-wide linkage studies have identified the 9q22 chromosomal region as linked with colorectal cancer (CRC) predisposition. A candidate gene in this region is transforming growth factor beta receptor 1 (TGFBR1). Investigation of TGFBR1 has focused on the common genetic variant rs11466445, a short exonic deletion of nine base pairs which results in truncation of a stretch of nine alanine residues to six alanine residues in the gene product. While the six alanine (*6A) allele has been reported to be associated with increased risk of CRC in some population based study groups this association remains the subject of robust debate. To date, reports have been limited to population-based case-control association studies, or case-control studies of CRC families selecting one affected individual per family. No study has yet taken advantage of all the genetic information provided by multiplex CRC families. Methods: We have tested for an association between rs11466445 and risk of CRC using several family-based statistical tests in a new study group comprising members of non-syndromic high risk CRC families sourced from three familial cancer centres, two in Australia and one in Spain. Results: We report a finding of a nominally significant result using the pedigree-based association test approach (PBAT; p = 0.028), while other family-based tests were non-significant, but with a p-value < 0.10 in each instance. These other tests included the Generalised Disequilibrium Test (GDT; p = 0.085), parent of origin GDT Generalised Disequilibrium Test (GDT-PO; p = 0.081) and empirical Family-Based Association Test (FBAT; p = 0.096, additive model). Related-person case-control testing using the 'More Powerful' Quasi-Likelihood Score Test did not provide any evidence for association (M-QL5; p = 0.41). Conclusions: After conservatively taking into account considerations for multiple hypothesis testing, we find little evidence for an association between the TGFBR1*6A allele and CRC risk in these families. The weak support for an increase in risk in CRC predisposed families is in agreement with recent meta-analyses of case-control studies, which estimate only a modest increase in sporadic CRC risk among 6*A allele carriers.

MAST3 : facteur de risque génétique aux maladies inflammatoires de l’intestin et modulateur d’inflammation

La maladie de Crohn (MC) et la colite ulcéreuse (CU) sont des maladies inflammatoires chroniques du tube digestif qu’on regroupe sous le terme maladies inflammatoires de l’intestin (MII). Les mécanismes moléculaires menant au développement des MII ne sont pas entièrement connus, mais des études génétiques et fonctionnelles ont permis de mettre en évidence des interactions entre des prédispositions génétiques et des facteurs environnementaux - notamment la flore intestinale – qui contribuent au développement d’une dérégulation de la réponse immunitaire menant à l’inflammation de la muqueuse intestinale. Des études d’association pangénomiques et ciblées ont permis d’identifier plusieurs gènes de susceptibilité aux MII mais les estimations de la contribution de ces gènes à l’héritabilité suggèrent que plusieurs gènes restent à découvrir. Certains d’entre eux peuvent se trouver dans les régions identifiées par des études de liaison génétique. L’objectif de mon projet de doctorat était d’identifier un ou des facteurs de risque génétique dans la région chromosomale 19p (identifiée comme région de liaison IBD6) et de le/les caractériser au niveau fonctionnel. Nous avons d’abord entrepris une cartographie d’association de la région 19p. À la suite du génotypage successif de deux cohortes indépendantes, nous avons identifié un SNP intronique et quatre SNP codants dont un non-synonyme, rs8108738, tous localisés dans le gène microtubule associated serine threonine kinase gene-3 (MAST3) et associés aux MII. Peu d’information fonctionnelle sur MAST3 était disponible. Par contre MAST2, une protéine encodée par un gène de la même famille, régule l’activité du facteur de transcription inflammatoire NF-kappaB. Nous avons confirmé l’implication de MAST3 dans l’activité de NF-kappaB via un knockdown de MAST3 et des essais gène-rapporteur. Pour poursuivre la caractérisation fonctionnelle de MAST3, nous avons choisi une approche non ciblée pour étudier les effets de la variation des niveaux d’expression de MAST3 sur la cellule. C’est-à-dire que nous avons créé un 1er modèle cellulaire de surexpression du gène MAST3 dans les cellules HEK293 et analysé l’expression pangénomique endogène. La validation de l’expression génique dans un 2e modèle cellulaire de knockdown et de type cellulaire différent (THP1), nous a permis d’identifier et de contrer les effets non-spécifiques dus aux niveaux non-physiologiques. Notre étude d’expression a mené à l’identification d’un groupe de gènes dont l’expression est régulée par MAST3. Ces gènes sont majoritairement impliqués dans des fonctions immunitaires (cytokines pro-inflammatoires, régulateurs de NF-kappaB, migration cellulaire, etc.) et une forte proportion est régulée par NF-kappaB. Nous avons évalué l’importance du groupe de gènes régulés par MAST3 dans la présentation clinique des MII à travers des études d’expression dans des biopsies intestinales de patients atteints de CU. Nous avons constaté que l’expression de ces gènes est significativement supérieure dans les régions enflammées par rapport aux régions saines de la muqueuse intestinale des patients atteints de CU. Globalement, les résultats de nos études suggèrent que le facteur de risque aux MII MAST3 agit via la voie du facteur de transcription NF-kappaB pour influencer l’expression d’un groupe de gènes impliqués dans l’inflammation intestinale typique des MII. Chaque étude génétique sur les MII a le potentiel d’orienter les recherches fonctionnelles vers de nouvelles voies biologiques causales. Le dévoilement des mécanismes moléculaires sous-jacents à ces voies permet d’augmenter les connaissances sur le développement de ces maladies vers une compréhension plus complète de la pathogenèse qui permettra d’optimiser le diagnostic et le traitement de ces maladies.

Progress in the genetics of common obesity and type 2 diabetes

The prevalence of obesity and diabetes, which are heritable traits that arise from the interactions of multiple genes and lifestyle factors, continues to rise worldwide, causing serious health problems and imposing a substantial economic burden on societies. For the past 15 years, candidate gene and genome-wide linkage studies have been the main genetic epidemiological approaches to identify genetic loci for obesity and diabetes, yet progress has been slow and success limited. The genome-wide association approach, which has become available in recent years, has dramatically changed the pace of gene discoveries. Genome-wide association is a hypothesis-generating approach that aims to identify new loci associated with the disease or trait of interest. So far, three waves of large-scale genome-wide association studies have identified 19 loci for common obesity and 18 for common type 2 diabetes. Although the combined contribution of these loci to the variation in obesity and diabetes risk is small and their predictive value is typically low, these recently identified loci are set to substantially improve our insights into the pathophysiology of obesity and diabetes. This will require integration of genetic epidemiological methods with functional genomics and proteomics. However, the use of these novel insights for genetic screening and personalised treatment lies some way off in the future.

Genetic and Functional Evidence Implicating DLL1 as the Gene That Influences Susceptibility to Visceral Leishmaniasis at Chromosome 6q27

Background. Visceral leishmaniasis (VL) is caused by Leishmania donovani and Leishmania infantum chagasi. Genome-wide linkage studies from Sudan and Brazil identified a putative susceptibility locus on chromosome 6q27. Methods. Twenty-two single-nucleotide polymorphisms (SNPs) at genes PHF10, C6orf70, DLL1, FAM120B, PSMB1, and TBP were genotyped in 193 VL cases from 85 Sudanese families, and 8 SNPs at genes PHF10, C6orf70, DLL1, PSMB1, and TBP were genotyped in 194 VL cases from 80 Brazilian families. Family-based association, haplotype, and linkage disequilibrium analyses were performed. Multispecies comparative sequence analysis was used to identify conserved noncoding sequences carrying putative regulatory elements. Quantitative reverse-transcription polymerase chain reaction measured expression of candidate genes in splenic aspirates from Indian patients with VL compared with that in the control spleen sample. Results. Positive associations were observed at PHF10, C6orf70, DLL1, PSMB1, and TBP in Sudan, but only at DLL1 in Brazil (combined P = 3 x 10(-4) at DLL1 across Sudan and Brazil). No functional coding region variants were observed in resequencing of 22 Sudanese VL cases. DLL1 expression was significantly (P = 2 x 10(-7)) reduced (mean fold change, 3.5 [SEM, 0.7]) in splenic aspirates from patients with VL, whereas other 6q27 genes showed higher levels (1.27 x 10(-6) < P < .01) than did the control spleen sample. A cluster of conserved noncoding sequences with putative regulatory variants was identified in the distal promoter of DLL1. Conclusions. DLL1, which encodes Delta-like 1, the ligand for Notch3, is strongly implicated as the chromosome 6q27 VL susceptibility gene.

Modelling dominance in a flexible intercross analysis

The aim of this paper is to develop a flexible model for analysis of quantitative trait loci (QTL) in outbred line crosses, which includes both additive and dominance effects. Our flexible intercross analysis (FIA) model accounts for QTL that are not fixed within founder lines and is based on the variance component framework. Genome scans with FIA are performed using a score statistic, which does not require variance component estimation. RESULTS: Simulations of a pedigree with 800 F2 individuals showed that the power of FIA including both additive and dominance effects was almost 50% for a QTL with equal allele frequencies in both lines with complete dominance and a moderate effect, whereas the power of a traditional regression model was equal to the chosen significance value of 5%. The power of FIA without dominance effects included in the model was close to those obtained for FIA with dominance for all simulated cases except for QTL with overdominant effects. A genome-wide linkage analysis of experimental data from an F2 intercross between Red Jungle Fowl and White Leghorn was performed with both additive and dominance effects included in FIA. The score values for chicken body weight at 200 days of age were similar to those obtained in FIA analysis without dominance. CONCLUSION: We have extended FIA to include QTL dominance effects. The power of FIA was superior, or similar, to standard regression methods for QTL effects with dominance. The difference in power for FIA with or without dominance is expected to be small as long as the QTL effects are not overdominant. We suggest that FIA with only additive effects should be the standard model to be used, especially since it is more computationally efficient.

Análise de ligação e associação no genoma com gagueira desenvolvimental persistente em famílias do Estado de São Paulo-Brasil

Pós-graduação em Ciências Biológicas (Genética) - IBB

QUANTITATIVE TRAIT LOCI FOR AGRONOMIC AND END-USE QUALITY PERFORMANCE AND THE EFFECT OF SOILBORNE WHEAT MOSAIC VIRUS IN A HARD WINTER WHEAT POPULATION IN NEBRASKA

To better understand agronomic and end-use quality in wheat (Triticum aestivum L.) we developed a population containing 154 F6:8 recombinant inbred lines (RILs) from the cross TAM107-R7/Arlin. The parental lines and RILs were phenotyped at six environments in Nebraska and differed for resistance to Wheat soilborne mosaic virus (WSBMV), morphological, agronomic, and end-use quality traits. Additionally, a 2300 cM genome-wide linkage map was created for quantitative trait loci (QTL) analysis. Based on our results across multiple environments, the best RILs could be used for cultivar improvement. The population and marker data are publicly available for interested researchers for future research. The population was used to determine the effect of WSBMV on agronomic and end-use quality and for the mapping of a resistance locus. Results from two infected environments showed that all but two agronomic traits were significantly affected by the disease. Specifically, the disease reduced grain yield by 30% of susceptible RILs and they flowered 5 d later and were 11 cm shorter. End-use quality traits were not negatively affected but flour protein content was increased in susceptible RILs. The resistance locus SbmTmr1 mapped to 27.1 cM near marker wPt-5870 on chromosome 5DL using ELISA data. Finally, we investigated how WSBMV affected QTL detection in the population. QTLs were mapped at two WSBMV infected environments, four uninfected environments, and in the resistant and susceptible RIL subpopulations in the infected environments. Fifty-two significant (LOD≥3) QTLs were mapped in RILs at uninfected environments. Many of the QTLs were pleiotropic or closely linked at 6 chromosomal regions. Forty-seven QTLs were mapped in RILs at WSBMV infected environments. Comparisons between uninfected and infected environments identified 20 common QTLs and 21 environmentally specific QTLs. Finally, 24 QTLs were determined to be affected by WSBMV by comparing the subpopulations in QTL analyses within the same environment. The comparisons were statistically validated using marker by disease interactions. These results showed that QTLs can be affected by WSBMV and careful interpretation of QTL results is needed where biotic stresses are present. Finally, beneficial QTLs not affected by WSBMV or the environment are candidates for marker-assisted selection.

Bipolar disorder susceptibility region on chromosome 3q29 not confirmed in a case-control association study

We identified a bipolar disorder (BPD) susceptibility region on chromosome 3q29 in a genome-wide linkage scan (Bailer et al. 2002 (Biol Psychiatry 52: 40), NPL-score 4.09) and follow-up linkage analysis (Schosser et al. 2004 (J Psychiatr Res 38(3): 357), NPL-scores >3 with five markers). These findings were supported by further fine-mapping of this region (Schosser et al. 2007 (Eur Neuropsychopharmacol 17(6-7): 501)), finding NPL-scores >3.9 with SNPs (single nucleotide polymorphisms) spanning a region of 3.46 Mbp in BPD families. Since genetic association studies are more powerful than linkage studies for detecting susceptibility genes of small effect size, we aimed to replicate these findings in an independent case-control sample collected in London (UK) and Vienna (Austria).

Study of a Swiss dopa-responsive dystonia family with a deletion in GCH1: redefining DYT14 as DYT5

OBJECTIVE: To report the study of a multigenerational Swiss family with dopa-responsive dystonia (DRD). METHODS: Clinical investigation was made of available family members, including historical and chart reviews. Subject examinations were video recorded. Genetic analysis included a genome-wide linkage study with microsatellite markers (STR), GTP cyclohydrolase I (GCH1) gene sequencing, and dosage analysis. RESULTS: We evaluated 32 individuals, of whom 6 were clinically diagnosed with DRD, with childhood-onset progressive foot dystonia, later generalizing, followed by parkinsonism in the two older patients. The response to levodopa was very good. Two additional patients had late onset dopa-responsive parkinsonism. Three other subjects had DRD symptoms on historical grounds. We found suggestive linkage to the previously reported DYT14 locus, which excluded GCH1. However, further study with more stringent criteria for disease status attribution showed linkage to a larger region, which included GCH1. No mutation was found in GCH1 by gene sequencing but dosage methods identified a novel heterozygous deletion of exons 3 to 6 of GCH1. The mutation was found in seven subjects. One of the patients with dystonia represented a phenocopy. CONCLUSIONS: This study rules out the previously reported DYT14 locus as a cause of disease, as a novel multiexonic deletion was identified in GCH1. This work highlights the necessity of an accurate clinical diagnosis in linkage studies as well as the need for appropriate allele frequencies, penetrance, and phenocopy estimates. Comprehensive sequencing and dosage analysis of known genes is recommended prior to genome-wide linkage analysis.

CHARACTERIZING A NOVEL GENETIC LOCUS ASSOCIATED WITH FAMILIAL CO-OCCURRENCE OF THORACIC AORTIC ANEURYSMS AND INTRACRANIAL ANEURYSMS

The Mendelian inheritance of genetic mutations can lead to adult-onset cardiovascular disease. Several genetic loci have been mapped for the familial form of Thoracic Aortic Aneurysms (TAA), and many causal mutations have been identified for this disease. Intracranial Aneurysms (ICA) also show linkage heterogeneity, but no mutations have been identified causing familial ICA alone. Here, we characterized a large family (TAA288) with an autosomal dominant pattern of inherited aneurysms. It is intriguing that female patients predominantly present with ICA and male patients predominantly with TAA in this family. To identify a causal mutation in this family, a genome-wide linkage analysis was previously performed on nine members of this family using the 50k GenChips Hind array from Affymetrix. This analysis eventually identified a single disease-segregating locus, on chromosome 5p15. We build upon this previous analysis in this study, hypothesizing that a genetic mutation inherited in this locus leads to the sex-specific phenotype of TAA and ICA in this family First we refined the boundaries of the 5p15 disease linked locus down to the genomic coordinates 5p15: 3,424,465- 6,312,925 (GRCh37/hg19 Assembly). This locus was named the TAA288 critical interval. Next, we sequenced candidate genes within the TAA288 critical interval. The selection of genes was simplified by the relatively small number of well-characterized genetic elements within the region. Seeking novel or rare disease-segregating variants, we initially observed a single point alteration in the metalloproteinase gene ADAMTS16 fulfilling this criteria. This variant was later classified as a low-frequency population polymorphism (rs72647757), but we continued to explore the potential role of the ADAMTS16 as the cause of disease in TAA288. We observed that fibroblasts cultured from TAA288 patients consistently upregulated the expression of this gene more strongly compared to matched control fibroblasts when treated with the cytokine TGF-β1, though there was some variation in the exact nature of this expression. We also observed evidence that this protein is expressed at elevated levels in aortic aneurysm tissue from patients with mutations in the gene TGFBR2 and Marfan syndrome, shown by immunohistochemical detection of this protein.

Nonlinear tests for genetic studies of complex disease

Linkage and association studies are major analytical tools to search for susceptibility genes for complex diseases. With the availability of large collection of single nucleotide polymorphisms (SNPs) and the rapid progresses for high throughput genotyping technologies, together with the ambitious goals of the International HapMap Project, genetic markers covering the whole genome will be available for genome-wide linkage and association studies. In order not to inflate the type I error rate in performing genome-wide linkage and association studies, multiple adjustment for the significant level for each independent linkage and/or association test is required, and this has led to the suggestion of genome-wide significant cut-off as low as 5 × 10 −7. Almost no linkage and/or association study can meet such a stringent threshold by the standard statistical methods. Developing new statistics with high power is urgently needed to tackle this problem. This dissertation proposes and explores a class of novel test statistics that can be used in both population-based and family-based genetic data by employing a completely new strategy, which uses nonlinear transformation of the sample means to construct test statistics for linkage and association studies. Extensive simulation studies are used to illustrate the properties of the nonlinear test statistics. Power calculations are performed using both analytical and empirical methods. Finally, real data sets are analyzed with the nonlinear test statistics. Results show that the nonlinear test statistics have correct type I error rates, and most of the studied nonlinear test statistics have higher power than the standard chi-square test. This dissertation introduces a new idea to design novel test statistics with high power and might open new ways to mapping susceptibility genes for complex diseases. ^

Gene by smoking interaction in blood pressure

Hypertension (HT) is mediated by the interaction of many genetic and environmental factors. Previous genome-wide linkage analysis studies have found many loci that show linkage to HT or blood pressure (BP) regulation, but the results were generally inconsistent. Gene by environment interaction is among the reasons that potentially explain these inconsistencies between studies. Here we investigate influences of gene by smoking (GxS) interaction on HT and BP in European American (EA), African American (AA) and Mexican American (MA) families from the GENOA study. A variance component-based method was utilized to perform genome-wide linkage analysis of systolic blood pressure (SBP), diastolic blood pressure (DBP), and HT status, as well as bivariate analysis for SBP and DBP for smokers, non-smokers, and combined groups. The most significant results were found for SBP in MA. The strongest signal was for chromosome 17q24 (LOD = 4.2), increased to (LOD = 4.7) in bivariate analysis but there was no evidence of GxS interaction at this locus (p = 0.48). Two signals were identified only in one group: on chromosome 15q26.2 (LOD = 3.37) in non-smokers and chromosome 7q21.11 (LOD = 1.4) in smokers, both of which had strong evidence for GxS interaction (p = 0.00039 and 0.009 respectively). There were also two other signals, one on chromosome 20q12 (LOD = 2.45) in smokers, which became much higher in the combined sample (LOD = 3.53), and one on chromosome 6p22.2 (LOD = 2.06) in non-smokers. Neither peak had very strong evidence for GxS interaction (p = 0.08 and 0.06 respectively). A fine mapping association study was performed using 200 SNPs in 30 genes located under the linkage signals on chromosomes 15 and 17. Under the chromosome 15 peak, the association analysis identified 6 SNPs accounting for a 7 mmHg increase in SBP in MA non-smokers. For the chromosome 17 linkage peak, the association analysis identified 3 SNPs accounting for a 6 mmHg increase in SBP in MA. However, none of these SNPs was significant after correcting for multiple testing, and accounting for them in the linkage analysis produced very small reductions in the linkage signal. ^ The linkage analysis of BP traits considering the smoking status produced very interesting signals for SBP in the MA population. The fine mapping association analysis gave some insight into the contribution of some SNPs to two of the identified signals, but since these SNPs did not remain significant after multiple testing correction and did not explain the linkage peaks, more work is needed to confirm these exploratory results and identify the culprit variations under these linkage peaks. ^

IDENTIFYING GENETIC VARIANTS AND CHARACTERIZING THEIR ROLE IN CLUBFOOT

Clubfoot is a common, complex birth defect affecting 4,000 newborns in the United States and 135,000 world-wide each year. The clubfoot deformity is characterized by inward and rigid downward displacement of one or both feet, along with persistent calf muscle hypoplasia. Despite strong evidence for a genetic liability, there is a limited understanding of the genetic and environmental factors contributing to the etiology of clubfoot. The studies described in this dissertation were performed to identify variants and/or genes associated with clubfoot. Genome-wide linkage scan performed on ten multiplex clubfoot families identified seven new chromosomal regions that provide new areas to search for clubfoot genes. Troponin C (TNNC2) the strongest candidate gene, located in 20q12-q13.11, is involved in muscle contraction. Exon sequencing of TNNC2 did not identify any novel coding variants. Interrogation of fifteen muscle contraction genes found strong associations with SNPs located in potential regulatory regions of TPM1 (rs4075583 and rs3805965), TPM2 (rs2025126 and rs2145925) and TNNC2 (rs383112 and rs437122). In previous studies, a strong association was found with rs3801776 located in the basal promoter of HOXA9, a gene also involved in muscle development and patterning. Altogether, this data suggests that SNPs located in potential regulatory regions of genes involved in muscle development and function could alter transcription factor binding leading to changes in gene expression. Functional analysis of 3801776/HOXA9, rs2025126/TPM2 and rs2145925/TPM2 showed altered protein binding, which significantly influenced promoter activity. Although the ancestral allele (G) of rs4075583/TPM1 creates a DNA-protein complex, it did not affect TPM1 promoter activity. However and importantly, in the context of a haplotype, rs4075583/G significantly decreased TPM1 promoter activity. These results suggest dysregulation of multiple skeletal muscle genes, TPM1, TPM2, TNNC2 and HOXA9, working in concert may contribute to clubfoot. However, specific allelic combinations involving these four regulatory SNPs did not confer a significantly higher risk for clubfoot. Other combinations of these variants are being evaluated. Moreover, these variants may interact with yet to be discovered variants in other genes to confer a higher clubfoot risk. Collectively, we show novel evidence for the role of skeletal muscle genes in clubfoot indicating that there are multiple genetic factors contributing to this complex birth defect.