The YTHDF proteins ECT2 and ECT3 bind largely overlapping target sets and influence target mRNA abundance, not alternative polyadenylation
Abstract
Gene regulation via N6-methyladenosine (m6A) in mRNA involves RNA-binding proteins that recognize m6A via a YT521-B homology (YTH) domain. The plant YTH domain proteins ECT2 and ECT3 act genetically redundantly in stimulating cell proliferation during organogenesis, but several fundamental questions regarding their mode of action remain unclear. Here, we use HyperTRIBE (targets of RNA-binding proteins identified by editing) to show that most ECT2 and ECT3 targets overlap, with only a few examples of preferential targeting by either of the two proteins. HyperTRIBE in different mutant backgrounds also provides direct views of redundant, ectopic, and specific target interactions of the two proteins. We also show that contrary to conclusions of previous reports, ECT2 does not accumulate in the nucleus. Accordingly, inactivation of ECT2, ECT3, and their surrogate ECT4 does not change patterns of polyadenylation site choice in ECT2/3 target mRNAs, but does lead to lower steady-state accumulation of target mRNAs. In addition, mRNA and microRNA expression profiles show indications of stress response activation in ect2/ect3/ect4 mutants, likely via indirect effects. Thus, previous suggestions of control of alternative polyadenylation by ECT2 are not supported by evidence, and ECT2 and ECT3 act largely redundantly to regulate target mRNA, including its abundance, in the cytoplasm.
Introduction
N6-methyladenosine (m6A) in mRNA is of fundamental importance in eukaryotic gene regulation (Zhao et al., 2017). It is deposited in the nucleus (Salditt-Georgieff et al., 1976; Ke et al., 2017; Huang et al., 2019) by the METTL3/METTL14-methyltransferase complex (MTA/MTB in plants) (Bokar et al., 1997; Zhong et al., 2008; Liu et al., 2014) and many functions of m6A involve RNA-binding proteins that recognize m6A in mRNA via a YT521-B homology (YTH) domain (Stoilov et al., 2002; Dominissini et al., 2012; Wang et al., 2014; Zaccara et al., 2019). The YTH-domain family is subdivided into two phylogenetic groups, YTHDF and YTHDC (Patil et al., 2018; Balacco and Soller, 2019), but the biochemistry used for m6A recognition is identical in both groups: an aromatic cage provides a hydrophobic environment for the N6-methyl group and stacking interactions with the adenine ring (Li et al., 2014b; Luo and Tong, 2014; Theler et al., 2014; Wang et al., 2014; Xu et al., 2014; Zhu et al., 2014). YTHDF proteins consist of an N-terminal intrinsically disordered region (IDR) followed by the YTH domain (Patil et al., 2018). Early reports seemed to indicate functional specialization of vertebrate YTHDF proteins for either translational activation or mRNA decay (Wang et al., 2014; Wang et al., 2015; Li et al., 2017; Shi et al., 2017), whereas recent studies support functional redundancy among the three YTHDFs in mammals and zebrafish (Kontur et al., 2020; Lasman et al., 2020; Zaccara and Jaffrey, 2020), similar to the functional overlap described earlier for a subgroup of plant YTHDFs (Arribas-Hernández et al., 2018).
In plants, the YTHDF family is greatly expanded, with 11 members in Arabidopsis, referred to as EVOLUTIONARILY CONSERVED C-TERMINAL REGION1-11 (ECT1-11), compared to 3 in humans (Li et al., 2014a; Scutenaire et al., 2018). To date, only ECT2, ECT3, and ECT4 have been studied using genetic approaches. These studies show that the three m6A readers have important functions in post-embryonic development, but appear to work largely redundantly, at least in formal genetic terms. This is because single knockouts of ECT2 or ECT3 produce only subtle phenotypes related to branching patterns of epidermal hairs and root growth directionality, while simultaneous knockout of ECT2 and ECT3 results in delayed organogenesis and defective morphology of leaves, roots, stems, flowers, and fruits; defects that are exacerbated by additional mutation of ECT4 in most cases (Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020). It remains unclear, however, which mRNAs are targeted by ECT2/3, and what the effects of ECT2/3 binding to them may be (Arribas-Hernández and Brodersen, 2020). In particular, it is not clear whether the formal genetic redundancy between ECT2 and ECT3 is reflected in an overlapping target set, as would be expected for truly redundant action, or whether ECT2 and ECT3 might bind separate targets in wild type plants, but are able to replace each other in the artificial situation created by gene knockouts. The fact that knockouts of ECT2 and ECT3 have opposite effects on root growth directionality (Arribas-Hernández et al., 2020) indicates that at least some level of functional specialization exists between them, despite the obvious genetic redundancy observed in the control of organogenesis. Thus, it is an open question of fundamental importance for understanding plant m6A-YTHDF axes whether, and to what degree, mRNA targets of ECT2 and ECT3 overlap.
ECT2 has previously been suggested to act in the nucleus to influence alternative polyadenylation of targets (Wei et al., 2018). This model implies that plant ECT2 would act fundamentally differently from metazoan YTHDF proteins that are thought to be exclusively cytoplasmic and act to control mRNA fate via accelerated mRNA decay, or translational status (Patil et al., 2018; Zaccara et al., 2019; Kan et al., 2021; Worpenberg et al., 2021), perhaps in some cases by influencing the ability of other RNA binding proteins to associate with specific mRNAs (Worpenberg et al., 2021). The evidence for nuclear localization of ECT2 is not unequivocal, however, because the ECT2 signal presumed to be nuclear has not been examined relative to a nuclear envelope marker. In contrast, all studies examining the subcellular localization of ECT2 (and ECT3 and ECT4) have clearly established their presence in the cytoplasm (Arribas-Hernández et al., 2018; Scutenaire et al., 2018; Wei et al., 2018; Arribas-Hernández et al., 2020). In addition, the model of ECT2-mediated gene regulation via alternative polyadenylation has not been tested by direct experimentation.
Here, we use the proximity-labeling method HyperTRIBE (targets of RNA binding proteins identified by editing) (McMahon et al., 2016; Xu et al., 2018) to identify mRNA targets of ECT3. HyperTRIBE uses fusion of an RNA-binding protein to a hyperactive mutant (E488Q) of the catalytic domain of the Drosophila melanogaster adenosine deaminase acting on RNA (DmADARcd) to obtain an A-G mutation profile specifically in mRNAs bound by the RNA-binding protein of interest in vivo. We combine the comparative analysis of this data set with the target identification of ECT2 by HyperTRIBE and iCLIP (individual nucleotide resolution crosslinking and immunoprecipitation) (König et al., 2010) reported in the accompanying paper (Arribas-Hernández et al., 2021), a series of transcriptomic analyses in ect2/ect3/ect4 triple knockout mutants, and super-resolution microscopy of ECT2 localization, to establish three fundamental properties of mRNA regulation by ECT2 and ECT3. (1) Most targets are shared between ECT2 and ECT3, and the two proteins act genuinely redundantly in vivo to bind to and regulate many targets, in agreement with their similar expression patterns and genetically redundant functions (Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020). (2) ECT2/3/4 do not appreciably influence alternative polyadenylation of target mRNAs, consistent with the absence of ECT2-mCherry from the nucleoplasm. (3) In ECT2-expressing cell populations, the abundance of the majority of ECT2/3-target mRNAs is reduced upon loss of ECT2/3/4 activity.
Results
Identification of ECT3 target mRNAs using HyperTRIBE
To identify target mRNAs of ECT3 transcriptome-wide, we chose HyperTRIBE, because of our demonstration in the accompanying paper (Arribas-Hernández et al., 2021) that it efficiently identifies ECT2 targets with little expression bias. We therefore proceeded in exactly the same way as described for ECT2, in this case using transgenic lines expressing AtECT3pro:AtECT3-FLAG-DmADARE488Qcd-AtECT3ter (henceforth ‘ECT3-FLAG-ADAR’) in ect3-1 single mutants, and of FLAG-ADAR under the control of the ECT3 promoter (henceforth simply ‘FLAG-ADAR’) in wild type background as negative control (Figure 1A). Complementation of the developmental phenotype of triple ect2-1/ect3-1/ect4-2 (te234) mutants (Arribas-Hernández et al., 2018) by expression of ECT3-FLAG-ADAR at comparable levels was also verified (Figure 1A, Figure 1—figure supplement 1A). Five lines of each type were used for mRNA-seq of dissected shoot and root apices, and the data was analyzed to identify differentially edited sites (Figure 1B, Figure 1—figure supplements 1 and 2). Despite the lower expression of ECT3 compared to ECT2 (Arribas-Hernández et al., 2018) and, consequently, generally lower editing proportions in ECT3-FLAG-ADAR lines compared to ECT2-FLAG-ADAR lines (Arribas-Hernández et al., 2021; Figure 1C), the implementation of the HyperTRIBER pipeline (Rennie et al., 2021) to call significant editing sites successfully identified 2448 targets in aerial tissues, and 3493 in roots (ECT3 HT-targets) (Figure 1B, Supplementary file 1). As seen for ECT2 (Arribas-Hernández et al., 2021), the ECT3 target genes in shoot and root apices largely overlapped, and the editing proportions of individual editing sites showed a strong correlation (Figure 1D and E). Accordingly, most aerial- or root-specific targets could be explained by differences in expression between tissues (Figure 1F). The identification of strongly overlapping target sets in roots and shoots is expected from the similar roles of ECT3 in promoting growth and cell division in the two tissues (Arribas-Hernández et al., 2020) and, therefore, constitutes an argument for robustness of ECT3 target identification by the HyperTRIBE method.
ECT2 and ECT3 bind to overlapping sets of targets
We next analyzed the degree to which ECT2 and ECT3 HT-targets overlap. The data sets are directly comparable, as growth conditions, tissue dissection, RNA extraction, library construction, and sequencing depth for target identification of ECT3 by HyperTRIBE were identical to those used for ECT2 (Arribas-Hernández et al., 2021; Figure 2—figure supplement 1). Remarkably, more than 94% of ECT3 HT-targets overlapped with the larger group of ECT2 HT-targets in both aerial and root tissues, and there was a clear correlation between the editing proportions of the common editing sites, albeit with higher editing by ECT2-FLAG-ADAR overall (Figure 2A and B). Indeed, the pattern of editing sites resulting from fusion of ADAR to ECT2 or ECT3 was similar for many targets, with a few more sites typically detected in the ECT2-HT data set (e.g., ATPQ, Figure 2C, left panel). Nevertheless, we also noticed examples with preferential targeting by ECT2-FLAG-ADAR (e.g., TUA4, Figure 2C, middle panel) or, interestingly, by the less abundant ECT3-FLAG-ADAR (e.g., UBQ6, Figure 2C, right panel), perhaps hinting to molecular explanations for the recently described non-redundant roles of ECT2 and ECT3 in determining root growth directionality (Arribas-Hernández et al., 2020). Overall, however, the overwhelming overlap between ECT2 and ECT3 HT-targets in both tissues suggests that binding to the same mRNA targets underlies their genetically redundant functions in leaf and root formation (Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020). We also observed that ECT2 and ECT3 mRNAs contain m6A sites in seedlings according to published data sets (Shen et al., 2016; Parker et al., 2020), and that their protein products target their own and each other’s transcripts (Figure 2—figure supplement 2), indicating that autoregulatory feedback may contribute to control their expression.
HyperTRIBE provides direct views of redundant target mRNA Interactions with ECT2 and ECT3
Although the demonstration that ECT2 and ECT3 bind to strongly overlapping target sets is consistent with largely redundant in vivo function, it does not constitute a direct proof. For example, the proteins may bind to the same targets, but in different cells such they act de facto non-redundantly. We reasoned that HyperTRIBE might provide a means to observe directly whether ECT2 and ECT3 act specifically or redundantly on shared targets, and whether one ECT protein acquires non-natural targets upon knockout of the other by comparison of editing proportions measured with ADAR fusions expressed in single versus triple mutant backgrounds. The single mutant background (e.g., ECT2-FLAG-ADAR in ect2-1) would mimic the wild type setting, while the triple mutant background (e.g., ECT2-FLAG-ADAR in te234) would probe target interactions in the absence of redundant or competing proteins, but still in plants with wild type growth rates (Figure 1A; Arribas-Hernández et al., 2021). Redundant target interactions would be expected to result in generally higher editing proportions of the same target sites in triple mutant than in single mutant backgrounds, especially for the least expressed protein, ECT3 (Figure 3—figure supplement 1). Conversely, specific interactions would cause one of two possible signatures. (i) In the case of cell type-specific interactions, no change in editing proportions between single and triple mutant backgrounds should be detectable. (ii) In the case of specific interactions within the same cells in wild type, but non-natural targeting in the absence of other ECT proteins, acquisition of non-natural targets in triple mutant backgrounds is expected. We observed widespread increases in editing proportions for ECT3 targets upon removal of ECT2/4, while such increases occurred only sporadically for ECT2 targets (Figure 3A–C, Supplementary file 2), supporting the idea of largely redundant target interactions. Furthermore, although more sites showed increased editing by ECT3-FLAG-ADAR in aerial tissues than in roots in the absence of ECT2/4 (Figure 3B), the net increase of editing proportions was higher in roots for both ECT2 and ECT3 (Figure 3C), consistent with a more dominant role of ECT2 over ECT3 in aerial tissues compared to roots (Figure 2A and B). Importantly, the higher editing proportions in the triple mutant background cannot be trivially explained by higher expression of the transgene in these lines, as the average expression was comparable or slightly lower (Figure 3—figure supplement 1). These observations directly support genuinely redundant interactions of ECT2 and ECT3 with the majority of their mRNA targets in vivo.
Small sets of specific ECT2/3 targets acquire unnatural ECT interactions in knockout backgrounds
Although the tendency of ECT2/3 to show redundant target mRNA interactions was widespread, we also looked for examples of specific interactions in the HyperTRIBE data in single and triple mutant backgrounds. A priori, we considered targets to be ECT2-specific if they were edited by ECT2-FLAG-ADAR, but not ECT3-FLAG-ADAR, in single mutants (strictly specific), or became edited by ECT3-FLAG-ADAR only in the triple mutant. The definition of ECT3-specific targets followed analogous criteria. However, because ECT2 expression is much higher than ECT3 expression (Figure 3—figure supplement 1), ECT2-specific targets identified in this way may simply be below the detection limit of the less expressed ECT3-FLAG-ADAR transgene. Hence, arguments for existence of bona fide specific targets must take detectability by ECT3-FLAG-ADAR into account. Consistent with the expectation from the different ECT2/ECT3 dosage, much larger numbers of strictly ECT2-specific transcripts were identified compared to ECT3: 2414 ECT2-specific and 93 ECT3-specific targets were identified in aerial tissues, while in roots, 1738 were ECT2-specific and 197 were ECT3-specific (Figure 3D and E, Figure 3—figure supplement 2). In addition, small sets of specific target mRNAs became targets of the other ECT protein upon knockout of its genuine interacting protein (110 and 24 for ECT2-specific targets in aerial and root tissues, respectively, and 2 for ECT3-specific targets in roots) (Figure 3D and E, Figure 3—figure supplement 2). These sets constitute outstanding candidates for ECT2/3-specific mRNA targets. Interestingly, mRNAs of four tubulin subunits (TUBA1, TUBB1, TUBB5, and TUBB9) were among the many ECT2-specific targets in roots. This finding may be biologically relevant, because ect2 and ect3 single mutants exhibit distinct root slanting phenotypes (Arribas-Hernández et al., 2020), and misregulation of tubulin subunits causes root slanting defects (Smyth, 2016). Curiously, a few transcripts (21 in aerial tissues that include four photosynthesis-related genes [LHCB1.2, LHB1B2, LHB1B1, and AT3G63540], and 9 in roots) were edited by either ECT2 or ECT3 only in the triple mutant background (Figure 3D and E, Figure 3—figure supplement 2). Because ECT4 is also knocked out in te234, these sets define putative ECT4-specific targets. In summary, our comparative analyses of ECT2/3 HyperTRIBE data obtained in single and triple mutant backgrounds indicate that redundant target interaction is pervasive, but they also identify small target sets with properties consistent with preferential interaction with only one ECT protein.
ECT2/3 targets tend to be co-expressed in proliferating cells and are enriched in functions related to basic metabolism and protein synthesis
We next combined the ECT3-HT target set described here with the ECT2 iCLIP and ECT2-HT data (Arribas-Hernández et al., 2021) to define three gene sets of particular interest for functional analysis of ECT2/3: The permissive target set (6528 genes) defined as genes with either ECT2 HT, ECT3 HT, or ECT2 iCLIP support, the stringent target set (1992 genes) defined as all ECT2 or ECT3 HT-targets that are also in the ECT2 iCLIP target set, and the non-target set (13,504 genes) defined as all expressed genes not contained in the permissive target set (Figure 4A, Figure 4—figure supplement 1, Supplementary file 3). As an initial check of consistency of the target sets with the biological context in which ECT2 and ECT3 function, we used single-cell transcriptome analysis of Arabidopsis roots (Denyer et al., 2019; Ma et al., 2020) to analyze the overlap of ECT2/3 expression with the enrichment of markers for different cell types in the permissive target set (roots only, Figure 4—figure supplement 1). This analysis showed reassuring congruence between predominant expression of ECT2/3 in meristem clusters and marker enrichment for these same clusters among targets (Figure 4B, Figure 4—figure supplement 2). We also analyzed the permissive target sets for groups of functionally related genes, and found that ECT2/3 targets are enriched in housekeeping genes, many related to basic metabolism and protein synthesis (Figure 4C). This result is in agreement with the functional enrichment analysis of m6A-containing transcripts in mature leaves (Anderson et al., 2018) and of genes differentially expressed in rosettes of mutants partially depleted of m6A (Bodi et al., 2012; Anderson et al., 2018). Together with the recurrent function of ECT2/3/4 in organogenesis (Arribas-Hernández et al., 2020), these results suggest that regulation of metabolism and growth is the main biological function of the m6A-ECT2/3/4 regulatory axis. These initial analyses also provide well-defined common ECT2/3 target sets for further functional studies.
Recovery of ECT2-expressing cell populations with and without ECT2/ECT3/ECT4 activity
ECT2, ECT3, and ECT4 expression is largely restricted to rapidly dividing cells of organ primordia (Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020), and since many ECT2/3 targets are broadly expressed housekeeping genes (Figure 4C), cell populations expressing ECT2–4 need to be isolated prior to transcriptome analyses to avoid confounding effects from cells that do not express these m6A readers. We therefore used the fact that ect2-1/ECT2-mCherry exhibits root growth rates similar to wild type while te234/ECT2W464A-mCherry exhibits clearly reduced root growth rates nearly identical to te234 triple knockouts (Arribas-Hernández et al., 2020), and applied fluorescence-associated cell sorting to select mCherry-expressing cell populations from root protoplasts of three independent transgenic lines for each of these two genetic backgrounds (Figure 5A). Because wild type and mutant fluorescent proteins have the same expression level, pattern, and intracellular localization (Figure 5B and C), this procedure yielded comparable ECT2-expressing cell populations (Figure 5D, Figure 5—figure supplement 1) with (ECT2-mCherry/ect2-1 henceforth ‘wild type’) or without (ECT2W464A-mCherry/ect2-1/ect3-1/ect4-2, henceforth ‘mutant’) ECT2/3/4 function. We therefore isolated mRNA and constructed Smart-Seq2 libraries for comparison of poly(A) sites (PASs) and abundance of ECT2/3 targets and non-targets in ECT2-expressing cells from plants of the two different genetic backgrounds. Compared to standard mRNA-seq, Smart-seq2 recovers more reads with untemplated As (beginning of poly(A) tails) in addition to gene-specific sequence and can, therefore, be used for PAS mapping. We note that the selection of ECT2-expressing cells from the root meristem of wild type and mutant lines also circumvents the trouble of preparing comparable samples from intact tissues of plants at different developmental stages despite having the same age.
ECT2/3/4 do not play a direct role in alternative polyadenylation of targets
We first addressed the conjecture on a nuclear role of ECT2 in PAS selection (Wei et al., 2018). In plants, PASs are not sharply defined but rather spread along localized regions and can be grouped into PAS clusters (PACs) for analysis (Wu et al., 2011; Sherstnev et al., 2012). Using a modification of the nanoPARE analysis pipeline (Schon et al., 2018) to map PASs from reads with ≥9 untemplated As, we identified a total of 14,667 PACs belonging to 12,662 genes after filtering possible false positives (see Materials and methods; Figure 5—figure supplement 2A,B, Supplementary file 4). We found no tendency for ECT2/3 target mRNAs to have more PACs than non-targeted genes (Figure 5—figure supplement 2C), suggesting that differential PAC location in ECT2/3 targets between mutant and wild type is not prevalent. Nevertheless, we specifically tested whether PASs could be affected by the loss of ECT2/3/4 function in two different ways: either a shift of the dominant PAC to an alternative PAC altogether, or a shift in the most common PAS within clusters. Sorting the 206 genes for which the dominant PAC differed between wild type and mutant samples (18.5% of the 1114 genes with more than one PAC) into the ECT2/3 target groups in roots (Figure 4—figure supplement 1B, Supplementary file 4) showed that both the permissive and stringent targets were significantly less likely than non-targets to have a different dominant PAC upon loss of ECT2/3/4 function (p=0.013 and p=1.21e−5 for strictly permissive and stringent targets, respectively; Fisher’s exact test) (Figure 5E, Figure 5—figure supplement 2C,D). This significant difference may be an effect of the higher expression of targets compared to non-targets (Arribas-Hernández et al., 2021), as accuracy of PAS detection increases with transcript abundance (see Figure 5—figure supplement 2E for details). The result indicates that the alternative polyadenylation observed upon loss of ECT2/3/4 function is not prevalent among ECT2/3 targets. Finally, we examined changes to the local distribution of PASs within clusters. We defined the most common PAS as the single position in all overlapping PACs with the most reads, and determined the distance between such dominant PASs in wild type and mutant samples. Comparison of the distances revealed that the most common PAS does not change by more than 10 bp in the majority of genes, and is not more likely to be different in ECT2/3 targets than in non-targets (Figure 5F). In fact, the most common PAS is more likely to be unchanged in targets than in non-targets (Figure 5G) (p=0.028 and p=2.2e−16 for strictly permissive and stringent targets, respectively; Fisher’s exact test). Taken together, these analyses show that neither the usage of alternative PACs nor the dominant PASs within clusters have any tendency to change in ECT2/3 targets upon loss of ECT2/3/4 function.
ECT2-mCherry does not localize to the nucleoplasmic side of the nuclear envelope
To further investigate whether ECT2 may have any nuclear functions, we revisited the evidence for localization of ECT2 in the nucleoplasm, which is based on confocal fluorescence microscopy of ECT2-GFP or YFP-ECT2 in DAPI-stained root cells of stable Arabidopsis lines (Scutenaire et al., 2018; Wei et al., 2018). Because (i) the localization of ECT2-mCherry in living root cells of our lines has a general sharp boundary with what we interpreted to be the nucleus (Figure 5B; Arribas-Hernández et al., 2018) and does not overlap with nucleoplasmic MTA-TFP (Arribas-Hernández et al., 2020), (ii) paraformaldehyde fixation routinely used to permeate DAPI inside plant tissues (used by Scutenaire et al. and not specified by Wei et al.) can introduce artifacts in the localization of fluorescent proteins (Man-Wah et al., 2015), and (iii) the RNA-binding properties of DAPI could yield signal from the RNA-rich rough endoplasmic reticulum surrounding the nucleus (Tanious et al., 1992), we decided to examine the localization of ECT2-mCherry relative to the nuclear envelope in living cells. We therefore crossed lines expressing functional ECT2-mCherry (Arribas-Hernández et al., 2020) with plants expressing the outer nuclear envelope and nuclear pore complex-associated protein WIP1 fused to GFP (Xu et al., 2007). Confocal fluorescence microscopy of intact roots showed that the sharp boundaries of the ECT2-mCherry expression domain were delimited by the GFP-WIP1 signal from the nuclear envelope (Figure 6A, Figure 6—figure supplement 1A). Importantly, the occasional points at which the ECT2-mCherry signal seemed to fuzzily spill into the nucleus (white arrows in Figure 6A, Figure 6—figure supplement 1A) overlapped with equally blurry GFP-WIP1 signal, probably due to lack of perpendicularity between the nuclear envelope and the optical section in these areas. In such cases, the cytoplasm, nucleus and nuclear envelope may be contained in the same region of the optical section and thus appear to be overlapping (Figure 6A, bottom panel). To verify this interpretation, we inspected our plants with the super-resolution confocal Airyscan detector (Huff, 2015) and, as expected, we did not observe ECT2-mCherry signal inside the GFP-WIP1-delimited nuclei in any instances (Figure 6B, Figure 6—figure supplement 1B,C). Based on these analyses, we conclude that ECT2 resides in the cytoplasm and its presence in the nucleus, if any, may be too transient to be detected by fluorescence microscopy. These results agree with the lack of evidence for a function of ECT2/3/4 in choice of PAS, and strongly suggest that the molecular basis for the importance of ECT2/3/4 should be sought in cytoplasmic properties of their mRNA targets.
ECT2/3 targets tend to show reduced abundance upon loss of ECT2/3/4
We next assessed the effect of loss of ECT2/3/4 function on target mRNA abundance, using the Smart-seq2 data from FACS-sorted root protoplasts described above. Principal component analysis showed that the three repeats of ‘wild type’ (ect2-1/ECT2-mCherry) were well separated from the three repeats of ‘mutant’ (te234/ECT2W464A-mCherry) along the first principal component (Figure 7—figure supplement 1A), indicating that the differential gene expression analysis was meaningful. We focused on stringent, permissive, and non-ECT2/3 targets in roots (Figure 4—figure supplement 1B, Supplementary file 3), and visualized their differential expression by scatter, volcano, and box plots (Figure 7A–C, Supplementary file 5). These approaches showed that stringent targets have a clear tendency toward downregulation upon loss of ECT2/3/4 function. This trend is maintained, but is less pronounced in permissive targets, and is reversed in non-targets (Figure 7A–C). Indeed, of the significantly differentially expressed stringent ECT2/3 targets, nearly all were downregulated in the mutant, while the majority of differentially expressed non-targets were upregulated compared to wild type (Figure 7D). Furthermore, ECT2/3 targets accounted for more than half of all significantly downregulated genes, but only about 15% of upregulated genes (Figure 7E). In contrast, highly upregulated genes tended to be non-targets (Figure 7B, right panel).
Functional groups of differentially expressed genes
To test if these differentially regulated gene sets represented subsets of functionally related genes within target and non-target groups, we analyzed their potential enrichment of gene ontology terms. This analysis revealed that downregulated ECT2/3 targets were particularly enriched in genes related to ribosome biogenesis and translation (Figure 7F), while upregulated non-targets were enriched in ‘abiotic stress responses’ with molecular function ‘transcription factor’ (Figure 7G). Because cell wall digestion required for protoplast isolation is a cellular stress, we tested the trivial possibility that loss of ECT2/3/4 function renders cells more susceptible to stress, and that such potential hyper-susceptibility underlies the observed differences of gene expression in ECT2-expressing root protoplasts. To this end, we isolated RNA from intact root apices of 4-day-old plants of Col-0 wild type and te234 mutants, and performed mRNA-seq analysis. These results recapitulated the trends of downregulation of stringent ECT2/3 targets and upregulation of stress-responsive non-targets, albeit with less pronounced differences than observed in the selected ECT2-expressing cell populations as expected (Figure 7—figure supplement 1). We also noticed that several stress-inducing and growth-restricting NF-YA-class transcription factors, all repressed by miR169defg (Song et al., 2019), were upregulated in root tips (Figure 8A), and used small RNA-seq to test if activation of the stress response was visible in the miRNA expression profile. Indeed, the miR169defg family was specifically repressed in te234 mutants, and the miR167 family, targeting growth-promoting auxin-response factors, was clearly upregulated (Figure 8B). In addition, the LTR-retrotransposon-targeting miR845a (Borges et al., 2018) was strongly upregulated. Thus, the stress response detected in te234 mutants comprises coherent changes of miRNA and transcription factor expression. Because it has been recently reported that processing of some miRNAs is affected by MTA-dependent methylation of their precursors in Arabidopsis (Bhat et al., 2020), we queried our iCLIP and HyperTRIBE data for evidence of direct binding of ECT2/3 to pri-miRNAs, possibly explaining the differential abundance of the miR169, miR167, and miR845 families in te234 mutants. However, we could not find evidence for ECT2/3 binding to these or any other pri-miRNA in our data, consistent with the cytoplasmic localization of ECT2 and the predominantly, though not exclusively, nuclear functions of pri-miRNA (Fang and Spector, 2007; Lauressergues et al., 2015). Thus, the differential abundance of mature miRNAs in roots of te234 seedlings compared to wild type is likely to be caused by the same indirect effects that lead to induction of stress-related mRNAs.
Altogether, the analysis of data from intact root tips confirms that the observed patterns of differential gene expression in selected protoplasts are genuine and biologically meaningful, and that the selection of ECT2-expressing cells ensures the most accurate description of differential gene expression resulting from loss of ECT2/3/4 function. We note that while the differential gene expression analysis suggests that ECT2/3/4 formally act to increase abundance of their mRNA targets, it does not allow conclusions to be drawn on how such activation is brought about: a direct stabilizing effect of ECT2/3/4 binding to their targets is consistent with the observed results, but indirect effects via transcriptional repression cannot be excluded, especially given the presence of stress-related transcription factors in the set of upregulated non-targets.
Discussion
Our identification and comparative analyses of mRNA targets of ECT2 (Arribas-Hernández et al., 2021) and ECT3, and the study of target behavior in terms of abundance and use of alternative polyadenylation in cell populations devoid of ECT2/3/4 activity, allow us to draw conclusions on (i) the redundancy, and (ii) the subcellular localization and functional impact of these YTHDF proteins in Arabidopsis.
Redundancy within the YTHDF protein family in plants
Combining the overlapping expression patterns of ECT2 and ECT3, their formal genetic redundancy in growth promotion (Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020), their overlapping target sets, and the signatures of redundant target interaction derived from HyperTRIBE in single and triple ect mutant backgrounds, we conclude that many target mRNAs can bind to either ECT2 or ECT3 with similar consequences. Thus, ECT2 and ECT3 can exhibit redundant function sensu stricto, not just the ability to replace function in the absence of the other protein. However, this does not mean that the ECT2 and ECT3 genes fulfill identical function under all conditions. For example, a recent study of the yeast stress-related transcription factors MSN2 and MSN4 shows that distributing fully redundant biochemical function over two genes expressed in the same cell can be beneficial, because strong induction by stress can be combined with steady low-noise expression under favorable conditions (Chapal et al., 2019); a combination difficult, if not impossible, to obtain with control of expression by a single promoter.
Looking at the complete ECT family beyond ECT2 and ECT3, our previous studies may offer additional insight into redundancy with ECT4. Because of the enhancement of developmental phenotypes of ect2/ect3 plants by mutation of ECT4, and the overlapping expression pattern of the three proteins (Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020), ECT4 may act redundantly with ECT2/3 on the same targets. The reduced impact of absence of ECT4 compared to that of ECT2/3 may simply be due to lower expression levels, but could also involve lower binding affinity, reduced activity, or a combination of those properties. However, such conclusions on potentially redundant function cannot be drawn for any of the eight additional members of the YTHDF family in Arabidopsis. Indeed, RNA-seq data suggests that at least some YTHDFs are expressed in different tissues. For example, ECT10 mRNA is highly abundant in pollen from which the rest of the family is absent (Arribas-Hernández et al., 2018). Although it is possible that specific expression patterns or induction by environmental cues are the main reasons for the expansion of the YTHDF family in plants (Li et al., 2014a; Scutenaire et al., 2018), functional specialization that may include different target-binding specificity (Arribas-Hernández et al., 2021) and effector function conferred by the different IDRs is also possible. For those ECTs with overlapping expression patterns, competition for targets leading to different molecular outcomes is an interesting possibility. Further genetic and biochemical studies on additional ECTs are necessary to understand redundant and specific functions of YTHDF proteins, a question that is currently under debate also in animal systems (Kontur et al., 2020; Lasman et al., 2020; Zaccara and Jaffrey, 2020).
Subcellular localization and functional impact on mRNA targets of plant YTHDF proteins
Our data show that ECT2 is not nuclear, and that ECT2/3/4 do not appreciably affect alternative polyadenylation in their direct mRNA targets. It is important to note that this conclusion on ECT2/3/4 does not extend to m6A altogether. Polyadenylation and transcription termination are clearly influenced by m6A in plants, as shown by studies of mutants in core N6-adenosine methyltransferase components VIR and FIP37. In vir-1 mutants, a tendency to use proximal alternative PASs in m6A-targets is prominent (Parker et al., 2020), and in fip37 mutants, several cases of defective transcription termination causing the production of chimeric transcripts were noted (Pontier et al., 2019). In the latter case, m6A recognition by the YTHDC-containing nuclear subunit of the Cleavage and Polyadenylation Specificity Factor CPSF30 was demonstrated to be required for regular transcription termination (Pontier et al., 2019), and subsequent studies of CPSF30, including of mutants specifically defective in m6A-binding, also indicated its role in mediating m6A-dependent alternative polyadenylation (Hou et al., 2021; Song et al., 2021). Thus, m6A does affect 3′-end formation of methylated pre-mRNAs in plants, but this process involves the nuclear CPSF30 which contains a YTHDC domain only in plants, rather than the cytoplasmic YTHDF protein ECT2.
The immediate implication of the conclusions that ECT2 does not act in the nucleus and does not influence alternative polyadenylation is that the conceptual framework for m6A-YTHDF action established mostly through studies in mammalian cell culture does extend to plants: m6A is installed by a nuclear methyltransferase complex (Zhong et al., 2008; Shen et al., 2016; Ružička et al., 2017), probably coupled to RNA Polymerase II transcription (Bhat et al., 2020), and YTHDF-mediated regulation of transcripts carrying the m6A mark takes place in the cytoplasm. Then, what is the molecular effect of ECT2/ECT3-binding to target mRNAs? Clearly, lack of ECT2/3/4 results in decreased abundance of direct targets. This result is in line with the lower accumulation of m6A-containing transcripts observed in plants partially depleted of m6A (Shen et al., 2016; Anderson et al., 2018; Parker et al., 2020) and in the ect2-1 knockout line (Wei et al., 2018) compared to wild type plants. However, differential expression analysis of ubiquitously expressed targets performed with RNA from entire plants is not easy to interpret when the regulatory process subject to study is located only in well-defined cell populations of meristematic tissues (Zhong et al., 2008; Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020). The problem is even more acute when the tissue composition of wild type and mutant individuals to be compared is different due to the developmental delay of the m6A-deficient lines. Finally, some of the previous studies of differential gene expression also revealed a number of m6A-containing transcripts with increased abundance in m6A-deficient mutants (Shen et al., 2016; Anderson et al., 2018), of which some were proposed to be of importance for the observed phenotype (Shen et al., 2016). Thus, the studies on differential gene expression in m6A-deficient mutants compared to wild type published thus far do not allow clear conclusions on the consequence of loss of m6A (or a reader protein) for target mRNA accumulation to be drawn (Arribas-Hernández and Brodersen, 2020). In contrast, our experimental setup, using a tissue with comparable cell-type composition in wild type and mutant plants and extracting RNA only from the cells in which target regulation by ECT2/3/4 takes place, establishes that targets generally have decreased abundance upon loss of ECT2/3/4.
The present study does not, however, elucidate the mechanisms involved in target regulation, because it does not directly measure target mRNA synthesis and degradation rates. It is possible that the reduced target mRNA accumulation in ect2/3/4 mutant cells is exclusively a direct consequence of ECT2/3/4 function at the post-transcriptional level. For example, ECT2/3/4 could cause mRNA stabilization by protection from endonucleolysis as previously suggested (Anderson et al., 2018). It is also possible that ECT2/3/4 act primarily to repress target mRNA translation, as suggested by the recent demonstration that although Drosophila Ythdf translationally represses the futsch target mRNA, futsch mRNA abundance is reduced in ythdf mutants (Worpenberg et al., 2021). We cannot at present exclude, however, that more indirect effects also play a role, perhaps related to transcriptional repression of ECT2/3/4 targets via stress responses activated upon loss of ECT2/3/4 function. We anticipate that clear answers to this question must await development of tools for conditional inactivation of ECT2/3 function, such that consequences for mRNA target stability and translatability can be studied immediately after loss of ECT2/3 binding. Finally, we note that the constitutive stress response activation is consistent with the stunted phenotype of mta and te234 mutants (Bodi et al., 2012; Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020), and that the earliest report of differentially expressed genes in an m6A-deficient organism, performed in plant leaves, also showed enrichment of stress-responsive genes, interpreted by the authors as possibly, ‘the consequence of a perceived stress due to reduced methylase activity’ (Bodi et al., 2012). Thus, disentangling the direct effects of m6A-ECT2/3/4 on growth via mRNA target regulation from possible indirect effects arising from stress response activation in knockout mutants will be of major importance in future studies.
Materials and methods
All data analyses were performed using TAIR10 and Araport11 as the reference genome and transcriptome, respectively. Unless otherwise stated, data analyses were performed in R (https://www.R-project.org/) and plots generated using either base R, IGV (for genomic data) (Robinson et al., 2011), or ggplot2 (https://ggplot2.tidyverse.org).
Definitions of experiment, biological replicates, and technical replicates
Request a detailed protocolWe use the terms ‘biological replicate’ (the only type of replicate present in this study), and ‘experiment’ in the same way as in the accompanying manuscript (Arribas-Hernández et al., 2021).
Plant material
Request a detailed protocolAll lines employed in this study are in the Arabidopsis thaliana Col-0 ecotype. The mutant alleles or their combinations: ect2-1 (SALK_002225) (Arribas-Hernández et al., 2018; Scutenaire et al., 2018; Wei et al., 2018), ect3-1 (SALK_077502), ect4-2 (GK_241H02), and ect2-1/ect3-1/ect4-2 (te234) (Arribas-Hernández et al., 2018) have been previously described. The transgenic lines ect2-1 ECT2-FLAG-ADAR and te234 ECT2-FLAG-ADAR (Arribas-Hernández et al., 2021), GFP:WIP1 (Xu et al., 2007) and those expressing ECT2pro:ECT2-mCherry-ECT2ter and ECT2pro:ECT2W464A-mCherry-ECT2ter (Arribas-Hernández et al., 2018; Arribas-Hernández et al., 2020) have also been described. Plants co-expressing ECT2-mCherry and GFP-WIP1 used for fluorescence microscopy were the F1 progeny of a genetic cross between GFP-WIP1 and ECT2-mCherry-expressing plants.
Growth conditions
Request a detailed protocolSeeds for HyperTRIBE, protoplast-FACS sorting and differential expression analyses of intact root tips were sterilized, stratified, and germinated, and seedlings were grown, in identical conditions than those described for HyperTRIBE and CLIP experiments in the accompanying manuscript (Arribas-Hernández et al., 2021). Similarly, we also used the growth conditions described in the corresponding section of the same article for phenotypic characterization of plants expressing ECT3-FLAG-ADAR.
Generation of transgenic lines for ECT3-HyperTRIBE
Request a detailed protocolWe generated lines expressing ECT3pro:ECT3-FLAG-DmADARE488Qcd-ECT3ter and ECT3pro:FLAG-DmADARE488Qcd-ECT3ter by USER cloning (Bitinaite and Nichols, 2009) and agrobacterium-mediated transformation in the same way as for the ECT2 equivalents (Arribas-Hernández et al., 2021). Primer sequences and their combination to obtain cloning fragments are detailed in Appendix 1. As for ECT2-HT, we selected five independent lines of each type based on segregation studies of sensitivity to glufosinate-ammonium (to isolate single T-DNA insertions), phenotypic complementation (in the te234 background), and transgene expression levels.
HyperTRIBE
Request a detailed protocolThe HyperTRIBE experiments were performed once, using five biological replicates (independent lines) for each of the groups (genotypes) used. Growth conditions and experimental procedures were identical for all the groups compared in this study. Root and aerial tissue were dissected from the same plants in all cases. Tissue dissection from of 10-day-old T2 seedlings, RNA extraction, and library preparation (by Novogene) were done as described for ECT2-HyperTRIBE (Arribas-Hernández et al., 2021).
Analysis of HyperTRIBE data
Request a detailed protocolSignificant differentially edited sites between ECT3-FLAG-ADAR (fusion) and FLAG-ADAR (control) samples were called according to our HyperTRIBER pipeline (Rennie et al., 2021) as described for ECT2-HyperTRIBE (Arribas-Hernández et al., 2021), without removal of any sample. Specific scripts for this analysis can be found at https://github.com/sarah-ku/targets_arabidopsis; Rennie, 2021.
For the analysis of editing sites by ECT2/3-FLAG-ADAR in triple (te234) versus single (either ect2-1 or ect3-1) mutant background, the HyperTRIBER pipeline (Rennie et al., 2021) was run between the two types of samples without taking into account the free ADAR controls, in order to detect positions edited preferentially in one or the other background. To account for low power as a result of high variance in editing proportions due to transgene expression differences across samples, scaled ADAR abundance was treated as an extra covariate in the model. This resulted in enriched sensitivity to specifically call A-to-G positions, which were subsequently considered as significant if they had an adjusted p-value<0.1 and an absolute log2(fold change)>0.25. We also required positions to be a significantly edited site in at least one of the single or triple mutant set ups against the free FLAG-ADAR control samples.
Comparison with root single-cell data
Request a detailed protocolThe expression matrix based on a total of 4727 individual cells from scRNA-seq in roots was downloaded from Denyer et al., 2019, together with extensive lists of marker genes associated with 15 clusters annotated to cell types in roots. To calculate the proportion of markers at target genes: for each of the 15 clusters, the proportion of marker genes that are ECT2 or ECT3 targets (based on ECT2-HT and ECT3-HT, respectively, in roots) was calculated. Proportions were then overlayed onto a t-SNE diagram (Denyer et al., 2019), according to relevant clusters of cells. Scripts for this analysis can be found at https://github.com/sarah-ku/targets_arabidopsis.
Preparation and sorting of protoplasts
Request a detailed protocolWe harvested roots from 5-day-old T4 seedlings grown on vertical square plates (20 plates with 4 rows of densely spotted seeds in each plate per line/replicate) to digest in 20 ml of protoplasting solution (20 mM MES, 0.4 M D-Mannitol, 20 mM KCl, 1.25% w/v Cellulase, 0.3 % Macerozyme, 0.1% w/v BSA, 10 mM CaCl2, and 5 mM β-mercaptoethanol; pH 5.7), following Benfey’s lab procedure (Birnbaum et al., 2005; Bargmann and Birnbaum, 2010). After 75 min of incubation at 27°C with gentle agitation, we filtered the cell-suspensions through a 40 μm strainer and pelleted cells by centrifugation at 500×g for 10 min at room temperature in a swinging-bucket centrifuge. Pellets were gently resuspended in 400 μl of protoplasting solution for direct sorting in a FACSAriaIII cytometer. The flow stream was adjusted to 20 psi sheath pressure with a 100 μm nozzle aperture. Sorted cells were collected into RTL buffer (QIAGEN) supplemented with 40 mM DTT (3.5 vol of buffer per volume of cell suspension) and lysed by vortexing. The protoplast extracts were flash-frozen on dry ice until extraction with the RNeasy Plus Micro Kit (QIAGEN) following the manufacturer’s instructions. The yield was ~300,000 cells in a volume of 1.5 ml (per sample). Samples were harvested, prepared, and sorted with a 15 min lapse between them to account for sorting time. In that way, every sample was processed in the same amount of time (~2 h from the start of harvesting to sorting). To prevent any possible bias, the samples of each genotype (3+3) were alternated during all the processing.
Smart-seq2
Request a detailed protocolSmart-seq2 libraries were generated according to Picelli et al., 2013 using the Illumina DNA Nextera Flex Kit (now called Illumina DNA Prep) from total RNA extracted with the RNeasy Plus Micro Kit (QIAGEN) from FACS-sorted root protoplasts (Birnbaum et al., 2005). The libraries were sequenced in PE75 mode on an Illumina NextSeq550 sequencer. Nextera transposase adapters were trimmed from all reads using cutadapt (Kechin et al., 2017).
Polyadenylation site analysis
Request a detailed protocolSmart-seq2 reads with at least nine 3′-terminal A nucleotides or 5′-terminal T nucleotides were labeled as putative poly(A)-containing reads and the oligo-A/T sequences were removed with a maximum allowed mismatch rate of 6%. All putative poly(A)-containing reads with a length >20 and a mean quality score >25 after trimming were retained along with their mate pair and mapped to the A. thaliana TAIR10 genome using STAR (Dobin et al., 2013) with the following parameters:
--alignIntronMax 5000 --alignMatesGapMax 5500 --outFilterMatchNmin 20
--alignSJDBoverhangMin 1 --outFilterMismatchNmax 5
--outFilterMismatchNoverLmax .05 --outFilterType BySJout
--outFilterIntronMotifs RemoveNoncanonicalUnannotated
Putative poly(A)-containing reads that mapped to the genome were filtered for false positives by examining the adjacent nucleotides in the genome: reads were removed if the putative poly(A) site was immediately upstream of a 15 nt region that is at least 80% purines (which are likely sites of oligo-dT mispriming). All putative poly(A)-containing reads not filtered in this way were retained as poly(A) sites and were counted for each position in the genome based on the most 3′ nucleotide of each read (allowing 3′-terminal mismatches).
Polyadenylation site clusters (PACs) were identified using a modification of the nanoPARE analysis pipeline (Schon et al., 2018) (specific scripts for this analysis are available at https://github.com/maschon0/ect_polyA_analysis; Schon, 2021). Briefly, reads from the samples above that did not contain untemplated poly(A) tails were mapped to the genome and used as a negative control of ‘gene body reads.’ Then, subtractive kernel density estimation was performed for each sample using endGraph.sh with default settings to produce a BED file for each sample containing poly(A) site clusters. As a final filter against oligo-dT mispriming events, the reads removed as false positives in the previous step that overlap with each cluster were counted. If the cluster contained more filtered signal than unfiltered signal, the entire cluster was considered a false positive cluster. Clusters were retained if an overlapping site was identified in at least two of the three replicates of both ect2-1/ECT2-mCherry and te234/ECT2W464A-mCherry genotypes. These two sets of clusters were merged using bedtools merge (Quinlan and Hall, 2010). Clusters mapping to the mitochondrial and chloroplast genome and the two rDNA loci were discarded, and the rest were retained for quantification with Salmon (Patro et al., 2017).
mRNA-seq and small RNA-seq from root tips
Request a detailed protocolTotal RNA purified from manually dissected root tips of 4-day-old plants (using the same growing conditions and methodology as for the HyperTRIBE lines) was used for the preparation of Illumina mRNA-seq (same methodology as for HyperTRIBE) and small RNA-seq libraries (NEBNext small RNA library prep set for Illumina). The experiment was performed once, using three biological replicates for mRNA-seq, and two for small RNA-seq.
Differential expression analysis (mRNA and miRNAs)
Request a detailed protocolDifferential gene expression analysis of mRNA was performed from processed (cutadapt Kechin et al., 2017), mapped (STAR, Dobin et al., 2013) and quantified (Salmon, Patro et al., 2017) Smart-seq2 or mRNA-seq data using DESeq2 (Love et al., 2014), for all genes with at least 1 TPM in all six samples (three biological replicates of the two types) and a total sum of at least 5 TPM. Significantly differentially expressed genes (FDR<0.05) were considered to be upregulated in the mutants if the fold change between mutant and wild type samples was higher than 1.5, or downregulated if lower than 1/1.5.
For small RNA-seq, raw reads were trimmed with cutadapt (Kechin et al., 2017) to lengths of 18–28 nt, and mapped to the Arabidopsis genome using STAR (Dobin et al., 2013) with genome indexes built on the Araport11_GTF_genes_transposons.Mar202021.gtf annotation. Mapped reads were counted using featureCounts (Liao et al., 2014). Genes with less than 1 RPM in all four samples (two biological replicates of the two types) were excluded from the analysis. Differential expression analysis was conducted using DESeq2 (Love et al., 2014) on the resulting matrix. Genes with adjusted p-value (FDR) lower or equal to 0.05 were considered significantly differentially expressed, and upregulated in te234 mutants if the fold change between te234 and Col-0 WT samples was higher than 1.5, or downregulated if lower than 1/1.5.
Gene ontology enrichment analysis
Request a detailed protocolThe functional enrichment analysis was carried out using the R package gprofiler2 (Raudvere et al., 2019).
Fluorescence microscopy
Request a detailed protocolEntire root tips growing inside MS-agar plates were imaged with a Leica MZ16 F stereomicroscope equipped with a Sony α6000 camera. Standard confocal fluorescence microscopy images of cells in root meristems were acquired with a Zeiss LSM700 confocal microscope as described in Arribas-Hernández et al., 2018 using ~7-day-old seedlings grown on MS-agar plates and freshly mounted in water. For super-resolution fluorescence microscopy, we used a Zeiss LSM900 equipped with the Airyscan detector (Huff, 2015). Fluorescence intensity plots were obtained with the tool ‘Plot Profile’ of the image-processing package ImageJ (Schindelin et al., 2012).
Appendix 1
Cloning | |
---|---|
Constructs | Primer pairs (fragments) |
ECT3pro:ECT3-FLAG-ADAR-ECT3ter (in pCAMBIA3300-U) | LA617-692 (ECT3pro:ECT3), LA693-642 (FLAG-ADAR), andLA626-622 (ECT3ter) |
ECT3pro:FLAG-ADAR-ECT3ter (in pCAMBIA3300-U) | LA617-627 (ECT3pro), LA694-642 (FLAG-ADAR),and LA626-622 (ECT3ter) |
Primers for USER-cloning | |
LA617.U-ECT3P.F | GGCTTAAUAGGCTTGGTTAGCAGAAGG |
LA622.ECT3T-U.R | GGTTTAAUTGGTTAACTCTATGGACTCATC |
LA626.ADAR/ECT3T.F | AACCAACCUTTGGTTTTAAGTGGGAAC |
LA627.ECT3P/FLAG.R | ATGGCTGCGUAGTGAGTGGCTTAG |
LA642.dADAR/ECT3T.R | AGGTTGGTUCATTCGGCAAGACCGAACTC |
LA692.ECT3/FLAG.R | AGGCAGUAGCAACAGCATTTTTCTCC |
LA693.ECT3/FLAG.F | ACTGCCUGCGATTACAAGGATGACGATGAC |
LA694.ECT3P/FLAG.F | ACGCAGCCAUGGATTACAAGGATGACGATG |
Data availability
Accession numbers: The raw and processed data for ECT3-HyperTRIBE, Smart-seq2 from root protoplasts and RNA-seq from root tips have been deposited in the European Nucleotide Archive (ENA) at EMBL-EBI under the accession number PRJEB44359.
Code availability: The scripts for the HyperTRIBE and scRNA-seq analyses can be found at
https://github.com/sarah-ku/targets_arabidopsis (copy archived at https://archive.softwareheritage.org/swh:1:rev:ad524fd57073b569320998bb79ecc66d433b37c7), and for the PAS analysis at https://github.com/maschon0/ect_polyA_analysis (copy archived at https://archive.softwareheritage.org/swh:1:rev:f469a02ac12ed83f442c9936b46c4f94336c10c1).
-
European Nucleotide ArchiveID PRJEB44359. Principles of mRNA targeting and regulation via Arabidopsis YTHDF proteins.
References
-
Fluorescence activated cell sorting of plant protoplastsJournal of Visualized Experiments 18:1673.https://doi.org/10.3791/1673
-
DNA cloning and engineering by uracil excisionCurrent Protocols in Molecular Biology 3:86.https://doi.org/10.1002/0471142727.mb0321s86
-
Purification and cDNA cloning of the AdoMet-binding subunit of the human mRNA (N6-adenosine)-methyltransferaseRNA 3:1233–1247.
-
Resolving noise-control conflict by gene duplicationPLOS Biology 17:e3000289.https://doi.org/10.1371/journal.pbio.3000289
-
Star: Ultrafast Universal RNA-SEQ alignerBioinformatics 29:15–21.https://doi.org/10.1093/bioinformatics/bts635
-
Cutprimers: A new tool for accurate cutting of primers from reads of targeted next generation sequencingJournal of Computational Biology 24:1138–1143.https://doi.org/10.1089/cmb.2017.0096
-
iCLIP reveals the function of HNRNP particles in splicing at individual nucleotide resolutionNature Structural & Molecular Biology 17:909–915.https://doi.org/10.1038/nsmb.1838
-
Context-dependent functional compensation between Ythdf M6A reader proteinsGenes & Development 34:1373–1391.https://doi.org/10.1101/gad.340695.120
-
Cytoplasmic m6A reader YTHDF3 promotes mRNA translationCell Research 27:444–447.https://doi.org/10.1038/cr.2017.10
-
A METTL3-METTL14 complex mediates mammalian nuclear RNA N6-adenosine methylationNature Chemical Biology 10:93–95.https://doi.org/10.1038/nchembio.1432
-
PscB: A Browser to Explore Plant Single Cell RNA-Sequencing Data SetsPlant Physiology 183:464–467.https://doi.org/10.1104/pp.20.00250
-
Reading m6A in the Transcriptome: m6A-Binding ProteinsTrends in Cell Biology 28:113–127.https://doi.org/10.1016/j.tcb.2017.10.001
-
Smart-seq2 for sensitive full-length transcriptome profiling in single cellsNature Methods 10:1096–1098.https://doi.org/10.1038/nmeth.2639
-
g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 updateNucleic Acids Research 47:W191–W198.https://doi.org/10.1093/nar/gkz369
-
Fiji: An open-source platform for biological-image analysisNature Methods 9:676–682.https://doi.org/10.1038/nmeth.2019
-
NanoPARE: parallel analysis of RNA 5’ ends from low-input RNAGenome Research 28:1931–1942.https://doi.org/10.1101/gr.239202.118
-
Direct sequencing of Arabidopsis thaliana RNA reveals patterns of cleavage and polyadenylationNature Structural & Molecular Biology 19:845–852.https://doi.org/10.1038/nsmb.2345
-
Helical growth in plant organs: Mechanisms and significanceDevelopment 143:3272–3282.https://doi.org/10.1242/dev.134064
-
MicroRNAs and Their Regulatory Roles in Plant–Environment InteractionsAnnual Review of Plant Biology 70:489–525.https://doi.org/10.1146/annurev-arplant-050718-100334
-
YTH: A new domain in nuclear proteinsTrends in Biochemical Sciences 27:495–497.https://doi.org/10.1016/s0968-0004(02)02189-8
-
Solution structure of the YTH domain in complex with N6-methyladenosine RNA: A reader of methylated RNANucleic Acids Research 42:13911–13919.https://doi.org/10.1093/nar/gku1116
-
Structural basis for selective binding of m6A RNA by the YTHDC1 YTH domainNature Chemical Biology 10:927–929.https://doi.org/10.1038/nchembio.1654
-
Reading, writing and erasing mRNA methylationNature Reviews. Molecular Cell Biology 20:608–624.https://doi.org/10.1038/s41580-019-0168-5
-
Post-transcriptional gene regulation by mRNA modificationsNature Reviews. Molecular Cell Biology 18:31–42.https://doi.org/10.1038/nrm.2016.132
Article and author information
Author details
Funding
H2020 European Research Council (PATHORISC ERC-2016-COG 726417)
- Peter Brodersen
Independent Research Fund Denmark (9040-00409B)
- Peter Brodersen
H2020 European Research Council (638173)
- Robin Andersson
Independent Research Fund Denmark (6108-00038B)
- Robin Andersson
H2020 European Research Council (63788)
- Michael D Nodine
The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.
Acknowledgements
The authors thank Lena Bjørn Johansson and Phillip Andersen for technical assistance in the construction of transgenic lines, Emilie Oksbjerg for preparing and sequencing small RNA libraries, Theo Bølsterli, René Hvidberg Petersen and their teams for plant care, and Anna Fossum and Rajesh Somasundaram for their assistance with FACS. Kim Rewitz and Kenneth Halberg are thanked for their assistance with Airyscan microscopy. The authors acknowledge Mathias Tankmar, Alexander JH Andersen, and Freja Asmussen for experimental support, and Tom Denyer and Marja Timmermans for their input and support in the analyses of ECT2/3 expression and target enrichment in their scRNAseq data. The authors are grateful to Norman R Groves and Iris Meier for the kind donation of GFP-WIP1 seeds.
Copyright
© 2021, Arribas-Hernández et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,806
- views
-
- 380
- downloads
-
- 41
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
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)
Further reading
-
- Chromosomes and Gene Expression
- Genetics and Genomics
Annotation of newly sequenced genomes frequently includes genes, but rarely covers important non-coding genomic features such as the cis-regulatory modules—e.g., enhancers and silencers—that regulate gene expression. Here, we begin to remedy this situation by developing a workflow for rapid initial annotation of insect regulatory sequences, and provide a searchable database resource with enhancer predictions for 33 genomes. Using our previously developed SCRMshaw computational enhancer prediction method, we predict over 2.8 million regulatory sequences along with the tissues where they are expected to be active, in a set of insect species ranging over 360 million years of evolution. Extensive analysis and validation of the data provides several lines of evidence suggesting that we achieve a high true-positive rate for enhancer prediction. One, we show that our predictions target specific loci, rather than random genomic locations. Two, we predict enhancers in orthologous loci across a diverged set of species to a significantly higher degree than random expectation would allow. Three, we demonstrate that our predictions are highly enriched for regions of accessible chromatin. Four, we achieve a validation rate in excess of 70% using in vivo reporter gene assays. As we continue to annotate both new tissues and new species, our regulatory annotation resource will provide a rich source of data for the research community and will have utility for both small-scale (single gene, single species) and large-scale (many genes, many species) studies of gene regulation. In particular, the ability to search for functionally related regulatory elements in orthologous loci should greatly facilitate studies of enhancer evolution even among distantly related species.
-
- Epidemiology and Global Health
- Genetics and Genomics
Alzheimer’s disease (AD) is a complex degenerative disease of the central nervous system, and elucidating its pathogenesis remains challenging. In this study, we used the inverse-variance weighted (IVW) model as the major analysis method to perform hypothesis-free Mendelian randomization (MR) analysis on the data from MRC IEU OpenGWAS (18,097 exposure traits and 16 AD outcome traits), and conducted sensitivity analysis with six models, to assess the robustness of the IVW results, to identify various classes of risk or protective factors for AD, early-onset AD, and late-onset AD. We generated 400,274 data entries in total, among which the major analysis method of the IVW model consists of 73,129 records with 4840 exposure traits, which fall into 10 categories: Disease, Medical laboratory science, Imaging, Anthropometric, Treatment, Molecular trait, Gut microbiota, Past history, Family history, and Lifestyle trait. More importantly, a freely accessed online platform called MRAD (https://gwasmrad.com/mrad/) has been developed using the Shiny package with MR analysis results. Additionally, novel potential AD therapeutic targets (CD33, TBCA, VPS29, GNAI3, PSME1) are identified, among which CD33 was positively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. TBCA and VPS29 were negatively associated with the main outcome traits of AD, as well as with both EOAD and LOAD. GNAI3 and PSME1 were negatively associated with the main outcome traits of AD, as well as with LOAD, but had no significant causal association with EOAD. The findings of our research advance our understanding of the etiology of AD.