Learning Enhancer-Gene associations from Bulk Transcriptomic and Epigenetic Sequencing Data with STITCHIT.

To reveal gene regulation mechanisms, it is essential to understand the role of regulatory elements, which are possibly distant from gene promoters. Integrative analysis of epigenetic and transcriptomic data can be used to gain insights into gene-expression regulation in specific phenotypes. Here, we discuss STITCHIT, an approach to dissect epigenetic variation in a gene-specific manner across many samples for the identification of regulatory elements without relying on peak calling algorithms. The obtained genomic regions are then further refined using a regularized linear model approach, which can also be used to predict gene expression. We illustrate the use of STITCHIT using H3k27ac ChIP-seq and RNA-seq data from the International Human Epigenome Consortium (IHEC).

Unveiling the Complexity of cis-Regulation Mechanisms in Kinases: A Comprehensive Analysis.

Protein cis-regulatory elements (CREs) are regions that modulate the activity of a protein through intramolecular interactions. Kinases, pivotal enzymes in numerous biological processes, often undergo regulatory control via inhibitory interactions in cis. This study delves into the mechanisms of cis regulation in kinases mediated by CREs, employing a combined structural and sequence analysis. To accomplish this, we curated an extensive dataset of kinases featuring annotated CREs, organized into homolog families through multiple sequence alignments. Key molecular attributes, including disorder and secondary structure content, active and ATP-binding sites, post-translational modifications, and disease-associated mutations, were systematically mapped onto all sequences. Additionally, we explored the potential for conformational changes between active and inactive states. Finally, we explored the presence of these kinases within membraneless organelles and elucidated their functional roles therein. CREs display a continuum of structures, ranging from short disordered stretches to fully folded domains. The adaptability demonstrated by CREs in achieving the common goal of kinase inhibition spans from direct autoinhibitory interaction with the active site within the kinase domain, to CREs binding to an alternative site, inducing allosteric regulation revealing distinct types of inhibitory mechanisms, which we exemplify by archetypical representative systems. While this study provides a systematic approach to comprehend kinase CREs, further experimental investigations are imperative to unravel the complexity within distinct kinase families. The insights gleaned from this research lay the foundation for future studies aiming to decipher the molecular basis of kinase dysregulation, and explore potential therapeutic interventions.

The role of structure in regulatory RNA elements.

Regulatory RNA elements fulfill functions such as translational regulation, control of transcript levels, and regulation of viral genome replication. Trans-acting factors (i.e. RNA-binding proteins) bind the so-called cis elements and confer functionality to the complex. The specificity during protein-RNA complex (RNP) formation often exploits the structural plasticity of RNA. Functional integrity of cis-trans pairs depends on the availability of properly folded RNA elements, and RNA conformational transitions can cause diseases. Knowledge of RNA structure and the conformational space is needed for understanding complex formation and deducing functional effects. However, structure determination of RNAs under in vivo conditions remains challenging. This review provides an overview of structured eukaryotic and viral RNA cis elements and discusses the effect of RNA structural equilibria on RNP formation. We showcase implications of RNA structural changes for diseases, outline strategies for RNA structure-based drug targeting, and summarize the methodological toolbox for deciphering RNA structures.

Regulatory elements in SEM1-DLX5-DLX6 (7q21.3) locus contribute to genetic control of coronal nonsyndromic craniosynostosis and bone density-related traits.

Purpose: The etiopathogenesis of coronal nonsyndromic craniosynostosis (cNCS), a congenital condition defined by premature fusion of 1 or both coronal sutures, remains largely unknown. Methods: We conducted the largest genome-wide association study of cNCS followed by replication, fine mapping, and functional validation of the most significant region using zebrafish animal model. Results: Genome-wide association study identified 6 independent genome-wide-significant risk alleles, 4 on chromosome 7q21.3 SEM1-DLX5-DLX6 locus, and their combination conferred over 7-fold increased risk of cNCS. The top variants were replicated in an independent cohort and showed pleiotropic effects on brain and facial morphology and bone mineral density. Fine mapping of 7q21.3 identified a craniofacial transcriptional enhancer (eDlx36) within the linkage region of the top variant (rs4727341; odds ratio [95% confidence interval], 0.48[0.39-0.59]; P = 1.2E-12) that was located in SEM1 intron and enriched in 4 rare risk variants. In zebrafish, the activity of the transfected human eDlx36 enhancer was observed in the frontonasal prominence and calvaria during skull development and was reduced when the 4 rare risk variants were introduced into the sequence. Conclusion: Our findings support a polygenic nature of cNCS risk and functional role of craniofacial enhancers in cNCS susceptibility with potential broader implications for bone health.

Poly(G)7 box: a functional element of mammalian 18S rRNA involved in translation.

In eukaryotes, the ribosomal small subunit (40S) is composed of 18S rRNA and 33 ribosomal proteins. 18S rRNA has a special secondary structure and is an indispensable part of the translation process. Herein, a special sequence located in mammalian 18S rRNA named Poly(G)7box, which is composed of seven guanines, was found. Poly(G)7 can form a special and stable secondary structure by binding to the translation elongation factor subunit eEF1D and the ribosomal protein RPL32. Poly(G)7box was transfected into cells, and the translation efficiency of cells was inhibited. We believe that Poly(G)7box is an important translation-related functional element located on mammalian 18S rRNA, meanwhile the Poly(G)7 located on mRNA 5' and 3' box does not affect mRNA translation.

Synthetic directed evolution for targeted engineering of plant traits.

Improving crop traits requires genetic diversity, which allows breeders to select advantageous alleles of key genes. In species or loci that lack sufficient genetic diversity, synthetic directed evolution (SDE) can supplement natural variation, thus expanding the possibilities for trait engineering. In this review, we explore recent advances and applications of SDE for crop improvement, highlighting potential targets (coding sequences and cis-regulatory elements) and computational tools to enhance crop resilience and performance across diverse environments. Recent advancements in SDE approaches have streamlined the generation of variants and the selection processes; by leveraging these advanced technologies and principles, we can minimize concerns about host fitness and unintended effects, thus opening promising avenues for effectively enhancing crop traits.

Beyond modular enhancers: new questions in cis-regulatory evolution.

Our understanding of how cis-regulatory elements work has advanced rapidly, outpacing our evolutionary models. In this review, we consider the implications of new mechanistic findings for evolutionary developmental biology. We focus on three different debates: whether evolutionary innovation occurs more often via the modification of old cis-regulatory elements or the emergence of new ones; the extent to which individual elements are specific and autonomous or multifunctional and interdependent; and how the robustness of cis-regulatory architectures influences the rate of trait evolution. These discussions lead us to propose new questions for the evo-devo of cis-regulation.

G-quadruplexes as pivotal components of cis-regulatory elements in the human genome.

BACKGROUND: Cis-regulatory elements (CREs) are crucial for regulating gene expression, and G-quadruplexes (G4s), as prototypal non-canonical DNA structures, may play a role in this regulation. However, the relationship between G4s and CREs, especially with non-promoter-like functional elements, requires further systematic investigation. We aimed to investigate the associations between G4s and human cCREs (candidate CREs) inferred from the Encyclopedia of DNA Elements (ENCODE) data. RESULTS: We found that G4s are prominently enriched in most types of cCREs, especially those with promoter-like signatures (PLS). The co-occurrence of CTCF signals with H3K4me3 or H3K27ac signals strengthens the association between cCREs and G4s. Genetic variants in G4s, particularly within their G-runs, exhibit higher regulatory potential and deleterious effects compared to cCREs. The G-runs within G4s near transcriptional start sites (TSSs) are more evolutionarily constrained compared to G-runs in cCREs, while those far from the TSS are relatively less conserved. The presence of G4s is often linked to a more favorable local chromatin environment for the activation and execution of regulatory function of cCREs, potentially attributable to the formation of G4 secondary structures. Finally, we discovered that G4-associated cCREs exhibit widespread activation in a variety of cancers. CONCLUSIONS: Our study suggests that G4s are integral components of human cis-regulatory elements, extending beyond their potential role in promoters. The G4 primary sequences are associated with the localization of CREs, while the G4 structures are linked to the activation of these elements. Therefore, we propose defining G4s as pivotal regulatory elements in the human genome.

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) in Macrophages.

Cleavage Under Targets and Release Using Nuclease (CUT&RUN) is a method to detect specific interactions between DNA and DNA-associated proteins. It is valuable for the characterization of the binding of transcription factors or co-regulators genome wide. Furthermore, it can be used for epigenetic profiling, chromatin accessibility assessment, and identification of regulatory elements. Compared to the more commonly used chromatin immunoprecipitation (ChIP), CUT&RUN has several advantages including an in situ approach as well as no need for sonication. However, the biggest advantage is the reduced cell amounts that are required for CUT&RUN, which makes it more attractive for experiments with limited cell numbers. In this chapter, we describe a reliable CUT&RUN protocol for macrophages that can be performed within 2 days and includes a library preparation so that the sample can be directly sequenced.

Alfin-like (AL) transcription factor family in Oryza sativa L.: Genome-wide analysis and expression profiling under different stresses.

Oryza sativa L. is the world's most essential and economically important food crop. Climate change and ecological imbalances make rice plants vulnerable to abiotic and biotic stresses, threatening global food security. The Alfin-like (AL) transcription factor family plays a crucial role in plant development and stress responses. This study comprehensively analyzed this gene family and their expression profiles in rice, revealing nine AL genes, classifying them into three distinct groups based on phylogenetic analysis and identifying four segmental duplication events. RNA-seq data analysis revealed high expression levels of OsALs in different tissues, growth stages, and their responsiveness to stresses. RT-qPCR data showed significant expression of OsALs in different abiotic stresses. Identification of potential cis-regulatory elements in promoter regions has also unveiled their involvement. Tertiary structures of the proteins were predicted. These findings would lay the groundwork for future research to reveal their molecular mechanism in stress tolerance and plant development.

Genome-Wide Profiling of Cis-regulatory Elements in Mammalian Skin.

The modulation of cis-regulatory elements (e.g., enhancers and promoters) is a major mechanism by which gene expression can be controlled in a temporal and spatially restricted manner. However, methods for both identifying these elements and inferring their activity are limited and often require a substantial investment of time, money, and resources. Here, using mammalian skin as a model, we demonstrate a streamlined protocol by which these hurdles can be overcome using a novel chromatin profiling technique (CUT&RUN) to map histone modifications genome-wide. This protocol can be used to map the location and activity of putative cis-regulatory elements, providing mechanistic insight into how differential gene expression is controlled in mammalian tissues.

Bioinformatic Prediction and High Throughput In Vivo Screening to Identify Cis-Regulatory Elements for the Development of Algal Synthetic Promoters.

Algae biotechnology holds immense promise for revolutionizing the bioeconomy through the sustainable and scalable production of various bioproducts. However, their development has been hindered by the lack of advanced genetic tools. This study introduces a synthetic biology approach to develop such tools, focusing on the construction and testing of synthetic promoters. By analyzing conserved DNA motifs within the promoter regions of highly expressed genes across six different algal species, we identified cis-regulatory elements (CREs) associated with high transcriptional activity. Combining the algorithms POWRS, STREME, and PhyloGibbs, we predicted 1511 CREs and inserted them into a minimal synthetic promoter sequence in 1, 2, or 3 copies, resulting in 4533 distinct synthetic promoters. These promoters were evaluated in vivo for their capacity to drive the expression of a transgene in a high-throughput manner through next-generation sequencing post antibiotic selection and fluorescence-activated cell sorting. To validate our approach, we sequenced hundreds of transgenic lines showing high levels of GFP expression. Further, we individually tested 14 identified promoters, revealing substantial increases in GFP expression─up to nine times higher than the baseline synthetic promoter, with five matching or even surpassing the performance of the native AR1 promoter. As a result of this study, we identified a catalog of CREs that can now be used to build superior synthetic algal promoters. More importantly, here we present a validated pipeline to generate building blocks for innovative synthetic genetic tools applicable to any algal species with a sequenced genome and transcriptome data set.

Distinct epigenomic landscapes underlie tissue-specific memory T cell differentiation.

The memory CD8+ T cell pool contains phenotypically and transcriptionally heterogeneous subsets with specialized functions and recirculation patterns. Here, we examined the epigenetic landscape of CD8+ T cells isolated from seven non-lymphoid organs across four distinct infection models, alongside their circulating T cell counterparts. Using single-cell transposase-accessible chromatin sequencing (scATAC-seq), we found that tissue-resident memory T (TRM) cells and circulating memory T (TCIRC) cells develop along distinct epigenetic trajectories. We identified organ-specific transcriptional regulators of TRM cell development, including FOSB, FOS, FOSL1, and BACH2, and defined an epigenetic signature common to TRM cells across organs. Finally, we found that although terminal TEX cells share accessible regulatory elements with TRM cells, they are defined by TEX-specific epigenetic features absent from TRM cells. Together, this comprehensive data resource shows that TRM cell development is accompanied by dynamic transcriptome alterations and chromatin accessibility changes that direct tissue-adapted and functionally distinct T cell states.

Machine learning-based classification reveals distinct clusters of non-coding genomic allelic variations associated with Erm-mediated antibiotic resistance.

The erythromycin resistance RNA methyltransferase (erm) confers cross-resistance to all therapeutically important macrolides, lincosamides, and streptogramins (MLS phenotype). The expression of erm is often induced by the macrolide-mediated ribosome stalling in the upstream co-transcribed leader sequence, thereby triggering a conformational switch of the intergenic RNA hairpins to allow the translational initiation of erm. We investigated the evolutionary emergence of the upstream erm regulatory elements and the impact of allelic variation on erm expression and the MLS phenotype. Through systematic profiling of the upstream regulatory sequences across all known erm operons, we observed that specific erm subfamilies, such as ermB and ermC, have independently evolved distinct configurations of small upstream ORFs and palindromic repeats. A population-wide genomic analysis of the upstream ermB regions revealed substantial non-random allelic variation at numerous positions. Utilizing machine learning-based classification coupled with RNA structure modeling, we found that many alleles cooperatively influence the stability of alternative RNA hairpin structures formed by the palindromic repeats, which, in turn, affects the inducibility of ermB expression and MLS phenotypes. Subsequent experimental validation of 11 randomly selected variants demonstrated an impressive 91% accuracy in predicting MLS phenotypes. Furthermore, we uncovered a mixed distribution of MLS-sensitive and MLS-resistant ermB loci within the evolutionary tree, indicating repeated and independent evolution of MLS resistance. Taken together, this study not only elucidates the evolutionary processes driving the emergence and development of MLS resistance but also highlights the potential of using non-coding genomic allele data to predict antibiotic resistance phenotypes. IMPORTANCE: Antibiotic resistance (AR) poses a global health threat as the efficacy of available antibiotics has rapidly eroded due to the widespread transmission of AR genes. Using Erm-dependent MLS resistance as a model, this study highlights the significance of non-coding genomic allelic variations. Through a comprehensive analysis of upstream regulatory elements within the erm family, we elucidated the evolutionary emergence and development of AR mechanisms. Leveraging population-wide machine learning (ML)-based genomic analysis, we transformed substantial non-random allelic variations into discernible clusters of elements, enabling precise prediction of MLS phenotypes from non-coding regions. These findings offer deeper insight into AR evolution and demonstrate the potential of harnessing non-coding genomic allele data for accurately predicting AR phenotypes.

Transcriptomic analyses of differentially expressed human genes, micro RNAs and long-non-coding RNAs in severe, symptomatic and asymptomatic malaria infection.

Malaria transmission and endemicity in Africa remains hugely disproportionate compared to the rest of the world. The complex life cycle of P. falciparum (Pf) between the vertebrate human host and the anopheline vector results in differential expression of genes within and between hosts. An in-depth understanding of Pf interaction with various human genes through regulatory elements will pave way for identification of newer tools in the arsenal for malaria control. Therefore, the regulatory elements (REs) involved in the over- or under-expression of various host immune genes hold the key to elucidating alternative control measures that can be applied for disease surveillance, prompt diagnosis and treatment. We carried out an RNAseq analysis to identify differentially expressed genes and network elucidation of non-coding RNAs and target genes associated with immune response in individuals with different clinical outcomes. Raw RNAseq datasets, retrieved for analyses include individuals with severe (Gambia-20), symptomatic (Burkina Faso-15), asymptomatic (Mali-16) malaria as well as uninfected controls (Tanzania-20; Mali-36). Of the total 107 datasets retrieved, we identified 5534 differentially expressed genes (DEGs) among disease and control groups. A peculiar pattern of DEGs was observed, with individuals presenting with severe/symptomatic malaria having the highest and most diverse upregulated genes, while a reverse phenomenon was recorded among asymptomatic and uninfected individuals. In addition, we identified 141 differentially expressed micro RNA (miRNA), of which 78 and 63 were upregulated and downregulated respectively. Interactome analysis revealed a moderate interaction between DEGs and miRNAs. Of all identified miRNA, five were unique (hsa-mir-32, hsa-mir-25, hsa-mir-221, hsa-mir-29 and hsa-mir-148) because of their connectivity to several genes, including hsa-mir-221 connected to 16 genes. Six-hundred and eight differentially expressed long non coding RNA (lncRNA) were also identified, including SLC7A11, LINC01524 among the upregulated ones. Our study provides important insight into host immune genes undergoing differential expression under different malaria conditions. It also identified unique miRNAs and lncRNAs that modify and/or regulate the expression of various immune genes. These regulatory elements we surmise, have the potential to serve a diagnostic purpose in discriminating between individuals with severe/symptomatic malaria and those with asymptomatic infection or uninfected, following further clinical validation from field isolates.

Epigenomic identification of vernalization cis-regulatory elements in winter wheat.

BACKGROUND: Winter wheat undergoes vernalization, a process activated by prolonged exposure to low temperatures. During this phase, flowering signals are generated and transported to the apical meristems, stimulating the transition to the inflorescence meristem while inhibiting tiller bud elongation. Although some vernalization genes have been identified, the key cis-regulatory elements and precise mechanisms governing this process in wheat remain largely unknown. RESULTS: In this study, we construct extensive epigenomic and transcriptomic profiling across multiple tissues-leaf, axillary bud, and shoot apex-during the vernalization of winter wheat. Epigenetic modifications play a crucial role in eliciting tissue-specific responses and sub-genome-divergent expressions during vernalization. Notably, we observe that H3K27me3 primarily regulates vernalization-induced genes and has limited influence on vernalization-repressed genes. The integration of these datasets enables the identification of 10,600 putative vernalization-related regulatory elements including distal accessible chromatin regions (ACRs) situated 30Kb upstream of VRN3, contributing to the construction of a comprehensive regulatory network. Furthermore, we discover that TaSPL7/15, integral components of the aging-related flowering pathway, interact with the VRN1 promoter and VRN3 distal regulatory elements. These interactions finely regulate their expressions, consequently impacting the vernalization process and flowering. CONCLUSIONS: Our study offers critical insights into wheat vernalization's epigenomic dynamics and identifies the putative regulatory elements crucial for developing wheat germplasm with varied vernalization characteristics. It also establishes a vernalization-related transcriptional network, and uncovers that TaSPL7/15 from the aging pathway participates in vernalization by directly binding to the VRN1 promoter and VRN3 distal regulatory elements.

Transitions in chromatin conformation shaped by fatty acids and the circadian clock underlie hepatic transcriptional reorganization in obese mice.

The circadian clock system coordinates metabolic, physiological, and behavioral functions across a 24-h cycle, crucial for adapting to environmental changes. Disruptions in circadian rhythms contribute to major metabolic pathologies like obesity and Type 2 diabetes. Understanding the regulatory mechanisms governing circadian control is vital for identifying therapeutic targets. It is well characterized that chromatin remodeling and 3D structure at genome regulatory elements contributes to circadian transcriptional cycles; yet the impact of rhythmic chromatin topology in metabolic disease is largely unexplored. In this study, we explore how the spatial configuration of the genome adapts to diet, rewiring circadian transcription and contributing to dysfunctional metabolism. We describe daily fluctuations in chromatin contacts between distal regulatory elements of metabolic control genes in livers from lean and obese mice and identify specific lipid-responsive regions recruiting the clock molecular machinery. Interestingly, under high-fat feeding, a distinct interactome for the clock-controlled gene Dbp strategically promotes the expression of distal metabolic genes including Fgf21. Alongside, new chromatin loops between regulatory elements from genes involved in lipid metabolism control contribute to their transcriptional activation. These enhancers are responsive to lipids through CEBPß, counteracting the circadian repressor REVERBa. Our findings highlight the intricate coupling of circadian gene expression to a dynamic nuclear environment under high-fat feeding, supporting a temporally regulated program of gene expression and transcriptional adaptation to diet.

Rational Design and Screening of Synthetic Promoters in Chlamydomonas reinhardtii.

Synthetic promoters are powerful tools to boost the biotechnological potential of microalgae as eco-sustainable industrial hosts. The increasing availability of transcriptome data on microalgae in a variety of environmental conditions allows to identify cis-regulatory elements (CREs) that are responsible for the transcriptional output. Furthermore, advanced cloning technologies, such as golden gate-based MoClo toolkits, enable the creation of modular constructs for testing multiple promoters and a range of reporter systems in a convenient manner. In this chapter, we will describe how to introduce in silico-identified CREs into promoter sequences, and how to clone the modified promoters into MoClo compatible vectors. We will then describe how these promoters can be evaluated and screened for transgene expression in an established microalgal model for genetic perturbation, i.e., Chlamydomonas reinhardtii.

The contribution of silencer variants to human diseases.

BACKGROUND: Although disease-causal genetic variants have been found within silencer sequences, we still lack a comprehensive analysis of the association of silencers with diseases. Here, we profiled GWAS variants in 2.8 million candidate silencers across 97 human samples derived from a diverse panel of tissues and developmental time points, using deep learning models. RESULTS: We show that candidate silencers exhibit strong enrichment in disease-associated variants, and several diseases display a much stronger association with silencer variants than enhancer variants. Close to 52% of candidate silencers cluster, forming silencer-rich loci, and, in the loci of Parkinson's-disease-hallmark genes TRIM31 and MAL, the associated SNPs densely populate clustered candidate silencers rather than enhancers displaying an overall twofold enrichment in silencers versus enhancers. The disruption of apoptosis in neuronal cells is associated with both schizophrenia and bipolar disorder and can largely be attributed to variants within candidate silencers. Our model permits a mechanistic explanation of causative SNP effects by identifying altered binding of tissue-specific repressors and activators, validated with a 70% of directional concordance using SNP-SELEX. Narrowing the focus of the analysis to individual silencer variants, experimental data confirms the role of the rs62055708 SNP in Parkinson's disease, rs2535629 in schizophrenia, and rs6207121 in type 1 diabetes. CONCLUSIONS: In summary, our results indicate that advances in deep learning models for the discovery of disease-causal variants within candidate silencers effectively "double" the number of functionally characterized GWAS variants. This provides a basis for explaining mechanisms of action and designing novel diagnostics and therapeutics.