Structure of the human SAGA coactivator complex

The SAGA complex is a regulatory hub involved in gene regulation, chromatin modification, DNA damage repair and signaling. While structures of yeast SAGA (ySAGA) have been reported, there are noteworthy functional and compositional differences for this complex in metazoans. Here we present the cryogenic-electron microscopy (cryo-EM) structure of human SAGA (hSAGA) and show how the arrangement of distinct structural elements results in a globally divergent organization from that of yeast, with a different interface tethering the core module to the TRRAP subunit, resulting in a dramatically altered geometry of functional elements and with the integration of a metazoan-specific splicing module. Our hSAGA structure reveals the presence of an inositol hexakisphosphate (InsP6) binding site in TRRAP and an unusual property of its pseudo-(Ψ)PIKK. Finally, we map human disease mutations, thus providing the needed framework for structure-guided drug design of this important therapeutic target for human developmental diseases and cancer.

Cis-regulatory architecture of human ESC-derived hypothalamic neuron differentiation aids in variant-to-gene mapping of relevant complex traits

The hypothalamus regulates metabolic homeostasis by influencing behavior and endocrine systems. Given its role governing key traits, such as body weight and reproductive timing, understanding the genetic regulation of hypothalamic development and function could yield insights into disease pathogenesis. However, given its inaccessibility, studying human hypothalamic gene regulation has proven challenging. To address this gap, we generate a high-resolution chromatin architecture atlas of an established embryonic stem cell derived hypothalamic-like neuron model across three stages of in vitro differentiation. We profile accessible chromatin and identify physical contacts between gene promoters and putative cis-regulatory elements to characterize global regulatory landscape changes during hypothalamic differentiation. Next, we integrate these data with GWAS loci for various complex traits, identifying multiple candidate effector genes. Our results reveal common target genes for these traits, potentially affecting core developmental pathways. Our atlas will enable future efforts to determine hypothalamic mechanisms influencing disease susceptibility.

The Role of Long Non-Coding RNAs in Endometriosis

Endometriosis is a chronic gynecological disorder affecting the quality of life and fertility of many women around the world. Heterogeneous and non-specific symptoms may lead to a delay in diagnosis, with treatment options limited to surgery and hormonal therapy. Hence, there is a need to better understand the pathogenesis of the disease to improve diagnosis and treatment. Long non-coding RNAs (lncRNAs) have been increasingly shown to be involved in gene regulation but remain relatively under investigated in endometriosis. Mutational and transcriptomic studies have implicated lncRNAs in the pathogenesis of endometriosis. Single-nucleotide polymorphisms (SNPs) in lncRNAs or their regulatory regions have been associated with endometriosis. Genome-wide transcriptomic studies have identified lncRNAs that show deregulated expression in endometriosis, some of which have been subjected to further experiments, which support a role in endometriosis. Mechanistic studies indicate that lncRNAs may regulate genes involved in endometriosis by acting as a molecular sponge for miRNAs, by directly targeting regulatory elements via interactions with chromatin or transcription factors or by affecting signaling pathways. Future studies should concentrate on determining the role of uncharacterized lncRNAs revealed by endometriosis transcriptome studies and the relevance of lncRNAs implicated in the disease by in vitro and animal model studies.

Enhancer RNA m6A methylation facilitates transcriptional condensate formation and gene activation

The mechanistic understanding of nascent RNAs in transcriptional control remains limited. Here, by a high sensitivity method methylation-inscribed nascent transcripts sequencing (MINT-seq), we characterized the landscapes of N6-methyladenosine (m6A) on nascent RNAs. We uncover heavy but selective m6A deposition on nascent RNAs produced by transcription regulatory elements, including promoter upstream antisense RNAs and enhancer RNAs (eRNAs), which positively correlates with their length, inclusion of m6A motif, and RNA abundances. m6A-eRNAs mark highly active enhancers, where they recruit nuclear m6A reader YTHDC1 to phase separate into liquid-like condensates, in a manner dependent on its C terminus intrinsically disordered region and arginine residues. The m6A-eRNA/YTHDC1 condensate co-mixes with and facilitates the formation of BRD4 coactivator condensate. Consequently, YTHDC1 depletion diminished BRD4 condensate and its recruitment to enhancers, resulting in inhibited enhancer and gene activation. We propose that chemical modifications of eRNAs together with reader proteins play broad roles in enhancer activation and gene transcriptional control.

Cell-type-specific effects of genetic variation on chromatin accessibility during human neuronal differentiation

Common genetic risk for neuropsychiatric disorders is enriched in regulatory elements active during cortical neurogenesis. However, it remains poorly understood as to how these variants influence gene regulation. To model the functional impact of common genetic variation on the noncoding genome during human cortical development, we performed the assay for transposase accessible chromatin using sequencing (ATAC-seq) and analyzed chromatin accessibility quantitative trait loci (QTL) in cultured human neural progenitor cells and their differentiated neuronal progeny from 87 donors. We identified significant genetic effects on 988/1,839 neuron/progenitor regulatory elements, with highly cell-type and temporally specific effects. A subset (roughly 30%) of chromatin accessibility-QTL were also associated with changes in gene expression. Motif-disrupting alleles of transcriptional activators generally led to decreases in chromatin accessibility, whereas motif-disrupting alleles of repressors led to increases in chromatin accessibility. By integrating cell-type-specific chromatin accessibility-QTL and brain-relevant genome-wide association data, we were able to fine-map and identify regulatory mechanisms underlying noncoding neuropsychiatric disorder risk loci.

Cis-regulatory evolution spotlights species differences in the adaptive potential of gene expression plasticity

Phenotypic plasticity is the variation in phenotype that a single genotype can produce in different environments and, as such, is an important component of individual fitness. However, whether the effect of new mutations, and hence evolution, depends on the direction of plasticity remains controversial. Here, we identify the cis-acting modifications that have reshaped gene expression in response to dehydration stress in three Arabidopsis species. Our study shows that the direction of effects of most cis-regulatory variants differentiating the response between A. thaliana and the sister species A. lyrata and A. halleri depends on the direction of pre-existing plasticity in gene expression. A comparison of the rate of cis-acting variant accumulation in each lineage indicates that the selective forces driving adaptive evolution in gene expression favors regulatory changes that magnify the stress response in A. lyrata. The evolutionary constraints measured on the amino-acid sequence of these genes support this interpretation. In contrast, regulatory changes that mitigate the plastic response to stress evolved more frequently in A. halleri. Our results demonstrate that pre-existing plasticity may be a stepping stone for adaptation, but its selective remodeling differs between lineages.

A cis-regulatory atlas in maize at single-cell resolution

cis-regulatory elements (CREs) encode the genomic blueprints of spatiotemporal gene expression programs enabling highly specialized cell functions. Using single-cell genomics in six maize organs, we determined the cis- and trans-regulatory factors defining diverse cell identities and coordinating chromatin organization by profiling transcription factor (TF) combinatorics, identifying TFs with non-cell-autonomous activity, and uncovering TFs underlying higher-order chromatin interactions. Cell-type-specific CREs were enriched for enhancer activity and within unmethylated long terminal repeat retrotransposons. Moreover, we found cell-type-specific CREs are hotspots for phenotype-associated genetic variants and were targeted by selection during modern maize breeding, highlighting the biological implications of this CRE atlas. Through comparison of maize and Arabidopsis thaliana developmental trajectories, we identified TFs and CREs with conserved and divergent chromatin dynamics, showcasing extensive evolution of gene regulatory networks. In addition to this rich dataset, we developed single-cell analysis software, Socrates, which can be used to understand cis-regulatory variation in any species.

Pioneering the developmental frontier

Coordinated changes in gene expression allow a single fertilized oocyte to develop into a complex multi-cellular organism. These changes in expression are controlled by transcription factors that gain access to discrete cis-regulatory elements in the genome, allowing them to activate gene expression. Although nucleosomes present barriers to transcription factor occupancy, pioneer transcription factors have unique properties that allow them to bind DNA in the context of nucleosomes, define cis-regulatory elements, and facilitate the subsequent binding of additional factors that determine gene expression. In this capacity, pioneer factors act at the top of gene-regulatory networks to control developmental transitions. Developmental context also influences pioneer factor binding and activity. Here we discuss the interplay between pioneer factors and development, their role in driving developmental transitions, and the influence of the cellular environment on pioneer factor binding and activity.

Cloning and functional characterization of seed-specific LEC1A promoter from peanut (Arachis hypogaea L.)

LEAFY COTYLEDON1 (LEC1) is a HAP3 subunit of CCAAT-binding transcription factor, which controls several aspects of embryo and postembryo development, including embryo morphogenesis, storage reserve accumulation and skotomorphogenesis. Herein, using the method of chromosomal walking, a 2707bp upstream sequence from the ATG initiation codon site of AhLEC1A which is a homolog of Arabidopsis LEC1 was isolated in peanut. Its transcriptional start site confirmed by 5’ RACE was located at 82 nt from 5’ upstream of ATG. The bioinformatics analysis revealed that there existed many tissue-specific elements and light responsive motifs in its promoter. To identify the functional region of the AhLEC1A promoter, seven plant expression vectors expressing the GUS (β-glucuronidase) gene, driven by 5’ terminal series deleted fragments of AhLEC1A promoter, were constructed and transformed into Arabidopsis. Results of GUS histochemical staining showed that the regulatory region containing 82bp of 5’ UTR and 2228bp promoter could facilitate GUS to express preferentially in the embryos at different development periods of Arabidopsis. Taken together, it was inferred that the expression of AhLEC1A during seed development of peanut might be controlled positively by several seed-specific regulatory elements, as well as negatively by some other regulatory elements inhibiting its expression in other organs. Moreover, the GUS expression pattern of transgenic seedlings in darkness and in light was relevant to the light-responsive elements scattered in AhLEC1A promoter segment, implying that these light-responsive elements harbored in the AhLEC1A promoter regulate skotomorphogenesis of peanut seeds, and AhLEC1A expression was inhibited after the germinated seedlings were transferred from darkness to light.

CRISPRi screens reveal a DNA methylation-mediated 3D genome dependent causal mechanism in prostate cancer

Prostate cancer (PCa) risk-associated SNPs are enriched in noncoding cis-regulatory elements (rCREs), yet their modi operandi and clinical impact remain elusive. Here, we perform CRISPRi screens of 260 rCREs in PCa cell lines. We find that rCREs harboring high risk SNPs are more essential for cell proliferation and H3K27ac occupancy is a strong indicator of essentiality. We also show that cell-line-specific essential rCREs are enriched in the 8q24.21 region, with the rs11986220-containing rCRE regulating MYC and PVT1 expression, cell proliferation and tumorigenesis in a cell-line-specific manner, depending on DNA methylation-orchestrated occupancy of a CTCF binding site in between this rCRE and the MYC promoter. We demonstrate that CTCF deposition at this site as measured by DNA methylation level is highly variable in prostate specimens, and observe the MYC eQTL in the 8q24.21 locus in individuals with low CTCF binding. Together our findings highlight a causal mechanism synergistically driven by a risk SNP and DNA methylation-mediated 3D genome architecture, advocating for the integration of genetics and epigenetics in assessing risks conferred by genetic predispositions.

Detection of gene cis-regulatory element perturbations in single-cell transcriptomes

We introduce poly-adenine CRISPR gRNA-based single-cell RNA-sequencing (pAC-Seq), a method that enables the direct observation of guide RNAs (gRNAs) in scRNA-seq. We use pAC-Seq to assess the phenotypic consequences of CRISPR/Cas9 based alterations of gene cis-regulatory regions. We show that pAC-Seq is able to detect cis-regulatory-induced alteration of target gene expression even when biallelic loss of target gene expression occurs in only ~5% of cells. This low rate of biallelic loss significantly increases the number of cells required to detect the consequences of changes to the regulatory genome, but can be ameliorated by transcript-targeted sequencing. Based on our experimental results we model the power to detect regulatory genome induced transcriptomic effects based on the rate of mono/biallelic loss, baseline gene expression, and the number of cells per target gRNA.

Fruitless decommissions regulatory elements to implement cell-type-specific neuronal masculinization

In the fruit fly Drosophila melanogaster, male-specific splicing and translation of the Fruitless transcription factor (Fru^M) alters the presence, anatomy, and/or connectivity of >60 types of central brain neurons that interconnect to generate male-typical behaviors. While the indispensable function of Fru^M in sex-specific behavior has been understood for decades, the molecular mechanisms underlying its activity remain unknown. Here, we take a genome-wide, brain-wide approach to identifying regulatory elements whose activity depends on the presence of Fru^M. We identify 436 high-confidence genomic regions differentially accessible in male fruitless neurons, validate candidate regions as bona fide, differentially regulated enhancers, and describe the particular cell types in which these enhancers are active. We find that individual enhancers are not activated universally but are dedicated to specific fru⁺ cell types. Aside from fru itself, genes are not dedicated to or common across the fru circuit; rather, Fru^M appears to masculinize each cell type differently, by tweaking expression of the same effector genes used in other circuits. Finally, we find Fru^M motifs enriched among regulatory elements that are open in the female but closed in the male. Together, these results suggest that Fru^M acts cell-type-specifically to decommission regulatory elements in male fruitless neurons.

Courtship behavior in male Drosophila melanogaster is controlled by a well-defined neural circuit that is labeled by the male-specific transcription factor Fruitless (Fru^M). While Fru^M is known to change the number, anatomy and connectivity of neurons which comprise the circuit and has been suggested to repress the expression of a few gene targets, the mechanism of how Fru^M regulates genes across many different kinds of neurons is unknown. Using an approach to identify gene regulatory elements based on their chromatin accessibility states (ATAC-seq), we identified a large set of chromatin accessibility changes downstream of Fruitless. By examining the activity of these regulatory elements in vivo, we found that their activity was 1) sexually dimorphic and 2) specific to a single class of Fru^M neurons, suggesting that Fru^M acts on different chromatin targets in different neuron classes comprising the courtship circuit. Further, we found a known Fru^M-regulated enhancer of the Fru^M-repressed gene Lgr3 to have closed chromatin specifically in Fru^M neurons. Combined with an enrichment of Fru^M motifs in regions which are closed in Fru^M neurons, we present a mechanism where Fru^M directs the decommissioning of sex-shared regulatory elements to masculinize neurons in a cell-type specific manner.

Inference of phenotype-relevant transcriptional regulatory networks elucidates cancer type-specific regulatory mechanisms in a pan-cancer study

Reconstruction of transcriptional regulatory networks (TRNs) is a powerful approach to unravel the gene expression programs involved in healthy and disease states of a cell. However, these networks are usually reconstructed independent of the phenotypic (or clinical) properties of the samples. Therefore, they may confound regulatory mechanisms that are specifically related to a phenotypic property with more general mechanisms underlying the full complement of the analyzed samples. In this study, we develop a method called InPheRNo to identify "phenotype-relevant" TRNs. This method is based on a probabilistic graphical model that models the simultaneous effects of multiple transcription factors (TFs) on their target genes and the statistical relationship between the target genes' expression and the phenotype. Extensive comparison of InPheRNo with related approaches using primary tumor samples of 18 cancer types from The Cancer Genome Atlas reveals that InPheRNo can accurately reconstruct cancer type-relevant TRNs and identify cancer driver TFs. In addition, survival analysis reveals that the activity level of TFs with many target genes could distinguish patients with poor prognosis from those with better prognosis.

HRT Atlas v1.0 database: redefining human and mouse housekeeping genes and candidate reference transcripts by mining massive RNA-seq datasets

Housekeeping (HK) genes are constitutively expressed genes that are required for the maintenance of basic cellular functions. Despite their importance in the calibration of gene expression, as well as the understanding of many genomic and evolutionary features, important discrepancies have been observed in studies that previously identified these genes. Here, we present Housekeeping and Reference Transcript Atlas (HRT Atlas v1.0, www.housekeeping.unicamp.br) a web-based database which addresses some of the previously observed limitations in the identification of these genes, and offers a more accurate database of human and mouse HK genes and transcripts. The database was generated by mining massive human and mouse RNA-seq data sets, including 11 281 and 507 high-quality RNA-seq samples from 52 human non-disease tissues/cells and 14 healthy tissues/cells of C57BL/6 wild type mouse, respectively. User can visualize the expression and download lists of 2158 human HK transcripts from 2176 HK genes and 3024 mouse HK transcripts from 3277 mouse HK genes. HRT Atlas also offers the most stable and suitable tissue selective candidate reference transcripts for normalization of qPCR experiments. Specific primers and predicted modifiers of gene expression for some of these HK transcripts are also proposed. HRT Atlas has also been integrated with a regulatory elements resource from Epiregio server.

Massively parallel techniques for cataloguing the regulome of the human brain

Complex brain disorders are highly heritable and arise from a complex polygenic risk architecture. Many disease-associated loci are found in non-coding regions that house regulatory elements. These elements influence the transcription of target genes-many of which demonstrate cell-type-specific expression patterns-and thereby affect phenotypically relevant molecular pathways. Thus, cell-type-specificity must be considered when prioritizing candidate risk loci, variants and target genes. This Review discusses the use of high-throughput assays in human induced pluripotent stem cell-based neurodevelopmental models to probe genetic risk in a cell-type- and patient-specific manner. The application of massively parallel reporter assays in human induced pluripotent stem cells can characterize the human regulome and test the transcriptional responses of putative regulatory elements. Parallel CRISPR-based screens can further functionally dissect this genetic regulatory architecture. The integration of these emerging technologies could decode genetic risk into medically actionable information, thereby improving genetic diagnosis and identifying novel points of therapeutic intervention.

Chromatin accessibility landscape and regulatory network of high-altitude hypoxia adaptation

High-altitude adaptation of Tibetans represents a remarkable case of natural selection during recent human evolution. Previous genome-wide scans found many non-coding variants under selection, suggesting a pressing need to understand the functional role of non-coding regulatory elements (REs). Here, we generate time courses of paired ATAC-seq and RNA-seq data on cultured HUVECs under hypoxic and normoxic conditions. We further develop a variant interpretation methodology (vPECA) to identify active selected REs (ASREs) and associated regulatory network. We discover three causal SNPs of EPAS1, the key adaptive gene for Tibetans. These SNPs decrease the accessibility of ASREs with weakened binding strength of relevant TFs, and cooperatively down-regulate EPAS1 expression. We further construct the downstream network of EPAS1, elucidating its roles in hypoxic response and angiogenesis. Collectively, we provide a systematic approach to interpret phenotype-associated noncoding variants in proper cell types and relevant dynamic conditions, to model their impact on gene regulation.

TENET 2.0: Identification of key transcriptional regulators and enhancers in lung adenocarcinoma

Lung cancer is the leading cause of cancer-related death and lung adenocarcinoma is its most common subtype. Although genetic alterations have been identified as drivers in subsets of lung adenocarcinoma, they do not fully explain tumor development. Epigenetic alterations have been implicated in the pathogenesis of tumors. To identify epigenetic alterations driving lung adenocarcinoma, we used an improved version of the Tracing Enhancer Networks using Epigenetic Traits method (TENET 2.0) in primary normal lung and lung adenocarcinoma cells. We found over 32,000 enhancers that appear differentially activated between normal lung and lung adenocarcinoma. Among the identified transcriptional regulators inactivated in lung adenocarcinoma vs. normal lung, NKX2-1 was linked to a large number of silenced enhancers. Among the activated transcriptional regulators identified, CENPA, FOXM1, and MYBL2 were linked to numerous cancer-specific enhancers. High expression of CENPA, FOXM1, and MYBL2 is particularly observed in a subgroup of lung adenocarcinomas and is associated with poor patient survival. Notably, CENPA, FOXM1, and MYBL2 are also key regulators of cancer-specific enhancers in breast adenocarcinoma of the basal subtype, but they are associated with distinct sets of activated enhancers. We identified individual lung adenocarcinoma enhancers linked to CENPA, FOXM1, or MYBL2 that were associated with poor patient survival. Knockdown experiments of FOXM1 and MYBL2 suggest that these factors regulate genes involved in controlling cell cycle progression and cell division. For example, we found that expression of TK1, a potential target gene of a MYBL2-linked enhancer, is associated with poor patient survival. Identification and characterization of key transcriptional regulators and associated enhancers in lung adenocarcinoma provides important insights into the deregulation of lung adenocarcinoma epigenomes, highlighting novel potential targets for clinical intervention.

Design of orthogonal regulatory systems for modulating gene expression in plants

Agricultural biotechnology strategies often require the precise regulation of multiple genes to effectively modify complex plant traits. However, most efforts are hindered by a lack of characterized tools that allow for reliable and targeted expression of transgenes. We have successfully engineered a library of synthetic transcriptional regulators that modulate expression strength in planta. By leveraging orthogonal regulatory systems from Saccharomyces spp., we have developed a strategy for the design of synthetic activators, synthetic repressors, and synthetic promoters and have validated their use in Nicotiana benthamiana and Arabidopsis thaliana. This characterization of contributing genetic elements that dictate gene expression represents a foundation for the rational design of refined synthetic regulators. Our findings demonstrate that these tools provide variation in transcriptional output while enabling the concerted expression of multiple genes in a tissue-specific and environmentally responsive manner, providing a basis for generating complex genetic circuits that process endogenous and environmental stimuli.

Transcriptional Regulators and Human-Specific/Primate-Specific Genes in Neocortical Neurogenesis

During development, starting from a pool of pluripotent stem cells, tissue-specific genetic programs help to shape and develop functional organs. To understand the development of an organ and its disorders, it is important to understand the spatio-temporal dynamics of the gene expression profiles that occur during its development. Modifications in existing genes, the de-novo appearance of new genes, or, occasionally, even the loss of genes, can greatly affect the gene expression profile of any given tissue and contribute to the evolution of organs or of parts of organs. The neocortex is evolutionarily the most recent part of the brain, it is unique to mammals, and is the seat of our higher cognitive abilities. Progenitors that give rise to this tissue undergo sequential waves of differentiation to produce the complete sets of neurons and glial cells that make up a functional neocortex. We will review herein our understanding of the transcriptional regulators that control the neural precursor cells (NPCs) during the generation of the most abundant class of neocortical neurons, the glutametergic neurons. In addition, we will discuss the roles of recently-identified human- and primate-specific genes in promoting neurogenesis, leading to neocortical expansion.

Investigating the effect of dependence between conditions with Bayesian Linear Mixed Models for motif activity analysis

MOTIVATION

Cellular identity and behavior is controlled by complex gene regulatory networks. Transcription factors (TFs) bind to specific DNA sequences to regulate the transcription of their target genes. On the basis of these TF motifs in cis-regulatory elements we can model the influence of TFs on gene expression. In such models of TF motif activity the data is usually modeled assuming a linear relationship between the motif activity and the gene expression level. A commonly used method to model motif influence is based on Ridge Regression. One important assumption of linear regression is the independence between samples. However, if samples are generated from the same cell line, tissue, or other biological source, this assumption may be invalid. This same assumption of independence is also applied to different yet similar experimental conditions, which may also be inappropriate. In theory, the independence assumption between samples could lead to loss in signal detection. Here we investigate whether a Bayesian model that allows for correlations results in more accurate inference of motif activities.

RESULTS

We extend the Ridge Regression to a Bayesian Linear Mixed Model, which allows us to model dependence between different samples. In a simulation study, we investigate the differences between the two model assumptions. We show that our Bayesian Linear Mixed Model implementation outperforms Ridge Regression in a simulation scenario where the noise, which is the signal that can not be explained by TF motifs, is uncorrelated. However, we demonstrate that there is no such gain in performance if the noise has a similar covariance structure over samples as the signal that can be explained by motifs. We give a mathematical explanation to why this is the case. Using four representative real datasets we show that at most ∼​40% of the signal is explained by motifs using the linear model. With these data there is no advantage to using the Bayesian Linear Mixed Model, due to the similarity of the covariance structure.

AVAILABILITY & IMPLEMENTATION

The project implementation is available at https://github.com/Sim19/SimGEXPwMotifs.

Effects of Certain cis-Regulatory Elements on Stage-Specific vitellogenin Expression in the Bombyx mori (Lepidoptera: Bombycidae)

Bombyx mori vitellogenin (BmVg) is highly upregulated during pupation, and the 20-hydroxyecdysone and amino acids may regulate stage-specific BmVg expression. However, previous studies showed that other factors may also affect stage-specific BmVg expression. Here, we characterized effective BmVg transcription factors by identifying the corresponding cis-regulatory elements (CREs). We prepared transgenic B. mori, in which DsRed was driven by various lengths of BmVg promoter. qRT-PCR analysis showed that DsRed expression driven by a 1.0-kb BmVg promoter (VgP1.0K) was consistent with endogenous BmVg. VgP1.0K specificity was closer to the endogenous BmVg promoter than that of VgP0.8K. These results suggest that CREs affecting stage-specific BmVg expression were localized to the 1.0-kb BmVg promoter. We investigated the effects of certain CREs that could influence the stage specificity of BmVg promoter on BmVg expression in transgenic B. mori. The relative DsRed expression was significantly reduced in transgenic female B. mori and the peak in DsRed expression was delayed after E-box CRE mutation. These results demonstrate that the E-box element enhanced BmVg expression and also affected stage-specific BmVg expression. Moreover, the relative DsRed expression was significantly increased in transgenic female of B. mori after 3×BD CRE mutation in BmVg promoter. However, the stage specificity of the mutated promoter was consistent with that of the endogenous BmVg promoter. The 3×BD element downregulated BmVg but had no effect on stage-specific BmVg expression. The present study promoted the process of elucidating the regulatory network for stage-specific BmVg expression and furnished a theoretical basis for the application of BmVg promoter.

Expression Analyses of Soybean VOZ Transcription Factors and the Role of GmVOZ1G in Drought and Salt Stress Tolerance

Vascular plant one-zinc-finger (VOZ) transcription factor, a plant specific one-zinc-finger-type transcriptional activator, is involved in regulating numerous biological processes such as floral induction and development, defense against pathogens, and response to multiple types of abiotic stress. Six VOZ transcription factor-encoding genes (GmVOZs) have been reported to exist in the soybean (Glycine max) genome. In spite of this, little information is currently available regarding GmVOZs. In this study, GmVOZs were cloned and characterized. GmVOZ genes encode proteins possessing transcriptional activation activity in yeast cells. GmVOZ1E, GmVOZ2B, and GmVOZ2D gene products were widely dispersed in the cytosol, while GmVOZ1G was primarily located in the nucleus. GmVOZs displayed a differential expression profile under dehydration, salt, and salicylic acid (SA) stress conditions. Among them, GmVOZ1G showed a significantly induced expression in response to all stress treatments. Overexpression of GmVOZ1G in soybean hairy roots resulted in a greater tolerance to drought and salt stress. In contrast, RNA interference (RNAi) soybean hairy roots suppressing GmVOZ1G were more sensitive to both of these stresses. Under drought treatment, soybean composite plants with an overexpression of hairy roots had higher relative water content (RWC). In response to drought and salt stress, lower malondialdehyde (MDA) accumulation and higher peroxidase (POD) and superoxide dismutase (SOD) activities were observed in soybean composite seedlings with an overexpression of hairy roots. The opposite results for each physiological parameter were obtained in RNAi lines. In conclusion, GmVOZ1G positively regulates drought and salt stress tolerance in soybean hairy roots. Our results will be valuable for the functional characterization of soybean VOZ transcription factors under abiotic stress.

Genome-wide identification and expression analysis of the phosphatase 2A family in rubber tree (Hevea brasiliensis)

The protein phosphatase 2As (PP2As) play a key role in manipulating protein phosphorylation. Although a number of proteins in the latex of laticifers are phosphorylated during latex regeneration in rubber tree, information about the PP2A family is limited. In the present study, 36 members of the HbPP2A family were genome-wide identified. They were clustered into five subgroups: the subgroup HbPP2AA (4), HbPP2AB' (14), HbPP2AB'' (6), HbPP2AB55 (4), and HbPP2AC (8). The members within the same subgroup shared highly conserved gene structures and protein motifs. Most of HbPP2As possessed ethylene- and wounding-responsive cis-acting elements. The transcripts of 29 genes could be detected in latex by using published high-throughput sequencing data. Of the 29 genes, seventeen genes were significantly down-regulated while HbPP2AA1-1 and HbPP2AB55α/Bα-1were up-regulated by tapping. Of the 17 genes, 14 genes were further significantly down-regulated by ethrel application. The down-regulated expression of a large number of HbPP2As may attribute to the enhanced phosphorylation of the proteins in latex from the tapped trees and the trees treated with ethrel application.

Genome-Wide Characterization and Analysis of CIPK Gene Family in Two Cultivated Allopolyploid Cotton Species: Sequence Variation, Association with Seed Oil Content, and the Role of GhCIPK6

Calcineurin B-like protein-interacting protein kinases (CIPKs), as key regulators, play an important role in plant growth and development and the response to various stresses. In the present study, we identified 80 and 78 CIPK genes in the Gossypium hirsutum and G. barbadense, respectively. The phylogenetic and gene structure analysis divided the cotton CIPK genes into five groups which were classified into an exon-rich clade and an exon-poor clade. A synteny analysis showed that segmental duplication contributed to the expansion of Gossypium CIPK gene family, and purifying selection played a major role in the evolution of the gene family in cotton. Analyses of expression profiles showed that GhCIPK genes had temporal and spatial specificity and could be induced by various abiotic stresses. Fourteen GhCIPK genes were found to contain 17 non-synonymous single nucleotide polymorphisms (SNPs) and co-localized with oil or protein content quantitative trait loci (QTLs). Additionally, five SNPs from four GhCIPKs were found to be significantly associated with oil content in one of the three field tests. Although most GhCIPK genes were not associated with natural variations in cotton oil content, the overexpression of the GhCIPK6 gene reduced the oil content and increased C18:1 and C18:1+C18:1d6 in transgenic cotton as compared to wild-type plants. In addition, we predicted the potential molecular regulatory mechanisms of the GhCIPK genes. In brief, these results enhance our understanding of the roles of CIPK genes in oil synthesis and stress responses.

Genomic Characterization of the Zinc Transcriptional Regulatory Element Reveals Potential Functional Roles of ZNF658

The zinc transcriptional regulatory element (ZTRE) is a newly reported binding motif for human zinc finger protein ZNF658, which alters gene expression in response to cellular zinc. The ZTRE has two nucleotide components-the palindromic flanking pairs and the bridging "N" bases between these flanks that range in number from 0 to 100. There are 12 pairs of ZTRE flanks (designated A-L). Three thousand five hundred twenty-five genes contain one or more ZTREs - 1000 to + 200 bp from their transcriptional start site (TSS). ZTRE-E is observed at a greater frequency, and ZTRE containing 25 bridging bases are less frequent, within - 200 bp from the TSS. The genes with ZTREs in this range are enriched in processes that may compensate zinc deficiency, while other genes with ZTREs outside this range are enriched in transcriptional activation processes. The division of ZTREs into two groups may imply a dual role of ZNF658, similar to the homologous yeast protein Zap1, via binding to low or high affinity sequences dependent upon cellular zinc. The KLF/Sp1-family binding motif is prevalent within the ZTRE "N" bridging bases, suggesting ZNF658 may compete with Sp1-like transactivators to suppress transcription.

Identification of gene specific cis-regulatory elements during differentiation of mouse embryonic stem cells: An integrative approach using high-throughput datasets

Gene expression governs cell fate, and is regulated via a complex interplay of transcription factors and molecules that change chromatin structure. Advances in sequencing-based assays have enabled investigation of these processes genome-wide, leading to large datasets that combine information on the dynamics of gene expression, transcription factor binding and chromatin structure as cells differentiate. While numerous studies focus on the effects of these features on broader gene regulation, less work has been done on the mechanisms of gene-specific transcriptional control. In this study, we have focussed on the latter by integrating gene expression data for the in vitro differentiation of murine ES cells to macrophages and cardiomyocytes, with dynamic data on chromatin structure, epigenetics and transcription factor binding. Combining a novel strategy to identify communities of related control elements with a penalized regression approach, we developed individual models to identify the potential control elements predictive of the expression of each gene. Our models were compared to an existing method and evaluated using the existing literature and new experimental data from embryonic stem cell differentiation reporter assays. Our method is able to identify transcriptional control elements in a gene specific manner that reflect known regulatory relationships and to generate useful hypotheses for further testing.

Direct prediction of regulatory elements from partial data without imputation

Genome segmentation approaches allow us to characterize regulatory states in a given cell type using combinatorial patterns of histone modifications and other regulatory signals. In order to analyze regulatory state differences across cell types, current genome segmentation approaches typically require that the same regulatory genomics assays have been performed in all analyzed cell types. This necessarily limits both the numbers of cell types that can be analyzed and the complexity of the resulting regulatory states, as only a small number of histone modifications have been profiled across many cell types. Data imputation approaches that aim to estimate missing regulatory signals have been applied before genome segmentation. However, this approach is computationally costly and propagates any errors in imputation to produce incorrect genome segmentation results downstream. We present an extension to the IDEAS genome segmentation platform which can perform genome segmentation on incomplete regulatory genomics dataset collections without using imputation. Instead of relying on imputed data, we use an expectation-maximization approach to estimate marginal density functions within each regulatory state. We demonstrate that our genome segmentation results compare favorably with approaches based on imputation or other strategies for handling missing data. We further show that our approach can accurately impute missing data after genome segmentation, reversing the typical order of imputation/genome segmentation pipelines. Finally, we present a new 2D genome segmentation analysis of 127 human cell types studied by the Roadmap Epigenomics Consortium. By using an expanded set of chromatin marks that have been profiled in subsets of these cell types, our new segmentation results capture a more complex picture of combinatorial regulatory patterns that appear on the human genome.

Distinct Genome Replication and Transcription Strategies within the Growing Filovirus Family

Research on filoviruses has historically focused on the highly pathogenic ebola- and marburgviruses. Indeed, until recently, these were the only two genera in the filovirus family. Recent advances in sequencing technologies have facilitated the discovery of not only a new ebolavirus, but also three new filovirus genera and a sixth proposed genus. While two of these new genera are similar to the ebola- and marburgviruses, the other two, discovered in saltwater fishes, are considerably more diverse. Nonetheless, these viruses retain a number of key features of the other filoviruses. Here, we review the key characteristics of filovirus replication and transcription, highlighting similarities and differences between the viruses. In particular, we focus on key regulatory elements in the genomes, replication and transcription strategies, and the conservation of protein domains and functions among the viruses. In addition, using computational analyses, we were able to identify potential homology and functions for some of the genes of the novel filoviruses with previously unknown functions. Although none of the newly discovered filoviruses have yet been isolated, initial studies of some of these viruses using minigenome systems have yielded insights into their mechanisms of replication and transcription. In general, the Cuevavirus and proposed Dianlovirus genera appear to follow the transcription and replication strategies employed by the ebola- and marburgviruses, respectively. While our knowledge of the fish filoviruses is currently limited to sequence analysis, the lack of certain conserved motifs and even entire genes necessitates that they have evolved distinct mechanisms of replication and transcription.

Identifying clusters of cis-regulatory elements underpinning TAD structures and lineage-specific regulatory networks

Cellular identity relies on cell-type-specific gene expression controlled at the transcriptional level by cis-regulatory elements (CREs). CREs are unevenly distributed across the genome, giving rise to individual CREs and clusters of CREs (COREs). Technical and biological features hinder CORE identification. We addressed these issues by developing an unsupervised machine learning approach termed clustering of genomic regions analysis method (CREAM). CREAM automates CORE detection from chromatin accessibility profiles that are enriched in CREs strongly bound by master transcription regulators, proximal to highly expressed and essential genes, and discriminating cell identity. Although COREs share similarities with super-enhancers, we highlight differences in terms of the genomic distribution and structure of these cis-regulatory units. We further show the enhanced value of COREs over super-enhancers to identify master transcription regulators, highly expressed and essential genes defining cell identity. COREs enrich at topologically associated domain (TAD) boundaries. They are also preferentially bound by the chromatin looping factors CTCF and cohesin, in contrast to super-enhancers, forming clusters of CTCF and cohesin binding regions and defining homotypic clusters of transcription regulator binding regions (HCTs). Finally, we show the clinical utility of CREAM to identify COREs across chromatin accessibility profiles to stratify more than 400 tumor samples according to their cancer type and to delineate cancer type-specific active biological pathways. Collectively, our results support the utility of CREAM to delineate COREs underlying, with greater accuracy than individual CREs or super-enhancers, the cell-type-specific biological underpinning across a wide range of normal and cancer cell types.

Regulatory interaction between the ZPBP2-ORMDL3/Zpbp2-Ormdl3 region and the circadian clock

Genome-wide association study (GWAS) loci for several immunity-mediated diseases (early onset asthma, inflammatory bowel disease (IBD), primary biliary cholangitis, and rheumatoid arthritis) map to chromosomal region 17q12-q21. The predominant view is that association between 17q12-q21 alleles and increased risk of developing asthma or IBD is due to regulatory variants. ORM sphingolipid biosynthesis regulator (ORMDL3) residing in this region is the most promising gene candidate for explaining association with disease. However, the relationship between 17q12-q21 alleles and disease is complex suggesting contributions from other factors, such as trans-acting genetic and environmental modifiers or circadian rhythms. Circadian rhythms regulate expression levels of thousands of genes and their dysregulation is implicated in the etiology of several common chronic inflammatory diseases. However, their role in the regulation of the 17q12-q21 genes has not been investigated. Moreover, the core clock gene nuclear receptor subfamily 1, group D, member 1 (NR1D1) resides about 200 kb distal to the GWAS region. We hypothesized that circadian rhythms influenced gene expression levels in 17q12-q21 region and conversely, regulatory elements in this region influenced transcription of the core clock gene NR1D1 in cis. To test these hypotheses, we examined the diurnal expression profiles of zona pellucida binding protein 2 (ZPBP2/Zpbp2), gasdermin B (GSDMB), and ORMDL3/Ormdl3 in human and mouse tissues and analyzed the impact of genetic variation in the ZPBP2/Zpbp2 region on NR1D1/Nr1d1 expression. We found that Ormdl3 and Zpbp2 were controlled by the circadian clock in a tissue-specific fashion. We also report that deletion of the Zpbp2 region altered the expression profile of Nr1d1 in lungs and ileum in a time-dependent manner. In liver, the deletion was associated with enhanced expression of Ormdl3. We provide the first evidence that disease-associated genes Zpbp2 and Ormdl3 are regulated by circadian rhythms and the Zpbp2 region influences expression of the core clock gene Nr1d1.

CRISPR/Cas9-based editing of a sensitive transcriptional regulatory element to achieve cell type-specific knockdown of the NEMO scaffold protein

The use of alternative promoters for the cell type-specific expression of a given mRNA/protein is a common cell strategy. NEMO is a scaffold protein required for canonical NF-κB signaling. Transcription of the NEMO gene is primarily controlled by two promoters: one (promoter B) drives NEMO transcription in most cell types and the second (promoter D) is largely responsible for NEMO transcription in liver cells. Herein, we have used a CRISPR/Cas9-based approach to disrupt a core sequence element of promoter B, and this genetic editing essentially eliminates expression of NEMO mRNA and protein in 293T human kidney cells. By cell subcloning, we have isolated targeted 293T cell lines that express no detectable NEMO protein, have defined genomic alterations at promoter B, and do not support activation of canonical NF-κB signaling in response to treatment with tumor necrosis factor. Nevertheless, non-canonical NF-κB signaling is intact in these NEMO-deficient cells. Expression of ectopic wild-type NEMO, but not certain human NEMO disease mutants, in the edited cells restores downstream NF-κB signaling in response to tumor necrosis factor. Targeting of the promoter B element does not substantially reduce NEMO expression (from promoter D) in the human SNU-423 liver cancer cell line. Thus, we have created a strategy for selectively eliminating cell type-specific expression from an alternative promoter and have generated 293T cell lines with a functional knockout of NEMO. The implications of these findings for further studies and for therapeutic approaches to target canonical NF-κB signaling are discussed.

What functional genomics has taught us about transcriptional regulation in malaria parasites

Malaria parasites are characterized by a complex life cycle that is accompanied by dynamic gene expression patterns. The factors and mechanisms that regulate gene expression in these parasites have been searched for even before the advent of next generation sequencing technologies. Functional genomics approaches have substantially boosted this area of research and have yielded significant insights into the interplay between epigenetic, transcriptional and post-transcriptional mechanisms. Recently, considerable progress has been made in identifying sequence-specific transcription factors and DNA-encoded regulatory elements. Here, we review the insights obtained from these efforts including the characterization of core promoters, the involvement of sequence-specific transcription factors in life cycle progression and the mapping of gene regulatory elements. Furthermore, we discuss recent developments in the field of functional genomics and how they might contribute to further characterization of this complex gene regulatory network.

RNAi-mediated depletion of the NSL complex subunits leads to abnormal chromosome segregation and defective centrosome duplication in Drosophila mitosis

The Drosophila Nonspecific Lethal (NSL) complex is a major transcriptional regulator of housekeeping genes. It contains at least seven subunits that are conserved in the human KANSL complex: Nsl1/Wah (KANSL1), Dgt1/Nsl2 (KANSL2), Rcd1/Nsl3 (KANSL3), Rcd5 (MCRS1), MBD-R2 (PHF20), Wds (WDR5) and Mof (MOF/KAT8). Previous studies have shown that Dgt1, Rcd1 and Rcd5 are implicated in centrosome maintenance. Here, we analyzed the mitotic phenotypes caused by RNAi-mediated depletion of Rcd1, Rcd5, MBD-R2 or Wds in greater detail. Depletion of any of these proteins in Drosophila S2 cells led to defects in chromosome segregation. Consistent with these findings, Rcd1, Rcd5 and MBD-R2 RNAi cells showed reduced levels of both Cid/CENP-A and the kinetochore component Ndc80. In addition, RNAi against any of the four genes negatively affected centriole duplication. In Wds-depleted cells, the mitotic phenotypes were similar but milder than those observed in Rcd1-, Rcd5- or MBD-R2-deficient cells. RT-qPCR experiments and interrogation of published datasets revealed that transcription of many genes encoding centromere/kinetochore proteins (e.g., cid, Mis12 and Nnf1b), or involved in centriole duplication (e.g., Sas-6, Sas-4 and asl) is substantially reduced in Rcd1, Rcd5 and MBD-R2 RNAi cells, and to a lesser extent in wds RNAi cells. During mitosis, both Rcd1-GFP and Rcd5-GFP accumulate at the centrosomes and the telophase midbody, MBD-R2-GFP is enriched only at the chromosomes, while Wds-GFP accumulates at the centrosomes, the kinetochores, the midbody, and on a specific chromosome region. Collectively, our results suggest that the mitotic phenotypes caused by Rcd1, Rcd5, MBD-R2 or Wds depletion are primarily due to reduced transcription of genes involved in kinetochore assembly and centriole duplication. The differences in the subcellular localizations of the NSL components may reflect direct mitotic functions that are difficult to detect at the phenotypic level, because they are masked by the transcription-dependent deficiency of kinetochore and centriolar proteins.

The Drosophila Nonspecific Lethal (NSL) complex is a conserved protein assembly that controls transcription of more than 4,000 housekeeping genes. We analyzed the mitotic functions of four genes, Rcd1, Rcd5, MBD-R2 and wds, encoding NSL subunits. Inactivation of these genes by RNA interference (RNAi) resulted in defects in both chromosome segregation and centrosome duplication. Our analyses indicate that RNAi against Rcd1, Rcd5 or MBD-R2 reduces transcription of genes involved in centromere/kinetochore assembly and centriole replication. During interphase, Rcd1, Rcd5, MBD-R2 and Wds are confined to the nucleus, as expected for transcription factors. However, during mitosis each of these proteins relocates to specific mitotic structures. Our results suggest that the four NSL components work together as a complex to stimulate transcription of genes encoding important mitotic determinants. However, the different localization of the proteins during mitosis suggests that they might have acquired secondary “moonlighting” functions that directly contribute to the mitotic process.

Mitigation of off-target toxicity in CRISPR-Cas9 screens for essential non-coding elements

Pooled CRISPR-Cas9 screens are a powerful method for functionally characterizing regulatory elements in the non-coding genome, but off-target effects in these experiments have not been systematically evaluated. Here, we investigate Cas9, dCas9, and CRISPRi/a off-target activity in screens for essential regulatory elements. The sgRNAs with the largest effects in genome-scale screens for essential CTCF loop anchors in K562 cells were not single guide RNAs (sgRNAs) that disrupted gene expression near the on-target CTCF anchor. Rather, these sgRNAs had high off-target activity that, while only weakly correlated with absolute off-target site number, could be predicted by the recently developed GuideScan specificity score. Screens conducted in parallel with CRISPRi/a, which do not induce double-stranded DNA breaks, revealed that a distinct set of off-targets also cause strong confounding fitness effects with these epigenome-editing tools. Promisingly, filtering of CRISPRi libraries using GuideScan specificity scores removed these confounded sgRNAs and enabled identification of essential regulatory elements.

Human evolved regulatory elements modulate genes involved in cortical expansion and neurodevelopmental disease susceptibility

Modern genetic studies indicate that human brain evolution is driven primarily by changes in gene regulation, which requires understanding the biological function of largely non-coding gene regulatory elements, many of which act in tissue specific manner. We leverage chromatin interaction profiles in human fetal and adult cortex to assign three classes of human-evolved elements to putative target genes. We find that human-evolved elements involving DNA sequence changes and those involving epigenetic changes are associated with human-specific gene regulation via effects on different classes of genes representing distinct biological pathways. However, both types of human-evolved elements converge on specific cell types and laminae involved in cerebral cortical expansion. Moreover, human evolved elements interact with neurodevelopmental disease risk genes, and genes with a high level of evolutionary constraint, highlighting a relationship between brain evolution and vulnerability to disorders affecting cognition and behavior. These results provide novel insights into gene regulatory mechanisms driving the evolution of human cognition and mechanisms of vulnerability to neuropsychiatric conditions.

Discovering Transcriptional Regulatory Elements From Run-On and Sequencing Data Using the Web-Based dREG Gateway

Transcription is a chromatin mark that can be used effectively to identify the location of active enhancers and promoters, collectively known as transcriptional regulatory elements (TREs). We recently introduced dREG, a tool for the identification of TREs using run-on and sequencing (RO-seq) assays, including global run-on and sequencing (GRO-seq), precision run-on and sequencing (PRO-seq), and chromatin run-on and sequencing (ChRO-seq). In this protocol, we present step-by-step instructions for running dREG on an arbitrary run-on and sequencing dataset. Users provide dREG with bigWig files (in which each read is represented by a single base) representing the location of RNA polymerase in a cell or tissue sample of interest, and dREG returns a list of genomic regions that are predicted to be active TREs. Finally, we demonstrate the use of dREG regions in discovering transcription factors controlling response to a stimulus and predicting their target genes. Together, this protocol provides detailed instructions for running dREG on arbitrary run-on and sequencing data. © 2018 by John Wiley & Sons, Inc.

Identification of the Regulatory Elements and Target Genes of Megakaryopoietic Transcription Factor MEF2C

Megakaryopoiesis produces specialized haematopoietic stem cells in the bone marrow that give rise to megakaryocytes which ultimately produce platelets. Defects in megakaryopoiesis can result in altered platelet counts and physiology, leading to dysfunctional haemostasis and thrombosis. Additionally, dysregulated megakaryopoiesis is also associated with myeloid pathologies. Transcription factors play critical roles in cell differentiation by regulating the temporal and spatial patterns of gene expression which ultimately decide cell fate. Several transcription factors have been described as regulating megakaryopoiesis including myocyte enhancer factor 2C (MEF2C); however, the genes regulated by MEF2C that influence megakaryopoiesis have not been reported. Using chromatin immunoprecipitation-sequencing and Gene Ontology data we identified five candidate genes that are bound by MEF2C and regulate megakaryopoiesis: MOV10, AGO3, HDAC1, RBBP5 and WASF2. To study expression of these genes, we silenced MEF2C gene expression in the Meg01 megakaryocytic cell line and in induced pluripotent stem cells by CRISPR/Cas9 editing. We also knocked down MEF2C expression in cord blood-derived haematopoietic stem cells by siRNA. We found that absent or reduced MEF2C expression resulted in defects in megakaryocytic differentiation and reduced levels of the candidate target genes. Luciferase assays confirmed that genomic sequences within the target genes are regulated by MEF2C levels. Finally, we demonstrate that small deletions linked to a platelet count-associated single nucleotide polymorphism alter transcriptional activity, suggesting a mechanism by which genetic variation in MEF2C alters platelet production. These data help elucidate the mechanism behind MEF2C regulation of megakaryopoiesis and genetic variation driving platelet production.

Evolutionary Origins of Pseudogenes and Their Association with Regulatory Sequences in Plants

Pseudogenes (Ψs), nonfunctional relatives of functional genes, form by duplication or retrotransposition, and loss of gene function by disabling mutations. Evolutionary analysis provides clues to Ψ origins and effects on gene regulation. However, few systematic studies of plant Ψs have been conducted, hampering comparative analyses. Here, we examined the origin, evolution, and expression patterns of Ψs and their relationships with noncoding sequences in seven angiosperm plants. We identified ∼250,000 Ψs, most of which are more lineage specific than protein-coding genes. The distribution of Ψs on the chromosome indicates that genome recombination may contribute to Ψ elimination. Most Ψs evolve rapidly in terms of sequence and expression levels, showing tissue- or stage-specific expression patterns. We found that a surprisingly large fraction of nontransposable element regulatory noncoding RNAs (microRNAs and long noncoding RNAs) originate from transcription of Ψ proximal upstream regions. We also found that transcription factor binding sites preferentially occur in putative Ψ proximal upstream regions compared with random intergenic regions, suggesting that Ψs have conditioned genome evolution by providing transcription factor binding sites that serve as promoters and enhancers. We therefore propose that rapid rewiring of Ψ transcriptional regulatory regions is a major mechanism driving the origin of novel regulatory modules.

Shared cis-regulatory architecture identified across defense response genes is associated with broad-spectrum quantitative resistance in rice

Plant disease resistance that is durable and effective against diverse pathogens (broad-spectrum) is essential to stabilize crop production. Such resistance is frequently controlled by Quantitative Trait Loci (QTL), and often involves differential regulation of Defense Response (DR) genes. In this study, we sought to understand how expression of DR genes is orchestrated, with the long-term goal of enabling genome-wide breeding for more effective and durable resistance. We identified short sequence motifs in rice promoters that are shared across Broad-Spectrum DR (BS-DR) genes co-expressed after challenge with three major rice pathogens (Magnaporthe oryzae, Rhizoctonia solani, and Xanthomonas oryzae pv. oryzae) and several chemical elicitors. Specific groupings of these BS-DR-associated motifs, called cis-Regulatory Modules (CRMs), are enriched in DR gene promoters, and the CRMs include cis-elements known to be involved in disease resistance. Polymorphisms in CRMs occur in promoters of genes in resistant relative to susceptible BS-DR haplotypes providing evidence that these CRMs have a predictive role in the contribution of other BS-DR genes to resistance. Therefore, we predict that a CRM signature within BS-DR gene promoters can be used as a marker for future breeding practices to enrich for the most responsive and effective BS-DR genes across the genome.

The cis-Regulatory Atlas of the Mouse Immune System

A complete chart of cis-regulatory elements and their dynamic activity is necessary to understand the transcriptional basis of differentiation and function of an organ system. We generated matched epigenome and transcriptome measurements in 86 primary cell types that span the mouse immune system and its differentiation cascades. This breadth of data enable variance components analysis that suggests that genes fall into two distinct classes, controlled by either enhancer- or promoter-driven logic, and multiple regression that connects genes to the enhancers that regulate them. Relating transcription factor (TF) expression to the genome-wide accessibility of their binding motifs classifies them as predominantly openers or closers of local chromatin accessibility, pinpointing specific cis-regulatory elements where binding of given TFs is likely functionally relevant, validated by chromatin immunoprecipitation sequencing (ChIP-seq). Overall, this cis-regulatory atlas provides a trove of information on transcriptional regulation through immune differentiation and a foundational scaffold to define key regulatory events throughout the immunological genome.

PlantPAN3.0: a new and updated resource for reconstructing transcriptional regulatory networks from ChIP-seq experiments in plants

The Plant Promoter Analysis Navigator (PlantPAN; http://PlantPAN.itps.ncku.edu.tw/) is an effective resource for predicting regulatory elements and reconstructing transcriptional regulatory networks for plant genes. In this release (PlantPAN 3.0), 17 230 TFs were collected from 78 plant species. To explore regulatory landscapes, genomic locations of TFBSs have been captured from 662 public ChIP-seq samples using standard data processing. A total of 1 233 999 regulatory linkages were identified from 99 regulatory factors (TFs, histones and other DNA-binding proteins) and their target genes across seven species. Additionally, this new version added 2449 matrices extracted from ChIP-seq peaks for cis-regulatory element prediction. In addition to integrated ChIP-seq data, four major improvements were provided for more comprehensive information of TF binding events, including (i) 1107 experimentally verified TF matrices from the literature, (ii) gene regulation network comparison between two species, (iii) 3D structures of TFs and TF-DNA complexes and (iv) condition-specific co-expression networks of TFs and their target genes extended to four species. The PlantPAN 3.0 can not only be efficiently used to investigate critical cis- and trans-regulatory elements in plant promoters, but also to reconstruct high-confidence relationships among TF-targets under specific conditions.

Reveal cell type-specific regulatory elements and their characterized histone code classes via a hidden Markov model

BACKGROUND

With the maturity of next generation sequencing technology, a huge amount of epigenomic data have been generated by several large consortia in the last decade. These plenty resources leave us the opportunity about sufficiently utilizing those data to explore biological problems.

RESULTS

Here we developed an integrative and comparative method, CsreHMM, which is based on a hidden Markov model, to systematically reveal cell type-specific regulatory elements (CSREs) along the whole genome, and simultaneously recognize the histone codes (mark combinations) charactering them. This method also reveals the subclasses of CSREs and explicitly label those shared by a few cell types. We applied this method to a data set of 9 cell types and 9 chromatin marks to demonstrate its effectiveness and found that the revealed CSREs relates to different kinds of functional regulatory regions significantly. Their proximal genes have consistent expression and are likely to participate in cell type-specific biological functions.

CONCLUSIONS

These results suggest CsreHMM has the potential to help understand cell identity and the diverse mechanisms of gene regulation.

Loss of hierarchical imprinting regulation at the Prader-Willi/Angelman syndrome locus in human iPSCs

The human chr15q11-q13 imprinted cluster is linked to several disorders, including Prader-Willi (PWS) and Angelman (AS) syndromes. Recently, disease modeling approaches based on induced pluripotent stem cells (iPSCs) have been used to study these syndromes. A concern regarding the use of these cells for imprinted disease modeling is the numerous imprinting defects found in many iPSCs. Here, by reprogramming skin fibroblasts from a control and AS individuals, we generated several iPSC lines and addressed the stability of imprinting status across the PWS/AS domain. We focused on three important regulatory DNA elements which are all differentially methylated regions (DMRs), methylated on the maternal allele: the PWS imprinting center (PWS-IC), which is a germline DMR and the somatic NDN and MKRN3 DMRs, hierarchically controlled by PWS-IC. Normal PWS-IC methylation pattern was maintained in most iPSC lines; however, loss of maternal methylation in one out of five control iPSC lines resulted in a monoallelic to biallelic switch for many imprinted genes in this domain. Surprisingly, MKRN3 DMR was found aberrantly hypermethylated in all control and AS iPSCs, regardless of the methylation status of the PWS-IC master regulator. This suggests a loss of hierarchical control of imprinting at PWS/AS region. We confirmed these results in established iPSC lines derived using different reprogramming procedures. Overall, we show that hierarchy of imprinting control in donor cells might not apply to iPSCs, accounting for their spectrum of imprinting alterations. Such differences in imprinting regulation should be taken into consideration for the use of iPSCs in disease modeling.

Towards a map of cis-regulatory sequences in the human genome

Accumulating evidence indicates that transcription factor (TF) binding sites, or cis-regulatory elements (CREs), and their clusters termed cis-regulatory modules (CRMs) play a more important role than do gene-coding sequences in specifying complex traits in humans, including the susceptibility to common complex diseases. To fully characterize their roles in deriving the complex traits/diseases, it is necessary to annotate all CREs and CRMs encoded in the human genome. However, the current annotations of CREs and CRMs in the human genome are still very limited and mostly coarse-grained, as they often lack the detailed information of CREs in CRMs. Here, we integrated 620 TF ChIP-seq datasets produced by the ENCODE project for 168 TFs in 79 different cell/tissue types and predicted an unprecedentedly completely map of CREs in CRMs in the human genome at single nucleotide resolution. The map includes 305 912 CRMs containing a total of 1 178 913 CREs belonging to 736 unique TF binding motifs. The predicted CREs and CRMs tend to be subject to either purifying selection or positive selection, thus are likely to be functional. Based on the results, we also examined the status of available ChIP-seq datasets for predicting the entire regulatory genome of humans.

DeFine: deep convolutional neural networks accurately quantify intensities of transcription factor-DNA binding and facilitate evaluation of functional non-coding variants

The complex system of gene expression is regulated by the cell type-specific binding of transcription factors (TFs) to regulatory elements. Identifying variants that disrupt TF binding and lead to human diseases remains a great challenge. To address this, we implement sequence-based deep learning models that accurately predict the TF binding intensities to given DNA sequences. In addition to accurately classifying TF-DNA binding or unbinding, our models are capable of accurately predicting real-valued TF binding intensities by leveraging large-scale TF ChIP-seq data. The changes in the TF binding intensities between the altered sequence and the reference sequence reflect the degree of functional impact for the variant. This enables us to develop the tool DeFine (Deep learning based Functional impact of non-coding variants evaluator, http://define.cbi.pku.edu.cn) with improved performance for assessing the functional impact of non-coding variants including SNPs and indels. DeFine accurately identifies the causal functional non-coding variants from disease-associated variants in GWAS. DeFine is an effective and easy-to-use tool that facilities systematic prioritization of functional non-coding variants.

Concomitant Regulation by a LacI-Type Transcriptional Repressor XylR on Genes Involved in Xylan and Xylose Metabolism and the Type III Secretion System in Rice Pathogen Xanthomonas oryzae pv. oryzae

The hypersensitive response and pathogenicity (hrp) genes of Xanthomonas oryzae pv. oryzae, the causal agent of bacterial leaf blight of rice, encode components of the type III secretion system and are essential for virulence. Expression of hrp genes is regulated by two key hrp regulators, HrpG and HrpX; HrpG regulates hrpX and hrpA, and HrpX regulates the other hrp genes on hrpB-hrpF operons. We previously reported the sugar-dependent quantitative regulation of HrpX; the regulator highly accumulates in the presence of xylose, followed by high hrp gene expression. Here, we found that, in a mutant lacking the LacI-type transcriptional regulator XylR, HrpX accumulation and hrp gene expression were high even in the medium without xylose, reaching the similar levels present in the wild type incubated in the xylose-containing medium. XylR also negatively regulated one of two xylose isomerase genes (xylA2 but not xylA1) by binding to the motif sequence in the upstream region of the gene. Xylose isomerase is an essential enzyme in xylose metabolism and interconverts between xylose and xylulose. Our results suggest that, in the presence of xylose, inactivation of XylR leads to greater xylan and xylose utilization and, simultaneously, to higher accumulation of HrpX, followed by higher hrp gene expression in the bacterium.

Diagnostic utility of DREAM gene mRNA levels in thyroid tumours

OBJECTIVE

The transcriptional repressor DREAM is involved in thyroid-specific gene expression, thyroid enlargement and nodular development, but its clinical utility is still uncertain. In this study we aimed to investigate whether DREAM mRNA levels differ in different thyroid tumors and how this possible difference would allow the use of DREAM gene expression as molecular marker for diagnostic and/or prognosis purpose.

MATERIALS AND METHODS

We quantified DREAM gene mRNA levels and investigated its mutational status, relating its expression and genetic changes to diagnostic and prognostic features of 200 thyroid tumors, being 101 malignant [99 papillary thyroid carcinomas (PTC) and 2 anaplastic thyroid carcinomas] and 99 benign thyroid lesions [49 goiter and 50 follicular adenomas (FA)].

RESULTS

Levels of mRNA of DREAM gene were higher in benign (0.7909 ± 0.6274 AU) than in malignant (0.3373 ± 0.6274 AU) thyroid lesions (p < 0.0001). DREAM gene expression was able to identify malignancy with 66.7% sensitivity, 85.4% specificity, 84.2% positive predictive value (PPV), 68.7% negative predictive value (NPV), and 75.3% accuracy. DREAM mRNA levels were also useful distinguishing the follicular lesions FA and FVPTC with 70.2% sensitivity, 73.5% specificity, 78.5% PPV, 64.1% NPV, and 71.6% accuracy. However, DREAM gene expression was neither associated with clinical features of tumor aggressiveness, nor with recurrence or survival. Six different genetic changes in non-coding regions of DREAM gene were also found, not related to DREAM gene expression or tumor features.

CONCLUSION

We suggest that DREAM gene expression may help diagnose thyroid nodules, identifying malignancy and characterizing follicular-patterned thyroid lesions; however, it is not useful as a prognostic marker.

Cis-regulatory landscapes of four cell types of the retina

The retina is composed of ∼50 cell-types with specific functions for the process of vision. Identification of the cis-regulatory elements active in retinal cell-types is key to elucidate the networks controlling this diversity. Here, we combined transcriptome and epigenome profiling to map the regulatory landscape of four cell-types isolated from mouse retinas including rod and cone photoreceptors as well as rare inter-neuron populations such as horizontal and starburst amacrine cells. Integration of this information reveals sequence determinants and candidate transcription factors for controlling cellular specialization. Additionally, we refined parallel reporter assays to enable studying the transcriptional activity of large collection of sequences in individual cell-types isolated from a tissue. We provide proof of concept for this approach and its scalability by characterizing the transcriptional capacity of several hundred putative regulatory sequences within individual retinal cell-types. This generates a catalogue of cis-regulatory regions active in retinal cell types and we further demonstrate their utility as potential resource for cellular tagging and manipulation.

DNA Sequence Constraints Define Functionally Active Steroid Nuclear Receptor Binding Sites in Chromatin

Gene regulatory programs are encoded in the sequence of the DNA. Since the completion of the Human Genome Project, millions of gene regulatory elements have been identified in the human genome. Understanding how each of those sites functionally contributes to gene regulation, however, remains a challenge for nearly every field of biology. Transcription factors influence cell function by interpreting information contained within cis-regulatory elements in chromatin. Whereas chromatin immunoprecipitation-sequencing has been used to identify and map transcription factor-DNA interactions, it has been difficult to assign functionality to the binding sites identified. Thus, in this study, we probed the transcriptional activity, DNA-binding competence, and functional activity of select nuclear receptor mutants in cellular and animal model systems and used this information to define the sequence constraints of functional steroid nuclear receptor cis-regulatory elements. Analysis of the architecture within sNR chromatin interacting sites revealed that only a small fraction of all sNR chromatin-interacting events is associated with transcriptional output and that this functionality is restricted to elements that vary from the consensus palindromic elements by one or two nucleotides. These findings define the transcriptional grammar necessary to predict functionality from regulatory sequences, with a multitude of future implications.

Functional analysis of the promoters of the small GTPases RhoA and RhoB in embryonic stem cells

Small GTPases of the Rho subfamily have been implicated in many physiological and pathological processes in various cell types including embryonic stem cells (ESCs). In the present study we performed a functional analysis of the promoters of the RhoA and the RhoB genes in order to identify regulatory elements that are important for their transcriptional control in ESCs. We show that RhoA mRNA levels were significantly higher compared with the RhoB mRNA levels in ESCs as well in various cancer cell lines and this difference could be accounted for by differences in the activities of the corresponding promoters. Deletion analysis of the RhoA and RhoB promoters in ESCs revealed that the proximal regions contain regulatory elements that are critical for their activity. Both proximal promoters contain CCAAT boxes and mutagenesis of these elements decreased significantly the activity of both promoters suggesting a coordinated regulation of the two genes by CCAAT box binding factors. Finally, we show that both genes are subjects to autoregulation in ESCs and in the case of RhoB, this autoregulation requires the GTPase activity of the Rho proteins. Understanding the mechanisms that control the transcription of Rho GTPases in ESCs may shed new light into the still unknown roles of these proteins in stem cell functions.