Transcription factor map alignment of promoter regions

We address the problem of comparing and characterizing the promoter regions of genes with similar expression patterns. This remains a challenging problem in sequence analysis, because often the promoter regions of co-expressed genes do not show discernible sequence conservation. In our approach, thus, we have not directly compared the nucleotide sequence of promoters. Instead, we have obtained predictions of transcription factor binding sites, annotated the predicted sites with the labels of the corresponding binding factors, and aligned the resulting sequences of labels—to which we refer here as transcription factor maps (TF-maps). To obtain the global pairwise alignment of two TF-maps, we have adapted an algorithm initially developed to align restriction enzyme maps. We have optimized the parameters of the algorithm in a small, but well-curated, collection of human–mouse orthologous gene pairs. Results in this dataset, as well as in an independent much larger dataset from the CISRED database, indicate that TF-map alignments are able to uncover conserved regulatory elements, which cannot be detected by the typical sequence alignments.

Assaying the regulatory potential of mammalian conserved non-coding sequences in human cells.

BACKGROUND: Conserved non-coding sequences in the human genome are approximately tenfold more abundant than known genes, and have been hypothesized to mark the locations of cis-regulatory elements. However, the global contribution of conserved non-coding sequences to the transcriptional regulation of human genes is currently unknown. Deeply conserved elements shared between humans and teleost fish predominantly flank genes active during morphogenesis and are enriched for positive transcriptional regulatory elements. However, such deeply conserved elements account for <1% of the conserved non-coding sequences in the human genome, which are predominantly mammalian. RESULTS: We explored the regulatory potential of a large sample of these 'common' conserved non-coding sequences using a variety of classic assays, including chromatin remodeling, and enhancer/repressor and promoter activity. When tested across diverse human model cell types, we find that the fraction of experimentally active conserved non-coding sequences within any given cell type is low (approximately 5%), and that this proportion increases only modestly when considered collectively across cell types. CONCLUSIONS: The results suggest that classic assays of cis-regulatory potential are unlikely to expose the functional potential of the substantial majority of mammalian conserved non-coding sequences in the human genome.

Functional interactions between the estrogen receptor and the transcription activator Sp1 regulate the estrogen-dependent transcriptional activity of the vitellogenin A1 io promoter.

Two distinct, TATA box-containing promoters regulate the transcriptional activity of the Xenopus vitellogenin A1 gene. These two promoters are of different strength and are separated by 1.8 kilobase pairs of untranslated sequence. Estrogen receptor (ER) and its ligand, 17beta-estradiol, induce the activity of both promoters. The estrogen response elements (EREs) are located proximal to the downstream i promoter while no ERE-like sequences have been identified in the vicinity of the upstream io promoter. We show here, that transcriptional activity of the upstream io promoter is Sp1-dependent. Moreover, we demonstrate that estrogen inducibility of the io promoter results from functional interactions between the io bound Sp1 and the ER bound at the proximity of i. Functional interactions between Sp1 and ER do not require the presence of a TATA box for transcriptional activation, as is demonstrated using the acyl-CoA oxidase promoter. The relative positions that ER and Sp1 occupy with respect to the initiation site determines whether these two transcription activators can synergize for transcription initiation.

Promoter recognition and activation by the global response regulator CbrB in Pseudomonas aeruginosa.

In Pseudomonas aeruginosa, the CbrA/CbrB two-component system is instrumental in the maintenance of the carbon-nitrogen balance and for growth on carbon sources that are energetically less favorable than the preferred dicarboxylate substrates. The CbrA/CbrB system drives the expression of the small RNA CrcZ, which antagonizes the repressing effects of the catabolite repression control protein Crc, an RNA-binding protein. Dicarboxylates appear to cause carbon catabolite repression by inhibiting the activity of the CbrA/CbrB system, resulting in reduced crcZ expression. Here we have identified a conserved palindromic nucleotide sequence that is present in upstream activating sequences (UASs) of promoters under positive control by CbrB and σ(54) RNA polymerase, especially in the UAS of the crcZ promoter. Evidence for recognition of this palindromic sequence by CbrB was obtained in vivo from mutational analysis of the crcZ promoter and in vitro from electrophoretic mobility shift assays using crcZ promoter fragments and purified CbrB protein truncated at the N terminus. Integration host factor (IHF) was required for crcZ expression. CbrB also activated the lipA (lipase) promoter, albeit less effectively, apparently by interacting with a similar but less conserved palindromic sequence in the UAS of lipA. As expected, succinate caused CbrB-dependent catabolite repression of the lipA promoter. Based on these results and previously published data, a consensus CbrB recognition sequence is proposed. This sequence has similarity to the consensus NtrC recognition sequence, which is relevant for nitrogen control.

The major transcription initiation site of the SV40 late promoter is a potent thyroid hormone response element.

Thyroid hormone receptors (TRs) are members of the nuclear hormone receptor superfamily, which act as transcription factors upon binding to specific DNA sequences called thyroid hormone (T3) response elements (TREs). Such elements are found in the upstream regulatory region of promoters as well as in intragenic sequences of T3-responsive genes. In this report, we demonstrate that SV40 late (SVL) promoter activity is strongly down-regulated by TR in the absence of ligand. Addition of T3 releases this repression, but does not further induce SVL promoter activity. Electrophoretic mobility shift analyses reveal a TR binding element that overlaps with the SV40 major late transcription initiation site. This element closely fits the consensus TRE, formed of two hexanucleotides organized in a tandem repeat separated by 4 nt, and is able to confer T3 responsiveness on a heterologous promoter. We further show that, although the presence of TR leads to quantitatively modified expression of an SVL-driven reporter gene, neither displacement of the site of transcription initiation nor modification of the splicing pattern of the primary transcripts occur.

Synergistic action of GA-binding protein and glucocorticoid receptor in transcription from the mouse mammary tumor virus promoter.

B lymphocytes are among the first cells to be infected by mouse mammary tumor virus (MMTV), and they play a crucial role in its life cycle. To study transcriptional regulation of MMTV in B cells, we have analyzed two areas of the long terminal repeat (LTR) next to the glucocorticoid receptor binding site, fp1 (at position -139 to -146 from the cap site) and fp2 (at -157 to -164). Both showed B-cell-specific protection in DNase I in vitro footprinting assays and contain binding sites for Ets transcription factors, a large family of proteins involved in cell proliferation and differentiation and oncogenic transformation. In gel retardation assays, fp1 and fp2 bound the heterodimeric Ets factor GA-binding protein (GABP) present in B-cell nuclear extracts, which was identified by various criteria: formation of dimers and tetramers, sensitivity to pro-oxidant conditions, inhibition of binding by specific antisera, and comigration of complexes with those formed by recombinant GABP. Mutations which prevented complex formation in vitro abolished glucocorticoid-stimulated transcription from an MMTV LTR linked to a reporter gene in transiently transfected B-cell lines, whereas they did not affect the basal level. Exogenously expressed GABP resulted in an increased level of hormone response of the LTR reporter plasmid and produced a synergistic effect with the coexpressed glucocorticoid receptor, indicating cooperation between the two. This is the first example of GABP cooperation with a steroid receptor, providing the opportunity for studying the integration of their intracellular signaling pathways.

Validation of real-time methylation-specific PCR to determine O6-methylguanine-DNA methyltransferase gene promoter methylation in glioma.

Epigenetic silencing of the DNA repair protein O(6)-methylguanine-DNA methyltransferase (MGMT) by promoter methylation predicts successful alkylating agent therapy, such as with temozolomide, in glioblastoma patients. Stratified therapy assignment of patients in prospective clinical trials according to tumor MGMT status requires a standardized diagnostic test, suitable for high-throughput analysis of small amounts of formalin-fixed, paraffin-embedded tumor tissue. A direct, real-time methylation-specific PCR (MSP) assay was developed to determine methylation status of the MGMT gene promoter. Assay specificity was obtained by selective amplification of methylated DNA sequences of sodium bisulfite-modified DNA. The copy number of the methylated MGMT promoter, normalized to the beta-actin gene, provides a quantitative test result. We analyzed 134 clinical glioma samples, comparing the new test with the previously validated nested gel-based MSP assay, which yields a binary readout. A cut-off value for the MGMT methylation status was suggested by fitting a bimodal normal mixture model to the real-time results, supporting the hypothesis that there are two distinct populations within the test samples. Comparison of the tests showed high concordance of the results (82/91 [90%]; Cohen's kappa = 0.80; 95% confidence interval, 0.82-0.95). The direct, real-time MSP assay was highly reproducible (Pearson correlation 0.996) and showed valid test results for 93% (125/134) of samples compared with 75% (94/125) for the nested, gel-based MSP assay. This high-throughput test provides an important pharmacogenomic tool for individualized management of alkylating agent chemotherapy.

Gene Promoter Evolution Targets the Center of the Human Protein Interaction Network

Assessing the contribution of promoters and coding sequences to gene evolution is an important step toward discovering the major genetic determinants of human evolution. Many specific examples have revealed the evolutionary importance of cis-regulatory regions. However, the relative contribution of regulatory and coding regions to the evolutionary process and whether systemic factors differentially influence their evolution remains unclear. To address these questions, we carried out an analysis at the genome scale to identify signatures of positive selection in human proximal promoters. Next, we examined whether genes with positively selected promoters (Prom+ genes) show systemic differences with respect to a set of genes with positively selected protein-coding regions (Cod+ genes). We found that the number of genes in each set was not significantly different (8.1% and 8.5%, respectively). Furthermore, a functional analysis showed that, in both cases, positive selection affects almost all biological processes and only a few genes of each group are located in enriched categories, indicating that promoters and coding regions are not evolutionarily specialized with respect to gene function. On the other hand, we show that the topology of the human protein network has a different influence on the molecular evolution of proximal promoters and coding regions. Notably, Prom+ genes have an unexpectedly high centrality when compared with a reference distribution (P = 0.008, for Eigenvalue centrality). Moreover, the frequency of Prom+ genes increases from the periphery to the center of the protein network (P = 0.02, for the logistic regression coefficient). This means that gene centrality does not constrain the evolution of proximal promoters, unlike the case with coding regions, and further indicates that the evolution of proximal promoters is more efficient in the center of the protein network than in the periphery. These results show that proximal promoters have had a systemic contribution to human evolution by increasing the participation of central genes in the evolutionary process.

Molecular interactions involved in the transactivation of the human T-cell leukemia virus type 1 promoter mediated by Tax and CREB-2 (ATF-4).

The human T-cell leukemia virus type 1 (HTLV-1) Tax protein activates viral transcription through three 21-bp repeats located in the U3 region of the HTLV-1 long terminal repeat and called Tax-responsive elements (TxREs). Each TxRE contains nucleotide sequences corresponding to imperfect cyclic AMP response elements (CRE). In this study, we demonstrate that the bZIP transcriptional factor CREB-2 is able to bind in vitro to the TxREs and that CREB-2 binding to each of the 21-bp motifs is enhanced by Tax. We also demonstrate that Tax can weakly interact with CREB-2 bound to a cellular palindromic CRE motif such as that found in the somatostatin promoter. Mutagenesis of Tax and CREB-2 demonstrates that both N- and C-terminal domains of Tax and the C-terminal region of CREB-2 are required for direct interaction between the two proteins. In addition, the Tax mutant M47, defective for HTLV-1 activation, is unable to form in vitro a ternary complex with CREB-2 and TxRE. In agreement with recent results suggesting that Tax can recruit the coactivator CREB-binding protein (CBP) on the HTLV-1 promoter, we provide evidence that Tax, CREB-2, and CBP are capable of cooperating to stimulate viral transcription. Taken together, our data highlight the major role played by CREB-2 in Tax-mediated transactivation.

Identification of the minimal promoter region of the mouse NKX5-3, a transcription factor implicated in eye development.

Early ocular development is controlled by a complex network of transcription factors, cell cycle regulators, and diffusible signalling molecules. Together, these molecules regulate cell proliferation and apoptosis, and specify retinal fate. NKX5-3 is a homeobox transcription factor implicated in eye development. The analysis of the 5'-flanking region of the mouse Nkx5-3 gene revealed a predicted TATA-less promoter sequence between -416 and -166 of the translation start site. To functionally characterise Nkx5-3 promoter activity, serial deletions of the promoter sequence were introduced in pGL-3 basic vector and promoter activity of these 5'- and 3'-deleted constructions was tested in HeLa and CHO cells. Transactivation assays identified a region between -350 and -296 exhibiting promoter-like activity. Combined analysis by deletions and point mutations showed that this sequence, containing multiple Sp1 binding sites was necessary to promote transcriptional activity. Binding of Sp1 to this region was confirmed by electrophoretic mobility shift assay (EMSA) and chromatin immunoprecipitation, using an antibody specific for Sp1. Altogether, these results demonstrated that the immediate upstream region of Nkx5-3 gene possessed a strong intrinsic promoter activity in vitro, suggesting a potential role in Nkx5-3 transcription in vivo.

A chromosomal insertion toolbox for promoter probing, mutant complementation and pathogenicity studies in Ralstonia solanacearum

We describe here the construction of a delivery system for stable and directed insertion of gene constructs in a permissive chromosomal site of the bacterial wilt pathogen Ralstonia solanacearum. The system consists of a collection of suicide vectors the Ralstonia chromosome (pRC) series that carry an integration element flanked by transcription terminators and two sequences of homology to the chromosome of strain GMI1000, where the integration element is inserted through a double recombination event. Unique restriction enzyme sites and a GATEWAY cassette enable cloning of any promoter::gene combination in the integration element. Variants endowed with different selectable antibiotic resistance genes and promoter::gene combinations are described. We show that the system can be readily used in GMI1000 and adapted to other R. solanacearum strains using an accessory plasmid. We prove that the pRC system can be employed to complement a deletion mutation with a single copy of the native gene, and to measure transcription of selected promoters in monocopy both in vitro and in planta. Finally, the system has been used to purify and study secretion type III effectors. These novel genetic tools will be particularly useful for the construction of recombinant bacteria that maintain inserted genes or reporter fusions in competitive situations (i.e., during plant infection).

Sequencing and promoter analysis of the nifENXorf3orf5fdxAnifQ operon from Azospirillum brasilense Sp7

A 40-kb DNA region containing the major cluster of nif genes has been isolated from the Azospirillum brasilense Sp7 genome. In this region three nif operons have been identified: nifHDKorf1Y, nifENXorf3orf5fdxAnifQ and orf2nifUSVorf4. The operons containing nifENX and nifUSV genes are separated from the structural nifHDKorf1Y operon by about 5 kb and 10 kb, respectively. The present study shows the sequence analysis of the 6045-bp DNA region containing the nifENX genes. The deduced amino acid sequences from the open reading frames were compared to the nif gene products of other diazotrophic bacteria and indicate the presence of seven ORFs, all reading in the same direction as that of the nifHDKorf1Y operon. Consensus sigma54 and NifA-binding sites are present only in the promoter region upstream of the nifE gene. This promoter is activated by NifA protein and is approximately two-times less active than the nifH promoter, as indicated by the ß-galactosidase assays. This result suggests the differential expression of the nif genes and their respective products in Azospirillum.

Control of the rat angiotensin I converting enzyme gene by CRE-like sequences

We characterized the role of potential cAMP-responsive elements (CRE) in basal and in induced angiotensin converting enzyme (ACE) gene promoter activity in order to shed light on the regulation of somatic ACE expression. We identified stimulators and repressors of basal expression between 122 and 288 bp and between 415 and 1303 bp upstream from the transcription start site, respectively, using a rabbit endothelial cell (REC) line. These regions also contained elements associated with the response to 8BrcAMP. When screening for CRE motifs we found pCRE, a proximal sequence between 209 and 222 bp. dCRE, a distal tandem of two CRE-like sequences conserved between rats, mice and humans, was detected between 834 and 846 bp. Gel retardation analysis of nuclear extracts of REC indicated that pCRE and dCRE bind to the same protein complexes as bound by a canonical CRE. Mutation of pCRE and dCRE in REC established the former as a positive element and the latter as a negative element. In 293 cells, a renal cell line, pCRE and dCRE are negative regulators. Co-transfection of ATF-2 or ATF-2 plus c-Jun repressed ACE promoter activity, suggesting that the ACE gene is controlled by cellular stress. Although mapping of cAMP responsiveness was consistent with roles for pCRE and dCRE, mutation analysis indicated that they were not required for cAMP responsiveness. We conclude that the basal activity of the somatic ACE promoter is controlled by proximal and distal CREs that can act as enhancers or repressors depending on the cell context.

Methylation status of ANAPC1, CDKN2A and TP53 promoter genes in individuals with gastric cancer

Gastric cancer is the forth most frequent malignancy and the second most common cause of cancer death worldwide. DNA methylation is the most studied epigenetic alteration, occurring through a methyl radical addition to the cytosine base adjacent to guanine. Many tumor genes are inactivated by DNA methylation in gastric cancer. We evaluated the DNA methylation status of ANAPC1, CDKN2A and TP53 by methylation-specific PCR in 20 diffuse- and 26 intestinal-type gastric cancer samples and 20 normal gastric mucosa in individuals from Northern Brazil. All gastric cancer samples were advanced stage adenocarcinomas. Gastric samples were surgically obtained at the João de Barros Barreto University Hospital, State of Pará, and were stored at -80°C before DNA extraction. Patients had never been submitted to chemotherapy or radiotherapy, nor did they have any other diagnosed cancer. None of the gastric cancer samples presented methylated DNA sequences for ANAPC1 and TP53. CDKN2A methylation was not detected in any normal gastric mucosa; however, the CDKN2A promoter was methylated in 30.4% of gastric cancer samples, with 35% methylation in diffuse-type and 26.9% in intestinal-type cancers. CDKN2A methylation was associated with the carcinogenesis process for ~30% diffuse-type and intestinal-type compared to non-neoplastic samples. Thus, ANAPC1 and TP53 methylation was probably not implicated in gastric carcinogenesis in our samples. CDKN2A can be implicated in the carcinogenesis process of only a subset of gastric neoplasias.

Promoter region sequence differences in the A and G gamma globin genes of Brazilian sickle cell anemia patients

Fetal hemoglobin (HbF), encoded by the HBG2 and HBG1 genes, is the best-known genetic modulator of sickle cell anemia, varying dramatically in concentration in the blood of these patients. This variation is partially associated with polymorphisms located in the promoter region of the HBG2 and HBG1 genes. In order to explore known and unknown polymorphisms in these genes, the sequences of their promoter regions were screened in sickle cell anemia patients and correlated with both their HbF levels and their βS-globin haplotypes. Additionally, the sequences were compared with genes from 2 healthy groups, a reference one (N = 104) and an Afro-descendant one (N = 98), to identify polymorphisms linked to the ethnic background.The reference group was composed by healthy individuals from the general population. Four polymorphisms were identified in the promoter region of HBG2 and 8 in the promoter region of HBG1 among the studied groups. Four novel single nucleotide polymorphisms (SNP) located at positions -324, -317, -309 and -307 were identified in the reference group. A deletion located between -396 and -391 in the HBG2 promoter region and the SNP -271 C→T in the HBG1 promoter region were associated with the Central African Republic βS-globin haplotype. In contrast, the -369 C→G and 309 A→G SNPs in the HBG2 promoter region were correlated to the Benin haplotype. The polymorphisms -396_-391 del HBG2, -369 SNP HBG2 and -271 SNP HBG1 correlated with HbF levels. Hence, we suggest an important role of HBG2 and HBG1 gene polymorphisms on the HbF synthesis.