Transcriptional Regulatory Motifs Research Articles

Cell-type-specific effects of autism-associated chromosome 15q11.2-13.1 duplications in human brain.

Recurrent copy number variation represents one of the most well-established genetic drivers in neurodevelopmental disorders, including autism spectrum disorder (ASD). Duplication of 15q11.2-13.1 (dup15q) is a well-described neurodevelopmental syndrome that increases the risk of ASD by over 40-fold. However, the effects of this duplication on gene expression and chromatin accessibility in specific cell types in the human brain remain unknown. To identify the cell-type-specific transcriptional and epigenetic effects of dup15q in the human frontal cortex we conducted single-nucleus RNA-sequencing and multi-omic sequencing on dup15q cases (n=6) as well as non-dup15q ASD (n=7) and neurotypical controls (n=7). Cell-type-specific differential expression analysis identified significantly regulated genes, critical biological pathways, and differentially accessible genomic regions. Although there was overall increased gene expression across the duplicated genomic region, cellular identity represented an important factor mediating gene expression changes. Neuronal subtypes, showed greater upregulation of gene expression across a critical region within the duplication as compared to other cell types. Genes within the duplicated region that had high baseline expression in control individuals showed only modest changes in dup15q, regardless of cell type. Of note, dup15q and ASD had largely distinct signatures of chromatin accessibility, but shared the majority of transcriptional regulatory motifs, suggesting convergent biological pathways. However, the transcriptional binding factor motifs implicated in each condition implicated distinct biological mechanisms; neuronal JUN/FOS networks in ASD vs. an inflammatory transcriptional network in dup15q microglia. This work provides a cell-type-specific analysis of how dup15q changes gene expression and chromatin accessibility in the human brain and finds evidence of marked cell-type-specific effects of this genetic driver. These findings have implications for guiding therapeutic development in dup15q syndrome, as well as understanding the functional effects CNVs more broadly in neurodevelopmental disorders.

A quantitative evaluation of topological motifs and their coupling in gene circuit state distributions

One of the major challenges in biology is to understand how gene interactions collaborate to determine overall functions of biological systems. Here, we present a new computational framework that enables systematic, high-throughput, and quantitative evaluation of how small transcriptional regulatory circuit motifs, and their coupling, contribute to functions of a dynamical biological system. We illustrate how this approach can be applied to identify four-node gene circuits, circuit motifs, and motif coupling responsible for various gene expression state distributions, including those derived from single-cell RNA sequencing data. We also identify seven major classes of four-node circuits from clustering analysis of state distributions. The method is applied to establish phenomenological models of gene circuits driving human neuron differentiation, revealing important biologically relevant regulatory interactions. Our study will shed light on a better understanding of gene regulatory mechanisms in creating and maintaining cellular states.

Gliomas display distinct sex-based differential methylation patterns based on molecular subtype.

BackgroundGliomas are the most common type of primary brain tumor and one of many cancers where males are diagnosed with greater frequency than females. However, little is known about the sex-based molecular differences in glioblastomas (GBMs) or lower grade glioma (non-GBM) subtypes. DNA methylation is an epigenetic mechanism involved in regulating gene transcription. In glioma and other cancers, hypermethylation of specific gene promoters downregulates transcription and may have a profound effect on patient outcome. The purpose of this study was to determine if sex-based methylation differences exist in different glioma subtypes.MethodsMolecular and clinical data from glioma patients were obtained from The Cancer Genome Atlas and grouped according to tumor grade and molecular subtype (IDH1/2 mutation and 1p/19q chromosomal deletion). Sex-specific differentially methylated probes (DMPs) were identified in each subtype and further analyzed to determine if they were part of differentially methylated regions (DMRs) or associated with differentially methylated DNA transcription regulatory binding motifs.ResultsAnalysis of methylation data in 4 glioma subtypes revealed unique sets of both sex-specific DMPs and DMRs in each subtype. Motif analysis based on DMP position also identified distinct sex-based sets of DNA-binding motifs that varied according to glioma subtype. Downstream targets of 2 of the GBM-specific transcription binding sites, NFAT5 and KLF6, showed differential gene expression consistent with increased methylation mediating downregulation.ConclusionDNA methylation differences between males and females in 4 glioma molecular subtypes suggest an important, sex-specific role for DNA methylation in epigenetic regulation of gliomagenesis.

Evolution of UCP1 Transcriptional Regulatory Elements Across the Mammalian Phylogeny.

Uncoupling protein 1 (UCP1) permits non-shivering thermogenesis (NST) when highly expressed in brown adipose tissue (BAT) mitochondria. Exclusive to placental mammals, BAT has commonly been regarded to be advantageous for thermoregulation in hibernators, small-bodied species, and the neonates of larger species. While numerous regulatory control motifs associated with UCP1 transcription have been proposed for murid rodents, it remains unclear whether these are conserved across the eutherian mammal phylogeny and hence essential for UCP1 expression. To address this shortcoming, we conducted a broad comparative survey of putative UCP1 transcriptional regulatory elements in 139 mammals (135 eutherians). We find no evidence for presence of a UCP1 enhancer in monotremes and marsupials, supporting the hypothesis that this control region evolved in a stem eutherian ancestor. We additionally reveal that several putative promoter elements (e.g., CRE-4, CCAAT) identified in murid rodents are not conserved among BAT-expressing eutherians, and together with the putative regulatory region (PRR) and CpG island do not appear to be crucial for UCP1 expression. The specificity and importance of the upTRE, dnTRE, URE1, CRE-2, RARE-2, NBRE, BRE-1, and BRE-2 enhancer elements first described from rats and mice are moreover uncertain as these motifs differ substantially—but generally remain highly conserved—in other BAT-expressing eutherians. Other UCP1 enhancer motifs (CRE-3, PPRE, and RARE-3) as well as the TATA box are also highly conserved in nearly all eutherian lineages with an intact UCP1. While these transcriptional regulatory motifs are generally also maintained in species where this gene is pseudogenized, the loss or degeneration of key basal promoter (e.g., TATA box) and enhancer elements in other UCP1-lacking lineages make it unlikely that the enhancer region is pleiotropic (i.e., co-regulates additional genes). Importantly, differential losses of (or mutations within) putative regulatory elements among the eutherian lineages with an intact UCP1 suggests that the transcriptional control of gene expression is not highly conserved in this mammalian clade.

Cas9/gRNA targeted excision of cystic fibrosis-causing deep-intronic splicing mutations restores normal splicing of CFTR mRNA.

Cystic Fibrosis is an autosomal recessive disorder caused by mutations in the CFTR gene. CRISPR mediated, template-dependent homology-directed gene editing has been used to correct the most common mutation, c.1521_1523delCTT / p.Phe508del (F508del) which affects ~70% of individuals, but the efficiency was relatively low. Here, we describe a high efficiency strategy for editing of three different rare CFTR mutations which together account for about 3% of individuals with Cystic Fibrosis. The mutations cause aberrant splicing of CFTR mRNA due to the creation of cryptic splice signals that result in the formation of pseudoexons containing premature stop codons c.1679+1634A>G (1811+1.6kbA>G) and c.3718-2477C>T (3849+10kbC>T), or an out-of-frame 5’ extension to an existing exon c.3140-26A>G (3272-26A>G). We designed pairs of Cas9 guide RNAs to create targeted double-stranded breaks in CFTR either side of each mutation which resulted in high efficiency excision of the target genomic regions via non-homologous end-joining repair. When evaluated in a mini-gene splicing assay, we showed that targeted excision restored normal splicing for all three mutations. This approach could be used to correct aberrant splicing signals or remove disruptive transcription regulatory motifs caused by deep-intronic mutations in a range of other genetic disorders.

Deep intronic mutations and human disease.

Next-generation sequencing has revolutionized clinical diagnostic testing. Yet, for a substantial proportion of patients, sequence information restricted to exons and exon-intron boundaries fails to identify the genetic cause of the disease. Here we review evidence from mRNA analysis and entire genomic sequencing indicating that pathogenic mutations can occur deep within the introns of over 75 disease-associated genes. Deleterious DNA variants located more than 100 base pairs away from exon-intron junctions most commonly lead to pseudo-exon inclusion due to activation of non-canonical splice sites or changes in splicing regulatory elements. Additionally, deep intronic mutations can disrupt transcription regulatory motifs and non-coding RNA genes. This review aims to highlight the importance of studying variation in deep intronic sequence as a cause of monogenic disorders as well as hereditary cancer syndromes.

Cluster analysis of Plasmodium RNA-seq time-course data identifies stage-specific co-regulated biological processes and regulatory elements

In this study, we interpreted RNA-seq time-course data of three developmental stages of Plasmodium species by clustering genes based on similarities in their expression profile without prior knowledge of the gene function. Functional enrichment of clusters of upregulated genes at specific time-points reveals potential targetable biological processes with information on their timings. We identified common consensus sequences that these clusters shared as potential points of coordinated transcriptional control. Five cluster groups showed upregulated profile patterns of biological interest. This included two clusters from the Intraerythrocytic Developmental Cycle (cluster 4 = 16 genes, and cluster 9 = 32 genes), one from the sexual development stage (cluster 2 = 851 genes), and two from the gamete-fertilization stage in the mosquito host (cluster 4 = 153 genes, and cluster 9 = 258 genes). The IDC expressed the least numbers of genes with only 1448 genes showing any significant activity of the 5020 genes (~29%) in the experiment. Gene ontology (GO) enrichment analysis of these clusters revealed a total of 671 uncharacterized genes implicated in 14 biological processes and components associated with these stages, some of which are currently being investigated as drug targets in on-going research. Five putative transcription regulatory binding motifs shared by members of each cluster were also identified, one of which was also identified in a previous study by separate researchers. Our study shows stage-specific genes and biological processes that may be important in antimalarial drug research efforts. In addition, timed-coordinated control of separate processes may explain the paucity of factors in parasites.

Cluster analysis of Plasmodium RNA-seq time-course data identifies stage-specific co-regulated biological processes and regulatory elements.

In this study, we interpreted RNA-seq time-course data of three developmental stages of Plasmodium species by clustering genes based on similarities in their expression profile without prior knowledge of the gene function. Functional enrichment of clusters of upregulated genes at specific time-points reveals potential targetable biological processes with information on their timings. We identified common consensus sequences that these clusters shared as potential points of coordinated transcriptional control. Five cluster groups showed upregulated profile patterns of biological interest. This included two clusters from the Intraerythrocytic Developmental Cycle (cluster 4 = 16 genes, and cluster 9 = 32 genes), one from the sexual development stage (cluster 2 = 851 genes), and two from the gamete-fertilization stage in the mosquito host (cluster 4 = 153 genes, and cluster 9 = 258 genes). The IDC expressed the least numbers of genes with only 1448 genes showing any significant activity of the 5020 genes (~29%) in the experiment. Gene ontology (GO) enrichment analysis of these clusters revealed a total of 671 uncharacterized genes implicated in 14 biological processes and components associated with these stages, some of which are currently being investigated as drug targets in on-going research. Five putative transcription regulatory binding motifs shared by members of each cluster were also identified, one of which was also identified in a previous study by separate researchers. Our study shows stage-specific genes and biological processes that may be important in antimalarial drug research efforts. In addition, timed-coordinated control of separate processes may explain the paucity of factors in parasites.

Isolation and Characterization of Dromedary Camel Coronavirus UAE-HKU23 from Dromedaries of the Middle East: Minimal Serological Cross-Reactivity between MERS Coronavirus and Dromedary Camel Coronavirus UAE-HKU23.

Recently, we reported the discovery of a dromedary camel coronavirus UAE-HKU23 (DcCoV UAE-HKU23) from dromedaries in the Middle East. In this study, DcCoV UAE-HKU23 was successfully isolated in two of the 14 dromedary fecal samples using HRT-18G cells, with cytopathic effects observed five days after inoculation. Northern blot analysis revealed at least seven distinct RNA species, corresponding to predicted subgenomic mRNAs and confirming the core sequence of transcription regulatory sequence motifs as 5′-UCUAAAC-3′ as we predicted previously. Antibodies against DcCoV UAE-HKU23 were detected in 58 (98.3%) and 59 (100%) of the 59 dromedary sera by immunofluorescence and neutralization antibody tests, respectively. There was significant correlation between the antibody titers determined by immunofluorescence and neutralization assays (Pearson coefficient = 0.525, p < 0.0001). Immunization of mice using recombinant N proteins of DcCoV UAE-HKU23 and Middle East respiratory syndrome coronavirus (MERS-CoV), respectively, and heat-inactivated DcCoV UAE-HKU23 showed minimal cross-antigenicity between DcCoV UAE-HKU23 and MERS-CoV by Western blot and neutralization antibody assays. Codon usage and genetic distance analysis of RdRp, S and N genes showed that the 14 strains of DcCoV UAE-HKU23 formed a distinct cluster, separated from those of other closely related members of Betacoronavirus 1, including alpaca CoV, confirming that DcCoV UAE-HKU23 is a novel member of Betacoronavirus 1.

Comprehensive Tissue-Specific Transcriptome Analysis Reveals Distinct Regulatory Programs during Early Tomato Fruit Development.

Fruit formation and early development involve a range of physiological and morphological transformations of the various constituent tissues of the ovary. These developmental changes vary considerably according to tissue type, but molecular analyses at an organ-wide level inevitably obscure many tissue-specific phenomena. We used laser-capture microdissection coupled to high-throughput RNA sequencing to analyze the transcriptome of ovaries and fruit tissues of the wild tomato species Solanum pimpinellifolium. This laser-capture microdissection-high-throughput RNA sequencing approach allowed quantitative global profiling of gene expression at previously unobtainable levels of spatial resolution, revealing numerous contrasting transcriptome profiles and uncovering rare and cell type-specific transcripts. Coexpressed gene clusters linked specific tissues and stages to major transcriptional changes underlying the ovary-to-fruit transition and provided evidence of regulatory modules related to cell division, photosynthesis, and auxin transport in internal fruit tissues, together with parallel specialization of the pericarp transcriptome in stress responses and secondary metabolism. Analysis of transcription factor expression and regulatory motifs indicated putative gene regulatory modules that may regulate the development of different tissues and hormonal processes. Major alterations in the expression of hormone metabolic and signaling components illustrate the complex hormonal control underpinning fruit formation, with intricate spatiotemporal variations suggesting separate regulatory programs.

In silico functional pathway annotation of 86 established prostate cancer risk variants.

Heritability is one of the strongest risk factors of prostate cancer, emphasizing the importance of the genetic contribution towards prostate cancer risk. To date, 86 established prostate cancer risk variants have been identified by genome-wide association studies (GWAS). To determine if these risk variants are located near genes that interact together in biological networks or pathways contributing to prostate cancer initiation or progression, we generated gene sets based on proximity to the 86 prostate cancer risk variants. We took two approaches to generate gene lists. The first strategy included all immediate flanking genes, up- and downstream of the risk variant, regardless of distance from the index variant, and the second strategy included genes closest to the index GWAS marker and to variants in high LD (r2 ≥0.8 in Europeans) with the index variant, within a 100 kb window up- and downstream. Pathway mapping of the two gene sets supported the importance of the androgen receptor-mediated signaling in prostate cancer biology. In addition, the hedgehog and Wnt/β-catenin signaling pathways were identified in pathway mapping for the flanking gene set. We also used the HaploReg resource to examine the 86 risk loci and variants high LD (r2 ≥0.8) for functional elements. We found that there was a 12.8 fold (p = 2.9 x 10-4) enrichment for enhancer motifs in a stem cell line and a 4.4 fold (p = 1.1 x 10-3) enrichment of DNase hypersensitivity in a prostate adenocarcinoma cell line, indicating that the risk and correlated variants are enriched for transcriptional regulatory motifs. Our pathway-based functional annotation of the prostate cancer risk variants highlights the potential regulatory function that GWAS risk markers, and their highly correlated variants, exert on genes. Our study also shows that these genes may function cooperatively in key signaling pathways in prostate cancer biology.

Rational Design for Enhanced Gene Therapy with DNA Transposons

Several genetic diseases caused by deficiencies in enzymes required throughout the body can be treated by enzyme replacement therapy at costs of upward of hundreds of thousands of dollars per year per patient. In the long run, such costs are prohibitive. Gene therapy is therefore an appealing option in that its aim is to treat patients after only a single or very few interventions. Clinical-grade viral vectors can efficiently deliver therapeutic expression cassettes, but such viruses can still be expensive to prepare. By contrast, clinical-grade nonviral nucleic acids are relatively inexpensive to produce but difficult to deliver into cells. In this issue of Molecular Therapy, two reports from the same group describe improvements to two limiting factors of nonviral gene therapy using DNA transposons delivered to the liver.1,2 The authors substantially increased efficiency of delivery of transgenes to the target organ and the amount of therapeutic protein exported per cell following transgene delivery.

Microarray-based Identification of Individual HERV Loci Expression: Application to Biomarker Discovery in Prostate Cancer

The prostate-specific antigen (PSA) is the main diagnostic biomarker for prostate cancer in clinical use, but it lacks specificity and sensitivity, particularly in low dosage values1​​. ‘How to use PSA' remains a current issue, either for diagnosis as a gray zone corresponding to a concentration in serum of 2.5-10 ng/ml which does not allow a clear differentiation to be made between cancer and noncancer2 or for patient follow-up as analysis of post-operative PSA kinetic parameters can pose considerable challenges for their practical application3,4. Alternatively, noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease, e.g. PCA3 in prostate cancer5,6 and to reveal uncharacterized aspects of tumor biology. Moreover, data from the ENCODE project published in 2012 showed that different RNA types cover about 62% of the genome. It also appears that the amount of transcriptional regulatory motifs is at least 4.5x higher than the one corresponding to protein-coding exons. Thus, long terminal repeats (LTRs) of human endogenous retroviruses (HERVs) constitute a wide range of putative/candidate transcriptional regulatory sequences, as it is their primary function in infectious retroviruses. HERVs, which are spread throughout the human genome, originate from ancestral and independent infections within the germ line, followed by copy-paste propagation processes and leading to multicopy families occupying 8% of the human genome (note that exons span 2% of our genome). Some HERV loci still express proteins that have been associated with several pathologies including cancer7-10. We have designed a high-density microarray, in Affymetrix format, aiming to optimally characterize individual HERV loci expression, in order to better understand whether they can be active, if they drive ncRNA transcription or modulate coding gene expression. This tool has been applied in the prostate cancer field (Figure 1).

Microarray-based Identification of Individual HERV Loci Expression: Application to Biomarker Discovery in Prostate Cancer

The prostate-specific antigen (PSA) is the main diagnostic biomarker for prostate cancer in clinical use, but it lacks specificity and sensitivity, particularly in low dosage values1​​. ‘How to use PSA' remains a current issue, either for diagnosis as a gray zone corresponding to a concentration in serum of 2.5-10 ng/ml which does not allow a clear differentiation to be made between cancer and noncancer2 or for patient follow-up as analysis of post-operative PSA kinetic parameters can pose considerable challenges for their practical application3,4. Alternatively, noncoding RNAs (ncRNAs) are emerging as key molecules in human cancer, with the potential to serve as novel markers of disease, e.g. PCA3 in prostate cancer5,6 and to reveal uncharacterized aspects of tumor biology. Moreover, data from the ENCODE project published in 2012 showed that different RNA types cover about 62% of the genome. It also appears that the amount of transcriptional regulatory motifs is at least 4.5x higher than the one corresponding to protein-coding exons. Thus, long terminal repeats (LTRs) of human endogenous retroviruses (HERVs) constitute a wide range of putative/candidate transcriptional regulatory sequences, as it is their primary function in infectious retroviruses. HERVs, which are spread throughout the human genome, originate from ancestral and independent infections within the germ line, followed by copy-paste propagation processes and leading to multicopy families occupying 8% of the human genome (note that exons span 2% of our genome). Some HERV loci still express proteins that have been associated with several pathologies including cancer7-10. We have designed a high-density microarray, in Affymetrix format, aiming to optimally characterize individual HERV loci expression, in order to better understand whether they can be active, if they drive ncRNA transcription or modulate coding gene expression. This tool has been applied in the prostate cancer field (Figure 1).

Statistical analysis of combinatorial transcriptional regulatory motifs in human intron-containing promoter sequences

The combinatorial regulation is an important feature of transcriptional regulation of genes. Although a number of researches have been done, the general properties of combinatorial transcriptional regulatory elements, including their patterns and distributions, are still not well understood. Specially, the understanding for laws of transcriptional regulation involving introns is challenging. Here, we performed statistical analyses on potentially interacting motif pairs in human housekeeping genes and tissue-specific genes of some tissues. We first identified transcriptional regulatory motifs in promoter sequences, including the first introns, using a sequence-based approach which incorporates several different sequence backgrounds. And then we predicted the potentially interacting motif pairs using Hypergeometric distribution and Mann–Whitney U-test. The validation of the predicted interacting motif pairs has been carried out via analyses of orthologous conservation and enrichment. The distribution characteristics of the identified motif pairs were examined using a probability model. We found that interacting motif pairs have some significant preferences in location, distance and orientation in housekeeping and tissue-specific genes. The results suggest that the features of transcriptional regulation in housekeeping genes are different from the tissue-specific genes. Besides, location preferences of the motif pairs provide clues to the way of the introns participating in transcriptional regulation in different tissue-specific genes.

Abstract 4206: Identification of transcriptional regulatory motifs that control gene expression in premenopausal women with previous history of breast cancer

Abstract Background: Age is an important risk factor for Breast Cancer. Older women have a 50 times greater risk for breast cancer as compared to younger women (under the age of 30 years). Aging is associated with activation of specific transcription factor networks. We hypothesized that premature activation of these programs may predispose to cancer. The aims of this study were 1) to identify active transcription factor networks controlling the expression of specific cancer-related genes in the normal breast and 2) to analyze whether premature activation of this network is associated with neoplastic transformation. Methods: Fresh frozen normal breast tissue of women, under the age of 50 years, with and without history of breast cancer was obtained (ncancer=8, nnon cancer=4). Total RNA was extracted, cDNA annealed and Whole-genome DASL assay (Illumina Corp) was conducted to identify genes that correlated with cancer outcome in premenopausal women. Differential expression analysis was performed using Partek® software program to identify the up or down regulated genes in women with cancer compared to those without. Transcription factors affecting the regulation of genes in patients with cancer were predicted using MotifModeler program (Liu et al). The resulting data was confirmed using qPCR analysis. Results: Whole genome DASL analysis showed differential expression of 34 genes (p-value &amp;lt;0.01). MotifModeler analysis predicted 18 transcription factors as major regulators of gene expression in premenopausal women with previous history of breast cancer. JUNB and FOS were the most significantly altered in this comparison. qPCR pre-validated the over expression of JUNB and FOS in association with cancer development. Further validation studies in an independent cohort are ongoing using Formalin-fixed Paraffin-embedded (FFPE) samples. Conclusion: Gene expression analysis identified transcriptional features of cancer outcome. Differential expression of the JUNB and FOS family of proteins, which are part of AP-1 transcription factor complex, may be correlated with neoplastic transformation in premenopausal women. Breast Cancer Program of the IU Simon Cancer Center funded this project. Citation Format: {Authors}. {Abstract title} [abstract]. In: Proceedings of the 103rd Annual Meeting of the American Association for Cancer Research; 2012 Mar 31-Apr 4; Chicago, IL. Philadelphia (PA): AACR; Cancer Res 2012;72(8 Suppl):Abstract nr 4206. doi:1538-7445.AM2012-4206

Promoter activation by ACGT in response to salicylic and abscisic acids is differentially regulated by the spacing between two copies of the motif

A variety of small sequence motifs, located upstream of minimal promoters influence gene expression. The function of the core sequence ACGT, present in a family of commonly occurring cis acting transcription regulatory motifs was investigated in the background of an artificially designed synthetic promoter sequence. The ACGT was placed in one or two copies and separated by different spacer lengths between the two copies, to study their affect on the expression of basal promoter in plant cells. The activation of transcription by the ACGT element was examined by transient and stable transformation in tobacco, using gusA as the reporter gene. The analysis shows that the expression of the reporter gene was influenced differently by spacing between two adjacent copies of the motif. Two copies of the ACGT element separated by 5 nucleotides gave highest activation. This configuration imparted salicylic acid inducibility to the basal promoter. However, two copies of ACGT separated by 25 nucleotides allowed the promoter to be induced by abscisic acid but not salicylic acid. Computational analysis of the Arabidopsis thaliana genome database showed the presence of these motifs in several genes associated with a variety of stress responses. The results on motif-related inducibility by salicylic acid and abscisic acid, as seen in a synthetic sequence background were validated by experiments on the expression of the native promoter of a protein phosphatase 2C-like gene of Arabidopsis in tobacco leaves. The study supports the importance of spacing between the ACGT sequence motifs on the elicitor specific modulation of gene expression and demonstrates the role of different ACGT motifs in the regulation of PP2C gene expression.

Identification of Potential E2F Target Genes Through cis-regulatory Modules Derived From Chromatin Immunoprecipitation Microarray Data

In the postgenome era, employment of high-throughput data via the integrated use of resources from various domains will lead to the generation of new knowledge and testable hypothesis. In this bioinformatic research, we utilized published chromatin immunoprecipitation microarray (ChIP-chip) results for the promoter binding by E2F1 or E2F4 proteins observed in the primary human WI-38 cells to infer the potential transcription factor binding sites (TFBS). We have compiled “gene vs . motif” and “motif vs . gene” tables from more than 2,700,000 computational predicted transcriptional regulatory motifs representing the regulatory potential for 230 transcription factors families within the proximal promoter sequence (1,200 nucleotides) of human genome. From this approach, for the first time, the transcription of 23 genes is predicted to be under the control of a cis -regulatory module containing four TFBS motifs (CREB, E2F, NF-Y and Nrf-1).

Analysis of transcriptional synergy between upstream regions and introns in ribosomal protein genes of yeast

Transcriptional regulation in eukaryotic genes generally requires combinatorial binding on DNA of multiple transcription factors. Though many analyses have been performed for identification of combinatorial patterns in promoter sequences, there are few studies concerned with introns of genes. Here our study focuses on the transcriptional synergistic (cooperative) regulation between upstream promoters and introns of ribosomal protein (RP) genes in Saccharomyces cerevisiae yeast. We first extract some potential transcriptional regulatory motifs based on a statistical comparative analysis. 98% of these motifs are accordance with experimental analyses. Then by pairing these motifs each other, we identify some potential synergistic motif pairs between upstream regions and introns of yeast RP genes (RPGs). Among 48 detected motif pairs, 44 match the binding sites for interacting transcriptional factors known from experiments or predictions. Checking the positions of these motif pairs in yeast RPGs, it is found that both motifs of the detected motif pairs are enriched in specific regions of upstream regions and introns, respectively. Some motif pairs present distance and orientation preferences, which may be favorable for transcription factors to bind simultaneously to DNA. These results will be helpful to understand the mechanism of synergistic regulation in yeast RPGs.

Spatial Epigenetic Control of Mono- and Bistable Gene Expression

Author SummaryIn the simplest scenario, a gene is expressed when an activator protein binds to its regulatory sequence, and is silenced when the regulatory sequence is bound by a repressor. Many genes are regulated by both activators and repressors, with the response determined by the combined influence of both factors. When the response is monostable graded, expression is finely tuned to a level that reflects the proportion of the bound activator to the bound repressor. Monostable graded systems allow cells to respond precisely to stimuli. If the response is bistable, the response of each cell depends on whether the activator or the repressor wins. Bistable regulation results in the same gene being expressed in some cells and silenced in others, an outcome that promotes cellular differentiation. It remains unclear, however, how different genetic regulatory structures code for monostable graded and bistable responses. We modeled mathematically the behavior of repressors as they bind to and spread their inhibitory effect along genes and found that the spatial distribution of the binding sites determines which response is chosen. If repressors bind both upstream and downstream of the coding sequence, the response is bistable. If they bind only to one side of the coding sequence, the response is monostable. We confirmed our theoretical findings using synthetic genetic constructs in yeast. These findings help to explain how variations in the location of regulatory elements can lead to cellular differentiation and adaption to varying environments.