Association of expression kinetics with mRNA stability states of TNF-inducible genes in RA FLS.

<p>For (a-b), two biological replicates of RA FLS (derived from two different RA patients) were exposed to a single dose of TNF (10 ng/ml) for 1-72h. Subsequently, Act D (10 μg/ml) was added for 3h and gene expression was measured by RNA sequencing. 386 genes were identified as highly induced (≥5-fold) by TNF at any time point and were clustered into 6 clusters with distinct kinetics of peak expression. (a), Heatmap illustrating the expression kinetics of the 6 clusters (red represents the maximum and blue the minimum expression level across the lane). (b), Stacked bar graphs illustrating the stability states of genes for Cluster 1, Clusters 2 &3, Cluster 4, and Clusters 5 & 6. For (c-f), RA FLS were exposed to a single dose of TNF (10 ng/ml) for 1–72 hours. Primers specific for the eighth intronic region of <i>MMP3</i> and for the first intronic region of <i>CCL5</i> were designed to capture primary transcripts (PT) of <i>MMP3</i> and <i>CCL5</i>. qPCR was used to measure the levels of PT and total mRNA of MMP3 (c-d) and CCL5 (e-f). Cumulative results from six independent experiments are shown. Values were normalized relative to mRNA for GAPDH and are presented as mean ±SEM.</p

TNF induces a temporal switch from an early program dominated by unstable transcripts to a late program with expansion of the pool of stable transcripts.

<p>(a-b), Density-plots illustrating the distribution of TNF-inducible genes (≥2-fold by TNF at 1 (a) and 72h (b)) based on their mRNA stability (shades of red represent areas with the highest numbers of genes (highest density)). The mRNA stability was calculated genome-wide in 2 biologic replicates as the ratio of RPKM levels at the TNF+Act D condition divided to the RPKM levels at the TNF condition. This ratio ranges from 0 (very unstable transcripts) to 1 (very stable transcripts). (c), RA FLS were stimulated with a single dose of TNF (10 ng/ml) for 72 hours and then Act D (10 μg/ml) was added for 1 hour. Subsequently cells were washed and fresh serum-free medium + Act D was replenished. Supernatants were collected 6 hours later and the protein levels of IL-6, CXCL1, CCL2, IL-8, RANTES and IP-10 were measured by magnetic bead-based multiplex immunoassay. Values are the mean ±SEM of three independent experiments. <i>P</i> values were calculated by one-way ANOVA and Tukey post-test analysis (* = <i>p</i><0.05, ** = <i>p</i><0.01, and *** = <i>p</i><0.001).</p

Genome-wide identification of transcripts stabilized by TNF in RA FLS.

<p>Two biologic replicates of RA FLS (derived from two different RA patients) were exposed to a single dose of TNF (10 ng/ml) for 1 or 72h. Subsequently, Act D was added for 3h and gene expression was measured by RNA sequencing. The degree of TNF-induced mRNA stabilization was calculated as the log₂ difference of TNF+Act D/TNF ratio between 1 and 72h of TNF stimulation and the adjusted p values of TNF-induced stabilization were calculated by RiboDiff. (a), Scatter-plot of the genes displaying TNF-induced mRNA stabilization comparing the degree of mRNA stabilization (y axis) to the adjusted p values of the stabilizing effect of TNF (x-axis). (b), The top 40 genes displaying the highest TNF-induced mRNA stabilization ranked by the degree of stabilization. (c), Enriched biological processes identified by GSEA/MSigDB pathway analysis of the top 10% of the genes (n = 593) displaying the highest degree of TNF-induced mRNA stabilization.</p

TNF induces late stabilization of IL-6 mRNA in RA FLS.

<p>RA FLS were exposed to a single dose of TNF (10 ng/ml) for 1-72h. Subsequently, actinomycin D (Act D; 10 μg/ml) or triptolide (1 μM), or flavopiridol (Flav; 0.5 μM) was added for 1 or 3h and qPCR was used to measure the mRNA levels of IL-6 (a-c), the primary transcripts (PT) of <i>IL6</i> (d-e), and the mRNA levels of CCL5 (f). For (b-c), the remaining expression of IL-6 after exposure to inhibitors of transcription (Act D, triptolide, and flavopiridol) was calculated as % of the IL-6 mRNA expressed in the absence of inhibitor at the corresponding TNF-stimulated condition. For (a-b) and (f), cumulative results from 7 independent experiments are shown. For (d-e), FLS were exposed to a single dose of TNF (10 ng/ml) for 72 hours and then inhibitors (Act D or flavopiridol) were added for the indicated time points to block active transcription. Primers specific for the fourth intronic region of <i>IL6</i> were designed to capture primary transcripts of <i>IL6</i>.Values were normalized relative to GAPDH mRNA and presented as mean ±SEM. GAPDH was considered an appropriate internal control for normalizing qPCR results since TNF stimulation had no impact on expression levels and stability status of GAPDH mRNA (<a href="http://www.plosone.org/article/info:doi/10.1371/journal.pone.0179762#pone.0179762.s001" target="_blank">S1 Table</a>). <i>P</i> values were calculated by one-way ANOVA and Tukey post-test analysis (* = <i>p</i><0.05, ** = <i>p</i><0.01, *** = <i>p</i><0.001, ns = not significant, and ND = not detected).</p

TNF induces expression of mRNA-stabilizing pathways and mRNA stabilization is MAPK-dependent.

<p>(a), RNA sequencing was performed in 2 biological replicates (derived from two different RA patients) of TNF-stimulated RA FLS and Panther-Gene Ontology was used to evaluate their enrichment for the biological process “Regulation of RNA stability” (GO:0043487 or GO:0043488). F.E = fold enrichment and ns = not significant. (b-h), RA FLS were exposed to a single dose of TNF (10 ng/ml) for 72h and then Act D (10 μg/ml) was added for 20 mins to block active transcription. Subsequently, the cells were treated for 4h with SB202190 (p38 inhibitor) alone or in various combinations with U0126 (MEK inhibitor) and SP600125 (JNK inhibitor). qPCR was used to measure the mRNA levels of CCL5 (b), IL-6 (c), IL-8 (d), CXCL3 (e), CCL2 (f), PTGS2 (g), and CXCL1 (h). Cumulative results from 4 independent experiments are shown. Values were normalized relative to GAPDH mRNA and presented as mean ±SEM. The mRNA expression at the TNF+Act D condition was set to 100 and the mRNA expression at all the other conditions was calculated as % of the TNF+Act D condition. <i>P</i> values were calculated by one-way ANOVA and Tukey post-test analysis (* = <i>p</i><0.05, ** = <i>p</i><0.01, *** = <i>p</i><0.001, and ns = not significant).</p

Genome-wide evaluation of mRNA stability states of expressed genes in RA FLS.

<p>(a-c), Gene tracks showing sequencing reads from RNA sequencing mapped to <i>CCL20</i> (a), <i>JUN</i> (b) and <i>IRF1</i> (c) genes. The sequencing reads after TNF stimulation for 1 hour without (blue) or with Act D (orange) are shown. (d), Stacked bar graphs illustrating the mRNA stability states of genes expressed in unstimulated (Control) and TNF-stimulated FLS (1, 3, 24 and 72 hours of TNF stimulation). The mRNA stability status was calculated as the ratio of expression levels at the TNF+Act D condition divided to the expression levels at the TNF condition. This ratio ranges from 0 to 1 and classifies genes to a spectrum from very unstable to very stable transcripts. The expressed genes were classified into five groups with distinct stability states and the size of each group is represented as % of total number of expressed genes for each condition.</p

Genome-wide association of the gene expression levels with the mRNA stability states in rheumatoid arthritis (RA) fibroblast-like synoviocytes (FLS).

<p>(a), Table of the top 10% of highly expressed genes at 1, 3, 24 and 72 hours of TNF stimulation grouped by their mRNA stability status. The number of genes included in each stability group is presented. (b), Table of the highly expressed genes with very unstable mRNAs. Their ranking by level of expression is presented in parentheses. (c), Kinetics of expression and TNF-induced dynamics of mRNA stability for IL-6, IL-8 and CCL2. Shades of blue represent unstable status and shades of red represent stable status.</p

Scatterplots comparing the expression levels to the mRNA stability states of the expressed genes in RA FLS.

<p>Two biological replicates of RA FLS (derived from two different RA patients) were exposed to a single dose of TNF (10 ng/ml) for 1, 3, 24, or 72 hours. Subsequently, actinomycin D (Act D, 10μg/ml) was added for 3 hours to block active transcription and gene expression was measured by RNA sequencing. RPKM values were generated using CuffDiff2. The mRNA stability status was calculated genome-wide as the ratio of RPKM levels at the TNF+Act D condition divided to the RPKM levels at the TNF condition. This ratio ranges from 0 to 1 and classifies genes to a spectrum from very unstable to very stable transcripts. The genes expressed at 1 (a), 3 (b), 24 (c), and 72 (d) hours of TNF stimulation were plotted based on their expression levels and the mRNA stability states. Shades of blue represent the region of unstable genes, and shades of red represent the zone of stable genes.</p