Regulation of XMyoDa transcription by the binding of XMEF2A to the TATA box

MEF2 is a $\underline{\rm m}$yocyte-specific $\underline{\rm e}$nhancer-binding $\underline{\rm f}$actor that binds a conserved DNA sequence, CTA(A/T)$\sb4$TAG. A MEF2 binding site in the XMyoDa promoter overlaps with the TATA box and is required for muscle specific expression. To examine the potential role of MEF2 in the regulation of MyoD transcription during early development, the appearance of MEF2 binding activity in developing Xenopus embryos was analyzed with the electrophoretic mobility shift assay. Two genes were isolated from a X. Laevis stage 24 cDNA library that encode factors that bind the XMyoDa TFIID/MEF2 site. Both genes are highly homologous to each other, belong to the MADS ($\underline{\rm M}$CM1-$\underline{\rm A}$rg80-agamous-$\underline{\rm d}$eficiens-$\underline{\rm S}$RF) protein family, and most highly related to the mammalian MEF2A gene, hence they are designated as XMEF2A1 and XMEF2A2. Proteins encoded by both cDNAs form specific complexes with the MEF2 binding site and show the same binding specificity as the endogenous MEF2 binding activity. XMEF2A transcripts accumulate preferentially in developing somites after the appearance of XMyoD transcripts. XMEF2 protein begins to accumulate in somites at tailbud stages. Transcriptional activation of XMyoD promoter by XMEF2A required only the MADS box and MEF2-specific domain when XMEF2A is bound at the TATA box. However, a different downstream transactivation domain was required when XMEF2A activates transcription through binding to multiple upstream sites. These results suggest that different activation mechanisms are involved, depending on where the factor is bound. Mutations in several basic amino acid clusters in the MADS box inhibit DNA binding suggesting these amino acids are essential for DNA binding. Mutation of Thr-20 and Ser-36 to the negatively charged amino acid residue, aspartic acid, abolish DNA binding. XMEF2A activity may be regulated by phosphorylation of these amino acids. A dominant negative mutant was made by mutating one of the basic amino acid clusters and deleting the downstream transactivation domain. In vivo roles of MEF2 in the regulation of MyoD transcription were investigated by overexpression of wild type MEF2 and dominant negative mutant of XMEF2A in animal caps and assaying for the effects on the level of expression of MyoD genes. Overexpression of MEF2 activates the transcription of endogenous MyoD gene family while expression of a dominant negative mutant reduces the level of transcription of XMRF4 and myogenin genes. These results suggest that MEF2 is downstream of MyoD and Myf5 and that MEF2 is involved in maintaining and amplifying expression of MyoD and Myf5. MEF2 is upstream of MRF4 and myogenin and plays a role in activating their expression. ^

Functional characterization of transcription factor USF

USF, Upstream Stimulatory Factor, is a family of ubiquitous transcription factors that contain highly conserved basic helix-loop-helix leucine zipper DNA binding domains and recognize the core DNA sequence CACGTG. In human and mouse, two members of the USF family, USF1 and USF2, encoded by two different genes, contribute to the USF activity. In order to gain insights into the mechanisms by which USFs function as transcriptional activators, different approaches were used to map the domains of USF2 responsible for nuclear localization and transcriptional activation. Two stretches of amino acids, one in the basic region of the DNA binding domain, the other in a highly conserved N-terminal region, were found to direct nuclear localization independently of one another. Two distinct activation domains were also identified. The first one, located in the conserved N-terminal region that overlaps the C-terminal nuclear localization signal, functioned only in the presence of an initiator element in the promoter of the reporter. The second, in a nonconserved region, activated transcription in the absence of an initiator element or when fused to a heterologous DNA binding domain. These results suggest that USF2 functions in different promoter contexts by selectively utilizing different activation domains.^ The deletion analysis of USF2 also identified two dominant negative mutants of USF, one lacking the activation domain, the other lacking the basic domain. The latter proved useful for testing the direct involvement of USFs in the transcriptional activation mediated by the viral protein IE62.^ To investigate the biological function of USFs, foci and colony formation assays were used to study the growth regulation by USFs. It was found that USFs had a strong antagonistic effect on cellular transformation mediated by the bHLH/LZ protein Myc. This effect required the DNA binding activity of either USF 1 or USF2. Moreover, USF2, but not USF1 or other mutants of USFs, was also found to have strong inhibitory effect on the cellular transformation by E1a and on the growth of HeLa cells. These results demonstrate that USFs could potentially regulate growth through two mechanisms, one by antagonizing the function of Myc in cellular transformation, the other by mediating a more general growth inhibitory effect. ^

Transcriptional regulation of the human {\it PAX6\/} gene

PAX6, a member of the paired-type homeobox gene family, is expressed in a partially and temporally restricted pattern in the developing central nervous system, and its mutation is responsible for human aniridia (AN) and mouse small eye (Sey). The objective of this study was to characterize the PAX6 gene regulation at the transcriptional level, and thereby gain a better understanding of the molecular basis of the dynamic expression pattern and the diversified function of the human PAX6 gene.^ Initially, we examined the transcriptional regulation of the PAX6 gene by transient transfection assays and identified multiple cis-regulatory elements that function differently in different cell lines. The transcriptional initiation site was identified by RNase protection and primer extension assays. Examination of the genomic DNA sequence indicated that the PAX6 promoter has a TATA like-box (ATATTTT) at $-$26 bp, and two CCAAT-boxes are located at positions $-$70 and $-$100 bp. A 38 bp ply (CA) sequence was located 992 bp upstream from the initiation site. Transient transfection assays in glioblastoma cells and leukemia cells indicate that a 92 bp region was required for basal level PAX6 promoter activity. Gel retardation assays showed that this 92 bp sequence can form four DNA-protein complexes which can be specifically competed by a 31-mer oligonucleotide containing a PAX6 TATA-like sequence or an adenovirus TATA box. The activation of the promoter is positively correlated with the expression of PAX6 transcripts in cells tested.^ Based on the results obtained from the in vitro transfection assays, we did further dissection assay and functional analysis in both cell-culture and transgenic mice. We found that a 5 kb upstream promoter sequence is required for the tissue specific expression in the forebrain region which is consistent with that of the endogenous PAX6 gene. A 267 bp cell-type specific repressor located within the 5 kb fragment was identified and shown to direct forebrain specific expression. The cell-type specific repressor element has been narrowed to a 30 bp region which contains a consensus E-box by in vitro transfection assays. The third regulatory element identified was contained in a 162 bp sequence (+167 to +328) which functions as a midbrain repressor, and it appeared to be required for establishing the normal expression pattern of the PAX6 gene. Finally, a highly conserved 216 bp sequence identified in intron 4 exhibited as a spinal cord specific enhancer. And this 216 bp cis-regulatory element can be used as a marker to trace the differentiation and migration of progenitor cells in the developing spinal cord. These studies show that the concerted action of multiple cis-acting regulatory elements located upstream and downstream of the transcription initiation site determines the tissue specific expression of PAX6 gene. ^

myogenin: Human candidate gene and knockout mouse effect

The myogenin gene encodes an evolutionarily conserved basic helix-loop-helix transcription factor that regulates the expression of skeletal muscle-specific genes and its homozygous deletion results in mice who die of respiratory failure at birth. The histology of skeletal muscle in the myogenin null mice is reminiscent of that found in some severe congenital myopathy patients, many of whom also die of respiratory complications and provides the rationale that an aberrant human myogenin (myf4) coding region could be associated with some congenital myopathy conditions.^ With PCR, we found similarly sized amplimers for the three exons of the myogenin gene in 37 patient and 40 control samples. In contrast to the GeneBank sequence for human myogenin, we report several differences in flanking and coding regions plus an additional 659 and 498 bps in the first and second introns, respectively, in all patients and controls. We also find a novel (CA)-dinucleotide repeat in the second intron. No causative mutations were detected in the myogenin coding regions of genomic DNA from patients with severe congenital myopathy.^ Severe congenital myopathies in humans are often associated with respiratory complications and pulmonary hypoplasia. We have employed the myogenin null mouse, which lacks normal development of skeletal muscle fibers as a genetically defined severe congenital myopathy mouse model to evaluate the effect of absent fetal breathing movement on pulmonary development.^ Significant differences are observed at embryonic days E14, E17 and E20 of lung:body weight, total DNA and histologically, suggesting that the myogenin null lungs are hypoplastic. RT-PCR, in-situ immunofluorescence and EM reveal pneumocyte type II differentiation in both null and wild lungs as early as E14. However, at E14, myogenin null lungs have decreased BrdU incorporation while E17 through term, augmented cell death is detected in the myogenin null lungs, not seen in wild littermates. Absent mechanical forces appear to impair normal growth, but not maturation, of the developing lungs in myogenin null mouse.^ These investigations provide the basis for delineating the DNA sequence of the myogenin gene and and highlight the importance of skeletal muscle development in utero for normal lung organogenesis. My observation of no mutations within the coding regions of the human myogenin gene in DNA from patients with severe congenital myopathy do not support any association with this condition. ^

Two prognostic factors in breast cancer: Novel functions of EGFR and beta-catenin

In this study, we demonstrated the novel functions of two important prognostic markers in breast cancer, EGFR and b -catenin in proliferation and/or other transformation phenotype. ^ First we demonstrated that EGFR could be detected in the nucleus in highly proliferating tissues, including primary breast cancer samples and a breast cancer cell line. We found that EGFR contained a strong transactivation domain, complexed with an AT-rich consensus DNA sequence and activated promoters containing this sequence, including cyclin D1 promoter. Therefore, EGFR may function as a transcription factor to activate genes required for highly proliferating activity such as cyclin D1 in breast cancer. ^ In the second part of this study, we identified b -catenin as an important prognostic factor in breast cancer. We found that cyclin D1 was one of the genes regulated by b -catenin in breast cancer cells. The transactivation activity of b -catenin correlated significantly with cyclin D1 expression in both breast cancer cell lines and in breast cancer patient samples, in which high b -catenin activity correlated with poor prognosis of the patients. Moreover, blockage of b -catenin activity significantly inhibited transformation phenotypes in breast cancer cells. Therefore, our results indicate that b -catenin can be involved in breast cancer formation and/or progression and may serve as a target for breast cancer therapy. ^

New algorithms for automated phylogenetic reconstruction using artificial intelligence and data mining techniques

Academic and industrial research in the late 90s have brought about an exponential explosion of DNA sequence data. Automated expert systems are being created to help biologists to extract patterns, trends and links from this ever-deepening ocean of information. Two such systems aimed on retrieving and subsequently utilizing phylogenetically relevant information have been developed in this dissertation, the major objective of which was to automate the often difficult and confusing phylogenetic reconstruction process. ^ Popular phylogenetic reconstruction methods, such as distance-based methods, attempt to find an optimal tree topology (that reflects the relationships among related sequences and their evolutionary history) by searching through the topology space. Various compromises between the fast (but incomplete) and exhaustive (but computationally prohibitive) search heuristics have been suggested. An intelligent compromise algorithm that relies on a flexible “beam” search principle from the Artificial Intelligence domain and uses the pre-computed local topology reliability information to adjust the beam search space continuously is described in the second chapter of this dissertation. ^ However, sometimes even a (virtually) complete distance-based method is inferior to the significantly more elaborate (and computationally expensive) maximum likelihood (ML) method. In fact, depending on the nature of the sequence data in question either method might prove to be superior. Therefore, it is difficult (even for an expert) to tell a priori which phylogenetic reconstruction method—distance-based, ML or maybe maximum parsimony (MP)—should be chosen for any particular data set. ^ A number of factors, often hidden, influence the performance of a method. For example, it is generally understood that for a phylogenetically “difficult” data set more sophisticated methods (e.g., ML) tend to be more effective and thus should be chosen. However, it is the interplay of many factors that one needs to consider in order to avoid choosing an inferior method (potentially a costly mistake, both in terms of computational expenses and in terms of reconstruction accuracy.) ^ Chapter III of this dissertation details a phylogenetic reconstruction expert system that selects a superior proper method automatically. It uses a classifier (a Decision Tree-inducing algorithm) to map a new data set to the proper phylogenetic reconstruction method. ^

Candidate gene identification following linkage: Search for hypertension susceptibility genes

Following up genetic linkage studies to identify the underlying susceptibility gene(s) for complex disease traits is an arduous yet biologically and clinically important task. Complex traits, such as hypertension, are considered polygenic with many genes influencing risk, each with small effects. Chromosome 2 has been consistently identified as a genomic region with genetic linkage evidence suggesting that one or more loci contribute to blood pressure levels and hypertension status. Using combined positional candidate gene methods, the Family Blood Pressure Program has concentrated efforts in investigating this region of chromosome 2 in an effort to identify underlying candidate hypertension susceptibility gene(s). Initial informatics efforts identified the boundaries of the region and the known genes within it. A total of 82 polymorphic sites in eight positional candidate genes were genotyped in a large hypothesis-generating sample consisting of 1640 African Americans, 1339 whites, and 1616 Mexican Americans. To adjust for multiple comparisons, resampling-based false discovery adjustment was applied, extending traditional resampling methods to sibship samples. Following this adjustment for multiple comparisons, SLC4A5, a sodium bicarbonate transporter, was identified as a primary candidate gene for hypertension. Polymorphisms in SLC4A5 were subsequently genotyped and analyzed for validation in two populations of African Americans (N = 461; N = 778) and two of whites (N = 550; N = 967). Again, SNPs within SLC4A5 were significantly associated with blood pressure levels and hypertension status. While not identifying a single causal DNA sequence variation that is significantly associated with blood pressure levels and hypertension status across all samples, the results further implicate SLC4A5 as a candidate hypertension susceptibility gene, validating previous evidence for one or more genes on chromosome 2 that influence hypertension related phenotypes in the population-at-large. The methodology and results reported provide a case study of one approach for following up the results of genetic linkage analyses to identify genes influencing complex traits. ^

Systematic interpretation of molecular beacon polymerase chain reaction for identifying rpoB mutations in Mycobacterium tuberculosis isolates with mixed resistant and susceptible bacteria

Detection of multidrug-resistant tuberculosis (MDR-TB), a frequent cause of treatment failure, takes 2 or more weeks to identify by culture. RIF-resistance is a hallmark of MDR-TB, and detection of mutations in the rpoB gene of Mycobacterium tuberculosis using molecular beacon probes with real-time quantitative polymerase chain reaction (qPCR) is a novel approach that takes ≤2 days. However, qPCR identification of resistant isolates, particularly for isolates with mixed RIF-susceptible and RIF-resistant bacteria, is reader dependent and limits its clinical use. The aim of this study was to develop an objective, reader-independent method to define rpoB mutants using beacon qPCR. This would facilitate the transition from a research protocol to the clinical setting, where high-throughput methods with objective interpretation are required. For this, DNAs from 107 M. tuberculosis clinical isolates with known susceptibility to RIF by culture-based methods were obtained from 2 regions where isolates have not previously been subjected to evaluation using molecular beacon qPCR: the Texas–Mexico border and Colombia. Using coded DNA specimens, mutations within an 81-bp hot spot region of rpoB were established by qPCR with 5 beacons spanning this region. Visual and mathematical approaches were used to establish whether the qPCR cycle threshold of the experimental isolate was significantly higher (mutant) compared to a reference wild-type isolate. Visual classification of the beacon qPCR required reader training for strains with a mixture of RIF-susceptible and RIF-resistant bacteria. Only then had the visual interpretation by an experienced reader had 100% sensitivity and 94.6% specificity versus RIF-resistance by culture phenotype and 98.1% sensitivity and 100% specificity versus mutations based on DNA sequence. The mathematical approach was 98% sensitive and 94.5% specific versus culture and 96.2% sensitive and 100% specific versus DNA sequence. Our findings indicate the mathematical approach has advantages over the visual reading, in that it uses a Microsoft Excel template to eliminate reader bias or inexperience, and allows objective interpretation from high-throughput analyses even in the presence of a mixture of RIF-resistant and RIF-susceptible isolates without the need for reader training.^

Structure and biosynthesis of a pyruvate -dependent enzyme---phosphatidylserine decarboxylase from {\it Escherichia coli\/}

Phosphatidylserine decarboxylase of E. coli, a cytoplasmic membrane protein, catalyzes the formation of phosphatidylethanolamine, the principal phospholipid of the organism. The activity of the enzyme is dependent on a covalently bound pyruvate (Satre and Kennedy (1978) J. Biol. Chem. 253, 479-483). This study shows that the enzyme consists of two nonidentical subunits, $\alpha$ (Mr = 7,332) and $\beta$ (Mr = 28,579), with the pyruvate prosthetic group in amide linkage to the amino-terminus of the $\alpha$ subunit. Partial protein sequence and DNA sequence analysis reveal that the two subunits are derived from a proenzyme ($\pi$ subunit, Mr = 35,893) through a post-translational event. During the conversion of the proenzyme to the $\alpha$ and $\beta$ subunits, the peptide bond between Gly253-Ser254 is cleaved, and Ser254 is converted to the pyruvate prosthetic group at the amino-terminus of the $\alpha$ subunit (Li and Dowhan (1988) J. Biol. Chem. 263, 11516-11522).^ The proenzyme cannot be detected in cells carrying either single or multiple copies of the gene (psd), but can be observed in a T7 RNA polymerase/promoter and transcription-translation system. The cleavage of the wild-type proenzyme occurs rapidly with a half-time on the order of 2 min. Changing of the Ser254 to cysteine (S254C) or threonine (S254T) slows the cleavage rate dramatically and results in mutants with a half-time for processing of around 2-4 h. Change of the Ser254 to alanine (S254A) blocks the cleavage of the proenzyme. The reduced processing rate with the mutations of the proenzyme is consistent with less of the functional enzyme being made. Mutants S254C and S254T produce $\sim$15% and $\sim$1%, respectively, of the activity of the wild-type allele, but can still complement a temperature-sensitive mutant of the psd locus. Neither detectable activity nor complementation is observed by mutant S254A. These results are consistent with the hydroxyl-group of the Ser254 playing a critical role in the cleavage of the peptide bond Gly253-Ser254 of the pro-phosphatidylserine decarboxylase, and support the mechanism proposed by Snell and co-workers (Recsei and Snell (1984) Annu. Rev. Biochem. 53, 357-387) for the formation of the prosthetic group of pyruvate-dependent decarboxylases. ^

ATXA -dependent gene expression in Bacillus anthracis

The Bacillus anthracis toxin genes, cya, lef , and pag, can be viewed as a regulon, in which transcription of all three genes is activated in trans by the same regulatory gene, atxA, in response to the same signal, CO2. I determined that several phenotypes are associated with the atxA gene. In addition to being toxin-deficient, an atxA -null mutant grows poorly on minimal media and sporulates early compared to the parent strain. Furthermore, an atxA-null mutant has an altered 2-D gel protein profile. I used a genetic approach to find additional atxA-regulated genes. Random transcriptional lacZ fusions were generated in B. anthracis using transposon Tn 917-LTV3. Transposon-insertion libraries were screened for mutants expressing increased β-galactosidase activity in 5% CO2. Introduction of an atxA-null mutation in these mutants revealed that 79% of the CO2-regulated fusions were also atxA-dependent. DNA sequence analysis of transposon insertion sites in mutants carrying CO 2/atxA-regulated fusions revealed ten mutants harboring transposon insertions in loci distinct from the toxin genes. The majority of the tcr (toxin co-regulated) loci mapped within the pXO1 pathogenicity island. These results indicate a clear association of atxA with CO2-enhanced gene expression in B. anthracis and provide evidence that atxA regulates genes other than the structural genes for the anthrax toxin proteins. ^ Characterization of one tcr locus revealed a new regulatory gene, pagR. The pagR gene (300 nt) is located downstream of pag. pagR is cotranscribed with pag and is responsible for autogenous control of the operon. pagR also represses expression of cya and lef. Repression of toxin gene expression by pagR may be mediated by atxA. The steady state level of atxA mRNA is increased in a pagR mutant. Recombinant PagR protein purified from Escherichia coli did not specifically bind the promoter regions of pagA or atxA. An unidentified factor in B. anthracis crude extracts, however, was able to bind the atxA promoter in the absence of PagR or AtxA. These investigations increase our knowledge of virulence regulation in B. anthracis and ultimately will lead to a better understanding of anthrax disease. ^

Expression regulation of AP -2 gamma gene

AP-2γ is a member of the AP-2 transcription factor family, is highly enriched in the trophoblast cell lineage, and is essential for placenta development. In an effort to identify factors regulating AP-2γ gene expression we isolated and characterized the promoter and 5′ flanking region of the mouse and human AP-2γ genes. The transcription start site of the mouse AP-2γ gene was mapped by primer extension and 5′ RACE. Transient gene transfer studies showed that basal promoter activity resides within a highly conserved ∼200 by DNA sequence located immediately upstream of the transcription start site. The conserved region is highly GC-rich and lacks typical TATA or CCAAT boxes. Multiple potential Sp and AP-2 binding sites are clustered within this region. Electrophoretic mobility shift assays demonstrated that Sp1 and Sp3 bind to three sites in the promoter region of the mouse AP-2γ gene. Combined mutation of the three putative Sp sites reduced promoter activity by 80% in trophoblast and non-trophoblast cells, demonstrating the functional importance of these sites in AP-2γ gene expression. ^ Mutational analysis of the 5′-flanking region revealed a 117-bp positive regulatory region of the mouse AP-2γ gene located between −5700 and −5583 upstream of the transcription start site. This 117-bp positive regulatory element provided approximately 7-fold enhancement of reporter gene expression in cultured trophoblast cells. A C/EBP-Sp1 transcription factor-binding module is located in this DNA sequence. Electrophoretic mobility shift assays demonstrated that transcription factors Sp1, Sp3 and C/EBP bind to the enhancer element. Mutation of each protein-binding site reduced the enhanced expression significantly. Mutagenesis assays showed that two other protein-binding sites also contribute to the enhancer activity. In summary, we have shown that Sp1 and Sp3 bind to cis-regulatory elements located in the promoter region and contribute to basal promoter activity. We have identified a 117-bp positive regulatory element of AP-2γ gene, and we have shown that Sp and C/EBP proteins bind to the cis -regulatory elements and contribute to the enhanced gene expression. ^

REPSA: A method for determining the sequence specificity of DNA binding ligands

We developed a novel combinatorial method termed restriction endonuclease protection selection and amplification (REPSA) to identify consensus binding sites of DNA-binding ligands. REPSA uses a unique enzymatic selection based on the inhibition of cleavage by a type IIS restriction endonuclease, an enzyme that cleaves DNA at a site distal from its recognition sequence. Sequences bound by a ligand are protected from cleavage while unprotected sequences are cleaved. This enzymatic selection occurs in solution under mild conditions and is dependant only on the DNA-binding ability of the ligand. Thus, REPSA is useful for a broad range of ligands including all classes of DNA-binding ligands, weakly binding ligands, mixed populations of ligands, and unknown ligands. Here I describe REPSA and the application of this method to select the consensus DNA-binding sequences of three representative DNA-binding ligands; a nucleic acid (triplex-forming single-stranded DNA), a protein (the TATA-binding protein), and a small molecule (Distamycin A). These studies generated new information regarding the specificity of these ligands in addition to establishing their DNA-binding sequences. ^