Search

HNF4α, a member of the nuclear receptor superfamily, regulates the genes involved in lipid and glucose metabolism. The expression of the RARβ gene in the liver of HNF4α knock-out mice was higher versus wildtype controls, whereas oppositely, RARβ promoter activity was 50% reduced by the overexpression of HNF4α in HepG2 cells, and treatment with retinoic acid (RA), a major metabolite of vitamin A, increased RARβ promoter activity 15-fold. The human RARβ2 promoter contains two DR5 and one DR8 binding motifs, as RA response elements (RARE) proximal to the transcription start site. While DR5 RARE1 was previously reported to be responsive to RARs but not to other nuclear receptors, we show here that mutation in DR5 RARE2 suppresses the promoter response to HNF4α and RARα/RXRα. Mutational analysis of ligand-binding pocket amino acids shown to be critical for fatty acid (FA) binding indicated that RA may interfere with interactions of FA carboxylic acid headgroups with side chains of S190 and R235, and the aliphatic group with I355. These results could explain the partial suppression of HNF4α transcriptional activation toward gene promoters that lack RARE, including APOC3 and CYP2C9, while conversely, HNF4α may bind to RARE sequences in the promoter of the genes such as CYP26A1 and RARβ, activating these genes in the presence of RA. Thus, RA could act as either an antagonist towards HNF4α in genes lacking RAREs, or as an agonist for RARE-containing genes. Overall, RA may interfere with the function of HNF4α and deregulate HNF4α targets genes, including the genes important for lipid and glucose metabolism.

Background Palustrin-2CE2 and brevinin-2CE3 are antimicrobial peptides from Rana chensinensis. In R. chensinensis tadpoles, the expression of prepropalustrin-2CE2 and preprobrevinin-2CE3 increased with the developmental stage. In addition, the expression of the two genes was dramatically upregulated with stimulation by Escherichia coli, Staphylococcus aureus, and the chemical lipopolysaccharide (LPS). The genomic organization of the two antimicrobial peptide genes was confirmed. Both prepropalustrin-2CE2 and preprobrevinin-2CE3 contain three exons separated by two large introns. Additionally, several presumed transcription factor binding sites were identified in the promoter sequence. Functional analysis of the promoter was performed using a luciferase reporter system, and further confirmed by yeast one-hybrid experiment and EMSA assay. The results indicated that the transcription factors NF-κB and RelA are involved in regulating the expression of prepropalustrin-2CE2 and preprobrevinin-2CE3. As amphibian populations decline globally, this study provides new data demonstrating how frogs defend against pathogens from the environment by regulating AMP expression. For amphibians, antimicrobial peptides are innate immune molecules that resist adverse external environmental stimuli. However, the regulation mechanism of antimicrobial peptide gene expression in frogs is still unclear. Objective The two antimicrobial peptides, palustrin-2CE2 and brevinin-2CE3, are produced under external stimulation in Rana chensinensis. Using this model, we analyzed the gene structure and regulatory elements of the two antimicrobial peptide genes and explored the regulatory effects of related transcription factors on the two genes. Method Different stimuli such as E. coli, S. aureus, and chemical substance lipopolysaccharide (LPS) were applied to Rana chensinensis tadpoles at different developmental stages, and antimicrobial peptide expression levels were detected by RT-PCR. Bioinformatics analysis and 5'-RACE and genome walking technologies were employed to analyze the genome structure and promoter region of the antimicrobial peptide genes. With dual-luciferase reporter gene assays, yeast one-hybrid experiment and EMSA assays, we assessed the regulatory effect of the endogenous regulators of the cell on the antimicrobial peptide promoter. Results The transcription levels of prepropalustrin-2CE2 and preprobrevinin-2CE3 were significantly upregulated after different stimulations. Genomic structure analysis showed that both genes contained three exons and two introns. Promoter analysis indicated that there are binding sites for regulatory factors of the NF-κB family in the promoter region, and experiments showed that endogenous NF-κB family regulatory factors in frog cells activate the promoters of the antimicrobial peptide genes. Yeast one-hybrid experiment and EMSA assay demonstrated that RelA and NF-κB1 might interact with specific motifs in the prepropalustrin-2CE2 promoter. Conclusion In this paper, we found that the gene expression levels of the antimicrobial peptides, palustrin-2CE2 and brevinin-2CE3, in R. chensinensis will increase under environmental stimuli, and we verified that the changes in gene expression levels are affected by the transcription factors RelA and NF-κB1. The yeast one-hybrid experiment and EMSA assay confirmed that RelA and NF-κB1 could directly interact with the frog antimicrobial peptide gene promoter, providing new data for the regulatory mechanism of antimicrobial peptides in response to environmental stimuli.

Summary Chromatin immunoprecipitation-sequencing is widely used to find transcription factor binding sites, but suffers from various sources of noise. Knocking out the target factor mitigates noise by acting as a negative control. Paired wild-type and knockout (KO) experiments can generate improved motifs but require optimal differential analysis. We introduce peaKO—a computational method to automatically optimize motif analyses with KO controls, which we compare to two other methods. PeaKO often improves elucidation of the target factor and highlights the benefits of KO controls, which far outperform input controls. Availability and implementation PeaKO is freely available at https://peako.hoffmanlab.org. Contact michael.hoffman@utoronto.ca

DNA methylation is a well-characterized epigenetic modification involved in numerous molecular and cellular functions. Methylation patterns have also been associated with aging mechanisms. However, how DNA methylation patterns change within key brain regions involved in memory formation in an age- and sex-specific manner remains unclear. Here, we performed reduced representation bisulfite sequencing (RRBS) from mouse dorsal hippocampus - which is necessary for the formation and consolidation of specific types of memories - in young and aging mice of both sexes. Overall, our findings demonstrate that methylation levels within the dorsal hippocampus are divergent between sexes during aging in genomic features correlating to mRNA functionality, transcription factor binding sites, and gene regulatory elements. These results define age-related changes in the methylome across genomic features and build a foundation for investigating potential target genes regulated by DNA methylation in an age- and sex-specific manner.

PRODORIC is worldwide one of the largest collections of prokaryotic transcription factor binding sites from multiple bacterial sources with corresponding interpretation and visualization tools. With the introduction of PRODORIC2 in 2017, the transition to a modern web interface and maintainable backend was started. With this latest PRODORIC release the database backend is now fully API-based and provides programmatical access to the complete PRODORIC data. The visualization tools Genome Browser and ProdoNet from the original PRODORIC have been reintroduced and were integrated into the PRODORIC website. Missing input and output options from the original Virtual Footprint were added again for position weight matrix pattern-based searches. The whole PRODORIC dataset was reannotated. Every transcription factor binding site was re-evaluated to increase the overall database quality. During this process, additional parameters, like bound effectors, regulation type and different types of experimental evidence have been added for every transcription factor. Additionally, 109 new transcription factors and 6 new organisms have been added. PRODORIC is publicly available at https://www.prodoric.de.

Pathogens alter defense gene responses, genome stability, and DNA repair machinery of plants. The regulation of defense gene expression relies much on pathogenesis-related transcription factors (PRTFs). Binding sites for PRTFs belonging to the families such as MYB, WRKY, NAC, and bZIP are present in the promoters of defense genes such as PATHOGENESIS-RELATED (PR) genes. It could be possible that these PRTFs also regulate the expression of genes unrelated to defense, such as DNA repair genes. Toward this end, this work aims at analyzing the existence and distribution of the putative binding sites of 26 PRTFs in the promoters of 122 DNA repair genes, belonging to 15 categories, of Arabidopsis thaliana. All the DNA repair genes analyzed had at least one putative PRTF binding site in their promoters. Distribution of these sites in the promoters of PR and DNA repair genes revealed much similarity. An analysis of previous transcriptome datasets of these DNA repair genes under various biotic stresses revealed that all of them showed altered expression under any one biotic stress condition. These findings justify the reason for the presence of PRTFs in the DNA repair gene promoters as well as altered genome instability in plants due to pathogen attack.

Transcriptome profiling of maize inbred and hybrid seedlings subjected to heat or cold stress was used to identify key cis-regulatory elements and develop models to predict gene expression responses., Changes in gene expression are important for responses to abiotic stress. Transcriptome profiling of heat- or cold-stressed maize genotypes identifies many changes in transcript abundance. We used comparisons of expression responses in multiple genotypes to identify alleles with variable responses to heat or cold stress and to distinguish examples of cis- or trans-regulatory variation for stress-responsive expression changes. We used motifs enriched near the transcription start sites (TSSs) for thermal stress-responsive genes to develop predictive models of gene expression responses. Prediction accuracies can be improved by focusing only on motifs within unmethylated regions near the TSS and vary for genes with different dynamic responses to stress. Models trained on expression responses in a single genotype and promoter sequences provided lower performance when applied to other genotypes but this could be improved by using models trained on data from all three genotypes tested. The analysis of genes with cis-regulatory variation provides evidence for structural variants that result in presence/absence of transcription factor binding sites in creating variable responses. This study provides insights into cis-regulatory motifs for heat- and cold-responsive gene expression and defines a framework for developing models to predict expression responses across multiple genotypes.

Inositol and its’ metabolic derivatives such as galactinol, raffinose, methylated inositols etc. are well known significant osmoprotectants produced in response to various abiotic stresses in plants. The present work focuses on the importance of phylogenetic analysis in search of a presumptive ‘master switch’ involved in inducing a set of stress-responsive coregulated genes in the pathways involved in inositol synthesis and its further metabolism. It can be reasonably assumed that the entire network of genes induced under a particular stress is governed by a few common stress-responsive transcription factor binding sites or cis elements in their respective promoters. A phylogenetic analysis of such regulatory regions of the coregulated genes participating in entire inositol metabolic pathway throughout evolution, can form a novel approach in recognizing the conserved stress-regulatory signature(s) or ‘master switch(es)’ regulating this important pathway under the stress impact.

El texto completo de este trabajo no está disponible en el Repositorio Académico UPC por restricciones de la casa editorial donde ha sido publicado. During the life cycle of enteric bacterium Escherichia coli, it encounters a wide spectrum of pH changes. The asymmetric dimer of the cAMP receptor protein, CRP, plays a key role in regulating the expressions of genes and the survival of E. coli. To elucidate the pH effects on the mechanism of signal transmission, we present a combination of results derived from ITC, crystallography, and computation. CRP responds to a pH change by inducing a differential effect on the affinity for the binding events to the two cAMP molecules, ensuing in a reversible conversion between positive and negative cooperativity at high and low pH, respectively. The structures of four crystals at pH ranging from 7.8 to 6.5 show that CRP responds by inducing a differential effect on the structures of the two subunits, particularly in the DNA binding domain. Employing the COREX/BEST algorithm, computational analysis shows the change in the stability of residues at each pH. The change in residue stability alters the connectivity between residues including those in cAMP and DNA binding sites. Consequently, the differential impact on the topology of the connectivity surface among residues in adjacent subunits is the main reason for differential change in affinity; that is, the pH-induced differential change in residue stability is the biothermodynamic basis for the change in allosteric behavior. Furthermore, the structural asymmetry of this homodimer amplifies the differential impact of any perturbations. Hence, these results demonstrate that the combination of these approaches can provide insights into the underlying mechanism of an apparent complex allostery signal and transmission in CRP. National Institutes of Health Revisión por pares

Hemophilia B is a monogenic X-linked recessive disorder associated with blood clotting reduction caused by mutations in the clotting factor IX gene (F9). An unusual form of this disorder, hemophilia B Leyden, caused by mutations in the promoter region of the F9 gene, is described. This form of hemophilia B is characterized by a significant, within norm, increase in factor IX in the patient’s plasma after puberty. Mutations in the promoter are grouped in discrete regions of transcription factor binding sites and alter the F9 gene expression during the ontogeny in different ways. This disorder demonstrates how single-nucleotide pathogenic variants in cis-regulatory regions can have a dramatic impact on gene expression. In this article, we present the current data on the F9 gene promoter structure and the results of examining three patients with hemophilia B Leyden who underwent genetic testing in the DNA-diagnostics laboratory of the Bochkov Research Centre for Medical Genetics.

To more sustainably mitigate the impact of crop diseases on plant health and productivity, there is a need for broader spectrum, long-lasting resistance traits. Defense response (DR) genes, located throughout the genome, participate in cellular and system-wide defense mechanisms to stave off infection by diverse pathogens. This multigenic resistance avoids rapid evolution of a pathogen to overcome host resistance. DR genes reside within resistance-associated quantitative trait loci (QTL), and alleles of DR genes in resistant varieties are more active during pathogen attack relative to susceptible haplotypes. Differential expression of DR genes results from polymorphisms in their regulatory regions, that includes cis-regulatory elements such as transcription factor binding sites as well as features that influence epigenetic structural changes to modulate chromatin accessibility during infection. Many of these elements are found in clusters, known as cis-regulatory modules (CRMs), which are distributed throughout the host genome. Regulatory regions involved in plant-pathogen interactions may also contain pathogen effector binding elements that regulate DR gene expression, and that, when mutated, result in a change in the plants' response. We posit that CRMs and the multiple regulatory elements that comprise them are potential targets for marker-assisted breeding for broad-spectrum, durable disease resistance.

Background We observe the emergence of several promoter-variant viral strains in India during recent years. The variant viral promoters contain additional copies of transcription factor binding sites present in the viral modulatory region or enhancer, including RBEIII, LEF-1, Ap-1 and/or NF-κB. These sites are crucial for governing viral gene expression and latency. Here, we infer that one variant viral promoter R2N3-LTR containing two copies of RBF-2 binding sites (an RBEIII site duplication) and three copies of NF-κB motifs may demonstrate low levels of gene expression noise as compared to the canonical RN3-LTR or a different variant R2N4-LTR (a duplication of an RBEIII site and an NF-κB motif). To demonstrate this, we constructed a panel of sub-genomic viral vectors of promoter-variant LTRs co-expressing two reporter proteins (mScarlet and Gaussia luciferase) under the dual-control of Tat and Rev. We established stable pools of CEM.NKR-CCR5 cells (CEM-CCR5RL reporter cells) and evaluated reporter gene expression under different conditions of cell activation. Results The R2N3-LTR established stringent latency that was highly resistant to reversal by potent cell activators such as TNF-α or PMA, or even to a cocktail of activators, compared to the canonical RN3- or the variant R2N4-LTR. The R2N3-LTR exhibited low-level basal gene expression in the absence of cell activation that enhanced marginally but significantly when activated. In the presence of Tat and Rev, trans-complemented in the form of an infectious virus, the R2N3-LTR demonstrated gene expression at levels comparable to the wild-type viral promoter. The R2N3-LTR is responsive to Tat and Rev factors derived from viral strains representing diverse genetic subtypes. Conclusion With extremely low-level transcriptional noise, the R2N3-LTR can serve as an excellent model to examine the establishment, maintenance, and reversal of HIV-1 latency. The R2N3-LTR would also be an ideal viral promoter to develop high-throughput screening assays to identify potent latency-reversing agents since the LTR is not affected by the usual background noise of the cell.

The rapid development of high-throughput sequencing technology provides unique opportunities for studying of transcription factor binding sites, but also brings new computational challenges. Recently, a series of discriminative motif discovery (DMD) methods have been proposed and offer promising solutions for addressing these challenges. However, because of the huge computation cost, most of them have to choose approximate schemes that either sacrifice the accuracy of motif representation or tune motif parameter indirectly. In this paper, we propose a bag-based classifier combined with a multi-fold learning scheme (BCMF) to discover motifs from ChIP-seq datasets. First, BCMF formulates input sequences as a labeled bag naturally. Then, a bag-based classifier, combining with a bag feature extracting strategy, is applied to construct the objective function, and a multi-fold learning scheme is used to solve it. Compared with the existing DMD tools, BCMF features three improvements: 1) Learning position weight matrix (PWM) directly in a continuous space; 2) Proposing to represent a positive bag with a feature fused by its k “most positive” patterns. 3) Applying a more advanced learning scheme. The experimental results on 134 ChIP-seq datasets show that BCMF substantially outperforms existing DMD methods (including DREME, HOMER, XXmotif, motifRG, EDCOD and our previous work).

A data set of promoter and 5′UTR sequences of homoeo‐alleles of 459 wheat genes that contribute to agriculturally important traits in 95 ancestral and commercial wheat cultivars is presented here. The high‐stringency myBaits technology used made individual capture of homoeo‐allele promoters possible, which is reported here for the first time. Promoters of most genes are remarkably conserved across the 83 hexaploid cultivars used with

BACKGROUND: MERTK encodes a receptor tyrosine kinase that regulates immune homeostasis via phagocytosis of apoptotic cells and cytokine-mediated immunosuppression. MERTK is highly expressed in the central nervous system (CNS), specifically in myeloid derived innate immune cells and its dysregulation is implicated in CNS pathologies including the autoimmune disease multiple sclerosis (MS). OBJECTIVE: While the cell types and tissues that express MERTK have been well described, the genetic elements that define the gene’s promoter and regulate specific transcription domains remain unknown. The primary objective of this study was to define and characterise the human MERTK promoter region. METHODS: We cloned and characterized the 5’ upstream region of MERTK to identify cis-acting DNA elements that promote gene transcription in luciferase reporter assays. In addition, promoter regions were tested for sensitivity to the anti-inflammatory glucocorticoid dexamethasone. RESULTS: This study identified identified both proximal and distal-acting DNA elements that promote transcription. The strongest promoter activity was identified in an ∼850 bp region situated 3 kb upstream of the MERTK transcription start site. Serial deletions of this putative enhancer revealed that the entire region is essential for expression activity. Using in silico analysis, we identified several candidate transcription factor binding sites. Despite a well-established upregulation of MERTK in response to anti-inflammatory glucocorticoids, no DNA region within the 5 kb putative promoter was found to directly respond to dexamethasone treatment. CONCLUSIONS: Elucidating the genetic mechanisms that regulate MERTK expression gives insights into gene regulation during homeostasis and disease, providing potential targets for therapeutic modulation of MERTK transcription.

To investigate how exogenous DNA concatemerizes to form episomal artificial chromosomes (ACs), acquire equal segregation ability and maintain stable holocentromeres, we injected DNA sequences with different features, including sequences that are repetitive or complex, and sequences with different AT-contents, into the gonad of Caenorhabditis elegans to form ACs in embryos, and monitored AC mitotic segregation. We demonstrated that AT-poor sequences (26% AT-content) delayed the acquisition of segregation competency of newly formed ACs. We also co-injected fragmented Saccharomyces cerevisiae genomic DNA, differentially expressed fluorescent markers and ubiquitously expressed selectable marker to construct a less repetitive, more complex AC. We sequenced the whole genome of a strain which propagates this AC through multiple generations, and de novo assembled the AC sequences. We discovered CENP-AHCP-3 domains/peaks are distributed along the AC, as in endogenous chromosomes, suggesting a holocentric architecture. We found that CENP-AHCP-3 binds to the unexpressed marker genes and many fragmented yeast sequences, but is excluded in the yeast extremely high-AT-content centromeric and mitochondrial DNA (> 83% AT-content) on the AC. We identified A-rich motifs in CENP-AHCP-3 domains/peaks on the AC and on endogenous chromosomes, which have some similarity with each other and similarity to some non-germline transcription factor binding sites.

Mutation in the germline is the ultimate source of genetic variation, but little is known about the influence of germline chromatin structure on mutational processes. Using ATAC-seq, we profile the open chromatin landscape of human spermatogonia, the most proliferative cell type of the germline, identifying transcription factor binding sites (TFBSs) and PRDM9 binding sites, a subset of which will initiate meiotic recombination. We observe an increase in rare structural variant (SV) breakpoints at PRDM9-bound sites, implicating meiotic recombination in the generation of structural variation. Many germline TFBSs, such as NRF1, are also associated with increased rates of SV breakpoints, apparently independent of recombination. Singleton short insertions (≥5 bp) are highly enriched at TFBSs, particularly at sites bound by testis active TFs, and their rates correlate with those of structural variant breakpoints. Short insertions often duplicate the TFBS motif, leading to clustering of motif sites near regulatory regions in this male-driven evolutionary process. Increased mutation loads at germline TFBSs disproportionately affect neural enhancers with activity in spermatogonia, potentially altering neurodevelopmental regulatory architecture. Local chromatin structure in spermatogonia is thus pervasive in shaping both evolution and disease.

Background Sleep-disordered breathing is a common disorder associated with significant morbidity. The genetic architecture of sleep-disordered breathing remains poorly understood. Through the NHLBI Trans-Omics for Precision Medicine (TOPMed) program, we performed the first whole-genome sequence analysis of sleep-disordered breathing. Methods The study sample was comprised of 7988 individuals of diverse ancestry. Common-variant and pathway analyses included an additional 13,257 individuals. We examined five complementary traits describing different aspects of sleep-disordered breathing: the apnea-hypopnea index, average oxyhemoglobin desaturation per event, average and minimum oxyhemoglobin saturation across the sleep episode, and the percentage of sleep with oxyhemoglobin saturation Results We identified a multi-ethnic set-based rare-variant association (p = 3.48 × 10−8) on chromosome X with ARMCX3. Additional rare-variant associations include ARMCX3-AS1, MRPS33, and C16orf90. Novel common-variant loci were identified in the NRG1 and SLC45A2 regions, and previously associated loci in the IL18RAP and ATP2B4 regions were associated with novel phenotypes. Transcription factor binding site enrichment identified associations with genes implicated with respiratory and craniofacial traits. Additional analyses identified significantly associated pathways. Conclusions We have identified the first gene-based rare-variant associations with objectively measured sleep-disordered breathing traits. Our results increase the understanding of the genetic architecture of sleep-disordered breathing and highlight associations in genes that modulate lung development, inflammation, respiratory rhythmogenesis, and HIF1A-mediated hypoxic response.

Identification of transcription factor binding sites is one of the most important steps in understanding the function of transcription factors and regulatory networks in organisms. The assay for transposase accessible chromatin sequencing (ATAC-seq) is a simple protocol for detection of open chromatin that could be a powerful tool to advance studies of protein-DNA interactions. Although ATAC-seq has been used in systematic identification of cis-regulatory regions in animal and plant genomes, this method has been rarely applied in fungi. Here, we describe a valuable ATAC-seq resource in the genome of an economically important phytopathogen, the rice false smut fungus Ustilaginoidea virens. The ATAC-seq data of U. virens mycelia collected from potato sucrose broth (PSB) and PSB supplied with rice spikelet extract were both generated. This is the first genome-wide profiling of open chromatin and transcription factor binding sites in U. virens. [Formula: see text] Copyright © 2021 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license .

Motivation Increasing evidence suggests that post-transcriptional ribonucleic acid (RNA) modifications regulate essential biomolecular functions and are related to the pathogenesis of various diseases. Precise identification of RNA modification sites is essential for understanding the regulatory mechanisms of RNAs. To date, many computational approaches for predicting RNA modifications have been developed, most of which were based on strong supervision enabled by base-resolution epitranscriptome data. However, high-resolution data may not be available. Results We propose WeakRM, the first weakly supervised learning framework for predicting RNA modifications from low-resolution epitranscriptome datasets, such as those generated from acRIP-seq and hMeRIP-seq. Evaluations on three independent datasets (corresponding to three different RNA modification types and their respective sequencing technologies) demonstrated the effectiveness of our approach in predicting RNA modifications from low-resolution data. WeakRM outperformed state-of-the-art multi-instance learning methods for genomic sequences, such as WSCNN, which was originally designed for transcription factor binding site prediction. Additionally, our approach captured motifs that are consistent with existing knowledge, and visualization of the predicted modification-containing regions unveiled the potentials of detecting RNA modifications with improved resolution. Availability implementation The source code for the WeakRM algorithm, along with the datasets used, are freely accessible at: https://github.com/daiyun02211/WeakRM Supplementary information Supplementary data are available at Bioinformatics online.

Objective The objective of this study was identification of the transcription factor binding sites (TFBS) in the promoter of HOX genes and elucidation of the comprehensive interaction of transcription factors (TFs)/genes with HOX. Methodology Promoter sequences of HOXA3, HOXA5, HOXA9, HOXA10, HOXA13, HOXB5, HOXC10, HOXC12, and HOXD10 were analyzed to predict the TFBS and their targets using TRANSFAC, TRRUST, and Harmonizome. Functional analysis of the processed data sets was carried out using DAVID and GATHER gene annotation tools. A network of regulatory interactions was constructed using NetworkAnalyst and a comprehensive illustration of the TF-gene network was constructed with HOX as a central hub using the Encyclopedia of DNA Elements chromatin immunoprecipitation sequencing data. Further, the enriched network was constructed to elucidate the roles of these genes in the various pathways. Results Binding sites for E2F1, HNF3α, SP3, and KLF6 were common to promoter regions of all of the HOX genes. The functional annotation and pathway analysis elucidated the regulatory activity of a distinct set of TF-genes in interaction with HOX. A P value ≤.05 and false discovery rate ≤0.01 were considered statistically significant. Conclusion We have confirmed that the predicted TFBSs in the HOX gene promoters function in transcriptional regulation by modulating target gene activity. TF-gene interactions are crucial to understanding oral carcinogenesis.

Single-cell RNA sequencing (scRNA-seq) has emerged as a powerful experimental approach to study cellular heterogeneity. One of the challenges in scRNA-seq data analysis is integrating different types of biological data to consistently recognize discrete biological functions and regulatory mechanisms of cells, such as transcription factor activities and gene regulatory networks in distinct cell populations. We have developed an approach to infer transcription factor activities from scRNA-seq data that leverages existing biological data on transcription factor binding sites. The Bayesian inference transcription factor activity model (BITFAM) integrates ChIP-seq transcription factor binding information into scRNA-seq data analysis. We show that the inferred transcription factor activities for key cell types identify regulatory transcription factors that are known to mechanistically control cell function and cell fate. The BITFAM approach not only identifies biologically meaningful transcription factor activities, but also provides valuable insights into underlying transcription factor regulatory mechanisms.

Targeted engineering of plant gene expression holds great promise for ensuring food security and for producing biopharmaceuticals in plants. However, this engineering requires thorough knowledge of cis-regulatory elements to precisely control either endogenous or introduced genes. To generate this knowledge, we used a massively parallel reporter assay to measure the activity of nearly complete sets of promoters from Arabidopsis, maize and sorghum. We demonstrate that core promoter elements—notably the TATA box—as well as promoter GC content and promoter-proximal transcription factor binding sites influence promoter strength. By performing the experiments in two assay systems, leaves of the dicot tobacco and protoplasts of the monocot maize, we detect species-specific differences in the contributions of GC content and transcription factors to promoter strength. Using these observations, we built computational models to predict promoter strength in both assay systems, allowing us to design highly active promoters comparable in activity to the viral 35S minimal promoter. Our results establish a promising experimental approach to optimize native promoter elements and generate synthetic ones with desirable features. A massively parallel reporter assay was used to measure the activity of nearly complete sets of promoters from Arabidopsis, maize and sorghum in two assay systems, uncovering the sequence features affecting promoter strength and facilitating promoter strength prediction and synthetic design.

Transcription factors (TFs) bind DNA in a sequence-specific manner and thereby regulate target gene expression. TF binding and its regulatory activity is highly context dependent, and is not only determined by specific cell types or differentiation stages but also relies on other regulatory mechanisms, such as DNA and chromatin modifications. Interactions between TFs and their DNA binding sites are critical mediators of phenotypic variation and play important roles in the onset of disease. A continuously growing number of studies therefore attempts to elucidate TF:DNA interactions to gain knowledge about regulatory mechanisms and disease-causing variants. Here we summarize how TF-binding characteristics and the impact of variants can be investigated, how bioinformatic tools can be used to analyze and predict TF:DNA binding, and what additional information can be obtained from the TF protein structure.

MotivationAberrant DNA methylation in transcription factor binding sites has been shown to lead to anomalous gene regulation that is strongly associated with human disease. However, the majority of methylation-sensitive positions within transcription factor binding sites remain unknown. Here we introduce SEMplMe, a computational tool to generate predictions of the effect of methylation on transcription factor binding strength in every position within a transcription factor’s motif.ResultsSEMplMe uses ChIP-seq and whole genome bisulfite sequencing to predict effects of methylation within binding sites. SEMplMe validates known methylation sensitive and insensitive positions within a binding motif, identifies cell type specific transcription factor binding driven by methylation, and outperforms SELEX-based predictions for CTCF. These predictions can be used to identify aberrant sites of DNA methylation contributing to human disease.Availability and ImplementationSEMplMe is available fromhttps://github.com/Boyle-Lab/SEMplMe.Contactapboyle@umich.eduSupplementary InformationSupplementary data are available atBioinformaticsonline.

Genome sequencing of tens of thousands of human individuals has recently enabled the measurement of large selective effects for mutations to protein-coding genes. Here we describe a new method, called ExtRaINSIGHT, for measuring similar selective effects at individual sites in noncoding as well as in coding regions of the human genome. ExtRaINSIGHT estimates the prevalance of strong purifying selection, or "ultraselection" ({lambda}s), as the fractional depletion of rare single-nucleotide variants (minor allele frequency < 0.1%) in a target set of genomic sites relative to matched sites that are putatively neutrally evolving, in a manner that controls for local variation and neighbor-dependence in mutation rate. We show using simulations that, above an appropriate threshold,{lambda} s is closely related to the average site-specific selection coefficient against heterozygous point mutations, as predicted at mutation-selection balance. Applying ExtRaINSIGHT to 71,702 whole genome sequences from gnomAD v3, we find particularly strong evidence of ultraselection in evolutionarily ancient miRNAs and neuronal protein-coding genes, as well as at splice sites. Moreover, our estimated selection coefficient against heterozygous amino-acid replacements across the genome (at 1.4%) is substantially larger than previous estimates based on smaller sample sizes. By contrast, we find weak evidence of ultraselection in other noncoding RNAs and transcription factor binding sites, and only modest evidence in ultraconserved elements and human accelerated regions. We estimate that [~]0.3-0.5% of the human genome is ultraselected, with one third to one half of ultraselected sites falling in coding regions. These estimates suggest [~]0.3-0.4 lethal or nearly lethal de novo mutations per potential human zygote, together with [~]2 de novo mutations that are more weakly deleterious. Overall, our study sheds new light on the genome-wide distribution of fitness effects for new point mutations by combining deep new sequencing data sets and classical theory from population genetics.

Quorum sensing gene expression in vibrios is regulated by the LuxR/HapR family of transcriptional factors, which includes Vibrio vulnificus SmcR. The consensus binding site of Vibrio LuxR/HapR/SmcR proteins is palindromic but highly degenerate with sequence variations at each promoter. To examine the mechanism by which SmcR recognizes diverse DNA sites, we generated SmcR separation-of-function mutants that either repress or activate transcription but not both. SmcR N55I is restricted in recognition of single base-pair variations in DNA binding site sequences and thus is defective at transcription activation but retains interaction with RNA polymerase (RNAP) alpha. SmcR S76A, L139R and N142D substitutions disrupt the interaction with RNAP alpha but retain functional DNA binding activity. X-ray crystallography and small angle X-ray scattering data show that the SmcR DNA binding domain exists in two conformations (wide and narrow), and the protein complex forms a mixture of dimers and tetramers in solution. The three RNAP interaction-deficient variants also have two DNA binding domain conformations, whereas SmcR N55I exhibits only the wide conformation. These data support a model in which two mechanisms drive SmcR transcriptional activation: interaction with RNAP and a multi-conformational DNA binding domain that permits recognition of variable DNA sites.

Chromatin accessibility of a promoter is fundamental in regulating transcriptional activity. The histone variant H2A.Z has been shown to contribute to this regulation, but its role has remained poorly understood. Here, we prepare high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for all human Pol II promoters in epithelial, mesenchymal and isogenic cancer cell lines. We find that, in contrast to the prevailing model, many different types of active and inactive promoter structures are observed that differ in their nucleosome organization and sensitivity to MNase digestion. Key aspects of an active chromatin structure include positioned H2A.Z MNase resistant nucleosomes upstream or downstream of the TSS, and a MNase sensitive nucleosome at the TSS. Furthermore, the loss of H2A.Z leads to a dramatic increase in the accessibility of transcription factor binding sites. Collectively, these results suggest that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter., Histone variant H2A.Z has been suggested to contribute to the regulation of promoter accessibility. Here, the authors present high-depth maps of the position and accessibility of H2A.Z-containing nucleosomes for human Pol II promoters and provide evidence that H2A.Z has multiple and distinct roles in regulating gene expression dependent upon its location in a promoter.

Enhancers represent non-coding regulatory regions of the genome located distantly from their target genes. They regulate gene expression programs in a context-specific manner via interacting with promoters of one or more target genes and are generally associated with transcription factor binding sites and epi(genomic)/chromatin features, such as regions of chromatin accessibility and histone modifications. The enhancers are difficult to identify due to the modularity of their associated features. Although enhancers have been studied extensively in human and animals, only a handful of them has been identified in few plant species till date due to non-availability of plant-specific experimental and computational approaches for enhancer identification. Being an important regulatory component of the genome, enhancers represent potential targets for engineering agronomic traits, including salinity stress tolerance in plants. Here, we provide a review of the available experimental and computational approaches along with the associated sequence and chromatin/epigenetic features for the discovery of enhancers in plants. In addition, we provide insights into the challenges and future prospects of enhancer research in plant biology with emphasis on potential applications in engineering salinity stress tolerance in crop plants.

The glucocorticoid (GR) and androgen (AR) receptors execute unique functions in vivo, yet have nearly identical DNA binding specificities. To identify mechanisms that facilitate functional diversification among these transcription factor paralogs, we studied AR and GR in an equivalent cellular context. Analysis of chromatin and sequence features suggest that divergent binding, and corresponding gene regulation, are driven by different abilities of AR and GR to interact with relatively inaccessible chromatin. Divergent genomic binding patterns can also be the results of subtle differences in DNA binding preference between AR and GR. Furthermore, the sequence composition of large regions (>10 kb) surrounding selectively occupied binding sites differs significantly, indicating a role for the sequence environment in selectively guiding AR and GR to distinct binding sites. The comparison of binding sites that are shared between AR and GR shows that the specificity paradox can also be resolved by differences in the events that occur downstream of receptor binding. Specifically, we find that shared binding sites display receptor-specific enhancer activity, cofactor recruitment and changes in histone modifications. Genomic deletion of shared binding sites demonstrates their contribution to directing receptor-specific gene regulation. Together, these data suggest that differences in genomic occupancy as well as divergence in the events that occur downstream of receptor binding direct functional diversification among transcription factor paralogs.

Background: Overlap of pathways enriched by single nucleotide polymorphisms and DNA-methylation underlying chronic postsurgical pain (CPSP), prompted pilot study of CPSP-associated methylation quantitative trait loci (meQTL). Materials & methods: Children undergoing spine-fusion were recruited prospectively. Logistic-regression for genome- and epigenome-wide CPSP association and DNA-methylation-single nucleotide polymorphism association/mediation analyses to identify meQTLs were followed by functional genomics analyses. Results: CPSP (n = 20/58) and non-CPSP groups differed in pain-measures. Of 2753 meQTLs, DNA-methylation at 127 cytosine–guanine dinucleotides mediated association of 470 meQTLs with CPSP (p Clinical trial registration: NCT01839461 , NCT01731873 (ClinicalTrials.gov).

Enhancers located in introns are abundant and play a major role in the regulation of gene expression in mammalian species. By contrast, the functions of intronic enhancers in plants have largely been unexplored and only a handful of plant intronic enhancers have been reported. We performed a genome-wide prediction of intronic enhancers in Arabidopsis thaliana using open chromatin signatures based on DNase I sequencing. We identified 941 candidate intronic enhancers associated with 806 genes in seedling tissue and 1,271 intronic enhancers associated with 1,069 genes in floral tissue. We validated the function of 15 of 21 (71%) of the predicted intronic enhancers in transgenic assays using a reporter gene. We also created deletion lines of three intronic enhancers associated with two different genes using CRISPR/Cas. Deletion of these enhancers, which span key transcription factor binding sites, did not abolish gene expression but caused varying levels of transcriptional repression of their cognate genes. Remarkably, the transcriptional repression of the deletion lines occurred at specific developmental stages and resulted in distinct phenotypic effects on plant morphology and development. Clearly, these three intronic enhancers are important in fine-tuning tissue- and development-specific expression of their cognate genes.

The most popular model for the search of ChIP-seq data for transcription factor binding sites (TFBS) is the positional weight matrix (PWM). However, this model does not take into account dependencies between nucleotide occurrences in different site positions. Currently, two recently proposed models, BaMM and InMoDe, can do as much. However, application of these models was usually limited only to comparing their recognition accuracies with that of PWMs, while none of the analyses of the co-prediction and relative positioning of hits of different models in peaks has yet been performed. To close this gap, we propose the pipeline called MultiDeNA. This pipeline includes stages of model training, assessing their recognition accuracy, scanning ChIP-seq peaks and their classification based on scan results. We applied our pipeline to 22 ChIP-seq datasets of TF FOXA2 and considered PWM, dinucleotide PWM (diPWM), BaMM and InMoDe models. The combination of these four models allowed a significant increase in the fraction of recognized peaks compared to that for the sole PWM model: the increase was 26.3 %. The BaMM model provided the main contribution to the recognition of sites. Although the major fraction of predicted peaks contained TFBS of different models with coincided positions, the medians of the fraction of peaks containing the predictions of sole models were 1.08, 0.49, 4.15 and 1.73 % for PWM, diPWM, BaMM and InMoDe, respectively. Thus, FOXA2 BSs were not fully described by only a sole model, which indicates theirs heterogeneity. We assume that the BaMM model is the most successful in describing the structure of the FOXA2 BS in ChIP-seq datasets under study.. В настоящее время самой распространенной моделью поиска сайтов связывания транскрипционных факторов (ССТФ) в пиках ChIP-seq является позиционная весовая матрица (position weight matrix, PWM). Но эта модель не учитывает взаимосвязи между частотами встреч нуклеотидов в разных позициях ССТФ, поэтому не способна гарантировать определение всех возможных структурных вариантов ССТФ. На сегодняшний день уже предложены альтернативные модели, например BaMM и InMoDe, которые учитывают такие взаимосвязи. Однако применение этих моделей обычно сводилось к сравнению их точности с точностью традиционной модели PWM, тогда как анализ совместной встречаемости и относительного расположения ССТФ разных моделей в пиках не производился. В нашей работе мы предлагаем конвейер программ MultiDeNA, позволяющий сочетать разные модели de novo поиска ССТФ для выявления структурной гетерогенности ССТФ в данных ChIP-seq. Разработанный конвейер включает этапы построения моделей на основе заданного набора пиков, оценки точности распознавания моделей с помощью перекрестных тестов, выбора порогов, сканирования пиков ChIP-seq и классификацию пиков по результатам сканирования. С применением конвейера нами проведен анализ 22 экспериментов ChIP-seq для ТФ FOXA2 с помощью четырех моделей: PWM, diPWM, BaMM и InMoDe. Показано, что сочетание моделей позволяет существенно увеличить общее количество распознанных пиков (на 26.3 %) по сравнению с применением только PWM; при этом основной вклад в распознавание внесла модель BaMM. В значительной доле пиков разные модели распознают совпадающие ССТФ; однако для моделей PWM, diPWM, BaMM и InMoDe медианы доли пиков, которые содержали ССТФ только одной модели, составили 1.08, 0.49, 4.15 и 1.73 % соответственно. Таким образом, совокупность ССТФ FOXA2 не описывается полностью только одной моделью, что свидетельствует о наличии структурной гетерогенности в ССТФ у FOXA2.

Transposable element (TE) invasions have shaped vertebrate genomes over the course of evolution. They have contributed an extra layer of species-specific gene regulation by providing novel transcription factor binding sites. In humans, SINE-VNTR-Alu (SVA) elements are one of three still active TE families; approximately 2800 SVA insertions exist in the human genome, half of which are human-specific. TEs are often silenced by KRAB zinc finger (KZNF) proteins recruiting corepressor proteins that establish a repressive chromatin state. A number of KZNFs have been reported to bind SVAs, but their individual contribution to repressing SVAs and their roles in suppressing SVA-mediated gene-regulatory effects remains elusive. We analyzed the genome-wide binding profile for ZNF91 in human cells and found that ZNF91 interacts with the VNTR region of SVAs. Through CRISPR-Cas9-mediated deletion of ZNF91 in human embryonic stem cells, we established that loss of ZNF91 results in increased transcriptional activity of SVAs. In contrast, SVA activation was not observed upon genetic deletion of the ZNF611 gene encoding another strong SVA interactor. Epigenetic profiling confirmed the loss of SVA repression in the absence of ZNF91 and revealed that mainly evolutionary young SVAs gain gene activation-associated epigenetic modifications. Genes close to activated SVAs showed a mild up-regulation, indicating SVAs adopt properties of cis-regulatory elements in the absence of repression. Notably, genome-wide derepression of SVAs elicited the communal up-regulation of KZNFs that reside in KZNF clusters. This phenomenon may provide new insights into the potential mechanisms used by the host genome to sense and counteract TE invasions.

TEOSINTE BRANCHED 1/CYCLOIDEA/PROLIFERATING CELL FACTOR 1 (TCP) transcription factors control multiple aspects of growth and development in various plant species. However, few genes were reported to be directly targeted and regulated by them through their specific binding sites, and then uncover their functions in plants. A consensus DNA-binding site motif of TCP2 was identified by random binding site selection (RBSS). DNA recognized by TCP2 contained the motif G(G/T)GGNCC(A/C), which showed high consistency with motifs bound by other TCP domain proteins. Consequently, this motif was regarded as the specific DNA-binding sites of TCP2. Circadian clock associated 1 (CCA1) and EARLY FLOWERING 3 (ELF3) were subsequently considered as potential target genes owing to the containing of the similar TCP2 binding sites or core binding sites GGNCC and found to be positively regulated by TCP2 via DNA binding. Phenotype analysis results showed that mutation and over-expression of TCP2 resulted in variations in leaf morphogenesis, especially the double or triple mutations of TCP2, 4 and 10. Mutations in TCPs caused late flowering. Finally, TCP2 was shown to influence hypocotyl elongation by mediating the jasmonate signaling pathway. Overall, these results provide a basis for future studies aimed at distinguishing the target genes of TCP2 and elucidating the important roles of TCP2 in plant growth and development.

The zinc finger-containing transcription factor Gli3 is a key mediator of Hedgehog (Hh) signaling pathway. In vertebrates, Gli3 has widespread expression pattern during early embryonic development. Along the anteroposterior axes of the central nervous system (CNS), dorsoventral neural pattern elaboration is achieved through Hh mediated spatio-temporal deployment of Gli3 transcripts. Previously, we and others uncovered a set of enhancers that mediate many of the known aspects of Gli3 expression during neurogenesis. However, the potential role of Gli3 associated enhancers in trait evolution has not yet received any significant attention. Here, we investigate the evolutionary patterns of Gli3 associated CNS-specific enhancers that have been reported so far. A subset of these enhancers has undergone an accelerated rate of molecular evolution in the human lineage in comparison to other primates/mammals. These fast-evolving enhancers have acquired human-specific changes in transcription factor binding sites (TFBSs). These human-unique changes within subset of Gli3 associated CNS-specific enhancers were further validated as single nucleotide polymorphisms through 1000 Genome Project Phase 3 data. This work not only infers the molecular evolutionary patterns of Gli3 associated enhancers but also provides clues for putative genetic basis of the population-specificity of gene expression regulation.

Although Saccharomyces cerevisiae and S. cerevisiae var. boulardii share more than 95% genome sequence homology, only S. cerevisiae var. boulardii displays probiotic activity. In this study, the transcriptomic differences exhibited by S. cerevisiae and S. cerevisiae var. boulardii in intestinal like medium were evaluated. S. cerevisiae was found to display stress response overexpression, consistent with higher ability of S. cerevisiae var. boulardii to survive within the human host, while S. cerevisiae var. boulardii exhibited transcriptional patterns associated with probiotic activity, suggesting increased acetate biosynthesis. Resorting to the creation of a S. cerevisiae var. boulardii genomic database within Yeastract+, a possible correlation between loss or gain of transcription factor binding sites in S. cerevisiae var. boulardii promoters and the transcriptomic pattern is discussed. This study suggests that S. cerevisiae var. boulardii probiotic activity, when compared to S. cerevisiae, relies, at least partially, on differential expression regulation, based on promoter variability.

Whole genome and whole exome sequencing technologies play a very important role in the studies of the genetic aspects of the pathogenesis of various diseases. The ample use of genome-wide and exome-wide association study methodology (GWAS and EWAS) made it possible to identify a large number of genetic variants associated with diseases. This information is accumulated in the databases like GWAS central, GWAS catalog, OMIM, ClinVar, etc. Most of the variants identified by the GWAS technique are located in the noncoding regions of the human genome. According to the ENCODE project, the fraction of regions in the human genome potentially involved in transcriptional control is many times greater than the fraction of coding regions. Thus, genetic variation in noncoding regions of the genome can increase the susceptibility to diseases by disrupting various regulatory elements (promoters, enhancers, silencers, insulator regions, etc.). However, identification of the mechanisms of influence of pathogenic genetic variants on the diseases risk is difficult due to a wide variety of regulatory elements. The present review focuses on the molecular genetic mechanisms by which pathogenic genetic variants affect gene expression. At the same time, attention is concentrated on the transcriptional level of regulation as an initial step in the expression of any gene. A triggering event mediating the effect of a pathogenic genetic variant on the level of gene expression can be, for example, a change in the functional activity of transcription factor binding sites (TFBSs) or DNA methylation change, which, in turn, affects the functional activity of promoters or enhancers. Dissecting the regulatory roles of polymorphic loci have been impossible without close integration of modern experimental approaches with computer analysis of a growing wealth of genetic and biological data obtained using omics technologies. The review provides a brief description of a number of the most well-known public genomic information resources containing data obtained using omics technologies, including (1) resources that accumulate data on the chromatin states and the regions of transcription factor binding derived from ChIP-seq experiments; (2) resources containing data on genomic loci, for which allele-specific transcription factor binding was revealed based on ChIP-seq technology; (3) resources containing in silico predicted data on the potential impact of genetic variants on the transcription factor binding sites.Полногеномные и полноэкзомные технологии секвенирования играют важную роль в исследованиях генетических аспектов патогенеза различных заболеваний. Широкое применение методов полногеномного и полноэкзомного анализа ассоциаций позволяет идентифицировать множество вариантов геномной изменчивости (ГИ), ассоциированных с заболеваниями. Эта информация накапливается в базах данных GWAS central, GWAS catalog, OMIM, ClinVar и др. Большинство вариантов, идентифицированных методикой полногеномного анализа ассоциаций, располагается в некодирующих областях генома человека. По данным проекта ENCODE, доля участков в геноме человека, потенциально задействованных в регуляции транскрипции, во много раз превышает долю кодирующих областей. Таким образом, геномная изменчивость в некодирующих областях генома может повышать предрасположенность к заболеваниям, нарушая функционирование различных регуляторных элементов (промоторов, энхансеров, участков, определяющих 3D структуру хроматина и т. д.). Однако идентификация механизмов влияния патогенных вариантов ГИ на риск развития заболеваний затруднена ввиду большого разнообразия регуляторных элементов. В обзоре рассмотрены молекулярно-генетические механизмы влияния патогенных вариантов ГИ на экспрессию генов. При этом внимание сосредоточено на транскрипционном уровне регуляции как ключевой стадии, запускающей последовательность этапов экспрессии любого гена. Пусковым событием, опосредующим влияние патогенного варианта ГИ на уровень экспрессии гена, может быть, например, изменение функциональной активности сайтов связывания транскрипционных факторов или уровня метилирования ДНК, что, в свою очередь, отражается на функциональной активности промоторов или энхансеров. Выявление регуляторных эффектов полиморфных локусов невозможно без тесной интеграции современных экспериментальных подходов с компьютерным анализом больших массивов генетических данных, получаемых на основе омиксных технологий. В обзоре кратко описаны наиболее известные открытые полногеномные информационные ресурсы, содержащие данные, полученные на основе омиксных технологий, в том числе: ресурсы, накапливающие сведения о состоянии хроматина и участках его связывания с транскрипционными факторами, выявленными с помощью технологии ChIP-seq; ресурсы по геномным локусам, для которых на основе данных ChIP-seq выявлено аллель-специфичное связывание с транскрипционными факторами; а также ресурсы, содержащие предсказанные in silico данные о потенциальном влиянии геномной изменчивости на сайты связывания транскрипционных факторов.

Perturbation and mechanistic studies have shown that the DNA-binding transcription factor Zeb2 controls cell fate decision, differentiation and/or maturation in multiple cell lineages in embryos and after birth. In cultured embryonic stem cells (ESCs) Zeb2s strong upregulation is necessary for their exit from primed pluripotency and entering neural and general differentiation. We engineered mouse ESCs to produce epitope-tagged Zeb2 from one of its two endogenous alleles. Using crosslinking ChIP-sequencing, we mapped for the first time 2,432 DNA-binding sites of Zeb2 in ESC-derived neuroprogenitor cells (NPCs). A new, major site maps promoter-proximal to Zeb2 itself, and its homozygous removal demonstrates that Zeb2 autoregulation is necessary to elicit proper Zeb2 effects in ESC[->]NPC differentiation. We then cross-referenced all Zeb2 DNA-binding sites with transcriptome data from Zeb2 perturbations in ESCs, ventral forebrain in mouse embryos, and adult neurogenesis. While the characteristics of these neurodevelopmental systems differ, we find interesting overlaps. Collectively, these new results obtained in ESC-derived NPCs significantly add to Zeb2s role as neurodevelopmental regulator as well as the causal gene in Mowat-Wilson Syndrome. Also, Zeb2 was found to map to loci mutated in other human congenital syndromes, making variant or disturbed levels of ZEB2 a candidate modifier principle in these. Significance statementZeb2 is needed for exit from primed pluripotency and differentiation of ESCs and being studied in various cell lineages and tissues/organs in knockouts and adult mice. Phenotyping has mainly documented transcriptomes of Zeb2-mutant affected cell types, but ChIP-seq data are still not available to document Zeb2s direct target loci as repressor or activator of transcription. This study maps for the first time the Zeb2 genome-wide binding sites in ESC-derived neuroprogenitors, and discovered Zeb2 autoregulation. Cross-referencing these ChIP-seq data with transcriptome profiles of Zeb2 perturbation further explains Mowat-Wilson Syndrome defects and proposes ZEB2 as candidate modifier gene in other syndromes.

Aims: Robust and more accurate method for identifying transcription factor binding sites (TFBS) for gene expression. Background: Deep neural networks (DNNs) have shown promising growth in solving complex machine learning problems. Conventional techniques are comfortably replaced by DNNs in computer vision, signal processing, healthcare, and genomics. Understanding DNA sequences is always a crucial task in healthcare and regulatory genomics. For DNA motif prediction, choosing the right dataset with a sufficient number of input sequences is crucial in order to design an effective model. Objective: Designing a new algorithm which works on different dataset while an improved performance for TFBS prediction. Methods: With the help of Layerwise Relevance Propagation, the proposed algorithm identifies the invariant features with adaptive noise patterns. Results: The performance is compared by calculating various metrics on standard as well as recent methods and significant improvement is noted. Conclusion: By identifying the invariant and robust features in the DNA sequences, the classification performance can be increased.

Tilapia is a hardy fish which can survive in water bodies polluted with heavy metals. Metal resistance is conferred by higher expression of metallothionein gene (mt) in many organisms. Level, time and tissue-specificity of gene expression is regulated through transcription factor binding sites (TFBS) which may be present in the upstream, downstream, or even in the introns of a gene. So, as a candidate regulatory region, the 5’upstream sequence of mt gene in three tilapia species, Oreochromis aureus, O. niloticus and O. mossambicus was studied. The targeted region was PCR-amplified and then sequenced using a pair of custom-designed primer. A total of only 2.7% variation was found in the sequenced genomic region among the three species. Metal-related TFBS were predicted from these sequences. A total of twenty eight TFBS were found in O. aureus and twenty nine in O. mossambicus and O. niloticus. The number of metalrelated TFBS predicted in the targeted sequence was significantly higher compared to that found in randomly selected other genomic regions of same size from O. niloticus genome. Thus, the results suggest the presence of putative regulatory elements in the targeted upstream region which might have important role in the regulation of mt gene function. Dhaka Univ. J. Biol. Sci. 30(1): 95-103, 2021 (January)

Objective: Vanin1 (VNN1) is a pantetheinase that can catalyze the hydrolysis of pantetheine to produce pantothenic acid and cysteamine. Our previous studies showed that VNN1 is specifically expressed in chicken liver. In this study, we aimed to investigate the roles of peroxisome proliferators activated receptor α (PPARα) and miRNA-181a-5p in regulating VNN1 gene expression in chicken liver.Methods: 5′-RACE was performed to identify the transcription start site of chicken VNN1. JASPAR and TFSEARCH were used to analyze the potential transcription factor binding sites in the promoter region of chicken VNN1 and miRanda was used to search miRNA binding sites in 3′ untranslated region (3′UTR) of chicken VNN1. We used a knock-down strategy to manipulate PPARα (or miRNA-181a-5p) expression levels in vitro to further investigate its effect on VNN1 gene transcription. Luciferase reporter assays were used to explore the specific regions of VNN1 targeted by PPARα and miRNA-181a-5p.Results: Sequence analysis of the VNN1 promoter region revealed several transcription factor-binding sites, including hepatocyte nuclear factor 1α (HNF1α), PPARα, and CCAAT/enhancer binding protein α. GW7647 (a specific agonist of PPARα) increased the expression level of VNN1 mRNA in chicken primary hepatocytes, whereas knockdown of PPARα with siRNA increased VNN1 mRNA expression. Moreover, the predicted PPARα-binding site was confirmed to be necessary for PPARα regulation of VNN1 gene expression. In addition, the VNN1 3′UTR contains a sequence that is completely complementary to nucleotides 1 to 7 of miRNA-181a-5p. Overexpression of miR-181a-5p significantly decreased the expression level of VNN1 mRNA.Conclusion: This study demonstrates that PPARα is an important transcriptional activator of VNN1 gene expression and that miRNA-181a-5p acts as a negative regulator of VNN1 expression in chicken hepatocytes.