Biochimica et Biophysica Acta (BBA) - Gene Regulatory Mechanisms

Chromatin remodeler imitation switch (ISWI) plays an important role in regulating chromatin structure through sliding and spacing nucleosomes. Despite the enormous progress in regulatory elements and mechanisms of the activity of ISWI in recent years, there are still some unclear structures and mechanisms in different species. Here, we studied the ATPase activity and nucleosome binding affinity of Chaetomium thermophilum ISWI (hereafter referred to as CtISWIWT) and several mutants, further proving the importance of these mutated residues in the inhibition of AutoN. We also analyzed the effects of dsDNA and ssDNA on ATPase activity of CtISWI, suggesting the potential interaction between HSS and ATPase domain. Notably, we provided a predicted structural model based on the sequence of CtISWIWT, proposing a two-step activating mechanism of conformation change and activity regulation. Taken together, our findings elucidate a different model of ISWI self-maintenance and action, providing a new mechanism of regulation supporting chromatin remodeling. HighlightsO_LIStructural modeling of ISWI: a chromatin remodeler that couples to ATP hydrolysis to slide and space composition of nucleosome. C_LIO_LIThe ATPase activity of ISWI can be stimulated by exogenous DNA, with opposite promoting effects by dsDNA and ssDNA. C_LIO_LIThe mechanism by which ISWI is activated upon binding with nucleosome has been the subject of debate, and a more comprehensive mechanism for regulation of ISWI activity. C_LI SignificanceIn the past decades, a variety of regulatory mechanisms of the activity of chromatin remodeling factor ISWI have been proposed. Based on the hypothesis of nucleosome complex structure analysis, these studies attempted to explore the mechanism of chromatin remodeling, a gene expression regulation activity. However, previous studies have basically focused on the binding and regulation mechanism of ISWI ATPase domain and nucleosomes, without mentioning the activity mode of full-length ISWI. Therefore, our study mainly focuses on the nearly full-length ISWI containing HSS domain, exploring the mechanism of the active state transition of ISWI in remodeling activities from this perspective. It enriches and supplements the research on chromatin remodeling, an important physiological activity.

The chromatin remodeling protein, dIno80 (Drosophila Ino80) regulates homeotic genes. We show that Ino80, along with Trx and ETP (Enhancer of Trithorax and Polycomb) proteins, interacts with two Polycomb/Trithorax Responsive Elements (PRE/TRE), iab-7 and bxd PRE in flies and the larval imaginal discs. In S2 cells, dIno80 localizes to the endogenous iab-7 and bxd-PREs. The localization of Ino80 and Pleiohomeotic (Pho) at the PRE is sensitive to the cellular abundance of each other; when levels of Ino80 are limiting, there is increased Pho enrichment, and Pho knock-down leads to increased enrichment of Ino80. We demonstrate that over-expression of dIno80 rescues the pupal lethality in pleiohomeotic (pho) deficient flies, which suggests that dIno80 has a role in cellular memory. The apparent competition between Pho and Ino80 for binding at the PRE indicates that Ino80 may act as a potential recruiter of the regulatory complex in addition to being a chromatin remodeler.\n\nAuthor SummaryThe null mutants of Pho and dIno80 show lethality at different stages of development in the fly, implying that they may function independent of each other. The observation that Pho-lethality can be rescued by overexpression of dIno80 with significant penetrance and that Ino80 has its own DNA binding domain, led us to predict that Ino80 may have Pho-independent functions, perhaps through non-canonical complexes. In the current study, we show that dIno80 interacts with bxd and iab-7 PRE in cooperation with Polycomb and Trithorax proteins and regulate the homeotic genes. The effect of knock-down or mutation of dIno80 results in altered phenotype in adult flies and rescue of Lac-Z expression in imaginal discs, in parallel with similar effect of Pho mutation or knock-down. We provide evidence of direct interaction of dIno80 with iab7- and bxd-PRE using chromatin immunoprecipitation. The dIno80 localization in and around the PRE sequence was enhanced in the absence of Pho, indicating competition between Pho and dIno80 for binding at the PRE.

Sonic hedgehog (SHH) signaling from the frontonasal ectodermal zone (FEZ) is a key regulator of craniofacial morphogenesis. Along with SHH, pre-B-cell leukemia homeobox (PBX) transcription factors regulate midfacial development. PBXs act in the epithelium during fusion of facial primordia, but their specific interactions with SHH have not been fully investigated. We hypothesized that PBX1/3 regulate SHH expression in the FEZ by activating or repressing transcription. The hypothesis was tested by manipulating PBX1/3 expression in chick embryos and profiling epigenomic landscapes at early developmental stages. PBX1/3 expression was perturbed in the chick face beginning at stage 10 (HH10) using RCAS viruses, and the resulting SHH expression was assessed at HH22. Overexpressing PBX1 expanded SHH expression, while overexpressing PBX3 decreased SHH expression. Conversely, reducing PBX1 expression decreased SHH expression, but reducing PBX3 induced ectopic SHH expression. We performed ATAC-seq and mapped binding of PBX1 and PBX3 with ChIP-seq on the FEZ at HH22 to assess direct interactions of PBX1/3 with the SHH locus. These multi-omics approaches uncovered a 400 bp PBX1-enriched element within intron 1 of SHH (chr2:8,173,222-8,173,621). Enhancer activity of this element was demonstrated by electroporation of reporter constructs in ovo and luciferase reporter assays in vitro. When bound by PBX1, this element upregulates transcription, while it downregulates transcription when bound by PBX3. The present study identifies a cis-regulatory element, named SFE1, that interacts with PBX1/3 to modulate SHH expression in the FEZ and establishes that PBX1 and PBX3 play complementary roles in SHH regulation during embryonic development.

Histone modifications are one of key mechanisms to regulate gene expression. Ash1 is a histone H3K36 methyltransferase and involved in gene activation. Ash1 forms a large complex with Mrg15 and Caf1/p55/Nurf55/RbAp48 (AMC complex). Ash1 subunit alone has very low activity due to the auto-inhibition and the binding of Mrg15 releases the auto-inhibition. Caf1 is a scaffolding protein commonly found in several chromatin modifying complexes. Caf1 has an ability to sense unmodified histone H3K4 residue. However, the role of Caf1 in AMC complex has not been investigated. Here, we dissected the interaction among the AMC complex subunits, revealing that Caf1 uses the histone H4 binding pocket to interact with Ash1 near the histone binding module cluster. Furthermore, we show that H3K4 methylation inhibits AMC HMTase activity via Caf1 sensing unmodified histone H3K4 to regulate the activity in an inter-nucleosomal manner, suggesting that there is a crosstalk between H3K4 and H3K36 methylations. Our work reveals a delicate regulatory mechanism of AMC histone H3K36 methyltransferase complex.

NifA is a {sigma}54 activator that turns on bacterial nitrogen fixation under reducing conditions and when fixed cellular nitrogen levels are low. The redox sensing mechanism in -proteobacterial NifA is poorly understood. In this work, we examine if a Cys pair that is part of a C(X)5C motif and located immediately upstream of NifAs DNA binding domain is involved in redox sensing in NifA from the -proteobacterium Gluconacetobacter diazotrophicus (Gd). We hypothesize that the Cys residues redox state may directly influence the DNA binding domains DNA binding affinity and/or alter the proteins oligomeric sate. Two DNA binding domain constructs were generated, a longer construct (2C-DBD), consisting of the DNA binding domain with the upstream Cys pair, and a shorter construct (NC-DBD) that lacks the Cys pair. The Kd of NC-DBD for its cognate DNA sequence (nifH-UAS) is equal to 20.0 M. The Kd of 2C-DBD for nifH-UAS when the Cys pair is oxidized is 34.5 M. Reduction of the disulfide bond does not change the DNA binding affinity. Additional experiments indicate that the redox state of the Cys residues does not influence the secondary structure or oligomerization state of the NifA DNA binding domain. Together, these results demonstrate that the Cys pair upstream of the DNA binding domain of Gd-NifA does not regulate DNA binding or domain dimerization in a redox dependent manner. This suggests that other Cys residues in NifA, such as those located in the central AAA+ domain, are responsible for redox sensing.

Most known modification-dependent restriction endonucleases target 5-methylcytosine, only a few N6-methyladenine (6mA)-dependent restriction endonucleases have been well-characterized, and the majority of them recognize the G6mATC motif (e.g., DpnI, HHPV4I). Here, we report the identification of a novel 6mA-dependent DNA-binding protein from Vibrio cholerae, VchI, which specifically recognizes the G6mAG motif. VchI contains a winged helix (wH) domain that is homologous to the wH domain in DpnI. However, several key residues involved in 6mA recognition differ between VchI and DpnI, which may contribute to the discrepancy in their recognition specificities. These findings advance our understanding of prokaryotic 6mA modification diversity and the 6mA recognition mechanism of the wH domain, while simultaneously providing an innovative tool for epigenetic research.

The m6A epitranscriptomic mark is the most abundant and widespread internal RNA chemical modification, which through the control of RNA acts as an important actor of eukaryote reproduction, growth, morphogenesis and stress response. The main m6A readers constitute a super family of proteins with hundreds of members that share a so-called YTH RNA binding domain. The majority of YTH proteins carry no obvious additional domain except for an Intrinsically Disordered Region (IDR). In Arabidopsis thaliana IDRs are important for the functional specialization among the different YTH proteins, known as Evolutionarily Conserved C-Terminal region, ECT 1 to 12. Here by studying the ECT2 protein and using an in vitro biochemical characterization, we show that full length ECT2 and its YTH domain alone have a distinct ability to bind m6A, conversely to previously characterized YTH readers. We identify peptide regions outside of ECT2 YTH domain, in the N-terminal IDR, that regulate its binding to m6A-methylated RNA. Furthermore, we show that the selectivity of ECT2 binding for m6A is enhanced by a high uridine content within its neighboring sequence, where ECT2 N-terminal IDR is believed to contact the target RNA in vivo. Finally, we also identify small structural elements, located next to ECT2 YTH domain and conserved in a large set of YTH proteins, that enhance its binding to m6A-methylated RNA. We propose from these findings that some of these regulatory regions are not limited to ECT2 or YTH readers of the flowering plants but may be widespread among the eukaryotic YTH readers.

Adipose tissue is a key metabolic organ for carbohydrate and lipid metabolism. Multiple nutrient sensing pathways and metabolic enzymes operate in adipocytes to control energy production and lipid mobilization in response to physiological conditions. Important regulators of glycolytic and lipid metabolism in many species are members of the estrogen-related receptor (ERR) family of nuclear receptors. We previously showed that ERR is required for transcriptional regulation of glycolytic and pentose phosphate pathway enzymes in adult Drosophila females, consistent with what is observed in ERR mutant males and larvae. However, the cell-type specific targets of ERR in adipose and other tissues have not been fully elucidated. Here, using a modified targeted DamID approach (NanoDam), we identified ERR occupancy specifically in adult adipocytes at the loci that encode genes involved in glycolysis, the pentose phosphate pathway, and fatty acid metabolism. Overall, our results predict that ERR centrally functions in adipose tissue to regulate multiple metabolic pathways.

CHT7 is a regulator of quiescence repression and TAG degradation between the nitrogen deprived and the nitrogen replenished states in Chlamydomonas reinhardtii. Initially it was thought that the CHT7s repression activity is managed by its DNA binding CXC domain which is a tandem repeat of two cysteine rich subdomains. Later, it was found that the CXC (CHT7_CXC) domain is effectively dispensable for CHT7s activities. Rather, CHT7s predicted protein binding domains are proposed to be involved in gene regulation activities by binding through other repressors in the cell. Yet, it remains unclear why and how CHT7 manages to refrain its own CXC domain from participating in any transcriptional activities. The question becomes more intriguing, because CXC binding regions are available in promoter regions of some of the misregulated genes in the CHT7 mutant (cht7). Through the combination of biophysical experiments and molecular dynamics approaches, we have studied the DNA recognition behavior of CHT7_CXC. The results show that CHT7_CXC domain is highly selective towards DNA sequences and this selectivity is imparted due to the differential binding abilities of the CXC subdomains. Further, to understand if the case is - that CXC looses its DNA binding capabilities in the vicinity of other repressor molecules, we carried out CHT7_CXCs DNA binding stability test by simulating the spatial constraint conditions using the AsLOV2- CXC fusion. Our test results show limited ability of CHT7_CXC to withstand steric forces and provide insights to why and how algal cells may hold back CHT7_CXCs indulgence in quiescence repression. SignificanceMicroalgae, under nutrient rich conditions, provide biomass. Whereas, nutrient deprivation leads to accumulation of biofuel feedstock, but cells enter quiescence. Net enhancement in feedstock, therefore relies on the precision of the quiescence regulator. In Chlamydomonas reinhardtii, CHT7 is a central regulator of quiescence. Surprisingly, rather than using its own DNA binding domain (DBD) for the regulatory activities, CHT7 recruits external transcriptional regulators using its non DBDs. To ensure smooth functioning, CHT7s DBD must rapidly switch to inactive form. Modifications in DNA binding profiles of DBDs due to non DBDs are seen in transcription factors of many organisms. The switching mechanism discussed could therefore be a generic approach of timely regulation of individual components of the complex transcriptional machineries.

Chromosomes in Drosophila melanogaster are organized into distinct topologically associated domains delimited by boundaries bound by insulator proteins. The insulator protein BEAF (Boundary Element-Associated Factor of 32kDa) plays roles in both chromatin organization and transcriptional regulation, yet its precise molecular mechanisms remain elusive. To examine the role of BEAF in insulator function we used yeast 2-hybrid assays, pull-down assays with bacterially expressed proteins, and bimolecular fluorescence complementation assays in S2 cells to characterize interactions between BEAF and three co-insulator proteins: CP190, Pzg, and Chro. Our studies pinpoint minimal regions of CP190, Pzg, and Chro that directly interact with BEAF, as well as parts of BEAF that are crucial for its interactions with these co-insulator proteins. Functional analyses in transfected S2 cells revealed distinct regulatory roles for BEAF in association with CP190, Pzg, and Chro. CP190 showed a weak interaction with BEAF, and CP190 bound 2.3 kb upstream of promoter-proximal BEAF could not loop out the intervening DNA to effectively communicate with BEAF for luciferase reporter gene activation. In contrast, more robust interactions were observed between BEAF and both Pzg and Chro. Both could also effectively interact with BEAF from a distance to activate the reporter gene. It is likely that the role of CP190 in long-range insulator interactions is mediated by insulator proteins other than BEAF, although BEAF could help stabilize interactions. On the other hand, we propose that BEAF can directly collaborate with Pzg and Chro to mediate long-range chromatin interactions.

The spatial organization of bacterial chromosomes is crucial for cellular functions. It remains unclear how bacterial chromosomes adapt to high temperature stress. This study delves into the 3D genome architecture and transcriptomic responses of Escherichia coli under heat stress condition to unravel the intricate interplay between chromosome structure and environmental cues. By examining the role of macrodomains, chromosome interaction domains (CIDs), and nucleoid-associated proteins (NAPs), this work unveils the dynamic changes in chromosome conformation and gene expression patterns induced by high temperature stress. It was observed that under heat stress, short-range interaction frequency of chromosome decreased, while the long-range interaction frequency of the Ter macrodomain increased. Furthermore, two metrics, namely, Global Compactness (GC) and Local Compactness (LC), were devised to measure and compare the compactness of chromosomes based on their 3D structure models. The findings in this work shed light on the molecular mechanisms underlying thermal adaptation and chromosomal organization in bacterial cells, offering valuable insights into the complex interrelationships between environmental stimuli and genomic responses.

The mechanisms of gene expression regulation are essential for cell identity and function, and their disruption usually leads to human disease. The 3 untranslated region (UTR) of mRNA contains important regulatory elements of gene expression, including upstream sequence elements (USEs) that are cis-regulatory sequences localized upstream of polyA signals (PAS). One of the best functionally characterised USEs is located in the 3UTR of the Drosophilas polo gene, which disruption leads to critical phenotypes in adult flies. In this work we found that the USE of the Drosophilas polo gene (DplUSE) is also found in in the 3UTR of vertebrate genes, including zebrafish, mouse and human genes, showing higher levels of conservation than the whole 3UTR sequence. Using reporter assays, we show that DplUSE is able to increase gene expression in vitro in human cell lines and in vivo in zebrafish embryos. Importantly, in humans, the DplUSE containing genes are enriched for genes associated to serious diseases such as Congenital abnormalities and Malignant neoplasms, illustrating the potential of this sequence to modulate genes with relevant biological functions and related with human health. Concomitantly, when sequestering the molecular machinery that operates at the DplUSE using a dominant negative strategy, we show that this is enough to dysregulate DplUSE containing genes in human cells and disrupt proper embryo development in zebrafish. Aiming to understand the molecular mechanism operating at the DplUSE, we identified three RNA binding proteins (RBP) that specifically bind to the DplUSE in vertebrates. Importantly, one of such RBPs is PTBP1, the vertebrate orthologue of the fruit flys RBP Heph, that was demonstrated to be required for the DplUSE function in Drosophila. To test if PTBP1is essential for the DplUSE function, as observed in Drosophila, we depleted PTBP1 from human cells and observed a downregulation of the expression of DplUSE containing genes, demonstrating that the molecular mechanisms that operate at DplUSE are ultra-conserved. Finally, we explored if variants in DplUSE consensus could be associated to human disease. We found a reported single nucleotide polymorphism (SNP; rs3087967) that is associated with malignant tumor of colon and generates an ectopic consensus of DplUSE in the 3 UTR of the tumorigenic POU2AF2/C11orf53 gene. We further show that this ectopic DplUSE motif causes a gain-of-function in vivo in zebrafish gut cells, suggesting its involvement in colon cancer development. These results show that a short motif present in the 3UTR of genes from phylogenetically distant bilaterians, from fruit flies to humans, control genes expression through an ultra-conserved mechanism involving RBPs binding and its dysregulation might impact in human disease.

Transcription factors (TFs) regulate gene expression by specifically binding to DNA targets. Many factors have been revealed to influence TF-DNA binding specificity. Coevolution of residues in proteins occurs due to a common evolutionary history. However, it is unclear how coevolving residues in TFs contribute to DNA binding specificity. Here, we systematically analyzed TF-DNA interactions from high-throughput experiments for seven TF families, including Homeobox, HLH, bZIP_1, Ets, HMG_box, zf-C4 and Zn_clus TFs. Based on TF-DNA interactions, we detected TF subclass determining sites (TSDSs) defining the heterogeneity of DNA binding preference for each TF family. We showed that the TSDSs were more likely to be coevolving with TSDSs than with non-TSDSs, particularly for Homeobox, HLH, Ets, bZIP_1 and HMG_box TF families. Mutation of the highly coevolving residues could significantly reduce the stability of TF-DNA complex. The distant residues from the DNA interface also contributed to TF-DNA binding activity. Overall, our study gave evidence of the functional importance of coevolved residues in refining transcriptional regulation and provided clues to the application of engineered DNA-binding domains and protein.

Protein-DNA interactions are crucial for cellular processes, but current quantification methods lack sensitivity. We developed the Real-time PCR-based DNA Binding Assay (RP-DBA) to detect and quantify these interactions. The target protein, expressed as a Strep-tag II fusion, is purified and incubated with double-stranded DNA probes containing 4 bp 3 overhangs. Protein-DNA complexes are immobilized on Strep-Tactin beads, washed, and eluted. A complementary single-stranded DNA amplification arm is added, extended by Taq polymerase, and quantified via qPCR with SYBR Green. RP-DBA enables real-time kinetic analysis and is 4- to 10-fold more sensitive than EMSA, depending on amplification arm length (30-90 bp). Its simplicity, speed, accuracy, and high-throughput potential make it a valuable tool for advancing molecular biology research.

Transcription factors (TFs) regulate gene expression by binding cis-regulatory DNA elements, yet how trans-regulatory characteristics such as TF affinity, concentration, and localization interact with cis-regulatory elements remains largely unclear. We systematically analyzed TF affinity mutants across abundance, and localization states and found that promoter binding-site strength most readily modulated expression levels, followed by TF localization and concentration, while affinity variations were mainly buffered. We further uncover performance trade-offs between TF abundance, localization, and affinity. Together, these results reveal how trans and cis factors collectively shape gene-regulatory output.

Limited oxygen supply or hypoxia can impair fetal development and lead to developmental disorders, but the molecular mechanism underlying this phenomenon remains poorly understood. It is also well known that hypoxia results in transcriptomic alterations and epigenetic reprogramming. Drosophila melanogaster (fruit fly) has been used for decades as a powerful model to dissect the molecular mechanisms regulating development. To better understand the role of early hypoxic stress on development, we performed single-cell joint analysis of chromatin accessibility and transcriptome to characterize the influence of hypoxia on Drosophila embryonic development. We identified hypoxia-induced alterations in both gene expression and chromatin accessibility across 22 cell groups, especially in the genes regulating organogenesis and development of neuronal, tracheal, and muscular systems, including a reduction of germ cells, suggesting a long-lasting influence of hypoxic stress at an early embryonic stage on development and reproduction. In summary, this study demonstrates that early embryonic hypoxia induces cell type- and dose-dependent changes in chromatin accessibility and gene expression, leading to distinct developmental phenotypic responses, such as reduced number of germ cells under both 3% and 5% O2. We further conclude that the tramtrack (ttk) gene is critical in germ cell development and reproduction in Drosophila melanogaster.

Brasilicardin A, BraA, is a secondary metabolite produced by the bacterium Nocardia terpenica, and a promising drug due to its potent immunosuppressive activity and low cytotoxicity. Currently, a semisynthetic approach confers production of a complete compound but suffers from insufficient heterologous biosynthesis of BraA intermediates used in the chemical semi-synthesis steps leading to only lab scale quantities of the compound. A better understanding of the involved gene expression regulatory pathways within the brasilicardin biosynthetic gene cluster, Bra-BGC, is a prerequisite to further improve production titers. However, the transcriptional regulation of the Bra-BGC has only been superficially analyzed, till now. In this study, we comprehensively analyze the functions of several unstudied transcriptional regulators, KstR, SdpR and OmpR, encoded within the close vicinity of the Bra-BGC, and delve into the role of the previously described cluster-situated activator Bra12. We present, that Bra12 and the novel regulator SdpR, bind several DNA sequences located in the promoter regions of the genes essential for BraA biosynthesis. Subsequently, we demonstrate the complex regulatory network through which both regulators are capable of controlling activity of those gene promoters and thus gene expression in Bra-BGC. Furthermore, using the heterologous producer strain Amycolatopsis japonicum, we present, that Bra12 and SdpR regulators play opposite roles in brasilicardin congener biosynthesis. Finally, we propose a comprehensive model of multilevel gene expression regulation in Bra-BGC and propose the roles of locally encoded transcriptional regulators.

The accumulation of manganese ions is crucial for scavenging reactive oxygen species (ROS) and protecting the proteome of Deinococcus radiodurans (Dr). However, metal homeostasis still needs to be tightly regulated to avoid toxicity. DR2539, a dimeric transcription regulator, plays a key role in Dr manganese homeostasis. Despite comprising three well-conserved domains: a DNA binding domain, a dimerization domain, and an ancillary domain, both the metal ion activation mechanism and the DNA recognition mechanism remain elusive. In this study, we present biophysical analyses and the structure of the dimerization and DNA binding domains of DR2539 in its holo form and in complex with the 21 bp pseudo-palindromic repeat of the dr1709 promotor region. These findings shed light into the activation and recognition mechanisms. The dimer presents eight manganese binding sites that induce structural conformations essential for DNA binding. The analysis of the protein-DNA interfaces elucidates the significance of Tyr59 and helix H3 sequence in the interaction with the DNA. Finally, the structure in solution as determined by small angle X-ray scattering experiments and supported by AlphaFold modelling provides a model illustrating the conformational changes induced upon metal binding.

Members of the R2R3-MYB4 subgroup are well-known negative regulatory transcription factors of phenylpropane and lignin pathways. In this study, we found that transgenic tobacco plants overexpressing a R2R3-MYB4 subgroup gene from Camellia sinensis (CsMYB4a) showed inhibited growth that was not regulated by phenylpropane and lignin pathways, and these plants exhibited altered sensitivity to synthetic auxin 1-naphthaleneacetic acid (-NAA) treatment. An auxin/indole-3-acetic acid 4 (AUX/IAA4) gene from Camellia sinensis (CsIAA4) participating in the regulation of the auxin signal transduction pathway was screened from the yeast two-hybrid library with CsMYB4a as the bait protein, and tobacco plants overexpressing this gene showed a series of auxin-deficiency phenotypes, such as dwarfism, small leaves, reduced lateral roots, and a shorter primary root. CsIAA4 transgenic tobacco plants were less sensitive to exogenous -NAA than control plants, which was consistent with the findings for CsMYB4a transgenic tobacco plants. The knockout of the endogenous NtIAA4 gene (a CsIAA4 homologous gene) in tobacco plants alleviated growth inhibition in CsMYB4a transgenic tobacco plants. Furthermore, protein-protein interaction experiments proved that domain II of CsIAA4 is the key motif for the interaction between CsIAA4 and CsMYB4a and that the degradation of CsIAA4 is prevented when CsMYB4a interacts with CsIAA4. In summary, our results suggest that CsMYB4a is a multifunctional transcription factor that regulates the auxin signaling pathway, phenylpropane and lignin pathways. This study provides new insights into the multiple functions of R2R3-MYB4 subgroup members as a group of well-known negative regulatory transcription factors. One-sentence summaryCsMYB4a act as multifunctional transcription factor that regulates the auxin signaling pathway, phenylpropane and lignin pathways.

Nickel is an occupational and environmental toxicant associated with a number of diseases in humans including pulmonary fibrosis, bronchitis and lung and nasal cancers. Our earlier studies showed that the nickel-exposure-induced genome-wide transcriptional changes, which persist even after the termination of exposure may underlie nickel pathogenesis. However, the mechanisms that drive nickel-induced persistent changes to the transcriptome remain elusive. To elucidate the mechanisms that underlie nickel induced long-term transcriptional changes, in this study, we examined the transcriptome and the epigenome of human lung epithelial cells during nickel exposure and after the termination of exposure. We identified two categories of persistently differentially expressed genes based on the timing of expression changes: i) the genes that were differentially expressed during nickel exposure; and ii) the genes that were differentially expressed only after the termination of nickel exposure. Interestingly, the majority of nickel-induced transcriptional changes occurred only after the termination of exposure. We found robust genome-wide alterations to the activating histone modification, H3K4me3, after the termination of nickel exposure, which coincided with the post-exposure gene expression changes. In addition, we found significant post-exposure alterations to the repressive histone modification, H3K27me3. By uncovering a new category of transcriptional and epigenetic changes, which occur only after the termination of exposure, this study sheds new light on the post-exposure effects of nickel and provides a novel understanding of the long-term deleterious consequences of nickel exposure on human health.