Versatile applications of transcriptional pulsing to study mRNA turnover in mammalian cells.

Development of transcriptional pulsing approaches using the c-fos and Tet-off promoter systems greatly facilitated studies of mRNA turnover in mammalian cells. However, optimal protocols for these approaches vary for different cell types and/or physiological conditions, limiting their widespread application. In this study, we have further optimized transcriptional pulsing systems for different cell lines and developed new protocols to facilitate investigation of various aspects of mRNA turnover. We apply the Tet-off transcriptional pulsing strategy to investigate ARE-mediated mRNA decay in human erythroleukemic K562 cells arrested at various phases of the cell cycle by pharmacological inhibitors. This application facilitates studies of the role of mRNA stability in control of cell-cycle dependent gene expression. To advance the investigation of factors involved in mRNA turnover and its regulation, we have also incorporated recently developed transfection and siRNA reagents into the transcriptional pulsing approach. Using these protocols, siRNA and DNA plasmids can be effectively cotransfected into mouse NIH3T3 cells to obtain high knockdown efficiency. Moreover, we have established a tTA-harboring stable line using human bronchial epithelial BEAS-2B cells and applied the transcriptional pulsing approach to monitor mRNA deadenylation and decay kinetics in this cell system. This broadens the application of the transcriptional pulsing system to investigate the regulation of mRNA turnover related to allergic inflammation. Critical factors that need to be considered when employing these approaches are characterized and discussed.

The molecular and genetic characterization of the MADS box transcription factor {\it Drosophila\/} MEF2

The myocyte enhancer factor (MEF)-2 family of transcription factors has been implicated in the regulation of muscle transcription in vertebrates, but the precise position of these regulators within the genetic hierarchy leading to myogenesis is unclear. The MEF2 proteins bind to a conserved A/T-rich DNA sequence present in numerous muscle-specific genes, and they are expressed in the cells of the developing somites and in the embryonic heart at the onset of muscle formation in mammals. The MEF2 genes belong to the MADS box family of transcription factors, which control specific programs of gene expression in species ranging from yeast to humans. Each MEF2 family member contains two highly conserved protein motifs, the MADS domain and the MEF2-specific domain, which together provide the MEF2 factors with their unique DNA binding and dimerization properties. In an effort to further define the function of the MEF2 proteins, and to evaluate the degree of conservation shared among these factors and the phylogenetic pathways that they regulate, we sought to identify MEF2 family members in other species. In Drosophila, a homolog of the vertebrate MEF2 genes was identified and termed D-mef2. The D-MEF2 protein binds to the consensus MEF2 element and can activate transcription through tandem copies of that site. During Drosophila embryogenesis, D-MEF2 is specific to the mesoderm germ layer of the developing embryo and becomes expressed in all muscle cell types within the embryo. The role of D-mef2 in Drosophila embryogenesis was examined by generating a loss-of-function mutation in the D-mef2 gene. In embryos homozygous for this mutant allele, somatic, cardiac, and visceral muscles fail to differentiate, but precursors of these myogenic lineages are normally specified and positioned. These results demonstrate that different muscle cell types share a common myogenic differentiation program controlled by MEF2 and suggest that this program has been conserved from Drosophila to mammals. ^

GENETIC ANALYSIS OF THE FUNCTION OF THE DROSOPHILA DOUBLESEX-RELATED FACTOR dmrt93B

DMRT (Doublesex and Mab-3 related transcription factor) proteins generally associated with sexual differentiation in many organisms share a common DNA binding domain and are often expressed in reproductive tissues. Aside from doublesex, which is a central factor in the regulation of sex determination, Drosophila possesses three different dmrt genes that are of unknown function. Because the association with sexual differentiation and reproduction is not universal and some DMRT proteins have been found to play other developmental roles we chose to further characterize one of these Drosophila genes. We carried out genetic analysis of dmrt93B, which was previously found to be expressed sex-specifically in the developing somatic gonad and to affect testis morphogenesis in RNAi knockdowns. In order to disrupt this gene, the GAL4 yeast transcriptional activator followed by a polyadenylation signal was inserted after the dmrt93B start codon and introduced into the genome by homologous recombination. Analysis of the knock-in mutation as well as a small deletion removing all dmrt93B sequence demonstrate that loss of function causes partial lethality at the late pupal stage. Surprisingly, these mutations have no significant effect on gonad formation or male fertility. Analysis of GAL4-driven GFP reporter expression indicates that the dmrt93B promoter activity is highly specific to neurons in the suboesophageal and proventricular ganglion in larva and adult of both sexes suggesting a possible role in digestive tract function. Using the Capillary Feeder (CAFÉ) assay to measure daily food intake we find that reduction in this gene’s function leads to an increase in food consumption. These results suggest dmrt93 plays an important role in the formation or maintenance of neurons that affect feeding and support the idea that dmrt genes may not be restricted to roles in sexual differentiation.

Genetic networks controlling retinal injury.

PURPOSE: The present study defines genomic loci underlying coordinate changes in gene expression following retinal injury. METHODS: A group of acute phase genes expressed in diverse nervous system tissues was defined by combining microarray results from injury studies from rat retina, brain, and spinal cord. Genomic loci regulating the brain expression of acute phase genes were identified using a panel of BXD recombinant inbred (RI) mouse strains. Candidate upstream regulators within a locus were defined using single nucleotide polymorphism databases and promoter motif databases. RESULTS: The acute phase response of rat retina, brain, and spinal cord was dominated by transcription factors. Three genomic loci control transcript expression of acute phase genes in brains of BXD RI mouse strains. One locus was identified on chromosome 12 and was highly correlated with the expression of classic acute phase genes. Within the locus we identified the inhibitor of DNA binding 2 (Id2) as a candidate upstream regulator. Id2 was upregulated as an acute phase transcript in injury models of rat retina, brain, and spinal cord. CONCLUSIONS: We defined a group of transcriptional changes associated with the retinal acute injury response. Using genetic linkage analysis of natural transcript variation, we identified regulatory loci and candidate regulators that control transcript levels of acute phase genes.

A Cbfa1-dependent genetic pathway controls bone formation beyond embryonic development.

The molecular mechanisms controlling bone extracellular matrix (ECM) deposition by differentiated osteoblasts in postnatal life, called hereafter bone formation, are unknown. This contrasts with the growing knowledge about the genetic control of osteoblast differentiation during embryonic development. Cbfa1, a transcriptional activator of osteoblast differentiation during embryonic development, is also expressed in differentiated osteoblasts postnatally. The perinatal lethality occurring in Cbfa1-deficient mice has prevented so far the study of its function after birth. To determine if Cbfa1 plays a role during bone formation we generated transgenic mice overexpressing Cbfa1 DNA-binding domain (DeltaCbfa1) in differentiated osteoblasts only postnatally. DeltaCbfa1 has a higher affinity for DNA than Cbfa1 itself, has no transcriptional activity on its own, and can act in a dominant-negative manner in DNA cotransfection assays. DeltaCbfa1-expressing mice have a normal skeleton at birth but develop an osteopenic phenotype thereafter. Dynamic histomorphometric studies show that this phenotype is caused by a major decrease in the bone formation rate in the face of a normal number of osteoblasts thus indicating that once osteoblasts are differentiated Cbfa1 regulates their function. Molecular analyses reveal that the expression of the genes expressed in osteoblasts and encoding bone ECM proteins is nearly abolished in transgenic mice, and ex vivo assays demonstrated that DeltaCbfa1-expressing osteoblasts were less active than wild-type osteoblasts. We also show that Cbfa1 regulates positively the activity of its own promoter, which has the highest affinity Cbfa1-binding sites characterized. This study demonstrates that beyond its differentiation function Cbfa1 is the first transcriptional activator of bone formation identified to date and illustrates that developmentally important genes control physiological processes postnatally.

Genetic and physical analysis of the bacteriophage P22 tail protein

A plasmid based genetic system was developed for the tail protein of the Salmonella typhimurium bacteriophage P22 and used to isolate and characterize tail protein mutants. The tail protein is a trimeric structural protein of the phage and an endorhamnosidase whose activity is essential for infection. The gene for the tail protein has previously been cloned into a plasmid expression vector and sequenced. A plate complementation assay for tail protein produced from the cloned gene was developed and used to isolate 27 tail protein mutants following mutagenesis of the cloned gene. These mutations were mapped into 12 deletion intervals using deletions which were made on plasmids in vitro and crossed onto P22. The base substitutions were determined by DNA sequencing. The majority of mutants had missense or nonsense mutations in the protein coding portion of the gene; however four of the mutants were in the putative transcription terminator. The oligomeric state of tail protein from the 15 missense mutants was investigated using SDS and nondenaturing polyacrylamide gel electrophoresis of cell lysates. Wild-type tail protein retains its trimeric structure in SDS gels at room temperature. Two of the mutant proteins also migrated as trimers in SDS gels, yet one of these had a considerably faster mobility than wild-type trimer. Its migration was the same as wild-type in a nondenaturing gel, so it is thought to be a trimer which is partially denatured by SDS. Four of the mutants produced proteins which migrate at the position of a monomer in an SDS gel but cannot be seen on a nondenaturing gel. These proteins are thought to be either monomers or soluble aggregates which cannot enter the nondenaturing gel. The remainder of mutants produce protein which is degraded. The mutant tail protein which had normal trimeric mobility on SDS and nondenaturing gels was purified. This protein has essentially wild-type ability to attach to phage capsids, but its endorhamnosidase activity is only 4% of wild-type. ^

STRUCTURAL ANALYSIS OF THE TRANSCRIPTION PRODUCTS OF A TEMPERATURE-SENSITIVE PHENOTYPIC MUTANT OF MOLONEY SARCOMA VIRUS--MUSV TS110

Cells infected with a temperature sensitive phenotypic mutant of Moloney sarcoma virus (MuSVts110) exhibit a transformed phenotype at 33('(DEGREES)) and synthesize two virus specific proteins, p85('gag-mos), a gag-mos fusion protein and p58('gag), a truncated gag precursor protein (the gag gene codes for viral structural proteins and mos is the MuSV transforming gene). At 39('(DEGREES)) only p58('gag) is synthesized and the morphology of the cells is similar to uninfected NRK parental cells. Two MuSVts110 specific RNAs are made in MuSVts110-infected cells, one of 4.0 kb in length, the other of 3.5 kb. Previous work indicated that each of these RNAs arose by a single central deletion of parental MuSV genetic material, and that p58('gag) was made by the 4.0 kb RNA and p85('gag-mos) from the 3.5 kb RNA. The objective of my dissertation research was to map precisely the deletion boundaries of both of the MuSVts110 RNAs, and to determine the proper reading frame across both deletion borders. This work succeeded in arriving at the following conclusions: (a) Using S-1 nuclease analysis and primer extension sequencing, it was found that the 4.0 kb MuSVts110 RNA arose by a 1488 base deletion of 5.2 kb parental MuSV genomic RNA. This deletion resulted in an out of frame fusion of the gag and mos genes that resulted in the formation of a "stop" codon which causes termination of translation just beyond the c-terminus of the gag region. Thus, this RNA can only be translated into the truncated gag protein p58('gag). (b) S-1 analysis of RNA from cells cultivated at different temperatures demonstrated that the 4.0 kb RNA was synthesized at all temperatures but that synthesis of the 3.5 kb RNA was temperature sensitive. These observations supported the data derived from blot hybridization experiments the interpretation of which argued for the existence of a single provirus in MuSVts110 infected cells, and hence only a single primary transcript (the 4.0 kb RNA). (c) Analyses similar to those described in (a) above showed that the 3.5 kb RNA was derived from the 4.0 kb MuSVts110 RNA by a further deletion of 431 bases, fusing the gag and mos genes into a continuous reading frame capable of directing synthesis of the p85('gag-mos) protein. These sequence data and the presence of only one MuSVts110-specific provirus, indicate that a splice mechanism is employed to generate the 3.5 kb RNA since the gag and mos genes are observed to be fused in frame in this RNA. . . . (Author's abstract exceeds stipulated maximum length. Discontinued here with permission of author.) UMI ^

GENETIC ORGANIZATION OF THE STRUCTURAL GENES FOR NITRATE REDUCTASE (TN5)

Nitrate reductase in Escherichia coli is a membrane-bound anaerobic enzyme that is repressed by oxygen and induced by nitrate. The genetic organization of the structural genes for the two larger subunits of nitrate reductase ((alpha) and (beta)) was determined by immunoprecipitation analysis of the formation of these proteins in nitrate reductase-deficient mutants resulting from transposon Tn5 mutagenesis. The results suggested that the genes encoding the (alpha) and (beta) subunits (narG and H) were arranged in an operon with transcription in the direction promoter(--->)(alpha)(--->)(beta). Segments of the chromosome containing the Tn5 inserts from several of the mutants were cloned into plasmid pBR322 and the positions of the transposons determined by restriction mapping. The Tn5 insertion sites were localized on two contiguous EcoRI fragments spanning about 6.6 kilobases of DNA. The narI gene (proposed to encode the (gamma) subunit) was positioned immediately downstream from the (beta)-gene (narH) by Southern analysis of Tn10 insertions into the narI locus. A Tn10 insertion into the narK locus, proposed to encode a nitrate-sensitive repressor of other anaerobic enzymes, was located about 1.5 kilobases upstream from the narGHI operon promoter. The narL locus, proposed to encode a nitrate-sensitive positive regulator of the narGHI operon and known to be genetically linked to the other nar genes, was demonstrated to lie outside a 19.3-kilobase region of the chromosome which encompasses the other nar genes. The physical limit of the narGHI promoter was defined by studying the effect of Tn5 insertions into a hybrid plasmid containing the functional operon. The points of origin of the coding regions for the (alpha) and (beta) genes were deduced by alignment of the chromosomal map of Tn5 insertion sites with the sizes of (alpha) and (beta) subunit fragments produced by plasmids carrying these Tn5 inserts in the nar operon. The coding region for the (alpha) subunit (143,000 daltons) begins about 250 nucleotides downstream from the deduced limit of the promoter region and includes about 4.0 kilobases of DNA; the region encoding (beta) (60,000 daltons) lies immediately downstream from the (alpha)-gene and is approximately 1.6 kilobases in length. The adjacent region encoding the (gamma) subunit (19,000 daltons) is approximately 0.5 kilobase in length. ^

In vivo definition of genetic and cellular aspects of mammalian sexual differentiation

The differentiation of the reproductive organs is an essential developmental process required for the proper transmission of the genetic material. Müllerian inhibiting substance (MIS) is produced by testes and is necessary for the regression of the Müllerian ducts: the anlagen of the uterus, fallopian tubes and cervix. In vitro and standard transgenic mouse studies indicate that the nuclear hormone receptor Steroidogenic factor 1 (SF-1) and the transcription factor SOX9 play an essential role in the regulation of Mis. To test this hypothesis, mutations in the endogenous SF-1 and SOX9 binding sites in the mouse Mis promoter were introduced by gene targeting in embryonic stem (ES) cells. In disagreement with cell culture and transgenic mouse studies, male mice homozygous for the mutant SF-1 binding site correctly initiated Mis transcription in the fetal testes, although at significantly reduced levels. Surprisingly, sufficient Mis was produced for complete elimination of the Müllerian duct system. However, when the SF-1 binding site mutation was combined with an Mis -null allele, the further decrease in Mis levels led to a partial retention of uterine tissue, but only at a distance from the testes. In contrast, males homozygous for the mutant SOX9 binding site did not initiate Mis transcription, resulting in pseudohermaphrodites with a uterus and oviducts. These studies suggest an essential role for SOX9 in the initiation of Mis transcription, whereas SF-1 appears to act as a quantitative regulator of Mis transcript levels perhaps for influencing non-Müllerian duct tissues. ^ The Mis type II receptor, a member of the TGF- b superfamily, is also required for the proper regression of the Müllerian ducts. Mis type II receptor-deficient human males and their murine counterparts develop as pseudohermaphrodites. A lacZ reporter cassette was introduced into the mouse Mis type II receptor gene, by homologous recombination in ES cells. Expression studies, based on b -galactosidase activity, show marked expression of the MIS type II receptor in the postnatal Sertoli cells of the testis as well as in the prenatal and postnatal granulosa cells of the ovary. Expression is also seen in the mesenchymal cells surrounding the Müllerian duct and in the longitudinal muscle layer of the uterus. ^

Biochemical and genetic analysis of two protein kinases JNKK2 and MEKK3

Mitogen-activated protein kinase (MAPK) cascades are conserved eukaryotic signaling modules consisting of a MAPK, a MAPKK and a MAP3K. MAPK cascades are involved in many cellular responses including proliferation, differentiation, apoptosis, stress and immune responses. ^ The first part of this thesis describes the cloning and biochemical analysis of JNKK2, a member of MAPKK gene family. Our results demonstrate that JNKK2 is a specific JNK activator and activates the JNK-dependent signal transduction pathway in vivo by inducing c-Jun and ATF2-mediated gene expression. We also found that JNKK2 is specifically activated by a MAP3K MEKK2 through formation of MEKK2-JNKK2-JNK1 triple complex module. JNKK2 is likely to mediate specific upstream signals to activate JNK cascade. ^ The second part of this thesis describes biochemical and gene disruption analysis of MEKK3, a member of MAP3K gene family. We showed that overexpression of MEKK3 strongly activates both JNK and p38 MAPKs but only weakly activates ERK. MEKK−/− embryos die at about embryonic day (E) 11. MEKK3−/− embryos displayed defects in blood vessel development in the yolk sacs, and in the myocardium and endocardium development at E9.5. The angiogenesis in the head, intersomitic region and placenta was also abnormal. These results demonstrate that MEKK3, a member of MAP3K MEKK/STE11 subgene family, is essential for early embryonic cardiovascular development. Furthermore, it was found that disruption of MEKK3 did not alter the expression of vascular endothelial growth factor-1 (VEGF-1), angiopoietin-1, -2 and their respective receptors Flt-1, Flk-1, Tie-1, Tie-2. Finally, MEKK3 was shown to activate myocyte-specific enhancer factor 2C (MEF2C), a crucial transcription factor for early embryonic cardiovascular development through the p38 MAPK cascade, suggesting that MEF2C is one of the key targets of the MEEKK3 signaling pathway during early embryonic cardiovascular development. ^

Genetic mechanisms of Na+, K+-ATPase regulation in hypertension

The spontaneously hypertensive rat (SHR) is a model of essential hypertension. During the early development of hypertension, the SHR demonstrates increased proximal tubule (PT) Na+ reabsorption. I hypothesized that the increased PT Na+ reabsorption exhibited by the young SHR was due to altered sub-cellular distribution of Na+, K +-ATPase compared to the normotensive Wistar Kyoto (WKY). The hypothesis is supported, herein, by observations of greater Na+, K +-ATPase α 1 abundance in PT plasma membrane and lower abundance in late endosomes of 4wk SHR despite no difference in total PT α 1 abundance. There is a greater amount of Ser-18 unphosphorylated α 1 in the 4wk SHR PT. Total PT Na+, K+-ATPase γ abundance is greater in SHR at 4wk and 16wk but γ abundance in plasma membrane is greater only at 4wk. The phosphatase, calcineurin, was chosen for study because it is involved in the stimulation of Na+, K +-ATPase. No difference in calcineurin coding sequence, expression, or activity was observed in SHR. Gene expression arrays were next used to find candidate genes involved in the regulation of Na+, K +-ATPase. The first candidate analyzed was soluble epoxide hydrolase (sEH). The gene encoding sEH (EPHX2) showed lower expression in SHR. There was also a reduction in sEH protein abundance but there was no correlation between protein abundance and blood pressure in F2 progeny. Two EPHX2 alleles were identified, an ancestral allele and a variant allele containing four polymorphisms. sEH activity was greater in animals carrying the variant allele but the inheritance of the variant allele did not correlate with blood pressure. Gene expression arrays also led to the examination of genes involved in redox balance/Na+, K+-ATPase regulation. A pattern of lower expression of genes involved in reactive radical detoxification in SHR was discerned. Six transcription factor binding sites were identified that occurred more often in these genes. Three transcription factors that bind to the HNF1 site were expressed at lower levels in SHR. This points to the HNF1 transcriptional complex as an important trans-acting regulator of a wide range of genes involved in altered redox balance in SHR. ^

Molecular and genetic analyses of Fgf signaling during cardiac outflow tract development

Epithelial-mesenchymal tissue interactions regulate the development of derivatives of the caudal pharyngeal arches (PAs) to govern the ultimate morphogenesis of the aortic arch and outflow tract (OFT) of the heart. Disruption of these signaling pathways is thought to contribute to the pathology of a significant proportion of congenital cardiovascular defects in humans. In this study, I tested whether Fibroblast Growth Factor 15 (Fgf15), a secreted signaling molecule expressed within the PAs, is an extracellular mediator of tissue interactions during PA and OFT development. Analyses of Fgf15−/− mouse embryonic hearts revealed abnormalities primarily localized to the OFT, correlating with aberrant cardiac neural crest cell behavior. The T-box-containing transcription factor Tbx1 has been implicated in the cardiovascular defects associated with the human 22q11 Deletion Syndromes, and regulates the expression of other Fgf family members within the mouse PAs. However, expression and genetic interaction studies incorporating mice deficient for Tbx1, its upstream regulator, Sonic Hedgehog (Shh), or its putative downstream effector, Fgf8, indicated that Fgf15 functions during OFT development in a manner independent of these factors. Rather, analyses of compound mutant mice indicated that Fgf15 and Fgf9, an additional Fgf family member expressed within the PAs, genetically interact, providing insight into the factors acting in conjunction with Fgf15 during OFT development. Finally, in an effort to further characterize this Fgf15-mediated developmental pathway, promoter deletion analyses were employed to isolate a 415bp sequence 7.1Kb 5′ to the Fgf15 transcription start site both necessary and sufficient to drive reporter gene expression within the epithelium of the PAs. Sequence comparisons among multiple mammalian species facilitated the identification of evolutionarily conserved potential trans-acting factor binding sites within this fragment. Subsequent studies will investigate the molecular pathway(s) through which Fgf15 functions via identification of factors that bind to this element to govern Fgf15 gene expression. Furthermore, targeted deletion of this element will establish the developmental requirement for pharyngeal epithelium-derived Fgf15 signaling function. Taken as a whole, these data demonstrate that Fgf15 is a component of a novel, Tbx1-independent molecular pathway, functioning within the PAs in a manner cooperative with Fgf9, required for proper development of the cardiac OFT. ^

The eukaryotic cellular stress response: Biochemical and genetic analyses in Saccharomyces cerevisiae

All cells must have the ability to deal with a variety of environmental stresses. Failure to correctly adapt to and/or protect against adverse stress conditions can lead to cell death. In humans, stress response defects have been linked to a number of neurodegenerative diseases and cancer, underscoring the importance of developing a fundamental understanding of the eukaryotic stress response.^ In an effort to characterize cellular response to high temperature stress, I identified and described one member of a novel gene family— RTR1. I show that the RTR1 gene and its protein product genetically and biochemically interact with core subunits of the RNA polymerase II enzyme. Appropriately, loss of RTR1 results in defective transcription from multiple promoters. These data provide evidence that Rtr1, which is essential under stress conditions, acts as a key regulator of transcription.^ In addition to transcriptional regulation, cells deal with many stressors by inducing molecular chaperones. Molecular chaperones are ubiquitous in all living cells and bind unfolded or damaged proteins and catalyze refolding or degradation. Hsp90 is a unique chaperone because it targets specific clients—typically signaling proteins—for maturation. While it has been shown that Sse1, the yeast Hsp110, is a critical regulator of the Hsp90 chaperone cycle, this work describes the molecular basis for that regulation. I show that Sse1 modulates Hsp90 function through regulation of Hsp70 nucleotide exchange. Further, Hsp110-type nucleotide exchange factors (NEFs) appear to have a specific role in modulating Hsp90 function in this manner. Finally, in addition to Hsp110, the eukaryotic cytosol contains two other types of Hsp70 NEF: Snl1 (BAG-domain protein) and Fes1 (HspBP1-like protein). I investigated the cellular roles of these NEFs to better understand the reason that eukaryotic cells contain three distinct protein families that perform the same biochemical function. I show that while cytsolic Hsp70 NEFs have some degree of functional overlap, they also exhibit striking divergence. Taken together, the work presented in this dissertation provides a more detailed understanding of the eukaryotic stress response. ^

Molecular epidemiological studies on genetic predispositions to squamous cell carcinoma of the head and neck

Two molecular epidemiological studies were conducted to examine associations between genetic variation and risk of squamous cell carcinoma of the head and neck (SCCHN). In the first study, we hypothesized that genetic variation in p53 response elements (REs) may play roles in the etiology of SCCHN. We selected and genotyped five polymorphic p53 REs as well as a most frequently studied p53 codon 72 (Arg72Pro, rs1042522) polymorphism in 1,100 non-Hispanic White SCCHN patients and 1,122 age-and sex-matched cancer-free controls recruited at The University of Texas M. D. Anderson Cancer Center. In multivariate logistic regression analysis with adjustment for age, sex, smoking and drinking status, marital status and education level, we observed that the EOMES rs3806624 CC genotype had a significant effect of protection against SCCHN risk (adjusted odds ratio= 0.79, 95% confidence interval =0.64–0.98), compared with the -838TT+CT genotypes. Moreover, a significantly increased risk associated with the combined genotypes of p53 codon 72CC and EOMES -838TT+CT was observed, especially in the subgroup of non-oropharyneal cancer patients. The values of false-positive report probability were also calculated for significant findings. In the second study, we assessed the association between SCCHN risk and four potential regulatory single nucleotide polymorphisms (SNPs) of DEC1 (deleted in esophageal cancer 1) gene, a candidate tumor suppressor gene for esophageal cancer. After adjustment for age, sex, and smoking and drinking status, the variant -606CC (i.e., -249CC) homozygotes had a significantly reduced SCCHN risk (adjusted odds ratio = 0.71, 95% confidence interval = 0.52–0.99), compared with the -606TT homozygotes. Stratification analyses showed that a reduced risk associated with the -606CC genotype was more pronounced in subgroups of non-smokers, non-drinkers, younger subjects (defined as ≤ 57 years), carriers of TP53 Arg/Arg (rs1042522) genotype, patients with oropharyngeal cancer or late-stage SCCHN. Further in silico analysis revealed that the -249 T-to-C change led to a gain of a transcription factor binding site. Additional functional analysis showed that the -249T-to-C change significantly enhanced transcriptional activity of the DEC1 promoter and the DNA-protein binding activity. We conclude that the DEC1 promoter -249 T>C (rs2012775) polymorphism is functional, modulating susceptibility to SCCHN among non-Hispanic Whites. Additional large-scale, preferably population-based studies are needed to validate our findings.^

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.