Drosophila DAxud1: A New Element in Transcriptional Pausing Complex Stabilization

Drosophila DAxud1: A New Element in Transcriptional Pausing Complex Stabilization | Research Square window.SnipcartSettings = { analytics: { enabled: false } }; (function() { var accessVector = localStorage.getItem('access_vector') || ''; window.dataLayer = window.dataLayer || []; if (accessVector) { window.dataLayer.push({ user: { profile: { profileInfo: { snid: accessVector } } } }); } })(); (function(w,d,s,l,i){w[l]=w[l]||[];w[l].push({'gtm.start':new Date().getTime(),event:'gtm.js'});var f=d.getElementsByTagName(s)[0],j=d.createElement(s),dl=l!='dataLayer'?'&l='+l:'';j.async=true;j.src='https://www.googletagmanager.com/gtm.js?id='+i+dl;f.parentNode.insertBefore(j,f);})(window,document,'script','dataLayer','GTM-K279D39R'); Browse Preprints In Review Journals COVID-19 Preprints AJE Video Bytes Research Tools Research Promotion AJE Professional Editing AJE Rubriq About Preprint Platform In Review Editorial Policies Our Team Advisory Board Help Center Sign In Submit a Preprint Cite Share Download PDF Research Drosophila DAxud1: A New Element in Transcriptional Pausing Complex Stabilization JM Zuñiga-Hernández, C Meneses, M Bastías, ML Allende, Alvaro Glavic This is a preprint; it has not been peer reviewed by a journal. https://doi.org/ 10.21203/rs.3.rs-1040684/v1 This work is licensed under a CC BY 4.0 License Status: Posted Version 1 posted You are reading this latest preprint version Abstract Background A rapid transcriptional response under an acute stimulus is common in all cellular systems and is an adaptation that allows tolerance to environmental changes. A gene group that has been studied because of its fast response and cytoprotective effects are the hsp genes (encodingHeat Shock Proteins(HSPs), conserved chaperones).. Under normal conditions, the mRNA and protein levels of the main hsp genes are low but they increase rapidly upon heat shock (HS). This is achieved due to the presence of an RNA Polymerase II pausing complex located +30-50 bp from TSS. This complex maintains a partially synthesized RNA strand of said length, poised to resume synthesis, and undergoes subsequent transcriptional inactivation to restore transcript levels after environmental stabilization. Methods The Gal4/UAS system was used to modify dAxud1 expression in a tissue specific manner. A DAxud1-GFP fusion was expressed in salivary glands to perform polytene chromosome immunofluorescence and chromatin immunoprecipitation. DAxud1 genome occupancy data was achieved expressing Dam-DAxud1 in imaginal wing discs using Gal4/UAS (TaDa-seq). Results Using TaDa-seq, we demonstrate that DAxud1 protein is present mainly near the TSS of significant occupied genes, most frequently in the first intron. This results also revealed DAxud1 is present in hsp genes, mainly in promoter zone. Following these results, we found that, under dAxud knockdown, larvae and adults flies have a diminished thermotolerance, despite showing an increase in hsp transcripts in larval tissues. We performed polytene chromosome immunofluorescence for DAxud1-GFP, revealing extensive, but dynamic localization on chromatin in hsp70 loci . This was confirmed with chromatin immunoprecipitation. We also found that DAxud1 overexpression leads to an enrichment of RNA Polymerase II at the 5’ end of the hsp70 gene, with a decrease in its transcripts. Importantly, we show interaction of DAxud1 with NELF-B, a component of the transcriptional pausing complex, and knockdown of both genes individually has similar effects on hsp70 transcription. Conclusion DAxud1 protein is a component of chromatin, that relocates under stress conditions such as heat shock, playing a role in maintaining RNA Polymerase II stalled at the 5’ of hsp70, possibly through a pausing mechanism based on its interaction with NELF-B. Epigenetics & Genomics heat shock hsp70 transcription pausing complex DAxud1 NELF Figures Figure 1 Figure 2 Figure 3 Figure 4 Introduction In natural environments, organisms are exposed to different stresses that affect their survival. During evolution, different mechanisms have emerged to respond to stress, some of them widely conserved among various kingdoms, but with differences in their temporal regulation and cell types involved. One of the conserved mechanisms for responding to stress is rapid biosynthesis of Heat Shock Proteins (HSPs), in which their mRNA levels can rise several fold in a few minutes, mainly due to transcriptional promoting elements preloaded at the core promoter and within the 5’ end of the gene body [1, 2]. HSPs were initially associated with heat stress responses. Further research revealed that these proteins also behave as important cytoprotectors under other stress conditions, including oxidative stress, accumulation of protein aggregates, an acidic environment, and osmolarity changes [3]. In animals, hsp transcriptional induction is very fast; the hsp70 promoter recruits specific factors in response to stress conditions [4], resulting in significative rise of hsp70 mRNA in less than 5 minutes [5]. This recruitment during heat shock is commanded mainly by HSF (Heat Shock Factor), a transcriptional activator, one of five known chromatin proteins involved in the heat stress transcriptional response, that induces Hsp70 and other small HSPs (Hsp27, Hsp26) after acute thermal stress. However, the repertoire of chromatin associated proteins varies for other stress conditions and the differences are still poorly understood [6–8]. We will focus mainly on the hsp70 paralogs because there is significant knowledge implicating the chromatin state of their promoters and the pausing phenomenon [9–14], as well as additional factors regulating their efficient transcriptional induction. In Drosophila melanogaster , the control of hsp70 transcription depends on at least 12 proteins that relocate to the hsp70 loci after two minutes of exposure to 37°C in 3 rd instar larvae [4, 13–16]. These proteins participate regulating chromatin condensation, recruiting transcription factors like HSF, engaging in transcriptional pausing, and controlling the initiation/elongation of hsp70 transcription. In hsp70 genes, transcriptional pausing is crucial for a quick transcriptional response upon stress or extracellular signaling. When pausing, RNA Polymerase II elongates only 30-50 bp from the transcription start site (TSS) and remains in that position until induced, where it resumes elongation and provides a transcriptional burst of the stress response genes [14, 17]. Transcriptional pausing is mediated by the NELF (NELF A, B, C/D, E) and DSIF (Spt5, Spt6) complexes, which retain RNA Polymerase II. After induction by heat shock, HSF recruits the p-TEFb kinase that phosphorylates NELF, DSIF and RNA Polymerase II, causing complete dissociation of NELF and the and the escape of RNA Polymerase, with DSIF acting as an elongation factor [4, 18]. In this work, we present evidence that introduces DAxud1 as an additional regulatory factor that participates in hsp transcription regulation, through its positioning at the transcription start site (TSS) and by stabilizing the pausing complex via its interaction with NELF and DSIF complexes. The dAxud1 gene is the only Drosophila ortholog of the human and mice CSRNP protein family (Cysteine Serine Rich Nuclear Proteins). Members of this protein family, which are conserved only in metazoans, have been proposed to be putative transcription factors with specific roles in apoptosis and neural development [19–21]. Initial studies proposed a physiological function for one of these proteins as a tumor suppressor in human tissues [19]. Further characterization of its ortholog in Drosophila suggested that its proapoptotic function relies on the activation of the JNK pathway [21]. Other studies indicate a strong association of the CSRNP family with stress responses as their mRNAs rise in cells exposed to different kinds of stressful stimuli such as cold [22], oxidative stress [23, 24], bacterial infection [25], pressure [26], and acute stress by sprint running in skeletal muscle [27]. These observations prompted us to analyze the molecular bases of the adaptive and non-adaptive roles of DAxud1 during the stress response. Here, in a screen using whole-genome occupancy analysis (TaDa-seq), we reveal how DAxud1 has a connection to the general stress response by analyzing its occupancy of the hsp70 locus . In this gene, DAxud1 associates with the Pol II pausing complex through its interaction with the NELF complex and thereby directly influences hsp transcription during the stress response. Methods Fly stocks Loss of function was achieved by expressing dAxud1 RNAi through the GAL4/UAS system using the Vienna stock V26479. For the experiment for polytene chromosomes, we used flies with the genotype nub > Gal4/Y ; UAS-dAxud1::GFP/+. For DamID-seq experiments (TaDa-seq), the genotypes were nub >Gal4/Y; UAS-mCherry.Dam and nub > Gal4 ; UAS-mCherry.Dam::dAxud1 , both generated from the plasmid pUAST-attb-LT3-Dam donated by the Andrea Brand Laboratory [41]. The Dam fusion plasmids were injected with pBS130 that express the phi -131 integrase. The genotype of injected flies was P{y[+t7.7]=CaryP}attP2 (Bloomington stock 8622), in accordance with methods in the previous reference [18]. For HA-NELF-B experiments in imaginal wing discs protein extract, stock F003904 from FlyORF was used for expression under the UAS/Gal4 system. Lifespan assay For every replicate, between 80-100 flies were collected no longer than 48 hours from pupae eclosion and separated into males and females. Three replicates for every condition were used, and dead and live flies were counted every day. The control condition was at 29°C, with the flies changed from the tubes every two days. The heat shock condition was at 37°C for 30 minutes every day. Immunofluorescence in polytene chromosomes The squash protocol and antibody staining for polytene chromosomes were performed according to the Johanssen protocol [79]. The GFP antibody used was Cell Signaling (code D5.1), dilution 1:250, and the RNA Polymerase IIo antibody (full phosphorylated CTD) used was Abcam (ab5408), dilution 1:1000. DAPI stain was used at 1:10000. The glands were dissected from larvae with nub>Gal4; UAS-dAxud1::GFP genotype. TaDa-seq In total, 100 to 120 imaginal wing discs per sample were dissected from third instar larvae ( nub>Gal4 ; UAS-mCherry.Dam & nub>Gal4 ; UAS-mCherry.Dam::dAxud1). The larvae grew for 7 days at 17°C and were then placed in 29°C for 24 hours, to the third instar stage. The discs were dissected and placed in PBS 1X, and the DNA was extracted according to the Southall protocol [41]. Accordingly, two biological replicates for every condition (Dam and Dam-DAxud1) were used to generate libraries, sequenced in Illumina Hiseq-2500, single end 100 bp. The data were processed with the pipeline designed and published by Marshall [80] using the dm6 Drosophila genome version (release 6.22). The replicates were merged and data was generated in bedgraph format, according with the pipeline . The GFF files with the coordinates of significant peaks were visualized with the WashU epigenome browser. These data were processed with HOMER [35] for motif discovery and finding. The annotation was done with T-Gene from Meme-Suite [81]. The metagene and heatmap profiles were generated with the deeptools toolbox platform [82]. RNA extraction and qPCR RNA extraction was carried out with the TriZol reagent following manufacturer instructions. The RNA was resuspended in water and the RT reaction was performed with the iScript® kit from BioRad. The qPCR mix reactions were done with Brilliant II SYBR® reaction solution with 10 minutes of initial denaturation and 40 cycles in the following sequence: 30 seconds at 90°C; 30 seconds at 60°C; and 30 seconds at 72°C. The primers used were as follows: dAxud1 Fw 5’- AGGGGACCACCAGCCTAAC-3’, Rv 5’-GGTTCGCTCTGATTATCCTTGTG-3’; hsp26 Fw 5’-ATGCCCACGATCTGTTCCATC-3’, Rv 5’-GTACGCGAATAACGACGAC; hsp67Bc Fw 5’-GACTCCCCGGACTCCATGTA, Rv 5’-GACCCAGGGTGTGCAAATCAA-3’; hsp70 for B paralogs (for Ba, Bb and Bc) Fw 5’- AAGAACCTCAAGGGTGAGCG-3’, Rv 5’-CGAACAGAGATCCCTCGTCG-3’. For hsp70B 5’-UTR Fw 5’- GCTAAGCAAATAAACAAGCGCAG-3’, Rv 5’-CAGTTGATTTACTTGGTTGCTGGT-3’ ; hsp70B 3’-UTR Fw 5’-GAGGATTTGGCGGCTACTCT-3’, Rv 5’-TTTAAAAACTTAAGCCAGGAACTGA-3’; actin-42A (as normalizer) Fw 5’-GCGTCGGTCAATTCAATCTT-3’, Rv 5’-AAGCTGCAACCTCTTCGTCA-3’ Chip-PCR protocol Chromatin immunoprecipitation was performed from salivary glands expressing DAxud1-GFP or dAxud1 RNAi in third instar larvae growth at 25°C using the driver nubbin >Gal4 (X). Larvae from the heat shock condition were placed at 37°C for twenty minutes. The glands were dissected and processed according to Ghosh et al.’s protocol [54] for three biological replicates per condition. Ten pairs of salivary glands from third-instar larvae were incubated for five minutes on ice in 100 μl of 1% formaldehyde in PBS 1X and then at RT for 7 minutes. The cross-linking reaction was quenched by adding 2.5 M glycine to a final concentration of 125 mM; then the glands were placed on ice for two minutes. They were centrifuged at 900 × g for two minutes at 4°C, and the supernatant was removed. A 100-μl volume of sonication buffer (20 mM Tris [pH 8.0], 0.5% SDS, 2 mM EDTA, 0.5 mM EGTA, 0.5 mM phenylmethylsulphonyl fluoride [PMSF]) and 1 μl Halt protease inhibitor cocktail (Thermo Fisher cat. 78430) were added to the glands to incubate at room temperature for ten minutes and then on ice for 10 minutes. The glands were vigorously shaken for ten minutes and homogenized with a small pestle. Lysates were sonicated at 4°C in Omniruptor-4000 at 100% power, 90% pulse for fifteen minutes to shear the DNA to an average fragment size of 400 bp. The lysate was clarified by centrifugation at 14,000 × g for 7 minutes. In total, 45 μl per lysate were used for each immunoprecipitation assay (mock and IP), with 10 μl for input. Immunoprecipitation was performed using 1 μl of anti-GFP antibody (Abcam ab290, 5µg for chromatin from teen pairs of salivary glands) or 30 μl of RNA Polymerase II antibody (Santa Cruz, 8WG16, 6µg for chromatin from ten pairs of salivary glands), diluted half in IP Buffer (50mM Tris-HCl pH 8; 100mM NaCl; 2mM EDTA; 1mM EDTA, 1% NP40, 1X Thermo Fisher protease cocktail), and immunoprecipitated with 30 μl of protein-A Dynabeads® solution, with the same quantity for the mock assay. Beads were washed twice in sequence with CHIP 1 Buffer (IP buffer + 1% Na-deoxycholate), CHIP 2 Buffer (IP Elution buffer + 1% Na-deoxycholate + 500mM NaCl), and CHIP 3 Buffer (IP Elution buffer + 1% Na-deoxycholate + 270mM LiCl), then washed twice in Tris-HCl 10 mM pH 8. Finally, elution was performed with elution buffer (NaHCO 3 pH 8.8 0.1M; SDS 1%), added to the precipitated samples until 200 μl was reached (same as input). Decrosslinking was performed by adding 10 μl of NaCl 4M and 0.5 μl of fungal proteinase K (Thermo Fisher), with six hours of incubation at 65°C. The samples were purified and eluted with MicroChip Diapure Columns. Samples were analyzed with qPCR. The primers used were as follows: hsp70Aa TSS (-800,-680) Fw 5’-AACAACAAATTCCAAGTTTGCAC-3’, Rv 5’-CACGGTTTTTGCATATGCTTT-3’; hsp70B paralogs 5’-UTR Fw 5’- GCTAAGCAAATAAACAAGCGCAG-3’, Rv 5’-CAGTTGATTTACTTGGTTGCTGGT-3’; hsp70 common AGAGTG motif Fw 5’-TTCTCTGGCCGTTATTCGTT-3’, Rv 5’-TCGAACCAACGAGAGCAGTA-3’; hsp70 TSS Fw 5’-CGACATACTGCTCTCGTTGG-3’, Rv 5’-CAGCTGCGCTTGTTTATTTG-3’; hsp70B 3’ end Fw 5’-AATGGAATCCTGAACGTCAGC-3’, Rv 5’-CACATTGAAGACGTAGCTCTCC-3’. Additional primers for 5’ and 3’ were extracted from reference [16] and were named CHIP-A and CHIP-B, respectively. Co-immunoprecipitation Coexpression of DAxud1-GFP and NELF-B-HA on the imaginal wing disc was achieved with driver nub>Gal4. In total, 150 wing discs were dissected on PBS 1X, then precipitated and resuspended in RIPA buffer. Protein extracts were immunoprecipitated with GFP antibody (Abcam ab290) with O.N incubation at 4°C, then isolated with Dynabeads® Protein-A. Western blot assay was performed in cells with INPUT (total protein), NB (non bound fraction, supernatant post Dynabeads isolation), and protein precipitated with Dynabeads. Page ruler 2166 (thermoscientific) was used as a weight marker. Results Distribution of DAxud1 across the genome According to results from previous reports, vertebrate DAxud1 orthologs possess transcription factor-like features, including DNA binding, acidic and transcription transactivation domains [20, 21, 28]. They have also been detected bound to specific promoters related with neural crest differentiation [28]. Further, vertebrate and invertebrate orthologs of DAxud1 are related to the stress response and cancer through induction of apoptosis [21, 23], although, until now, this process has not been linked to its putative transcriptional function. Considering this information, we asked which groups or categories of genes are transcriptionally regulated and/or bind DAxud1 in Drosophila , using this information as a proxy to get a better perspective on the function of dAxud1 . For this purpose, we performed a TaDa-Seq (DamID-seq) experiment using the method described by Southall et al. [29], expressing the Dam-DAxud1 fusion protein to explore the loci in which DAxud1 could reside. This was performed in imaginal wing disc tissue using the nub >Gal4 driver. Two replicates were generated per condition (third instar larvae imaginal wing disc, with Dam as a control or Dam-DAxud1). The aligned sequences were analyzed with the findpeaks script, generated by Marshall Owen [30]. Using this method, we identified 1811 significant peaks across the Drosophila genome, representing zones where Dam-DAxud1 has a stronger, more stable positioning than the Dam signal. As shown in Fig. 1A, the Dam-DAxud1 distribution was mainly intronic, with 52.3% of peaks localizing to these gene regions. The first intron accounted for more than half of all instances of intron localization, representing 25.7% of overall genome occupancy. The other main significant locations were the intergenic and the proximal promoter regions (-1000 to +1). These regions are also the most enriched in transcription factor binding sites [31, 32]. To better understand the distribution of DAxud1 across the genome and in specific genes, we used the information from significant peaks to create a metagene profile, using gene bodies as a reference scaled to 1000 bp. The resulting metagene profile (Fig. 1B) shows a robust average signal at the promoter zone and a maximum within the 5’-half of gene bodies, confirming the previously described distribution (Fig. 1A) and suggesting a role in activating or poised gene transcription. Several reports describe DAxud1 vertebrate orthologs as stress response genes [24, 26, 33, 34]. Therefore, we searched for significant DAxud1 binding on stress response gene loci (Additional Table S1) and found heat shock response genes as recurrent targets. Further, we generated alignment charts to visualize the distribution of peaks within this group of genes (Fig. 1C). These graphics clearly reveal that the distribution of Dam-DAxud1 specifically on hsp genes has a highly similar pattern to that observed in the metagene profile. The intergenic, promoter, and gene-body (intron and exon) peaks were sorted and annotated for nearby genes using HOMER [35]. Additional Table S2 shows the most relevant enrichment analysis results. Within the “biological processes” category, the two major enriched groups are related to neural development, but no groups show a significative score for stress response or tolerance. One pathway category identified involves Wnt signaling, concordant with reports that DAxud1 orthologs in mammals are linked to this signaling pathway [19, 36, 37] and apoptosis [21, 38]. Potential DAxud1 DNA binding motifs After having found genes in which Dam-DAxud1 is enriched, we searched for recurrent binding motifs within the Dam-DAxud1 peaks. The sequence peaks had a size between 50-5000 bp and were extracted for motif find analysis from HOMER [39]. As shown in Fig. 1D, the most recurrent DNA motif (TACATACATA), present in 1019 sequence peaks out of 1811, is novel compared with a previous report on Axud1 vertebrate orthologs [28], possibly due to the wide-range of sequences that the TaDa-seq experiment provides [40–42], but this information could provide insights from the chromatin context surrounding DAxud1. Using the matrix sequence, we conducted another analysis using this data and the GOMO platform, which provides information about the gene ontology (GO) of the promoter regions, using a range between -1000 and +200 bp from the TSS for each gene in the Drosophila melanogaster genome. For the TACATACATA motif, the most related GO found is “Heat shock mediated polytene chromosome puffing” (GO:0035080). This motif also found a match in the Drosophila Topoisomerase 2 gene, which encodes a protein that removes supercoils in chromatin to facilitate transcription, with an essential role in RNA polymerase pausing-release in hsps genes in the fast stress response [43, 44]. These analyses suggest that DAxud1 might play a role in the transcription of heat shock and other stress-induced genes. To confirm this last hypothesis, we performed chromatin immunoprecipitation on the hsp70 promoter with DAxud1-GFP as bait, followed by a CHIP-PCR with primers flanking the TACATACATA motif. The result confirms the presence of DAxud1 at this locus (Additional Fig. S9), reinforcing the idea that it regulates the stress response by modulating hsp gene expression and also that the regions discovered using TaDa-seq are bona fide DAxud1 binding sites representing potentially functional interactions. DAxud1 regulates thermotolerance and lifespan in Drosophila Our results reveal the presence of DAxud1 bound to hsp genes (Fig. 1; Additional Fig. S9), suggesting a possible role in stress tolerance, a hypothesis supported by studies on Axud1 orthologs in mammals [23, 45]. Although “stress response” is not the major GO category of DAxud1 genome occupancy (Additional Table S2), we tested whether flies with a global knockdown of this gene exhibit an altered stress response or aberrant hsps expression, given its presence in this class of genes. For this, we used flies expressing a DAxud1 RNAi construct (Vienna stock 26479, UAS-IR-DAxud1) in all tissues using the Tubulin-Gal4 driver. Flies overexpressing DAxud1 cannot be assessed in this type of experiment since ubiquitous expression of this protein is lethal [21]. We established adult flies at 29°C for an optimal expression of the UAS RNAi construct; then, control and experimental adult flies were exposed to a heat shock of 37°C for a half-hour every day as the thermal stress condition, and their survival was measured. Knockdown of DAxud1 leads to a diminished lifespan of adult flies compared to the control genotype under the same stress condition (Fig. 2A-2B). Curiously, under control conditions (no heat shock), DAxud1 knockdown animals have a longer lifespan compared to control animals for both males and females (Figs. 2A and 2B). To confirm whether there is a fluctuation in hsp gene expression, and to relate this to the observed phenotypes, we performed qPCR analysis for different hsp genes in imaginal wing discs and salivary glands, expressing the RNAi construct or DAxud1-GFP version using the nubbin -Gal4 driver. This allowed us to study the effects of diminished DAxud1 levels in salivary glands and imaginal wing discs, tissues that have cells in an endo replication state (salivary glands) and in a mitotic state (imaginal wing discs). Figs. 2C, 2D, and 2E show the qPCR results for three hsp genes ( hsp70, hsp26, hsp67 ). In these experiments, we can appreciate that DAxud1 knockdown increases the expression of hsp genes in control conditions, but does not potentiate the transcriptional heat shock response in either tissue. Notably, the rise of hsp mRNA expression at the control temperature can explain the longer lifespan since there is evidence that hsp overexpression generates this effect [46]. On the other hand, DAxud1 overexpression reduces hsp basal expression at the control temperature, perhaps not inhibiting the hsp stress response but generating a milder induction (Figs. 2C-E). These results, especially those observed after DAxud1 knockdown, are in conflict with our first hypothesis in which DAxud1 could function as a transcription factor that facilitates hsp expression. Therefore, we then set out to further clarify the role of DAxud1 on the hsp70 locus . DAxud1 exhibits a widespread presence across the genome, whereas heat shock induces its recruitment to the hsp70 locus Further examination of the results of hsp gene expression in DAxud1 knockdown conditions (Fig. 2), suggests that DAxud1 may participate directly in hsp transcription, possibly in a repressive manner. Previous studies have reported that DAxud1 orthologs (CSRNP in mouse) have transcription factor features, including a C-terminal trans-activating domain and a DNA binding motif [47]. However, there is no evidence indicating which type of regulatory factor it is, and its dynamics within the nuclear structure. Drosophila melanogaster’s polytene chromosomes represent a suitable model to answer these questions. As there is no available antibody against DAxud1, we used the GFP-tagged version of DAxud1 [21] (DAxud1-GFP) to perform immunofluorescence using an anti-GFP antibody. RNA Polymerase IIo (Hyper-phosphorylated) was used as a positive control for a chromosome attached protein and for detecting transcriptionally active loci. As seen in Fig. 3A, DAxud1-GFP on the polytene chromosome exhibits a pattern with widespread RNA Pol IIo co-distribution, which means the specific function on hsp genes may be distinct from other functions of DAxud1 in gene expression. The hsp70 locus (Fig. 3B) has been extensively analyzed with regards to its chromatin rearrangements under stress conditions. Its documented that during stress conditions such as heat shock, cytological zones known as chromatin puffs appear in hsp70 region, clearly visible with confocal microscopy, a reflection of chromatin opening for transcription and transcription factor recruitment [48, 49]. We evaluated chromatin and DAxud1-GFP dynamics on this locus (at both 17°C and 37°C, Fig. 3C) and observed that DAxud1-GFP localizes to the 87A-87B locus together with RNA Polymerase IIo after heat shock treatment. Transcription factors that relocate in this way are identified mostly as activators, for instance HSF, p-TEFb [2], and DSIF [16]. From these results, we reasoned that DAxud1 is a chromatin element that relocates to the hsp70 locus during heat stress, but it does not interfere with the transcriptional induction since hsps mRNA levels rise by the same magnitude as in control, under a DAxud1 overexpression condition (Fig. 2C-E). To confirm these results, we performed chromatin immunoprecipitation for DAxud1-GFP (Fig. 3D) to detect the presence of this factor on the hsp70 paralogs hsp70Ab and hsp70Bb/Bc , or loci 87A and 87C, respectively. According to those results, DAxud1-GFP has a clear presence on hsp70 promoters as do other heat shock transcription factors [4]. Interestingly, DAxud1-GFP did not exhibit enrichment after oxidative stress induction by peroxide (Fig. 3D), a condition that increases hsp mRNA levels [50], indicating that the involvement of DAxud1 in the heat shock response is different to the mechanism evoked during oxidative stress. Together, the data obtained here, have described the genome-wide distribution of DAxud1, its influence on hsps transcription, and its effects on lifespan. However, how does DAxud1 exert its function and what is its relevance for cell physiology? In the hsp70 gene, there is a special chromatin configuration, known as the pausing complex, in which the activated RNA Polymerase II pauses the transcription process at positions +30 to +50 bp of the gene, and remains stalled as a “poised polymerase” [44]. This poised polymerase forms a complex with two main components, NELF and DSIF, whose role is to maintain the polymerase in the paused state. Transcription resumes when the p-TEFb complex phosphorylates NELF and DSIF, releasing the active polymerase for resumption of elongation [44]. hsp70 genes have been used as transcriptional pausing models because the resumption of transcription occurs quickly when heat shock or stress signals reach the cell [44]. Considering the data on the position and behavior of DAxud1 on the hsp70 gene after heat shock and that elements of the pausing complex on hsp70 exhibit a similar pattern of reorganization on chromatin during heat shock [4, 18], we conjectured that DAxud1 might interact with the pausing complex. This hypothesis is also supported by RNA-seq evidence (Additional Figs. S3 and S4), in which DAxud1 overexpression increases hsp70 RNA synthesis only at the 5’ end of the gene, with an incomplete synthesis of hsp70 mRNAs, specifically on hsp70B paralogs, which display high reads only between +1 and +100 region. This effect was not detected in qPCRs performed in Figs. 2C-E because the primers were designed to detect the 3’ side. In addition, this role in pausing is also supported by the evidence of double-hybrid experiments with Drosophila proteins, which demonstrated a physical interaction between DAxud1 and Spt5 [51], a member of the DSIF complex, with a pivotal role in the RNA Polymerase transcriptional pausing and elongation complex [52]. Notably, the DAxud1 overexpression wing phenotype is reversed in a heterozygous spt5 mutant background (Additional Fig. S6). Another component of the pausing machinery is the NELF complex (NELF-A, NELF-B, NELF-C/D, NELF-E), a stabilizer factor of the stalled polymerase in basal conditions. These components dissociate from the polymerase when transcription is resumed after heat shock [18], and the loss of function of some of its components generates higher levels of hsp70 [53], which is similar to our observations in the case of DAxud1 knockdown (Fig. 2C). Considering this information, we tested whether DAxud1 influences RNA Polymerase II location within the hsp70 gene body [53]. Chromatin immunoprecipitation (CHIP) for RNA Pol II was performed under different DAxud1 levels of expression (Fig. 4A). In the RNA Pol II CHIP, we find that, in the case of DAxud1 knockdown, its levels decrease uniformly in the hsp70 gene locus , and there is no change in the 5’/3’ rate of RNA Pol II compared to the control and heat shock conditions (Fig. 4B). On the other hand, DAxud1 overexpression generates an increase of this ratio, more so than in the control condition, due to an increase of RNA Pol II on 5’ end of the hsp70 gene body (Fig 4A). To find out whether DAxud1 interacts with pausing factors, we performed a co-immunoprecipitation assay in imaginal wing disc expressing NELF-B-HA and DAxud1-GFP. In this assay (Fig. 4C), a fraction of NELF-B-HA was co-immunoprecipitated with DAxud1-GFP, indicating a physical interaction between these two proteins, suggesting they could act in the same process in hsp70 expression. The information provided by the physical and genetic interaction of DAxud1 with NELF-B and Spt5, respectively, strongly suggests that it could be influencing RNA Polymerase II dynamics on the hsp70 gene. To confirm this, we performed a CHIP for RNA Polymerase on the 5’ and 3’ regions of the hsp70 gene, under DAxud1 overexpression and knockdown conditions. The results, presented in Fig. 4, reveal that DAxud1 overexpression generates an apparent pausing effect on RNA Polymerase II distribution (Fig. 4), coincident with extensive upregulation of the 5’ RNA levels of hsp70 (Additional Figs. S3 and S4). Promoter-proximal pausing in hsp genes allows a rapid expression after heat shock induction due to pausing elements as DSIF/NELF complexes. Also, this rapid expression has a shut-off mechanism to restore hsp70 mRNA levels after the heat shock condition is reverted (heat shock recovery), and there is evidence that supports the role of NELF complex in the shut-off mechanism [54]. We tested the physical interaction between NELF-B and DAxud1 (Fig. 4C) as well as its genetic interaction (Additional Fig. S7), allowing us to further clarify the role of DAxud1 in hsp70 shut-down levels during heat shock recovery. To evaluate the effect of DAxud1 in this process, we performed qPCR for hsp70 mRNA 40’ after a 20’ heat shock (37°), as shown in Fig. 4D, in DAxud1 knockdown salivary glands. It is possible to observe that hsp70 induction has an apparently slower turn-over than in control animals (Fig. 4D), suggesting that DAxud1 plays a role in shutting-down hsp70 locus induction, and confirming DAxud1 has a synergistic function with NELF-B in the transcriptional pausing of hsp70 as well as the complex stabilization for further heat shock stimulus. Discussion Role of DAxud1 in tissue homeostasis through hsp regulation In the first part of this study, we searched for the most frequent sites where DAxud1 locates in the Drosophila genome, and upon finding hsp -type genes in these sites we performed thermotolerance analysis under DAxud1 knockdown in order to test the physiological relevance of this gene. We further studied hsp expression during DAxud1 knockdown to mechanistically link these cellular players. In lifespan assays with DAxud1 knockdown (Figs. 2A-B), we found it generates adults with extended lifespan in control conditions, but reduced thermal resistance compared with control genotype animals, with no other change in developmental timing or apparent alterations (Fig. 1; Additional Fig. S2). Interestingly, in DAxud1 knockdown larvae raised at 29°C (no heat shock), the organisms exhibit an increase in the level of hsp70 expression in salivary glands and imaginal wing discs, compared to the control animals subjected to the same temperature. In Drosophila , the effect of hsps overexpression is well documented, in that one of its effects is precisely an extended lifespan due to its cytoprotective and anti-apoptotic effect [46, 55], explaining why DAxud1 knockdown extends lifespan in control temperature, compared with animals only expressing the Gal4 driver (Figs. 2A-B). On the other hand, when animals with DAxud1 knockdown are exposed to daily thermal stress, lifespan shortens, situation that did not occur with control animals in which lifespan remained within the same range, as in control temperature, during daily thermal stress (Figs. 2A-B). This seems paradoxical since animals with higher levels of Hsps would be expected to maintain enhanced thermotolerance. However, the diminished tolerance to stress in DAxud1 knockdown could be explained alternatively by the reported pro-apoptotic DAxud1 function through an activating role on JNK signaling in a DAxud1 overexpression background [21, 26]. Thus, in the context of DAxud1 knockdown, the JNK signaling pathway could be less strongly activated, allowing damaged cells to survive though interfering with development and regenerative processes, which are known consequences of impaired JNK signaling [56, 57]. This last suggestion is supported by the occupancy of Dam-DAxud1 on the loci of genes associated with the GO- KEEG apoptosis pathway, related to JNK (Additional Tables S2 and S5). Regarding this point, the pro-apoptotic activity of DAxud1 is associated with control of tumorigenesis [19, 21] and cancers with poor prognosis [58], which may be due not only to decreased pro-apoptotic activity resulting from lower Axud1 levels, but also to increased hsp gene expression. hsp overexpression strongly correlates with cancer cell progression and poor prognosis [59, 60], so this condition of certain cancer cells may be due to low levels of Axud1, a condition also associated with aggressive cancers, as previously described [19, 21, 58]. This function in tissue homeostasis through transcription modulation could be one of the reasons DAxud1 is conserved among metazoans. However, the reduced thermotolerance phenotype would not be directly related to the expression of hsps , despite the fact that DAxud1 exhibits recurrence in these genes (Fig. 1) and at the same time the strong effect on their expression generated by the knockdown of DAxud1. These observations prompted us to investigate the relationship between hsps , widely conserved in all kingdoms, and DAxud1, which only has orthologs in metazoans, probably participating in a type of modulation on hsps expression exclusive to metazoans. This led us to focus on its effects at the chromatin level. DAxud1 function in chromatin as pausing factor on the hsp70 gene Although dAxud1 mRNA levels do not change significantly during heat shock (Additional Fig. S8), the DAxud1 protein relocates to the hsp70 loci . The localization of DAxud1 to chromatin, observed in Fig. 3C, can mediated by putative phosphorylation sites [21] or by the DAxud1 cysteine-rich region [20, 21]. These cysteine-rich regions have been shown to act as modulators of protein conformation, in which heat or oxidative stress can alter the disulfide bonds, changing the conformational state of the protein. This phenomenon is well documented [61, 62] in proteins with cysteine-rich regions, including the DAxud1 CSRNP family [20]. Thus, changes in DAxud1 protein conformation induced by post-translational modifications could be key in the regulation of its localization and its effect on gene expression. Transcriptional pausing is a state in which the RNA Polymerase initiates transcription but remains stalled on the first 30-50 base pairs of the gene, forming a complex with the DSIF (Spt5/Stp4 proteins) and NELF complexes [16, 48, 63]. The pausing is released when internal or external cell signals activate p-TEFb, which phosphorylates RNA Polymerase II itself as well as the DSIF, and NELF complexes which, in turn, disengage from the RNA Polymerase [18, 64, 65]. Also, there is evidence that pausing relies on the chromatin state, in which NELF exerts it pausing effect by depleting H3K4me3 histones and, therefore, when NELF dissociates from the pausing complex, the polymerase elongates in a favorable transcriptional context [53]. This mechanism is found in fast-response genes like hsp genes or in innate immune response genes [66]. Strikingly, overexpression of DAxud1 generates an increase in hsp70B transcription only on its 5’ end of the transcript sequence (Additional Figs. S3 and S4). This increment of 5’ sequences of the hsp70B genes indicates that there is an increase of mRNA synthesis, from approximately +1 to +80 bp. We interpret that DAxud1 retains the pausing complex completely stalled, generating an enrichment of 5’ hsp70 mRNA. Whatever the explanation is at this level, the effect of DAxud1 on this particular hsp gene expression could involve the transcriptional pausing mechanism, in which DAxud1 might interact with the pausing complex. The possible role of DAxud1 in the pausing complex is also supported by the co-immunoprecipitation of DAxud1-GFP with NELF-B (Fig. 4C). NELF-B acts as a component of the pausing complex, stabilizing and pausing RNA Polymerase II in a hypo-phosphorylated CTD state, and dissociates from the complex under heat shock or another signal [54]. Additionally, data concerning Spt5, another component of the pausing complex (DSIF complex), further supports this observation: Biogrid [51], a database that documents physical interactions between proteins, reported an interaction of DAxud1 with Spt5 [4] and, importantly, the DAxud1 overexpression phenotype is partially reverted in a spt5 heterozygous mutant background (Additional Fig. S6). Furthermore, Spt5 interacts with NELF-B in the pausing complex [67], and stabilizes it synergistically with the NELF complex, maintaining the RNA Polymerase stalled. This pausing complex keeps stable until p-TEFb kinase phosphorylates both complexes, with dissociation of NELF complex [68] and DSIF (Stp5/6) remaining as a component of RNA Polymerase holoenzyme, acting as an elongation factor. In this context, DAxud1 may play a role in pausing the stabilization complex, since DAxud1 knockdown causes an increase in the transcription rate of hsp70 mRNA, while its overexpression causes stalled transcription at the 5’ region, perhaps instigating the formation of the pausing complex (Fig. 4) in the same way as NELF does [53], but not interfering with RNA Polymerase II pausing release in heat shock conditions. This interaction could be conserved in other instances, such as in the regulation of the expression of the MMP1 gene in humans, which requires the ortholog of DAxud1 (CSRNP1) for its transcription upon cytokine activation [69]. More importantly, there is detailed evidence that MMP1 is a transcriptionally paused gene, in which the NELF/DSIF pausing complex is necessary for its appropriate expression in response to immune-activating signals [66]. Also, in myeloid cells, the NELF complex and Spt5 are essential for maintaining the repression of pro-apoptotic genes during myeloid development in zebrafish. In this case, they rely on their role in RNA Pol II pausing of pro-apoptotic genes [70] and, accordingly, Spt5’s deficiency generates a loss of the myeloid line, due to early apoptosis of hematopoietic stem cells, the same phenotype observed after Csrnp1 reduction in morphant zebrafish larvae [71]. Curiously, in imaginal wing discs and, subsequently, in the adult wing, the overexpression of NELF-B reverts the DAxud1 overexpression phenotype, while NELF-B knockdown exacerbates the DAxud1 overexpression phenotype (Additional Fig. S7). These results suggest an antagonistic role of these two proteins on gene expression, though not necessarily related to transcription of hsp genes, because the knockdown of both genes separately provokes a similar effect in the increase of hsp transcription. Instead, such an effect could be found in genes related to apoptosis, in which NELF-B functions as a repressive element, as previously described [53]. This is also similar to the interaction between DSIF and NELF, both are necessary for RNA Pol II pausing, but after p-TEFb activation (through heat or morphogen signals) their roles are quite opposite, in which the DSIF elements are necessary for RNA Pol II elongation while NELF dissociates completely from the transcription complex [18]. The relation between pausing and DAxud1 motif founded with TaDa-seq In support of the pausing-release hypothesis, we found that the most enriched motif found in TaDa-seq peaks is “TACATACATACA.” This sequence strongly matches the binding site of Topoisomerase 2 (Top2), a chromatin modifying protein (Fig. 1D). There is evidence that Top2 plays a role in pausing-release and in the hsp70 fast transcriptional response in Drosophila Kc cells [72]. In addition, GOMO analysis shows that this motif is present in promoters of genes that belong to the “heat shock-mediated polytene chromosome puff” (GO:0035080) and “response to hypoxia” (GO:0001666) GO categories. Additionally, there is evidence that Topoisomerase 2, due to its structural chromatin function, plays a role in the cell cycle, in which inhibition or missteps in the arrangement of this protein lead to cell cycle arrest in G2-M [73]. The same cellular phenotype was previously reported when DAxud1 was overexpressed [21], raising the possibility that these proteins interact at the chromatin topology level, likely interfering with Topoisomerase 2 positioning before the G2-M checkpoint. However, more analysis is required to support this conjecture. In previous studies, the vertebrate orthologues of DAxud1 were associated with the binding motif AGAGTS [28, 69, 74]. Curiously, with HOMER ( denovo motif discovery function), no similar sequence was found. We thus used HOMER to search for the occurrence of the AGAGTS motif among the peak sequences; we found AGATGS 4657 times in 1407 peaks among the 1811 significant sequences. In total, 895 genes were identified (Additional Table S4). Again, pathways like Wnt and Apoptosis are present in this group (Additional Tables S2 and S5), most of those belonging to the apoptosis group are part of JNK signaling, including msn , puc , rpr , and hid , a finding which is in agreement with bsk and puc upregulation upon DAxud1 overexpression (Additional Fig. S5, [21]). This result indicates that DAxud1 acts at different levels on transcriptional regulation to balance cell physiology. With regard to its physiological role, low levels of human Axud1 expression are related to tumor development [19]. Considering the role of Axud1 in the stress response, this condition could predispose cells to undergo more damage in the first stage of tumor development due to the loss of hsp gene regulation, leading to higher levels of Hsp proteins, thereby buffering damage signals and likely inhibiting the apoptotic process and extending the lifespan of impaired cells. Later on, in the next stage, Axud1 could be necessary at higher expression levels to initiate apoptosis through JNK signaling [21], though this Axud1 upregulation does not occur in tumorigenesis [19]. Developmental role of Axud1 A previous study showed that the DAxud1 vertebrate ortholog, Axud1, acts as a transcription factor, located only in the promoter region and related to the positive transcriptional activity of the Pax7 and Msx1 transcription factors, driving the development of neural crest by Wnt signaling [28]. According to this study, the presence of Axud1 on these promoters depends on the AGAGTS binding site, previously described for human and mouse orthologs [11], from in vitro experiments. The AGAGTS motif was found in the TaDa-seq data, with a frequency of 4657 times out of 1407 unique peak sequences. The AGAGTS sequence is present 88,167 times in the whole Drosophila genome and, therefore, the Dam-DAxud1 peak sequences comprise 5.28% of the whole genome motif’s abundance. Out of those 1400 peaks, gene ontology analysis mainly points towards morphogenesis, axon guidance, and Wnt signaling (Additional Tables S2 and S5). This information is in line with previous work on vertebrate Axud1 [19, 28, 75] and our data in which Axin and Wnt are downregulated upon DAxud1 overexpression (Additional Fig. S5), which could indicate that DAxud1 is part of a feedback loop in the Wnt pathway. Mainly, the data from other studies supports the role of Axud1 in neural development as previously described, in which vertebrate homologs of DAxud1 could exert a pro-neural function in both CNS and neural crest development, and this depends on wnt signaling [28, 76]. Therefore, this regulatory network is conserved in animals with different roles in development. This function in neural structures is a hallmark of DAxud1 orthologs because our previous data showed that zebrafish Csrnp1a is highly expressed in the growing brain [77] and in the ventral nerve cord in Drosophila during germ band elongation stage [21], making this protein a regulator of a specific process more than a general factor of animal development, and not classifiable in a binary category of transcription activator/repressor. Overall, our study raises new insights about DAxud1 function, demonstrating new chromatin features related to the stress response. We propose this is mainly supported on the stabilization of the pausing complex, maintaining a chromatin configuration for rapid hsp expression under stress conditions in the first stage of the thermal stress response. This function might not represent the complete chromatin regulatory function of DAxud1, which could play a role in JNK signaling in long-term stress to promote apoptosis, and an additional role in neurogenesis as a mediator of Wnt signaling, as described in vertebrates [78]. Declarations Ethics approval: Our laboratory and the manage of animals were approved to perform the experiments exposed in this article. This was approved by Bioethics committee of the Faculty of Sciences (Universidad de Chile), chaired by Dr. Marco Méndez from 2013. Competing interests: No potential competing interest was reported by the authors. Author’s contributions: JM Zúñiga-Hernández was the main executor, writer of the original draft, designer of the experiments, as well as the data analysis, including NGS data. Also, JM Zúñiga-Hernández is the corresponding author. Meneses C contributed to facilitate the access to Illumina sequencing platform, and support experimental design. Bastías M was the technician in charge to process the libraries for the Illumina sequencing. Allende ML contributed to the editing and review of the manuscript, funding and provider of resources. Glavic A was the supervisor and reviewer of the manuscript, funding provider, as well as the manager of the initial idea for this research. Author information: Jorge Zúñiga Hernández, PhD (Zúñiga-Hernández JM): Corresponding author, works as Postdoctoral researcher in Center for Genome Regulation and Laboratory of Developmental Biology, Department of Biology, Faculty of Sciences, University of Chile. Santiago, Chile. Mail: [email protected] Claudio Meneses, PhD (Meneses C): Investigator in Plant Biotechnology Center, University Andres Bello. Santiago, Chile. Mail: [email protected] Macarena Bastías (Bastías M): Technician in in Plant Biotechnology Center, University Andres Bello. Santiago, Chile. Mail: [email protected] Miguel Allende, PhD (Allende ML): Director of Center for Genome Regulation and Laboratory of Developmental Biology, Department of Biology, Faculty of Sciences, University of Chile. Santiago, Chile. Mail: [email protected] Alvaro Glavic, PhD (Glavic A): Co investigator in Center for Genome Regulation and Laboratory of Developmental Biology, Department of Biology, Faculty of Sciences, University of Chile. Santiago, Chile. Mail: [email protected] Consent of participation: No humans were part of the experiments exposed in this article. Consent of publication: No humans were part of the experiments exposed in this article. Funding: Anillo ACT1401 / Center for Genome Regulation (FONDAP 15200002) Availability of data: TaDa-seq are available in Sequence Read Archive (SRA), ID: PRJNA776616 Acknowledgements: We thank to Andrea Brand for provide plasmid pUAST-attb-LT3-Dam. Also acknowledgments to FONDECYT (ANID) for Doctoral Fellowship 21110721. References Richter K, Haslbeck M, Buchner J. The Heat Shock Response: Life on the Verge of Death. Mol Cell. 2010;40:253–66. doi: 10.1016/j.molcel.2010.10.006 . Boehm AK, Saunders A, Werner J, Lis JT. Transcription Factor and Polymerase Recruitment, Modification, and Movement on dhsp70 In Vivo in the Minutes following Heat Shock. Mol Cell Biol. 2003;23:7628–37. Mayer MP, Bukau B. Hsp70 chaperones: Cellular functions and molecular mechanism. Cell Mol Life Sci. 2005;62:670–84. Andrulis ED. High-resolution localization of Drosophila Spt5 and Spt6 at heat shock genes in vivo: roles in promoter proximal pausing and transcription elongation. Genes Dev. 2000;14:2635–49. doi: 10.1101/gad.844200 . O’Brien T, Lis JT. Rapid changes in Drosophila transcription after an instantaneous heat shock. Mol Cell Biol. 1993;13:3456–63. San Gil R, Ooi L, Yerbury JJ, Ecroyd H. The heat shock response in neurons and astroglia and its role in neurodegenerative diseases. Mol Neurodegener. 2017;12:1–20. Chen B, Feder ME, Kang L. Evolution of heat-shock protein expression underlying adaptive responses to environmental stress. Mol Ecol. 2018;27:3040–54. Lang BJ, Guerrero ME, Prince TL, Okusha Y, Bonorino C, Calderwood SK. The functions and regulation of heat shock proteins; key orchestrators of proteostasis and the heat shock response. Arch Toxicol. 2021;95:1943–70. doi: 10.1007/s00204-021-03070-8 . Donovan MR, Marr MT. DFOXO activates large and small heat shock protein genes in response to oxidative stress to maintain proteostasis in drosophila. J Biol Chem. 2016;291:19042–50. Li Z, Srivastava P. Heat-shock proteins. Curr Protoc Immunol. 2004;Appendix 1:Appendix 1T. doi: 10.1002/0471142735.ima01ts58 . King AM, MacRae TH. Insect heat shock proteins during stress and diapause. Annu Rev Entomol. 2015;60:59–75. doi: 10.1146/annurev-ento-011613-162107 . Kwak H, Fuda NJ, Core LJ, Lis JT. Precise maps of RNA polymerase reveal how promoters direct initiation and pausing. Science. 2013;339:950–3. doi: 10.1126/science.1229386 . Dong B, Jaeger AM, Thiele DJ. Inhibiting Heat Shock Factor 1 in Cancer: A Unique Therapeutic Opportunity. Trends Pharmacol Sci. 2019;40:986–1005. Vihervaara A, Duarte FM, Lis JT. Molecular mechanisms driving transcriptional stress responses. Nat Rev Genet. 2018;19:385–97. doi: 10.1038/s41576-018-0001-6 . Jennings BH, Shah S, Yamaguchi Y, Seki M, Phillips RG, Handa H, et al. Locus-Specific Requirements for Spt5 in Transcriptional Activation and Repression in Drosophila. 2004;14:1680–4. Wu CH, Yamaguchi Y, Benjamin LR, Horvat-Gordon M, Washinsky J, Enerly E, et al. NELF and DSIF cause promoter proximal pausing on the hsp70 promoter in Drosophila. Genes Dev. 2003;17:1402–14. O’Brien T, Lis JT. RNA polymerase II pauses at the 5’ end of the transcriptionally induced Drosophila hsp70 gene. Mol Cell Biol. 1991;11:5285–90. doi: 10.1128/MCB.11.10.5285 . Yamaguchi Y, Shibata H, Handa H. Transcription elongation factors DSIF and NELF: Promoter-proximal pausing and beyond. Biochim Biophys Acta - Gene Regul Mech. 2013;1829:98–104. doi: 10.1016/j.bbagrm.2012.11.007 . Ishiguro H, Tsunoda T, Tanaka T, Fujii Y, Nakamura Y, Furukawa Y. Identification of AXUD1, a novel human gene induced by AXIN1 and its reduced expression in human carcinomas of the lung, liver, colon and kidney. Oncogene. 2001;20:5062–6. doi: 10.1038/sj.onc.1204603 . Gingras S, Pelletier S, Boyd K, Ihle JN. Characterization of a family of novel cysteine- serine-rich nuclear proteins (CSRNP). PLoS One. 2007;2:e808. doi: 10.1371/journal.pone.0000808 . Glavic A, Molnar C, Cotoras D, de Celis JF. Drosophila Axud1 is involved in the control of proliferation and displays pro-apoptotic activity. Mech Dev. 2008;126:184–97. doi: 10.1016/j.mod.2008.11.005 . Pasovic L, Eidet JR, Olstad OK, Chen DF, Lyberg T, Utheim TP. Impact of Storage Temperature on the Expression of Cell Survival Genes in Cultured ARPE-19 Cells. Curr Eye Res. 2016;3683 June:1–11. doi: 10.3109/02713683.2016.1145236 . Yi X, Bekeredjian R, DeFilippis NJ, Siddiquee Z, Fernandez E, Shohet R V. Transcriptional analysis of doxorubicin-induced cardiotoxicity. Am J Physiol Heart Circ Physiol. 2006;290:H1098-102. doi: 10.1152/ajpheart.00832.2005 . Cheng Z, Zhao H, Ze Y, Su J, Li B, Sheng L, et al. Gene-expression changes in cerium chloride-induced injury of mouse hippocampus. PLoS One. 2013;8:e60092. doi: 10.1371/journal.pone.0060092 . McDermott JE, Vartanian KB, Mitchell H, Stevens SL, Sanfilippo A, Stenzel-Poore MP. Identification and validation of ifit1 as an important innate immune bottleneck. PLoS One. 2012;7:e36465. doi: 10.1371/journal.pone.0036465 . Diercke K, Kohl A, Lux CJ, Erber R. Compression of human primary cementoblasts leads to apoptosis: A possible cause of dental root resorption? J Orofac Orthop = Fortschritte der Kieferorthopadie Organ/official J Dtsch Gesellschaft fur Kieferorthopadiedg. 2014;:1–16. doi: 10.1007/s00056-014-0237-5 . Rundqvist HC, Montelius A, Osterlund T, Norman B, Esbjornsson M, Jansson E. Acute sprint exercise transcriptome in human skeletal muscle. PLoS One. 2019;14:1–24. doi: 10.1371/journal.pone.0223024 . Simões-costa M, Stone M, Bronner ME. Axud1 integrates Wnt signaling and transcriptional inputs to drive neural crest formation. Southall TD, Gold KS, Egger B, Davidson CM, Caygill EE, Marshall OJ, et al. Cell-type-specific profiling of gene expression and chromatin binding without cell isolation: assaying RNA Pol II occupancy in neural stem cells. Dev Cell. 2013;26:101–12. doi: 10.1016/j.devcel.2013.05.020 . Marshall OJ, Brand AH. Damidseq-pipeline: An automated pipeline for processing DamID sequencing datasets. Bioinformatics. 2015;31:3371–3. Hernandez-Garcia CM, Finer JJ. Identification and validation of promoters and cis-acting regulatory elements. Plant Sci. 2014;217–218:109–19. doi: 10.1016/j.plantsci.2013.12.007 . Jo B-S, Choi SS. Introns: The Functional Benefits of Introns in Genomes. Genomics Inform. 2015;13:112–8. doi: 10.5808/GI.2015.13.4.112 . Yi X, Bekeredjian R, DeFilippis NJ, Siddiquee Z, Fernandez E, Shohet R V. Transcriptional analysis of doxorubicin-induced cardiotoxicity. Am J Physiol Heart Circ Physiol. 2006;290:H1098-102. doi: 10.1152/ajpheart.00832.2005 . Gui S, Sang X, Zheng L, Ze Y, Zhao X, Sheng L, et al. Intragastric exposure to titanium dioxide nanoparticles induced nephrotoxicity in mice, assessed by physiological and gene expression modifications. Part Fibre Toxicol. 2013;10:4. doi: 10.1186/1743-8977-10-4 . Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol Cell. 2010;38:576–89. doi: 10.1016/j.molcel.2010.05.004 . Fang D, Li Soungyu L, Xu Wanfu, KE zhiyong ZF. CSRNP. J South Med Univ. 2013;33:1122–6. Contreras O, Cruz-Soca M, Theret M, Soliman H, Tung LW, Groppa E, et al. Cross-talk between TGF-β and PDGFRα signaling pathways regulates the fate of stromal fibro-adipogenic progenitors. J Cell Sci. 2019;132. Diercke K, Kohl A, Lux CJ, Erber R. Compression of human primary cementoblasts leads to apoptosis A possible cause of dental root resorption ? Kompression führt in primären humanen Zementoblasten zur Apoptose Eine mögliche Ursache für Wurzelresorptionen. 2014;:430–45. Heinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol Cell. 2010;38:576–89. doi: 10.1016/j.molcel.2010.05.004 . Germann S, Juul-Jensen T, Letarnec B, Gaudin V. DamID, a new tool for studying plant chromatin profiling in vivo, and its use to identify putative LHP1 target loci. Plant J. 2006;48:153–63. Southall TD, Gold KS, Egger B, Davidson CM, Caygill EE, Marshall OJ, et al. Cell-type-specific profiling of gene expression and chromatin binding without cell isolation: assaying RNA Pol II occupancy in neural stem cells. Dev Cell. 2013;26:101–12. doi: 10.1016/j.devcel.2013.05.020 . Aughey GN, Southall TD. Dam it’s good! DamID profiling of protein-DNA interactions. Wiley Interdiscip Rev Dev Biol. 2016;5:25–37. Johnson DG, Walker CL. Cyclins and cell cycle checkpoints. Annu Rev Pharmacol Toxicol. 1999;39:295–312. doi: 10.1146/annurev.pharmtox.39.1.295 . Bunch H. RNA polymerase II pausing and transcriptional regulation of the HSP70 expression. Eur J Cell Biol. 2017;96:739–45. doi: 10.1016/j.ejcb.2017.09.003 . Korb K, Katsikogianni E, Zingler S, Daum E, Lux CJ, Hohenstein A, et al. Inhibition of AXUD1 attenuates compression-dependent apoptosis of cementoblasts. Clin Oral Investig. 2016;20:2333–41. doi: 10.1007/s00784-016-1740-4 . Zhao Y, Sun H, Lu J, Li X, Chen X, Tao D, et al. Lifespan extension and elevated hsp gene expression in Drosophila caused by histone deacetylase inhibitors. J Exp Biol. 2005;208:697–705. Gingras S, Pelletier S, Boyd K, Ihle JN. Characterization of a family of novel cysteine- serine-rich nuclear proteins (CSRNP). PLoS One. 2007;2:e808. doi: 10.1371/journal.pone.0000808 . Missra A, Gilmour DS. Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation complex. Proc Natl Acad Sci U S A. 2010;107:11301–6. doi: 10.1073/pnas.1000681107 . Zhimulev IF, Belyaeva ES, Semeshin VF, Koryakov DE, Demakov S a, Demakova O V, et al. Polytene chromosomes: 70 years of genetic research. Int Rev Cytol. 2004;241:203–75. doi: 10.1016/S0074-7696(04)41004-3 . Courgeon A-M, Rollet E, BEcket J, Maisonhaute C, Best-Belpomme M. Hydrogen peroxide (H2O2) induces actin and some heat-shock proteins in Drosophila cells. Eur J Biochem. 1988;171:163–70. doi: 10.1111/j.1432-1033.1988.tb13772.x . Chatr-Aryamontri A, Breitkreutz BJ, Oughtred R, Boucher L, Heinicke S, Chen D, et al. The BioGRID interaction database: 2015 update. Nucleic Acids Res. 2015;43:D470–8. Saunders A, Werner J, Andrulis ED, Nakayama T, Hirose S, Reinberg D, et al. Tracking FACT and the RNA polymerase II elongation complex through chromatin in vivo. Science. 2003;301:1094–6. Gilchrist DA, Nechaev S, Lee C, Ghosh SKB, Collins JB, Li L, et al. NELF-mediated stalling of Pol II can enhance gene expression by blocking promoter-proximal nucleosome assembly. Genes Dev. 2008;22:1921–33. Ghosh SKB, Missra A, Gilmour DS. Negative Elongation Factor Accelerates the Rate at Which Heat Shock Genes Are Shut off by Facilitating Dissociation of Heat Shock Factor. Mol Cell Biol. 2011;31:4232–43. Sørensen JG, Loeschcke V. Larval crowding in Drosophila melanogaster induces Hsp70 expression, and leads to increased adult longevity and adult thermal stress resistance. J Insect Physiol. 2001;47:1301–7. Bosch M, Serras F, Martín-Blanco E, Baguñà J. JNK signaling pathway required for wound healing in regenerating Drosophila wing imaginal discs. Dev Biol. 2005;280:73–86. doi: 10.1016/j.ydbio.2005.01.002 . Zhu S, Chen R, Soba P, Jan YN. JNK signaling coordinates with ecdysone signaling to promote pruning of Drosophila sensory neuron dendrites. Dev. 2019;146:1–10. Zhang H, Qiu X, Yang G. The CSRNP Gene Family Serves as a Prognostic Biomarker in Clear Cell Renal Cell Carcinoma. 2021;11 March:1–12. Sherman MY, Gabai VL. Hsp70 in cancer: back to the future. 2015; October 2014:4153–61. Jagadish N, Agarwal S, Gupta N, Fatima R, Devi S, Kumar V, et al. Heat shock protein 70-2 (HSP70-2) overexpression in breast cancer. J Exp Clin Cancer Res. 2016;2:1–14. McDuffee AT, Senisterra G, Huntley S, Lepock JR, Sekhar KR, Meredith MJ, et al. Proteins containing non-native disulfide fonds generated by oxidative stress can act as signals for the induction of the heat shock response. J Cell Physiol. 1997;171:143–51. Mosaddegh B, Takalloo Z, Sajedi RH, Shirin Shahangian S, Hassani L, Rasti B. An inter-subunit disulfide bond of artemin acts as a redox switch for its chaperone-like activity. Cell Stress Chaperones. 2018;23:685–93. Wu CH, Lee C, Fan R, Smith MJ, Yamaguchi Y, Handa H, et al. Molecular characterization of Drosophila NELF. Nucleic Acids Res. 2005;33:1269–79. Sims RJ, Mandal SS, Reinberg D. Recent highlights of RNA-polymerase-II-mediated transcription. Curr Opin Cell Biol. 2004;16:263–71. Lis JT, Mason P, Peng J, Price DH, Werner J. P-TEFb kinase recruitment and function at heat shock loci P-TEFb kinase recruitment and function at heat shock loci. 2000;:792–803. Gilchrist DA, Fromm G, dos Santos G, Pham LN, Mcdaniel IE, Burkholder A, et al. Regulating the regulators: The pervasive effects of Pol II pausing on stimulus-responsive gene networks. Genes Dev. 2012;26:933–44. Missra A, Gilmour DS. Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation complex. Proc Natl Acad Sci U S A. 2010;107:11301–6. doi: 10.1073/pnas.1000681107 . Missra A, Gilmour DS. Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation complex. Proc Natl Acad Sci U S A. 2010;107:11301–6. doi: 10.1073/pnas.1000681107 . Macdonald CD, Falconer AMD, Chan CM, Wilkinson DJ, Skelton A, Reynard L, et al. Cytokine-induced cysteine- serine-rich nuclear protein-1 (CSRNP1) selectively contributes to MMP1 expression in human chondrocytes. PLoS One. 2018;13:1–16. Yang Q, Liu X, Zhou T, Cook J, Nguyen K, Bai X. RNA polymerase II pausing modulates hematopoietic stem cell emergence in zebrafish. Blood. 2016;128:1701–10. Solı C, Espina J, Feijo CG. csrnp1a Is Necessary for the Development of Primitive Hematopoiesis Progenitors in Zebrafish. 2013;8. Kroeger PE, Rowe TC. Analysis of Topoisomerase I and II Cleavage Sites on the Drosophila Actin and Hsp70 Heat Shock Genes. Biochemistry. 1992;31:2492–501. Downes CS, Clarke DJ, Mullinger AM, Giménez-Abián JF, Creighton AM, Johnson RT. A topoisomerase II-dependent G2 cycle checkpoint in mammalian cells. Nature. 1994;372:467–70. Gingras S, Pelletier S, Boyd K, Ihle JN. Characterization of a family of novel cysteine- serine-rich nuclear proteins (CSRNP). PLoS One. 2007;2:e808. doi: 10.1371/journal.pone.0000808 . Yamada K, Akiyama N, Yamada S, Tanaka H, Saito S, Hiraoka M, et al. Taip2 is a novel cell death-related gene expressed in the brain during development. Biochem Biophys Res Commun. 2008;369:426–31. doi: 10.1016/j.bbrc.2008.02.041 . Feijóo CG, Sarrazin AF, Allende ML, Glavic A. Cystein-serine-rich nuclear protein 1, Axud1/Csrnp1, is essential for cephalic neural progenitor proliferation and survival in zebrafish. Dev Dyn. 2009;238:2034–43. doi: 10.1002/dvdy.22006 . Feijóo CG, Sarrazin AF, Allende ML, Glavic A. Cystein-serine-rich nuclear protein 1, Axud1/Csrnp1, is essential for cephalic neural progenitor proliferation and survival in zebrafish. Dev Dyn. 2009;238:2034–43. doi: 10.1002/dvdy.22006 . Azambuja AP, Simoes-Costa M. A regulatory sub-circuit downstream of Wnt signaling controls developmental transitions in neural crest formation. PLoS Genet. 2021;17:1–23. doi: 10.1371/journal.pgen.1009296 . Johansen KM, Cai W, Deng H, Bao X, Zhang W, Girton J, et al. Polytene chromosome squash methods for studying transcription and epigenetic chromatin modification in Drosophila using antibodies. Methods. 2009;48:387–97. doi: 10.1016/j.ymeth.2009.02.019 . Marshall OJ, Brand AH. Damidseq-pipeline: An automated pipeline for processing DamID sequencing datasets. Bioinformatics. 2015;31:3371–3. Bailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME Suite: Tools for motif discovery and searching. Nucleic Acids Res. 2009;37 SUPPL. 2:202–8. Ramírez F, Dündar F, Diehl S, Grüning BA, Manke T. DeepTools: A flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42:187–91. Supplementary Files Additional.TableS1GeneAnnotationAllPeaks.xltx Additional.TableS2GeneOntologyAll.xlsx Additional.TableS3AGAGTGannotation.xlsx Additional.TableS4AGAGTGPeaksGenes.xlsx Additional.TableS5GeneOntologyAGAGTG.xlsx AdditionalFiguresBMCsubmission.pdf Cite Share Download PDF Status: Posted Version 1 posted You are reading this latest preprint version Research Square lets you share your work early, gain feedback from the community, and start making changes to your manuscript prior to peer review in a journal. As a division of Research Square Company, we’re committed to making research communication faster, fairer, and more useful. We do this by developing innovative software and high quality services for the global research community. Our growing team is made up of researchers and industry professionals working together to solve the most critical problems facing scientific publishing. Also discoverable on Platform About Our Team In Review Editorial Policies Advisory Board Help Center Resources Author Services Accessibility API Access RSS feed Manage Cookie Preferences © Research Square 2026 | ISSN 2693-5015 (online) Privacy Policy Terms of Service Do Not Sell My Personal Information {"props":{"pageProps":{"initialData":{"identity":"rs-1040684","acceptedTermsAndConditions":true,"allowDirectSubmit":true,"archivedVersions":[],"articleType":"Research","associatedPublications":[],"authors":[{"id":61780229,"identity":"c2980e58-e1bd-4e8a-884b-43e988578340","order_by":0,"name":"JM Zuñiga-Hernández","email":"data:image/png;base64,iVBORw0KGgoAAAANSUhEUgAAAZAAAAAyAQMAAABI0h/eAAAABlBMVEX///8AAABVwtN+AAAACXBIWXMAAA7EAAAOxAGVKw4bAAAA9klEQVRIiWNgGAWjYBACxgYwJcHAwA5hJbCxg6gCYrTwHIBqYQZRBsTYJ5EA0cJASAtze3fyZ94dFnb9M98Yv2Bss8vjY2Y+wPADjxbGnrMbjHnPSCTPuJ1jZsHYllzMxsyWwNiDT8uM3A3JvG0SyQxALQaMbQcS25h5DJjxOQyk5TBIi/zNMzAt/B8IadnYDNRiZ3CDx/gB1BYG/Fp6zm5mnNsmkWB4Jq2MIeFcMlALm8FBfH4xbO/d/OFtW5293PHDmz98KLNLnN/e/PDBjwo8WhogdCKQZoNGDQPDAdwaGBjkobQ9EDN/wKdyFIyCUTAKRi4AAEG/TMcongaGAAAAAElFTkSuQmCC","orcid":"https://orcid.org/0000-0001-5821-4384","institution":"University of Chile Faculty of Sciences: Universidad de Chile Facultad de Ciencias","correspondingAuthor":true,"submittingAuthor":false,"prefix":"","firstName":"JM","middleName":"","lastName":"Zuñiga-Hernández","suffix":""},{"id":61780230,"identity":"9e55701c-4943-42dd-9e45-b5b030844d40","order_by":1,"name":"C Meneses","email":"","orcid":"","institution":"Universidad Andrés Bello: Universidad Andres Bello","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"C","middleName":"","lastName":"Meneses","suffix":""},{"id":61780231,"identity":"f81b8462-93cf-4b60-9c7c-83113a2256b2","order_by":2,"name":"M Bastías","email":"","orcid":"","institution":"Universidad Andrés Bello: Universidad Andres Bello","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"M","middleName":"","lastName":"Bastías","suffix":""},{"id":61780232,"identity":"69c9c739-ed6f-4889-85f4-905bc2738c8d","order_by":3,"name":"ML Allende","email":"","orcid":"","institution":"Universidad de Chile Facultad de Ciencias","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"ML","middleName":"","lastName":"Allende","suffix":""},{"id":61780233,"identity":"aacc5ce9-aa25-4883-bf08-4556d8a70aa8","order_by":4,"name":"Alvaro Glavic","email":"","orcid":"","institution":"Universidad de Chile Facultad de Ciencias","correspondingAuthor":false,"submittingAuthor":false,"prefix":"","firstName":"Alvaro","middleName":"","lastName":"Glavic","suffix":""}],"badges":[],"createdAt":"2021-11-01 17:49:25","currentVersionCode":1,"declarations":"","doi":"10.21203/rs.3.rs-1040684/v1","doiUrl":"https://doi.org/10.21203/rs.3.rs-1040684/v1","draftVersion":[],"editorialEvents":[],"editorialNote":"","failedWorkflow":false,"files":[{"id":15441290,"identity":"3cced9fd-9fc1-49dc-9832-fbc76fc9be04","added_by":"auto","created_at":"2021-11-11 15:22:45","extension":"png","order_by":1,"title":"Figure 1","display":"","copyAsset":false,"role":"figure","size":97816,"visible":true,"origin":"","legend":"TaDa-seq from Dam-DAxud1/Dam shows the genome-wide distribution and main motif in significative binding regions. 1811 peaks were identified as significative over the background from the subtraction of signals (Dam-DAxud1/Dam (no fusion) using the Dam/TaDa pipeline and annotated with HOMER. (A) Annotation of significative peaks reveals that Dam-DAxud1 is prevalent in the first intron, intergenic, and promoter zones. (B) The metagene profile and heatmap show the main signals on the TSS surrounding zone with a low signal at the gene body and TTS (transcription termination sites). Gene bodies of all genes were scaled to 1000 bp. (C) Gene tracks for hsp genes confirm the presence of Dam-DAxud1 in these sites, concordantly with the genome distribution in (A). Positive average scores are in blue and negative are in orange. (D) The main motif identified with HOMER using the denovo application, matching with the Topoisomerase 2 binding site motif, whose search using GOMO (Gene Ontology of Motif using promoter zones) matches the biological process related to heat stress and hypoxia. ","description":"","filename":"Onlinefloatimage1.png","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/275259024840943743696f69.png"},{"id":15441660,"identity":"60d259c6-46c7-47a8-8a8d-06f37594f815","added_by":"auto","created_at":"2021-11-11 15:25:44","extension":"png","order_by":2,"title":"Figure 2","display":"","copyAsset":false,"role":"figure","size":41062,"visible":true,"origin":"","legend":"Effect of dAxud1 knockdown in pupation and lifespan. (A) and (B) represent the survival chart from dAxud1 knockdown (yw; UAS-IRdAxud1/+; Tub\u003eGal4/+) genotypes under heat shock episodes of 37° for 30 minutes every day (DHS) and compared with flies in stable temperature at 29°C (NHS). p\u003c 0.0001 with Log-rank statistical analysis. Every assay was performed with 25-35 flies with 3 replicates. Figs. C, D, and E represent the gene expression of two different tissues for three heat shock genes, hsp70, hsp26, and hsp67, respectively. The expression was tested on 5 different conditions and compared to control conditions, including dAxud knockdown (nubbin Gal4/+; UAS-IRdAxud1) and dAxud overexpression (nubbin Gal4/+; UAS-dAxud1-GFP). The control genotype was nubbin Gal4/+ (Driver only).","description":"","filename":"Onlinefloatimage2.png","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/c9d8fef93385affe8cff20bb.png"},{"id":15440199,"identity":"e4569ddf-7c21-4086-be00-1938a5411dc5","added_by":"auto","created_at":"2021-11-11 15:19:44","extension":"png","order_by":3,"title":"Figure 3","display":"","copyAsset":false,"role":"figure","size":102484,"visible":true,"origin":"","legend":"DAxud1-GFP has an extensive distribution on polytene chromosomes. In accordance with Johanssen’s protocol [30], immunostaining was performed for RNA polymerase IIo and GFP for DAxud1-GFP protein localization; DAPI was added to visualize condensed DNA. A) 20X picture from 17°C growth third instar larvae. Bar 20 µm. B) Scheme with hsp70 genes cytolocation at 87A-87C. C) Changes in hsp70 locus during heat shock (37°C at 20 minutes) reveal the relocation of DAxud1-GFP in the same way as RNA Pol IIo. Bar: 5 µm. D) CHIP analysis from third instar larvae salivary glands for DAxud1-GFP on hsp70 locus for 5’ and 3’ gene zones. Control conditions were defined as 25°C. The oxidative stress condition was induced by feeding larvae H2O2 0.5% v/v for 24 hours (late 2nd to 3rd instar larvae, then dissected). Heat shock was performed by exposing 3rd instar larvae to 37°C for 20 minutes; then the salivary glands were dissected for chromatin extraction and immunoprecipitation.","description":"","filename":"Onlinefloatimage3.png","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/94ad0d3755c614db8c8e0693.png"},{"id":15441287,"identity":"20cb4318-6995-4553-8c23-2703c91d88a4","added_by":"auto","created_at":"2021-11-11 15:22:44","extension":"png","order_by":4,"title":"Figure 4","display":"","copyAsset":false,"role":"figure","size":48714,"visible":true,"origin":"","legend":"CHIP analysis of RNA Polymerase IIo occupancy on hsp70 locus under DAxud1 overexpression and DAxud1 knockdown. A) CHIP analysis RNA Pol II (RPB1 subunit) occupancy on hsp70A locus. Each condition was examined in triplicate B) Comparison of the 5’/3’ occupancy ratio between conditions. Statistical analysis was performed as a Tukey test, * p \u003c 0.05. D) Western blot from the co-immunoprecipitation of DAxud1-GFP with NELF-B-HA. D) qPCR for hsp70 from salivary glands (3rd instar) expressing drive nubbin\u003eGal4 and Dicer2 (to improve RNAi silencing activity). The hsp70 levels were measured at 20’ heat shock exposure (37°C) and heat shock recovery (40’ at 25°C after 20’ heat shock) to see if there were differences between levels of hsp70 at the same recovery time. Statistical analysis was performed as a t-test, * p \u003c 0.05. ","description":"","filename":"Onlinefloatimage4.png","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/0fdc7f48b23afdfc6750608e.png"},{"id":16504420,"identity":"42db7d40-f949-4ae6-a552-602568038613","added_by":"auto","created_at":"2021-12-16 06:56:37","extension":"pdf","order_by":0,"title":"","display":"","copyAsset":false,"role":"manuscript-pdf","size":1801982,"visible":true,"origin":"","legend":"","description":"","filename":"manuscript.pdf","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/f79c8b11-4969-4bbf-9ccd-a709b1be6d6c.pdf"},{"id":15440198,"identity":"2c6b283f-6536-48df-a1a0-e969a27921f3","added_by":"auto","created_at":"2021-11-11 15:19:44","extension":"xltx","order_by":1,"title":"","display":"","copyAsset":false,"role":"supplement","size":62095,"visible":true,"origin":"","legend":"","description":"","filename":"Additional.TableS1GeneAnnotationAllPeaks.xltx","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/ef8f41b26e6713757bb52104.xltx"},{"id":15440206,"identity":"d5783a06-9524-43ab-8f13-0fda53f4bb39","added_by":"auto","created_at":"2021-11-11 15:19:45","extension":"xlsx","order_by":2,"title":"","display":"","copyAsset":false,"role":"supplement","size":29795,"visible":true,"origin":"","legend":"","description":"","filename":"Additional.TableS2GeneOntologyAll.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/0bda7e8af16e43821cc4e867.xlsx"},{"id":15440205,"identity":"805fca9e-0dc7-44e4-95dd-c9ff3b20bfd6","added_by":"auto","created_at":"2021-11-11 15:19:45","extension":"xlsx","order_by":3,"title":"","display":"","copyAsset":false,"role":"supplement","size":76296,"visible":true,"origin":"","legend":"","description":"","filename":"Additional.TableS3AGAGTGannotation.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/2518b69c9cb8892b09290e8a.xlsx"},{"id":15440203,"identity":"877ca579-fc6f-4ce3-8627-6ee1b7de6767","added_by":"auto","created_at":"2021-11-11 15:19:45","extension":"xlsx","order_by":4,"title":"","display":"","copyAsset":false,"role":"supplement","size":718248,"visible":true,"origin":"","legend":"","description":"","filename":"Additional.TableS4AGAGTGPeaksGenes.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/d4c764dd9e78cddc490e0e36.xlsx"},{"id":15440202,"identity":"5e557f48-9356-4afb-839b-f9170dd78648","added_by":"auto","created_at":"2021-11-11 15:19:45","extension":"xlsx","order_by":5,"title":"","display":"","copyAsset":false,"role":"supplement","size":14408,"visible":true,"origin":"","legend":"","description":"","filename":"Additional.TableS5GeneOntologyAGAGTG.xlsx","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/37994ee05d83fa77290dd909.xlsx"},{"id":15441289,"identity":"dbdc57c9-9068-4ad0-b894-be4fbac22736","added_by":"auto","created_at":"2021-11-11 15:22:45","extension":"pdf","order_by":6,"title":"","display":"","copyAsset":false,"role":"supplement","size":890159,"visible":true,"origin":"","legend":"","description":"","filename":"AdditionalFiguresBMCsubmission.pdf","url":"https://assets-eu.researchsquare.com/files/rs-1040684/v1/48d56e923777ba9839928687.pdf"}],"financialInterests":"","formattedTitle":"\u003cp\u003eDrosophila DAxud1: A New Element in Transcriptional Pausing Complex Stabilization\u003c/p\u003e","fulltext":[{"header":"Introduction","content":"\u003cp\u003eIn natural environments, organisms are exposed to different stresses that affect their survival. During evolution, different mechanisms have emerged to respond to stress, some of them widely conserved among various kingdoms, but with differences in their temporal regulation and cell types involved. One of the conserved mechanisms for responding to stress is rapid biosynthesis of Heat Shock Proteins (HSPs), in which their mRNA levels can rise several fold in a few minutes, mainly due to transcriptional promoting elements preloaded at the core promoter and within the 5\u0026rsquo; end of the gene body\u0026nbsp;[1, 2]. HSPs were initially associated with heat stress responses. Further research revealed that these proteins also behave as important cytoprotectors under other stress conditions, including oxidative stress, accumulation of protein aggregates, an acidic environment, and osmolarity changes\u0026nbsp;[3]. In animals, \u003cem\u003ehsp\u003c/em\u003e transcriptional induction is very fast; the \u003cem\u003ehsp70\u003c/em\u003e promoter recruits specific factors in response to stress conditions\u0026nbsp;[4], resulting in significative rise of \u003cem\u003ehsp70\u003c/em\u003e mRNA in less than 5 minutes\u0026nbsp;[5]. This recruitment during heat shock is commanded mainly by HSF (Heat Shock Factor), a transcriptional activator, one of five known chromatin proteins involved in the heat stress transcriptional response, that induces Hsp70 and other small HSPs (Hsp27, Hsp26) after acute thermal stress. However, the repertoire of chromatin associated proteins varies for other stress conditions and the differences are still poorly understood\u0026nbsp;[6\u0026ndash;8]. We will focus mainly on the \u003cem\u003ehsp70\u003c/em\u003e paralogs because there is significant knowledge implicating the chromatin state of their promoters and the pausing phenomenon\u0026nbsp;[9\u0026ndash;14], as well as additional factors regulating their efficient transcriptional induction.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;In \u003cem\u003eDrosophila melanogaster\u003c/em\u003e, the control of \u003cem\u003ehsp70\u003c/em\u003e transcription depends on at least 12 proteins that relocate to the \u003cem\u003ehsp70\u003c/em\u003e \u003cem\u003eloci\u0026nbsp;\u003c/em\u003eafter two minutes of exposure to 37\u0026deg;C in 3\u003csup\u003erd\u003c/sup\u003e instar larvae\u0026nbsp;[4, 13\u0026ndash;16]. These proteins participate regulating chromatin condensation, recruiting transcription factors like HSF, engaging in transcriptional pausing, and controlling the initiation/elongation of \u003cem\u003ehsp70\u003c/em\u003e transcription. In \u003cem\u003ehsp70\u003c/em\u003e genes, transcriptional pausing is crucial for a quick transcriptional response upon stress or extracellular signaling. When pausing, RNA Polymerase II elongates only 30-50 bp from the transcription start site (TSS) and remains in that position until induced, where it resumes elongation and provides a transcriptional burst of the stress response genes\u0026nbsp;[14, 17]. Transcriptional pausing is mediated by the NELF (NELF A, B, C/D, E) and DSIF (Spt5, Spt6) complexes, which retain RNA Polymerase II. After induction by heat shock, HSF recruits the p-TEFb kinase that phosphorylates NELF, DSIF and RNA Polymerase II, causing complete dissociation of NELF and the and the escape of RNA Polymerase, with DSIF acting as an elongation factor\u0026nbsp;[4, 18].\u0026nbsp;In this work, we present evidence that introduces DAxud1 as an additional regulatory factor that participates in \u003cem\u003ehsp\u003c/em\u003e transcription regulation, through its positioning at the transcription start site (TSS) and by stabilizing the pausing complex via its interaction with NELF and DSIF complexes.\u003c/p\u003e\n\u003cp\u003e\u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp; \u0026nbsp;\u0026nbsp;The \u003cem\u003edAxud1\u003c/em\u003e gene is the only \u003cem\u003eDrosophila\u003c/em\u003e ortholog of the human and mice CSRNP protein family (Cysteine Serine Rich Nuclear Proteins). Members of this protein family, which are conserved only in metazoans, have been proposed to be putative transcription factors with specific roles in apoptosis and neural development\u0026nbsp;[19\u0026ndash;21]. Initial studies proposed a physiological function for one of these proteins as a tumor suppressor in human tissues\u0026nbsp;[19]. Further characterization of its ortholog in \u003cem\u003eDrosophila\u003c/em\u003e suggested that its proapoptotic function relies on the activation of the JNK pathway\u0026nbsp;[21]. Other studies indicate a strong association of the CSRNP family with stress responses as their mRNAs rise in cells exposed to different kinds of stressful stimuli such as cold\u0026nbsp;[22], oxidative stress\u0026nbsp;[23, 24], bacterial infection\u0026nbsp;[25], pressure\u0026nbsp;[26], and acute stress by sprint running in skeletal muscle\u0026nbsp;[27].\u0026nbsp;These observations prompted us to analyze the molecular bases of the adaptive and non-adaptive roles of DAxud1 during the stress response. Here, in a screen using whole-genome occupancy analysis (TaDa-seq), we reveal how DAxud1 has a connection to the general stress response by analyzing its occupancy of the \u003cem\u003ehsp70 locus\u003c/em\u003e. In this gene, DAxud1 associates with the Pol II pausing complex through its interaction with the NELF complex and thereby directly influences \u003cem\u003ehsp\u003c/em\u003e transcription during the stress response.\u0026nbsp;\u003c/p\u003e"},{"header":"Methods","content":"\u003cp\u003e\u003cstrong\u003eFly stocks\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eLoss of function was achieved by expressing \u003cem\u003edAxud1\u003c/em\u003e RNAi through the GAL4/UAS system using the Vienna stock V26479. For the experiment for polytene chromosomes, we used flies with the genotype \u003cem\u003enub\u003c/em\u003e\u0026gt;\u003cem\u003eGal4/Y\u003c/em\u003e; \u003cem\u003eUAS-dAxud1::GFP/+.\u003c/em\u003e For DamID-seq experiments (TaDa-seq), the genotypes were \u003cem\u003enub\u003c/em\u003e\u0026gt;Gal4/Y; \u003cem\u003eUAS-mCherry.Dam\u0026nbsp;\u003c/em\u003eand \u003cem\u003enub\u003c/em\u003e\u0026gt;\u003cem\u003eGal4\u003c/em\u003e; \u003cem\u003eUAS-mCherry.Dam::dAxud1\u003c/em\u003e, both generated from the plasmid pUAST-attb-LT3-Dam donated by the Andrea Brand Laboratory\u0026nbsp;[41]. The Dam fusion plasmids were injected with pBS130 that express the \u003cem\u003ephi\u003c/em\u003e-131 integrase. The genotype of injected flies was \u003cem\u003eP{y[+t7.7]=CaryP}attP2\u003c/em\u003e (Bloomington stock 8622), in accordance with methods in the previous reference [18]. For HA-NELF-B experiments in imaginal wing discs protein extract, stock F003904 from FlyORF was used for expression under the UAS/Gal4 system.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eLifespan assay\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFor every replicate, between 80-100 flies were collected no longer than 48 hours from pupae eclosion and separated into males and females. Three replicates for every condition were used, and dead and live flies were counted every day. The control condition was at 29\u0026deg;C, with the flies changed from the tubes every two days. The heat shock condition was at 37\u0026deg;C for 30 minutes every day.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eImmunofluorescence in polytene chromosomes\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eThe squash protocol and antibody staining for polytene chromosomes were performed according to the Johanssen protocol\u0026nbsp;[79]. The GFP antibody used was Cell Signaling (code D5.1), dilution 1:250, and the RNA Polymerase IIo antibody (full phosphorylated CTD) used was Abcam (ab5408), dilution 1:1000. DAPI stain was used at 1:10000. The glands were dissected from larvae with \u003cem\u003enub\u0026gt;Gal4; UAS-dAxud1::GFP\u003c/em\u003e genotype.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eTaDa-seq\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn total, 100 to 120 imaginal wing discs per sample were dissected from third instar larvae (\u003cem\u003enub\u0026gt;Gal4\u003c/em\u003e; \u003cem\u003eUAS-mCherry.Dam\u003c/em\u003e \u0026amp; \u003cem\u003enub\u0026gt;Gal4\u003c/em\u003e; \u003cem\u003eUAS-mCherry.Dam::dAxud1).\u0026nbsp;\u003c/em\u003eThe larvae grew for 7 days at 17\u0026deg;C and were then placed in 29\u0026deg;C for 24 hours, to the third instar stage. The discs were dissected and placed in PBS 1X, and the DNA was extracted according to the Southall protocol\u0026nbsp;[41]. Accordingly, two biological replicates for every condition (Dam and Dam-DAxud1) were used to generate libraries, sequenced in Illumina Hiseq-2500, single end 100 bp. The data were processed with the pipeline designed and published by Marshall\u0026nbsp;[80]\u0026nbsp;using the dm6 \u003cem\u003eDrosophila\u003c/em\u003e genome version (release 6.22). The replicates were merged and data was generated in bedgraph format, according with the pipeline . The GFF files with the coordinates of significant peaks were visualized with the WashU epigenome browser. These data were processed with HOMER\u0026nbsp;[35]\u0026nbsp;for motif discovery and finding. The annotation was done with T-Gene from Meme-Suite\u0026nbsp;[81]. The metagene and heatmap profiles were generated with the deeptools toolbox platform\u0026nbsp;[82].\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eRNA extraction and qPCR\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eRNA extraction was carried out with the TriZol reagent following manufacturer instructions. The RNA was resuspended in water and the RT reaction was performed with the iScript\u0026reg; kit from BioRad. The qPCR mix reactions were done with Brilliant II SYBR\u0026reg; reaction solution with 10 minutes of initial denaturation and 40 cycles in the following sequence: 30 seconds at 90\u0026deg;C; 30 seconds at 60\u0026deg;C; and 30 seconds at 72\u0026deg;C. The primers used were as follows: \u003cem\u003edAxud1\u003c/em\u003e Fw 5\u0026rsquo;- AGGGGACCACCAGCCTAAC-3\u0026rsquo;, Rv 5\u0026rsquo;-GGTTCGCTCTGATTATCCTTGTG-3\u0026rsquo;; \u003cem\u003ehsp26\u0026nbsp;\u003c/em\u003eFw 5\u0026rsquo;-ATGCCCACGATCTGTTCCATC-3\u0026rsquo;, Rv 5\u0026rsquo;-GTACGCGAATAACGACGAC; \u003cem\u003ehsp67Bc\u0026nbsp;\u003c/em\u003eFw 5\u0026rsquo;-GACTCCCCGGACTCCATGTA, Rv 5\u0026rsquo;-GACCCAGGGTGTGCAAATCAA-3\u0026rsquo;; \u003cem\u003ehsp70\u003c/em\u003e for B paralogs (for Ba, Bb and Bc) Fw 5\u0026rsquo;- AAGAACCTCAAGGGTGAGCG-3\u0026rsquo;, Rv 5\u0026rsquo;-CGAACAGAGATCCCTCGTCG-3\u0026rsquo;. For \u003cem\u003ehsp70B\u0026nbsp;\u003c/em\u003e5\u0026rsquo;-UTR Fw 5\u0026rsquo;- GCTAAGCAAATAAACAAGCGCAG-3\u0026rsquo;, Rv 5\u0026rsquo;-CAGTTGATTTACTTGGTTGCTGGT-3\u0026rsquo; ; \u003cem\u003ehsp70B\u0026nbsp;\u003c/em\u003e3\u0026rsquo;-UTR Fw 5\u0026rsquo;-GAGGATTTGGCGGCTACTCT-3\u0026rsquo;, Rv 5\u0026rsquo;-TTTAAAAACTTAAGCCAGGAACTGA-3\u0026rsquo;; \u003cem\u003eactin-42A\u003c/em\u003e (as normalizer) Fw 5\u0026rsquo;-GCGTCGGTCAATTCAATCTT-3\u0026rsquo;, Rv 5\u0026rsquo;-AAGCTGCAACCTCTTCGTCA-3\u0026rsquo;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eChip-PCR protocol\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eChromatin immunoprecipitation was performed from salivary glands expressing DAxud1-GFP or \u003cem\u003edAxud1\u003c/em\u003e RNAi in third instar larvae growth at 25\u0026deg;C using the driver \u003cem\u003enubbin\u003c/em\u003e\u0026gt;Gal4 (X). Larvae from the heat shock condition were placed at 37\u0026deg;C for twenty minutes. The glands were dissected and processed according to Ghosh et al.\u0026rsquo;s protocol\u0026nbsp;[54]\u0026nbsp;for three biological replicates per condition.\u0026nbsp;Ten pairs of salivary glands from third-instar larvae were incubated for five minutes on ice in 100 \u0026mu;l of 1% formaldehyde in PBS 1X and then at RT for 7 minutes. The cross-linking reaction was quenched by adding 2.5 M glycine to a final concentration of 125 mM; then the glands were placed on ice for two minutes. They were centrifuged at 900 \u0026times; \u003cem\u003eg\u003c/em\u003e for two minutes at 4\u0026deg;C, and the supernatant was removed. A 100-\u0026mu;l volume of sonication buffer (20 mM Tris [pH 8.0], 0.5% SDS, 2 mM EDTA, 0.5 mM EGTA, 0.5 mM phenylmethylsulphonyl fluoride [PMSF]) and 1 \u0026mu;l Halt protease inhibitor cocktail (Thermo Fisher cat. 78430) were added to the glands to incubate at room temperature for ten minutes and then on ice for 10 minutes. The glands were vigorously shaken for ten minutes and homogenized with a small pestle. Lysates were sonicated at 4\u0026deg;C in Omniruptor-4000 at 100% power, 90% pulse for fifteen minutes to shear the DNA to an average fragment size of 400 bp. The lysate was clarified by centrifugation at 14,000 \u0026times; \u003cem\u003eg\u003c/em\u003e for 7 minutes. In total, 45 \u0026mu;l per lysate were used for each immunoprecipitation assay (mock and IP), with 10 \u0026mu;l for input. Immunoprecipitation was performed using 1 \u0026mu;l of anti-GFP antibody (Abcam ab290, 5\u0026micro;g for chromatin from teen pairs of salivary glands) or 30 \u0026mu;l of RNA Polymerase II antibody (Santa Cruz, 8WG16, \u0026nbsp;6\u0026micro;g for chromatin from ten pairs of salivary glands), diluted half in IP Buffer (50mM Tris-HCl pH 8; 100mM NaCl; 2mM EDTA; 1mM EDTA, 1% NP40, 1X Thermo Fisher protease cocktail), and immunoprecipitated with 30 \u0026mu;l of protein-A Dynabeads\u0026reg; solution, with the same quantity for the mock assay. Beads were washed twice in sequence with CHIP 1 Buffer (IP buffer + 1% Na-deoxycholate), CHIP 2 Buffer (IP Elution buffer + 1% Na-deoxycholate + 500mM NaCl), and CHIP 3 Buffer (IP Elution buffer + 1% Na-deoxycholate + 270mM LiCl), then washed twice in Tris-HCl 10 mM pH 8. Finally, elution was performed with elution buffer (NaHCO\u003csub\u003e3\u003c/sub\u003e pH 8.8 0.1M; SDS 1%), added to the precipitated samples until 200 \u0026mu;l was reached (same as input).\u003c/p\u003e\n\u003cp\u003eDecrosslinking was performed by adding 10 \u0026mu;l of NaCl 4M and 0.5 \u0026mu;l of fungal proteinase K (Thermo Fisher), with six hours of incubation at 65\u0026deg;C. The samples were purified and eluted with MicroChip Diapure Columns. Samples were analyzed with qPCR.\u0026nbsp;The primers used were as follows: \u003cem\u003ehsp70Aa\u0026nbsp;\u003c/em\u003eTSS\u003cem\u003e\u0026nbsp;\u003c/em\u003e(-800,-680)\u003cem\u003e\u0026nbsp;\u003c/em\u003eFw 5\u0026rsquo;-AACAACAAATTCCAAGTTTGCAC-3\u0026rsquo;, Rv 5\u0026rsquo;-CACGGTTTTTGCATATGCTTT-3\u0026rsquo;; \u003cem\u003ehsp70B\u003c/em\u003e paralogs\u003cem\u003e\u0026nbsp;\u003c/em\u003e5\u0026rsquo;-UTR Fw 5\u0026rsquo;- GCTAAGCAAATAAACAAGCGCAG-3\u0026rsquo;, Rv 5\u0026rsquo;-CAGTTGATTTACTTGGTTGCTGGT-3\u0026rsquo;; \u003cem\u003ehsp70\u003c/em\u003e common AGAGTG motif Fw 5\u0026rsquo;-TTCTCTGGCCGTTATTCGTT-3\u0026rsquo;, Rv 5\u0026rsquo;-TCGAACCAACGAGAGCAGTA-3\u0026rsquo;; hsp70 TSS Fw 5\u0026rsquo;-CGACATACTGCTCTCGTTGG-3\u0026rsquo;, Rv 5\u0026rsquo;-CAGCTGCGCTTGTTTATTTG-3\u0026rsquo;; \u003cem\u003ehsp70B\u003c/em\u003e 3\u0026rsquo; end Fw 5\u0026rsquo;-AATGGAATCCTGAACGTCAGC-3\u0026rsquo;, Rv 5\u0026rsquo;-CACATTGAAGACGTAGCTCTCC-3\u0026rsquo;. Additional primers for 5\u0026rsquo; and 3\u0026rsquo; were extracted from reference\u0026nbsp;[16]\u0026nbsp;and were named CHIP-A and CHIP-B, respectively.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eCo-immunoprecipitation\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eCoexpression of DAxud1-GFP and NELF-B-HA on the imaginal wing disc was achieved with driver nub\u0026gt;Gal4. In total, 150 wing discs were dissected on PBS 1X, then precipitated and resuspended in RIPA buffer. Protein extracts were immunoprecipitated with GFP antibody (Abcam ab290) with O.N incubation at 4\u0026deg;C, then isolated with Dynabeads\u0026reg; Protein-A. Western blot assay was performed in cells with INPUT (total protein), NB (non bound fraction, supernatant post Dynabeads isolation), and protein precipitated with Dynabeads. Page ruler 2166 (thermoscientific) was used as a weight marker.\u003c/p\u003e"},{"header":"Results","content":"\u003cp\u003e\u003cstrong\u003eDistribution of DAxud1 across the genome\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAccording to results from previous reports, vertebrate DAxud1 orthologs possess transcription factor-like features, including DNA binding, acidic and transcription transactivation domains [20, 21, 28]. They have also been detected bound to specific promoters related with neural crest differentiation [28]. Further, vertebrate and invertebrate orthologs of DAxud1 are related to the stress response and cancer through induction of apoptosis [21, 23], although, until now, this process has not been linked to its putative transcriptional function. Considering this information, we asked which groups or categories of genes are transcriptionally regulated and/or bind DAxud1 in \u003cem\u003eDrosophila\u003c/em\u003e, using this information as a proxy to get a better perspective on the function of \u003cem\u003edAxud1\u003c/em\u003e. For this purpose, we performed a TaDa-Seq (DamID-seq) experiment using the method described by Southall et al. [29], expressing the Dam-DAxud1 fusion protein to explore the \u003cem\u003eloci\u003c/em\u003e in which DAxud1 could reside. This was performed in imaginal wing disc tissue using the \u003cem\u003enub\u003c/em\u003e\u0026gt;Gal4 driver. Two replicates were generated per condition (third instar larvae imaginal wing disc, with Dam as a control or Dam-DAxud1). The aligned sequences were analyzed with the \u003cem\u003efindpeaks\u003c/em\u003e script, generated by Marshall Owen [30]. Using this method, we identified 1811 significant peaks across the \u003cem\u003eDrosophila\u0026nbsp;\u003c/em\u003egenome, representing zones where Dam-DAxud1 has a stronger, more stable positioning than the Dam signal. As shown in Fig. 1A, the Dam-DAxud1 distribution was mainly intronic, with 52.3% of peaks localizing to these gene regions. The first intron accounted for more than half of all instances of intron localization, representing 25.7% of overall genome occupancy. The other main significant locations were the intergenic and the proximal promoter regions (-1000 to +1). These regions are also the most enriched in transcription factor binding sites [31, 32]. To better understand the distribution of DAxud1 across the genome and in specific genes, we used the information from significant peaks to create a metagene profile, using gene bodies as a reference scaled to 1000 bp. The resulting metagene profile (Fig. 1B) shows a robust average signal at the promoter zone and a maximum within the 5\u0026rsquo;-half of gene bodies, confirming the previously described distribution (Fig. 1A) and suggesting a role in activating or poised gene transcription.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eSeveral reports describe DAxud1 vertebrate orthologs as stress response genes [24, 26, 33, 34]. Therefore, we searched for significant DAxud1 binding on stress response gene \u003cem\u003eloci\u003c/em\u003e (Additional Table S1) and found heat shock response genes as recurrent targets. Further, we generated alignment charts to visualize the distribution of peaks within this group of genes (Fig. 1C). These graphics clearly reveal that the distribution of Dam-DAxud1 specifically on \u003cem\u003ehsp\u003c/em\u003e genes has a highly similar pattern to that observed in the metagene profile.\u003c/p\u003e\n\u003cp\u003eThe intergenic, promoter, and gene-body (intron and exon) peaks were sorted and annotated for nearby genes using HOMER\u0026nbsp;[35]. Additional Table S2 shows the most relevant enrichment analysis results. Within the \u0026ldquo;biological processes\u0026rdquo; category, the two major enriched groups are related to neural development, but no groups show a significative score for stress response or tolerance. One pathway category identified involves Wnt signaling, concordant with reports that DAxud1 orthologs in mammals are linked to this signaling pathway\u0026nbsp;[19, 36, 37]\u0026nbsp;and apoptosis\u0026nbsp;[21, 38].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003ePotential DAxud1 DNA binding motifs\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAfter having found genes in which Dam-DAxud1 is enriched, we searched for recurrent binding motifs within the Dam-DAxud1 peaks. The sequence peaks had a size between 50-5000 bp and were extracted for \u003cem\u003emotif find analysis\u003c/em\u003e from HOMER [39]. As shown in Fig. 1D, the most recurrent DNA motif (TACATACATA), present in 1019 sequence peaks out of 1811, is novel compared with a previous report on Axud1 vertebrate orthologs [28], possibly due to the wide-range of sequences that the TaDa-seq experiment provides [40\u0026ndash;42], but this information could provide insights from the chromatin context surrounding DAxud1. Using the matrix sequence, we conducted another analysis using this data and the GOMO platform, which provides information about the gene ontology (GO) of the promoter regions, using a range between -1000 and +200 bp from the TSS for each gene in the \u003cem\u003eDrosophila melanogaster\u003c/em\u003e genome. For the TACATACATA motif, the most related GO found is \u0026ldquo;Heat shock mediated polytene chromosome puffing\u0026rdquo; (GO:0035080). This motif also found a match in the \u003cem\u003eDrosophila\u003c/em\u003e Topoisomerase 2 gene, which encodes a protein that removes supercoils in chromatin to facilitate transcription, with an essential role in RNA polymerase pausing-release in \u003cem\u003ehsps\u003c/em\u003e genes in the fast stress response [43, 44]. These analyses suggest that DAxud1 might play a role in the transcription of heat shock and other stress-induced genes. To confirm this last hypothesis, we performed chromatin immunoprecipitation on the \u003cem\u003ehsp70\u003c/em\u003e promoter with DAxud1-GFP as bait, followed by a CHIP-PCR with primers flanking the TACATACATA motif. The result confirms the presence of DAxud1 at this \u003cem\u003elocus\u003c/em\u003e (Additional Fig. S9), reinforcing the idea that it regulates the stress response by modulating \u003cem\u003ehsp\u003c/em\u003e gene expression and also that the regions discovered using TaDa-seq are \u003cem\u003ebona fide\u003c/em\u003e DAxud1 binding sites representing potentially functional interactions.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDAxud1 regulates thermotolerance and lifespan in Drosophila\u0026nbsp;\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eOur results reveal the presence of DAxud1 bound to \u003cem\u003ehsp\u003c/em\u003e genes (Fig. 1; Additional Fig. S9), suggesting a possible role in stress tolerance, a hypothesis supported by studies on Axud1 orthologs in mammals [23, 45]. \u0026nbsp; Although \u0026ldquo;stress response\u0026rdquo; is not the major GO category of DAxud1 genome occupancy (Additional Table S2), we tested whether flies with a global knockdown of this gene exhibit an altered stress response or aberrant \u003cem\u003ehsps\u003c/em\u003e expression, given its presence in this class of genes. For this, we used flies expressing a DAxud1 RNAi construct (Vienna stock 26479, UAS-IR-DAxud1) in all tissues using the Tubulin-Gal4 driver. Flies overexpressing DAxud1 cannot be assessed in this type of experiment since ubiquitous expression of this protein is lethal [21]. We established adult flies at 29\u0026deg;C for an optimal expression of the UAS RNAi construct; then, control and experimental adult flies were exposed to a heat shock of 37\u0026deg;C for a half-hour every day as the thermal stress condition, and their survival was measured. Knockdown of DAxud1 leads to a diminished lifespan of adult flies compared to the control genotype under the same stress condition (Fig. 2A-2B). Curiously, under control conditions (no heat shock), DAxud1 knockdown animals have a longer lifespan compared to control animals for both males and females (Figs. 2A and 2B). To confirm whether there is a fluctuation in \u003cem\u003ehsp\u003c/em\u003e gene expression, and to relate this to the observed phenotypes, we performed qPCR analysis for different\u003cem\u003e\u0026nbsp;hsp\u0026nbsp;\u003c/em\u003egenes in imaginal wing discs and salivary glands, expressing the RNAi construct or DAxud1-GFP version using the \u003cem\u003enubbin\u003c/em\u003e-Gal4 driver. This allowed us to study the effects of diminished DAxud1 levels in salivary glands and imaginal wing discs, tissues that have cells in an endo replication state (salivary glands) and in a mitotic state (imaginal wing discs). Figs. 2C, 2D, and 2E show the qPCR results for three \u003cem\u003ehsp\u0026nbsp;\u003c/em\u003egenes (\u003cem\u003ehsp70, hsp26, hsp67\u003c/em\u003e). In these experiments, we can appreciate that DAxud1 knockdown increases the expression of \u003cem\u003ehsp\u003c/em\u003e genes in control conditions, but does not potentiate the transcriptional heat shock response in either tissue. Notably, the rise of \u003cem\u003ehsp\u003c/em\u003e mRNA expression at the control temperature can explain the longer lifespan since there is evidence that \u003cem\u003ehsp\u003c/em\u003e overexpression generates this effect [46]. On the other hand, DAxud1\u003cem\u003e\u0026nbsp;\u003c/em\u003eoverexpression reduces \u003cem\u003ehsp\u003c/em\u003e basal expression at the control temperature, perhaps not inhibiting the \u003cem\u003ehsp\u003c/em\u003e stress response but generating a milder induction (Figs. 2C-E). These results, especially those observed after DAxud1 knockdown, are in conflict with our first hypothesis in which DAxud1 could function as a transcription factor that facilitates \u003cem\u003ehsp\u003c/em\u003e expression. Therefore, we then set out to further clarify the role of DAxud1 on the \u003cem\u003ehsp70 locus\u003c/em\u003e.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDAxud1 exhibits a widespread presence across the genome, whereas heat shock induces its recruitment to the \u003cem\u003ehsp70 locus\u003c/em\u003e\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eFurther examination of the results of \u003cem\u003ehsp\u003c/em\u003e gene expression in DAxud1 knockdown conditions (Fig. 2), suggests that DAxud1 may participate directly in \u003cem\u003ehsp\u003c/em\u003e transcription, possibly in a repressive manner. Previous studies have reported that DAxud1 orthologs (CSRNP in mouse) have transcription factor features, including a C-terminal trans-activating domain and a DNA binding motif [47]. However, there is no evidence indicating which type of regulatory factor it is, and its dynamics within the nuclear structure. \u003cem\u003eDrosophila melanogaster\u0026rsquo;s\u0026nbsp;\u003c/em\u003epolytene chromosomes represent a suitable model to answer these questions.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAs there is no available antibody against DAxud1, we used the GFP-tagged version of DAxud1 [21] (DAxud1-GFP) to perform immunofluorescence using an anti-GFP antibody. RNA Polymerase IIo (Hyper-phosphorylated) was used as a positive control for a chromosome attached protein and for detecting transcriptionally active \u003cem\u003eloci.\u0026nbsp;\u003c/em\u003eAs seen in Fig. 3A, DAxud1-GFP on the polytene chromosome exhibits a pattern with widespread RNA Pol IIo co-distribution, which means the specific function on \u003cem\u003ehsp\u003c/em\u003e genes may be distinct from other functions of DAxud1 in gene expression. The \u003cem\u003ehsp70\u003c/em\u003e \u003cem\u003elocus\u003c/em\u003e (Fig. 3B) has been extensively analyzed with regards to its chromatin rearrangements under stress conditions. Its documented that during stress conditions such as heat shock, cytological zones known as \u003cem\u003echromatin puffs\u0026nbsp;\u003c/em\u003eappear in \u003cem\u003ehsp70\u003c/em\u003e region,\u003cem\u003e\u0026nbsp;\u003c/em\u003eclearly\u003cem\u003e\u0026nbsp;\u003c/em\u003evisible\u003cem\u003e\u0026nbsp;\u003c/em\u003ewith confocal microscopy, a reflection of chromatin opening for transcription and transcription factor recruitment [48, 49]. We evaluated chromatin and DAxud1-GFP dynamics on this locus (at both 17\u0026deg;C and 37\u0026deg;C, Fig. 3C) and observed that DAxud1-GFP localizes to the 87A-87B \u003cem\u003elocus\u003c/em\u003e together with RNA Polymerase IIo after heat shock treatment. Transcription factors that relocate in this way are identified mostly as activators, for instance HSF, p-TEFb [2], and DSIF [16]. From these results, we reasoned that DAxud1 is a chromatin element that relocates to the \u003cem\u003ehsp70 locus\u003c/em\u003e during heat stress, but it does not interfere with the transcriptional induction since \u003cem\u003ehsps\u0026nbsp;\u003c/em\u003emRNA levels rise by the same magnitude as in control, under a DAxud1 overexpression condition (Fig. 2C-E).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTo confirm these results, we performed chromatin immunoprecipitation for DAxud1-GFP (Fig. 3D) to detect the presence of this factor on the \u003cem\u003ehsp70\u003c/em\u003e paralogs \u003cem\u003ehsp70Ab\u003c/em\u003e and \u003cem\u003ehsp70Bb/Bc\u003c/em\u003e, or \u003cem\u003eloci\u0026nbsp;\u003c/em\u003e87A and 87C, respectively. According to those results, DAxud1-GFP has a clear presence on \u003cem\u003ehsp70\u003c/em\u003e promoters as do other heat shock transcription factors [4]. Interestingly, DAxud1-GFP did not exhibit enrichment after oxidative stress induction by peroxide (Fig. 3D), a condition that increases \u003cem\u003ehsp\u003c/em\u003e mRNA levels [50], indicating that the involvement of DAxud1 in the heat shock response is different to the mechanism evoked during oxidative stress.\u003c/p\u003e\n\u003cp\u003eTogether, the data obtained here, have described the genome-wide distribution of DAxud1, its influence on \u003cem\u003ehsps\u003c/em\u003e transcription, and its effects on lifespan. However, how does DAxud1 exert its function and what is its relevance for cell physiology? In the \u003cem\u003ehsp70\u0026nbsp;\u003c/em\u003egene, there is a special chromatin configuration, known as the pausing complex, in which the activated RNA Polymerase II pauses the transcription process at positions +30 to +50 bp of the gene, and remains stalled as a \u0026ldquo;poised polymerase\u0026rdquo; [44]. This poised polymerase forms a complex with two main components, NELF and DSIF, whose role is to maintain the polymerase in the paused state. Transcription resumes when the p-TEFb complex phosphorylates NELF and DSIF, releasing the active polymerase for resumption of elongation [44]. \u003cem\u003ehsp70\u003c/em\u003e genes have been used as transcriptional pausing models because the resumption of transcription occurs quickly when heat shock or stress signals reach the cell [44]. Considering the data on the position and behavior of DAxud1 on the \u003cem\u003ehsp70\u003c/em\u003e gene after heat shock and that elements of the pausing complex on \u003cem\u003ehsp70\u003c/em\u003e\u0026nbsp; exhibit a similar pattern of reorganization on chromatin during heat shock [4, 18], we conjectured that DAxud1 might interact with the pausing complex.\u003c/p\u003e\n\u003cp\u003eThis hypothesis is also supported by RNA-seq evidence (Additional Figs. S3 and S4), in which DAxud1 overexpression increases \u003cem\u003ehsp70\u003c/em\u003e RNA synthesis only at the 5\u0026rsquo; end of the gene, with an incomplete synthesis of \u003cem\u003ehsp70\u003c/em\u003e mRNAs, specifically on \u003cem\u003ehsp70B\u003c/em\u003e paralogs, which display high reads only between +1 and +100 region. This effect was not detected in qPCRs performed in Figs. 2C-E because the primers were designed to detect the 3\u0026rsquo; side. In addition, this role in pausing is also supported by the evidence of double-hybrid experiments with \u003cem\u003eDrosophila\u003c/em\u003e proteins, which demonstrated a physical interaction between DAxud1 and Spt5 [51], a member of the DSIF complex, with a pivotal role in the RNA Polymerase transcriptional pausing and elongation complex [52]. Notably, the DAxud1 overexpression wing phenotype is reversed in a heterozygous \u003cem\u003espt5\u003c/em\u003e mutant background (Additional Fig. S6). Another component of the pausing machinery is the NELF complex (NELF-A, NELF-B, NELF-C/D, NELF-E), a stabilizer factor of the stalled polymerase in basal conditions. These components dissociate from the polymerase when transcription is resumed after heat shock [18], and the loss of function of some of its components generates higher levels of \u003cem\u003ehsp70\u003c/em\u003e [53], which is similar to our observations in the case of DAxud1 knockdown (Fig. 2C). Considering this information, we tested whether DAxud1 influences RNA Polymerase II location within the \u003cem\u003ehsp70\u003c/em\u003e gene body \u0026nbsp;[53]. Chromatin immunoprecipitation (CHIP) for RNA Pol II was performed under different DAxud1 levels of expression (Fig. 4A). In the RNA Pol II CHIP, we find that, in the case of DAxud1 knockdown, its levels decrease uniformly in the \u003cem\u003ehsp70\u003c/em\u003e gene \u003cem\u003elocus\u003c/em\u003e, and there is no change in the 5\u0026rsquo;/3\u0026rsquo; rate of RNA Pol II compared to the control and heat shock conditions (Fig. 4B). On the other hand, DAxud1 overexpression generates an increase of this ratio, more so than in the control condition, due to an increase of RNA Pol II on 5\u0026rsquo; end of the \u003cem\u003ehsp70\u003c/em\u003e gene body (Fig 4A).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eTo find out whether DAxud1 interacts with pausing factors, we performed a co-immunoprecipitation assay in imaginal wing disc expressing NELF-B-HA and DAxud1-GFP. In this assay (Fig. 4C), a fraction of NELF-B-HA was co-immunoprecipitated with DAxud1-GFP, indicating a physical interaction between these two proteins, suggesting they could act in the same process in \u003cem\u003ehsp70\u003c/em\u003e expression.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eThe information provided by the physical and genetic interaction of DAxud1 with NELF-B and Spt5, respectively, strongly suggests that it could be influencing RNA Polymerase II dynamics on the \u003cem\u003ehsp70\u003c/em\u003e gene. To confirm this, we performed a CHIP for RNA Polymerase on the 5\u0026rsquo; and 3\u0026rsquo; regions of the \u003cem\u003ehsp70\u003c/em\u003e gene, under DAxud1 overexpression and knockdown conditions. The results, presented in Fig. 4, reveal that DAxud1 overexpression generates an apparent pausing effect on RNA Polymerase II distribution (Fig. 4), coincident with extensive upregulation of the 5\u0026rsquo; RNA levels of \u003cem\u003ehsp70\u0026nbsp;\u003c/em\u003e(Additional Figs. S3 and S4).\u0026nbsp;\u003c/p\u003e\n\u003cp\u003ePromoter-proximal pausing in \u003cem\u003ehsp\u003c/em\u003e genes allows a rapid expression after heat shock induction due to pausing elements as DSIF/NELF complexes. Also, this rapid expression has a shut-off mechanism to restore \u003cem\u003ehsp70\u003c/em\u003e mRNA levels after the heat shock condition is reverted (heat shock recovery), and there is evidence that supports the role of NELF complex in the shut-off mechanism [54]. We tested the physical interaction between NELF-B and DAxud1 (Fig. 4C) as well as its genetic interaction (Additional Fig. S7), allowing us to further clarify the role of DAxud1 in \u003cem\u003ehsp70\u003c/em\u003e shut-down levels during heat shock recovery. To evaluate the effect of DAxud1 in this process, we performed qPCR for \u003cem\u003ehsp70\u003c/em\u003e mRNA 40\u0026rsquo; after a 20\u0026rsquo; heat shock (37\u0026deg;), as shown in Fig. 4D, in DAxud1 knockdown salivary glands. It is possible to observe that \u003cem\u003ehsp70\u003c/em\u003e induction has an apparently slower turn-over than in control animals (Fig. 4D), suggesting that DAxud1 plays a role in shutting-down \u003cem\u003ehsp70\u003c/em\u003e \u003cem\u003elocus\u0026nbsp;\u003c/em\u003einduction, and confirming DAxud1 has a synergistic function with NELF-B in the transcriptional pausing of \u003cem\u003ehsp70\u003c/em\u003e as well as the complex stabilization for further heat shock stimulus.\u003c/p\u003e"},{"header":"Discussion","content":"\u003cp\u003e\u003cstrong\u003eRole of DAxud1 in tissue homeostasis through \u003cem\u003ehsp\u003c/em\u003e regulation\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn the first part of this study, we searched for the most frequent sites where DAxud1 locates in the \u003cem\u003eDrosophila\u003c/em\u003e genome, and upon finding \u003cem\u003ehsp\u003c/em\u003e-type genes in these sites we performed thermotolerance analysis under DAxud1 knockdown in order to test the physiological relevance of this gene. We further studied \u003cem\u003ehsp\u003c/em\u003e expression during DAxud1 knockdown to mechanistically link these cellular players. In lifespan assays with DAxud1 knockdown (Figs. 2A-B), we found it generates adults with extended lifespan in control conditions, but reduced thermal resistance compared with control genotype animals, with no other change in developmental timing or apparent alterations (Fig. 1; Additional Fig. S2). Interestingly, in DAxud1 knockdown larvae raised at 29\u0026deg;C (no heat shock), the organisms exhibit an increase in the level of \u003cem\u003ehsp70\u003c/em\u003e expression in salivary glands and imaginal wing discs, compared to the control animals subjected to the same temperature. In \u003cem\u003eDrosophila\u003c/em\u003e, the effect of \u003cem\u003ehsps\u003c/em\u003e overexpression is well documented, in that one of its effects is precisely an extended lifespan due to its cytoprotective and anti-apoptotic effect\u0026nbsp;[46, 55], explaining why DAxud1 knockdown extends lifespan in control temperature, compared with animals only expressing the Gal4 driver (Figs. 2A-B). On the other hand, when animals with DAxud1 knockdown are exposed to daily thermal stress, lifespan shortens, situation that did not occur with control animals in which lifespan remained within the same range, as in control temperature, during daily thermal stress (Figs. 2A-B). This seems paradoxical since animals with higher levels of Hsps would be expected to maintain enhanced thermotolerance. However, the diminished tolerance to stress in DAxud1 knockdown could be explained alternatively by the reported pro-apoptotic DAxud1 function through an activating role on JNK signaling in a DAxud1 overexpression background\u0026nbsp;[21, 26]. Thus, in the context of DAxud1 knockdown, the JNK signaling pathway could be less strongly activated, allowing damaged cells to survive though interfering with development and regenerative processes, which are known consequences of impaired JNK signaling\u0026nbsp;[56, 57]. This last suggestion is supported by the occupancy of Dam-DAxud1 on the \u003cem\u003eloci\u003c/em\u003e of genes associated with the GO- KEEG apoptosis pathway, related to JNK (Additional Tables S2 and S5). Regarding this point, the pro-apoptotic activity of DAxud1 is associated with control of tumorigenesis\u0026nbsp;[19, 21]\u0026nbsp;and cancers with poor prognosis\u0026nbsp;[58], which may be due not only to decreased pro-apoptotic activity resulting from lower Axud1 levels, but also to increased \u003cem\u003ehsp\u003c/em\u003e gene expression. \u003cem\u003ehsp\u003c/em\u003e overexpression strongly correlates with cancer cell progression and poor prognosis\u0026nbsp;[59, 60], \u0026nbsp;so this condition of certain cancer cells may be due to low levels of Axud1, a condition also associated with aggressive cancers, as previously described\u0026nbsp;[19, 21, 58]. This function in tissue homeostasis through transcription modulation could be one of the reasons DAxud1 is conserved among metazoans. However, the reduced thermotolerance phenotype would not be directly related to the expression of \u003cem\u003ehsps\u003c/em\u003e, despite the fact that DAxud1 exhibits recurrence in these genes (Fig. 1) and at the same time the strong effect on their expression generated by the knockdown of DAxud1. These observations prompted us to investigate the relationship between \u003cem\u003ehsps\u003c/em\u003e, widely conserved in all kingdoms, and DAxud1, which only has orthologs in metazoans, probably participating in a type of modulation on \u003cem\u003ehsps\u0026nbsp;\u003c/em\u003eexpression exclusive to metazoans. This led us to focus on its effects at the chromatin level.\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDAxud1 function in chromatin as pausing factor on the \u003cem\u003ehsp70\u003c/em\u003e gene\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eAlthough \u003cem\u003edAxud1\u003c/em\u003e mRNA levels do not change significantly during heat shock (Additional Fig. S8), the DAxud1 protein relocates to the \u003cem\u003ehsp70\u003c/em\u003e \u003cem\u003eloci\u003c/em\u003e. The localization of DAxud1 to chromatin, observed in Fig. 3C, can mediated by putative phosphorylation sites\u0026nbsp;[21]\u0026nbsp;or by the DAxud1 cysteine-rich region\u0026nbsp;[20, 21]. These cysteine-rich regions have been shown to act as modulators of protein conformation, in which heat or oxidative stress can alter the disulfide bonds, changing the conformational state of the protein. This phenomenon is well documented\u0026nbsp;[61, 62]\u0026nbsp;in proteins with cysteine-rich regions, including the DAxud1 CSRNP family\u0026nbsp;[20]. Thus, changes in DAxud1 protein conformation induced by post-translational modifications could be key in the regulation of its localization and its effect on gene expression.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;Transcriptional pausing is a state in which the RNA Polymerase initiates transcription but remains stalled on the first 30-50 base pairs of the gene, forming a complex with the DSIF (Spt5/Stp4 proteins) and NELF complexes\u0026nbsp;[16, 48, 63]. The pausing is released when internal or external cell signals activate p-TEFb, which phosphorylates RNA Polymerase II itself as well as the DSIF, and NELF complexes which, in turn, disengage from the RNA Polymerase\u0026nbsp;[18, 64, 65]. Also, there is evidence that pausing relies on the chromatin state, in which NELF exerts it pausing effect by depleting H3K4me3 histones and, therefore, when NELF dissociates from the pausing complex, the polymerase elongates in a favorable transcriptional context \u0026nbsp;[53]. This mechanism is found in fast-response genes like \u003cem\u003ehsp\u0026nbsp;\u003c/em\u003egenes or in innate immune response genes\u0026nbsp;[66]. Strikingly, overexpression of DAxud1 generates an increase in \u003cem\u003ehsp70B\u003c/em\u003e transcription only on its 5\u0026rsquo; end of the transcript sequence (Additional Figs. S3 and S4).\u0026nbsp;This increment of 5\u0026rsquo; sequences of the \u003cem\u003ehsp70B\u003c/em\u003e genes indicates that there is an increase of mRNA synthesis, from approximately +1 to +80 bp. We interpret that DAxud1 retains the pausing complex completely stalled, generating an enrichment of 5\u0026rsquo; \u003cem\u003ehsp70\u003c/em\u003e mRNA. Whatever the explanation is at this level, the effect of DAxud1 on this particular \u003cem\u003ehsp\u003c/em\u003e gene expression could involve the transcriptional pausing mechanism, in which DAxud1 might interact with the pausing complex.\u003c/p\u003e\n\u003cp\u003eThe possible role of DAxud1 in the pausing complex is also supported by the co-immunoprecipitation of DAxud1-GFP with NELF-B (Fig. 4C). NELF-B acts as a component of\u0026nbsp;the\u0026nbsp;pausing complex, stabilizing and pausing RNA Polymerase II in a hypo-phosphorylated CTD state, and dissociates from the complex under heat shock or another signal \u0026nbsp;[54]. Additionally, data concerning Spt5, another component of\u0026nbsp;the pausing complex (DSIF complex), further supports this observation: Biogrid\u0026nbsp;[51], a database that documents physical interactions between proteins, reported an interaction of DAxud1 with Spt5\u0026nbsp;[4]\u0026nbsp;and, importantly,\u0026nbsp;the DAxud1 overexpression phenotype is partially reverted in a \u003cem\u003espt5\u003c/em\u003e heterozygous mutant background (Additional Fig. S6). Furthermore, Spt5 interacts with NELF-B in\u0026nbsp;the pausing complex\u0026nbsp;[67], and stabilizes it synergistically with the NELF complex, maintaining the RNA Polymerase stalled. This pausing complex keeps stable until p-TEFb kinase phosphorylates both complexes, with dissociation of NELF complex\u0026nbsp;[68]\u0026nbsp;and DSIF (Stp5/6) remaining as a component of RNA Polymerase holoenzyme, acting as an elongation factor. In this context, DAxud1 may play a role in pausing the stabilization complex, \u0026nbsp;since DAxud1\u003cem\u003e\u0026nbsp;\u003c/em\u003eknockdown causes an increase in the transcription rate of \u003cem\u003ehsp70\u003c/em\u003e mRNA, while its overexpression causes stalled transcription at the 5\u0026rsquo; region, perhaps instigating\u0026nbsp;the formation of\u0026nbsp;the pausing complex (Fig. 4) in\u0026nbsp;the same way as NELF does\u0026nbsp;[53], but not interfering with RNA Polymerase II pausing release in heat shock conditions. This interaction could be conserved in other instances, such as in the regulation of the expression of the MMP1 gene in humans, which requires the ortholog of DAxud1 (CSRNP1) for its transcription upon cytokine activation\u0026nbsp;[69]. More importantly, there is detailed evidence that MMP1 is a transcriptionally paused gene, in which the NELF/DSIF pausing complex is necessary for\u0026nbsp;its appropriate expression in response to immune-activating signals\u0026nbsp;[66]. Also, in myeloid cells, the NELF complex and Spt5 are essential for maintaining the repression of pro-apoptotic genes during myeloid development in zebrafish. In this case, they rely on their role in RNA Pol II pausing of pro-apoptotic genes\u0026nbsp;[70]\u0026nbsp;and, accordingly, Spt5\u0026rsquo;s deficiency generates a loss of\u0026nbsp;the myeloid line, due to early apoptosis of hematopoietic stem cells, the same phenotype observed after Csrnp1 reduction in morphant zebrafish larvae\u0026nbsp;[71].\u0026nbsp;Curiously, in imaginal wing discs and, subsequently, in the adult wing, the overexpression of NELF-B reverts the DAxud1 overexpression phenotype, while NELF-B knockdown exacerbates the DAxud1 overexpression phenotype (Additional Fig. S7). These results suggest an antagonistic role of these two proteins on gene expression, though not necessarily related to transcription of \u003cem\u003ehsp\u003c/em\u003e genes, because the knockdown of both genes separately provokes a similar effect in the increase of \u003cem\u003ehsp\u003c/em\u003e transcription. Instead, such an effect could be found in genes related to apoptosis, in which NELF-B functions as a repressive element, as previously described\u0026nbsp;[53]. This is also similar to the interaction between DSIF and NELF, both are necessary for RNA Pol II pausing, but after p-TEFb activation (through heat or morphogen signals) their roles are quite opposite, in which the DSIF elements are necessary for RNA Pol II elongation while NELF dissociates completely from the transcription complex\u0026nbsp;[18].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eThe relation between pausing and DAxud1 motif founded with TaDa-seq\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eIn support of the pausing-release hypothesis, we found that the most enriched motif found in TaDa-seq peaks is \u0026ldquo;TACATACATACA.\u0026rdquo; This sequence strongly matches the binding site of Topoisomerase 2 (Top2), a chromatin modifying protein (Fig. 1D). There is evidence that Top2 plays a role in pausing-release and in the \u003cem\u003ehsp70\u003c/em\u003e fast transcriptional response in \u003cem\u003eDrosophila\u003c/em\u003e Kc cells\u0026nbsp;[72]. In addition, GOMO analysis shows that this motif is present in promoters of genes that belong to the \u0026ldquo;heat shock-mediated polytene chromosome puff\u0026rdquo; (GO:0035080) and \u0026ldquo;response to hypoxia\u0026rdquo; (GO:0001666) GO categories. Additionally, there is evidence that Topoisomerase 2, due to its structural chromatin function, plays a role in the cell cycle, in which inhibition or missteps in the arrangement of this protein lead to cell cycle arrest in G2-M\u0026nbsp;[73]. The same cellular phenotype was previously reported when DAxud1 was overexpressed\u0026nbsp;[21], raising the possibility that these proteins interact at the chromatin topology level, likely interfering with Topoisomerase 2 positioning before the G2-M checkpoint. However, more analysis is required to support this conjecture.\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eIn previous studies, the vertebrate orthologues of DAxud1 were associated with the binding motif AGAGTS\u0026nbsp;[28, 69, 74]. Curiously, with HOMER (\u003cem\u003edenovo\u003c/em\u003e \u003cem\u003emotif discovery\u003c/em\u003e function), no similar sequence was found. We thus used HOMER to search for the occurrence of the AGAGTS motif among the peak sequences; we found AGATGS 4657 times in 1407 peaks among the 1811 significant sequences. In total, 895 genes were identified (Additional Table S4). Again, pathways like Wnt and Apoptosis are present in this group (Additional Tables S2 and S5), most of those belonging to the apoptosis group are part of JNK signaling, including \u003cem\u003emsn\u003c/em\u003e, \u003cem\u003epuc\u003c/em\u003e,\u003cem\u003e\u0026nbsp;rpr\u003c/em\u003e, and \u003cem\u003ehid\u003c/em\u003e, a finding\u003cem\u003e\u0026nbsp;\u003c/em\u003ewhich is in agreement with \u003cem\u003ebsk\u0026nbsp;\u003c/em\u003eand\u003cem\u003e\u0026nbsp;puc\u0026nbsp;\u003c/em\u003eupregulation upon DAxud1 overexpression (Additional Fig. S5,\u0026nbsp;[21]). This result indicates that DAxud1 acts at different levels on transcriptional regulation to balance cell physiology.\u003c/p\u003e\n\u003cp\u003eWith regard to its physiological role, low levels of human Axud1 expression are related to tumor development\u0026nbsp;[19]. Considering the role of Axud1 in the stress response, this condition could predispose cells to undergo more damage in the first stage of tumor development due to the loss of \u003cem\u003ehsp\u003c/em\u003e gene regulation, leading to higher levels of Hsp proteins, thereby buffering damage signals and likely inhibiting the apoptotic process and extending the lifespan of impaired cells. Later on, in the next stage, Axud1 could be necessary at higher expression levels to initiate apoptosis through JNK signaling\u0026nbsp;[21], though this Axud1 upregulation does not occur in tumorigenesis\u0026nbsp;[19].\u003c/p\u003e\n\u003cp\u003e\u003cstrong\u003eDevelopmental role of Axud1\u003c/strong\u003e\u003c/p\u003e\n\u003cp\u003eA previous study showed that the DAxud1 vertebrate ortholog, Axud1, acts as a transcription factor, located only in the promoter region and related to the positive transcriptional activity of the Pax7 and Msx1 transcription factors, driving the development of neural crest by Wnt signaling\u0026nbsp;[28]. According to this study, the presence of Axud1 on these promoters depends on\u0026nbsp;the AGAGTS binding site, previously described for human and mouse orthologs [11], from \u003cem\u003ein vitro\u003c/em\u003e experiments. The AGAGTS motif was found in the TaDa-seq data, with a frequency of 4657 times out of 1407 unique peak sequences. The AGAGTS sequence is present 88,167 times in the whole \u003cem\u003eDrosophila\u003c/em\u003e genome and, therefore, the Dam-DAxud1 peak sequences comprise 5.28% of the whole genome motif\u0026rsquo;s abundance. Out of those 1400 peaks, gene ontology analysis mainly points towards morphogenesis, axon guidance, and Wnt signaling (Additional Tables S2 and S5). This information is in line with previous work on vertebrate Axud1\u0026nbsp;[19, 28, 75]\u0026nbsp;and our data in which Axin and Wnt are downregulated upon DAxud1 overexpression (Additional Fig. S5), which could indicate that DAxud1 is part of a feedback loop in the Wnt pathway. Mainly, the data from other studies supports the role of Axud1 in neural development as previously described, in which vertebrate homologs of DAxud1 could exert a pro-neural function in both CNS and neural crest development, and this depends on wnt signaling\u0026nbsp;[28, 76]. Therefore, this regulatory network is conserved in animals with different roles in development. This function in neural structures is a hallmark of DAxud1 orthologs because our previous data showed that zebrafish Csrnp1a is highly expressed in the growing brain\u0026nbsp;[77]\u0026nbsp;and in the ventral nerve cord in Drosophila during germ band elongation stage\u0026nbsp;[21], making this protein a regulator of a specific process more than a general factor of animal development, and not classifiable in a binary category of transcription activator/repressor.\u003c/p\u003e\n\u003cp\u003eOverall, our study raises new insights about DAxud1 function, demonstrating new chromatin features related to the stress response. We propose this is mainly supported on the stabilization of\u0026nbsp;the pausing complex, maintaining a chromatin configuration for rapid \u003cem\u003ehsp\u003c/em\u003e expression under stress conditions in the first stage of the thermal stress response. This function might not represent the complete chromatin regulatory function of DAxud1, which could play a role in JNK signaling in long-term stress to promote apoptosis, and an additional role in neurogenesis as a mediator of Wnt signaling, as described in vertebrates [78].\u003c/p\u003e"},{"header":"Declarations","content":"\u003cp\u003eEthics approval:\u003c/p\u003e\n\u003cp\u003eOur laboratory and the manage of animals were approved to perform the experiments exposed in this article. This was approved by Bioethics committee of the Faculty of Sciences (Universidad de Chile), chaired by Dr. Marco M\u0026eacute;ndez from 2013.\u003c/p\u003e\n\u003cp\u003eCompeting interests:\u003c/p\u003e\n\u003cp\u003eNo potential competing interest was reported by the authors.\u003c/p\u003e\n\u003cp\u003eAuthor\u0026rsquo;s contributions:\u003c/p\u003e\n\u003cp\u003eJM Z\u0026uacute;\u0026ntilde;iga-Hern\u0026aacute;ndez was the main executor, writer of the original draft, designer of the experiments, as well as the data analysis, including NGS data. Also, JM Z\u0026uacute;\u0026ntilde;iga-Hern\u0026aacute;ndez is the corresponding author. Meneses C contributed to facilitate the access to Illumina sequencing platform, and support experimental design. Bast\u0026iacute;as M was the technician in charge to process the libraries for the Illumina sequencing. Allende ML contributed to the editing and review of the manuscript, funding and provider of resources. Glavic A was the supervisor and reviewer of the manuscript, funding provider, as well as the manager of the initial idea for this research.\u003c/p\u003e\n\u003cp\u003eAuthor information:\u003c/p\u003e\n\u003cp\u003eJorge Z\u0026uacute;\u0026ntilde;iga Hern\u0026aacute;ndez, PhD (Z\u0026uacute;\u0026ntilde;iga-Hern\u0026aacute;ndez JM):\u0026nbsp;Corresponding author, works as Postdoctoral researcher in\u0026nbsp;Center for Genome Regulation and Laboratory of Developmental Biology, Department of Biology, Faculty of Sciences, University of Chile. Santiago, Chile. Mail: [email protected]\u003c/p\u003e\n\u003cp\u003eClaudio Meneses, PhD (Meneses C): Investigator in Plant Biotechnology Center, University Andres Bello. Santiago, Chile. Mail: [email protected]\u003c/p\u003e\n\u003cp\u003eMacarena Bast\u0026iacute;as (Bast\u0026iacute;as M): Technician in in Plant Biotechnology Center, University Andres Bello. Santiago, Chile. Mail:\u0026nbsp;[email protected]\u003c/p\u003e\n\u003cp\u003eMiguel Allende, PhD (Allende ML): Director of Center for Genome Regulation and Laboratory of Developmental Biology, Department of Biology, Faculty of Sciences, University of Chile. Santiago, Chile. Mail: [email protected]\u003c/p\u003e\n\u003cp\u003e\u0026nbsp;\u003c/p\u003e\n\u003cp\u003eAlvaro Glavic, PhD (Glavic A): Co investigator in Center for Genome Regulation and Laboratory of Developmental Biology, Department of Biology, Faculty of Sciences, University of Chile. Santiago, Chile. Mail: [email protected]\u003c/p\u003e\n\u003cp\u003eConsent of participation:\u003c/p\u003e\n\u003cp\u003eNo humans were part of the experiments exposed in this article. Consent of publication:\u003c/p\u003e\n\u003cp\u003eNo humans were part of the experiments exposed in this article.\u003c/p\u003e\n\u003cp\u003eFunding:\u003c/p\u003e\n\u003cp\u003eAnillo ACT1401 / Center for Genome Regulation (FONDAP 15200002)\u003c/p\u003e\n\u003cp\u003eAvailability of data:\u003c/p\u003e\n\u003cp\u003eTaDa-seq are available in\u0026nbsp;Sequence Read Archive (SRA), ID: PRJNA776616\u003c/p\u003e\n\u003cp\u003eAcknowledgements:\u003c/p\u003e\n\u003cp\u003eWe thank to Andrea Brand for provide plasmid pUAST-attb-LT3-Dam. Also acknowledgments to FONDECYT (ANID) for Doctoral Fellowship 21110721.\u003c/p\u003e"},{"header":"References","content":"\u003col\u003e\u003cli\u003e\u003cspan\u003eRichter K, Haslbeck M, Buchner J. The Heat Shock Response: Life on the Verge of Death. Mol Cell. 2010;40:253\u0026ndash;66. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.molcel.2010.10.006\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBoehm AK, Saunders A, Werner J, Lis JT. Transcription Factor and Polymerase Recruitment, Modification, and Movement on dhsp70 In Vivo in the Minutes following Heat Shock. Mol Cell Biol. 2003;23:7628\u0026ndash;37.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMayer MP, Bukau B. Hsp70 chaperones: Cellular functions and molecular mechanism. Cell Mol Life Sci. 2005;62:670\u0026ndash;84.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAndrulis ED. High-resolution localization of Drosophila Spt5 and Spt6 at heat shock genes in vivo: roles in promoter proximal pausing and transcription elongation. Genes Dev. 2000;14:2635\u0026ndash;49. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1101/gad.844200\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eO\u0026rsquo;Brien T, Lis JT. Rapid changes in Drosophila transcription after an instantaneous heat shock. Mol Cell Biol. 1993;13:3456\u0026ndash;63.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSan Gil R, Ooi L, Yerbury JJ, Ecroyd H. The heat shock response in neurons and astroglia and its role in neurodegenerative diseases. Mol Neurodegener. 2017;12:1\u0026ndash;20.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChen B, Feder ME, Kang L. Evolution of heat-shock protein expression underlying adaptive responses to environmental stress. Mol Ecol. 2018;27:3040\u0026ndash;54.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLang BJ, Guerrero ME, Prince TL, Okusha Y, Bonorino C, Calderwood SK. The functions and regulation of heat shock proteins; key orchestrators of proteostasis and the heat shock response. Arch Toxicol. 2021;95:1943\u0026ndash;70. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s00204-021-03070-8\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDonovan MR, Marr MT. DFOXO activates large and small heat shock protein genes in response to oxidative stress to maintain proteostasis in drosophila. J Biol Chem. 2016;291:19042\u0026ndash;50.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLi Z, Srivastava P. Heat-shock proteins. Curr Protoc Immunol. 2004;Appendix 1:Appendix 1T. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1002/0471142735.ima01ts58\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKing AM, MacRae TH. Insect heat shock proteins during stress and diapause. Annu Rev Entomol. 2015;60:59\u0026ndash;75. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1146/annurev-ento-011613-162107\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKwak H, Fuda NJ, Core LJ, Lis JT. Precise maps of RNA polymerase reveal how promoters direct initiation and pausing. Science. 2013;339:950\u0026ndash;3. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1126/science.1229386\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDong B, Jaeger AM, Thiele DJ. Inhibiting Heat Shock Factor 1 in Cancer: A Unique Therapeutic Opportunity. Trends Pharmacol Sci. 2019;40:986\u0026ndash;1005.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eVihervaara A, Duarte FM, Lis JT. Molecular mechanisms driving transcriptional stress responses. Nat Rev Genet. 2018;19:385\u0026ndash;97. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/s41576-018-0001-6\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJennings BH, Shah S, Yamaguchi Y, Seki M, Phillips RG, Handa H, et al. Locus-Specific Requirements for Spt5 in Transcriptional Activation and Repression in Drosophila. 2004;14:1680\u0026ndash;4.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWu CH, Yamaguchi Y, Benjamin LR, Horvat-Gordon M, Washinsky J, Enerly E, et al. NELF and DSIF cause promoter proximal pausing on the hsp70 promoter in Drosophila. Genes Dev. 2003;17:1402\u0026ndash;14.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eO\u0026rsquo;Brien T, Lis JT. RNA polymerase II pauses at the 5\u0026rsquo; end of the transcriptionally induced Drosophila hsp70 gene. Mol Cell Biol. 1991;11:5285\u0026ndash;90. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1128/MCB.11.10.5285\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYamaguchi Y, Shibata H, Handa H. Transcription elongation factors DSIF and NELF: Promoter-proximal pausing and beyond. Biochim Biophys Acta - Gene Regul Mech. 2013;1829:98\u0026ndash;104. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.bbagrm.2012.11.007\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eIshiguro H, Tsunoda T, Tanaka T, Fujii Y, Nakamura Y, Furukawa Y. Identification of AXUD1, a novel human gene induced by AXIN1 and its reduced expression in human carcinomas of the lung, liver, colon and kidney. Oncogene. 2001;20:5062\u0026ndash;6. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1038/sj.onc.1204603\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGingras S, Pelletier S, Boyd K, Ihle JN. Characterization of a family of novel cysteine- serine-rich nuclear proteins (CSRNP). PLoS One. 2007;2:e808. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0000808\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGlavic A, Molnar C, Cotoras D, de Celis JF. Drosophila Axud1 is involved in the control of proliferation and displays pro-apoptotic activity. Mech Dev. 2008;126:184\u0026ndash;97. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.mod.2008.11.005\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003ePasovic L, Eidet JR, Olstad OK, Chen DF, Lyberg T, Utheim TP. Impact of Storage Temperature on the Expression of Cell Survival Genes in Cultured ARPE-19 Cells. Curr Eye Res. 2016;3683 June:1\u0026ndash;11. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.3109/02713683.2016.1145236\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYi X, Bekeredjian R, DeFilippis NJ, Siddiquee Z, Fernandez E, Shohet R V. Transcriptional analysis of doxorubicin-induced cardiotoxicity. Am J Physiol Heart Circ Physiol. 2006;290:H1098-102. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1152/ajpheart.00832.2005\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCheng Z, Zhao H, Ze Y, Su J, Li B, Sheng L, et al. Gene-expression changes in cerium chloride-induced injury of mouse hippocampus. PLoS One. 2013;8:e60092. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0060092\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcDermott JE, Vartanian KB, Mitchell H, Stevens SL, Sanfilippo A, Stenzel-Poore MP. Identification and validation of ifit1 as an important innate immune bottleneck. PLoS One. 2012;7:e36465. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0036465\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDiercke K, Kohl A, Lux CJ, Erber R. Compression of human primary cementoblasts leads to apoptosis: A possible cause of dental root resorption? J Orofac Orthop = Fortschritte der Kieferorthopadie Organ/official J Dtsch Gesellschaft fur Kieferorthopadiedg. 2014;:1\u0026ndash;16. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s00056-014-0237-5\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRundqvist HC, Montelius A, Osterlund T, Norman B, Esbjornsson M, Jansson E. Acute sprint exercise transcriptome in human skeletal muscle. PLoS One. 2019;14:1\u0026ndash;24. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0223024\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSim\u0026otilde;es-costa M, Stone M, Bronner ME. Axud1 integrates Wnt signaling and transcriptional inputs to drive neural crest formation.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSouthall TD, Gold KS, Egger B, Davidson CM, Caygill EE, Marshall OJ, et al. Cell-type-specific profiling of gene expression and chromatin binding without cell isolation: assaying RNA Pol II occupancy in neural stem cells. Dev Cell. 2013;26:101\u0026ndash;12. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.devcel.2013.05.020\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMarshall OJ, Brand AH. Damidseq-pipeline: An automated pipeline for processing DamID sequencing datasets. Bioinformatics. 2015;31:3371\u0026ndash;3.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHernandez-Garcia CM, Finer JJ. Identification and validation of promoters and cis-acting regulatory elements. Plant Sci. 2014;217\u0026ndash;218:109\u0026ndash;19. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.plantsci.2013.12.007\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJo B-S, Choi SS. Introns: The Functional Benefits of Introns in Genomes. Genomics Inform. 2015;13:112\u0026ndash;8. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.5808/GI.2015.13.4.112\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYi X, Bekeredjian R, DeFilippis NJ, Siddiquee Z, Fernandez E, Shohet R V. Transcriptional analysis of doxorubicin-induced cardiotoxicity. Am J Physiol Heart Circ Physiol. 2006;290:H1098-102. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1152/ajpheart.00832.2005\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGui S, Sang X, Zheng L, Ze Y, Zhao X, Sheng L, et al. Intragastric exposure to titanium dioxide nanoparticles induced nephrotoxicity in mice, assessed by physiological and gene expression modifications. Part Fibre Toxicol. 2013;10:4. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1186/1743-8977-10-4\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHeinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol Cell. 2010;38:576\u0026ndash;89. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.molcel.2010.05.004\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFang D, Li Soungyu L, Xu Wanfu, KE zhiyong ZF. CSRNP. J South Med Univ. 2013;33:1122\u0026ndash;6.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eContreras O, Cruz-Soca M, Theret M, Soliman H, Tung LW, Groppa E, et al. Cross-talk between TGF-β and PDGFRα signaling pathways regulates the fate of stromal fibro-adipogenic progenitors. J Cell Sci. 2019;132.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDiercke K, Kohl A, Lux CJ, Erber R. Compression of human primary cementoblasts leads to apoptosis A possible cause of dental root resorption ? Kompression f\u0026uuml;hrt in prim\u0026auml;ren humanen Zementoblasten zur Apoptose Eine m\u0026ouml;gliche Ursache f\u0026uuml;r Wurzelresorptionen. 2014;:430\u0026ndash;45.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eHeinz S, Benner C, Spann N, Bertolino E, Lin YC, Laslo P, et al. Simple Combinations of Lineage-Determining Transcription Factors Prime cis-Regulatory Elements Required for Macrophage and B Cell Identities. Mol Cell. 2010;38:576\u0026ndash;89. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.molcel.2010.05.004\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGermann S, Juul-Jensen T, Letarnec B, Gaudin V. DamID, a new tool for studying plant chromatin profiling in vivo, and its use to identify putative LHP1 target loci. Plant J. 2006;48:153\u0026ndash;63.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSouthall TD, Gold KS, Egger B, Davidson CM, Caygill EE, Marshall OJ, et al. Cell-type-specific profiling of gene expression and chromatin binding without cell isolation: assaying RNA Pol II occupancy in neural stem cells. Dev Cell. 2013;26:101\u0026ndash;12. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.devcel.2013.05.020\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAughey GN, Southall TD. Dam it\u0026rsquo;s good! DamID profiling of protein-DNA interactions. Wiley Interdiscip Rev Dev Biol. 2016;5:25\u0026ndash;37.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJohnson DG, Walker CL. Cyclins and cell cycle checkpoints. Annu Rev Pharmacol Toxicol. 1999;39:295\u0026ndash;312. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1146/annurev.pharmtox.39.1.295\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBunch H. RNA polymerase II pausing and transcriptional regulation of the HSP70 expression. Eur J Cell Biol. 2017;96:739\u0026ndash;45. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.ejcb.2017.09.003\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKorb K, Katsikogianni E, Zingler S, Daum E, Lux CJ, Hohenstein A, et al. Inhibition of AXUD1 attenuates compression-dependent apoptosis of cementoblasts. Clin Oral Investig. 2016;20:2333\u0026ndash;41. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1007/s00784-016-1740-4\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhao Y, Sun H, Lu J, Li X, Chen X, Tao D, et al. Lifespan extension and elevated hsp gene expression in Drosophila caused by histone deacetylase inhibitors. J Exp Biol. 2005;208:697\u0026ndash;705.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGingras S, Pelletier S, Boyd K, Ihle JN. Characterization of a family of novel cysteine- serine-rich nuclear proteins (CSRNP). PLoS One. 2007;2:e808. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0000808\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMissra A, Gilmour DS. Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation complex. Proc Natl Acad Sci U S A. 2010;107:11301\u0026ndash;6. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1073/pnas.1000681107\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhimulev IF, Belyaeva ES, Semeshin VF, Koryakov DE, Demakov S a, Demakova O V, et al. Polytene chromosomes: 70 years of genetic research. Int Rev Cytol. 2004;241:203\u0026ndash;75. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/S0074-7696(04)41004-3\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eCourgeon A-M, Rollet E, BEcket J, Maisonhaute C, Best-Belpomme M. Hydrogen peroxide (H2O2) induces actin and some heat-shock proteins in Drosophila cells. Eur J Biochem. 1988;171:163\u0026ndash;70. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1111/j.1432-1033.1988.tb13772.x\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eChatr-Aryamontri A, Breitkreutz BJ, Oughtred R, Boucher L, Heinicke S, Chen D, et al. The BioGRID interaction database: 2015 update. Nucleic Acids Res. 2015;43:D470\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSaunders A, Werner J, Andrulis ED, Nakayama T, Hirose S, Reinberg D, et al. Tracking FACT and the RNA polymerase II elongation complex through chromatin in vivo. Science. 2003;301:1094\u0026ndash;6.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGilchrist DA, Nechaev S, Lee C, Ghosh SKB, Collins JB, Li L, et al. NELF-mediated stalling of Pol II can enhance gene expression by blocking promoter-proximal nucleosome assembly. Genes Dev. 2008;22:1921\u0026ndash;33.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGhosh SKB, Missra A, Gilmour DS. Negative Elongation Factor Accelerates the Rate at Which Heat Shock Genes Are Shut off by Facilitating Dissociation of Heat Shock Factor. Mol Cell Biol. 2011;31:4232\u0026ndash;43.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eS\u0026oslash;rensen JG, Loeschcke V. Larval crowding in Drosophila melanogaster induces Hsp70 expression, and leads to increased adult longevity and adult thermal stress resistance. J Insect Physiol. 2001;47:1301\u0026ndash;7.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBosch M, Serras F, Mart\u0026iacute;n-Blanco E, Bagu\u0026ntilde;\u0026agrave; J. JNK signaling pathway required for wound healing in regenerating Drosophila wing imaginal discs. Dev Biol. 2005;280:73\u0026ndash;86. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.ydbio.2005.01.002\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhu S, Chen R, Soba P, Jan YN. JNK signaling coordinates with ecdysone signaling to promote pruning of Drosophila sensory neuron dendrites. Dev. 2019;146:1\u0026ndash;10.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eZhang H, Qiu X, Yang G. The CSRNP Gene Family Serves as a Prognostic Biomarker in Clear Cell Renal Cell Carcinoma. 2021;11 March:1\u0026ndash;12.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSherman MY, Gabai VL. Hsp70 in cancer: back to the future. 2015; October 2014:4153\u0026ndash;61.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJagadish N, Agarwal S, Gupta N, Fatima R, Devi S, Kumar V, et al. Heat shock protein 70-2 (HSP70-2) overexpression in breast cancer. J Exp Clin Cancer Res. 2016;2:1\u0026ndash;14.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMcDuffee AT, Senisterra G, Huntley S, Lepock JR, Sekhar KR, Meredith MJ, et al. Proteins containing non-native disulfide fonds generated by oxidative stress can act as signals for the induction of the heat shock response. J Cell Physiol. 1997;171:143\u0026ndash;51.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMosaddegh B, Takalloo Z, Sajedi RH, Shirin Shahangian S, Hassani L, Rasti B. An inter-subunit disulfide bond of artemin acts as a redox switch for its chaperone-like activity. Cell Stress Chaperones. 2018;23:685\u0026ndash;93.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eWu CH, Lee C, Fan R, Smith MJ, Yamaguchi Y, Handa H, et al. Molecular characterization of Drosophila NELF. Nucleic Acids Res. 2005;33:1269\u0026ndash;79.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSims RJ, Mandal SS, Reinberg D. Recent highlights of RNA-polymerase-II-mediated transcription. Curr Opin Cell Biol. 2004;16:263\u0026ndash;71.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eLis JT, Mason P, Peng J, Price DH, Werner J. P-TEFb kinase recruitment and function at heat shock loci P-TEFb kinase recruitment and function at heat shock loci. 2000;:792\u0026ndash;803.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGilchrist DA, Fromm G, dos Santos G, Pham LN, Mcdaniel IE, Burkholder A, et al. Regulating the regulators: The pervasive effects of Pol II pausing on stimulus-responsive gene networks. Genes Dev. 2012;26:933\u0026ndash;44.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMissra A, Gilmour DS. Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation complex. Proc Natl Acad Sci U S A. 2010;107:11301\u0026ndash;6. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1073/pnas.1000681107\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMissra A, Gilmour DS. Interactions between DSIF (DRB sensitivity inducing factor), NELF (negative elongation factor), and the Drosophila RNA polymerase II transcription elongation complex. Proc Natl Acad Sci U S A. 2010;107:11301\u0026ndash;6. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1073/pnas.1000681107\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMacdonald CD, Falconer AMD, Chan CM, Wilkinson DJ, Skelton A, Reynard L, et al. Cytokine-induced cysteine- serine-rich nuclear protein-1 (CSRNP1) selectively contributes to MMP1 expression in human chondrocytes. PLoS One. 2018;13:1\u0026ndash;16.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYang Q, Liu X, Zhou T, Cook J, Nguyen K, Bai X. RNA polymerase II pausing modulates hematopoietic stem cell emergence in zebrafish. Blood. 2016;128:1701\u0026ndash;10.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eSolı C, Espina J, Feijo CG. csrnp1a Is Necessary for the Development of Primitive Hematopoiesis Progenitors in Zebrafish. 2013;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eKroeger PE, Rowe TC. Analysis of Topoisomerase I and II Cleavage Sites on the Drosophila Actin and Hsp70 Heat Shock Genes. Biochemistry. 1992;31:2492\u0026ndash;501.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eDownes CS, Clarke DJ, Mullinger AM, Gim\u0026eacute;nez-Abi\u0026aacute;n JF, Creighton AM, Johnson RT. A topoisomerase II-dependent G2 cycle checkpoint in mammalian cells. Nature. 1994;372:467\u0026ndash;70.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eGingras S, Pelletier S, Boyd K, Ihle JN. Characterization of a family of novel cysteine- serine-rich nuclear proteins (CSRNP). PLoS One. 2007;2:e808. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pone.0000808\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eYamada K, Akiyama N, Yamada S, Tanaka H, Saito S, Hiraoka M, et al. Taip2 is a novel cell death-related gene expressed in the brain during development. Biochem Biophys Res Commun. 2008;369:426\u0026ndash;31. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.bbrc.2008.02.041\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFeij\u0026oacute;o CG, Sarrazin AF, Allende ML, Glavic A. Cystein-serine-rich nuclear protein 1, Axud1/Csrnp1, is essential for cephalic neural progenitor proliferation and survival in zebrafish. Dev Dyn. 2009;238:2034\u0026ndash;43. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1002/dvdy.22006\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eFeij\u0026oacute;o CG, Sarrazin AF, Allende ML, Glavic A. Cystein-serine-rich nuclear protein 1, Axud1/Csrnp1, is essential for cephalic neural progenitor proliferation and survival in zebrafish. Dev Dyn. 2009;238:2034\u0026ndash;43. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1002/dvdy.22006\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eAzambuja AP, Simoes-Costa M. A regulatory sub-circuit downstream of Wnt signaling controls developmental transitions in neural crest formation. PLoS Genet. 2021;17:1\u0026ndash;23. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1371/journal.pgen.1009296\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eJohansen KM, Cai W, Deng H, Bao X, Zhang W, Girton J, et al. Polytene chromosome squash methods for studying transcription and epigenetic chromatin modification in Drosophila using antibodies. Methods. 2009;48:387\u0026ndash;97. doi:\u003cspan class=\"ExternalRef\"\u003e\u003cspan class=\"RefSource\"\u003e10.1016/j.ymeth.2009.02.019\u003c/span\u003e\u003c/span\u003e.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eMarshall OJ, Brand AH. Damidseq-pipeline: An automated pipeline for processing DamID sequencing datasets. Bioinformatics. 2015;31:3371\u0026ndash;3.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eBailey TL, Boden M, Buske FA, Frith M, Grant CE, Clementi L, et al. MEME Suite: Tools for motif discovery and searching. Nucleic Acids Res. 2009;37 SUPPL. 2:202\u0026ndash;8.\u003c/span\u003e\u003c/li\u003e \u003cli\u003e\u003cspan\u003eRam\u0026iacute;rez F, D\u0026uuml;ndar F, Diehl S, Gr\u0026uuml;ning BA, Manke T. DeepTools: A flexible platform for exploring deep-sequencing data. Nucleic Acids Res. 2014;42:187\u0026ndash;91.\u003c/span\u003e\u003c/li\u003e\u003c/ol\u003e"}],"fulltextSource":"","fullText":"","funders":[],"hasAdminPriorityOnWorkflow":false,"hasManuscriptDocX":true,"hasOptedInToPreprint":true,"hasPassedJournalQc":"","hasAnyPriority":false,"hideJournal":true,"highlight":"","institution":"","isAcceptedByJournal":false,"isAuthorSuppliedPdf":false,"isDeskRejected":"","isHiddenFromSearch":false,"isInQc":false,"isInWorkflow":false,"isPdf":false,"isPdfUpToDate":true,"isWithdrawnOrRetracted":false,"journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true},"keywords":"heat shock, hsp70 transcription, pausing complex, DAxud1, NELF","lastPublishedDoi":"10.21203/rs.3.rs-1040684/v1","lastPublishedDoiUrl":"https://doi.org/10.21203/rs.3.rs-1040684/v1","license":{"name":"CC BY 4.0","url":"https://creativecommons.org/licenses/by/4.0/"},"manuscriptAbstract":"\u003cp\u003e\u003cstrong\u003eBackground\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eA rapid transcriptional response under an acute stimulus is common in all cellular systems and is an adaptation that allows tolerance to environmental changes. A gene group that has been studied because of its fast response and cytoprotective effects are the\u0026nbsp;\u003cem\u003ehsp\u003c/em\u003e\u0026nbsp;genes (encodingHeat Shock Proteins(HSPs), conserved chaperones).. Under normal conditions, the mRNA and protein levels of the main \u003cem\u003ehsp\u003c/em\u003e genes are low but they increase rapidly upon heat shock (HS). This is achieved due to the presence of an RNA Polymerase II pausing complex located +30-50 bp from TSS. This complex maintains a partially synthesized RNA strand of said length, poised to resume synthesis, and undergoes subsequent transcriptional inactivation to restore transcript levels after environmental stabilization.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eMethods\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eThe Gal4/UAS system was used to modify \u003cem\u003edAxud1\u003c/em\u003e expression in a tissue specific manner. A DAxud1-GFP fusion was expressed in salivary glands to perform polytene chromosome immunofluorescence and chromatin immunoprecipitation. DAxud1 genome occupancy data was achieved expressing Dam-DAxud1 in imaginal wing discs using Gal4/UAS (TaDa-seq).\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eResults\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eUsing TaDa-seq, we demonstrate that DAxud1 protein is present mainly near the TSS of significant occupied genes, most frequently in the first intron. This results also revealed DAxud1 is present in \u003cem\u003ehsp\u003c/em\u003e genes, mainly in promoter zone. Following these results, we found that, under \u003cem\u003edAxud\u003c/em\u003e knockdown, larvae and adults flies have a diminished thermotolerance, despite showing an increase in \u003cem\u003ehsp\u003c/em\u003e transcripts in larval tissues. We performed polytene chromosome immunofluorescence for DAxud1-GFP, revealing extensive, but dynamic localization on chromatin in \u003cem\u003ehsp70 loci\u003c/em\u003e. This \u0026nbsp;was confirmed with chromatin immunoprecipitation. We also found that DAxud1 overexpression leads to an enrichment of RNA Polymerase II at the 5’ end of the \u003cem\u003ehsp70 \u003c/em\u003egene, with a decrease in its transcripts. Importantly, we show interaction of DAxud1 with NELF-B, a component of the transcriptional pausing complex, and knockdown of both genes individually has similar effects on \u003cem\u003ehsp70\u003c/em\u003e transcription.\u003c/p\u003e\u003cp\u003e\u003cstrong\u003eConclusion\u003c/strong\u003e\u003c/p\u003e\u003cp\u003eDAxud1 protein is a component of chromatin, that relocates under stress conditions such as heat shock, playing a role in maintaining RNA Polymerase II stalled at the 5’ of \u003cem\u003ehsp70, \u003c/em\u003epossibly\u003cem\u003e \u003c/em\u003ethrough a pausing mechanism based on its interaction with NELF-B.\u003c/p\u003e","manuscriptTitle":"Drosophila DAxud1: A New Element in Transcriptional Pausing Complex Stabilization","msid":"","msnumber":"","nonDraftVersions":[{"code":1,"date":"2021-11-11 15:19:43","doi":"10.21203/rs.3.rs-1040684/v1","editorialEvents":[{"type":"communityComments","content":0}],"status":"published","journal":{"display":true,"email":"[email protected]","identity":"researchsquare","isNatureJournal":false,"hasQc":true,"allowDirectSubmit":true,"externalIdentity":"","sideBox":"","snPcode":"","submissionUrl":"/submission","title":"Research Square","twitterHandle":"researchsquare","acdcEnabled":true,"dfaEnabled":false,"editorialSystem":"","reportingPortfolio":"","inReviewEnabled":false,"inReviewRevisionsEnabled":true}}],"origin":"","ownerIdentity":"b4396af8-5172-4dde-889a-677b88ea987f","owner":[],"postedDate":"November 11th, 2021","published":true,"recentEditorialEvents":[],"rejectedJournal":[],"revision":"","amendment":"","status":"posted","subjectAreas":[{"id":8450309,"name":"Epigenetics \u0026 Genomics"}],"tags":[],"updatedAt":"2021-12-16T06:56:28+00:00","versionOfRecord":[],"versionCreatedAt":"2021-11-11 15:19:43","video":"","vorDoi":"","vorDoiUrl":"","workflowStages":[]},"version":"v1","identity":"rs-1040684","journalConfig":"researchsquare"},"__N_SSP":true},"page":"/article/[identity]/[[...version]]","query":{"redirect":"/article/rs-1040684","identity":"rs-1040684","version":["v1"]},"buildId":"7rjqhiLT3MXkJMwkYKINL","isFallback":false,"isExperimentalCompile":false,"dynamicIds":[84888],"gssp":true,"scriptLoader":[]}