Metazoan promoters: emerging characteristics and insights into transcriptional regulation

Are complex traits underpinned by polygenic molecular traits? A reflection on the complexity of gene expression

Variation in complex traits is controlled by multiple genes. The prevailing assumption is that such polygenic complex traits are underpinned by variation in elementary molecular traits, such as gene expression, which themselves have a simple genetic basis. Here, we review recent advances that reveal the captivating complexity of gene regulation: the cell type, time point, and magnitude of gene expression are not merely dependent on a couple of regulators; rather, they result from a probabilistic process shaped by cis- and trans-regulatory elements collaboratively integrating internal and external cues with the tightly regulated dynamics of DNA. In addition, the finding that genetic variants linked to complex diseases in humans often do not co-localize with quantitative trait loci modulating gene expression, along with the role of nonfunctional transcription factor (TF) binding sites, suggests that some of the genetic effects influencing gene expression variation may be indirect. If the number of genomic positions responsible for TF binding, TF binding site search time, DNA conformation and accessibility as well as regulation of all trans-acting factors is indeed vast, is it plausible that the complexity of elementary molecular traits approaches the complexity of higher-level organismal traits? Although it is hard to know the answer to this question, we motivate it by reviewing the complexity of the molecular machinery further.

DNA hypermethylation preceded by H3K27 trimethylation is linked to downregulation of gene expression in disuse muscle atrophy in male mice

Disuse muscle atrophy can result in downregulated gene expression vital to muscle integrity, yet the mechanisms driving this downregulation remain unclear. Epigenetic alterations regulate transcriptional potential, with repressive changes suppressing gene expression. This study explored epigenetic mechanisms of gene downregulation during disuse muscle atrophy. Male C57BL/6J mice underwent hindlimb suspension for 3 or 7 days. The vastus intermedius (VI) muscle was analyzed, showing unchanged mass on day 3, but on day 7, decreased mass and reduced fiber size were assessed via immunohistochemistry. Corresponding to this atrophy timing, qPCR analysis revealed nine downregulated genes on day 7, which were selected for epigenetic analysis; collectively, they showed no downregulation on day 3. Among the nine genes, methylated DNA immunoprecipitation revealed significantly elevated DNA methylation (hypermethylation) in the upstream regions of transcription start sites (TSS) on day 7, which overall negatively correlated with gene expression. Histone marks (H3K27me3, H3K4me3, H3.3, and total H3) were also assessed using chromatin immunoprecipitation, revealing that the repressive histone mark H3K27me3 increased in the regions on day 3 but decreased on day 7. These findings suggest that DNA hypermethylation in the upstream regions preceded by H3K27me3 enrichment contributes to the downregulation of gene expression during disuse muscle atrophy.

True length of diverse capped RNA sequencing (TLDR-seq): 5'-3'-end sequencing of capped RNAs regardless of 3'-end status

Analysis of transcript function is greatly aided by knowledge of the full-length RNA sequence. New long-read sequencing enabled by Oxford Nanopore and PacBio devices have the potential to provide full-length transcript information; however, standard methods still lack the ability to capture true RNA 5' ends and select for polyadenylated (pA+) transcripts only. Here, we present a method that, by utilizing cap trapping and 3'-end adapter ligation, sequences transcripts between their exact 5' and 3' ends regardless of polyadenylation status and without the need for ribosomal RNA depletion, with the ability to characterize polyadenylation length of RNAs, if any. The method shows high reproducibility, can faithfully detect 5' ends, 3' ends and splice junctions, and produces gene-expression estimates that are highly correlated to those of short-read sequencing techniques. We also demonstrate that the method can detect and sequence full-length nonadenylated (pA-) RNAs, including long noncoding RNAs, promoter upstream transcripts, and enhancer RNAs, and present cases where pA+ and pA- RNAs show preferences for different but closely located transcription start sites. Our method is therefore useful for the characterization of diverse capped RNA species and analysis of relationships between transcription initiation, termination, and RNA processing.

Simultaneous measurement of intrinsic promoter and enhancer potential reveals principles of functional duality and regulatory reciprocity

Growing evidence indicates that transcriptional regulatory elements can exert both promoter and enhancer activity; however, the relationship and determinants of this dual functionality remain poorly understood. We developed a massively parallel dual reporter assay that enables simultaneous assessment of the intrinsic promoter and enhancer potential exerted by the same sequence. Parallel quantification for thousands of elements reveals that canonical human promoters and enhancers can act as both promoters and enhancers under the same contexts, and that promoter activity may be necessary but not sufficient for enhancer function. We find that regulatory potential is intrinsic to element sequences, irrespective of downstream features typically associated with distinct element classes. Perturbations to element transcription factor binding motifs lead to disruptions in both activities, implicating a shared syntax for the two regulatory functions. Combinations of elements with different minimal promoters reveal reciprocal activity modulation between associated elements and a strong positive correlation between promoter and enhancer functions imply a bidirectional feedback loop used to maintain environments of high transcriptional activity. Finally, our results indicate that the magnitude and balance between promoter and enhancer functions are shaped by both intrinsic sequence properties and contextual regulatory influences, suggesting a degree of plasticity in regulatory action. Our approach provides a new lens for understanding fundamental principles of regulatory element biology.

Predicting gene expression responses to environment in Arabidopsis thaliana using natural variation in DNA sequence

The evolution of gene expression responses is a critical component of adaptation to variable environments. Predicting how DNA sequence influences expression is challenging because the genotype to phenotype map is not well resolved for cis regulatory elements, transcription factor binding, regulatory interactions, and epigenetic features, not to mention how these factors respond to environment. We tested if flexible machine learning models could learn some of the underlying cis-regulatory genotype to phenotype map. We tested this approach using cold-responsive transcriptome profiles in 5 diverse Arabidopsis thaliana accessions. We first tested for evidence that cis regulation plays a role in environmental response, finding 14 and 15 motifs that were significantly enriched within the up- and down-stream regions of cold-responsive differentially regulated genes (DEGs). We next applied convolutional neural networks (CNNs), which learn de novo cis-regulatory motifs in DNA sequences to predict expression response to environment. We found that CNNs predicted differential expression with moderate accuracy, with evidence that predictions were hindered by biological complexity of regulation and the large potential regulatory code. Overall, approaches to predict DEGs between specific environments based only on proximate DNA sequences require further development, and additional information may be required.

RAEPI: Predicting Enhancer-Promoter Interactions Based on Restricted Attention Mechanism

Enhancer-promoter interactions (EPIs) are crucial in gene transcription regulation and cell differentiation. Traditional biological experiments are costly and time-consuming, motivating the development of computational prediction methods. However, existing EPI prediction methods inadequately capture the intricate direct interactions between enhancer and promoter sequences, which limits their prediction performance to some extent. In this work, we propose an innovative attention-based approach RAEPI, which uses convolutional neural networks to extract initial features of enhancers and promoters, combined with a specially designed Restricted Attention mechanism with Query-Key-Value constrained to simulate the interactions between them for further feature extraction. To improve cross-cell line prediction, we employ a transfer learning strategy for pre-training. Furthermore, we extracted sequence motifs to evaluate the RAEPI's effectiveness from a visualization perspective. Experimental results show that RAEPI achieves competitive prediction performance to existing methods on the benchmark dataset.

The combinatorial binding syntax of transcription factors in forebrain-specific enhancers

Tissue-specific gene regulation in mammals involves the coordinated binding of multiple transcription factors (TFs). Using the forebrain as a model, we investigated the syntax of TF occupancy to determine tissue-specific enhancer regions. We analyzed forebrain-exclusive enhancers from the VISTA Enhancer Browser and a curated set of 23 TFs relevant to forebrain development and disease. Our findings revealed multiple distinct patterns of combinatorial TF binding, with the HES5-FOXP2-GATA3 triad being the most frequent in forebrain-specific enhancers. This syntactic structure was detected in 2614 enhancers from a genome-wide catalog of 25,000 predicted human forebrain enhancers. Notably, this catalog represents a computationally predicted dataset, distinct from the in vivo validated set of enhancers obtained from the VISTA Enhancer Browser. The shortlisted 2614 enhancers were further analyzed using genome-wide epigenetic data and evaluated for evolutionary conservation and disease relevance. Our findings highlight the value of these 2614 enhancers in forebrain-specific gene regulation and provide a framework for discovering tissue-specific enhancers, enhancing the understanding of enhancer function.

Universal prediction of vertebrate species age at maturity

Animal age at maturity can be used as a universal and simple predictor of species extinction risk. At present, methods to estimate age at maturity are typically species-specific, limiting comparisons among species, or are infeasible due to practical constraints. To overcome this, here we develop a universal predictor of species-level age at maturity for vertebrates. We show that modelling the frequency of 'CG' sequences (CpG sites) in gene promoter regions yields rapid predictions of vertebrate age at maturity. Our models predict age at maturity with remarkable accuracy and generalisability, with median error rates of 30% (less than 1 year) and are robust to genome assemblies of varying quality. We generate predictions for 1912 vertebrate species for which age at maturity estimates were previously absent from public databases. The predictions can be used to help to inform management decisions for the many species for which more detailed population information is currently unavailable.

Oncogenic Enhancers in Leukemia

Although the study of leukemogenesis has traditionally focused on protein-coding genes, the role of enhancer dysregulation is becoming increasingly recognized. The advent of high-throughput sequencing, together with a better understanding of enhancer biology, has revealed how various genetic and epigenetic lesions produce oncogenic enhancers that drive transformation. These aberrations include translocations that lead to enhancer hijacking, point mutations that modulate enhancer activity, and copy number alterations that modify enhancer dosage. In this review, we describe these mechanisms in the context of leukemia and discuss potential therapeutic avenues to target these regulatory elements. Significance: Large-scale sequencing projects have uncovered recurrent gene mutations in leukemia, but the picture remains incomplete: some patients harbor no such aberrations, whereas others carry only a few that are insufficient to bring about transformation on their own. One of the missing pieces is enhancer dysfunction, which only recently has emerged as a critical driver of leukemogenesis. Knowledge of the various mechanisms of enhancer dysregulation is thus key for a complete understanding of leukemia and its causes, as well as the development of targeted therapies in the era of precision medicine.

Cooperative insulation of regulatory domains by CTCF-dependent physical insulation and promoter competition

The specificity of gene expression during development requires the insulation of regulatory domains to avoid inappropriate enhancer-gene interactions. In vertebrates, this insulator function is mostly attributed to clusters of CTCF sites located at topologically associating domain (TAD) boundaries. However, TAD boundaries allow some physical crosstalk across regulatory domains, which is at odds with the specific and precise expression of developmental genes. Here we show that developmental genes and nearby clusters of CTCF sites cooperatively foster the robust insulation of regulatory domains. By genetically dissecting a couple of representative loci in mouse embryonic stem cells, we show that CTCF sites prevent undesirable enhancer-gene contacts (i.e. physical insulation), while developmental genes preferentially contribute to regulatory insulation through non-structural mechanisms involving promoter competition rather than enhancer blocking. Overall, our work provides important insights into the insulation of regulatory domains, which in turn might help interpreting the pathological consequences of certain structural variants.

Genome-Wide Association Study for Weight Loss at the End of Dry-Curing of Hams Produced from Purebred Heavy Pigs

Dissecting the genetics of production traits in livestock is of outmost importance, both to understand biological mechanisms underlying those traits and to facilitate the design of selection programs incorporating that information. For the pig industry, traits related to curing are key for protected designation of origin productions. In particular, appropriate ham weight loss after dry-curing ensures high quality of the final product and avoids economic losses. In this study, we analyzed data (N = 410) of ham weight loss after approximately 20 months of dry-curing. The animals used for ham production were purebred pigs belonging to a commercial line. A genome-wide association study (GWAS) of 29,844 SNP markers revealed the polygenic nature of the trait: 221 loci explaining a small percentage of the variance (0.3-1.65%) were identified on almost all Sus scrofa chromosomes. Post-GWAS analyses revealed 32 windows located within regulatory regions and 94 windows located in intronic regions of specific genes. In total, 30 candidate genes encoding receptors and enzymes associated with ham weight loss (MTHFD1L, DUSP8), proteolysis (SPARCL1, MYH8), drip loss (TNNI2), growth (CDCA3, LSP1, CSMD1, AP2A2, TSPAN4), and fat metabolism (AGPAT4, IGF2R, PTDSS2, HRAS, TALDO1, BRSK2, TNNI2, SYT8, GTF2I, GTF2IRD1, LPCAT3, ATN1, GNB3, CMIP, SORCS2, CCSER1, SPP1) were detected.

Alternative TSS use is widespread in Cryptococcus fungi in response to environmental cues and regulated genome-wide by the transcription factor Tur1

Alternative transcription start site (TSS) usage regulation has been identified as a major means of gene expression regulation in metazoans. However, in fungi, its impact remains elusive as its study has thus far been restricted to model yeasts. Here, we first re-analyzed TSS-seq data to define genuine TSS clusters in 2 species of pathogenic Cryptococcus. We identified 2 types of TSS clusters associated with specific DNA sequence motifs. Our analysis also revealed that alternative TSS usage regulation in response to environmental cues is widespread in Cryptococcus, altering gene expression and protein targeting. Importantly, we performed a forward genetic screen to identify a unique transcription factor (TF) named Tur1, which regulates alternative TSS (altTSS) usage genome-wide when cells switch from exponential phase to stationary phase. ChiP-Seq and DamID-Seq analyses suggest that at some loci, the role of Tur1 might be direct. Tur1 has been previously shown to be essential for virulence in C. neoformans. We demonstrated here that a tur1Δ mutant strain is more sensitive to superoxide stress and phagocytosed more efficiently by macrophages than the wild-type (WT) strain.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically shows phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption-the root of the pathogenesis-is similar in the different disease-relevant cell types. This is possible in principle, because all these cell types are subject to effects from the same causative gene, which has the same kind of function (e.g., methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2 and find that the chromatin accessibility changes in neurons are mostly distinct from changes in B or T cells. This is not because the neuronal accessibility changes occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that regulatory elements disrupted in B/T cells do show chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators and suggest that blood-derived episignatures, although useful diagnostically, may not be well suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Integrative analysis of transcriptomic and epigenomic data reveals distinct patterns for developmental and housekeeping gene regulation

BACKGROUND

Regulation of transcription is central to the emergence of new cell types during development, and it often involves activation of genes via proximal and distal regulatory regions. The activity of regulatory elements is determined by transcription factors (TFs) and epigenetic marks, but despite extensive mapping of such patterns, the extraction of regulatory principles remains challenging.

RESULTS

Here we study differentially and similarly expressed genes along with their associated epigenomic profiles, chromatin accessibility and DNA methylation, during lineage specification at gastrulation in mice. Comparison of the three lineages allows us to identify genomic and epigenomic features that distinguish the two classes of genes. We show that differentially expressed genes are primarily regulated by distal elements, while similarly expressed genes are controlled by proximal housekeeping regulatory programs. Differentially expressed genes are relatively isolated within topologically associated domains, while similarly expressed genes tend to be located in gene clusters. Transcription of differentially expressed genes is associated with differentially open chromatin at distal elements including enhancers, while that of similarly expressed genes is associated with ubiquitously accessible chromatin at promoters.

CONCLUSION

Based on these associations of (linearly) distal genes' transcription start sites (TSSs) and putative enhancers for developmental genes, our findings allow us to link putative enhancers to their target promoters and to infer lineage-specific repertoires of putative driver transcription factors, within which we define subgroups of pioneers and co-operators.

Developmental and Housekeeping Genes: Two Types of Genetic Organization in the Drosophila Genome

We developed a procedure for locating genes on Drosophila melanogaster polytene chromosomes and described three types of chromosome structures (gray bands, black bands, and interbands), which differed markedly in morphological and genetic properties. This was reached through the use of our original methods of molecular and genetic analysis, electron microscopy, and bioinformatics data processing. Analysis of the genome-wide distribution of these properties led us to a bioinformatics model of the Drosophila genome organization, in which the genome was divided into two groups of genes. One was constituted by 65, in which the genome was divided into two groups, 62 genes that are expressed in most cell types during life cycle and perform basic cellular functions (the so-called "housekeeping genes"). The other one was made up of 3162 genes that are expressed only at particular stages of development ("developmental genes"). These two groups of genes are so different that we may state that the genome has two types of genetic organization. Different are the timings of their expression, chromatin packaging levels, the composition of activating and deactivating proteins, the sizes of these genes, the lengths of their introns, the organization of the promoter regions of the genes, the locations of origin recognition complexes (ORCs), and DNA replication timings.

Transcription factor B-H2 regulates CYP9J34 expression conveying deltamethrin resistance in Culex pipiens pallens

BACKGROUND

Mosquitoes are vectors of various diseases, posing significant health threats worldwide. Chemical pesticides, particularly pyrethroids like deltamethrin, are commonly used for mosquito control, but the emergence of resistant mosquito populations has become a concern. In the deltamethrin-resistant (DR) strain of Culex pipiens pallens, the highly expressed cytochrome P450 9 J34 (CYP9J34) gene is believed to play a role in resistance, yet the underlying mechanism remains unclear.

RESULTS

Quantitative polymerase chain reaction with reverse transcription (qRT-PCR) analysis revealed that the expression of CYP9J34 was 14.6-fold higher in DR strains than in deltamethrin-susceptible (DS) strains. The recombinant production of CYP9J34 protein of Cx. pipiens pallens showed that the protein could directly metabolize deltamethrin, yielding the major metabolite 4'-OH deltamethrin. Through dual luciferase reporter assays and RNA interference, the transcription factor homeobox protein B-H2-like (B-H2) was identified to modulate the expression of the CYP9J34 gene, contributing to mosquito resistance to deltamethrin.

CONCLUSIONS

Our findings demonstrate that the CYP9J34 protein could directly degrade deltamethrin, and the transcription factor B-H2 could regulate CYP9J34 expression, influencing the resistance of mosquitoes to deltamethrin. © 2023 Society of Chemical Industry.

SWI/SNF-dependent genes are defined by their chromatin landscape

SWI/SNF complexes are evolutionarily conserved, ATP-dependent chromatin remodeling machines. Here, we characterize the features of SWI/SNF-dependent genes using BRM014, an inhibitor of the ATPase activity of the complexes. We find that SWI/SNF activity is required to maintain chromatin accessibility and nucleosome occupancy for most enhancers but not for most promoters. SWI/SNF activity is needed for expression of genes with low to medium levels of expression that have promoters with (1) low chromatin accessibility, (2) low levels of active histone marks, (3) high H3K4me1/H3K4me3 ratio, (4) low nucleosomal phasing, and (5) enrichment in TATA-box motifs. These promoters are mostly occupied by the canonical Brahma-related gene 1/Brahma-associated factor (BAF) complex. These genes are surrounded by SWI/SNF-dependent enhancers and mainly encode signal transduction, developmental, and cell identity genes (with almost no housekeeping genes). Machine-learning models trained with different chromatin characteristics of promoters and their surrounding regulatory regions indicate that the chromatin landscape is a determinant for establishing SWI/SNF dependency.

Identification of teleost tnnc1a enhancers for specific pan-cardiac transcription

Troponin C regulates muscle contraction by forming the troponin complex with troponin I and troponin T. Different muscle types express different troponin C genes. The mechanisms of such differential transcription are not fully understood. The Zebrafish tnnc1a gene is restrictively expressed in cardiac muscles. We here identify the enhancers and promoters of the zebrafish and medaka tnnc1a genes, including intronic enhancers in zebrafish and medaka and an upstream enhancer in the medaka. The intronic and upstream enhancers are likely functionally redundant. The GFP transgenic reporter driven by these enhancers is expressed more strongly in the ventricle than in the atrium, recapitulating the expression pattern of the endogenous zebrafish tnnc1a gene. Our study identifies a new set of enhancers for cardiac-specific transgenic expression in zebrafish. These enhancers can serve as tools for future identification of transcription factor networks that drive cardiac-specific gene transcription.

Dynamics of Chromatin Opening across Larval Development in the Urochordate Ascidian Ciona savignyi

Ascidian larvae undergo tail elongation and notochord lumenogenesis, making them an ideal model for investigating tissue morphogenesis in embryogenesis. The cellular and mechanical mechanisms of these processes have been studied; however, the underlying molecular regulatory mechanism remains to be elucidated. In this study, assays for transposase-accessible chromatin using sequencing (ATAC-seq) and RNA sequencing (RNA-seq) were applied to investigate potential regulators of the development of ascidian Ciona savignyi larvae. Our results revealed 351 and 138 differentially accessible region genes through comparisons of ATAC-seq data between stages 21 and 24 and between stages 24 and 25, respectively. A joint analysis of RNA-seq and ATAC-seq data revealed a correlation between chromatin accessibility and gene transcription. We further verified the tissue expression patterns of 12 different genes. Among them, Cs-matrix metalloproteinase 24 (MMP24) and Cs-krüppel-like factor 5 (KLF5) were highly expressed in notochord cells. Functional assay results demonstrated that both genes are necessary for notochord lumen formation and expansion. Finally, we performed motif enrichment analysis of the differentially accessible regions in different tailbud stages and summarized the potential roles of these motif-bearing transcription factors in larval development. Overall, our study found a correlation between gene expression and chromatin accessibility and provided a vital resource for understanding the mechanisms of the development of ascidian embryos.

Participation of the TBP-associated factors (TAFs) in cell differentiation

The understanding of the mechanisms that regulate gene expression to establish differentiation programs and determine cell lineages, is one of the major challenges in Developmental Biology. Besides the participation of tissue-specific transcription factors and epigenetic processes, the role of general transcription factors has been ignored. Only in recent years, there have been scarce studies that address this issue. Here, we review the studies on the biological activity of some TATA-box binding protein (TBP)-associated factors (TAFs) during the proliferation of stem/progenitor cells and their involvement in cell differentiation. Particularly, the accumulated evidence suggests that TAF4, TAF4b, TAF7L, TAF8, TAF9, and TAF10, among others, participate in nervous system development, adipogenesis, myogenesis, and epidermal differentiation; while TAF1, TAF7, TAF15 may be involved in the regulation of stem cell proliferative abilities and cell cycle progression. On the other hand, evidence suggests that TBP variants such as TBPL1 and TBPL2 might be regulating some developmental processes such as germ cell maturation and differentiation, myogenesis, or ventral specification during development. Our analysis shows that it is necessary to study in greater depth the biological function of these factors and its participation in the assembly of specific transcription complexes that contribute to the differential gene expression that gives rise to the great diversity of cell types existing in an organism. The understanding of TAFs' regulation might lead to the development of new therapies for patients which suffer from mutations, alterations, and dysregulation of these essential elements of the transcriptional machinery.

Cis-Regulatory Elements in Mammals

In cis-regulatory elements, enhancers and promoters with complex molecular interactions are used to coordinate gene transcription through physical proximity and chemical modifications. These processes subsequently influence the phenotypic characteristics of an organism. An in-depth exploration of enhancers and promoters can substantially enhance our understanding of gene regulatory networks, shedding new light on mammalian development, evolution and disease pathways. In this review, we provide a comprehensive overview of the intrinsic structural attributes, detection methodologies as well as the operational mechanisms of enhancers and promoters, coupled with the relevant novel and innovative investigative techniques used to explore their actions. We further elucidated the state-of-the-art research on the roles of enhancers and promoters in the realms of mammalian development, evolution and disease, and we conclude with forward-looking insights into prospective research avenues.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

Regulatory Mechanisms for Transcriptional Bursting Revealed by an Event-Based Model

Gene transcription often occurs in discrete bursts, and it can be difficult to deduce the underlying regulatory mechanisms for transcriptional bursting with limited experimental data. Here, we categorize numerous states of single eukaryotic genes and identify 6 essential transcriptional events, each comprising a series of state transitions; transcriptional bursting is characterized as a sequence of 4 events, capable of being organized in various configurations, in addition to the beginning and ending events. By associating transcriptional kinetics with mean durations and recurrence probabilities of the events, we unravel how transcriptional bursting is modulated by various regulators including transcription factors. Through analytical derivation and numerical simulation, this study reveals key state transitions contributing to transcriptional sensitivity and specificity, typical characteristics of burst profiles, global constraints on intrinsic transcriptional noise, major regulatory modes in individual genes and across the genome, and requirements for fast gene induction upon stimulation. It is illustrated how biochemical reactions on different time scales are modulated to separately shape the durations and ordering of the events. Our results suggest that transcriptional patterns are essentially controlled by a shared set of transcriptional events occurring under specific promoter architectures and regulatory modes, the number of which is actually limited.

Understanding blood development and leukemia using sequencing-based technologies and human cell systems

Our current understanding of human hematopoiesis has undergone significant transformation throughout the years, challenging conventional views. The evolution of high-throughput technologies has enabled the accumulation of diverse data types, offering new avenues for investigating key regulatory processes in blood cell production and disease. In this review, we will explore the opportunities presented by these advancements for unraveling the molecular mechanisms underlying normal and abnormal hematopoiesis. Specifically, we will focus on the importance of enhancer-associated regulatory networks and highlight the crucial role of enhancer-derived transcription regulation. Additionally, we will discuss the unprecedented power of single-cell methods and the progression in using in vitro human blood differentiation system, in particular induced pluripotent stem cell models, in dissecting hematopoietic processes. Furthermore, we will explore the potential of ever more nuanced patient profiling to allow precision medicine approaches. Ultimately, we advocate for a multiparameter, regulatory network-based approach for providing a more holistic understanding of normal hematopoiesis and blood disorders.

Transcription factor IID parks and drives preinitiation complexes at sharp or broad promoters

Core promoters are sites where transcriptional regulatory inputs of a gene are integrated to direct the assembly of the preinitiation complex (PIC) and RNA polymerase II (Pol II) transcription output. Until now, core promoter functions have been investigated by distinct methods, including Pol II transcription initiation site mappings and structural characterization of PICs on distinct promoters. Here, we bring together these previously unconnected observations and hypothesize how, on metazoan TATA promoters, the precisely structured building up of transcription factor (TF) IID-based PICs results in sharp transcription start site (TSS) selection; or, in contrast, how the less strictly controlled positioning of the TATA-less promoter DNA relative to TFIID-core PIC components results in alternative broad TSS selections by Pol II.

Transposon accumulation at xenobiotic gene family loci in aphids

The evolution of resistance is a major challenge for the sustainable control of pests and pathogens. Thus, a deeper understanding of the evolutionary and genomic mechanisms underpinning resistance evolution is required to safeguard health and food production. Several studies have implicated transposable elements (TEs) in xenobiotic-resistance evolution in insects. However, analyses are generally restricted to one insect species and/or one or a few xenobiotic gene families (XGFs). We examine evidence for TE accumulation at XGFs by performing a comparative genomic analysis across 20 aphid genomes, considering major subsets of XGFs involved in metabolic resistance to insecticides: cytochrome P450s, glutathione S-transferases, esterases, UDP-glucuronosyltransferases, and ABC transporters. We find that TEs are significantly enriched at XGFs compared with other genes. XGFs show similar levels of TE enrichment to those of housekeeping genes. But unlike housekeeping genes, XGFs are not constitutively expressed in germline cells, supporting the selective enrichment of TEs at XGFs rather than enrichment owing to chromatin availability. Hotspots of extreme TE enrichment occur around certain XGFs. We find, in aphids of agricultural importance, particular enrichment of TEs around cytochrome P450 genes with known functions in the detoxification of synthetic insecticides. Our results provide evidence supporting a general role for TEs as a source of genomic variation at host XGFs and highlight the existence of considerable variability in TE content across XGFs and host species. These findings show the need for detailed functional verification analyses to clarify the significance of individual TE insertions and elucidate underlying mechanisms at TE-XGF hotspots.

Different Protein Groups Involved in Transcription Regulation in Development and Housekeeping Genes in Drosophila

Antibodies to histone modifications and an insulator protein involved in the processes of transcription initiation and elongation are mapped in Drosophila polytene chromosomes. The CHRIZ protein (chromatin insulator) and H3K36me3 histone modification (RNA elongation) are detected only in the localization of housekeeping genes (interbands and gray bands of polytene chromosomes) and never in the regions of developmental genes (black bands and large puffs arising from them). Antibodies to H3S10P histone modification, which is associated with the initial elongation of the RNA strand during transcription, are found exclusively in small puffs, but not in housekeeping gene localization sites or large ecdysone-induced puffs, where housekeeping genes are localized. Antibodies to H4R3me2 histone modification (a co-repressor of the ecdysone receptor) are detected only in large ecdysone-induced puffs.

Combining TSS-MPRA and sensitive TSS profile dissimilarity scoring to study the sequence determinants of transcription initiation

Cis-regulatory elements (CREs) can be classified by the shapes of their transcription start site (TSS) profiles, which are indicative of distinct regulatory mechanisms. Massively parallel reporter assays (MPRAs) are increasingly being used to study CRE regulatory mechanisms, yet the degree to which MPRAs replicate individual endogenous TSS profiles has not been determined. Here, we present a new low-input MPRA protocol (TSS-MPRA) that enables measuring TSS profiles of episomal reporters as well as after lentiviral reporter chromatinization. To sensitively compare MPRA and endogenous TSS profiles, we developed a novel dissimilarity scoring algorithm (WIP score) that outperforms the frequently used earth mover's distance on experimental data. Using TSS-MPRA and WIP scoring on 500 unique reporter inserts, we found that short (153 bp) MPRA promoter inserts replicate the endogenous TSS patterns of ∼60% of promoters. Lentiviral reporter chromatinization did not improve fidelity of TSS-MPRA initiation patterns, and increasing insert size frequently led to activation of extraneous TSS in the MPRA that are not active in vivo. We discuss the implications of our findings, which highlight important caveats when using MPRAs to study transcription mechanisms. Finally, we illustrate how TSS-MPRA and WIP scoring can provide novel insights into the impact of transcription factor motif mutations and genetic variants on TSS patterns and transcription levels.

Comparison of antiviral responses in two bat species reveals conserved and divergent innate immune pathways

Bats host a range of disease-causing viruses without displaying clinical symptoms. The mechanisms behind this are a continuous source of interest. Here, we studied the antiviral response in the Egyptian fruit bat and Kuhl's pipistrelle, representing two subordinal clades. We profiled the antiviral response in fibroblasts using RNA sequencing and compared bat with primate and rodent responses. Both bats upregulate similar genes; however, a subset of these genes is transcriptionally divergent between them. These divergent genes also evolve rapidly in sequence, have specific promoter architectures, and are associated with programs underlying tolerance and resistance. Finally, we characterized antiviral genes that expanded in bats, with duplicates diverging in sequence and expression. Our study reveals a largely conserved antiviral program across bats and points to a set of genes that rapidly evolve through multiple mechanisms. These can contribute to bat adaptation to viral infection and provide directions to understanding the mechanisms behind it.

Neuron-specific chromatin disruption at CpG islands and aging-related regions in Kabuki syndrome mice

Many Mendelian developmental disorders caused by coding variants in epigenetic regulators have now been discovered. Epigenetic regulators are broadly expressed, and each of these disorders typically exhibits phenotypic manifestations from many different organ systems. An open question is whether the chromatin disruption - the root of the pathogenesis - is similar in the different disease-relevant cell types. This is possible in principle, since all these cell-types are subject to effects from the same causative gene, that has the same kind of function (e.g. methylates histones) and is disrupted by the same germline variant. We focus on mouse models for Kabuki syndrome types 1 and 2, and find that the chromatin accessibility abnormalities in neurons are mostly distinct from those in B or T cells. This is not because the neuronal abnormalities occur at regulatory elements that are only active in neurons. Neurons, but not B or T cells, show preferential chromatin disruption at CpG islands and at regulatory elements linked to aging. A sensitive analysis reveals that the regions disrupted in B/T cells do exhibit chromatin accessibility changes in neurons, but these are very subtle and of uncertain functional significance. Finally, we are able to identify a small set of regulatory elements disrupted in all three cell types. Our findings reveal the cellular-context-specific effect of variants in epigenetic regulators, and suggest that blood-derived "episignatures" may not be well-suited for understanding the mechanistic basis of neurodevelopment in Mendelian disorders of the epigenetic machinery.

Three-dimensional genome structure and function

Linear DNA undergoes a series of compression and folding events, forming various three-dimensional (3D) structural units in mammalian cells, including chromosomal territory, compartment, topologically associating domain, and chromatin loop. These structures play crucial roles in regulating gene expression, cell differentiation, and disease progression. Deciphering the principles underlying 3D genome folding and the molecular mechanisms governing cell fate determination remains a challenge. With advancements in high-throughput sequencing and imaging techniques, the hierarchical organization and functional roles of higher-order chromatin structures have been gradually illuminated. This review systematically discussed the structural hierarchy of the 3D genome, the effects and mechanisms of cis-regulatory elements interaction in the 3D genome for regulating spatiotemporally specific gene expression, the roles and mechanisms of dynamic changes in 3D chromatin conformation during embryonic development, and the pathological mechanisms of diseases such as congenital developmental abnormalities and cancer, which are attributed to alterations in 3D genome organization and aberrations in key structural proteins. Finally, prospects were made for the research about 3D genome structure, function, and genetic intervention, and the roles in disease development, prevention, and treatment, which may offer some clues for precise diagnosis and treatment of related diseases.

Multi-landmark alignment of genomic signals reveals conserved expression patterns across transcription start sites

The prevalent one-dimensional alignment of genomic signals to a reference landmark is a cornerstone of current methods to study transcription and its DNA-dependent processes but it is prone to mask potential relations among multiple DNA elements. We developed a systematic approach to align genomic signals to multiple locations simultaneously by expanding the dimensionality of the genomic-coordinate space. We analyzed transcription in human and uncovered a complex dependence on the relative position of neighboring transcription start sites (TSSs) that is consistently conserved among cell types. The dependence ranges from enhancement to suppression of transcription depending on the relative distances to the TSSs, their intragenic position, and the transcriptional activity of the gene. Our results reveal a conserved hierarchy of alternative TSS usage within a previously unrecognized level of genomic organization and provide a general methodology to analyze complex functional relationships among multiple types of DNA elements.

Pre- and Post-Transcriptional Control of HBV Gene Expression: The Road Traveled towards the New Paradigm of HBx, Its Isoforms, and Their Diverse Functions

Hepatitis B virus (HBV) is an enveloped DNA human virus belonging to the Hepadnaviridae family. Perhaps its main distinguishable characteristic is the replication of its genome through a reverse transcription process. The HBV circular genome encodes only four overlapping reading frames, encoding for the main canonical proteins named core, P, surface, and X (or HBx protein). However, pre- and post-transcriptional gene regulation diversifies the full HBV proteome into diverse isoform proteins. In line with this, hepatitis B virus X protein (HBx) is a viral multifunctional and regulatory protein of 16.5 kDa, whose canonical reading frame presents two phylogenetically conserved internal in-frame translational initiation codons, and which results as well in the expression of two divergent N-terminal smaller isoforms of 8.6 and 5.8 kDa, during translation. The canonical HBx, as well as the smaller isoform proteins, displays different roles during viral replication and subcellular localizations. In this article, we reviewed the different mechanisms of pre- and post-transcriptional regulation of protein expression that take place during viral replication. We also investigated all the past and recent evidence about HBV HBx gene regulation and its divergent N-terminal isoform proteins. Evidence has been collected for over 30 years. The accumulated evidence simply strengthens the concept of a new paradigm of the canonical HBx, and its smaller divergent N-terminal isoform proteins, not only during viral replication, but also throughout cell pathogenesis.

Mutational scans reveal differential evolvability of Drosophila promoters and enhancers

Rapid enhancer and slow promoter evolution have been demonstrated through comparative genomics. However, it is not clear how this information is encoded genetically and if this can be used to place evolution in a predictive context. Part of the challenge is that our understanding of the potential for regulatory evolution is biased primarily toward natural variation or limited experimental perturbations. Here, to explore the evolutionary capacity of promoter variation, we surveyed an unbiased mutation library for three promoters in Drosophila melanogaster. We found that mutations in promoters had limited to no effect on spatial patterns of gene expression. Compared to developmental enhancers, promoters are more robust to mutations and have more access to mutations that can increase gene expression, suggesting that their low activity might be a result of selection. Consistent with these observations, increasing the promoter activity at the endogenous locus of shavenbaby led to increased transcription yet limited phenotypic changes. Taken together, developmental promoters may encode robust transcriptional outputs allowing evolvability through the integration of diverse developmental enhancers.

This article is part of the theme issue ‘Interdisciplinary approaches to predicting evolutionary biology’.

Alternative promoters in CpG depleted regions are prevalently associated with epigenetic misregulation of liver cancer transcriptomes

Transcriptional regulation is commonly governed by alternative promoters. However, the regulatory architecture in alternative and reference promoters, and how they differ, remains elusive. In 100 CAGE-seq libraries from hepatocellular carcinoma patients, here we annotate 4083 alternative promoters in 2926 multi-promoter genes, which are largely undetected in normal livers. These genes are enriched in oncogenic processes and predominantly show association with overall survival. Alternative promoters are narrow nucleosome depleted regions, CpG island depleted, and enriched for tissue-specific transcription factors. Globally tumors lose DNA methylation. We show hierarchical retention of intragenic DNA methylation with CG-poor regions rapidly losing methylation, while CG-rich regions retain it, a process mediated by differential SETD2, H3K36me3, DNMT3B, and TET1 binding. This mechanism is validated in SETD2 knockdown cells and SETD2-mutated patients. Selective DNA methylation loss in CG-poor regions makes the chromatin accessible for alternative transcription. We show alternative promoters can control tumor transcriptomes and their regulatory architecture.

Promoter evolution of mammalian gene duplicates

BACKGROUND

Gene duplication is thought to be a central process in evolution to gain new functions. The factors that dictate gene retention following duplication as well paralog gene divergence in sequence, expression and function have been extensively studied. However, relatively little is known about the evolution of promoter regions of gene duplicates and how they influence gene duplicate divergence. Here, we focus on promoters of paralog genes, comparing their similarity in sequence, in the sets of transcription factors (TFs) that bind them, and in their overall promoter architecture.

RESULTS

We observe that promoters of recent duplications display higher sequence similarity between them and that sequence similarity rapidly declines between promoters of more ancient paralogs. In contrast, similarity in cis-regulation, as measured by the set of TFs that bind promoters of both paralogs, does not simply decrease with time from duplication and is instead related to promoter architecture-paralogs with CpG Islands (CGIs) in their promoters share a greater fraction of TFs, while CGI-less paralogs are more divergent in their TF binding set. Focusing on recent duplication events and partitioning them by their duplication mechanism enables us to uncover promoter properties associated with gene retention, as well as to characterize the evolution of promoters of newly born genes: In recent retrotransposition-mediated duplications, we observe asymmetry in cis-regulation of paralog pairs: Retrocopy genes are lowly expressed and their promoters are bound by fewer TFs and are depleted of CGIs, in comparison with the original gene copy. Furthermore, looking at recent segmental duplication regions in primates enable us to compare successful retentions versus loss of duplicates, showing that duplicate retention is associated with fewer TFs and with CGI-less promoter architecture.

CONCLUSIONS

In this work, we profiled promoters of gene duplicates and their inter-paralog divergence. We also studied how their characteristics are associated with duplication time and duplication mechanism, as well as with the fate of these duplicates. These results underline the importance of cis-regulatory mechanisms in shaping the evolution of new genes and their fate following duplication.

Neutral Models of De Novo Gene Emergence Suggest that Gene Evolution has a Preferred Trajectory

New protein coding genes can emerge from genomic regions that previously did not contain any genes, via a process called de novo gene emergence. To synthesize a protein, DNA must be transcribed as well as translated. Both processes need certain DNA sequence features. Stable transcription requires promoters and a polyadenylation signal, while translation requires at least an open reading frame. We develop mathematical models based on mutation probabilities, and the assumption of neutral evolution, to find out how quickly genes emerge and are lost. We also investigate the effect of the order by which DNA features evolve, and if sequence composition is biased by mutation rate. We rationalize how genes are lost much more rapidly than they emerge, and how they preferentially arise in regions that are already transcribed. Our study not only answers some fundamental questions on the topic of de novo emergence but also provides a modeling framework for future studies.

Evolution of homology: From archetype towards a holistic concept of cell type

The concept of homology lies in the heart of comparative biological science. The distinction between homology as structure and analogy as function has shaped the evolutionary paradigm for a century and formed the axis of comparative anatomy and embryology, which accept the identity of structure as a ground measure of relatedness. The advent of single-cell genomics overturned the classical view of cell homology by establishing a backbone regulatory identity of cell types, the basic biological units bridging the molecular and phenotypic dimensions, to reveal that the cell is the most flexible unit of living matter and that many approaches of classical biology need to be revised to understand evolution and diversity at the cellular level. The emerging theory of cell types explicitly decouples cell identity from phenotype, essentially allowing for the divergence of evolutionarily related morphotypes beyond recognition, as well as it decouples ontogenetic cell lineage from cell-type phylogeny, whereby explicating that cell types can share common descent regardless of their structure, function or developmental origin. The article succinctly summarizes current progress and opinion in this field and formulates a more generalistic view of biological cell types as avatars, transient or terminal cell states deployed in a continuum of states by the developmental programme of one and the same omnipotent cell, capable of changing or combining identities with distinct evolutionary histories or inventing ad hoc identities that never existed in evolution or development. It highlights how the new logic grounded in the regulatory nature of cell identity transforms the concepts of cell homology and phenotypic stability, suggesting that cellular evolution is inherently and massively network-like, with one-to-one homologies being rather uncommon and restricted to shallower levels of the animal tree of life.

Fish species lifespan prediction from promoter cytosine-phosphate-guanine density

Lifespan is a key attribute of a species' life cycle and varies extensively among major lineages of animals. In fish, lifespan varies by several orders of magnitude, with reported values ranging from less than 1 year to approximately 400 years. Lifespan information is particularly useful for species management, as it can be used to estimate invasion potential, extinction risk and sustainable harvest rates. Despite its utility, lifespan is unknown for most fish species. This is due to the difficulties associated with accurately identifying the oldest individual(s) of a given species, and/or deriving lifespan estimates that are representative for an entire species. Recently it has been shown that CpG density in gene promoter regions can be used to predict lifespan in mammals and other vertebrates, with variable accuracy across taxa. To improve accuracy of lifespan prediction in a non-mammalian vertebrate group, here we develop a fish-specific genomic lifespan predictor. Our new model includes more than eight times the number of fish species included in the previous vertebrate model (n = 442) and uses fish-specific gene promoters as reference sequences. The model predicts fish lifespan from genomic CpG density alone (measured as CpG observed/expected ratio), explaining 64% of the variance between known and predicted lifespans. The predictions are highly robust to variation in genome quality and are applicable to all classes of fish; a taxonomically diverse and speciose group. The results demonstrate the value of promoter CpG density as a universal predictor of fish lifespan that can applied where empirical data are unavailable, or impracticable to obtain.

The miR-430 locus with extreme promoter density forms a transcription body during the minor wave of zygotic genome activation

In anamniote embryos, the major wave of zygotic genome activation starts during the mid-blastula transition. However, some genes escape global genome repression, are activated substantially earlier, and contribute to the minor wave of genome activation. The mechanisms underlying the minor wave of genome activation are little understood. We explored the genomic organization and cis-regulatory mechanisms of a transcription body, in which the minor wave of genome activation is first detected in zebrafish. We identified the miR-430 cluster as having excessive copy number and the highest density of Pol-II-transcribed promoters in the genome, and this is required for forming the transcription body. However, this transcription body is not essential for, nor does it encompasse, minor wave transcription globally. Instead, distinct minor-wave-specific promoter architecture suggests that promoter-autonomous mechanisms regulate the minor wave of genome activation. The minor-wave-specific features also suggest distinct transcription initiation mechanisms between the minor and major waves of genome activation.

The Role of PARP1 and PAR in ATP-Independent Nucleosome Reorganisation during the DNA Damage Response

The functioning of the eukaryotic cell genome is mediated by sophisticated protein-nucleic-acid complexes, whose minimal structural unit is the nucleosome. After the damage to genomic DNA, repair proteins need to gain access directly to the lesion; therefore, the initiation of the DNA damage response inevitably leads to local chromatin reorganisation. This review focuses on the possible involvement of PARP1, as well as proteins acting nucleosome compaction, linker histone H1 and non-histone chromatin protein HMGB1. The polymer of ADP-ribose is considered the main regulator during the development of the DNA damage response and in the course of assembly of the correct repair complex.

Loss of the Ash2l subunit of histone H3K4 methyltransferase complexes reduces chromatin accessibility at promoters

Changes in gene expression programs are intimately linked to cell fate decisions. Post-translational modifications of core histones contribute to control gene expression. Methylation of lysine 4 of histone H3 (H3K4) correlates with active promoters and gene transcription. This modification is catalyzed by KMT2 methyltransferases, which require interaction with 4 core subunits, WDR5, RBBP5, ASH2L and DPY30, for catalytic activity. Ash2l is necessary for organismal development and for tissue homeostasis. In mouse embryo fibroblasts (MEFs), Ash2l loss results in gene repression, provoking a senescence phenotype. We now find that upon knockout of Ash2l both H3K4 mono- and tri-methylation (H3K4me1 and me3, respectively) were deregulated. In particular, loss of H3K4me3 at promoters correlated with gene repression, especially at CpG island promoters. Ash2l loss resulted in increased loading of histone H3 and reduced chromatin accessibility at promoters, accompanied by an increase of repressing and a decrease of activating histone marks. Moreover, we observed altered binding of CTCF upon Ash2l loss. Lost and gained binding was noticed at promoter-associated and intergenic sites, respectively. Thus, Ash2l loss and reduction of H3K4me3 correlate with altered chromatin accessibility and transcription factor binding. These findings contribute to a more detailed understanding of mechanistic consequences of H3K4me3 loss and associated repression of gene transcription and thus of the observed cellular consequences.

Identification and Functional Characterization of the Transcription Factors AhR/ARNT in Dendroctonus armandi

The aryl hydrocarbon receptor (AhR) and aryl hydrocarbon receptor nuclear translocator (ARNT) belong to the bHLH-PAS (basic Helix-Loop-Helix-Period/ARNT/Single-minded) family of transcription factors, which participate in the sensing and transmitting stimuli of exogenous and endogenous chemical substances, and subsequently activates genes transcription involved in various detoxification and physiological functions. However, they have not been identified in Dendroctonus armandi, and their roles in the detoxification metabolism are unclear. In the present study, AhR and ARNT of D. armandi were characterized. Spatiotemporal expression profiling indicated that DaAhR and DaARNT were highly expressed in the adult and larval stages of D. armandi and mainly expressed in the midgut and Malpighian tubules of adults. Additionally, the expression of DaAhR and DaARNT significantly increased after exposure to (-)-𝛽-pinene, (+)-3-carene, and (±)-limonene. Silencing DaAhR and DaARNT increased the susceptibility of D. armandi to (-)-𝛽-pinene, (+)-3-carene, and (±)-limonene, and the activities of detoxification enzyme were also remarkably reduced. Moreover, DaCYP6DF1 and DaGSTs2 were significantly down-regulated after injections of dsAhR and dsARNT in the male and female adults, with the expression of DaCYP6DF1 decreasing by higher than 70%. The present study revealed that the transcription factors AhR and ARNT of D. armandi were induced by terpenoids and participated in the regulation of DaCYP6DF1 expression, which was associated with D. armandi's susceptibility to (-)-𝛽-pinene and (±)-limonene. These results may provide a theoretical basis for the integrated control of D. armandi and improve our comprehension of insect toxicology.

EPIHC: Improving Enhancer-Promoter Interaction Prediction by Using Hybrid Features and Communicative Learning

Enhancer-promoter interactions (EPIs) regulate the expression of specific genes in cells, which help facilitate understanding of gene regulation, cell differentiation and disease mechanisms. EPI identification approaches through wet experiments are often costly and time-consuming, leading to the design of high-efficiency computational methods is in demand. In this paper, we propose a deep neural network-based method named EPIHC to predict Enhancer-Promoter Interactions with Hybrid features and Communicative learning. EPIHC extracts enhancer and promoter sequence-derived features using convolutional neural networks (CNN), and then we design a communicative learning module to capture the communicative information between enhancer and promoter sequences. Besides, EPIHC takes the genomic features of enhancers and promoters into account, incorporating with the sequence-derived features to predict EPIs. The computational experiments show that EPIHC outperforms the existing state-of-the-art EPI prediction methods on the benchmark datasets and chromosome-split datasets, and the study reveals that the communicative learning module can bring explicit information about EPIs, which is ignored by CNN, and provide explainability about EPIs to some degree. Moreover, we consider two strategies to improve the performances of EPIHC in the cross-cell line prediction, and experimental results show that EPIHC constructed on some cell lines can exhibit good performances for other cell lines. The codes and data are available at https://github.com/BioMedicalBigDataMiningLab/EPIHC.

iPromoter-CLA: Identifying promoters and their strength by deep capsule networks with bidirectional long short-term memory

BACKGROUND AND OBJECTIVE

The promoter is a fragment of DNA and a specific sequence with transcriptional regulation function in DNA. Promoters are located upstream at the transcription start site, which is used to initiate downstream gene expression. So far, promoter identification is mainly achieved by biological methods, which often require more effort. It has become a more effective classification and prediction method to identify promoter types through computational methods.

METHODS

In this study, we proposed a new capsule network and recurrent neural network hybrid model to identify promoters and predict their strength. Firstly, we used one-hot to encode DNA sequence. Secondly, we used three one-dimensional convolutional layers, a one-dimensional convolutional capsule layer and digit capsule layer to learn local features. Thirdly, a bidirectional long short-time memory was utilized to extract global features. Finally, we adopted the self-attention mechanism to improve the contribution of relatively important features, which further enhances the performance of the model.

RESULTS

Our model attains a cross-validation accuracy of 86% and 73.46% in prokaryotic promoter recognition and their strength prediction, which showcases a better performance compared with the existing approaches in both the first layer promoter identification and the second layer promoter's strength prediction.

CONCLUSIONS

our model not only combines convolutional neural network and capsule layer but also uses a self-attention mechanism to better capture hidden information features from the perspective of sequence. Thus, we hope that our model can be widely applied to other components.

Exploiting Genomic Features to Improve the Prediction of Transcription Factor-Binding Sites in Plants

The identification of transcription factor (TF) target genes is central in biology. A popular approach is based on the location by pattern matching of potential cis-regulatory elements (CREs). During the last few years, tools integrating next-generation sequencing data have been developed to improve the performance of pattern matching. However, such tools have not yet been comprehensively evaluated in plants. Hence, we developed a new streamlined method aiming at predicting CREs and target genes of plant TFs in specific organs or conditions. Our approach implements a supervised machine learning strategy, which allows decision rule models to be learnt using TF ChIP-chip/seq experimental data. Different layers of genomic features were integrated in predictive models: the position on the gene, the DNA sequence conservation, the chromatin state and various CRE footprints. Among the tested features, the chromatin features were crucial for improving the accuracy of the method. Furthermore, we evaluated the transferability of predictive models across TFs, organs and species. Finally, we validated our method by correctly inferring the target genes of key TFs controlling metabolite biosynthesis at the organ level in Arabidopsis. We developed a tool-Wimtrap-to reproduce our approach in plant species and conditions/organs for which ChIP-chip/seq data are available. Wimtrap is a user-friendly R package that supports an R Shiny web interface and is provided with pre-built models that can be used to quickly get predictions of CREs and TF gene targets in different organs or conditions in Arabidopsis thaliana, Solanum lycopersicum, Oryza sativa and Zea mays.

Performance of abiotic stress-inducible synthetic promoters in genetically engineered hybrid poplar (Populus tremula × Populus alba)

Abiotic stresses can cause significant damage to plants. For sustainable bioenergy crop production, it is critical to generate resistant crops to such stress. Engineering promoters to control the precise expression of stress resistance genes is a very effective way to address the problem. Here we developed stably transformed Populus tremula × Populus alba hybrid poplar (INRA 717-1B4) containing one-of-six synthetic drought stress-inducible promoters (SDs; SD9-1, SD9-2, SD9-3, SD13-1, SD18-1, and SD18-3) identified previously by transient transformation assays. We screened green fluorescent protein (GFP) induction in poplar under osmotic stress conditions. Of six transgenic lines containing synthetic promoter, three lines (SD18-1, 9-2, and 9-3) had significant GFP expression in both salt and osmotic stress treatments. Each synthetic promoter employed heptamerized repeats of specific and short cis-regulatory elements (7 repeats of 7-8 bases). To verify whether the repeats of longer sequences can improve osmotic stress responsiveness, a transgenic poplar containing the synthetic promoter of the heptamerized entire SD9 motif (20 bases, containing all partial SD9 motifs) was generated and measured for GFP induction under osmotic stress. The heptamerized entire SD9 motif did not result in higher GFP expression than the shorter promoters consisting of heptamerized SD9-1, 9-2, and 9-3 (partial SD9) motifs. This result indicates that shorter synthetic promoters (~50 bp) can be used for versatile control of gene expression in transgenic poplar. These synthetic promoters will be useful tools to engineer stress-resilient bioenergy tree crops in the future.

Developmental and housekeeping transcriptional programs in Drosophila require distinct chromatin remodelers

Gene transcription is a highly regulated process in all animals. In Drosophila, two major transcriptional programs, housekeeping and developmental, have promoters with distinct regulatory compatibilities and nucleosome organization. However, it remains unclear how the differences in chromatin structure relate to the distinct regulatory properties and which chromatin remodelers are required for these programs. Using rapid degradation of core remodeler subunits in Drosophila melanogaster S2 cells, we demonstrate that developmental gene transcription requires SWI/SNF-type complexes, primarily to maintain distal enhancer accessibility. In contrast, wild-type-level housekeeping gene transcription requires the Iswi and Ino80 remodelers to maintain nucleosome positioning and phasing at promoters. These differential remodeler dependencies relate to different DNA-sequence-intrinsic nucleosome affinities, which favor a default ON state for housekeeping but a default OFF state for developmental gene transcription. Overall, our results demonstrate how different transcription-regulatory strategies are implemented by DNA sequence, chromatin structure, and remodeler activity.

Alternative Transcription Start Site Usage and Functional Implications in Pathogenic Fungi

Pathogenic fungi require delicate gene regulation mechanisms to adapt to diverse living environments and escape host immune systems. Recent advances in sequencing technology have exposed the complexity of the fungal genome, thus allowing the gradual disentanglement of multiple layers of gene expression control. Alternative transcription start site (aTSS) usage, previously reported to be prominent in mammals and to play important roles in physiopathology, is also present in fungi to fine-tune gene expression. Depending on the alteration in their sequences, RNA isoforms arising from aTSSs acquire different characteristics that significantly alter their stability and translational capacity as well as the properties and biologic functions of the resulting proteins. Disrupted control of aTSS usage has been reported to severely impair growth, virulence, and the infectious capacity of pathogenic fungi. Here, we discuss principle concepts, mechanisms, and the functional implication of aTSS usage in fungi.

Regulation of arsenic methylation: identification of the transcriptional region of the human AS3MT gene

The human enzyme As(III) S-adenosylmethionine methyltransferase (AS3MT) catalyzes arsenic biotransformations and is considered to contribute to arsenic-related diseases. AS3MT is expressed in various tissues and cell types including liver, brain, adrenal gland, and peripheral blood mononuclear cells but not in human keratinocytes, urothelial, or brain microvascular endothelial cells. This indicates that AS3MT expression is regulated in a tissue/cell type-specific manner, but the mechanism of transcriptional regulation of expression of the AS3MT gene is not known. In this study, we define the DNA sequence of the core promoter region of the human AS3MT gene. We identify a GC box in the promoter to which the stress-related transcription factor Sp1 binds, indicating involvement of regulatory elements in AS3MT gene expression.