The ORB2 RNA-binding protein represses translation of its target transcripts during the Drosophila maternal-to-zygotic transition via its functionally conserved Zinc-binding 'ZZ' domain

Abstract

RNA-binding proteins (RBPs) are key components of the post-transcriptional regulatory machinery. We show that the ORB2 RBP, the Drosophila ortholog of human Cytoplasmic Polyadenylation Element Binding Protein (hCPEB) 2-4 protein family, binds to hundreds of maternally provided mRNAs in early embryos, identify a U-rich motif enriched in ORB2's targets, and show that this motif confers ORB2 binding and repression to a luciferase reporter mRNA in S2 tissue culture cells. ORB2's target transcripts are translationally repressed and unstable during the maternal-to-zygotic transition (MZT), a developmental phase during which a large proportion of maternally provided mRNAs are repressed and cleared. We show that, when tethered to a luciferase reporter, ORB2 and hCPEB2 (but not ORB and hCPEB1) repress translation and that the C-terminal Zinc-binding ('ZZ') domain of ORB2 is necessary and sufficient for repression. ORB2 interacts with a suite of post-transcriptional regulators in early embryos; a subset of these interactions is lost upon deletion of the ZZ domain, notably with the Cup repressive complex. Analysis of the early embryo's translatome in the presence or absence of the endogenous ZZ domain, shows that ORB2's targets move onto polysomes upon ZZ domain deletion, indicating that this domain mediates translational repression of ORB2's targets during the MZT. Together, our results assign a function to the ZZ domain and support a significant role for ORB2 in post-transcriptional regulation of maternal mRNAs during the Drosophila MZT.

Keywords: Cup; Cytoplasmic Polyadenylation Element Binding Protein (CPEB); ORB2; RNA-binding protein (RBP); S2 cell; embryo; maternal-to-zygotic transition (MZT); post-transcriptional regulation; pre-initiation complex (PIC); translational regulation.

Figure 1.. RIP-seq identifies over 450 ORB2 targets in 0-3 hr embryos.

(a) Domain structure of ORB2B and location of the antigen used to screen for phage-displayed synthetic antibodies. (b) Scatterplot depicting the Log₂(Fold Enrichment)/Input Abundance (TPM) for all genes identified in the RIP-seq analysis; each point represents a single gene. Values represent the average enrichment from three independent biological replicates. ORB2 targets (red) are defined as having fold enrichment >1.5 in the anti-ORB2 RIP-seq compared to input and a ratio of ORB2/Control > 1.5, with an FDR adjusted P-value < 0.05; genes enriched in both ORB2 and control RIP-seq (green) are defined as having fold enrichment >1.5 in both anti-ORB2 and control RIP-seq, with a ratio of ORB2/Control < 1.5 and an FDR adjusted P-value < 0.05; Control enriched transcripts (blue) are defined as having fold enrichment <−1.5 in the control RIP-seq, fold enrichment < −1.5 in the anti-ORB2 RIP-seq, and an FDR adjusted P-value < 0.05. Dotted lines indicate ±1.5 fold-change thresholds. (c) Violin and boxplots depicting the distribution of transcript abundance as measured by average TPM from three biological replicates based on bound status as described in (A). Wilcoxon Rank Sum Test with Bonferroni multiple-test correction: ** = P < 10⁻¹⁰; *** = P < 10⁻¹⁵. (d) Validation of RIP-seq targets by RT-qPCR; bar graphs depict the log₂(fold enrichment) in an anti-ORB2 RIP-RT-qPCR compared to control RIP-RT-qPCR; P-values given are derived from Wilcoxon Signed Rank Test comparing ORB2 RIP-RT-qPCR to control RIP-RT-qPCR.

Figure 2.. ORB2-bound transcripts are enriched in translationally repressed, unstable transcripts during the MZT.

Violin and boxplots depicting the relationship between ORB2 association and translational status, change in poly(A) tail length, and change in transcript abundance. Unbound genes include all expressed genes that were not defined as ORB2 targets in the RIP-seq analysis. (a) Translation efficiency, calculated from ribosome footprinting performed on 0-1 hour embryo extracts (Eichhorn et al. 2016). (b) Translation efficiency, calculated from ribosome footprinting performed on 2-3 hour embryo extracts (Eichhorn et al. 2016). (c) Translation efficiency as measured by ribosome profiling performed on 0-2 hour embryo extracts (Dunn et al. 2013). (d) Translation index, calculated from polysome gradients performed on 0-3 hour embryo extracts (Chen et al. 2014). (e) Change in transcript poly(A) tail length, calculated by comparing 0-1 hour and 2-3 hour embryo extracts, as measured by PAL-Seq (Eichhorn et al. 2016). (f) Change in abundance, calculated by comparing 0-1 hour and 2-3 hour embryo extracts, as measured by RNA-Seq (Eichhorn et al. 2016). Wilcoxon Rank Sum Test÷ ** = P < 10⁻¹⁰; *** = P < 10⁻¹⁵.

Figure 3.. ORB2 binding motifs predicted from RIP-seq confer translational repression on a reporter mRNA in S2 cells.

(a) ORB2 recognition element (ORE) enriched in ORB2 targets identified by de novo motif discovery using #ATS (Li et al. 2010) on results from the anti-ORB2 RIP-seq in this study. (b) Diagram depicting the reporter Firefly Luciferase mRNA with six copies of different ORB2-binding motifs in the 3′UTR. (c) Bar graph depicting the enrichment of the luciferase reporter mRNA in an anti-ORB2 RIP-qPCR compared to Renilla Luciferase mRNA, normalized to a reporter without ORE motifs. Mutation of either the OREs or deletion of the RNA-binding domain (RRM) abrogate binding to the mRNA. (d) Bar graph depicting the ratio of luminescence between Firefly luciferase and Renilla luciferase. Mutation of either the OREs or deletion of the RNA-binding domain (RRM) abrogate repression. One-way ANOVA with Tukey multiple test correction: * = P < 0.05.

Figure 4.. The ORB2 and hCPEB2 ZZ domain but not the ORB and hCPEB1 ZZ domain repress a luciferase reporter.

(a) Diagram depicting the tethering assay in which a polypeptide of interest fused to the BIV Tat domain is tethered to a Firefly luciferase reporter mRNA carrying six tandem repeat TAR elements in the 3′UTR. (b) Bar graph depicting the expression of the reporter firefly luciferase (normalized to co-transfected Renilla luciferase and normalized to steady state RNA levels) when tethering whole ORB2 or ORB2 fragments compared to GFP control. The ORB2 ZZ domain is necessary and sufficient for repression. (c) Bar graph depicting the expression of the reporter firefly luciferase mRNA measured by RT-qPCR (normalized to co-transfected Renilla luciferase mRNA) when tethering whole ORB2 or ORB2 fragments compared to GFP control. Reporter mRNA levels are unaffected by ORB2 ZZ or ORB2ΔZZ. (d) amino acid sequence alignment of the ZZ domains of ORB2, ORB, and their human orthologs, hCPEB1 and hCPEB2. The arrowhead denotes the junction between the RBD and ZZ domains. (e) Bar graph depicting the expression of the reporter firefly luciferase (normalized to co-transfected Renilla luciferase) when tethering different CPEB ZZ domains. Only ORB2 and hPEB2 repress. (f) Bar graph depicting the expression of the reporter firefly luciferase when ORB ZZ chimeric domains harboring substitutions from the ORB2 ZZ domain, are tethered. The ORB2 ‘A’ and ‘C’ regions are each sufficient to convert the ORB ZZ domain into a repressor. One-way ANOVA with Hold-Sidak multiple test correction: * = P < 0.05; ** = P < 0.001; *** = P < 10⁻⁵.

Figure 5.. Anti-ORB2 IP-MS from early embryo extracts identifies known translation repressors and the PIC as binding partners.

(a) Interaction network depicting the connections between components of the Cup repressive complex and other translational repressors identified as ORB2 interactors by SAINT analysis. (b) Interaction network depicting the connections between components of eukaryotic PIC identified as ORB2 interactors by SAINT analysis. Interaction networks were produced by the GeneMANIA app for Cytoscape (Franz et al. 2018); size of node indicates average number of spectral counts obtained in the IP-MS; thickness of edges indicates strength of interaction assigned by GeneMANIA.

Figure 6.. Deletion of the ORB2 ZZ domain results in derepression of target RNAs in the early embryo.

(a) Violin and boxplots comparing target and non-target abundance as determined by RNA-Seq in 0-1, 1-2 and 2-3 hour old embryos laid by w¹¹¹⁸ and orb2^ΔZZ females. No significant difference was detected for either genotype at any of the time-points. (b) Violin and boxplots comparing target and non-target translation index as determined by polysome gradients followed by RNA-Seq in 0-1, 1-2 and 2-3 hour old embryos laid by w¹¹¹⁸ and orb2^ΔZZ females. ORB2 targets were significantly derepressed in 0-1 hour old orb2^ΔZZ embryos relative to w¹¹¹⁸. Log2(Ratio of mean normalized counts of pool 3 /mean normalized counts of pool 1 mean) was calculated for each transcript. Wilcoxon Rank Sum Test: ** = P < 0.006.