Genome-wide analysis of core promoter elements from conserved human and mouse orthologous pairs

Physical Peculiarity of Two Sites in Human Promoters: Universality and Diverse Usage in Gene Function

Since the discovery of physical peculiarities around transcription start sites (TSSs) and a site corresponding to the TATA box, research has revealed only the average features of these sites. Unsettled enigmas include the individual genes with these features and whether they relate to gene function. Herein, using 10 physical properties of DNA, including duplex DNA free energy, base stacking energy, protein-induced deformability, and stabilizing energy of Z-DNA, we clarified for the first time that approximately 97% of the promoters of 21,056 human protein-coding genes have distinctive physical properties around the TSS and/or position -27; of these, nearly 65% exhibited such properties at both sites. Furthermore, about 55% of the 21,056 genes had a minimum value of regional duplex DNA free energy within TSS-centered ±300 bp regions. Notably, distinctive physical properties within the promoters and free energies of the surrounding regions separated human protein-coding genes into five groups; each contained specific gene ontology (GO) terms. The group represented by immune response genes differed distinctly from the other four regarding the parameter of the free energies of the surrounding regions. A vital suggestion from this study is that physical-feature-based analyses of genomes may reveal new aspects of the organization and regulation of genes.

Possibilities of skin coat color-dependent risks and risk factors of squamous cell carcinoma and deafness of domestic cats inferred via RNA-seq data

The transcriptome data of skin cells from domestic cats with brown, orange, and white coats were analyzed using a public database to investigate the possible relationship between coat color-related gene expression and squamous cell carcinoma risk, as well as the mechanism of deafness in white cats. We found that the ratio of the expression level of genes suppressing squamous cell carcinoma to that of genes promoting squamous cell carcinoma might be considerably lower than the theoretical estimation in skin cells with orange and white coats in white-spotted cat. We also found the possibility of the frequent production of KIT lacking the first exon (d1KIT) in skin cells with white coats, and d1KIT production exhibited a substantial negative correlation with the expression of SOX10, which is essential for melanocyte formation and adjustment of hearing function. Additionally, the production of d1KIT was expected to be due to the insulating activity of the feline endogenous retrovirus 1 (FERV1) LTR in the first intron of KIT by its CTCF binding sequence repeat. These results contribute to basic veterinary research to understand the relationship between cat skin coat and disease risk, as well as the underlying mechanism.

Investigating the sequence landscape in the Drosophila initiator core promoter element using an enhanced MARZ algorithm

The core promoter elements are important DNA sequences for the regulation of RNA polymerase II transcription in eukaryotic cells. Despite the broad evolutionary conservation of these elements, there is extensive variation in the nucleotide composition of the actual sequences. In this study, we aim to improve our understanding of the complexity of this sequence variation in the TATA box and initiator core promoter elements in Drosophila melanogaster. Using computational approaches, including an enhanced version of our previously developed MARZ algorithm that utilizes gapped nucleotide matrices, several sequence landscape features are uncovered, including an interdependency between the nucleotides in position 2 and 5 in the initiator. Incorporating this information in an expanded MARZ algorithm improves predictive performance for the identification of the initiator element. Overall our results demonstrate the need to carefully consider detailed sequence composition features in core promoter elements in order to make more robust and accurate bioinformatic predictions.

CFA: An explainable deep learning model for annotating the transcriptional roles of cis-regulatory modules based on epigenetic codes

Metazoa gene expression is controlled by modular DNA segments called cis-regulatory modules (CRMs). CRMs can convey promoter/enhancer/insulator roles, generating additional regulation layers in transcription. Experiments for understanding CRM roles are low-throughput and costly. Large-scale CRM function investigation still depends on computational methods. However, existing in silico tools only recognize enhancers or promoters exclusively, thus accumulating errors when considering CRM promoter/enhancer/insulator roles altogether. Currently, no algorithm can concurrently consider these CRM roles. In this research, we developed the CRM Function Annotator (CFA) model. CFA provides complete CRM transcriptional role labeling based on epigenetic profiling interpretation. We demonstrated that CFA achieves high performance (test macro auROC/auPRC = 94.1%/90.3%) and outperforms existing tools in promoter/enhancer/insulator identification. CFA is also inspected to recognize explainable epigenetic codes consistent with previous findings when labeling CRM roles. By considering the higher-order combinations of the epigenetic codes, CFA significantly reduces false-positive rates in CRM transcriptional role annotation. CFA is available at https://github.com/cobisLab/CFA/.

Loss of TAF8 causes TFIID dysfunction and p53-mediated apoptotic neuronal cell death

Mutations in genes encoding general transcription factors cause neurological disorders. Despite clinical prominence, the consequences of defects in the basal transcription machinery during brain development are unclear. We found that loss of the TATA-box binding protein-associated factor TAF8, a component of the general transcription factor TFIID, in the developing central nervous system affected the expression of many, but notably not all genes. Taf8 deletion caused apoptosis, unexpectedly restricted to forebrain regions. Nuclear levels of the transcription factor p53 were elevated in the absence of TAF8, as were the mRNAs of the pro-apoptotic p53 target genes Noxa, Puma and Bax. The cell death in Taf8 forebrain regions was completely rescued by additional loss of p53, but Taf8 and p53 brains failed to initiate a neuronal expression program. Taf8 deletion caused aberrant transcription of promoter regions and splicing anomalies. We propose that TAF8 supports the directionality of transcription and co-transcriptional splicing, and that failure of these processes causes p53-induced apoptosis of neuronal cells in the developing mouse embryo.

Sequence determinants of human gene regulatory elements

DNA can determine where and when genes are expressed, but the full set of sequence determinants that control gene expression is unknown. Here, we measured the transcriptional activity of DNA sequences that represent an ~100 times larger sequence space than the human genome using massively parallel reporter assays (MPRAs). Machine learning models revealed that transcription factors (TFs) generally act in an additive manner with weak grammar and that most enhancers increase expression from a promoter by a mechanism that does not appear to involve specific TF–TF interactions. The enhancers themselves can be classified into three types: classical, closed chromatin and chromatin dependent. We also show that few TFs are strongly active in a cell, with most activities being similar between cell types. Individual TFs can have multiple gene regulatory activities, including chromatin opening and enhancing, promoting and determining transcription start site (TSS) activity, consistent with the view that the TF binding motif is the key atomic unit of gene expression.

Analysis of massively parallel reporter assays measuring the transcriptional activity of DNA sequences indicates that most transcription factor (TF) activity is additive and does not rely on specific TF–TF interactions. Individual TFs can have different gene regulatory activities.

Distinctive regulatory architectures of germline-active and somatic genes in C. elegans

RNA profiling has provided increasingly detailed knowledge of gene expression patterns, yet the different regulatory architectures that drive them are not well understood. To address this, we profiled and compared transcriptional and regulatory element activities across five tissues of Caenorhabditis elegans, covering ∼90% of cells. We find that the majority of promoters and enhancers have tissue-specific accessibility, and we discover regulatory grammars associated with ubiquitous, germline, and somatic tissue–specific gene expression patterns. In addition, we find that germline-active and soma-specific promoters have distinct features. Germline-active promoters have well-positioned +1 and −1 nucleosomes associated with a periodic 10-bp WW signal (W = A/T). Somatic tissue–specific promoters lack positioned nucleosomes and this signal, have wide nucleosome-depleted regions, and are more enriched for core promoter elements, which largely differ between tissues. We observe the 10-bp periodic WW signal at ubiquitous promoters in other animals, suggesting it is an ancient conserved signal. Our results show fundamental differences in regulatory architectures of germline and somatic tissue–specific genes, uncover regulatory rules for generating diverse gene expression patterns, and provide a tissue-specific resource for future studies.

Control of the neuroprotective Lipocalin Apolipoprotein D expression by alternative promoter regions and differentially expressed mRNA 5' UTR variants

The Lipocalin Apolipoprotein D (ApoD) is one of the few genes consistently overexpressed in the aging brain, and in most neurodegenerative and psychiatric diseases. Its functions include metabolism regulation, myelin management, neuroprotection, and longevity regulation. Knowledge of endogenous regulatory mechanisms controlling brain disease-triggered ApoD expression is relevant if we want to boost pharmacologically its neuroprotecting potential. In addition to classical transcriptional control, Lipocalins have a remarkable variability in mRNA 5’UTR-dependent translation efficiency. Using bioinformatic analyses, we uncover strong selective pressures preserving ApoD 5’UTR properties, indicating unexpected functional conservation. PCR amplifications demonstrate the production of five 5’UTR variants (A-E) in mouse ApoD, with diverse expression levels across tissues and developmental stages. Importantly, Variant E is specifically expressed in the oxidative stress-challenged brain. Predictive analyses of 5’UTR secondary structures and enrichment in elements restraining translation, point to Variant E as a tight regulator of ApoD expression. We find two genomic regions conserved in human and mouse ApoD: a canonical (α) promoter region and a previously unknown region upstream of Variant E that could function as an alternative mouse promoter (β). Luciferase assays demonstrate that both α and β promoter regions can drive expression in cultured mouse astrocytes, and that Promoter β activity responds proportionally to incremental doses of the oxidative stress generator Paraquat. We postulate that Promoter β works in association with Variant E 5’UTR as a regulatory tandem that organizes ApoD gene expression in the nervous system in response to oxidative stress, the most common factor in aging and neurodegeneration.

Identification of regulatory elements recapitulating early expression of L-plastin in the zebrafish enveloping layer and embryonic periderm

We have cloned and characterized an intronic fragment of zebrafish lymphocyte cytosolic protein 1 (lcp1, also called L-plastin) that drives expression to the zebrafish enveloping layer (EVL). L-plastin is a calcium-dependent actin-bundling protein belonging to the plastin/fimbrin family of proteins, and is necessary for the proper migration and attachment of several adult cell types, including leukocytes and osteoclasts. However, in zebrafish lcp1 is abundantly expressed much earlier, during differentiation of the EVL. The cells of this epithelial layer migrate collectively, spreading vegetally over the yolk. L-plastin expression persists into the larval periderm, a transient epithelial tissue that forms the first larval skin. This finding establishes that L-plastin is activated in two different embryonic waves, with a distinct regulatory switch between the early EVL and the later leukocyte. To better study L-plastin expressing cells we attempted CRISPR/Cas9 homology-driven recombination (HDR) to insert a self-cleaving peptide (Cre-P2A-EGFP-CAAX) downstream of the native lcp1 promoter. This produced a stable zebrafish line expressing Cre recombinase in EVL nuclei and green fluorescence in EVL cell membranes. In vivo tracking of these labeled cells provided enhanced views of EVL migration behavior, membrane extensions, and mitotic events. Finally, we experimentally dissected key elements of the targeted lcp1 locus, discovering a ∼300 bp intronic sequence sufficient to drive EVL expression. The lcp1: Cre-P2A-EGFP-CAAX zebrafish should be useful for studying enveloping layer specification, gastrulation movements and periderm development in this widely used vertebrate model. In addition, the conserved regulatory sequences we have isolated predict that L-plastin orthologs may have a similar early expression pattern in other vertebrate embryos.

Reporter-ChIP-nexus reveals strong contribution of the Drosophila initiator sequence to RNA polymerase pausing

RNA polymerase II (Pol II) pausing is a general regulatory step in transcription, yet the stability of paused Pol II varies widely between genes. Although paused Pol II stability correlates with core promoter elements, the contribution of individual sequences remains unclear, in part because no rapid assay is available for measuring the changes in Pol II pausing as a result of altered promoter sequences. Here, we overcome this hurdle by showing that ChIP-nexus captures the endogenous Pol II pausing on transfected plasmids. Using this reporter-ChIP-nexus assay in Drosophila cells, we show that the pausing stability is influenced by downstream promoter sequences, but that the strongest contribution to Pol II pausing comes from the initiator sequence, in which a single nucleotide, a G at the +2 position, is critical for stable Pol II pausing. These results establish reporter-ChIP-nexus as a valuable tool to analyze Pol II pausing.

PARP14 Controls the Nuclear Accumulation of a Subset of Type I IFN-Inducible Proteins

The enzymes of the poly-ADP-ribose polymerase (PARP) superfamily control many relevant cellular processes, but a precise understanding of their activities in different physiological or disease contexts is largely incomplete. We found that transcription of several Parp genes was dynamically regulated upon murine macrophage activation by endotoxin. PARP14 was strongly induced by several inflammatory stimuli and translocated into the nucleus of stimulated cells. Quantitative mass spectrometry analysis showed that PARP14 bound to a group of IFN-stimulated gene (ISG)-encoded proteins, most with an unknown function, and it was required for their nuclear accumulation. Moreover, PARP14 depletion attenuated transcription of primary antiviral response genes regulated by the IFN regulatory transcription factor 3, including Ifnb1, thus reducing IFN-β production and activation of ISGs involved in the secondary antiviral response. In agreement with the above-mentioned data, PARP14 hindered Salmonella typhimurium proliferation in murine macrophages. Overall, these data hint at a role of PARP14 in the control of antimicrobial responses and specifically in nuclear activities of a subgroup of ISG-encoded proteins.

Independent regulation of gene expression level and noise by histone modifications

The inherent stochasticity generates substantial gene expression variation among isogenic cells under identical conditions, which is frequently referred to as gene expression noise or cell-to-cell expression variability. Similar to (average) expression level, expression noise is also subject to natural selection. Yet it has been observed that noise is negatively correlated with expression level, which manifests as a potential constraint for simultaneous optimization of both. Here, we studied expression noise in human embryonic cells with computational analysis on single-cell RNA-seq data and in yeast with flow cytometry experiments. We showed that this coupling is overcome, to a certain degree, by a histone modification strategy in multiple embryonic developmental stages in human, as well as in yeast. Importantly, this epigenetic strategy could fit into a burst-like gene expression model: promoter-localized histone modifications (such as H3K4 methylation) are associated with both burst size and burst frequency, which together influence expression level, while gene-body-localized ones (such as H3K79 methylation) are more associated with burst frequency, which influences both expression level and noise. We further knocked out the only “writer” of H3K79 methylation in yeast, and observed that expression noise is indeed increased. Consistently, dosage sensitive genes, such as genes in the Wnt signaling pathway, tend to be marked with gene-body-localized histone modifications, while stress responding genes, such as genes regulating autophagy, tend to be marked with promoter-localized ones. Our findings elucidate that the “division of labor” among histone modifications facilitates the independent regulation of expression level and noise, extend the “histone code” hypothesis to include expression noise, and shed light on the optimization of transcriptome in evolution.

Gene expression noise, or cell-to-cell expression variability, has been a topic of intense interest for more than a decade. The prevailing model of “burst-like transcription” mediated by the promoter transitions between on and off states explains the formation of noise in eukaryotes. Albeit widely accepted, the cis- elements that determine the burst frequency and burst size remain largely unknown. Here we systematically examined the relationship between transcriptional burst frequency/size and all major histone modifications in various cell types, including human embryonic cells, mouse embryonic stem cells, and yeast, and found that histone markers can be divided into two groups based on their associations with burst frequency/ size. Coincidently, promoter-localized histone markers are associated with both burst size and burst frequency whereas gene-body-localized ones are more associated with burst frequency. We further knocked out a gene that is responsible for “writing” a gene-body histone mark in yeast, and found that burst frequency is indeed reduced. Our findings reveal a new mechanism of transcriptional burst regulation and shed light on the simultaneous optimization of gene expression level and noise in evolution.

Core promoter-specific gene regulation: TATA box selectivity and Initiator-dependent bi-directionality of serum response factor-activated transcription

Gene-specific activation by enhancers involves their communication with the basal RNA polymerase II transcription machinery at the core promoter. Core promoters are diverse and may contain a variety of sequence elements such as the TATA box, the Initiator (INR), and the downstream promoter element (DPE) recognized, respectively, by the TATA-binding protein (TBP) and TBP-associated factors of the TFIID complex. Core promoter elements contribute to the gene selectivity of enhancers, and INR/DPE-specific enhancers and activators have been identified. Here, we identify a TATA box-selective activating sequence upstream of the human β-actin (ACTB) gene that mediates serum response factor (SRF)-induced transcription from TATA-dependent but not INR-dependent promoters and requires the TATA-binding/bending activity of TBP, which is otherwise dispensable for transcription from a TATA-less promoter. The SRF-dependent ACTB sequence is stereospecific on TATA promoters but activates in an orientation-independent manner a composite TATA/INR-containing promoter. More generally, we show that SRF-regulated genes of the actin/cytoskeleton/contractile family tend to have a TATA box. These results suggest distinct TATA-dependent and INR-dependent mechanisms of TFIID-mediated transcription in mammalian cells that are compatible with only certain stereospecific combinations of activators, and that a TBP-TATA binding mechanism is important for SRF activation of the actin/cytoskeleton-related gene family.

Focused transcription from the human CR2/CD21 core promoter is regulated by synergistic activity of TATA and Initiator elements in mature B cells

Complement receptor 2 (CR2/CD21) is predominantly expressed on the surface of mature B cells where it forms part of a coreceptor complex that functions, in part, to modulate B-cell receptor signal strength. CR2/CD21 expression is tightly regulated throughout B-cell development such that CR2/CD21 cannot be detected on pre-B or terminally differentiated plasma cells. CR2/CD21 expression is upregulated at B-cell maturation and can be induced by IL-4 and CD40 signaling pathways. We have previously characterized elements in the proximal promoter and first intron of CR2/CD21 that are involved in regulating basal and tissue-specific expression. We now extend these analyses to the CR2/CD21 core promoter. We show that in mature B cells, CR2/CD21 transcription proceeds from a focused TSS regulated by a non-consensus TATA box, an initiator element and a downstream promoter element. Furthermore, occupancy of the general transcriptional machinery in pre-B versus mature B-cell lines correlate with CR2/CD21 expression level and indicate that promoter accessibility must switch from inactive to active during the transitional B-cell window.

TSS seq based core promoter architecture in blood feeding Tsetse fly (Glossina morsitans morsitans) vector of Trypanosomiasis

Background

Transcription initiation regulation is mediated by sequence-specific interactions between DNA-binding proteins (transcription factors) and cis-elements, where BRE, TATA, INR, DPE and MTE motifs constitute canonical core motifs for basal transcription initiation of genes. Accurate identification of transcription start site (TSS) and their corresponding promoter regions is critical for delineation of these motifs. To this end, the genome scale analysis of core promoter architecture in insects has been confined to Drosophila. The recently sequenced Tsetse fly genome provides a unique opportunity to analyze transcription initiation regulation machinery in blood-feeding insects.

Results

A computational method for identification of TSS in newly sequenced Tsetse fly genome was evaluated, using TSS seq tags sampled from two developmental stages namely; larvae and pupae. There were 3134 tag clusters among which 45.4 % (1424) of the tag clusters mapped to first coding exons or their proximal predicted 5′UTR regions and 1.0 % (31) tag clusters mapping to transposons, within a threshold of 100 tags per cluster. These 1393 non transposon-derived core promoters had propensity for AT nucleotides. The −1/+1 and 1/+1 positions in D. melanogaster, and G. m. morsitans had propensity for CA and AA dinucleotides respectively. The 1393 tag clusters comprised narrow promoters (5 %), broad with peak promoters (23 %) and broad without peak promoters (72 %). Two-way motif co-occurrence analysis showed that the MTE-DPE pair is over-represented in broad core promoters. The frequently occurring triplet motifs in all promoter classes are the INR-MTE-DPE, TATA-MTE-DPE and TATA-INR-DPE. Promoters without the TATA motif had higher frequency of the MTE and INR motifs than those observed in Drosophila, where the DPE motif occur more frequently in promoters without TATA motif. Gene ontology terms associated with developmental processes were overrepresented in the narrow and broad with peak promoters.

Conclusions

The study has identified different motif combinations associated with broad promoters in a blood-feeding insect. In the case of TATA-less core promoters, G.m. morsitans uses the MTE to compensate for the lack of a TATA motif. The increasing availability of TSS seq data allows for revision of existing gene annotation datasets with the potential of identifying new transcriptional units.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1921-6) contains supplementary material, which is available to authorized users.

Discovery of transcription start sites in the Chinese hamster genome by next-generation RNA sequencing

Chinese hamster ovary (CHO) cell lines are one of the major production tools for monoclonal antibodies, recombinant proteins, and therapeutics. Although many efforts have significantly improved the availability of sequence information for CHO cells in the last years, forthcoming draft genomes still lack the information depth known from the mouse or human genomes. Many genes annotated for CHO cells and the Chinese hamster reference genome still are in silico predictions, only insufficiently verified by biological experiments. The correct annotation of transcription start sites (TSSs) is of special interest for CHO cells, as these directly define the location of the eukaryotic core promoter. Our study aims to elucidate these largely unexplored regions, trying to shed light on promoter landscapes in the Chinese hamster genome. Based on a 5' enriched dual library RNA sequencing approach 6547 TSSs were identified, of which over 90% were assigned to known genes. These TSSs were used to perform extensive promoter studies using a novel, modular bioinformatics pipeline, incorporating analyses of important regulatory elements of the eukaryotic core promoter on per-gene level and on genomic scale.

Loss of SMAD4 from colorectal cancer cells promotes CCL15 expression to recruit CCR1+ myeloid cells and facilitate liver metastasis

BACKGROUND & AIMS

Loss of the tumor suppressor SMAD4 correlates with progression of colorectal cancer (CRC). In mice, colon tumors that express CCL9 recruit CCR1(+) myeloid cells, which facilitate tumor invasion and metastasis by secreting matrix metalloproteinase 9.

METHODS

We used human CRC cell lines to investigate the ability of SMAD4 to regulate expression of CCL15, a human ortholog of mouse CCL9. We used immunohistochemistry to compare levels of CCL15 and other proteins in 141 samples of human liver metastases.

RESULTS

In human CRC cell lines, knockdown of SMAD4 increased CCL15 expression, and overexpression of SMAD4 decreased it. SMAD4 bound directly to the promoter region of the CCL15 gene to negatively regulate its expression; transforming growth factor-β increased binding of SMAD4 to the CCL15 promoter and transcriptional repression. In livers of nude mice, SMAD4-deficient human CRC cells up-regulated CCL15 to recruit CCR1(+) cells and promote metastasis. In human tumor samples, there was a strong inverse correlation between levels of CCL15 and SMAD4; metastases that expressed CCL15 contained 3-fold more CCR1(+) cells than those without CCL15. Patients with CCL15-expressing metastases had significantly shorter times of disease-free survival than those with CCL15-negative metastases. CCR1(+) cells in the metastases expressed the myeloid cell markers CD11b and myeloperoxidase, and also matrix metalloproteinase 9.

CONCLUSIONS

In human CRC cells, loss of SMAD4 leads to up-regulation of CCL15 expression. Human liver metastases that express CCL15 contain higher numbers CCR1(+) cells; patients with these metastases have shorter times of disease-free survival. Reagents designed to block CCL15 recruitment of CCR1(+) cells could prevent metastasis of CRC to liver.

Genome-wide computational prediction and analysis of core promoter elements across plant monocots and dicots

Transcription initiation, essential to gene expression regulation, involves recruitment of basal transcription factors to the core promoter elements (CPEs). The distribution of currently known CPEs across plant genomes is largely unknown. This is the first large scale genome-wide report on the computational prediction of CPEs across eight plant genomes to help better understand the transcription initiation complex assembly. The distribution of thirteen known CPEs across four monocots (Brachypodium distachyon, Oryza sativa ssp. japonica, Sorghum bicolor, Zea mays) and four dicots (Arabidopsis thaliana, Populus trichocarpa, Vitis vinifera, Glycine max) reveals the structural organization of the core promoter in relation to the TATA-box as well as with respect to other CPEs. The distribution of known CPE motifs with respect to transcription start site (TSS) exhibited positional conservation within monocots and dicots with slight differences across all eight genomes. Further, a more refined subset of annotated genes based on orthologs of the model monocot (O. sativa ssp. japonica) and dicot (A. thaliana) genomes supported the positional distribution of these thirteen known CPEs. DNA free energy profiles provided evidence that the structural properties of promoter regions are distinctly different from that of the non-regulatory genome sequence. It also showed that monocot core promoters have lower DNA free energy than dicot core promoters. The comparison of monocot and dicot promoter sequences highlights both the similarities and differences in the core promoter architecture irrespective of the species-specific nucleotide bias. This study will be useful for future work related to genome annotation projects and can inspire research efforts aimed to better understand regulatory mechanisms of transcription.

Genome-wide organization of eukaryotic preinitiation complex is influenced by nonconsensus protein-DNA binding

Genome-wide binding preferences of the key components of eukaryotic preinitiation complex (PIC) have been recently measured at high resolution in Saccharomyces cerevisiae by Rhee and Pugh. However, the rules determining the PIC binding specificity remain poorly understood. In this study, we show that nonconsensus protein-DNA binding significantly influences PIC binding preferences. We estimate that such nonconsensus binding contributes statistically at least 2-3 kcal/mol (on average) of additional attractive free energy per protein per core-promoter region. The predicted attractive effect is particularly strong at repeated poly(dA:dT) and poly(dC:dG) tracts. Overall, the computed free-energy landscape of nonconsensus protein-DNA binding shows strong correlation with the measured genome-wide PIC occupancy. Remarkably, statistical PIC preferences of binding to both TFIID-dominated and SAGA-dominated genes correlate with the nonconsensus free-energy landscape, yet these two groups of genes are distinguishable based on the average free-energy profiles. We suggest that the predicted nonconsensus binding mechanism provides a genome-wide background for specific promoter elements, such as transcription-factor binding sites, TATA-like elements, and specific binding of the PIC components to nucleosomes. We also show that nonconsensus binding has genome-wide influence on transcriptional frequency.

Novel core promoter elements in the oomycete pathogen Phytophthora infestans and their influence on expression detected by genome-wide analysis

Background

The core promoter is the region flanking the transcription start site (TSS) that directs formation of the pre-initiation complex. Core promoters have been studied intensively in mammals and yeast, but not in more diverse eukaryotes. Here we investigate core promoters in oomycetes, a group within the Stramenopile kingdom that includes important plant and animal pathogens. Prior studies of a small collection of genes proposed that oomycete core promoters contain a 16 to 19 nt motif bearing an Initiator-like sequence (INR) flanked by a novel sequence named FPR, but this has not been extended to whole-genome analysis.

Results

We used expectation maximization to find over-represented motifs near TSSs of Phytophthora infestans, the potato blight pathogen. The motifs corresponded to INR, FPR, and a new element found about 25 nt downstream of the TSS called DPEP. TATA boxes were not detected. Assays of DPEP function by mutagenesis were consistent with its role as a core motif. Genome-wide searches found a well-conserved combined INR+FPR in only about 13% of genes after correcting for false discovery, which contradicted prior reports that INR and FPR are found together in most genes. INR or FPR were found alone near TSSs in 18% and 7% of genes, respectively. Promoters lacking the motifs had pyrimidine-rich regions near the TSS. The combined INR+FPR motif was linked to higher than average mRNA levels, developmentally-regulated transcription, and functions related to plant infection, while DPEP and FPR were over-represented in constitutively-expressed genes. The INR, FPR, and combined INR+FPR motifs were detected in other oomycetes including Hyaloperonospora arabidopsidis, Phytophthora sojae, Pythium ultimum, and Saprolegnia parasitica, while DPEP was found in all but S. parasitica. Only INR seemed present in a non-oomycete stramenopile.

Conclusions

The absence of a TATA box and presence of novel motifs show that the oomycete core promoter is diverged from that of model systems, and likely explains the lack of activity of non-oomycete promoters in Phytophthora transformants. The association of the INR+FPR motif with developmentally-regulated genes shows that oomycete core elements influence stage-specific transcription in addition to regulating formation of the pre-initiation complex.

Genome-wide analysis of Pax8 binding provides new insights into thyroid functions

Background

The transcription factor Pax8 is essential for the differentiation of thyroid cells. However, there are few data on genes transcriptionally regulated by Pax8 other than thyroid-related genes. To better understand the role of Pax8 in the biology of thyroid cells, we obtained transcriptional profiles of Pax8-silenced PCCl3 thyroid cells using whole genome expression arrays and integrated these signals with global cis-regulatory sequencing studies performed by ChIP-Seq analysis

Results

Exhaustive analysis of Pax8 immunoprecipitated peaks demonstrated preferential binding to intragenic regions and CpG-enriched islands, which suggests a role of Pax8 in transcriptional regulation of orphan CpG regions. In addition, ChIP-Seq allowed us to identify Pax8 partners, including proteins involved in tertiary DNA structure (CTCF) and chromatin remodeling (Sp1), and these direct transcriptional interactions were confirmed in vivo. Moreover, both factors modulate Pax8-dependent transcriptional activation of the sodium iodide symporter (Nis) gene promoter. We ultimately combined putative and novel Pax8 binding sites with actual target gene expression regulation to define Pax8-dependent genes. Functional classification suggests that Pax8-regulated genes may be directly involved in important processes of thyroid cell function such as cell proliferation and differentiation, apoptosis, cell polarity, motion and adhesion, and a plethora of DNA/protein-related processes.

Conclusion

Our study provides novel insights into the role of Pax8 in thyroid biology, exerted through transcriptional regulation of important genes involved in critical thyrocyte processes. In addition, we found new transcriptional partners of Pax8, which functionally cooperate with Pax8 in the regulation of thyroid gene transcription. Besides, our data demonstrate preferential location of Pax8 in non-promoter CpG regions. These data point to an orphan CpG island-mediated mechanism that represents a novel role of Pax8 in the transcriptional output of the thyrocyte.

Core promoter-selective function of HMGA1 and Mediator in Initiator-dependent transcription

The factors and mechanisms underlying the differential activity and regulation of eukaryotic RNA polymerase II on different types of core promoters have remained elusive. Here we show that the architectural factor HMGA1 and the Mediator coregulator complex cooperate to enhance basal transcription from core promoters containing both a TATA box and an Initiator (INR) element but not from "TATA-only" core promoters. INR-dependent activation by HMGA1 and Mediator requires the TATA-binding protein (TBP)-associated factors (TAFs) within the TFIID complex and counteracts negative regulators of TBP/TATA-dependent transcription such as NC2 and Topoisomerase I. HMGA1 interacts with TFIID and Mediator and is required for the synergy of TATA and INR elements in mammalian cells. Accordingly, natural HMGA1-activated genes in embryonic stem cells tend to have both TATA and INR elements in a synergistic configuration. Our results suggest a core promoter-specific regulation of Mediator and the basal transcription machinery by HMGA1.

Cell-type specific and combinatorial usage of diverse transcription factors revealed by genome-wide binding studies in multiple human cells

Cell-type diversity is governed in part by differential gene expression programs mediated by transcription factor (TF) binding. However, there are few systematic studies of the genomic binding of different types of TFs across a wide range of human cell types, especially in relation to gene expression. In the ENCODE Project, we have identified the genomic binding locations across 11 different human cell types of CTCF, RNA Pol II (RNAPII), and MYC, three TFs with diverse roles. Our data and analysis revealed how these factors bind in relation to genomic features and shape gene expression and cell-type specificity. CTCF bound predominantly in intergenic regions while RNAPII and MYC preferentially bound to core promoter regions. CTCF sites were relatively invariant across diverse cell types, while MYC showed the greatest cell-type specificity. MYC and RNAPII co-localized at many of their binding sites and putative target genes. Cell-type specific binding sites, in particular for MYC and RNAPII, were associated with cell-type specific functions. Patterns of binding in relation to gene features were generally conserved across different cell types. RNAPII occupancy was higher over exons than adjacent introns, likely reflecting a link between transcriptional elongation and splicing. TF binding was positively correlated with the expression levels of their putative target genes, but combinatorial binding, in particular of MYC and RNAPII, was even more strongly associated with higher gene expression. These data illuminate how combinatorial binding of transcription factors in diverse cell types is associated with gene expression and cell-type specific biology.

An efficient design strategy for a whole-cell biosensor based on engineered ribosome binding sequences

In prokaryotes, the ribosome binding sequence (RBS), located in the 5' untranslated region (5' UTR) of an mRNA, plays a critical role in enhancing mRNA translation and stability. To evaluate the effect of the RBS on the sensitivity and signal intensity of an environmental whole-cell biosensor, three Escherichia coli-based biosensors that respond to benzene, toluene, ethylbenzene, and the xylenes (BTEX) were constructed; the three biosensors have the same Pu promoter and xylR regulator from the Pseudomonas putida TOL plasmid but differ in the engineered RBS in their reporter genes. The results from time and dose-dependent induction of luminescence activity by 2-chlorotoluene showed that the BTEX-SE and BTEX-SD biosensors with engineered RBS had signal intensities approximately 10-35 times higher than the primary BTEX-W biosensor. The limits of detection (LOD) of the BTEX-SE and BTEX-SD biosensors were also significantly lower than the LOD of the BTEX-W biosensor (20 ± 5 μmol L(-1) and 25 ± 5 μmol L(-1) vs. 120 ± 10 μmol L(-1)). Moreover, the BTEX-SE and BTEX-SD biosensors responded three times more rapidly to the analytes. These results suggest that rationally designed RBS in the 5' UTR of a reporter gene may be a promising strategy for increasing the sensitivity, signal intensity, and response speed of whole-cell biosensors.

Comprehensive analysis of transcript start sites in ly49 genes reveals an unexpected relationship with gene function and a lack of upstream promoters

Comprehensive analysis of the transcription start sites of the Ly49 genes of C57BL/6 mice using the oligo-capping 5′-RACE technique revealed that the genes encoding the “missing self” inhibitory receptors, Ly49A, C, G, and I, were transcribed from multiple broad regions in exon 1, in the intron1/exon2 region, and upstream of exon -1b. Ly49E was also transcribed in this manner, and uniquely showed a transcriptional shift from exon1 to exon 2 when NK cells were activated in vitro with IL2. Remarkably, a large proportion of Ly49E transcripts was then initiated from downstream of the translational start codon. By contrast, the genes encoding Ly49B and Q in myeloid cells, the activating Ly49D and H receptors in NK cells, and Ly49F in activated T cells, were predominantly transcribed from a conserved site in a pyrimidine-rich region upstream of exon 1. An ∼200 bp fragment from upstream of the Ly49B start site displayed tissue-specific promoter activity in dendritic cell lines, but the corresponding upstream fragments from all other Ly49 genes lacked detectable tissue-specific promoter activity. In particular, none displayed any significant activity in a newly developed adult NK cell line that expressed multiple Ly49 receptors. Similarly, no promoter activity could be found in fragments upstream of intron1/exon2. Collectively, these findings reveal a previously unrecognized relationship between the pattern of transcription and the expression/function of Ly49 receptors, and indicate that transcription of the Ly49 genes expressed in lymphoid cells is achieved in a manner that does not require classical upstream promoters.

Ensemble approach combining multiple methods improves human transcription start site prediction

Background

The computational prediction of transcription start sites is an important unsolved problem. Some recent progress has been made, but many promoters, particularly those not associated with CpG islands, are still difficult to locate using current methods. These methods use different features and training sets, along with a variety of machine learning techniques and result in different prediction sets.

Results

We demonstrate the heterogeneity of current prediction sets, and take advantage of this heterogeneity to construct a two-level classifier ('Profisi Ensemble') using predictions from 7 programs, along with 2 other data sources. Support vector machines using 'full' and 'reduced' data sets are combined in an either/or approach. We achieve a 14% increase in performance over the current state-of-the-art, as benchmarked by a third-party tool.

Conclusions

Supervised learning methods are a useful way to combine predictions from diverse sources.

Genome-wide mapping of RNA Pol-II promoter usage in mouse tissues by ChIP-seq

Alternative promoters that are differentially used in various cellular contexts and tissue types add to the transcriptional complexity in mammalian genome. Identification of alternative promoters and the annotation of their activity in different tissues is one of the major challenges in understanding the transcriptional regulation of the mammalian genes and their isoforms. To determine the use of alternative promoters in different tissues, we performed ChIP-seq experiments using antibody against RNA Pol-II, in five adult mouse tissues (brain, liver, lung, spleen and kidney). Our analysis identified 38 639 Pol-II promoters, including 12 270 novel promoters, for both protein coding and non-coding mouse genes. Of these, 6384 promoters are tissue specific which are CpG poor and we find that only 34% of the novel promoters are located in CpG-rich regions, suggesting that novel promoters are mostly tissue specific. By identifying the Pol-II bound promoter(s) of each annotated gene in a given tissue, we found that 37% of the protein coding genes use alternative promoters in the five mouse tissues. The promoter annotations and ChIP-seq data presented here will aid ongoing efforts of characterizing gene regulatory regions in mammalian genomes.

Highly redundant function of multiple AT-rich sequences as core promoter elements in the TATA-less RPS5 promoter of Saccharomyces cerevisiae

In eukaryotes, protein-coding genes are transcribed by RNA polymerase II (pol II) together with general transcription factors (GTFs). TFIID, the largest GTF composed of TATA element-binding protein (TBP) and 14 TBP-associated factors (TAFs), plays a critical role in transcription from TATA-less promoters. In metazoans, several core promoter elements other than the TATA element are thought to be recognition sites for TFIID. However, it is unclear whether functionally homologous elements also exist in TATA-less promoters in Saccharomyces cerevisiae. Here, we identify the cis-elements required to support normal levels of transcription and accurate initiation from sites within the TATA-less and TFIID-dependent RPS5 core promoter. Systematic mutational analyses show that multiple AT-rich sequences are required for these activities and appear to function as recognition sites for TFIID. A single copy of these sequences can support accurate initiation from the endogenous promoter, indicating that they carry highly redundant functions. These results show a novel architecture of yeast TATA-less promoters and support a model in which pol II scans DNA downstream from a recruited site, while searching for appropriate initiation site(s).

Characterization of the human PANK2 promoter

Pantothenate kinase 2 (PANK2) is an essential regulatory enzyme in coenzyme A biosynthesis. PANK2 mutations cause pantothenate kinase-associated neurodegeneration (PKAN), which leads to pigmentary retinopathy, progressive dystonia and other abnormalities. Two nearly identical PANK2 isoforms have been described: short PANK2 and mature PANK2, which are processed from a precursor isoform. Since the biological relevance of these isoforms remains unclear, we sought to explore their transcriptional regulation. Here we show that their regulation is distinct and describe a promoter for the short isoform of PANK2. Moreover, we identify potential regulators of PANK2 expression, including NF-Y, FOXN4 and the human heterogeneous nuclear ribonucleoprotein A/B family. These findings validate expression of the short PANK2 isoform and enable predictions about potentially deleterious sequence variants in the regulatory region of this human disease gene.

SCS: signal, context, and structure features for genome-wide human promoter recognition

This paper integrates the signal, context, and structure features for genome-wide human promoter recognition, which is important in improving genome annotation and analyzing transcriptional regulation without experimental supports of ESTs, cDNAs, or mRNAs. First, CpG islands are salient biological signals associated with approximately 50 percent of mammalian promoters. Second, the genomic context of promoters may have biological significance, which is based on n-mers (sequences of n bases long) and their statistics estimated from training samples. Third, sequence-dependent DNA flexibility originates from DNA 3D structures and plays an important role in guiding transcription factors to the target site in promoters. Employing decision trees, we combine above signal, context, and structure features to build a hierarchical promoter recognition system called SCS. Experimental results on controlled data sets and the entire human genome demonstrate that SCS is significantly superior in terms of sensitivity and specificity as compared to other state-of-the-art methods. The SCS promoter recognition system is available online as supplemental materials for academic use and can be found on the Computer Society Digital Library at http://doi.ieeecomputersociety.org/10.1109/TCBB.2008.95.

High-resolution characterization of Toxoplasma gondii transcriptome with a massive parallel sequencing method

For the last couple of years, a method that permits the collection of precise positional information of transcriptional start sites (TSSs) together with digital information of the gene-expression levels in a high-throughput manner was established. We applied this novel method, ‘tss-seq’, to elucidate the transcriptome of tachyzoites of the Toxoplasma gondii, which resulted in the identification of 124 000 TSSs, and they were clustered into 10 000 transcription regions (TRs) with a statistics-based analysis. The TRs and annotated ORFs were paired, resulting in the identification of 30% of the TRs and 40% of the ORFs without their counterparts, which predicted undiscovered genes and stage-specific transcriptions, respectively. The massive data for TSSs make it possible to execute the first systematic analysis of the T. gondii core promoter structure, and the information showed that T. gondii utilized an initiator-like motif for their transcription in the major and novel motif, the downstream thymidine cluster, which was similar to the Y patch observed in plants. This encyclopaedic analysis also suggested that the TATA box, and the other well-known core promoter elements were hardly utilized.

Functional analysis and identification of cis-regulatory elements of human chromosome 21 gene promoters

Given the inherent limitations of in silico studies relying solely on DNA sequence analysis, the functional characterization of mammalian promoters and associated cis-regulatory elements requires experimental support, which demands cloning and analysis of putative promoter regions. Focusing on human chromosome 21, we cloned 182 gene promoters of 2500 bp in length and conducted reporter gene assays on transfected-cell arrays. We found 56 promoters that were active in HEK293 cells, while another 49 promoters could be activated by treatment of cells with Trichostatin A or depletion of serum. We observed high correlations between promoter activities and endogenous transcript levels, RNA polymerase II occupancy, CpG islands and core promoter elements. Truncation of a subset of 62 promoters to ∼500 bp revealed that truncation rarely resulted in loss of activity, but rather in loss of responses to external stimuli, suggesting the presence of cis-regulatory response elements within distal promoter regions. In these regions, we found a strong enrichment of transcription factor binding sites that could potentially activate gene expression in the presence of stimuli. This study illustrates the modular functional architecture of chromosome 21 promoters and helps to reveal the complex mechanisms governing transcriptional regulation.

Sugarcane for bioenergy production: an assessment of yield and regulation of sucrose content

An increasing number of plant scientists, including breeders, agronomists, physiologists and molecular biologists, are working towards the development of new and improved energy crops. Research is increasingly focused on how to design crops specifically for bioenergy production and increased biomass generation for biofuel purposes. The most important biofuel to date is bioethanol produced from sugars (sucrose and starch). Second generation bioethanol is also being targeted for studies to allow the use of the cell wall (lignocellulose) as a source of carbon. If a crop is to be used for bioenergy production, the crop should be high yielding, fast growing, low lignin content and requiring relatively small energy inputs for its growth and harvest. Obtaining high yields in nonprime agricultural land is a key for energy crop development to allow sustainability and avoid competition with food production. Sugarcane is the most efficient bioenergy crop of tropical and subtropical regions, and biotechnological tools for the improvement of this crop are advancing rapidly. We focus this review on the studies of sugarcane genes associated with sucrose content, biomass and cell wall metabolism and the preliminary physiological characterization of cultivars that contrast for sugar and biomass yield.

High-resolution analysis of DNA regulatory elements by synthetic saturation mutagenesis

We present a method that harnesses massively parallel DNA synthesis and sequencing for the high-throughput functional analysis of regulatory sequences at single-nucleotide resolution. As a proof of concept, we quantitatively assayed the effects of all possible single-nucleotide mutations for three bacteriophage promoters and three mammalian core promoters in a single experiment per promoter. The method may also serve as a rapid screening tool for regulatory element engineering in synthetic biology.

Promoting developmental transcription

Animal growth and development depend on the precise control of gene expression at the level of transcription. A central role in the regulation of developmental transcription is attributed to transcription factors that bind DNA enhancer elements, which are often located far from gene transcription start sites. Here, we review recent studies that have uncovered significant regulatory functions in developmental transcription for the TFIID basal transcription factors and for the DNA core promoter elements that are located close to transcription start sites.

Coding region structural heterogeneity and turnover of transcription start sites contribute to divergence in expression between duplicate genes

BACKGROUND

Gene expression divergence is one manifestation of functional differences between duplicate genes. Although rapid accumulation of expression divergence between duplicate gene copies has been observed, the driving mechanisms behind this phenomenon have not been explored in detail.

RESULTS

We examine which factors influence expression divergence between human duplicate genes, utilizing the latest genome-wide data sets. We conclude that the turnover of transcription start sites between duplicate genes occurs rapidly after gene duplication and that gene pairs with shared transcription start sites have significantly higher expression similarity than those without shared transcription start sites. Moreover, we find that most (55%) duplicate gene pairs do not retain the same coding sequence structure between the two duplicate copies and this also contributes to divergence in their expression. Furthermore, the proportion of aligned sequences in cis-regulatory regions between the two copies is positively correlated with expression similarity. Surprisingly, we find no effect of copy-specific transposable element insertions on the divergence of duplicate gene expression.

CONCLUSIONS

Our results suggest that turnover of transcription start sites, structural heterogeneity of coding sequences, and divergence of cis-regulatory regions between copies play a pivotal role in determining the expression divergence of duplicate genes.

Pol II promoter prediction using characteristic 4-mer motifs: a machine learning approach

Background

Eukaryotic promoter prediction using computational analysis techniques is one of the most difficult jobs in computational genomics that is essential for constructing and understanding genetic regulatory networks. The increased availability of sequence data for various eukaryotic organisms in recent years has necessitated for better tools and techniques for the prediction and analysis of promoters in eukaryotic sequences. Many promoter prediction methods and tools have been developed to date but they have yet to provide acceptable predictive performance. One obvious criteria to improve on current methods is to devise a better system for selecting appropriate features of promoters that distinguish them from non-promoters. Secondly improved performance can be achieved by enhancing the predictive ability of the machine learning algorithms used.

Results

In this paper, a novel approach is presented in which 128 4-mer motifs in conjunction with a non-linear machine-learning algorithm utilising a Support Vector Machine (SVM) are used to distinguish between promoter and non-promoter DNA sequences. By applying this approach to plant, Drosophila, human, mouse and rat sequences, the classification model has showed 7-fold cross-validation percentage accuracies of 83.81%, 94.82%, 91.25%, 90.77% and 82.35% respectively. The high sensitivity and specificity value of 0.86 and 0.90 for plant; 0.96 and 0.92 for Drosophila; 0.88 and 0.92 for human; 0.78 and 0.84 for mouse and 0.82 and 0.80 for rat demonstrate that this technique is less prone to false positive results and exhibits better performance than many other tools. Moreover, this model successfully identifies location of promoter using TATA weight matrix.

Conclusion

The high sensitivity and specificity indicate that 4-mer frequencies in conjunction with supervised machine-learning methods can be beneficial in the identification of RNA pol II promoters comparative to other methods. This approach can be extended to identify promoters in sequences for other eukaryotic genomes.

Characterization of the human alpha1 beta1 soluble guanylyl cyclase promoter: key role for NF-kappaB(p50) and CCAAT-binding factors in regulating expression of the nitric oxide receptor

Soluble guanylyl cyclase (sGC) is the principal receptor for NO and plays a ubiquitous role in regulating cellular function. This is exemplified in the cardiovascular system where sGC governs smooth muscle tone and growth, vascular permeability, leukocyte flux, and platelet aggregation. As a consequence, aberrant NO-sGC signaling has been linked to diseases including hypertension, atherosclerosis, and stroke. Despite these key (patho)physiological roles, little is known about the expressional regulation of sGC. To address this deficit, we have characterized the promoter activity of human α₁ and β₁ sGC genes in a cell type relevant to cardiovascular (patho)physiology, primary human aortic smooth muscle cells. Luciferase reporter constructs revealed that the 0.3- and 0.5-kb regions upstream of the transcription start sites were optimal for α₁ and β₁ sGC promoter activity, respectively. Deletion of consensus sites for c-Myb, GAGA, NFAT, NF-κB(p50), and CCAAT-binding factor(s) (CCAAT-BF) revealed that these are the principal transcription factors regulating basal sGC expression. In addition, under pro-inflammatory conditions, the effects of the strongest α₁ and β₁ sGC repressors were enhanced, and enzyme expression and activity were reduced; in particular, NF-κB(p50) is pivotal in regulating enzyme expression under such conditions. NO itself also elicited a cGMP-independent negative feedback effect on sGC promoter activity that is mediated, in part, via CCAAT-BF activity. In sum, these data provide a systematic characterization of the promoter activity of human sGC α₁ and β₁ subunits and identify key transcription factors that govern subunit expression under basal and pro-inflammatory (i.e. atherogenic) conditions and in the presence of ligand NO.

The RNA polymerase II core promoter - the gateway to transcription

The RNA polymerase II core promoter is generally defined to be the sequence that directs the initiation of transcription. This simple definition belies a diverse and complex transcriptional module. There are two major types of core promoters - focused and dispersed. Focused promoters contain either a single transcription start site or a distinct cluster of start sites over several nucleotides, whereas dispersed promoters contain several start sites over 50-100 nucleotides and are typically found in CpG islands in vertebrates. Focused promoters are more ancient and widespread throughout nature than dispersed promoters; however, in vertebrates, dispersed promoters are more common than focused promoters. In addition, core promoters may contain many different sequence motifs, such as the TATA box, BRE, Inr, MTE, DPE, DCE, and XCPE1, that specify different mechanisms of transcription and responses to enhancers. Thus, the core promoter is a sophisticated gateway to transcription that determines which signals will lead to transcription initiation.

Weak correlation between sequence conservation in promoter regions and in protein-coding regions of human-mouse orthologous gene pairs

Background

Interspecies sequence comparison is a powerful tool to extract functional or evolutionary information from the genomes of organisms. A number of studies have compared protein sequences or promoter sequences between mammals, which provided many insights into genomics. However, the correlation between protein conservation and promoter conservation remains controversial.

Results

We examined promoter conservation as well as protein conservation for 6,901 human and mouse orthologous genes, and observed a very weak correlation between them. We further investigated their relationship by decomposing it based on functional categories, and identified categories with significant tendencies. Remarkably, the 'ribosome' category showed significantly low promoter conservation, despite its high protein conservation, and the 'extracellular matrix' category showed significantly high promoter conservation, in spite of its low protein conservation.

Conclusion

Our results show the relation of gene function to protein conservation and promoter conservation, and revealed that there seem to be nonparallel components between protein and promoter sequence evolution.

Taip2 is a novel cell death-related gene expressed in the brain during development

TAIP2 was isolated as one of the homologous genes of TAIP3 (TGF-beta-up-regulated apoptosis-inducing-protein chromosome 3). The transcript of the mouse counterpart of TAIP2, designated mTaip2, was detected in several tissue specimens from embryos to adults, while mTaip2 was dominantly expressed in the embryonic brain. The overexpression of the full-length mTaip2 induced cell death in various cell lines. An analysis of mTaip2 deletion mutants revealed that the N-terminal half of mTaip2, but not the C-terminal half, had nuclear localization and cell death-inducing activities. The results indicate that mTaip2 is a novel cell death-related gene dominantly expressed in the embryonic brain, thus suggesting that mTaip2 may play a role in development of the brain.

Generic eukaryotic core promoter prediction using structural features of DNA

Despite many recent efforts, in silico identification of promoter regions is still in its infancy. However, the accurate identification and delineation of promoter regions is important for several reasons, such as improving genome annotation and devising experiments to study and understand transcriptional regulation. Current methods to identify the core region of promoters require large amounts of high-quality training data and often behave like black box models that output predictions that are difficult to interpret. Here, we present a novel approach for predicting promoters in whole-genome sequences by using large-scale structural properties of DNA. Our technique requires no training, is applicable to many eukaryotic genomes, and performs extremely well in comparison with the best available promoter prediction programs. Moreover, it is fast, simple in design, and has no size constraints, and the results are easily interpretable. We compared our approach with 14 current state-of-the-art implementations using human gene and transcription start site data and analyzed the ENCODE region in more detail. We also validated our method on 12 additional eukaryotic genomes, including vertebrates, invertebrates, plants, fungi, and protists.

Transcription of histone gene cluster by differential core-promoter factors

The 100 copies of tandemly arrayed Drosophila linker (H1) and core (H2A/B and H3/H4) histone gene cluster are coordinately regulated during the cell cycle. However, the molecular mechanisms that must allow differential transcription of linker versus core histones prevalent during development remain elusive. Here, we used fluorescence imaging, biochemistry, and genetics to show that TBP (TATA-box-binding protein)-related factor 2 (TRF2) selectively regulates the TATA-less Histone H1 gene promoter, while TBP/TFIID targets core histone transcription. Importantly, TRF2-depleted polytene chromosomes display severe chromosomal structural defects. This selective usage of TRF2 and TBP provides a novel mechanism to differentially direct transcription within the histone cluster. Moreover, genome-wide chromatin immunoprecipitation (ChIP)-on-chip analyses coupled with RNA interference (RNAi)-mediated functional studies revealed that TRF2 targets several classes of TATA-less promoters of >1000 genes including those driving transcription of essential chromatin organization and protein synthesis genes. Our studies establish that TRF2 promoter recognition complexes play a significantly more central role in governing metazoan transcription than previously appreciated.

Co-regulation of Gremlin and Notch signalling in diabetic nephropathy

Diabetic nephropathy is currently the leading cause of end-stage renal disease worldwide, and occurs in approximately one third of all diabetic patients. The molecular pathogenesis of diabetic nephropathy has not been fully characterized and novel mediators and drivers of the disease are still being described. Previous data from our laboratory has identified the developmentally regulated gene Gremlin as a novel target implicated in diabetic nephropathy in vitro and in vivo. We used bioinformatic analysis to examine whether Gremlin gene sequence and structure could be used to identify other genes implicated in diabetic nephropathy. The Notch ligand Jagged1 and its downstream effector, hairy enhancer of split-1 (Hes1), were identified as genes with significant similarity to Gremlin in terms of promoter structure and predicted microRNA binding elements. This led us to discover that transforming growth factor-beta (TGFbeta1), a primary driver of cellular changes in the kidney during nephropathy, increased Gremlin, Jagged1 and Hes1 expression in human kidney epithelial cells. Elevated levels of Gremlin, Jagged1 and Hes1 were also detected in extracts from renal biopsies from diabetic nephropathy patients, but not in control living donors. In situ hybridization identified specific upregulation and co-expression of Gremlin, Jagged1 and Hes1 in the same tubuli of kidneys from diabetic nephropathy patients, but not controls. Finally, Notch pathway gene clustering showed that samples from diabetic nephropathy patients grouped together, distinct from both control living donors and patients with minimal change disease. Together, these data suggest that Notch pathway gene expression is elevated in diabetic nephropathy, co-incident with Gremlin, and may contribute to the pathogenesis of this disease.

A short ultraconserved sequence drives transcription from an alternate FBN1 promoter

FBN1, the gene mutated in Marfan syndrome, encodes fibrillin-1, a large glycoprotein component of the extracellular microfibrils. Human FBN1 has three untranslated upstream exons, and homologous sequences can be identified in a number of mammalian species. In this work, we have used functional assays to characterize the FBN1 upstream region. Sequences upstream of exon 1 and at least two of the upstream untranslated exons were shown to possess promoter activity in vitro. The strongest activity in luciferase assays was shown for sequences upstream of the untranslated exon A. Sequence analysis of the sequences in and upstream of exon A in humans and six other mammalian species demonstrated several highly conserved potential cis-acting sequences as well as a 66-basepair (bp) ultraconserved sequence with nearly perfect conservation in the seven species. The ultraconserved sequence contains an initiator element (Inr), a downstream promoter element (DPE), and a 10-bp palindromic element. Mutational assays showed that both the Inr and the DPE are critical for full promoter activity. A mutation of the 10-bp palindromic element completely abolished basal promoter activity. The element was shown to bind specifically to an unknown nuclear protein by electrophoretic mobility shift assay. Ultraconservation within an alternate promoter has not been previously reported. We suggest that the ultraconservation may reflect the importance of finely tuned regulation of alternate transcription of FBN1 and that the sequences involved have been under negative selective pressure for at least the last 180 million years of mammalian evolution.

Computational analyses of eukaryotic promoters

Computational analysis of eukaryotic promoters is one of the most difficult problems in computational genomics and is essential for understanding gene expression profiles and reverse-engineering gene regulation network circuits. Here I give a basic introduction of the problem and recent update on both experimental and computational approaches. More details may be found in the extended references. This review is based on a summer lecture given at Max Planck Institute at Berlin in 2005.