Single-cell RNA-Seq analysis identifies a noncoding interleukin 4 (IL-4) RNA that post-transcriptionally up-regulates IL-4 production in T helper cells

High-throughput sequencing has revealed a tremendous complexity of cellular transcriptomes, which is partly due to the generation of multiple alternative transcripts from a single gene locus. Because alternative transcripts often have low abundance in bulk cells, the functions of most of these transcripts and their relationship with their canonical counterparts remain unclear. Here we applied single-cell RNA-Seq to analyze the transcriptome complexity of in vitro–differentiated, murine type 2 T helper (Th2) cells. We found that cytokine gene transcripts contribute most of the intercellular heterogeneity, with a group of universal cytokines, including interleukins 1a, 2, 3, and 16, being bimodally expressed. At the single-cell level, use of alternative promoters prevalently generated alternative transcripts. For instance, although undetectable in bulk cells, a noncoding RNA isoform of IL-4 (IL4nc), which was driven by an intronic promoter in the IL-4 locus, was predominantly expressed in a subset of Th2 cells. IL4nc displayed distinct temporal expression patterns compared with the canonical IL-4 mRNA and post-transcriptionally promoted the production of IL-4 protein in Th2 cells. In conclusion, our findings reveal a mechanism whereby minor noncanonical transcripts post-transcriptionally regulate expression of their cognate canonical genes.

Upon antigen stimulation, naïve helper T cells (Th) 2 differentiate into various effector subsets (i.e. Th1, Th2, Th17, and regulatory T cells (Treg)) to orchestrate a proper immune response against foreign pathogens (1). Each subset of T helper cells expresses their characteristic cytokines (i.e. Th1 cells express IFN-␥ but not IL-4, whereas Th2 cells express IL-4 but not IFN-␥) and master transcription factors (i.e. Tbet for Th1 cells and GATA-3 for Th2 cells) (2)(3)(4). In addition to these lineage-specific factors, effector Th cells also express common cytokines shared among virtually all effector T cells (i.e. IL-2 and tumor necrosis factor ␣) to carry out their physiological functions. The variable combinations of the unique and common factors further increase the diversity of effector Th cells. Although the lineage decision and functions of Th subsets have been intensively studied (2,5,6), the intra-subset diversity and its origin remain to be fully elucidated.
As a prototype of the Th subset, Th2 cells are induced by extracellular pathogens to initiate the humoral immune response. Th2 cells can be in vitro differentiated from naïve CD4ϩ T cells by antigen and IL-4 stimulation (7)(8)(9)(10). Effector activities of Th2 cells depend on an array of secreted proteins, including IL-4, IL-5, IL-9, and IL-13 (collectively called type 2 cytokines) (11,12) that have shared and unique functions. IL-4 plays a key role in both humoral and cellular immunity, including promoting proliferation and survival of activated B cells, and inducing B cell class switching to IgG1 and IgE (13)(14)(15). IL-13, on the other hand, promotes mucus production by goblet cells and airway hyperresponsiveness, a hallmark of allergic asthma (2,(11)(12)(13). In addition to type 2 cytokines, Th2 cells also secrete universal cytokines shared with other Th subsets, such as IL-2, IL-1␣, and IL-3, although their regulation and function in Th2 cells are largely undetermined.
Recent progress in high-throughput sequencing has revealed the great complexity of cellular transcriptomes, partly contributed by previously underestimated noncoding transcripts. For instance, long noncoding RNAs can be generated from either noncoding DNA or from coding genes through alternative splicing or alternative promoter usage (16,17). Furthermore, many enhancers are actively transcribed into enhancer RNAs, which are often unspliced, nonpolyadenylated, and associated with active gene expression (18). However, the functions of most noncoding transcripts and their relation to the coding counterparts remain poorly understood, partially because of their low abundance at the bulk cell level and the difficulty to separate them from transcriptional noise.
In this study, we analyzed the transcriptome complexity of in vitro differentiated Th2 cells using single-cell RNA-Seq.

cro ARTICLE
Our analysis revealed that cytokines were the most variably expressed genes among Th2 cells, with universal cytokines, rather than Th2 specific cytokines, bimodally expressed among single cells. Furthermore, we identified a noncoding RNA abundantly transcribed from the Il4 locus in a subset of Th2 cells, which is driven by an enhancer within the second intron of the mouse IL4 gene. We demonstrated that induction of this noncoding isoform preceded the canonical IL-4 mRNA and promoted IL-4 translation following T cell antigen receptor (TCR) restimulation. Thus, our study charted the intra-subset diversity of Th2 cells and revealed a novel mechanism whereby translation of an immune effector cytokine is regulated by an alternative noncoding transcript from the same locus.

Single-cell RNA-Seq reveals transcriptome heterogeneity of Th2 cells
To explore transcriptome diversification during T cell differentiation, we analyzed in vitro differentiated Th2 cells using single-cell RNA-Seq (scRNA-Seq). Purified naïve CD4 T cells (CD4 ϩ CD62L hi CD44 lo CD25 Ϫ ) were differentiated into the Th2 lineage in vitro by anti-CD3, anti-CD28, IL-4, and anti-IFN-␥ stimulation (7,8). Five days after differentiation, flow cytometric analysis identified that ϳ25% cells expressed IL-4 protein following TCR restimulation (Fig. S1A), confirming successful differentiation of Th2 cells (7,8). 56 single cells were then captured from the restimulated Th2 cells for scRNA-Seq analysis using a Fluidigm C1 platform (19,20). As a control, pooled 10,000 Th2 cells were subjected to a standard RNA-Seq library preparation protocol. Each library was subsequently sequenced to an average depth of 24 million reads with a singleend 100-bp read length (Fig. S1B), which is well above the sequence depth required to cover the transcriptomes of single Th2 cells (Fig. S2). To validate that our scRNA-Seq assay truly reflected cellular heterogeneity, we compared global gene expression between single cells and pooled cells. The average global gene expression of single cells was very similar to the pooled sample (Pearson correlation, r ϭ 0.85) (Fig. 1A, top panel). The transcriptomes of any two single cells, however, differed greatly (average Pearson correlation, r Ͻ 0.6), suggesting a large degree of cellular heterogeneity of in vitro differentiated Th2 cells (Fig. 1, A center panel, and B). The pooled sample also correlated poorly with almost all single cells, indicating that traditional RNA-Seq analysis was not able to truly portray transcriptomes of single cells (Fig. 1A, bottom panel, and B). Further hierarchical cluster analysis revealed two distinct populations within the Th2 cells and one outlier, which was excluded from the subsequent analysis ( Fig. 1B and Fig. S1C). The scRNA-Seq analysis revealed cellular heterogeneity of in vitro differentiated Th2 cells, which was masked in bulk cell RNA-Seq analysis.
To quantitate variation of gene expression among single cells, we calculated coefficient variation (CV) of 10,565 genes detected in Th2 cells (TPM Ն 10 in at least one single cell or pooled sample). Consistent with a previous scRNA-Seq analysis in immune cells (21), many housekeeping genes were expressed abundantly (average TPM () Ͼ 8) and invariably among single Th2 cells (CV Ͻ 0.2). In contrast, although a number of cytokines were also abundantly expressed ( Ͼ 8), their expression varied greatly among individual cells (CV Ͼ 0.28) (Figs. 1C and 2A). Indeed, gene ontology analysis revealed that cytokine genes were enriched in a class of the most variably expressed genes ( Ͼ 8, CV Ͼ 0.28) in Th2 cells (Fig. 1D). As expected, cellular housekeeping functions, including RNA binding and ribosome components, were enriched in the most consistently expressed genes. Together, these data indicate that cytokine genes are major contributors of cellular heterogeneity of differentiated Th2 cells.

Universal cytokines are bimodally expressed in Th2 cells
Given the highly variable expression and the importance of cytokines in physiological functions of Th2 cells, we further analyzed the expression and distribution of each cytokine gene among single cells. To our surprise, some cytokines exhibited a bimodal expression pattern in which the gene was highly expressed in a subpopulation of cells but minimally expressed
To quantitate the bimodal distribution of gene expression among single cells, we calculated a bimodal index (BI) for every cytokine gene (23), which uses a two-component mixture model to measure the separation of two normal distributions. Setting BI Ͼ 1.8 as a threshold, over 50% of universal cytokines were bimodally expressed, in contrast to 25% of the Th2-specific cytokines (Fig. 2B). The average BI of universal cytokines (BI ϭ 1.46) was also significantly higher than that of the Th2-specific cytokines (BI ϭ 0.92) (Fig. 2C). Together, these data demonstrate that the universal cytokines more likely exhibit a unique bimodal expression pattern among single cells and contribute most to the cell-tocell variation in Th2 cells.

Alternative promoter use contributes to cellular transcriptome heterogeneity
Although our scRNA-Seq analysis revealed that IL-4 RNA was uniformly expressed among Th2 cells ( Fig. 2A), IL-4 protein was only detectable in ϳ25% cells, as determined by flow cytometric analysis (Fig. S1A). To reconcile the discrepancy between IL-4 RNA and protein expression among single cells, we examined transcripts that originated from the Il4 locus in depth. Although over 80% of cells (45 of 55) had robust transcription from the Il4 locus (TPM Ͼ 100), only ϳ60% (34 of 55) abundantly expressed the full-length IL-4 mRNA (IL4 fl ) (Fig.  3A). The remaining cells predominantly expressed an IL-4 noncoding isoform (IL4 nc ) that utilized an alternative transcription start site between canonical exon 2 and exon 3 in the IL4 gene (Fig. 3A). The anti-IL-4 antibody was only able to detect IL4 flexpressing K562 cells but not IL4 nc -expressing cells, confirming that IL4 nc RNA is noncoding and undetectable by the anti-IL-4 antibody (Fig. S3, A and B). Sequence analysis further confirmed that, except for the first exon, the remaining sequence of IL4 nc was identical to IL4 fl . Indeed, six cells only expressed IL4 nc but not IL4 fl (IL4 nc /IL4 fl Ͼ ϭ 7). Thus, transcription of a noncoding RNA from the Il4 locus at least partially contributes to the fewer IL-4 protein-expressing cells compared with the RNA-expressing cells during in vitro Th2 cell differentiation.
Similar to IL4, the CAPG gene, which is specifically up-regulated during Th2 response to parasite infection (24), also utilized an alternative transcription start site in most (90%) cells (Fig. 3B). Analysis of alternative first exon (AFE) usage globally identified 1620 AFE events (percent spliced in, ⌿ Ͼ 0.1 in at least 20% of single cells). Many of these alternative first exons were considerably utilized together with their canonical coun- universal cytokines (second row), Th2-specific cytokines (third row), and Th2-specific transcription factors (bottom row) among 55 single cells. The BI is indicated above the violin plots for each gene. Th2-specific genes were defined as their microarray expression value in Th2 cells being at least 4-fold higher than in either Th1 or Th17 cells. Universal cytokines were defined as being abundantly expressed in all Th subsets (microarray expression value Ͼ 16 in Th1, Th2, and Th17 cells). The microarray data of gene expression of Th1, Th2, and Th17 cells were derived from GSE11924 (40). B and C, universal cytokines preferentially exhibit a bimodal expression pattern compared with Th2-specific cytokines. B, percentages of cytokine genes with BI Ͼ 1.8 were calculated in universal cytokines or Th2-specific cytokines. C, BI of the Th2-specific or universal cytokines were plotted. The p value was calculated using a Mann-Whitney test. *, p Ͻ 0.05. Error bars stand for the standard deviation of the mean.

IL4 non-coding RNA promotes IL-4 mRNA translation
terparts (⌿ around 0.5) (Fig. 3C), indicating that utilization of alternative promoters is a general phenomenon for transcriptome diversification at the single-cell level.

IL-4 nc RNA promotes IL-4 production in Th2 cells
The first exon of the IL4 nc RNA falls into an intragenic enhancer of the Il4 locus (HSII) that regulates IL-4 but not IL-5 or IL-13 expression in T cells (25). Enhancer-derived RNAs have been widely reported and are usually associated with target gene activation (18, 26 -28), although most of their functions remain unclear. The unique expression pattern of the IL4 nc RNA among single cells prompted us to explore its function in depth.
We first developed a TaqMan real-time PCR assay that was able to distinguish the IL4 nc RNA from the pre-mRNA of the IL-4 coding isoform (IL4 fl ) by designing probes covering unique junctions between their first and the second exons, respectively (Fig. S3C). Reverse transcription using oligo(dT) primers indicated that the IL4 nc RNA was polyadenylated (Fig. S3D). Among all Th cell subsets, IL4 nc RNA was detectable in Th0 cells and Th2 cells but not in Th1, Th17, or Treg cells (Fig. 4A), consistent with reduced Il4 locus accessibility in the latter lineages (29). Fractionation further showed that IL4 nc RNA predominately localized in the cytosol rather than in the nucleus in either resting or restimulated Th2 cells (Fig. 4B). Thus, unlike most enhancer RNAs, IL4 nc RNA is a spliced and polyadenylated RNA mainly located in the cytosol.
One of the key characteristics of effector Th cells is their capacity to quickly produce a large amount of effector cytokines following antigen receptor stimulation. Prior to TCR restimulation, we found that the IL4 nc RNA was the predominant transcript from the Il4 locus, which was expressed, on average, 4-fold higher than the IL4 fl RNA (Fig. 4C). Following TCR stimulation, induction of the IL4 nc RNA (peaked 1 h post-stimulation) also preceded the IL4 fl RNA (peaked 3 h post-stimulation). Despite quicker induction, IL4 nc RNA was only induced around 7-fold after TCR restimulation, much less than the magnitude of IL4 fl RNA (Ͼ200-fold induction) ( Fig. 4D and Fig.  S3E), rendering the IL4 nc isoform constituting 2% of total IL-4 transcripts 3 h after the TCR restimulation. These data together show that IL4 nc RNA is a primary transcript from the Il4 locus in resting Th2 cells, whereas transcription shifts to IL4 fl RNA after TCR restimulation.
To examine connections between IL4 nc RNA and IL-4 production in Th2 cells, we transduced IL4 nc RNA in primary Th2 cells (Fig. 5A). IL4 nc promoted the production of IL-4 protein in Th2 cells over 40% compared with the control group, as determined by either flow cytometry (Fig. 5, B and C) or ELISA (Fig.  5D). Expression of other cytokines, including IL-13, IL-2, IL-10, and tumor necrosis factor ␣, was not significantly changed (Fig.  5, B, and C, and Fig. S3, F-H), indicating that IL4 nc RNA specifically regulates IL-4 production in Th2 cells. We further mapped sequence elements required for the IL4 nc RNA to promote IL-4 production. Only a mutant IL4 nc RNA missing the last exon, but not mutants missing either the first or the second exon, completely lost its capacity to facilitate IL-4 production (Fig. 5, E-G). Further truncation of the last exon showed that the sequence overlapping the 3Ј UTR of the coding isoform was necessary and sufficient for enhancing IL-4 production (Fig.  S4). Interestingly, despite promoting IL-4 protein expression, the noncoding isoform did not affect the IL-4 mRNA level or its expression kinetics following TCR stimulation (Fig. 5, H and I), indicating that IL4 nc directly enhances translation of IL4 fl RNA.

IL4 non-coding RNA promotes IL-4 mRNA translation
Conversely, abrogating IL4 nc RNA inhibited IL-4 protein production in Th2 cells. We designed two independent shRNAs targeting the unique region of IL4 nc RNA (its first exon) and transduced the two shRNAs into differentiating Th2 cells. Either shRNA decreased ϳ50% of IL-4 noncoding RNA without affecting the expression of IL-4 -coding mRNA (Fig. 6D), confirming the specificity of these shRNAs. Despite an unaltered IL4 fl mRNA level, production of the IL-4 protein was reduced over 40% in shRNA-transduced cells compared with control cells, as determined by either flow cytometry (Fig. 6,  A and B) or ELISA (Fig. 6C). Similarly, in murine Th2 clone D10 cells (30,31), knockdown of IL4 nc suppressed IL-4 protein production without changing IL4 fl mRNA expression (Fig. S5, A-D).
Finally, we examined the translation efficiencies of IL4 fl mRNA in IL4 nc -deficient Th2 cells using polysomal profiling assays. Compared with control cells, IL-4 mRNA in polysomal fractions was decreased in the IL4 nc -deficient cells (Fig. 6E). Conversely, overexpression of the IL4 nc transcript in Th2 cells increased the ratios of the IL-4 mRNA in polysomal fractions (Fig. 6F). Distribution of ribosome RNAs was not affected by either IL4 nc knockdown or overexpression (Fig. S6). These data together indicate that IL4 nc RNA promotes production of IL-4 by enhancing its mRNA translation, consistent with its predominant cytoplasmic localization.

Discussion
In this study, we applied scRNA-seq to characterize the transcriptome complexity of Th2 cells. Our results revealed that variable cytokine expression contributed most to the cell-tocell heterogeneity of in vitro differentiated Th2 cells. We also found that alternative promoter use was an underlying mechanism of cellular heterogeneity. Moreover, we identified that an intragenic enhancer drives expression of an IL-4 noncoding Th2 cells were differentiated in vitro as described in Fig. 1 and transduced with a retrovirus containing either IL4 nc RNA or empty vector (Mock). 5 days after stimulation, the cells were restimulated with PMA (50 ng/ml) and ionomycin (500 ng/ml) for 5 h. Expression of IL-4 was determined by intracellular staining. Data shown are representative (B) or summary (C) of IL-4 ϩ % or IL-13 ϩ % cells from three independent experiments. D, Th2 cells were differentiated and stimulated as in B, and secreted IL-4 in the supernatant was determined by ELISA . Data are summary of three independent experiments. E-G, exon 3 of IL4 nc is required for promoting IL-4 production post-transcriptionally. As above, full-length IL4 nc (FL) or its mutants depleted the indicated exons (⌬exon1, ⌬exon2, and ⌬exon3, respectively) and were retrovirally expressed in Th2 cells. Production of IL-4 was analyzed by intracellular staining (E and F) or ELISA (G). H, as in E. Expression of IL-4 and GATA3 mRNAs 2 h after restimulation was analyzed by real-time quantitative PCR. I, as in B. Th2 cells were transduced with a retrovirus overexpressing IL4 nc RNA, and the cells were selected with puromycin for 3 days (3 g/ml). 5 days after differentiation, expression of IL4 fl was analyzed by real-time PCR after activation by anti-CD3 (5 g/ml) and anti-CD28 (1 g/ml) (left panel) or PMA (50 ng/ml) and ionomycin (500 ng/ml) (right panel) for the indicated time. Data are representative (B and E) or the sum (C, D, and F-I) of at least three independent experiments. Error bars stand for the standard deviation of the mean. *, p Ͻ 0.05; **, p Ͻ 0.01; paired Student's t test. ns, not significantly changed.

IL4 non-coding RNA promotes IL-4 mRNA translation
isoform that was constitutively expressed and promoted IL-4 mRNA translation. Antigen stimulation shifts transcription of the noncoding isoform to the canonical coding IL-4 mRNA, and the noncoding isoform promotes translation of the coding IL-4 mRNA. Our study thus exemplifies that an enhancer-derived noncoding RNA promotes the translation of coding RNA from the same locus.
We chose to analyze the in vitro differentiated Th2 cells using the Fluidigm platform. Compared with recently published scRNA-seq results of in vivo differentiated Th2 cells (32,33), our analysis likely reflects intrinsic cellular heterogeneity. One of the unexpected findings is that universal cytokines, but not Th2-specific cytokines, were bimodally expressed among single Th2 cells. Bimodal distribution among single cells often suggests distinct functions or response to external stimulation, suggesting that these cytokines respond digitally to the TCR stimulation. A similar digital response has been observed in signaling pathways downstream of TCR stimulation (i.e. phos-phorylation of extracellular signal-regulated kinase) (34), suggesting that these universal cytokines are directly downstream of TCR signaling events. Other processes, including stochastic effects, the cell cycle, or the frequencies of cell division from naïve CD4ϩ T cells, may also influence the expression pattern of these universal cytokines.
Although the Fluidigm platform was only able to capture a small number of cells, it enabled analysis of full-length RNA at the single-cell level, which is unachievable with droplet-based platforms for single-cell RNA-Seq (35)(36)(37). For instance, IL4 nc would be undistinguishable from IL-4 full-length mRNA if only the 3Ј end of RNA had been sequenced. Even with a relatively small number of cells, our data revealed heterogenicity of transcripts from the IL-4 locus. Given that over 30% of Th2 cells expressed the IL4 nc transcript, our current sequence depth and cell number sufficiently reflected the cellular heterogenicity of in vitro differentiated Th2 cells. The relatively abundant IL4 nc transcript was also confirmed by our TaqMan real-time PCR assays.
Previous studies have shown that locus accessibility contributes greatly to the variation of Th2 gene expression (38, 39). Our study provides another layer of transcriptional complexity. An alternative promoter-derived noncoding RNA is constitutively expressed and able to sequester negative regulatory elements of RNA translation, which may enable rapid production of IL-4 protein right after antigen stimulation.
Our data also suggest that the long-observed lower efficiency of Th2 cell differentiation compared with Th1 is probably caused by the transcription of IL4 nc in a subset of cells. Flow cytometry and intracellular staining of IL-4 protein have been the most widely used methods to detect IL-4-expressing cells, which, however miss IL4 nc -expressing cells. Moreover, it has long been observed that naïve CD4ϩ T cells express IL-4 mRNA but not IL-4 protein after TCR stimulation, and our analysis suggests that one reason for this is that naïve CD4 ϩ T cells only express IL4 nc . It will be interesting to explore whether a similar mechanism exists for other cellular differentiation processes.

Experimental procedures
Mice C57BL/6 mice were obtained from the Shanghai Laboratory Animal Center, Chinese Academy of Sciences. All mice were maintained under specific pathogen-free conditions and used at 8 -10 weeks of age. All animal studies were performed in compliance with the Guide for the Care and Use of Laboratory Animals and approved by the Institutional Biomedical Research Ethics Committee of the Shanghai Institutes for Biological Sciences.

Cell lines
Plat-E cells were obtained from the Cell Bank of the Type Culture Collection of the Chinese Academy of Sciences (Shanghai, China) and cultured in Dulbecco's modified Eagle's medium (Hyclone) with 10% FBS (Hyclone). D10.G4.1 cells were obtained from the ATCC and cultured in RPMI 1640 medium (Hyclone) with 10% FBS, 1ϫ antibiotics, 1ϫ nonessential amino acids, 50 M ␤-mercaptoethanol, 10 ng/ml IL-1␤, Figure 6. Knockdown of IL4 nc RNA represses IL-4 production. Th2 cells were differentiated in vitro as described in Fig. 1 and transduced with a retrovirus containing a scramble sequence (sh-Ctl) or two independent shRNAs against IL4 nc RNA. Five days after differentiation, the cells were restimulated with PMA (50 ng/ml) and ionomycin (500 ng/ml) for 5 h. A-C, production of IL-4 was determined by intracellular staining (A and B) and ELISA (C). D, as in A. The cells were restimulated with anti-CD3 and anti-CD28, and expression of IL4 fl (left panel) and IL4 nc (right panel) RNA was analyzed by real-time PCR at the indicated time after restimulation. E, as in A. 5 days after differentiation, the cells were restimulated with PMA (50 ng/ml) and ionomycin (500 ng/ml) for 3.5 h. The cells were then treated with cycloheximide (100 g/ml), lysed, and fractionated on a 10 -50% sucrose gradient. Distribution of IL4 mRNA at each fraction was analyzed by real-time PCR. Distribution of ribosomal RNA is shown in Fig. S6. F, as in E. In vitro differentiated Th2 cells were transduced with a retrovirus overexpressing IL4 nc (IL4 nc ) or an empty vector (Mock). Distribution of IL4 mRNA at each fraction after sucrose density gradient centrifugation was determined by real-time PCR. Data are representative (A) or the sum (B-F) of at least three independent experiments. Error bars stand for the standard deviation of the mean. *, p Ͻ 0.05; **, p Ͻ 0.01; paired Student's t test.

IL4 non-coding RNA promotes IL-4 mRNA translation
and 100 units/ml IL-2. K562 cells were obtained from the Cell Bank of the Type Culture Collection of the Chinese Academy of Sciences and cultured in Iscove's modified Dulbecco's medium (Gibco) with 10% FBS (Hyclone), 1ϫ glutamine, 1ϫ penicillin, 1ϫ streptomycin, and 50 M ␤-mercaptoethanol.

Single-cell and bulk RNA-Seq
Single cells were isolated with a Fluidigm C1 system. A total of 56 cells were collected and sequenced with an Illumina Hiseq system. For bulk RNA-Seq, 10,000 pooled cells were prepared using a standard protocol for the RNA-Seq library. All fastq files were deposited in the GEO database under accession number GSE103647. The microarray data of gene expression of Th1, Th2, and Th17 cells was derived from GSE11924.

Bimodal index analysis
The BI of gene expression was calculated as described previously (23). Briefly, if a gene is expected to have two populations in the RNA expression level by the model, then the mean expression values of the two populations ( 1 and 2 ), the common standard deviation (), and the percentage of one subset () are predicted, and then the standardized distance between two subsets(␦) and BI are defined as

RNA extraction, cDNA synthesis, and quantitative real-time PCR
Total RNA was isolated with the E.N.Z.A HP Total RNA Kit (Omega), and cDNA was synthesized using the PrimeScript TM RT Reagent Kit with gDNA Eraser (Takara). Quantitative PCR was performed using SYBR Premix Ex Taq (Takara) on a StepOnePlus real-time PCR system (Applied Biosystems). Gene expression was normalized to glyceraldehyde-3-phosphate dehydrogenase.

Intracellular cytokine staining and flow cytometry
Cells were stimulated with phorbol 12-myristate 13-acetate (PMA, Sigma, 50 ng/ml) and ionomycin (Sigma, 500 ng/ml) in the presence of Golgi-plug (BD Biosciences) for 5 h at 37°C. Cells were stained with the LIVE/DEAD Fixable Violet Dead Cell Stain Kit (Thermo Fisher), fixed/permeabilized with the Cytofix/Cytoperm Kit (BD Biosciences), and then stained with antibodies against the indicated cytokines. Fluorescence-labeled anti-IL-4 and anti-IL-13 antibodies were from eBioscience.

Retroviral infection
Retrovirus was produced by transfecting Plat-E cells with the indicated plasmids, and supernatant containing fresh virus was collected 48 h and 72 h after transfection. To infect K562 cells and D10.G4.1 cells, cells were suspended with fresh retrovirus, spin-infected for 1.5 h at 900 ϫ g at 37°C in the presence of Polybrene (8 g/ml), and cultured at 37°C for another 6 h before recovering to culture medium. Then puromycin (3 g/ml) was added for at least 3 days before analysis. To infect T H 2 cells, naive CD4ϩ T cells were first stimulated with anti-CD3, anti-CD28, anti-IFN-␥, and IL-4 as mentioned above. For the first infection, 80% of the stimulation medium was replaced with medium containing fresh retrovirus 18 h after stimulation. Then the CD4ϩ T cells were spin-infected for 1.5 h at 900 ϫ g at 37°C in the presence of Polybrene (8 g/ml) and cultured at 37°C for another 2 h before recovering to the stimulation medium. The second infection was carried out 40 h after stimulation the same as the first infection.

Plasmids
For knockdown of IL-4 noncoding RNA, shRNA sequences targeting mouse IL-4 noncoding RNA or a scramble sequence of the same length was inserted into the MSCV-LTRmiR30-PIG (LMP) vector. The LMP and pcl-eco plasmids were cotransfected into Plat-E to produce retrovirus. For the purpose of overexpression, the IL-4 coding sequence and IL-4 noncoding RNA and its truncated isoforms were inserted into the MSCV-PIG vector.

Polysomal profiling
After cycloheximide treatment (100 g/ml), 2.5 ϫ 107 Th2 cells were collected and lysed in 1.3 ml of lysis buffer (100 mM KCl, 5 mM MgCl 2 , 20 mM Tris-HCl (pH 7.5), 0.03% TritonX-100, 1 mM dithiothreitol, 40 units/ml RNase inhibitor (Takara), and protease inhibitor (Roche)) on ice for 30 min. Then cells were spun for 10 min at 12,000 rpm (13,800 ϫ g) at 4°C, and supernatant was collected. 10%, 20%, 30%, 40%, and 50% sucrose buffers were prepared in lysis buffer. To obtain a linear IL4 non-coding RNA promotes IL-4 mRNA translation gradient, 1.5 ml of each sucrose buffer was added into layers with the heavy sucrose solution at the bottom, and the gradient was preserved at 4°C overnight. The cell lysate was added to the top of the gradient. After ultracentrifugation (35,000 rpm for 2 h in a SW40Ti rotor at 4°C), 19 fractions from the sucrose gradient were collected with the same volume. The RNA extraction, cDNA synthesis, and quantitative real-time PCR were conducted as mentioned above.
Author contributions-W. Y. formal analysis; W. Y. and X. C. investigation; W. Y. and X. C. writing-original draft; W. Y. and X. C. project administration; Y. S. data curation; X. C. conceptualization; X. C. supervision; X. C. funding acquisition; X. C. writing-review and editing.