New method for quantification and statistical analysis of nociceptive reflex receptive fields in humans

A method for quantifying nociceptive withdrawal reflex receptive fields in human volunteers and patients is described. The reflex receptive field (RRF) for a specific muscle denotes the cutaneous area from which a muscle contraction can be evoked by a nociceptive stimulus. The method is based on random stimulations presented in a blinded sequence to 10 stimulation sites. The sensitivity map is derived by interpolating the reflex responses evoked from the 10 sites. A set of features describing the size and location of the RRF is presented based on statistical analysis of the sensitivity map within every subject. The features include RRF area, volume, peak location and center of gravity. The method was applied to 30 healthy volunteers. Electrical stimuli were applied to the sole of the foot evoking reflexes in the ankle flexor tibialis anterior. The RRF area covered a fraction of 0.57+/-0.06 (S.E.M.) of the foot and was located on the medial, distal part of the sole of the foot. An intramuscular injection into flexor digitorum brevis of capsaicin was performed in one spinal cord injured subject to attempt modulation of the reflex receptive field. The RRF area, RRF volume and location of the peak reflex response appear to be the most sensitive measures for detecting modulation of spinal nociceptive processing. This new method has important potential applications for exploring aspects of central plasticity in volunteers and patients. It may be utilized as a new diagnostic tool for central hypersensitivity and quantification of therapeutic interventions.

Inhibition of FcepsilonRI-dependent mediator release and calcium flux from human mast cells by sialic acid-binding immunoglobulin-like lectin 8 engagement

BACKGROUND: Sialic acid-binding immunoglobulin-like lectins (Siglecs) are a family of glycan-binding inhibitory receptors, and among them, Siglec-8 is selectively expressed on human eosinophils, basophils, and mast cells. On eosinophils, Siglec-8 engagement induces apoptosis, but its function on mast cells is unknown. OBJECTIVE: We sought to study the effect of Siglec-8 engagement on human mast cell survival and mediator release responses. METHODS: Human mast cells were generated from CD34+ precursors. Apoptosis was studied by using flow cytometry. Mast cell mediator release or human lung airway smooth muscle contraction was initiated by FcepsilonRI cross-linking with or without preincubation with Siglec-8 or control antibodies, and release of mediators was analyzed along with Ca++ flux. RBL-2H3 cells transfected with normal and mutated forms of Siglec-8 were used to study how Siglec-8 engagement alters mediator release. RESULTS: Siglec-8 engagement failed to induce human mast cell apoptosis. However, preincubation with Siglec-8 mAbs significantly (P < .05) inhibited FcepsilonRI-dependent histamine and prostaglandin D(2) release, Ca++ flux, and anti-IgE-evoked contractions of human bronchial rings. In contrast, release of IL-8 was not inhibited. Siglec-8 ligation was also shown to inhibit beta-hexosaminidase release and Ca++ flux triggered through FcepsilonRI in RBL-2H3 cells transfected with full-length human Siglec-8 but not in cells transfected with Siglec-8 containing a tyrosine to phenylalanine point mutation in the membrane-proximal immunoreceptor tyrosine-based inhibitory motif domain. CONCLUSION: These data represent the first reported inhibitory effects of Siglec engagement on human mast cells.

Mechanical signals regulating extracellular matrix gene expression in fibroblasts

Mechanical forces are essential for connective tissue homeostasis. The extracellular matrix (ECM) plays a key role in the transmission of forces generated by the organism (e.g. muscle contraction) and externally applied (e.g. gravity). The expression of specific ECM proteins such as collagens and tenascin-C, as well as of matrix metalloproteinases, involved in their turnover, is influenced by mechanical stimuli. The precise mechanisms by which mechanical strains are translated into chemical signals and lead to differential gene expression are however not fully understood. Cell-matrix adhesion sites are good candidates for hosting a "mechanosensory switch", as they transmit forces from the ECM to the cytoskeleton and vice versa by physically linking the cytoskeleton to the ECM. Integrins, transmembrane proteins located to these adhesion sites, have been shown to trigger a set of internal signaling cascades after mechanical stimulation. We have shown that the expression level of tenascin-C directly correlates with externally applied mechanical stress, as well as with RhoA/RhoA-dependent kinase-mediated cytoskeletal tension. Presumably other genes are regulated in a similar manner. The changes in ECM composition and mechanical properties derived from mechanical stress are relevant in medical intervention after ligament and tendon injury.

Quantitative assessment of the nociceptive withdrawal reflex in healthy, non-medicated experimental sheep

This study aimed to characterize the nociceptive withdrawal reflex (NWR) and to define the nociceptive threshold in 25 healthy, non-medicated experimental sheep in standing posture. Electrical stimulation of the dorsal lateral digital nerves of the right thoracic and the pelvic limb was performed and surface-electromyography (EMG) from the deltoid (all animals) and the femoral biceps (18 animals) or the peroneus tertius muscles (7 animals) was recorded. The behavioural reaction following each stimulation was scored on a scale from 0 (no reaction) to 5 (strong whole body reaction). A train-of-five 1 ms constant-current pulse was used and current intensity was stepwise increased until NWR threshold intensity was reached. The NWR threshold intensity (It) was defined as the minimal stimulus intensity able to evoke a reflex with a minimal Root-Mean-Square amplitude (RMSA) of 20 μV, a minimal duration of 10 ms and a minimal reaction score of 1 (slight muscle contraction of the stimulated limb) within the time window of 20 to 130 ms post-stimulation. Based on this value, further stimulations were performed below (0.9It) and above threshold (1.5It and 2It). The stimulus-response curve was described. Data are reported as medians and interquartile ranges. At the deltoid muscle It was 4.4 mA (2.9–5.7) with an RMSA of 62 μV (30–102). At the biceps femoris muscle It was 7.0 mA (4.0–10.0) with an RMSA of 43 μV (34–50) and at the peroneus tertius muscle It was 3.4 mA (3.1–4.4) with an RMSA of 38 μV (32–46). Above threshold, RMSA was significantly increased at all muscles. Below threshold, RMSA was only significantly smaller than at It for the peroneus tertius muscle but not for the other muscles.

Hydrostatic intestinal edema induced signaling pathways: potential role of mechanical forces.

BACKGROUND: Hydrostatic intestinal edema initiates a signal transduction cascade that results in smooth muscle contractile dysfunction. Given the rapid and concurrent alterations in the mechanical properties of edematous intestine observed with the development of edema, we hypothesize that mechanical forces may serve as a stimulus for the activation of certain signaling cascades. We sought to examine whether isolated similar magnitude mechanical forces induced the same signal transduction cascades associated with edema. METHODS: The distal intestine from adult male Sprague Dawley rats was stretched longitudinally for 2 h to 123% its original length, which correlates with the interstitial stress found with edema. We compared wet-to-dry ratios, myeloperoxidase activity, nuclear signal transduction and activator of transcription (STAT)-3 and nuclear factor (NF)-kappa B DNA binding, STAT-3 phosphorylation, myosin light chain phosphorylation, baseline and maximally stimulated intestinal contractile strength, and inducible nitric oxide synthase (iNOS) and sodium hydrogen exchanger 1-3 messenger RNA (mRNA) in stretched and adjacent control segments of intestine. RESULTS: Mechanical stretch did not induce intestinal edema or an increase in myeloperoxidase activity. Nuclear STAT-3 DNA binding, STAT-3 phosphorylation, and nuclear NF-kappa B DNA binding were significantly increased in stretched seromuscular samples. Increased expression of sodium hydrogen exchanger 1 was found but not an increase in iNOS expression. Myosin light chain phosphorylation was significantly decreased in stretched intestine as was baseline and maximally stimulated intestinal contractile strength. CONCLUSION: Intestinal stretch, in the absence of edema/inflammatory/ischemic changes, leads to the activation of signaling pathways known to be altered in intestinal edema. Edema may initiate a mechanotransductive cascade that is responsible for the subsequent activation of various signaling cascades known to induce contractile dysfunction.

Denaturation and unfolding of human anaphylatoxin C3a: an unusually low covalent stability of its native disulfide bonds.

The complement C3a anaphylatoxin is a major molecular mediator of innate immunity. It is a potent activator of mast cells, basophils and eosinophils and causes smooth muscle contraction. Structurally, C3a is a relatively small protein (77 amino acids) comprising a N-terminal domain connected by 3 native disulfide bonds and a helical C-terminal segment. The structural stability of C3a has been investigated here using three different methods: Disulfide scrambling; Differential CD spectroscopy; and Reductive unfolding. Two uncommon features regarding the stability of C3a and the structure of denatured C3a have been observed in this study. (a) There is an unusual disconnection between the conformational stability of C3a and the covalent stability of its three native disulfide bonds that is not seen with other disulfide proteins. As measured by both methods of disulfide scrambling and differential CD spectroscopy, the native C3a exhibits a global conformational stability that is comparable to numerous proteins with similar size and disulfide content, all with mid-point denaturation of [GdmCl](1/2) at 3.4-5M. These proteins include hirudin, tick anticoagulant protein and leech carboxypeptidase inhibitor. However, the native disulfide bonds of C3a is 150-1000 fold less stable than those proteins as evaluated by the method of reductive unfolding. The 3 native disulfide bonds of C3a can be collectively and quantitatively reduced with as low as 1mM of dithiothreitol within 5 min. The fragility of the native disulfide bonds of C3a has not yet been observed with other native disulfide proteins. (b) Using the method of disulfide scrambling, denatured C3a was shown to consist of diverse isomers adopting varied extent of unfolding. Among them, the most extensively unfolded isomer of denatured C3a is found to assume beads-form disulfide pattern, comprising Cys(36)-Cys(49) and two disulfide bonds formed by two pair of consecutive cysteines, Cys(22)-Cys(23) and Cys(56)-Cys(57), a unique disulfide structure of polypeptide that has not been documented previously.

Polysensory interneuronal projections to foot contractile pedal neurons in Hermissenda.

A Pavlovian-conditioning procedure may produce modifications in multiple behavioral responses. As an example, conditioning may result in the elicitation of a specific somatomotor conditioned response (CR) and, in addition, other motor and visceral CRs. In the mollusk Hermissenda conditioning produces two conditioned responses: foot-shortening and decreased locomotion. The neural circuitry supporting ciliary locomotion is well characterized, although the neural circuit underlying foot-shortening is poorly understood. Here we describe efferent neurons in the pedal ganglion that produce contraction or extension of specific regions of the foot in semi-intact preparations. Synaptic connections between polysensory type Ib and type Is interneurons and identified foot contractile efferent neurons were examined. Type Ib and type Is interneurons receive synaptic input from the visual, graviceptive, and somatosensory systems. Depolarization of type Ib interneurons evoked spikes in identified tail and lateral foot contractile efferent neurons. Mechanical displacement of the statocyst evoked complex excitatory postsynaptic potentials (EPSPs) and spikes recorded from type Ib and type Is interneurons and complex EPSPs and spikes in identified foot contractile efferent neurons. Depolarization of type Ib interneurons in semi-intact preparations produced contraction and shortening along the rostrocaudal axis of the foot. Depolarization of Is interneurons in semi-intact preparations produced contraction of the anterior region of the foot. Taken collectively, the results suggest that type Ib and type Is polysensory interneurons may contribute to the neural circuit underlying the foot-shortening CR in Hermissenda.

Western diet changes cardiac acyl-CoA composition in obese rats: a potential role for hepatic lipogenesis.

The "lipotoxic footprint" of cardiac maladaptation in diet-induced obesity is poorly defined. We investigated how manipulation of dietary lipid and carbohydrate influenced potential lipotoxic species in the failing heart. In Wistar rats, contractile dysfunction develops at 48 weeks on a high-fat/high-carbohydrate "Western" diet, but not on low-fat/high-carbohydrate or high-fat diets. Cardiac content of the lipotoxic candidates--diacylglycerol, ceramide, lipid peroxide, and long-chain acyl-CoA species--was measured at different time points by high-performance liquid chromatography and biochemical assays, as was lipogenic capacity in the heart and liver by qRT-PCR and radiometric assays. Changes in membranes fluidity were also monitored using fluorescence polarization. We report that Western feeding induced a 40% decrease in myocardial palmitoleoyl-CoA content and a similar decrease in the unsaturated-to-saturated fatty acid ratio. These changes were associated with impaired cardiac mitochondrial membrane fluidity. At the same time, hepatic lipogenic capacity was increased in animals fed Western diet (+270% fatty acid elongase activity compared with high-fat diet), while fatty acid desaturase activity decreased over time. Our findings suggest that dysregulation of lipogenesis is a significant component of heart failure in diet-induced obesity.

Pharmacological characterization of nipecotic acid ethyl ester: A novel cholinergic muscarinic agonist

The aim of this dissertation was to examine the hypothesis that (R)-nipecotic acid ethyl ester ((R)-NAEE) is a cholinergic agonist that is selective for a particular subclass (M$\sb1$ or M$\sb2$) of muscarinic receptors.^ Ligand binding studies indicated that like cholinergic agonists (R)-NAEE selectively interacts with rat heart (M$\sb2$) and brain (M$\sb1$) muscarinic binding sites. Physiological studies revealed that unlike cholinergic agonists (R)-NAEE stimulated only those responses coupled to M$\sb2$ muscarinic receptors (acid secretion, negative inotropic response, smooth muscle contraction). Moreover, in rat brain (R)-NAEE differentiated between M$\sb2$ receptors negatively coupled to adenylate cyclase activity and M$\sb1$ receptors mediating PI turnover, being a weak competitive antagonist at these latter sites. In isolated rat gastric mucosal cells (R)-NAEE also differentiated between two M$\sb2$ coupled responses where it potentiated acid secretion but could not stimulate PI turnover. Atropine, a selective antimuscarinic agent, competitively antagonized all agonist effects of (R)-NAEE.^ Unlike (R)-NAEE, the muscarinic agonist arecoline, which is structurally similar to (R)-NAEE, stimulates both M$\sb1$ and M$\sb2$ receptors. Structure activity studies revealed that saturation of the piperidine ring and the length of the ester side chain of (R)-NAEE are the most important determinants for both M$\sb2$ efficacy and selectivity.^ The results of this dissertation establish that (R)-NAEE is a cholinergic muscarinic receptor agonist that displays greater efficacy at M$\sb2$ than at M$\sb1$ receptors, being a weak antagonist at the M$\sb1$ site. With such selectivity, (R)-NAEE may be regarded as a prototype for a unique class of cholinergic muscarinic M$\sb2$ receptor agonists. Because of these unique properties, (R)-NAEE should be useful in the further characterization of muscarinic receptors, and could lead to the development of a new class of therapeutic agents. ^

Studies to localize and characterize the isoform specific functional differences between cardiac and skeletal troponin C

Thin filament regulation of muscle contraction is a calcium dependent process mediated by the Tn complex. Calcium is released into the sarcomere and is bound by TnC. The subsequent conformation change in TnC is thought to begin a cascade of events that result in the activation of the actin-myosin ATPase. While the general events of this cascade are known, the molecular mechanisms of this signal transduction event are not. Recombinant DNA techniques, physiological and biochemical studies have been used to localize and characterize the structural domains of TnC that play a role in the calcium dependent signal transduction event that serves to trigger muscle contraction. The strategy exploited the observed functional differences between the isoforms of TnC to map regions of functional significance to the proteins. Chimeric cardiac-skeletal TnC proteins were generated to localize the domains of TnC that are required for maximal function in the myofibrilar ATPase assay. Characterization of these regions has yielded information concerning the molecular mechanism of muscle contraction. ^

Molecular mechanisms of Ca$\sp{2+}$ signal transduction from cardiac troponin C to troponin I

$\rm Ca\sp{2+}$-dependent exposure of an N-terminal hydrophobic region in troponin C (TnC) is thought to be important for the regulation of contraction in striated muscle. To study these conformational changes in cardiac troponin (cTnC), the $\varepsilon$C and $\varepsilon$H chemical shifts for all 10 Met residues in cTnC were sequence-specific assigned on NMR spectra using a combination of two dimensional NMR techniques and site-directed mutagenesis. The assigned methyl-Met chemical shifts were used as structural markers to monitor conformational changes induced by $\rm Ca\sp{2+}.$ The results showed that binding of $\rm Ca\sp{2+}$ to the regulatory site in the N-domain induced large changes in the $\varepsilon$H and $\varepsilon$C chemical shifts of Met 45, Met 80, Met 81 in the predicted N-terminal hydrophobic region, but had no effect on the chemical shifts of Met residues located in the C-domain. These results suggest that the $\rm Ca\sp{2+}$-dependent functions of cTnC are mainly through N-terminal domain of cTnC.^ To further define the molecular mechanism by which TnC regulates muscle contraction, single Cys residues were engineered at positions 45, 81, 84 or 85 in the N-terminal hydrophobic region of cTnC to provide sites for attachment of specific blocking groups. Blocking groups were coupled to these Cys residues in cTnC mutants and the covalent adducts were tested for activity in TnC-extracted myofibrils. Covalent modification of cTnC(C45) had no effect on maximal myofibril ATPase activity. Greatly decreased myofibril ATPase activity resulted when the peptide or biotin was conjugated to residue 81 in cTnC(C81), while less inhibition resulted from covalent modification of cTnC(C84) or cTnC(C85). The results suggest that limited sites of the N-terminal hydrophobic region in cTnC are important for transducing the $\rm Ca\sp{2+}$ signal to troponin I (TnI) and are sensitive to modification, while other regions are less important or can adapt to steric hindrances introduced by bulky blocking groups.^ Although the exposed TnI interaction site in the N-terminal hydrophobic region of TnC is crucial for function of TnC, other regions in the N-domain of TnC may also participate in transducing the $\rm Ca\sp{2+}$ signal and conferring the maximal activation of actomyosin ATPase. The interactions between the B-/C-helices of cTnC and cTnI were characterized using a combination of site-directed mutagenesis, fluorescence and covalent modification. The results suggest that the $\rm Ca\sp{2+}$-dependent interactions of the B-/C-helices of cTnC with TnI may be required for the maximal activation of muscle contraction. ^

A novel animal model for external anal sphincter insufficiency

Abstract PURPOSE: Reliable animal models are essential to evaluate future therapeutic options like cell-based therapies for external anal sphincter insufficiency. The goal of our study was to describe the most reliable model for external sphincter muscle insufficiency by comparing three different methods to create sphincter muscle damage. METHODS: In an experimental animal study, female Lewis rats (200-250 g) were randomly assigned to three treatment groups (n = 5, each group). The external sphincter muscle was weakened in the left dorsal quadrant by microsurgical excision, cryosurgery, or electrocoagulation by diathermy. Functional evaluation included in vivo measurements of resting pressure, spontaneous muscle contraction, and contraction in response to electrical stimulation of the afferent nerve at baseline and at 2, 4, and 6 weeks after sphincter injury. Masson's trichrome staining and immunofluorescence for skeletal muscle markers was performed for morphological analysis. RESULTS: Peak contraction after electrical stimulation was significantly decreased after sphincter injury in all groups. Contraction forces recovered partially after cryosurgery and electrocoagulation but not after microsurgical excision. Morphological analysis revealed an incomplete destruction of the external sphincter muscle in the cryosurgery and electrocoagulation groups compared to the microsurgery group. CONCLUSIONS: For the first time, three different models of external sphincter muscle insufficiency were directly compared. The animal model using microsurgical sphincter destruction offers the highest level of consistency regarding tissue damage and sphincter insufficiency, and therefore represents the most reliable model to evaluate future therapeutic options. In addition, this study represents a novel model to specifically test the external sphincter muscle function.

Evolutionary Genomics and Conservation of the Endangered Przewalski's Horse.

Przewalski's horses (PHs, Equus ferus ssp. przewalskii) were discovered in the Asian steppes in the 1870s and represent the last remaining true wild horses. PHs became extinct in the wild in the 1960s but survived in captivity, thanks to major conservation efforts. The current population is still endangered, with just 2,109 individuals, one-quarter of which are in Chinese and Mongolian reintroduction reserves [1]. These horses descend from a founding population of 12 wild-caught PHs and possibly up to four domesticated individuals [2-4]. With a stocky build, an erect mane, and stripped and short legs, they are phenotypically and behaviorally distinct from domesticated horses (DHs, Equus caballus). Here, we sequenced the complete genomes of 11 PHs, representing all founding lineages, and five historical specimens dated to 1878-1929 CE, including the Holotype. These were compared to the hitherto-most-extensive genome dataset characterized for horses, comprising 21 new genomes. We found that loci showing the most genetic differentiation with DHs were enriched in genes involved in metabolism, cardiac disorders, muscle contraction, reproduction, behavior, and signaling pathways. We also show that DH and PH populations split ∼45,000 years ago and have remained connected by gene-flow thereafter. Finally, we monitor the genomic impact of ∼110 years of captivity, revealing reduced heterozygosity, increased inbreeding, and variable introgression of domestic alleles, ranging from non-detectable to as much as 31.1%. This, together with the identification of ancestry informative markers and corrections to the International Studbook, establishes a framework for evaluating the persistence of genetic variation in future reintroduced populations.

Exploration of genetic susceptibility to pancreatic cancer: A GWAS approach

Pancreatic cancer is the 4th most common cause for cancer death in the United States, accompanied by less than 5% five-year survival rate based on current treatments, particularly because it is usually detected at a late stage. Identifying a high-risk population to launch an effective preventive strategy and intervention to control this highly lethal disease is desperately needed. The genetic etiology of pancreatic cancer has not been well profiled. We hypothesized that unidentified genetic variants by previous genome-wide association study (GWAS) for pancreatic cancer, due to stringent statistical threshold or missing interaction analysis, may be unveiled using alternative approaches. To achieve this aim, we explored genetic susceptibility to pancreatic cancer in terms of marginal associations of pathway and genes, as well as their interactions with risk factors. We conducted pathway- and gene-based analysis using GWAS data from 3141 pancreatic cancer patients and 3367 controls with European ancestry. Using the gene set ridge regression in association studies (GRASS) method, we analyzed 197 pathways from the Kyoto Encyclopedia of Genes and Genomes (KEGG) database. Using the logistic kernel machine (LKM) test, we analyzed 17906 genes defined by University of California Santa Cruz (UCSC) database. Using the likelihood ratio test (LRT) in a logistic regression model, we analyzed 177 pathways and 17906 genes for interactions with risk factors in 2028 pancreatic cancer patients and 2109 controls with European ancestry. After adjusting for multiple comparisons, six pathways were marginally associated with risk of pancreatic cancer ( P < 0.00025): Fc epsilon RI signaling, maturity onset diabetes of the young, neuroactive ligand-receptor interaction, long-term depression (Ps < 0.0002), and the olfactory transduction and vascular smooth muscle contraction pathways (P = 0.0002; Nine genes were marginally associated with pancreatic cancer risk (P < 2.62 × 10−5), including five reported genes (ABO, HNF1A, CLPTM1L, SHH and MYC), as well as four novel genes (OR13C4, OR 13C3, KCNA6 and HNF4 G); three pathways significantly interacted with risk factors on modifying the risk of pancreatic cancer (P < 2.82 × 10−4): chemokine signaling pathway with obesity ( P < 1.43 × 10−4), calcium signaling pathway (P < 2.27 × 10−4) and MAPK signaling pathway with diabetes (P < 2.77 × 10−4). However, none of the 17906 genes tested for interactions survived the multiple comparisons corrections. In summary, our current GWAS study unveiled unidentified genetic susceptibility to pancreatic cancer using alternative methods. These novel findings provide new perspectives on genetic susceptibility to and molecular mechanisms of pancreatic cancer, once confirmed, will shed promising light on the prevention and treatment of this disease. ^

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.