Antagonistic roles for Ataxin-2 structured and disordered domains in RNP condensation

  1. Amanjot Singh
  2. Joern Hulsmeier
  3. Arvind Reddy Kandi
  4. Sai Shruti Pothapragada
  5. Jens Hillebrand
  6. Arnas Petrauskas
  7. Khushboo Agrawal
  8. Krishnan RT
  9. Devasena Thiagarajan
  10. Deepa Jayaprakashappa
  11. K VijayRaghavan
  12. Mani Ramaswami  Is a corresponding author
  13. Baskar Bakthavachalu  Is a corresponding author
  1. National Centre for Biological Sciences, India
  2. Trinity College Institute of Neuroscience, School of Genetics and Microbiology, Smurfit Institute of Genetics and School of Natural Sciences, Trinity College Dublin, Ireland
  3. Tata Institute for Genetics and Society Centre at inStem, Bellary Road, India
  4. School of Biotechnology, Amrita Vishwa Vidyapeetham University, India
  5. School of Basic Sciences, Indian Institute of Technology, India
6 figures, 2 tables and 3 additional files


Figure 1 with 4 supplements
Using Targets of RNA-Binding Proteins Identified by Editing (TRIBE) to identify Atx2-interacting mRNAs in the adult fly brain.

(A) Schematic and flow chart for TRIBE analysis: the Atx2-ADARcd fusion protein is expressed in the adult fly brain, total brain mRNA isolated, sequenced, and analyzed using a published TRIBE pipeline. (B) Heatmaps show edit percentages of individual transcript coordinates. Replicate experiments R1 and R2 identify largely overlapping edit sites and edited mRNAs. The common targets (intersect between R1 and R2) show almost reproducible edit levels. The inset (heatmaps on either sides of the common intersected list) indicates that several mRNAs are identified as ‘non-replicates’ between R1 or R2 because they do not cross quality control thresholds (edit percentages or read counts) and not because of robust differences between replicates. (C) Bar plot showing the number of edits with a significant number of genes edited at a single site. (D) The mRNA edits are independent of expression levels. A scatter plot shows expression differences of all the mRNAs expressed in fly brain. Red dots represent the edited mRNAs while gray dots represent the mRNAs of the brain transcriptome that are not edited. (E) Sanger sequencing confirms editing site identified by TRIBE analysis. Data shown for one target mRNA (Akh). Red bars show the edit percentages at the different modified nucleotides with respect to the total Akh mRNA shown in gray. The black bar indicates identified edits that are below 15% threshold (see also Figure 1—figure supplement 4).

Figure 1—source data 1

Data related to Figure 1B.

Percentage edits for Atx2-target mRNAs in RNAseq replicates.
Figure 1—source data 2

Data related to Figure 1C.

Atx2-target mRNAs along with the number of edits.
Figure 1—figure supplement 1
Atx2-ADARcd fusion protein constructs.

Schematic showing the Atx2-ADARcd construct.

Figure 1—figure supplement 2
UAS-Atx2-ADARcd expression was driven using elav-Gal4;tub-Gal80ts in adult fly brain for 5 days.

Brains dissected from flies maintained at 29℃ for 5 days were stained using anti-V5 (A, A') and anti-Atx2 (B, B') antibodies. (A-C') Confocal images of Oblique frontal view of the Drosophila brain. (A,B,C) (dorsal is more anterior) of a complete. (A',B'C') higher magnification of the respective overview images.

Figure 1—figure supplement 3
Unique mapping percentage for various RNAseq replicates.

Total RNA was isolated from 10 to 12 brains for each replicate. Poly(A)-enriched mRNA was sequenced using paired-end libraries on Illumina 2500. The reads obtained were of high quality, and mapping to Drosophila melanogaster genome was >80% for all the samples.

Figure 1—figure supplement 4
An example of a gene, 14-3-3 epsilon, showing multiple edits specifically in the exons.

Red lines show the edit percentages at the different modified nucleotides with respect to the total mRNA present at a given moment of time.

Figure 2 with 4 supplements
Atx2 preferentially edits AU-rich sequences in the 3′UTRs of the target mRNAs.

(A) Metaplot analysis showing Atx2 preferentially associates with 3′UTRs of the target mRNA. (B) Motif analyses involving ±100 bases around the edit site using MEME identifies AU-rich element sequences in the 3′UTRs of the target mRNAs.

Figure 2—figure supplement 1
Exon specificity of Atx2 mRNA interactions and sequence selectivity of observed Atx2–coding sequence (CDS) interactions.

Atx2-ADARcd fusion protein exclusively edits exons, indicating cytoplasmic binding to mature mRNAs.

Figure 2—figure supplement 2
Atx2-target mRNAs have longer-than-average 3′UTRs.

Mann–Whitney U test was used to determine the difference in the 3′UTR length between Atx2 targets and brain transcriptome. p-Value is <0.0001, and error bars represent SEM. Bootstraping was performed 1000 times with replacement. Statistics was performed using chi-square test with cutoff (p-value <0.05). Fisher's test was used to combine the p-values.

Figure 2—figure supplement 3
Coding sequence edits were co-existent with 3'UTR edits.

(A) Coding sequence edits often occurred together with a 3′UTR edit on the same mRNA consistent with an initial Atx2–3′UTR interaction enabling edits to a proximal coding sequence as well. (B) A GC-rich motif for Atx2–CDS association indicated by analyses of sequences flanking CDS edit sites.

Figure 2—figure supplement 4
Multiple-sequence alignment was performed for regions around Atx2 binding motifs in Akh and Vmat genes between Drosophila melanogaster and closely related Drosophila species.
Figure 3 with 6 supplements
Atx2–ARE interactions modulate mRNA stability.

(A) Schematic of strategy to induce RNAi expression specifically in adult fly brain. Total RNA extracted from brain post 5 days at 29°C was sequenced using Illumina 2500 and differential expression analyzed using DESeq2. (B) Normalized mRNA read counts showing Atx2 mRNA levels to be significantly reduced following RNAi expression compared to Gal4 control with a p-value of <0.0298 using Student’s t-test. (C) Effect on Atx2-target mRNAs following Atx2 knockdown: the majority are reduced in level, indicating a role for Atx2 is target stabilization. Bootstraping was performed 1000 times with replacement. Statistics was performed using chi-square test with cutoff (p-value <0.05). Fisher’s test was used to combine the p-values. (D) AREScore analysis showing AU-rich elements (AREs) to be enriched in Atx2-target mRNAs compared to the brain transcriptome. (E) Among Atx2-target mRNAs, higher ARE scores are seen in mRNAs whose levels are reduced following Atx2 knockdown.

Figure 3—source data 1

Data related to Figure 3C.

Fold changes in Atx2-target mRNAs upon Atx2 RNAi.
Figure 3—source data 2

Data related to Figure 3E.

Correlation of AREScore and target mRNA levels.
Figure 3—figure supplement 1
Atx2 RNAi knockdown reduces Atx2 expressions.

Ptc-Gal4 was used to drive Atx2 RNAi in Drosophila wing disc and immunostained with anti-beta-Gal (A) and anti-Atx2 (B) antibodies. Atx2 levels are reduced in the Gal4-expressing cells (C).

Figure 3—figure supplement 1—source data 1

Data related to RNAseq mapping percentage from Gal4 and Atx2 RNAi brains.
Figure 3—figure supplement 2
Correlation between AREScore and relative expression levels of Atx2-downregulated targets for three biological replicates (selected based on a p-value<0.05).

Heatmap indicates relative expression levels for each mRNA in control and experimental brain samples. Light gray bars represent Atx2-target mRNAs without significant AREScore.

Figure 3—figure supplement 3
Differential expression analysis of intron reads from control and Atx2 RNAi show no major reduction in levels of Atx2 Targets of RNA-Binding Proteins Identified by Editing nascent mRNA targets.
Figure 3—figure supplement 3—source data 1

Fold changes in the introns of Atx2-target mRNAs upon Atx2 RNAi.
Figure 3—figure supplement 4
Comparative analyses of Atx2 targets in Drosophila brain and S2 cells.

(A) Comparison of Atx2 S2 cell Targets of RNA-Binding Proteins Identified by Editing (TRIBE) targets with Drosophila brain TRIBE targets. (B) Metaplot showing Atx2 preferentially associates with 3′UTRs of target mRNAs in S2 cells.

Figure 3—figure supplement 4—source data 1

Percentage edits of Atx2-target mRNAs in S2 cells among RNAseq replicates.
Figure 3—figure supplement 5
Altering levels of Atx2 in S2 cells affects the levels of target mRNAs.

(A) Knockdown of Atx2 using dsRNA in S2 cells. Western blots with anti-Atx2 antibody and anti-tubulin show significant reduction in Atx2 levels. (B) Knockdown of Atx2 in S2 cells destabilizes majority of its Targets of RNA-Binding Proteins Identified by Editing targets, which are stabilized in overexpression, indicating a role for Atx2 is target stabilization.

Figure 3—figure supplement 5—source data 1

Fold changes in Atx2-target mRNAs when Atx2 is overexpressed or knocked down in S2 cells.
Figure 3—figure supplement 6
AREScore correlate directly with the stability of Atx2-target mRNAs in S2 cells.
Figure 4 with 3 supplements
Atx2 requires its cIDR domain for interaction with its target mRNA.

(A) An illustration of domains of Atx2. The structured (Like-Sm [Lsm] and Lsm-associated [LsmAD]) and disordered domains (mIDR and cIDR) of Atx2 protein were deleted one at a time to understand the domains necessary for its interaction with target RNA. The wild-type Atx2 and the deletions contained a c-terminal ADARcd and V5 tag. The deletions are shown by the dotted lines. (B) Deletion of LsmAD or mIDR domain reduced the ability of Atx2 to interact with its targets, while cIDR deletion almost entirely prevented Atx2 mRNA interactions. The Lsm-domain deletion showed an overall increase in mRNA edits. (C) More detailed heatmap view of individual target mRNAs showing target edit percentages in flies expressing control or domain-deleted forms Atx2 fused to ADARcd. (D) Most of the apparently new targets identified in Atx2∆Lsm Targets of RNA-Binding Proteins Identified by Editing analyses are also edited, albeit at a lower efficiency, in flies expressing full-length Atx2-ADAR fusions. This suggests that deletion of the Lsm domain increases the interaction of Atx2 with its native target mRNAs. (For C and D, 15% threshold edits identified for samples in column 1 [Atx2 WT in C and Atx2∆Lsm in D] was used to compare with all the edits of the rest of the columns to generate a heatmap.)

Figure 4—source data 1

Data related to Figure 4B.

mRNA targets identified for different Atx2 domain deletions.
Figure 4—source data 2

Data related to Figure 4C, D.

Edit percentages for Atx2-target mRNAs identified using different domain deletions.
Figure 4—figure supplement 1
Atx2 Like-Sm (Lsm), Lsm-associated (LsmAD), or mIDR domains are not essential for AU-rich motif preference.

MEME analysis showing Atx2 preferentially associates with AU-rich motif on target mRNA. Deletion of Lsm, LsmAD, or mIDR domains does not block the Atx2 AU-rich motif preference. Unique targets found in Lsm domain deletion have a GC-rich motif preference, which is also identified as the second most predominant motif that Atx2 WT bind.

Figure 4—figure supplement 1—source data 1

Mapping percentage for different Atx2 deletions.
Figure 4—figure supplement 2
Atx2∆Lsm unique targets were present at a reduced edit threshold in Atx2 WT.

Several of the Atx2∆Lsm edits found at 15% threshold were found in Atx2 WT edits at 5% threshold.

Figure 4—figure supplement 3
Equivalent read quality and ADAR expression across Atx2 deletion analyses.

Poly(A)-enriched mRNA from 10 to 12 brains was sequenced using paired-end libraries on Illumina 2500. (A) The reads obtained were comparable between genotypes and of high quality, and mapping to the Drosophila genome was >80% for all the samples (Figure 4—figure supplement 1—source data 1). (B) Quantification for ADAR mRNA as a read out for transgene expression showed expression above the endogenous control and comparable between Atx2 WT and the various mutant transgenic lines. (C) Western blots with anti-V5 antibody show levels of Atx2 transgene expression in Drosophila brain among transgenic lines.

Figure 5 with 1 supplement
Atx2 Like-Sm (Lsm) domain alters granule dynamics.

(A) Drosophila S2 cells expressing Atx2 protein with a C-terminal SNAP tag. Atx2∆Lsm form large cytoplasmic granules (i) that are much larger than WT Atx2 granules (ii). Atx2 with an additional Lsm domain in place of mIDR to create 2XLsm forms smaller granules compared to WT Atx2 in S2 cells (iii). Deletion of cIDR domain blocks Atx2 granule formation (iv). In the absence of cIDR, Lsm deletion does not rescue Atx2 granules (v). Domain map along with the granule phenotype is shown in (vi). The scale bar corresponds to 2 µm. (B) Radius of the stress granule (size) and (C) number of stress granules per cell were quantified and plotted. More than 80 cells per genotype were used for granule quantification.

Figure 5—figure supplement 1
Granules formed by overexpression of Atx2 and ∆Lsm resemble stress granules (SGs) that also contain SG protein Rasputin (Rin).

Atx2 WT forms cytoplasmic granules (i) that contain Rin (ii and iii). The deletion of cIDR blocks granule formation (iv, v, and vi). Atx2 granules formed upon Like-Sm (Lsm) domain deletion also resemble Atx2 WT granules and contain SG protein Rin (vii, viii, and ix). The deletion of both Lsm and cIDR domains blocks the ability of Atx2 to form granules (x, xi, and xii). Atx2-2XLsm, where mIDR was replaced with a second Lsm domain, reduces Atx2 granule size but resembles WT granules and contains Rin (xiii, xiv, and xv). The scale bars correspond to 2 µm.

Ataxin-2 domains in mRNA target regulation.

A model to explain the current and previous observations on Ataxin-2 as well as its function in mRNA regulation. (A) Ataxin-2 does not make direct contact with mRNAs in soluble mRNPs. Instead, it is recruited to RNA by other RBPs that bind to structured Like-Sm (Lsm) (or Lsm-associated [LsmAD]) domains. In these soluble mRNPs, the cIDR is buried and inaccessible. One class of mRNP (above) contains AU-rich elements (AREs); the second class below (e.g., per in Drosophila) does not. (B) Under specific signaling conditions (e.g., stress), RBP–Atx2 interactions are prevented. In these ‘remodeled’ mRNPs, the Atx2-cIDR is exposed. We speculate that a segment of this intrinsically disordered domain (IDR) directly binds to nearby AREs (or to ARE-binding proteins), while other segments of the IDR mediate multivalent interactions that contribute to mRNP condensation. (C) Ataxin-2 cIDR interactions enable mRNP assembly into granules, facilitated by RNA–RNA crosslinks and interactions mediated by other RBPs (e.g., G3BP/Rin). These RNP granules may include both ARE-containing mRNAs (edited) and mRNAs that do not. Note that additional RBPs (red) could determine not only Atx2 proximity to specific mRNAs (leading to editing) but also the effect of Atx2 on mRNA translation and/or stability. These and alternative models remain to be tested.


Table 1
Gene Ontology (GO) analysis of Atx2 brain targets.

GO analysis using PANTHER shows a large enrichment of neuronal mRNAs coding for neuropeptides and proteins involved in neuronal signaling pathways. FDR: false discovery rate.

Molecular functionFold enrichmentProbability (FDR)Targets (%)Genome (%)GO term
Signaling receptor binding7.754.47E-1114.461.87GO:0005102
Neuropeptide hormone activity27.672.52E-096.630.24GO:0005184
Hormone activity18.645.74E-086.630.36GO:0005179
Neuropeptide receptor binding38.742.08E-064.220.11GO:0071855
G protein-coupled receptor binding13.843.79E-066.030.44GO:0001664
Receptor regulator activity9.584.55E-067.230.76GO:0030545
Receptor ligand activity9.411.96E-056.630.71GO:0048018
Signaling receptor activator activity9.042.58E-056.630.74GO:0030546
Organic hydroxy compound transmembrane transporter activity83.023.38E-042.410.03GO:1901618
Monoamine transmembrane transporter activity66.415.52E-042.410.04GO:0008504
Phosphoric ester hydrolase activity4.936.43E-047.841.59GO:0042578
Anion:sodium symporter activity55.356.44E-042.410.05GO:0015373
Sodium:chloride symporter activity55.356.84E-042.410.05GO:0015378
Phosphoprotein phosphatase activity7.957.06E-045.430.69GO:0004721
Anion:cation symporter activity36.91.81E-032.410.07GO:0015296
Cation:chloride symporter activity36.91.91E-032.410.07GO:0015377
Phosphatase activity5.052.53E-036.631.32GO:0016791
Serotonin:sodium symporter activity83.024.17E-031.810.03GO:0005335
Neurotransmitter:sodium symporter activity14.825.04E-033.020.21GO:0005328
Protein tyrosine/serine/threonine phosphatase activity23.726.17E-032.410.11GO:0008138
Solute:sodium symporter activity11.531.26E-023.020.27GO:0015370
mRNA 3'UTR binding11.221.36E-023.020.27GO:0003730
Protein tyrosine phosphatase activity10.641.59E-023.020.29GO:0004725
Neurotransmitter transmembrane transporter activity10.641.65E-023.020.29GO:0005326
Adrenergic receptor activity24.913.98E-021.810.08GO:0004935
Key resources table
Reagent type
(species) or
DesignationSource or
Genetic reagent (Drosophila melanogaster)UAS-Atx2-WT-ADARcdThis paperN/ARelated to Figures 1, 2 and 4
Genetic reagent (Drosophila melanogaster)UAS-Atx2-ΔLsm-ADARcdThis paperN/ARelated to Figure 4
Genetic reagent (Drosophila melanogaster)UAS-Atx2- ΔLsm-AD-ADARcdThis paperN/ARelated to Figure 4
Genetic reagent (Drosophila melanogaster)(UAS-Atx2- ΔmIDR-ADARcd)This paperN/ARelated to Figure 4
Genetic reagent (Drosophila melanogaster)UAS-Atx2- ΔcIDR-ADARcdThis paperN/ARelated to Figure 4
Genetic reagent (Drosophila melanogaster)UAS-Atx2-only Lsm/Lsm-AD -ADARcdThis paperN/ARelated to Figure 4
Genetic reagent (Drosophila melanogaster)Elav-Gal4; Tub-Gal80tsBloomington Drosophila Stock centerRelated to Figure 3
Genetic reagent (Drosophila melanogaster)UAS-Atx2 RNAiVienna Drosophila RNAi Center stock collection34955Related to Figure 3
Cell line
(Drosophila melanogaster)
Recombinant DNA reagentpJFRC7-20XUAS-IVS-8_Atx2-ADARcd
This paperN/AConstruct to express WT Atx2-ADAR fusion protein
Recombinant DNA reagentpJFRC7-20XUAS-
This paperN/AConstruct to express Atx2∆mIDR-ADAR fusion protein
Recombinant DNA reagentpJFRC7-20XUAS-IVS-8_Atx2ΔcIDR -ADARcd
This paperN/AConstruct to express Atx2∆cIDR-ADAR fusion protein
Recombinant DNA reagentpJFRC7-20XUAS-IVS-8_Atx2ΔLsm -ADARcd
This paperN/AConstruct to express Atx2∆Lsm-ADAR fusion protein
Recombinant DNA reagentpJFRC7-20XUAS-IVS-8_Atx2ΔLsmAD-ADARcd
This paperN/AConstruct to express Atx2∆LsmAD-ADAR fusion protein
Recombinant DNA reagentpJFRC7-20XUAS-IVS-8_Atx2only-Lsm-LsmAD-ADARcd
This paperN/AConstruct to express only Lsm-LsmAD-ADAR fusion protein
AntibodyAnti-Atx2 (chicken polyclonal)Bakthavachalu et al., 2018IF (1:500)
WB (1:1000)
AntibodyAnti-Rasputin (rabbit polyclonal)Aguilera-Gomez et al., 2017IF (1:500)
AntibodyAnti-GFP (rabbit polyclonal)Molecular probesCat# A11122IF (1:500)
AntibodyAnti-GFP (chicken polyclonal)AbcamCat# mAb 13970IF (1:500)
AntibodyAnti-V5 (rabbit polyclonal)Santa Cruz BiotechnologyCat# sc83849-RIF (1:500)
WB (1:1000)
AntibodyAnti-nc82 (mouse monoclonal)Kittel, 2006IF (1:500)
AntibodyAnti-tubulin (mouse monoclonal)DSHBCat# E7CWB (1:2000)
AntibodyAlexa Fluor 555 (polyclonal goat anti-chicken IgG)InvitrogenCat# A21437IF (1:1000)
AntibodyAlexa Fluor 488 (polyclonal goat anti-chicken IgG)InvitrogenCat# A11039IF (1:1000)
AntibodyAlexa Fluor 647 (polyclonal goat anti-chicken IgG)InvitrogenCat# A21449IF (1:1000)
AntibodyAlexa Fluor 555 (polyclonal goat anti-rabbit IgG)InvitrogenCat# A21428IF (1:1000)
AntibodyAlexa Fluor 488 (polyclonal goat anti-rabbit IgG)InvitrogenCat# A11078IF (1:1000)
AntibodyAlexa Fluor 647 (polyclonal goat anti-rabbit IgG)InvitrogenCat# A21244IF (1:1000)
AntibodyAlexa Fluor 555 (polyclonal goat anti-mouse IgG)InvitrogenCat# A21422IF (1:1000)
AntibodyAlexa Fluor 488
(polyclonal goat anti-mouse IgG)
InvitrogenCat# A21121IF (1:1000)
AntibodyAlexa Fluor 647 (polyclonal goat anti-mouse IgG)InvitrogenCat# A21235IF (1:1000)
Chemical compoundVectashield Mounting MediumVector LaboratoriesCat# H-1000
Chemical compoundSNAP-Surface 549New England BiolabsCat# S9112SIF (1:500)
Software, algorithmTRIBEMcMahon et al., 2016
Software, algorithmSTAR v2.5.3Dobin et al., 2013
Software, algorithmHTSeq v0.11.2Anders et al., 2015
Software, algorithmDESeq2Love et al., 2014
Software, algorithmAREScoreSpasic et al., 2012
Software, algorithmGuitarCui et al., 2016
Software, algorithmBedtoolsQuinlan and Hall, 2010
Software, algorithmtwoBitToFa-
Software, algorithmMEME suiteBailey et al., 2009
Software, algorithmCellprofilerMcQuin et al., 2018
Software, algorithmImageJSchneider et al., 2012
Software, algorithmGgplot2Wilkinson, 2011
Software, algorithmPheatmap
Software, algorithmSnapDragon

Additional files

Supplementary file 1

Details of all the primers used in the study.
Supplementary file 2

All the edits identified (without thresholding) in Targets of RNA-Binding Proteins Identified by Editing experiments are consolidated as a single spreadsheet file.

Each genotype is presented as a sheet within this file, and the raw and processed sequencing data is available at GEO accession number GSE153985.
Transparent reporting form

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Amanjot Singh
  2. Joern Hulsmeier
  3. Arvind Reddy Kandi
  4. Sai Shruti Pothapragada
  5. Jens Hillebrand
  6. Arnas Petrauskas
  7. Khushboo Agrawal
  8. Krishnan RT
  9. Devasena Thiagarajan
  10. Deepa Jayaprakashappa
  11. K VijayRaghavan
  12. Mani Ramaswami
  13. Baskar Bakthavachalu
Antagonistic roles for Ataxin-2 structured and disordered domains in RNP condensation
eLife 10:e60326.