TiProD: the Tissue-specific Promoter Database

Engineering highly efficient backsplicing and translation of synthetic circRNAs

Circular RNAs (circRNAs) are highly stable RNA molecules that are attractive templates for expression of therapeutic proteins and non-coding RNAs. In eukaryotes, circRNAs are primarily generated by the spliceosome through backsplicing. Here, we interrogate different molecular elements including intron type and length, Alu repeats, internal ribosome entry sites (IRESs), and exon length essential for circRNA formation and exploit this information to engineer robust backsplicing and circRNA expression. Specifically, we leverage the finding that the downstream intron can tolerate large inserts without affecting splicing to achieve tandem expression of backspliced circRNAs and tRNA intronic circRNAs from the same template. Further, truncation of selected intronic regions markedly increased circRNA formation in different cell types in vitro as well as AAV-mediated circRNA expression in cardiac and skeletal muscle tissue in vivo. We also observed that different IRES elements and exon length influenced circRNA expression and translation, revealing an exonic contribution to splicing, as evidenced by different RNA species produced. Taken together, these data provide new insight into improving the design and expression of synthetic circRNAs. When combined with AAV capsid and promoter technologies, the backsplicing introns and IRES elements constituting this modular platform significantly expand the gene expression toolkit.

In vivo Application of the REMOTE-control System for the Manipulation of Endogenous Gene Expression

Here we describe a protocol for implementing the REMOTE-control system (Reversible Manipulation of Transcription at Endogenous loci), which allows for reversible and tunable expression control of an endogenous gene of interest in living model systems. The REMOTE-control system employs enhanced lac repression and tet activation systems to achieve down- or upregulation of a target gene within a single biological system. Tight repression can be achieved from repressor binding sites flexibly located far downstream of a transcription start site by inhibiting transcription elongation. Robust upregulation can be attained by enhancing the transcription of an endogenous gene by targeting tet transcriptional activators to the cognate promoter. This reversible and tunable expression control can be applied and withdrawn repeatedly in organisms. The potency and versatility of the system, as demonstrated for endogenous Dnmt1 here, will allow more precise in vivo functional analyses by enabling investigation of gene function at various expression levels and by testing the reversibility of a phenotype.

DJ1 and microRNA-214 act synergistically to rescue myoblast cells after ischemia/reperfusion injury

Ischemia/reperfusion injury is a tissue injury occurring post-reperfusion of tissues with pre-existing ischemia. A good blood supply to tissues aids in the survival of ischemic tissue, however, due to prolonged ischemia the levels of ATP decrease and pH declines leading to acidosis. Reduced ATP leads to an increase in the AMP/ATP ratio, causing cessation of intracellular calcium transport, hence calcium overload and cell death. In this study, we demonstrate the synergistic and antagonistic effect of DJ1 and microR-214 (miR-214) in rescuing myoblast C2C12 cells after ischemia/reperfusion in an in vitro model. Both DJ1 and miR-214 were cloned into a hypoxic inducible expression cassette and transfected into the C2C12 cells. We showed that DJ1 and miR-214 have synergistic effects in reducing intracellular lactate dehydrogenase and intracellular transient calcium levels after reoxygenation compared to control cells, in addition to reducing cell death via necrosis. Western blotting revealed a decrease in autophagosome formation in LC3II/I ratio and an increase in AKT expression in cells transfected with DJ1 and miR-214. Using quantitative real-time PCR, we demonstrated that DJ1 and miR-214 significantly reduced the expression of pro-apoptotic factors and autophagy compared to control. The results indicated DJ1 is an endogenous oxidative stress molecule and miR-214 is a potent inhibitor of the sodium calcium exchanger channel. DJ1 had the greatest effect to inhibiting mitochondrial cell death pathways by possibly acting as a modulator of autophagy. Additionally, we have concluded that miR-214 has an inhibitory effect on extrinsic cell death pathways such as necrosis and autophagy.

Multiplexed and programmable regulation of gene networks with an integrated RNA and CRISPR/Cas toolkit in human cells

RNA-based regulation and CRISPR/Cas transcription factors (CRISPR-TFs) have the potential to be integrated for the tunable modulation of gene networks. A major limitation of this methodology is that guide RNAs (gRNAs) for CRISPR-TFs can only be expressed from RNA polymerase III promoters in human cells, limiting their use for conditional gene regulation. We present new strategies that enable expression of functional gRNAs from RNA polymerase II promoters and multiplexed production of proteins and gRNAs from a single transcript in human cells. We use multiple RNA regulatory strategies, including RNA-triple-helix structures, introns, microRNAs, and ribozymes, with Cas9-based CRISPR-TFs and Cas6/Csy4-based RNA processing. Using these tools, we efficiently modulate endogenous promoters and implement tunable synthetic circuits, including multistage cascades and RNA-dependent networks that can be rewired with Csy4 to achieve complex behaviors. This toolkit can be used for programming scalable gene circuits and perturbing endogenous networks for biology, therapeutic, and synthetic biology applications.

Generation of a tumor- and tissue-specific episomal non-viral vector system

Background

A key issue for safe and reproducible gene therapy approaches is the autologous and tissue-specific expression of transgenes. Tissue-specific expression in vivo is either achieved by transfer vectors that deliver the gene of interest into a distinct cell type or by use of tissue-specific expression cassettes. Here we present the generation of non-viral, episomally replicating vectors that are able to replicate in a tissue specific manner thus allowing tissue specific transgene expression in combination with episomal replication. The episomal replication of the prototype vector pEPI-1 and its derivatives depends exclusively on a transcription unit starting from a constitutively active promoter extending into the scaffold/matrix attachment region (S/MAR).

Results

Here, we exchanged the constitutive promoter in the pEPI derivative pEPito by the tumor specific alpha fetoprotein (AFP) or the muscle specific smooth muscle 22 (SM22) promoter leading to specific transgene expression in AFP positive human hepatocellular carcinoma (HUH7) and in a SM22 positive cell line, respectively. The incorporation of the hCMV enhancer element into the expression cassette further boosted the expression levels with both promoters. Tissue specific-replication could be exemplary proven for the smooth muscle protein 22 (SM22) promoter in vitro. With the AFP promoter-driven pEPito vector hepatocellular carcinoma-specific expression could be achieved in vivo after systemic vector application together with polyethylenimine as transfection enhancer.

Conclusions

In this study we present an episomal plasmid system designed for tissue specific transgene expression and replication. The human AFP-promoter in combination with the hCMV enhancer element was demonstrated to be a valuable tissue-specific promoter for targeting hepatocellular carcinomas with non-viral gene delivery system, and tissue specific replication could be shown in vitro with the muscle specific SM22 promoter. In combination with appropriate delivery systems, the tissue specific pEPito vector system will allow higher tissue-specificity with less undesired side effects and is suitable for long term transgene expression in vivo within gene therapeutical approaches.

Sequence signatures extracted from proximal promoters can be used to predict distal enhancers

BACKGROUND

Gene expression is controlled by proximal promoters and distal regulatory elements such as enhancers. While the activity of some promoters can be invariant across tissues, enhancers tend to be highly tissue-specific.

RESULTS

We compiled sets of tissue-specific promoters based on gene expression profiles of 79 human tissues and cell types. Putative transcription factor binding sites within each set of sequences were used to train a support vector machine classifier capable of distinguishing tissue-specific promoters from control sequences. We obtained reliable classifiers for 92% of the tissues, with an area under the receiver operating characteristic curve between 60% (for subthalamic nucleus promoters) and 98% (for heart promoters). We next used these classifiers to identify tissue-specific enhancers, scanning distal non-coding sequences in the loci of the 200 most highly and lowly expressed genes. Thirty percent of reliable classifiers produced consistent enhancer predictions, with significantly higher densities in the loci of the most highly expressed compared to lowly expressed genes. Liver enhancer predictions were assessed in vivo using the hydrodynamic tail vein injection assay. Fifty-eight percent of the predictions yielded significant enhancer activity in the mouse liver, whereas a control set of five sequences was completely negative.

CONCLUSIONS

We conclude that promoters of tissue-specific genes often contain unambiguous tissue-specific signatures that can be learned and used for the de novo prediction of enhancers.

Promoter deletion analysis using a dual-luciferase reporter system

Promoter deletion analysis is a useful tool for identifying important regulatory regions involved in transcriptional control of gene expression. In this approach, a series of promoter deletion fragments are fused to a reporter gene, such as chloramphenicol acetyltransferase or luciferase gene in a vector, and then transfected into cells for induction. Screening the expression level of the reporter gene using either a qualitative or a quantitative assay, allows to identify the regulatory regions of interest (e.g., cis-acting elements or enhancer) in the promoter.Luciferase genes have been widely used as reporter genes for their sensitivity and efficiency. Firefly and Renilla luciferases are two commonly used reporters, which oxidize different substrates to generate quantifiable luminescence. Therefore, the enzymatic activities of firefly and Renilla luciferases can be sequentially measured in a single sample by controlling reaction conditions. Here, we describe a dual-luciferase reporter assay, where the promoter of interest is fused to a firefly luciferase reporter and is co-transfected into cells with an internal control vector (pRL-CMV) to express Renilla luciferase. Both the Firefly and Renilla luciferases are measured using a dual-luciferase reporter assay system which improves experimental accuracy.

Identification of key mechanisms controlling gene expression in Leishmania infected macrophages using genome-wide promoter analysis

The present study describes the in silico prediction of the regulatory network of Leishmania infected human macrophages. The construction of the gene regulatory network requires the identification of Transcription Factor Binding Sites (TFBSs) in the regulatory regions (promoters, enhancers) of genes that are regulated upon Leishmania infection. The promoters of human, mouse, rat, dog and chimpanzee genes were first identified in the whole genomes using available experimental data on full length cDNA sequences or deep CAGE tag data (DBTSS, FANTOM3, FANTOM4), mRNA models (ENSEMBL), or using hand annotated data (EPD, TRANSFAC). A phylogenetic footprinting analysis and a MATCH analysis of the promoter sequences were then performed to predict TFBS. Then, an SQL database that integrates all results of promoter analysis as well as other genome annotation information obtained from ENSEMBL, TRANSFAC, TRED (Transcription Regulatory Element Database), ORegAnno and the ENCODE project, was established. Finally publicly available expression data from human Leishmania infected macrophages were analyzed using the genome-wide information on predicted TFBS with the computer system ExPlain™. The gene regulatory network was constructed and activated signal transduction pathways were revealed. The Irak1 pathway was identified as a key pathway regulating gene expression changes in Leishmania infected macrophages.

Human-chimpanzee promoter comparisons: property-conserved evolution?

Identification of different functional elements and their properties is a fundamental need in biomedical research and phylogenetic comparisons of a growing number of sequenced genomes form a solid basis for this task. Most available phylogenetic approaches are focused on searching for individual sequence alterations, responsible for the observed phenotype, or statistically evaluate observed mutations to infer general trends. However, being applied to close genomes such methods suffer from poor statistics of rare mutations and give only (at its best) coarse results concerning the potential functional importance of the nucleotide differences. However, quantifying the changes in physical properties of DNA allows to see the strength of introduced mutations and hence to classify them for further investigations. In this work we present the comparative sequence analysis of two evolutionarily close species-human and chimpanzee. In contrast to previous studies we evaluate changes in melting enthalpy of DNA rather than count nucleotide mismatches. We find that nucleotide mismatches in promoters were apparently introduced in a correlated manner during the course of evolution, so that, for example, the DNA property "melting enthalpy" was retained. Such property conservation of promoters is significantly different from nucleotide conservation, shows significant positional and functional biases, and seems to represent a novel feature of gene regulation.

Creating PWMs of transcription factors using 3D structure-based computation of protein-DNA free binding energies

Background

Knowledge of transcription factor-DNA binding patterns is crucial for understanding gene transcription. Numerous DNA-binding proteins are annotated as transcription factors in the literature, however, for many of them the corresponding DNA-binding motifs remain uncharacterized.

Results

The position weight matrices (PWMs) of transcription factors from different structural classes have been determined using a knowledge-based statistical potential. The scoring function calibrated against crystallographic data on protein-DNA contacts recovered PWMs of various members of widely studied transcription factor families such as p53 and NF-κB. Where it was possible, extensive comparison to experimental binding affinity data and other physical models was made. Although the p50p50, p50RelB, and p50p65 dimers belong to the same family, particular differences in their PWMs were detected, thereby suggesting possibly different in vivo binding modes. The PWMs of p63 and p73 were computed on the basis of homology modeling and their performance was studied using upstream sequences of 85 p53/p73-regulated human genes. Interestingly, about half of the p63 and p73 hits reported by the Match algorithm in the altogether 126 promoters lay more than 2 kb upstream of the corresponding transcription start sites, which deviates from the common assumption that most regulatory sites are located more proximal to the TSS. The fact that in most of the cases the binding sites of p63 and p73 did not overlap with the p53 sites suggests that p63 and p73 could influence the p53 transcriptional activity cooperatively. The newly computed p50p50 PWM recovered 5 more experimental binding sites than the corresponding TRANSFAC matrix, while both PWMs showed comparable receiver operator characteristics.

Conclusions

A novel algorithm was developed to calculate position weight matrices from protein-DNA complex structures. The proposed algorithm was extensively validated against experimental data. The method was further combined with Homology Modeling to obtain PWMs of factors for which crystallographic complexes with DNA are not yet available. The performance of PWMs obtained in this work in comparison to traditionally constructed matrices demonstrates that the structure-based approach presents a promising alternative to experimental determination of transcription factor binding properties.

Finding human promoter groups based on DNA physical properties

DNA rigidity is an important physical property originating from the DNA three-dimensional structure. Although the general DNA rigidity patterns in human promoters have been investigated, their distinct roles in transcription are largely unknown. In this paper, we discover four highly distinct human promoter groups based on similarity of their rigidity profiles. First, we find that all promoter groups conserve relatively rigid DNAs at the canonical TATA box [a consensus TATA(A/T)A(A/T) sequence] position, which are important physical signals in binding transcription factors. Second, we find that the genes activated by each group of promoters share significant biological functions based on their gene ontology annotations. Finally, we find that these human promoter groups correlate with the tissue-specific gene expression.

Meta-coexpression conservation analysis of microarray data: a "subset" approach provides insight into brain-derived neurotrophic factor regulation

Background

Alterations in brain-derived neurotrophic factor (BDNF) gene expression contribute to serious pathologies such as depression, epilepsy, cancer, Alzheimer's, Huntington and Parkinson's disease. Therefore, exploring the mechanisms of BDNF regulation represents a great clinical importance. Studying BDNF expression remains difficult due to its multiple neural activity-dependent and tissue-specific promoters. Thus, microarray data could provide insight into the regulation of this complex gene. Conventional microarray co-expression analysis is usually carried out by merging the datasets or by confirming the re-occurrence of significant correlations across datasets. However, co-expression patterns can be different under various conditions that are represented by subsets in a dataset. Therefore, assessing co-expression by measuring correlation coefficient across merged samples of a dataset or by merging datasets might not capture all correlation patterns.

Results

In our study, we performed meta-coexpression analysis of publicly available microarray data using BDNF as a "guide-gene" introducing a "subset" approach. The key steps of the analysis included: dividing datasets into subsets with biologically meaningful sample content (e.g. tissue, gender or disease state subsets); analyzing co-expression with the BDNF gene in each subset separately; and confirming co- expression links across subsets. Finally, we analyzed conservation in co-expression with BDNF between human, mouse and rat, and sought for conserved over-represented TFBSs in BDNF and BDNF-correlated genes. Correlated genes discovered in this study regulate nervous system development, and are associated with various types of cancer and neurological disorders. Also, several transcription factor identified here have been reported to regulate BDNF expression in vitro and in vivo.

Conclusion

The study demonstrates the potential of the "subset" approach in co-expression conservation analysis for studying the regulation of single genes and proposes novel regulators of BDNF gene expression.

Gene therapy of benign gynecological diseases

Gene therapy is the introduction of genetic material into patient's cells to achieve therapeutic benefit. Advances in molecular biology techniques and better understanding of disease pathogenesis have validated the use of a variety of genes as potential molecular targets for gene therapy based approaches. Gene therapy strategies include: mutation compensation of dysregulated genes; replacement of defective tumor-suppressor genes; inactivation of oncogenes; introduction of suicide genes; immunogenic therapy and antiangiogenesis based approaches. Preclinical studies of gene therapy for various gynecological disorders have not only shown to be feasible, but also showed promising results in diseases such as uterine leiomyomas and endometriosis. In recent years, significant improvement in gene transfer technology has led to the development of targetable vectors, which have fewer side-effects without compromising their efficacy. This review provides an update on developing gene therapy approaches to treat common gynecological diseases such as uterine leiomyoma and endometriosis.

Towards accurate human promoter recognition: a review of currently used sequence features and classification methods

This review describes important advances that have been made during the past decade for genome-wide human promoter recognition. Interest in promoter recognition algorithms on a genome-wide scale is worldwide and touches on a number of practical systems that are important in analysis of gene regulation and in genome annotation without experimental support of ESTs, cDNAs or mRNAs. The main focus of this review is on feature extraction and model selection for accurate human promoter recognition, with descriptions of what they are, what has been accomplished, and what remains to be done.

AthaMap, integrating transcriptional and post-transcriptional data

The AthaMap database generates a map of predicted transcription factor binding sites (TFBS) for the whole Arabidopsis thaliana genome. AthaMap has now been extended to include data on post-transcriptional regulation. A total of 403 173 genomic positions of small RNAs have been mapped in the A. thaliana genome. These identify 5772 putative post-transcriptionally regulated target genes. AthaMap tools have been modified to improve the identification of common TFBS in co-regulated genes by subtracting post-transcriptionally regulated genes from such analyses. Furthermore, AthaMap was updated to the TAIR7 genome annotation, a graphic display of gene analysis results was implemented, and the TFBS data content was increased. AthaMap is freely available at http://www.athamap.de/.

Genome wide prediction of HNF4alpha functional binding sites by the use of local and global sequence context

An application of machine learning algorithms enables prediction of the functional context of transcription factor binding sites in the human genome.

We report an application of machine learning algorithms that enables prediction of the functional context of transcription factor binding sites in the human genome. We demonstrate that our method allowed de novo identification of hepatic nuclear factor (HNF)4α binding sites and significantly improved an overall recognition of faithful HNF4α targets. When applied to published findings, an unprecedented high number of false positives were identified. The technique can be applied to any transcription factor.

Evaluation of promoters for use in tissue-specific gene delivery

Vectors used in gene therapy require an expression cassette. The expression cassette consists of three important components: promoter, therapeutic gene and polyadenylation signal. The promoter is essential to control expression of the therapeutic gene. A tissue-specific promoter is a promoter that has activity in only certain cell types. Use of a tissue-specific promoter in the expression cassette can restrict unwanted transgene expression as well as facilitate persistent transgene expression. Therefore, choosing the correct promoter, especially a tissue-specific promoter, is a major step toward achieving successful therapeutic transgene expression. Ideally, the elements of the natural promoter region, necessary for obtaining the required level of the gene expression while retaining tissue-specificity, should be known. Also, it is important to understand whether interactions occur between the promoter region and the rest of the vector genome that could affect promoter activity and specificity. To assess this, it is helpful to select a suitable vector system that will be used in further gene therapy studies. Second, have one or several candidate tissue-specific promoters available for use. Third, ideally have an in vitro cell model suitable to evaluate tissue-specificity. Fourth, have a convenient in vivo animal model to use. Fifth, select a good reporter gene system. Next, using conventional recombinant DNA techniques create different promoter constructs with the selected vector system. Lastly, have a suitable transfection method to test the plasmid constructs in both the in vitro and the in vivo models.

Integrative content-driven concepts for bioinformatics “beyond the cell”

Bioinformatics has delivered great contributions to genome and genomics research, without which the world-wide success of this and other global ('omics') approaches would not have been possible. More recently, it has developed further towards the analysis of different kinds of networks thus laying the foundation for comprehensive description, analysis and manipulation of whole living systems in modern "systems biology". The next step which is necessary for developing a systems biology that deals with systemic phenomena is to expand the existing and develop new methodologies that are appropriate to characterize intercellular processes and interactions without omitting the causal underlying molecular mechanisms. Modelling the processes on the different levels of complexity involved requires a comprehensive integration of information on gene regulatory events, signal transduction pathways, protein interaction and metabolic networks as well as cellular functions in the respective tissues / organs.

Tissue-specific and cell type–specific RNA interference in vivo

RNA interference (RNAi) is an efficient method for silencing genes in cultured cells. Here we describe a simple RNAi approach for silencing genes in a cell type-specific and tissue-specific way in vivo. The approach, which mimics the means by which naturally occurring 'microRNA's are generated, uses a tissue-specific polymerase II promoter to drive the expression of a short hairpin RNA (shRNA) directed against the gene target. The shRNA is cleaved by ubiquitously expressed endonucleases to form an active small interfering RNA of about 22 nt. As a proof of principle, it has been shown that expression of a shRNA directed against the transcription factor Wilms tumor 1 in transgenic mice reduces that protein specifically in nurse cells in the testis. Our transgenic RNAi approach offers a cost-effective means of rapidly (within months) addressing the function(s) of genes of interest in a wide variety of specific cell types and tissues in mice in vivo.

Beyond microarrays: find key transcription factors controlling signal transduction pathways

Background

Massive gene expression changes in different cellular states measured by microarrays, in fact, reflect just an "echo" of real molecular processes in the cells. Transcription factors constitute a class of the regulatory molecules that typically require posttranscriptional modifications or ligand binding in order to exert their function. Therefore, such important functional changes of transcription factors are not directly visible in the microarray experiments.

Results

We developed a novel approach to find key transcription factors that may explain concerted expression changes of specific components of the signal transduction network. The approach aims at revealing evidence of positive feedback loops in the signal transduction circuits through activation of pathway-specific transcription factors. We demonstrate that promoters of genes encoding components of many known signal transduction pathways are enriched by binding sites of those transcription factors that are endpoints of the considered pathways. Application of the approach to the microarray gene expression data on TNF-alpha stimulated primary human endothelial cells helped to reveal novel key transcription factors potentially involved in the regulation of the signal transduction pathways of the cells.

Conclusion

We developed a novel computational approach for revealing key transcription factors by knowledge-based analysis of gene expression data with the help of databases on gene regulatory networks (TRANSFAC^® and TRANSPATH^®). The corresponding software and databases are available at .

Composite Module Analyst: a fitness-based tool for identification of transcription factor binding site combinations

MOTIVATION

Functionally related genes involved in the same molecular-genetic, biochemical or physiological process are often regulated coordinately. Such regulation is provided by precisely organized binding of a multiplicity of special proteins [transcription factors (TFs)] to their target sites (cis-elements) in regulatory regions of genes. Cis-element combinations provide a structural basis for the generation of unique patterns of gene expression.

RESULTS

Here we present a new approach for defining promoter models based on the composition of TF binding sites and their pairs. We utilize a multicomponent fitness function for selection of the promoter model that fits best to the observed gene expression profile. We demonstrate examples of successful application of the fitness function with the help of a genetic algorithm for the analysis of functionally related or co-expressed genes as well as testing on simulated and permutated data.

AVAILABILITY

The CMA program is freely available for non-commercial users. URL http://www.gene-regulation.com/pub/programs.html#CMAnalyst. It is also a part of the commercial system ExPlain (www.biobase.de) designed for causal analysis of gene expression data..

TRANSFAC and its module TRANSCompel: transcriptional gene regulation in eukaryotes

The TRANSFAC database on transcription factors, their binding sites, nucleotide distribution matrices and regulated genes as well as the complementing database TRANSCompel on composite elements have been further enhanced on various levels. A new web interface with different search options and integrated versions of Match and Patch provides increased functionality for TRANSFAC. The list of databases which are linked to the common GENE table of TRANSFAC and TRANSCompel has been extended by: Ensembl, UniGene, EntrezGene, HumanPSD and TRANSPRO. Standard gene names from HGNC, MGI and RGD, are included for human, mouse and rat genes, respectively. With the help of InterProScan, Pfam, SMART and PROSITE domains are assigned automatically to the protein sequences of the transcription factors. TRANSCompel contains now, in addition to the COMPEL table, a separate table for detailed information on the experimental EVIDENCE on which the composite elements are based. Finally, for TRANSFAC, in respect of data growth, in particular the gain of Drosophila transcription factor binding sites (by courtesy of the Drosophila DNase I footprint database) and of Arabidopsis factors (by courtesy of DATF, Database of Arabidopsis Transcription Factors) has to be stressed. The here described public releases, TRANSFAC 7.0 and TRANSCompel 7.0, are accessible under http://www.gene-regulation.com/pub/databases.html.