CoChat
Sign Up
Cited by 256RNA Research and SplicingRNA and protein synthesis mechanismsRNA regulation and disease

The DOI matches its Crossref record; 83 of the 83 references it cites were found in a registry or an index.

GC content shapes mRNA storage and decay in human cells

Open access
Courel, Maïté, Clément, Yves, Bossevain, Clémentine, Foretek, Dominika, Vidal Cruchez, Olivia, Yi, Zhou, Bénard, Marianne, Benassy, Marie-Noëlle, Kress, Michel, Vindry, Caroline, Ernoult-Lange, Michèle, Antoniewski, Christophe, Morillon, Antonin, Brest, Patrick, Hubstenberger, Arnaud, Roest Crollius, Hugues, Standart, Nancy, Weil, Dominique·2019·eLife
DOI 10.7554/elife.49708PMC PMC6944446

Published version of a preprint

Open this paper in CoChat — free account

Abstract

mRNA translation and decay appear often intimately linked although the rules of this interplay are poorly understood. In this study, we combined our recent P-body transcriptome with transcriptomes obtained following silencing of broadly acting mRNA decay and repression factors, and with available CLIP and related data. This revealed the central role of GC content in mRNA fate, in terms of P-body localization, mRNA translation and mRNA stability: P-bodies contain mostly AU-rich mRNAs, which have a particular codon usage associated with a low protein yield; AU-rich and GC-rich transcripts tend to follow distinct decay pathways; and the targets of sequence-specific RBPs and miRNAs are also biased in terms of GC content. Altogether, these results suggest an integrated view of post-transcriptional control in human cells where most translation regulation is dedicated to inefficiently translated AU-rich mRNAs, whereas control at the level of 5’ decay applies to optimally translated GC-rich mRNAs.

Overview

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Integrity signals

Registry

Crossref · 10.7554/elife.49708

Cites

83 of 83 references verified0 citation errors · 0 not checkable

Statistics

Figures

Figure 1. PB mRNAs are AU-rich and longer than average. (A) Long mRNAs are particularly enriched in PBs. Transcripts were subdivided into six classes depending on their length (from <1.5 kb to >10 kb). The boxplots represent the distribution of their respective enrichment in PBs. The boxes represent the 25–75 percentiles and the whiskers the 10–90 percentiles. rs, Spearman correlation coefficient. (B) AU-rich mRNAs are particularly enriched in PBs. Transcripts were subdivided into six classes depending on their GC content (from <40 to >60%) and analyzed as in (A). (C) PBs mostly contain the AU-rich fraction of the transcriptome. The human transcriptome was binned depending on its GC content (0.7% GC increments). The graph represents the number of PB-enriched (PB-in, p<0.05, n = 5200) and PB-excluded (PB-out, p<0.05, n = 4669) transcripts in each bin. The distribution of all transcripts is shown for comparison (n = 14443). The median GC value is indicated below for each group. (D) mRNA localization in PBs mostly depends on the GC content of their CDS and 3’UTR. The analysis was repeated as in (B) using the GC content of the 5’UTR, CDS or 3’UTR, as indicated. For 5’UTRs, the >60% class was subdivided into three classes to take into account their higher GC content compared to CDSs and 3’UTRs. −0.57 and −0.55 are not significantly different (p=0.17), while −0.22 and −0.55 are (p<0.0001) (E) GC content is lower in PB-enriched mRNAs than PB-excluded ones independently of their length. The GC content distribution of PB-enriched (PB-in, p<0.05) and PB-excluded (PB-out, p<0.05) mRNAs was analyzed as in (B).
Figure 1. PB mRNAs are AU-rich and longer than average. (A) Long mRNAs are particularly enriched in PBs. Transcripts were subdivided into six classes depending on their length (from <1.5 kb to >10 kb). The boxplots represent the distribution of their respective enrichment in PBs. The boxes represent the 25–75 percentiles and the whiskers the 10–90 percentiles. rs, Spearman correlation coefficient. (B) AU-rich mRNAs are particularly enriched in PBs. Transcripts were subdivided into six classes depending on their GC content (from <40 to >60%) and analyzed as in (A). (C) PBs mostly contain the AU-rich fraction of the transcriptome. The human transcriptome was binned depending on its GC content (0.7% GC increments). The graph represents the number of PB-enriched (PB-in, p<0.05, n = 5200) and PB-excluded (PB-out, p<0.05, n = 4669) transcripts in each bin. The distribution of all transcripts is shown for comparison (n = 14443). The median GC value is indicated below for each group. (D) mRNA localization in PBs mostly depends on the GC content of their CDS and 3’UTR. The analysis was repeated as in (B) using the GC content of the 5’UTR, CDS or 3’UTR, as indicated. For 5’UTRs, the >60% class was subdivided into three classes to take into account their higher GC content compared to CDSs and 3’UTRs. −0.57 and −0.55 are not significantly different (p=0.17), while −0.22 and −0.55 are (p<0.0001) (E) GC content is lower in PB-enriched mRNAs than PB-excluded ones independently of their length. The GC content distribution of PB-enriched (PB-in, p<0.05) and PB-excluded (PB-out, p<0.05) mRNAs was analyzed as in (B).CC BY · PubMed Central
Figure 1—figure supplement 1. PB-enriched mRNAs tend to be long and AU-rich. (A) The transcripts enriched in PBs tend to be long. The mRNA length was expressed as a function of PB enrichment. rs, Spearman correlation coefficient.(B) Number of transcripts in the different classes analyzed in Figure 1A,E. (C) PB mRNAs are longer than PB-excluded ones, but shorter than SG mRNAs. The boxplot represents the distribution of length of PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs, and SG-enriched (SG-in) and SG-depleted (SG-out) mRNAs, as taken from Wells et al. (2017). The length of PB-enriched and SG-enriched mRNAs was statistically different using a two-tail Mann-Whitney test. **, p<0.0001. (D) The transcripts enriched in PBs tend to be AU-rich. The mRNA GC content was expressed as a function of PB enrichment. rs, Spearman correlation coefficient. (E) Number of transcripts in the different classes analyzed in Figure 1B,D and Figure 1—figure supplement 1. (F) PBs mostly contain transcripts from AU-rich genes. The graph was constructed as in Figure 1C using the GC content of the genes. The median value is indicated below for each group. (G) Partial correlations between PB localization (PBs), mRNA GC content (GC.mRNA), mRNA expression (RPKM) and sequence conservation (dN/dS). For each line, individual cells represent the correlation between the two corresponding variables, controlling for the effects of all other variables. dN/dS is the ratio of non-synonymous (dN) to synonymous (dS) changes between two organisms, here computed between human and mouse. It is used as a measure of natural selection strength in coding sequence. (H) PB-enriched mRNAs are longer than PB-excluded ones independently of their GC content. The distribution of length of PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs was analyzed as in Figure 1E.
Figure 1—figure supplement 1. PB-enriched mRNAs tend to be long and AU-rich. (A) The transcripts enriched in PBs tend to be long. The mRNA length was expressed as a function of PB enrichment. rs, Spearman correlation coefficient.(B) Number of transcripts in the different classes analyzed in Figure 1A,E. (C) PB mRNAs are longer than PB-excluded ones, but shorter than SG mRNAs. The boxplot represents the distribution of length of PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs, and SG-enriched (SG-in) and SG-depleted (SG-out) mRNAs, as taken from Wells et al. (2017). The length of PB-enriched and SG-enriched mRNAs was statistically different using a two-tail Mann-Whitney test. **, p<0.0001. (D) The transcripts enriched in PBs tend to be AU-rich. The mRNA GC content was expressed as a function of PB enrichment. rs, Spearman correlation coefficient. (E) Number of transcripts in the different classes analyzed in Figure 1B,D and Figure 1—figure supplement 1. (F) PBs mostly contain transcripts from AU-rich genes. The graph was constructed as in Figure 1C using the GC content of the genes. The median value is indicated below for each group. (G) Partial correlations between PB localization (PBs), mRNA GC content (GC.mRNA), mRNA expression (RPKM) and sequence conservation (dN/dS). For each line, individual cells represent the correlation between the two corresponding variables, controlling for the effects of all other variables. dN/dS is the ratio of non-synonymous (dN) to synonymous (dS) changes between two organisms, here computed between human and mouse. It is used as a measure of natural selection strength in coding sequence. (H) PB-enriched mRNAs are longer than PB-excluded ones independently of their GC content. The distribution of length of PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs was analyzed as in Figure 1E.CC BY · PubMed Central
Figure 2. Codon usage is strongly biased in PBs. (A) PB mRNAs and PB-excluded mRNAs encode proteins with different amino acid usage. The graph represents the frequency of each amino acid in the proteins encoded by mRNAs enriched or excluded from PBs, using the indicated PB enrichment thresholds. (B) Codon usage bias in and out of PBs follows their GC content. The relative codon usage for each amino acid was calculated in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs, using a PB enrichment threshold of +/- 1 (in log2). The graph represents the log2 of their ratio (PB-in/PB-out) and was ranked by decreasing values for each amino acid. The GC content of each codon is gray-coded below, using the scale indicated on the right. (C) The usage of some codons is biased independently of their GC content. Two examples are shown encoding Leucine (L) and Glycine (G). (D) The frequency of low usage codons strongly correlates with the GC content of the CDS, independently of their PB localization. The frequency of low usage codons was calculated for mRNAs excluded (PB-out) and enriched (PB-in) in PBs using a PB enrichment threshold of +/- 1 (in log2). It was expressed as a function of the CDS GC content at position 3 (GC3). Note that the slopes of the tendency curves are similar for PB-enriched and PB-excluded transcripts. The difference between the Spearman correlation coefficients (rs) are nevertheless statistically significant (p<0.0001). (E) PB mRNAs have longer CDS than PB-excluded mRNAs. The analysis was performed as in Figure 1E. (F) The number of low usage codons per CDS is a good determinant of both protein yield and PB localization. The protein yield was expressed as a function of the number of low usage codons for PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. rs, Spearman correlation coefficient.
Figure 2. Codon usage is strongly biased in PBs. (A) PB mRNAs and PB-excluded mRNAs encode proteins with different amino acid usage. The graph represents the frequency of each amino acid in the proteins encoded by mRNAs enriched or excluded from PBs, using the indicated PB enrichment thresholds. (B) Codon usage bias in and out of PBs follows their GC content. The relative codon usage for each amino acid was calculated in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs, using a PB enrichment threshold of +/- 1 (in log2). The graph represents the log2 of their ratio (PB-in/PB-out) and was ranked by decreasing values for each amino acid. The GC content of each codon is gray-coded below, using the scale indicated on the right. (C) The usage of some codons is biased independently of their GC content. Two examples are shown encoding Leucine (L) and Glycine (G). (D) The frequency of low usage codons strongly correlates with the GC content of the CDS, independently of their PB localization. The frequency of low usage codons was calculated for mRNAs excluded (PB-out) and enriched (PB-in) in PBs using a PB enrichment threshold of +/- 1 (in log2). It was expressed as a function of the CDS GC content at position 3 (GC3). Note that the slopes of the tendency curves are similar for PB-enriched and PB-excluded transcripts. The difference between the Spearman correlation coefficients (rs) are nevertheless statistically significant (p<0.0001). (E) PB mRNAs have longer CDS than PB-excluded mRNAs. The analysis was performed as in Figure 1E. (F) The number of low usage codons per CDS is a good determinant of both protein yield and PB localization. The protein yield was expressed as a function of the number of low usage codons for PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. rs, Spearman correlation coefficient.CC BY · PubMed Central
Figure 2—figure supplement 1. Amino acid usage and codon usage biases in PBs. (A) PB mRNAs and PB-excluded mRNAs encode proteins with different amino acid usage. The boxplots represent the distribution of the frequency of each amino acid in the proteins encoded by mRNAs enriched (PB-in, n = 5500) or excluded from PBs (PB-out, n = 5104). The boxes represent the 25–75 percentiles and the whiskers the 10–90 percentiles. (B,C) Codon usage biases independent on the GC content. For each amino acid with 4 or six synonymous codons, the ratio of usage of NNU to NNA (B) or NNC to NNG (C) was calculated. The boxplots represent the distribution of these ratios in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs using a PB enrichment threshold of +/- 1(log2). All differences were statistically significant (*, p<0.05) except for the Serine NNT/NNA ratio. Note the change of scale between (B) and (C). (D) Low usage codons are more frequent in PB mRNAs. The ratio between the relative codon usage in PB-enriched and PB-excluded (PB-in/PB-out) mRNAs was expressed as a function of their global relative usage, after normalization by the number of synonymous codons. (E) The frequency of low usage codons correlates weakly with protein yield. The protein yield was expressed as a function of the frequency of low usage codons for PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. rs, Spearman correlation coefficient.
Figure 2—figure supplement 1. Amino acid usage and codon usage biases in PBs. (A) PB mRNAs and PB-excluded mRNAs encode proteins with different amino acid usage. The boxplots represent the distribution of the frequency of each amino acid in the proteins encoded by mRNAs enriched (PB-in, n = 5500) or excluded from PBs (PB-out, n = 5104). The boxes represent the 25–75 percentiles and the whiskers the 10–90 percentiles. (B,C) Codon usage biases independent on the GC content. For each amino acid with 4 or six synonymous codons, the ratio of usage of NNU to NNA (B) or NNC to NNG (C) was calculated. The boxplots represent the distribution of these ratios in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs using a PB enrichment threshold of +/- 1(log2). All differences were statistically significant (*, p<0.05) except for the Serine NNT/NNA ratio. Note the change of scale between (B) and (C). (D) Low usage codons are more frequent in PB mRNAs. The ratio between the relative codon usage in PB-enriched and PB-excluded (PB-in/PB-out) mRNAs was expressed as a function of their global relative usage, after normalization by the number of synonymous codons. (E) The frequency of low usage codons correlates weakly with protein yield. The protein yield was expressed as a function of the frequency of low usage codons for PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. rs, Spearman correlation coefficient.CC BY · PubMed Central
Figure 2—figure supplement 2. Codon usage biases and abundance of amino-acylated tRNA. Codon usage in PB-excluded mRNAs does not closely follow amino-acylated tRNA abundance. The abundance of charged tRNA isodecoders in HEK293 cells was taken from Evans et al. (24) and aggregated for each codon (gray bars). Codons with no cognate tRNA are indicated with a star, and their possible decoding by near-cognate tRNAs with arrows. Normalized codon usage in PB-enriched and PB-excluded mRNAs is on the right y axis.
Figure 2—figure supplement 2. Codon usage biases and abundance of amino-acylated tRNA. Codon usage in PB-excluded mRNAs does not closely follow amino-acylated tRNA abundance. The abundance of charged tRNA isodecoders in HEK293 cells was taken from Evans et al. (24) and aggregated for each codon (gray bars). Codons with no cognate tRNA are indicated with a star, and their possible decoding by near-cognate tRNAs with arrows. Normalized codon usage in PB-enriched and PB-excluded mRNAs is on the right y axis.CC BY · PubMed Central
Figure 3. DDX6 silencing has opposite effects on mRNA fate depending on their GC content. (A) mRNA stabilization after DDX6 silencing in HEK293 and K562 cells applies to GC-rich mRNAs. The fold-changes (FC) in mRNA accumulation were analyzed as in Figure 1B. (B) mRNA translation derepression after DDX6 silencing in HEK293 cells applies to AU-rich mRNAs. The fold-changes in translation rate (polysomal/total mRNA ratio) were analyzed as in (A). (C) GC-rich mRNAs are particularly enriched in the DDX6 CLIP experiment.
Figure 3. DDX6 silencing has opposite effects on mRNA fate depending on their GC content. (A) mRNA stabilization after DDX6 silencing in HEK293 and K562 cells applies to GC-rich mRNAs. The fold-changes (FC) in mRNA accumulation were analyzed as in Figure 1B. (B) mRNA translation derepression after DDX6 silencing in HEK293 cells applies to AU-rich mRNAs. The fold-changes in translation rate (polysomal/total mRNA ratio) were analyzed as in (A). (C) GC-rich mRNAs are particularly enriched in the DDX6 CLIP experiment.CC BY · PubMed Central
Figure 3—figure supplement 1. Polysome profiling following DDX6 silencing. (A) Western blot analysis of DDX6 in siRNA-transfected HEK293 cells. The DDX6 signal, normalized using the ribosomal protein S6, is indicated on the right panel (n = 3; p<0.001; t-test). Error bars, SD. (B) Representative polysome profiles of siRNA-transfected cells. (C,D) MA plots of mRNA fold changes between siRNA-transfected and control cells in total (C) and polysomal (D) mRNA. (E) Correlation between fold changes (FC) of randomly selected mRNAs measured by qPCR and calculated from RNA-Seq data, in total and polysomal fractions (n = 3). Error bars, SD. (F,G) Correlation between total (F) and polysomal (G) mRNA fold-changes and mRNA GC content. Note that correlation is positive for total mRNAs and negative for polysomal mRNAs. rs, Spearman correlation coefficient.
Figure 3—figure supplement 1. Polysome profiling following DDX6 silencing. (A) Western blot analysis of DDX6 in siRNA-transfected HEK293 cells. The DDX6 signal, normalized using the ribosomal protein S6, is indicated on the right panel (n = 3; p<0.001; t-test). Error bars, SD. (B) Representative polysome profiles of siRNA-transfected cells. (C,D) MA plots of mRNA fold changes between siRNA-transfected and control cells in total (C) and polysomal (D) mRNA. (E) Correlation between fold changes (FC) of randomly selected mRNAs measured by qPCR and calculated from RNA-Seq data, in total and polysomal fractions (n = 3). Error bars, SD. (F,G) Correlation between total (F) and polysomal (G) mRNA fold-changes and mRNA GC content. Note that correlation is positive for total mRNAs and negative for polysomal mRNAs. rs, Spearman correlation coefficient.CC BY · PubMed Central
Figure 3—figure supplement 2. Impact of DDX6 binding and mRNA length on DDX6 dependency. (A) Number of transcripts in the different classes analyzed in Figures 3 and 4A,B. (B) Dot plot representation of the changes in translation rate (Polysome/Total FC) as a function of changes in mRNA stability (Total FC) after DDX6 silencing. (C,D) DDX6-clipped targets tend to be particularly stabilized after DDX6 silencing in both HEK293 and K562 cells (C), but less affected than other mRNAs in terms of translation rate (D). The difference between DDX6-clipped mRNAs and all mRNAs was statistically different using a two-tail Mann-Whitney test. **, p<0.0001. (E) Number of transcripts in the different classes analyzed in (F) and (G). (F) Transcript length has a weak impact on mRNA abundance following DDX6 silencing. Transcripts were subdivided into six classes depending on their full length (from <1.5 to 10 kb, left panel), the length of their CDS (from <0.5 to 6 kb, middle panel) or 3’UTR (from <0.5 to 6 kb, right panel). The boxplots represent the distribution of their respective fold-changes in total mRNA. The difference between class 1 and 6 was statistically significant using a two tail Mann-Whitney test: p<0.0001 in all panels. rs, Spearman correlation coefficient. (G) The 3’UTR length has an impact on mRNA translation derepression following DDX6 silencing. The polysomal/total mRNA ratio was used as a proxy measurement of translation rate.The data are represented as in (F). Note that the difference between class 1 and 6 was statistically significant in left and right panels only.
Figure 3—figure supplement 2. Impact of DDX6 binding and mRNA length on DDX6 dependency. (A) Number of transcripts in the different classes analyzed in Figures 3 and 4A,B. (B) Dot plot representation of the changes in translation rate (Polysome/Total FC) as a function of changes in mRNA stability (Total FC) after DDX6 silencing. (C,D) DDX6-clipped targets tend to be particularly stabilized after DDX6 silencing in both HEK293 and K562 cells (C), but less affected than other mRNAs in terms of translation rate (D). The difference between DDX6-clipped mRNAs and all mRNAs was statistically different using a two-tail Mann-Whitney test. **, p<0.0001. (E) Number of transcripts in the different classes analyzed in (F) and (G). (F) Transcript length has a weak impact on mRNA abundance following DDX6 silencing. Transcripts were subdivided into six classes depending on their full length (from <1.5 to 10 kb, left panel), the length of their CDS (from <0.5 to 6 kb, middle panel) or 3’UTR (from <0.5 to 6 kb, right panel). The boxplots represent the distribution of their respective fold-changes in total mRNA. The difference between class 1 and 6 was statistically significant using a two tail Mann-Whitney test: p<0.0001 in all panels. rs, Spearman correlation coefficient. (G) The 3’UTR length has an impact on mRNA translation derepression following DDX6 silencing. The polysomal/total mRNA ratio was used as a proxy measurement of translation rate.The data are represented as in (F). Note that the difference between class 1 and 6 was statistically significant in left and right panels only.CC BY · PubMed Central
Figure 3—figure supplement 3. Impact of the GC content on DDX6-dependency. (A) Number of transcripts in the different classes analyzed in (B) and (C). (B,C) mRNA abundance (B) and translation derepression (C) following DDX6 silencing in HEK293 cells mostly depends on the GC content of their CDS and 3’ UTR. The analysis was performed as in Figure 1D, using the polysomal/total mRNA ratio as a proxy measurement of translation rate. In (B), the Spearman correlation coefficients (rs) 0.40, 0.42 and 0.18 are significantly different from each other (p<0.0001). In (C), −0.52 and −0.15 are significantly different (p<0.0001), while 0.53 and 0.52 are not (p=0.37).
Figure 3—figure supplement 3. Impact of the GC content on DDX6-dependency. (A) Number of transcripts in the different classes analyzed in (B) and (C). (B,C) mRNA abundance (B) and translation derepression (C) following DDX6 silencing in HEK293 cells mostly depends on the GC content of their CDS and 3’ UTR. The analysis was performed as in Figure 1D, using the polysomal/total mRNA ratio as a proxy measurement of translation rate. In (B), the Spearman correlation coefficients (rs) 0.40, 0.42 and 0.18 are significantly different from each other (p<0.0001). In (C), −0.52 and −0.15 are significantly different (p<0.0001), while 0.53 and 0.52 are not (p=0.37).CC BY · PubMed Central
Figure 4. XRN1 and PAT1B targets have distinct GC content. (A) mRNA stabilization after XRN1 silencing in HeLa and HCT116 cells applies to GC-rich mRNAs. The analysis was performed as in Figure 1B. The GC content distribution for all mRNAs is presented for comparison (in gray). (B) mRNA stabilization after PAT1B silencing in HEK293 cells applies to AU-rich mRNAs. The analysis was performed as in (A). (C) Read coverage of PAT1B targets (FC >0.7, n = 330, solid lines) and all mRNAs (n = 16000, dashed lines) in the siPAT1B dataset. The read coverage was analyzed in each duplicate experiment and normalized as described in Materials and methods. The average value in control cells (gray lines) and after PAT1B silencing (peach lines) was plotted, with the bars representing the duplicate values. An expanded view of the dashed box is presented on the right panel. (D) Read coverage of XRN1 targets (FC >0.8, n = 199, solid lines) and all mRNAs (n = 13760, dashed lines) in the siXRN1 dataset. The data were analyzed as in (C). (E) Clustering analysis of mRNAs depending on their GC content, their differential expression after silencing DDX6, XRN1 or PAT1B, and their enrichment in PBs. Raw GC content and log2 transformed ratio of the other datasets were used for the clustering of both transcripts (lines) and datasets (columns). The values were color-coded as indicated on the right scale, and the Spearman correlation matrix is presented below (all p<10−48). The heatmap highlights the distinct fate of GC-rich and AU-rich mRNAs.
Figure 4. XRN1 and PAT1B targets have distinct GC content. (A) mRNA stabilization after XRN1 silencing in HeLa and HCT116 cells applies to GC-rich mRNAs. The analysis was performed as in Figure 1B. The GC content distribution for all mRNAs is presented for comparison (in gray). (B) mRNA stabilization after PAT1B silencing in HEK293 cells applies to AU-rich mRNAs. The analysis was performed as in (A). (C) Read coverage of PAT1B targets (FC >0.7, n = 330, solid lines) and all mRNAs (n = 16000, dashed lines) in the siPAT1B dataset. The read coverage was analyzed in each duplicate experiment and normalized as described in Materials and methods. The average value in control cells (gray lines) and after PAT1B silencing (peach lines) was plotted, with the bars representing the duplicate values. An expanded view of the dashed box is presented on the right panel. (D) Read coverage of XRN1 targets (FC >0.8, n = 199, solid lines) and all mRNAs (n = 13760, dashed lines) in the siXRN1 dataset. The data were analyzed as in (C). (E) Clustering analysis of mRNAs depending on their GC content, their differential expression after silencing DDX6, XRN1 or PAT1B, and their enrichment in PBs. Raw GC content and log2 transformed ratio of the other datasets were used for the clustering of both transcripts (lines) and datasets (columns). The values were color-coded as indicated on the right scale, and the Spearman correlation matrix is presented below (all p<10−48). The heatmap highlights the distinct fate of GC-rich and AU-rich mRNAs.CC BY · PubMed Central
Figure 4—figure supplement 1. Transcriptome analysis following XRN1 and PAT1B silencing. (A) Western blot analysis of XRN1 in siRNA-transfected HeLa cells. The XRN1 signal, normalized using tubulin, is indicated below. (B) Western blot analysis of XRN1 in shRNA-transfected HCT116 cells. Cells transfected with the XRN1 or control shRNA were induced for 0 to 72 hr with Doxycycline (DOX). The XRN1 signal, normalized using RNA Pol II, is indicated below. The time point 48 hr was chosen for library preparation. (C) XRN1 targets tend to be excluded from PBs. The fold-changes after XRN1 silencing in HeLa cells were expressed as a function of PB enrichment. (D, E) PAT1B targets are regulated by DDX6 at the level of translation but not decay. The total (D) and polysomal/total (E) mRNA fold-changes (FC) after DDX6 silencing were expressed as a function of fold-changes after PAT1B silencing. (F) PAT1B targets are enriched in PBs. The fold-changes of total RNA after PAT1B silencing were expressed as a function of PB enrichment. rs, Spearman correlation coefficients.
Figure 4—figure supplement 1. Transcriptome analysis following XRN1 and PAT1B silencing. (A) Western blot analysis of XRN1 in siRNA-transfected HeLa cells. The XRN1 signal, normalized using tubulin, is indicated below. (B) Western blot analysis of XRN1 in shRNA-transfected HCT116 cells. Cells transfected with the XRN1 or control shRNA were induced for 0 to 72 hr with Doxycycline (DOX). The XRN1 signal, normalized using RNA Pol II, is indicated below. The time point 48 hr was chosen for library preparation. (C) XRN1 targets tend to be excluded from PBs. The fold-changes after XRN1 silencing in HeLa cells were expressed as a function of PB enrichment. (D, E) PAT1B targets are regulated by DDX6 at the level of translation but not decay. The total (D) and polysomal/total (E) mRNA fold-changes (FC) after DDX6 silencing were expressed as a function of fold-changes after PAT1B silencing. (F) PAT1B targets are enriched in PBs. The fold-changes of total RNA after PAT1B silencing were expressed as a function of PB enrichment. rs, Spearman correlation coefficients.CC BY · PubMed Central
Figure 5. GC biases in the targets of various RNA decay factors, translation regulators and miRNAs. (A) GC content biases in the targets of various RBPs. The targets of the indicated factors were defined using CLIP experiments or motif analysis (see Materials and methods). The boxplots represent the distribution of the GC content of their gene. The distribution for all mRNAs is presented for comparison (in gray) and the red dashed line indicates its median value. (B) Heatmap representation of the different factors depending on the behavior of their mRNA targets in the different datasets. The lines were ordered by increasing GC content, and the columns as in Figure 4E. (C) GC content biases in the targets of various AGO proteins and miRNAs. The AGO targets (in yellow) were defined using CLIP experiments and the miRNAs targets (in violet) using miRTarbase. The data are represented as in (A). (D) Heatmap representation of AGO and miRNAs depending on the behavior of their mRNA targets in the different datasets. The data were represented as in (B), using the same color code.
Figure 5. GC biases in the targets of various RNA decay factors, translation regulators and miRNAs. (A) GC content biases in the targets of various RBPs. The targets of the indicated factors were defined using CLIP experiments or motif analysis (see Materials and methods). The boxplots represent the distribution of the GC content of their gene. The distribution for all mRNAs is presented for comparison (in gray) and the red dashed line indicates its median value. (B) Heatmap representation of the different factors depending on the behavior of their mRNA targets in the different datasets. The lines were ordered by increasing GC content, and the columns as in Figure 4E. (C) GC content biases in the targets of various AGO proteins and miRNAs. The AGO targets (in yellow) were defined using CLIP experiments and the miRNAs targets (in violet) using miRTarbase. The data are represented as in (A). (D) Heatmap representation of AGO and miRNAs depending on the behavior of their mRNA targets in the different datasets. The data were represented as in (B), using the same color code.CC BY · PubMed Central
Figure 5—figure supplement 1. Targets of group I and II regulators (part I). (A) Number of transcripts in the mRNA subsets analyzed in (B–D) and in Figure 5A,B, Figure 5—figure supplement 2A–D and Figure 5—figure supplement 4. (B) Sensitivity of the different mRNA subsets to DDX6 silencing in HEK293 (upper panel) and K562 cells (lower panel). The behavior of all mRNAs is shown for comparison (in gray). Only the targets of group I factorsaccumulated following DDX6 depletion. While the presented analysis uses the YTHDF2 targets reported in Vindry et al. (2019), the same pattern was observed using the targets reported in Ozgur et al. (2010). (C) Sensitivity of the different mRNA subsets to XRN1 silencing in HeLa (upper panel) and HCT116 cells (lower panel). Only the targets of group I factors accumulated following XRN1 depletion in HCT116 cells, and only SMG6 targets did so in HeLa cells. (D) Sensitivity of the different mRNA subsets to PAT1B silencing in HEK293 cells. Only the targets of group II factors accumulated following PAT1B depletion, except ATXN2. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; ns, non-significant.
Figure 5—figure supplement 1. Targets of group I and II regulators (part I). (A) Number of transcripts in the mRNA subsets analyzed in (B–D) and in Figure 5A,B, Figure 5—figure supplement 2A–D and Figure 5—figure supplement 4. (B) Sensitivity of the different mRNA subsets to DDX6 silencing in HEK293 (upper panel) and K562 cells (lower panel). The behavior of all mRNAs is shown for comparison (in gray). Only the targets of group I factorsaccumulated following DDX6 depletion. While the presented analysis uses the YTHDF2 targets reported in Vindry et al. (2019), the same pattern was observed using the targets reported in Ozgur et al. (2010). (C) Sensitivity of the different mRNA subsets to XRN1 silencing in HeLa (upper panel) and HCT116 cells (lower panel). Only the targets of group I factors accumulated following XRN1 depletion in HCT116 cells, and only SMG6 targets did so in HeLa cells. (D) Sensitivity of the different mRNA subsets to PAT1B silencing in HEK293 cells. Only the targets of group II factors accumulated following PAT1B depletion, except ATXN2. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; ns, non-significant.CC BY · PubMed Central
Figure 5—figure supplement 2. Targets of group I and II regulators (part II). (A) PB localization of the different mRNA subsets. Only the targets of group II factors were enriched in PBs, except ATXN2. Note also that SMG6 targets and TOP mRNAs were particularly excluded from PBs. (B) Sensitivity of the different mRNA subsets to DDX6 silencing in HEK293 in terms of translation rate(polysomal/total mRNA ratio). Only the targets of group II factors were translationally derepressed, except ATXN2. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001. (C, D) The targets of group I regulators were stabilized after DDX6 (C) and XRN1 (D) silencing mostly like other PB-excluded mRNAs. The same analysis as in Figure 5—figure supplement 1B and C was conducted separately on PB-excluded (PB-out) and PB-enriched (PB-in) mRNAs. The distribution of all mRNAs (in gray), all PB-out and all PB-in mRNAs (dashed boxes) are shown for comparison. Two-tail Mann-Whitney test was performed with respect to cognate PB-out or PB-in mRNAs. **, p<0.0001; *, p=0.006 for (C) and 0.02 for (D); ns, non-significant. (E) DDX6-dependent mRNA decay does not correlate with polysome engagement. The fraction of mRNAs in polysomes (polysomal/total cpm) was calculated using the RNA-Seq data from control HEK293 cells. Transcripts were subdivided into 15 bins (1000 transcripts each) of increasing value. The boxplot represents their fold-change (FC) distribution in total RNA following DDX6 silencing.
Figure 5—figure supplement 2. Targets of group I and II regulators (part II). (A) PB localization of the different mRNA subsets. Only the targets of group II factors were enriched in PBs, except ATXN2. Note also that SMG6 targets and TOP mRNAs were particularly excluded from PBs. (B) Sensitivity of the different mRNA subsets to DDX6 silencing in HEK293 in terms of translation rate(polysomal/total mRNA ratio). Only the targets of group II factors were translationally derepressed, except ATXN2. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001. (C, D) The targets of group I regulators were stabilized after DDX6 (C) and XRN1 (D) silencing mostly like other PB-excluded mRNAs. The same analysis as in Figure 5—figure supplement 1B and C was conducted separately on PB-excluded (PB-out) and PB-enriched (PB-in) mRNAs. The distribution of all mRNAs (in gray), all PB-out and all PB-in mRNAs (dashed boxes) are shown for comparison. Two-tail Mann-Whitney test was performed with respect to cognate PB-out or PB-in mRNAs. **, p<0.0001; *, p=0.006 for (C) and 0.02 for (D); ns, non-significant. (E) DDX6-dependent mRNA decay does not correlate with polysome engagement. The fraction of mRNAs in polysomes (polysomal/total cpm) was calculated using the RNA-Seq data from control HEK293 cells. Transcripts were subdivided into 15 bins (1000 transcripts each) of increasing value. The boxplot represents their fold-change (FC) distribution in total RNA following DDX6 silencing.CC BY · PubMed Central
Figure 5—figure supplement 3. Targets of the miRNA pathway. (A,B) Number of AGO targets (A) and miRNA targets (B) analyzed in (C–H) and in Figure 5C,D. (C–E) AGO targets are enriched in PBs (C), translationally derepressed (D) and little or not stabilized (E) following DDX6 silencing. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; *, p=0.002, x, p=0.04; ns, non-significant. (F) AGO targets excluded from PBs accumulate following DDX6 silencing like other PB-excluded transcripts. The analysis was performed as in Figure 5—figure supplement 2C. Two-tail Mann-Whitney test was performed with respect to cognate PB-out or PB-in mRNAs. x, p=0.04; ns, non-significant. (G) AGO targets accumulate following PAT1B silencing. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; x, p=0.04. (H) The extent of PB enrichment of miRNA targets depends on the miRNA. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; ns, non-significant.
Figure 5—figure supplement 3. Targets of the miRNA pathway. (A,B) Number of AGO targets (A) and miRNA targets (B) analyzed in (C–H) and in Figure 5C,D. (C–E) AGO targets are enriched in PBs (C), translationally derepressed (D) and little or not stabilized (E) following DDX6 silencing. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; *, p=0.002, x, p=0.04; ns, non-significant. (F) AGO targets excluded from PBs accumulate following DDX6 silencing like other PB-excluded transcripts. The analysis was performed as in Figure 5—figure supplement 2C. Two-tail Mann-Whitney test was performed with respect to cognate PB-out or PB-in mRNAs. x, p=0.04; ns, non-significant. (G) AGO targets accumulate following PAT1B silencing. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; x, p=0.04. (H) The extent of PB enrichment of miRNA targets depends on the miRNA. Two-tail Mann-Whitney test was performed with respect to all mRNAs. **, p<0.0001; ns, non-significant.CC BY · PubMed Central
Figure 5—figure supplement 4. Targets of the group I and II regulators behave like mRNAs of similar GC content. (A–E) The whole transcriptome was binned depending on its GC content (bin size of 500 mRNAs). The median GC content of each bin was represented as a function of its median fold-change in the various datasets (in gray). The same analysis was performed for the targets of group I (in green) and II (in orange) regulators, and superimposed to the graphs. The dashed line indicates the median GC content for all mRNAs. Overall, in all datasets, the targets of the different regulators behave similarly to mRNAs of same GC content. For figure clarity, only the few regulators leading to distinct behaviors were indicated.
Figure 5—figure supplement 4. Targets of the group I and II regulators behave like mRNAs of similar GC content. (A–E) The whole transcriptome was binned depending on its GC content (bin size of 500 mRNAs). The median GC content of each bin was represented as a function of its median fold-change in the various datasets (in gray). The same analysis was performed for the targets of group I (in green) and II (in orange) regulators, and superimposed to the graphs. The dashed line indicates the median GC content for all mRNAs. Overall, in all datasets, the targets of the different regulators behave similarly to mRNAs of same GC content. For figure clarity, only the few regulators leading to distinct behaviors were indicated.CC BY · PubMed Central
Figure 5—figure supplement 5. miRNA targets behave like mRNAs of similar GC content. (A–E) The analysis was conducted as in Figure 5—figure supplement 4. Overall, the targets of the different miRNAs behave similarly to mRNAs of same GC content. (F) The GC content of miRNA targets correlates with the GC content of the miRNA itself. rp, Pearson correlation coefficient.
Figure 5—figure supplement 5. miRNA targets behave like mRNAs of similar GC content. (A–E) The analysis was conducted as in Figure 5—figure supplement 4. Overall, the targets of the different miRNAs behave similarly to mRNAs of same GC content. (F) The GC content of miRNA targets correlates with the GC content of the miRNA itself. rp, Pearson correlation coefficient.CC BY · PubMed Central
Figure 6. The GC content of the CDS and the 3’UTR both contribute to PB localization. (A) General importance of the CDS. Transcripts were subdivided into six classes depending on the GC content of their 3’UTR (from <40 to >55%). The boxplots represent the distribution of their CDS GC content at position 3 (GC3) in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. (B) Importance of the 3’UTR for GC-rich CDSs. Transcripts were subdivided into eight classes depending on their GC3 (from <40 to >70%). The boxplots represent the distribution of their 3’UTR GC content in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. (C) The transcripts of haplo-insufficiency genes are enriched in PBs. The haplo-insufficiency score is the probability that a gene is haplo-insufficient, as taken from the Huang et al. (2010) study. The analysis was performed for PB-enriched (PB-in, n = 4646, median score 0.26) and PB-excluded (PB-out, n = 4205, median score 0.17) mRNAs. The difference of distribution of haplo-insufficiency scores was statistically significant using a two tail Mann-Whitney test: p<0.0001. The results were similar using Steinberg et al. (2015) scores. (D) Protein yield is higher from a GC-rich than an AU-rich CDS. HEK293 cells were transfected with Rluc reporters differing by the GC content of their CDS, along with a control Fluc plasmid. After 24 hr, mRNA levels were measured by qPCR and protein levels by luciferase activity. The Rluc to Fluc ratio for the GC-rich reporter was set to 100 (n = 3). Error bars, SD. (E, F) Preferential localization of AU-rich transcripts in PBs. HEK293 cells expressing GFP-LSM14A were transfected with the AU-rich and GC-rich Rluc reporters and the localization of the Rluc transcripts (in red) was analyzed by smiFISH. Representative cells are shown in (E). Bar, 5 µm. Arrows indicate the PBs enlarged above. The experiment was performed in duplicate (exp. 1 and 2) and repeated in HEK293 cells where PBs were immunostained using DDX6 antibodies (exp. 3 and 4). The percentage of PBs containing clusters of Rluc transcripts in the four experiments is represented in (F). Exp.1: 56/75 PBs from 21/27 cells; exp.2: 87/75 PBs from 38/35 cells; exp.3: 31/32 PBs from 15/19 cells; exp.4: 72/83 PBs from 34/41 cells (G) Assembly of PB-like granules in cell-free extracts from HEK293 cells expressing GFP-LSM14A. The scheme recapitulates the main steps of the assay. Fluorescence microscopy images show that PBs in cells, PBs after cell lysis, and reconstituted PB-like granules have similar size. Bar, 10 µm. (H) AU-rich RNA favors the formation of PB-like granules. PB-like granules were assembled in cell-free extracts in the presence of AU-rich or GC-rich RNA, and counted by flow cytometry. Their number in the absence of added RNA was set to 100 (n = 3 experiments in duplicate, using two independent cell-free extracts and RNA preparations). Error bars, SD.
Figure 6. The GC content of the CDS and the 3’UTR both contribute to PB localization. (A) General importance of the CDS. Transcripts were subdivided into six classes depending on the GC content of their 3’UTR (from <40 to >55%). The boxplots represent the distribution of their CDS GC content at position 3 (GC3) in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. (B) Importance of the 3’UTR for GC-rich CDSs. Transcripts were subdivided into eight classes depending on their GC3 (from <40 to >70%). The boxplots represent the distribution of their 3’UTR GC content in PB-enriched (PB-in) and PB-excluded (PB-out) mRNAs. (C) The transcripts of haplo-insufficiency genes are enriched in PBs. The haplo-insufficiency score is the probability that a gene is haplo-insufficient, as taken from the Huang et al. (2010) study. The analysis was performed for PB-enriched (PB-in, n = 4646, median score 0.26) and PB-excluded (PB-out, n = 4205, median score 0.17) mRNAs. The difference of distribution of haplo-insufficiency scores was statistically significant using a two tail Mann-Whitney test: p<0.0001. The results were similar using Steinberg et al. (2015) scores. (D) Protein yield is higher from a GC-rich than an AU-rich CDS. HEK293 cells were transfected with Rluc reporters differing by the GC content of their CDS, along with a control Fluc plasmid. After 24 hr, mRNA levels were measured by qPCR and protein levels by luciferase activity. The Rluc to Fluc ratio for the GC-rich reporter was set to 100 (n = 3). Error bars, SD. (E, F) Preferential localization of AU-rich transcripts in PBs. HEK293 cells expressing GFP-LSM14A were transfected with the AU-rich and GC-rich Rluc reporters and the localization of the Rluc transcripts (in red) was analyzed by smiFISH. Representative cells are shown in (E). Bar, 5 µm. Arrows indicate the PBs enlarged above. The experiment was performed in duplicate (exp. 1 and 2) and repeated in HEK293 cells where PBs were immunostained using DDX6 antibodies (exp. 3 and 4). The percentage of PBs containing clusters of Rluc transcripts in the four experiments is represented in (F). Exp.1: 56/75 PBs from 21/27 cells; exp.2: 87/75 PBs from 38/35 cells; exp.3: 31/32 PBs from 15/19 cells; exp.4: 72/83 PBs from 34/41 cells (G) Assembly of PB-like granules in cell-free extracts from HEK293 cells expressing GFP-LSM14A. The scheme recapitulates the main steps of the assay. Fluorescence microscopy images show that PBs in cells, PBs after cell lysis, and reconstituted PB-like granules have similar size. Bar, 10 µm. (H) AU-rich RNA favors the formation of PB-like granules. PB-like granules were assembled in cell-free extracts in the presence of AU-rich or GC-rich RNA, and counted by flow cytometry. Their number in the absence of added RNA was set to 100 (n = 3 experiments in duplicate, using two independent cell-free extracts and RNA preparations). Error bars, SD.CC BY · PubMed Central
Figure 6—figure supplement 1. Role of PB localization in XRN1 and DDX6 sensitivity and importance of the coding property for PB localization. (A) AU-rich mRNAs are better XRN1 targets when excluded from PBs than when enriched in PBs. The analysis was conducted as in Figure 5—figure supplement 4B for all mRNAs (in gray) and for PB-excluded mRNAs (PB-out, in dark blue), except that mRNAs with GC values under 50% were binned by 200 to account for their low number. (B) AU-rich mRNAs are not better DDX6 decay targets when they are excluded from PBs. The dataset after DDX6 silencing in HEK293 cells was analyzed as in (A). (C) AU-rich mRNAs excluded from PBs preferentially encode proteins associated with membranous organelles or secreted. GO analysis was performed using the SLIM categories of cellular components to compare PB-excluded AU-rich mRNAs (GC <49.1%, 556 annotated) to all AU-rich mRNAs (GC <49.1%, 3896 annotated). The Volcano plot (left panel) represents the most enriched categories (FDR < 0.01). Both size and color of the dots are proportional to the size of the category within the query database. Transcripts corresponding to the mitochondria-related categories (green), other membranous organelles (red), secreted proteins (blue) and ribosomal proteins (yellow) were further subjected to a Venn diagram analysis (right panel). (D) LncRNAs poorly accumulate in PBs even when they are AU-rich. LncRNAs (n = 2589) were subdivided into 6 classes of increasing GC content and their accumulation in PBs was represented as for mRNAs in Figure 1B. rs, Spearman correlation coefficient. (E) Number of transcripts in the different classes analyzed in (D) and in Figure 6A,B.
Figure 6—figure supplement 1. Role of PB localization in XRN1 and DDX6 sensitivity and importance of the coding property for PB localization. (A) AU-rich mRNAs are better XRN1 targets when excluded from PBs than when enriched in PBs. The analysis was conducted as in Figure 5—figure supplement 4B for all mRNAs (in gray) and for PB-excluded mRNAs (PB-out, in dark blue), except that mRNAs with GC values under 50% were binned by 200 to account for their low number. (B) AU-rich mRNAs are not better DDX6 decay targets when they are excluded from PBs. The dataset after DDX6 silencing in HEK293 cells was analyzed as in (A). (C) AU-rich mRNAs excluded from PBs preferentially encode proteins associated with membranous organelles or secreted. GO analysis was performed using the SLIM categories of cellular components to compare PB-excluded AU-rich mRNAs (GC <49.1%, 556 annotated) to all AU-rich mRNAs (GC <49.1%, 3896 annotated). The Volcano plot (left panel) represents the most enriched categories (FDR < 0.01). Both size and color of the dots are proportional to the size of the category within the query database. Transcripts corresponding to the mitochondria-related categories (green), other membranous organelles (red), secreted proteins (blue) and ribosomal proteins (yellow) were further subjected to a Venn diagram analysis (right panel). (D) LncRNAs poorly accumulate in PBs even when they are AU-rich. LncRNAs (n = 2589) were subdivided into 6 classes of increasing GC content and their accumulation in PBs was represented as for mRNAs in Figure 1B. rs, Spearman correlation coefficient. (E) Number of transcripts in the different classes analyzed in (D) and in Figure 6A,B.CC BY · PubMed Central
Figure 6—figure supplement 2. The GC content of reporters RNAs is key for PB localization. (A) Specificity of the GC-rich and AU-rich Rluc probes. HEK293 cells were transfected or not with the GC-rich and AU-rich Rluc reporter constructs. After 24 hr, cells were analyzed by smiFISH using the GC-rich or AU-rich probes set, as indicated. Bar, 5 µm. (B) AU-rich transcripts localize preferentially in PBs. HEK293 cells transfected with the AU-rich or GC-rich Rluc reporters were subjected to smiFISH analysis, using specific AU-rich or GC-rich probes set to localize the Rluc transcripts (in red), and anti-DDX6 antibody to localize PBs (in blue). Representative cells are shown. Arrows indicate PBs enlarged above. Bar, 5 µm. (C) Recombinant DDX6 protein used in the cell-free assay. The double purified CBP-DDX6-HIS protein was migrated on a polyacrylamide gel along with a molecular weight marker (MW) before Coomassie staining. (D) PB-like granules formed in a cell-free extract can be counted by flow cytometry. The cell-free extract was prepared from HEK293 cells expressing GFP-LSM14A. Without introducing a trigger threshold, the detector was saturated with electronic noise, as illustrated using buffer only (upper panel), and the machine could not be used in the counting mode. After gating as indicated by the red lines (lower panel), PBs were readily detected in the resuspended lysate pellet (pelleted PBs), while almost absent in the supernatant corresponding to the cell free extract. Upon DDX6 addition to the cell-free extract, PB-like granules were detected at a size similar to PBs. Their fluorescence intensity was slightly lower, which likely reflected the limited amount of soluble GFP-LSM14A in the cell-free extract. (E) PB-like granules form in a DDX6-dependent manner. Increasing amount of recombinant DDX6 protein was added to the cell-free extract, and granules were counted by flow cytometry, as described in (D). (F) GC-rich and AU-rich mRNAs used in the cell-free extract. After in vitro transcription from GC-rich and AU-rich plasmids, the RNAs were migrated in an agarose gel to compare their size and concentration. Note that the conditions do not allow for an absolute size estimation, as the gel is non-denaturing and the molecular weight marker (MW) is DNA.
Figure 6—figure supplement 2. The GC content of reporters RNAs is key for PB localization. (A) Specificity of the GC-rich and AU-rich Rluc probes. HEK293 cells were transfected or not with the GC-rich and AU-rich Rluc reporter constructs. After 24 hr, cells were analyzed by smiFISH using the GC-rich or AU-rich probes set, as indicated. Bar, 5 µm. (B) AU-rich transcripts localize preferentially in PBs. HEK293 cells transfected with the AU-rich or GC-rich Rluc reporters were subjected to smiFISH analysis, using specific AU-rich or GC-rich probes set to localize the Rluc transcripts (in red), and anti-DDX6 antibody to localize PBs (in blue). Representative cells are shown. Arrows indicate PBs enlarged above. Bar, 5 µm. (C) Recombinant DDX6 protein used in the cell-free assay. The double purified CBP-DDX6-HIS protein was migrated on a polyacrylamide gel along with a molecular weight marker (MW) before Coomassie staining. (D) PB-like granules formed in a cell-free extract can be counted by flow cytometry. The cell-free extract was prepared from HEK293 cells expressing GFP-LSM14A. Without introducing a trigger threshold, the detector was saturated with electronic noise, as illustrated using buffer only (upper panel), and the machine could not be used in the counting mode. After gating as indicated by the red lines (lower panel), PBs were readily detected in the resuspended lysate pellet (pelleted PBs), while almost absent in the supernatant corresponding to the cell free extract. Upon DDX6 addition to the cell-free extract, PB-like granules were detected at a size similar to PBs. Their fluorescence intensity was slightly lower, which likely reflected the limited amount of soluble GFP-LSM14A in the cell-free extract. (E) PB-like granules form in a DDX6-dependent manner. Increasing amount of recombinant DDX6 protein was added to the cell-free extract, and granules were counted by flow cytometry, as described in (D). (F) GC-rich and AU-rich mRNAs used in the cell-free extract. After in vitro transcription from GC-rich and AU-rich plasmids, the RNAs were migrated in an agarose gel to compare their size and concentration. Note that the conditions do not allow for an absolute size estimation, as the gel is non-denaturing and the molecular weight marker (MW) is DNA.CC BY · PubMed Central
Figure 7. Schematic representation recapitulating the features of mRNA post-transcriptional regulation depending on their GC content.
Figure 7. Schematic representation recapitulating the features of mRNA post-transcriptional regulation depending on their GC content.CC BY · PubMed Central
Figure 7—figure supplement 1. Distribution of the gene GC content in various eukaryotic genomes. The graphs were constructed as in Figure 1—figure supplement 1F. The median value for each genome is indicated below.
Figure 7—figure supplement 1. Distribution of the gene GC content in various eukaryotic genomes. The graphs were constructed as in Figure 1—figure supplement 1F. The median value for each genome is indicated below.CC BY · PubMed Central

Verification

Claims the paper makes

The claims and statistics of this paper haven’t been checked yet.

Run the check — free, no account needed. It reads the open-access text and checks the sources it cites.

A few free checks a day · papers already verified open instantly · results are public

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

The self-check needs the paper’s full text.

Citations

Cites

Cited by

Sorted by citations

See all 256

Related

See all

Metadata

Authors
Courel, M., Clément, Y., Bossevain, C., Foretek, D., Vidal Cruchez, O., Yi, Z., Bénard, M., Benassy, M.-N., Kress, M., Vindry, C., Ernoult-Lange, M., Antoniewski, C., Morillon, A., Brest, P., Hubstenberger, A., Roest Crollius, H., Standart, N., & Weil, D.
Affiliations
  • Centre National de la Recherche Scientifique (FR)
  • École Normale Supérieure - PSL (FR)
  • Université Côte d'Azur, CNRS, INSERM, IRCAN, FHU-OncoAge, Nice, France
  • Université Côte d'Azur, CNRS, INSERM, iBV, Nice, France
  • University of Cambridge (GB)
  • Sorbonne Université, CNRS, Institut de Biologie Paris Seine (IBPS), Laboratoire de Biologie du Développement, Paris, France
  • Ecole Normale Supérieure, Institut de Biologie de l'ENS, IBENS, Paris, France
ORCID
Venue
eLife
Publisher
eLife Sciences Publications, Ltd
Published
2019-12-19
Volume
8
Pages
e49708
ISSN
2050-084X (ISSN-L)
DOI
10.7554/elife.49708
PMC
PMC6944446
OpenAlex
W2994774550
Cited by
250
References
83
Citations by year
2026: 38 · 2025: 40 · 2024: 44 · 2023: 35 · 2022: 25 · 2021: 44 · 2020: 17 · 2019: 7
FWCI
7.98
Type
journal article
Language
en
Open access
Licence
Locations
  • eLife (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • Europe PMC (PubMed Central) (publishedVersion, cc-by)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
  • HAL (submittedVersion, other-oa)
  • HAL (submittedVersion)
  • HAL (submittedVersion, cc-by)
  • DOAJ (DOAJ: Directory of Open Access Journals) (submittedVersion, cc-by-sa)
  • Apollo (University of Cambridge) (submittedVersion, cc-by)
  • Apollo (cc-by)
Funders
  • Association pour la Recherche sur le Cancer (Subvention Fixe)
  • Agence Nationale de la Recherche (ANR-14-CE09-0013-01, ANR-11-LABX-0028-01)
  • European Research Council (DARK consolidator grant)
  • Canceropôle PACA
  • Biotechnology and Biological Sciences Research Council
  • Isaac Newton Trust
  • Fondation Philippe Wiener - Maurice Anspach
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA Research and Splicing
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA and protein synthesis mechanisms
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › RNA regulation and disease
Keywords
GC content, mRNA stability, mRNA degradation, P-bodies, codon usage, translational regulation, post-transcriptional regulation, RNA-binding proteins, microRNA, human cells
MeSH terms
Base Composition, Gene Expression Regulation, Humans, RNA, Messenger, Protein Biosynthesis, RNA Stability, RNA, Messenger, Stored, MicroRNAs, Transcriptome

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review