RNA-seq Analysis of δ9-Tetrahydrocannabinol-treated T Cells Reveals Altered Gene Expression Profiles That Regulate Immune Response and Cell Proliferation*

Marijuana has drawn significant public attention and concern both for its medicinal and recreational use. Δ9-Tetrahydrocannabinol (THC), which is the main bioactive component in marijuana, has also been shown to possess potent anti-inflammatory properties by virtue of its ability to activate cannabinoid receptor-2 (CB-2) expressed on immune cells. In this study, we used RNA-seq to quantify the transcriptomes and transcript variants that are differentially regulated by THC in super antigen-activated lymph node cells and CD4+ T cells. We found that the expressions of many transcripts were altered by THC in both total lymph node cells and CD4+ T cells. Furthermore, the abundance of many miRNA precursors and long non-coding RNAs was dramatically altered in THC-treated mice. For example, the expression of miR-17/92 cluster and miR-374b/421 cluster was down-regulated by THC. On the other hand miR-146a, which has been shown to induce apoptosis, was up-regulated by THC. Long non-coding RNAs that are expressed from the opposite strand of CD27 and Appbp2 were induced by THC. In addition, THC treatment also caused alternative promoter usage and splicing. The functions of those altered transcripts were mainly related to immune response and cell proliferation.

Marijuana (Cannabis sativa) is the most widely used drug of abuse in the United States. Nonetheless, it is worth noting that its medicinal value has been known for many centuries and is getting increasing recognition in recent times. Thus, ⌬9-tetrahydrocannabinol (THC), 2 the major psychoactive principle in marijuana, is used for medicinal purposes such as for pain relief as well as an anti-emetic and appetite stimulant during cancer chemotherapy and in HIV/AIDS patients (1). In addition, marijuana is widely known anecdotally for its use to relieve pain, and prevent seizures in patients with epilepsy. Because of the potential therapeutic value of marijuana in a variety of clinical disorders, 23 states and the District of Columbia currently have laws legalizing marijuana for medicinal use. Four states and the District of Columbia have also legalized marijuana for recreational use. Clearly, legalization of marijuana for medicinal or recreational purposes raises significant questions on its long term use/abuse and the impact on human health.
THC has also been shown to modulate immune response primarily through activation of cannabinoid receptors (CB2) expressed primarily on cells of the immune system (2)(3)(4). In general, most studies have shown that THC has immunosuppressive properties mediated through multiple pathways (2,5,6). For example, it has been shown that THC suppresses Th1, whereas promoting Th2 cells (7). THC has also been shown to induce apoptosis in immune cells (8,9), and induce myeloidderived suppressor cells and Tregs, which are known to inhibit T cell proliferation (10 -12). Moreover, induction of myeloidderived suppressor cell by THC is associated with altered microRNA expression (13). Our previous studies have shown that epigenetic modifications such as histone methylation and acetylation play important roles in THC-mediated gene expression in immune cells (7).
In this study, we used RNA-seq to examine global gene expression in antigen-activated lymphocytes. RNA-seq is a next generation sequencing-based RNA quantification method (14). In this case, a population of RNA is isolated and converted to a library containing short cDNA fragments. Fragments with adapters to one or both ends are then sequenced using highthroughput sequencing technology. The fragments are then aligned to the genome to calculate the abundance of the transcript. Although RNA-seq is still under active development, it has shown many advantages over microarray (14). RNA-seq has a very low background signal because fragments can be unambiguously mapped to regions of the genome and unmapped sequences are filtered out. Because it is sequencing based, RNA-seq does not have an up or low limit for quantification. Therefore, it has a high sensitivity for genes expressed at very low or very high levels. The single-base resolution of RNA-seq can also detect the transcript boundaries that may lead to the discovery of alternative transcription start sites and splicing isoforms. Furthermore, novel transcripts including non-coding RNAs can be identified by RNA-seq.
Staphylococcal enterotoxin B (SEB), is a superantigen produced by Staphylococcus aureus and is responsible for causing food poisoning (15) and toxic shock (16). SEB binds to the nonpolymorphic regions of MHC II on antigen presenting cells and the specific V␤ regions of the T-cell receptor (TCR) such as murine V␤8, thereby activating a large population (up to 30%) of T cells that proliferate rapidly and trigger cytokine storm (17). Recent studies from our laboratory demonstrated that THC can attenuate SEB-mediated inflammation and toxicity to the lungs (18 -20). In this study we wanted to examine the effect of THC on global gene expression in activated lymphocytes. We used the same animal model and SEB was used as the reagent to induce inflammation. We then used RNA-seq to analyze the lymph node (LN) cells as well as purified CD4 ϩ T cells from mice immunized with SEB and treated with vehicle or THC to analyze RNA abundance. We found that the expression of many transcripts was altered by THC in both total lymph node cells and CD4 ϩ T cells. Functional analysis revealed that those altered genes were mainly related to immune response and cell proliferation. Furthermore, the abundance of many miRNA precursors and long noncoding RNA (lncRNA) was dramatically altered in THC-treated mice, suggesting that non-coding RNA expression might play an important role in THCmediated immune regulation.

THC-mediated Suppression of Lymphocyte Proliferation and
Inflammation-SEB is a potent superantigen that activates ϳ30% of T cells particularly those expressing V␤8 T cell receptor. Previous studies from our laboratory demonstrated that THC can attenuate SEB-induced inflammation (7,21). However, to ensure that this model was working in the immune cells under current investigation, we used SEB to induce inflammation and tested the effect of THC. To that end, we immunized mice with SEB into the footpads so that SEB-specific T cell response can be studied in the draining popliteal LN as described (22,23). To make sure that the THC effects were dependent on CB2 receptors, another group of mice received both THC and CB2 antagonist, SR144528 (24). The data shown FIGURE 1. Comparison of gene expression in vehicle-and THC-treated cells. C57BL/6J mice were treated with THC or vehicle as described under "Experimental Procedures." Two hours after the second THC treatment, 10 g of SEB was injected into each footpad. Three days after SEB challenge, draining popliteal lymph node cells (LN) of vehicle-or THC-treated mice (n ϭ 5) were harvested. CD4 ϩ cells were purified from the total LN cells. RNA was purified from those cells and the expression of each gene was quantified by RNA-seq. The expression level in total LN cells (a) and CD4 ϩ cells (b) under the vehicle and THC treatment is presented as Fragments Per Kilobase per Million mapped reads (FPKM). The difference in gene expression is also presented as fold-change (c and d). The FPKM value of each gene in THC-treated LN cells (c) and CD4 ϩ cells (d) was divided by its value in the vehicle-treated sample and expressed as log 2 .
in Fig. 1 demonstrated that THC treatment in SEB-immunized mice caused a significant decrease in total draining LN cellularity, CD4 ϩ T cells as well as Ifn-␥ ϩ CD4 ϩ T cells when compared with vehicle controls (supplemental Fig. S1). Furthermore, in mice that received THC ϩ SR144528, the immunosuppressive effects of THC were reversed (supplemental Fig. S1), thereby indicating that the THC effects in lymphocytes were mediated through CB2. These results were consistent with our previous report that THC suppresses cell proliferation and induces a shift in T cells from Th1 to Th2 (7).
THC-induced Differential Expression of Genes in Lymphocytes-To investigate the effect of THC treatment on SEBinduced gene expression, the abundance of each transcript was compared by RNA-seq between the SEB ϩ vehicle-and SEB ϩ THC-treated total lymph node cells as well as purified CD4 ϩ cells. The rationale for using both purified CD4 ϩ T cells and whole LN cells was as follows. 1) THC is well known to exert its effect on Th cell differentiation and therefore, studying CD4 ϩ T cells would provide useful clues. 2) The cytokine storm induced by SEB can also activate other immune cells that can be studied using whole LN cells. 3) Our previous ChIP-seq study was performed in the total lymph node cells (7). Thus, comparing gene expression profile and histone modification profile would allow us to determine whether THC-induced alteration in gene expression is mediated by histone modification. We obtained 30 to 40 million reads in each sample, and 70 -80% of reads were mapped to the mouse genome. In those mapped reads, 75-80% were uniquely mapped. Using Cufflinks and Cuffcompare, the expression level of each gene was calculated. Although most genes were expressed at a low level in both SEB ϩ vehicleand SEB ϩ THC-treated samples, there were many genes that were differentially expressed under these two treatments ( Fig.  1). In total LN cells, several genes were undetectable in SEB ϩ THC-treated sample, whereas highly expressed in the SEB ϩ vehicle-treated sample (Fig. 1a), thereby suggesting that THC was down-regulating the expression of genes induced by SEB. The differentially expressed genes were also presented as foldchanges ( Fig. 1, c and d).
Genes that showed more than 2-fold difference and with an "OK" status as determined by Cuffcompare were further analyzed. Compared with SEB ϩ THC group, 310 and 260 genes were highly expressed in the SEB ϩ vehicle-treated total LN cells and CD4 ϩ cells, respectively. Among those, 37 were common in total LN and CD4 ϩ cells (Fig. 2a). On the other hand, the expression of 328 genes in the total LN cells and 333 genes in CD4 ϩ cells were induced by THC ϩ SEB treatment. Among them, 56 were common in two cell populations (Fig. 2b). A list of those differentially expressed genes is presented in supplemental Data S1-S3. Interestingly, among those differentially expressed genes, many of them were miRNA precursors and lncRNAs. There were 122 and 128 such transcripts that were highly expressed in the SEB ϩ vehicle-treated total LN and CD4 ϩ cells, respectively. On the other hand, there were 111 and 100 non-coding RNAs that were highly expressed in SEB ϩ THC-treated LN and CD4 ϩ cells, respectively. However, few of them were common in total LN and CD4 ϩ cells (Fig. 2, c and d).
The majority of these miRNAs had unknown functions (see supplemental Data S1-S3 for lists of these miRNAs).

miRNA Precursors with Altered Expression after THC
Treatment-Because RNA-seq can only identify pri-miRNA, it was unclear whether the pri-miRNA level correlates with its mature miRNA expression. Suppression of miR-17/92 cluster by THC has been shown in another animal model (21). Indeed, the expression of precursor of this miRNA cluster was downregulated by THC in this study. We confirmed the down-regulation of mature miRNAs of this cluster by real-time PCR (data not shown). This RNA-seq study also identified another miRNA cluster, miR-374b/421, whose precursor was also down-regulated by THC in total LN cells. pri-miR-421 can be processed into 2 mature miRNAs, miR-421-1 and miR-421-2. They are located at the 5Ј and 3Ј of pri-miR-421, respectively. Real-time PCR quantification of mature miR-421-1 and miR-421-2 confirmed that their expression levels were lower in THC-treated samples (Fig. 3a). In our previous ChIP-seq study, we have found that the promoter of certain miRNAs is differentially associated with histone marks in the THC-treated lymphocytes (7). In the promoter of miR-374b/421 cluster, histone H3 trimethylation at Lys-4 and Lys-36 (H3K4me3 and H3K36me3), which are associated with transcription activation, were present in both SEB ϩ vehicle and SEB ϩ THC samples. However, histone H3 trimethylaiton at Lys-27 (H3K27me3), which is associated with transcription suppression, was only present in the promoter area of this miRNA cluster in SEB ϩ THC-treated sample, suggesting a possible histone methylation-mediated suppression (Fig. 3b). We also examined another miRNA, miR-146a. The expression of its precursor was induced in the LN cells by THC based on the RNA-seq data. Real-time PCR indeed showed that the mature miR-146a was significantly induced by THC (Fig. 3c). Again, we examined the association of histone marks under those two conditions. In this case, the association of H3K4me3, H3K4me27, and H3K36me3 did not differ significantly in the SEB ϩ vehicle and SEB ϩ THC-treated groups. However, his-FIGURE 2. THC-induced differentially expressed genes. Mice were treated with SEB and vehicle or THC, and LN cells or CD4 ϩ T cells were analyzed as described in the legend to Fig. 1. Genes with at least 2-fold differences in expression between the vehicle and THC-treated cells were identified. The Venn diagrams present the numbers of all genes (a), and miRNA precursors (c), whose expression was high in the vehicle-treated cells but suppressed by THC in the LN cells and CD4 ϩ cells, and the numbers of all genes (b), and miRNA precursors (d), whose expression was high in THC-treated cells but low in vehicle-treated cells.
tone H3 acetylation at Lys-9 (H3K9Ac) was only present in the promoter area of miR-146a in the THC-treated sample, suggesting histone acetylation might play a role in THC-induced expression of miR-146a (Fig. 3d). Fig. 4 is a Circos plot that shows genome-wide histone methylation, gene expression, and links between differentially expressed miRNAs and their target mRNAs. Two relatively well studied histone marks, H3K4me3 and H3K27me3, in total lymphocytes are presented as histograms. The differences in gene expression in both total lymphocytes and CD4 ϩ cells are presented as heat maps. Genes in red and brown were those induced by THC and genes in blue and green are those suppressed by THC. Genes in yellow did not show significant difference. The links show the genomic positions of the differentially expressed miRNAs and their target mRNAs that were also differentially expressed. Cold colors (green, blue etc.) represent miRNAs that were down-regulated by THC, and hot colors (red, brown etc.) represent THC-induced miRNAs. A list of those miRNAs and their potential target mRNA genes are in supplemental Data S1-S3.
Pathway Analysis of Differentially Expressed Genes-We further performed gene interaction analysis by Ingenuity pathway analysis on those differentially expressed genes. Most of those known miRNAs were involved in cell proliferation (Fig. 5). This was consistent with our previous study showing that THC promotes apoptosis in antigen-activated lymphocytes (7). For the protein-coding genes with more than 2-fold changes, THC decreased the expression of 188 genes in the total LN cells and 132 genes in CD4 ϩ T cells. Among those, there were 20 genes that overlapped (Fig. 6a). On the other hand, 217 and 233 genes were up-regulated by THC in total LN and CD4 ϩ cells, respectively, of which 39 genes overlapped (Fig. 6b). Functional analysis of those differentially expressed genes indicated that the major enriched pathway was inflammatory response (Fig. 6, c and d).
THC-induced Alteration in the Expression of Long Non-coding RNA-Besides alteration in miRNA and protein coding gene expression, THC also modulated the expression of several lncRNA (supplemental Data S1-S3). The functions of those lncRNAs are not yet clear. We selected 3 lncRNAs whose expression was induced by THC based on RNA-seq data. These lncRNAs are transcribed from the opposite strand of known genes, and may have a suppressive effect on those protein coding genes. 5730480H06Rik overlaps the 3Ј end of Kv channel interacting protein 4 (KCNIP4). Appbp2OS overlaps the transcription start site of amyloid ␤ precursor protein-binding protein 2. E130112N10Rik overlaps the whole transcript of Cd27 (Fig. 7a). To confirm that those non-coding RNAs were induced by THC through CB2 receptor, we examined their expression in an independent experiment by real-time PCR. First, we examined the marks of Th1 and Th2, which are known to be regulated by THC. Tbx-21, the transcription factor for Ifn-␥, and Il-2, the pro-inflammatory cytokine were suppressed, whereas Il-4, the anti-inflammatory cytokine was induced by THC (Fig. 7b). In the same set of samples, the expression of these non-coding RNAs was induced by THC, and the induction was reversed by SR144528 (Fig. 7c).
Transcripts Using Different Promoters or Having Different Splicing Forms after THC Treatment-RNA-seq not only can identify differentially expressed genes but also transcripts with different promoter usage and splicing form. Approximately 600 transcripts in both total LN cells and CD4 ϩ T cells used different promoters following SEB ϩ vehicle and SEB ϩ THC treatments. About half of those transcripts were common in both total LN and CD4 ϩ T cells (Fig. 8a). As for alternative splicing, THC induced about 2500 different splicing and nearly half of them were common in both types of cells (Fig. 8b). Functional analysis of those transcripts with different promoters or splicing forms showed that those transcripts were involved in vari-ous pathways, whereas cell death and survival was the top network (Tables 1 and 2 and supplemental Data S1-S3).

Discussion
Previous studies from our laboratory and elsewhere have shown that THC suppresses the expression of pro-inflammatory genes while inducing the expression of anti-inflammatory genes (6,7,25). To identify other genes whose expression might be regulated by THC, we used RNA-seq to compare the abundance in SEB-activated LN cells or CD4 ϩ T cells, with and without THC treatment. Some genes such as Il-4 identified in this study have already been shown in our previous studies and those of others (7,21). Interestingly, many genes such as Fpr2, which has been shown to have an anti-inflammatory role (26), but has never been shown previously to be induced by THC were also identified. However, the fact that only a small number of genes were either up-regulated or down-regulated by THC in total popliteal lymph node cells and CD4 ϩ cells, suggested that FIGURE 5. Top network of THC-regulated miRNAs. Mice were treated with SEB and vehicle or THC, and LN cells or CD4 ϩ T cells were analyzed as described in the legend to Fig. 1. Ingenuity pathway analysis indicated that cell death and survival was the most overrepresented network by miRNAs whose precursors were differentially expressed in the vehicle and THC-treated lymphocytes. THC up-regulated miRNAs are shown in red and THC down-regulated miRNAs are shown in green. JULY 22, 2016 • VOLUME 291 • NUMBER 30 THC may differentially regulate gene expression in different immune cells. This can be explained by the fact that although SEB activates primarily V␤8 ϩ T cells, the cytokine storm induced by this superantigen can also affect other immune cells. It is possible that THC may even have an opposite effect on the expression of certain genes in different subset of cells. It is, however, not clear whether this is the case. In the future, we will separate different types of cells from lymph nodes and compare their gene expression profiles. Despite few overlapping genes between total lymphocytes and CD4 ϩ cells, most of those altered were related to cell survival and inflammatory response. This is consistent with our flow cytometry result using the same animal model in our previous study (7), which showed that THC promotes apoptosis in SEB-activated lymphocytes, and inhibits production of Th1 cytokines while inducing the production of Th2 cytokines (7).

RNA-seq Analysis of THC-treated T Cells
A unique finding of this study is that the expression of many miRNA precursors was dramatically affected by THC. The majority of these miRNAs are newly identified with unknown functions. These data suggested that THC-mediated alteration in miRNA expression may play a critical role in cell prolifera-tion and immune response. These novel miRNAs may have important roles in immune response. In the future, we will first validate the expression of these miRNAs in their mature miRNA form and then investigate their potential roles in the regulation of immune functions. Among those known miRNAs, the miR-17/92 cluster has been shown to be induced by SEB in an acute lung injury mouse model by regulating PIK3/Akt pathway (21). Because THC suppresses the expression of the miR-17/92 cluster in different tissues (21), it is likely that this effect is not cell type specific. This miRNA cluster has been shown to have many functions and THC may exert a broad spectrum of functions through this miRNA cluster. However, it is not clear how THC modulates the expression of miR-17/92 cluster. We examined the association of 5 histone modification marks (H3K4me3, H3K27me3, H3K36me3, H3K9me3, and H3K9ac) in its gene. No significant difference was found between the vehicle and THC-treated samples (data not shown), suggesting that these histone marks may not play a key role in THC-mediated modulation of the expression of miR-17/92 cluster.
The miRNA cluster 23, including miR-374b to miR-421, was also down-regulated by THC. In human cells, miR-421 sup-FIGURE 6. Protein-coding genes that are differentially regulated by THC. Mice were treated with SEB and vehicle or THC, and LN cells or CD4 ϩ T cells were analyzed as described in the legend to Fig. 1. Numbers of protein-coding genes whose expression was inhibited by THC (a) or induced by THC (b) in total LN cells and CD4 ϩ cells are presented by Venn diagrams. Their differences under these two treatments were at least 2-fold. Functional analysis indicated that these differentially expressed genes, in both LN cells (c) and CD4 ϩ cells (d), are involved in inflammatory response. THC-induced genes are in red, whereas THC suppressed genes are in green.
presses the expression of Atm, a key kinase in the maintenance of genomic integrity, by targeting the 3Ј UTR of the Atm transcript, and the expression of miR-421 is induced by oncogene N-Myc (27). This raises the possibility that THC may regulate Atm-dependent DNA damage response. miR-421 has been shown to be overexpressed in various types of cancer (28 -31). A recent study suggested that miR-421 induces cell proliferation and apoptosis resistance by targeting Foxo4 (32). Down-regulation of miR-421 by THC is consistent with the notion that THC has a pro-apoptotic property. The miR-374b/421 cluster is located in the Ftx gene, which encodes a lncRNA. Another miRNA cluster, miR-545/374a, is also located in this lncRNA. This cluster is overexpressed in hepatitis B virus-induced hepatocellular carcinoma (33). It is unclear whether THC differentially regulates the expression of individual miRNAs in these clusters.
The expression of miR-146a is highly induced by THC at both the pri-miRNA level and mature miRNA level. This miRNA has been relatively well studied. Many studies indicated that miR-146a induces cell apoptosis and its abnormal expression has been found in various types of cancers (34 -37). miR-146a also plays an important role in immune response. A study in miR-146a knockout mice has shown that these mice produce an excessive amount of pro-inflammatory cytokines after bacterial challenge (38). A recent study showed that the expression of miR-146a and miR-146b is increased when monocytes differentiate into dendritic cells (39). Silencing of these two miRNAs enhances the expression of pro-inflammatory cytokines Il-12, Il-6, and Ifn-␥ (39). In our model, THC-induced expression of miR-146a correlates with the decreased expression of these cytokines. We also find that several histone markers are associated with the miR-146a gene, suggesting that the expression of this miRNA may be regulated by the environmental stimuli through epigenetic mechanisms, which is supported by the presence of H3K9ac in the promoter area of its genes in THC-treated cells (Fig. 3d).
In this study, we also identified several THC-regulated lncRNAs. The functions of most lncRNAs are not known and they are involved in a variety of biological processes through different mechanisms. lncRNAs can serve as signals for transcription, guides for recruiting chromatin-modifying enzymes, or scaffolds for protein complex assembly (40). The biggest challenge of studying lncRNA functions is that it is very difficult to predict by which mechanism that a particular lncRNA may exert its regulatory function. A large proportion of lncRNAs is antisense to protein-coding transcripts. Antisense RNAs can suppress the expression of their sense counterparts through inhibiting transcription and translation (41). Here we show 3

RNA-seq Analysis of THC-treated T Cells
examples of antisense lncRNA that overlap their sense proteincoding transcripts in different positions. Whether and how these lncRNAs suppress the expression of their counterparts needs further investigation. The purpose of this study is to screen for potentially interesting THC-regulated lncRNAs. KCNIP4 has been identified as an asthma gene (42). Marijuana has been linked to the cause as well as treatment for many respiratory diseases including asthma (9,43).
Another interesting lncRNA is Appbp2os. Appbp2 is a microtubule-interacting protein and is involved in the transporting of amyloid precursor protein. It is not known whether Appbp2os regulates the expression or function of Appbp2 protein. However, THC can inhibit amyloid ␤-peptide aggregation (44). Cannabis-based medicine has been shown to reduce multiple pathological processes in an animal model of Alzheimer disease (45). Recently, CB2 receptor has been linked to amyloid-␤ processing and shows a minor therapeutic property in a mouse model (46). How THC affects the pathology of Alzheimer disease is worth further investigation.
THC-induced lncRNA E130112N10Rik is much longer than its counterpart sense transcript, Cd27, and covers the whole transcript of Cd27. CD27 is required for immune cells to maintain long-term antigen-specific T cell immunity (47). A study also showed that Tregs inhibit CD27-dependent Th1 priming (48). So far there are no reports regarding THC regulating CD27-mediated immune response. However, the role of THC in inhibiting Th1 is well documented. This lncRNA could add another possible mechanism of THC-induced immune suppression.
Pathway analysis of THC-induced differentially expressed genes, both protein-coding and nonprotein-coding genes, has indicated that THC has a pleotropic effect, from cell proliferation to metabolism. Therefore, it is necessary to investigate the role of THC in modulating various biological processes to weigh the benefits and risks of marijuana use. This study also shows that THC induces alternative promoter usage and splicing. More than half of human genes contain multiple promoters and have different splicing forms (49). Alternative transcription start sites can yield different pre-mRNAs. Alternative splicing of these pre-mRNAs leads to multiple transcript variants. As a result, one gene can produce various gene products under different conditions. Increasing evidence indicates that many diseases are associated with altered transcript isoforms. For example, it has been shown that using alternative transcription start and termination sites is a major source of transcriptome diversity in developing mouse cerebella (50). Moreover, that study indicates that H3K4me3 defines the use of alternative promot- Functional analysis of genes with differential promoter usage ers that links a histone mark to promoter usage. The significance of THC-induced alternative promoter usage and splicing in immune cells is currently not clear.
In summary, our study has identified many protein coding and non-coding RNA transcripts whose expression is significantly altered by THC. This may provide information for future studies regarding the role of THC in various biological processes. It would also be interesting to investigate how histone modification, non-coding RNA, and alternative splicing coordinate each other for THC to exert its effects.

Experimental Procedures
Mice and Cell Isolation-Female C57BL/6J mice (6 -7 weeks old) received intraperitoneal (i.p.) injection of THC (Sigma, 20 g/kg of body weight) or the same amount of vehicle as described previously (51). The mice received the same treat-  JULY 22, 2016 • VOLUME 291 • NUMBER 30 ment 24 h later. Two hours after the second treatment, 10 g of SEB was injected in each hind footpads. To study the effect of CB2 antagonist, mice received SR144528 (Santa Cruz Biotechnology, 10 g/kg of body weight) by i.p. route, at the same time as THC. Popliteal draining lymph nodes from the vehicle and THC-treated mice were collected 72 h after SEB injection and single cells were prepared using micro-Biomaster (Seward Lab Systems, West Sussex, UK). Samples from 4 -7 mice in each group were pooled. Cells were filtered through a 70-m mesh and counted. CD4 ϩ cells were isolated using EasySep PE selection kit (StemCell Technologies, Vancouver, BC) according to the provided instruction. The antibody used was PE-antimouse CD4 (BioLegend, San Diego, CA).

RNA-seq Analysis of THC-treated T Cells
Staining and FACS Analysis-For CD3 and CD4 staining, lymph node cells were stained with PE/Cy7-CD3 and PE-CD4 immediately after harvest. For intracellular staining, lymph node cells were cultured in complete RPMI in the presence of 1 nM PMA (Sigma), 1 M calcium ionophore (Sigma), and 2 M protein transport inhibitor monensin (Biolegend, San Diego, CA) for 4 h. Cells were washed and resuspended in FACS buffer (PBS containing 2% FBS and 0.1% sodium azide). Fc receptors were blocked by adding anti-mouse CD16/CD32 (10 g/ml) followed by surface staining for PE-CD4 and intracellular staining with FITC-Ifn-␥. All antibodies were purchased from BD Biosciences. Cells were analyzed in a BC FC 500 flow cytometer.
RNA-seq-RNA-seq libraries were constructed using Illumina TruSeq RNA Sample Preparation kit. Briefly, total RNA was purified using Qiagen RNA easy kit. The oligo(dT) beads were added to 1 g of total RNA and the purified mRNA was fragmented. The RNA fragments were then reverse transcribed into double-stranded cDNA fragments. The DNA fragments were repaired to generate blunt ends using T4 DNA polymerase, Klenow polymerase, and T4 polynucleotide kinase. After DNA fragments were purified using a Qiagen PCR purification kit (Qiagen catalogue number 28004), an "A" base was added to the 3Ј end of the blunt DNA fragment by Klenow fragment. Sequencing adapters were ligated to the ends of DNA fragments using DNA ligase. The libraries were then amplified by limited PCR (15 cycles) using primers provided by the kit. The PCR products were then separated in 2% agarose gel and fragments with size range from 250 to 400 bp were excised and purified using QIAquick Gel Extraction Kit (Qiagen number 28704). The concentration and distribution of the library were determined by a NanoDrop spectrophotometer (Thermo Scientific, Wilmington, DE). The library was sequenced by illumina HiSeq 2000 at Tufts University Genomic core facility. Raw sequencing reads (50 bp single-end) were mapped to mouse genome build mm9 using Tophat 2 (52). The accepted hits were used for assembling transcripts and estimating their abundances using Cufflinks (53). The differentially expressed gene, promoter usage, and splicing form were determined by Cuffdiff and Cuffcompare (53).
Quantitative Real-time PCR-Mature miRNA was quantified by real-time PCR using miScript SYBR Green PCR kit and miScript Primer Assays (Qiagen). Quantitative real-time PCR was performed according to the protocol provided by the manufacturer. Snord96a and Gapdh were used as the internal con-trols for miRNAs and mRNA, respectively. The amount in the vehicle-treated sample was set as 1.
Functional Enrichment-Ingenuity pathway analysis (Qiagen) was used for functional analysis. Ingenuity pathway analysis has tools to identify the canonical pathways and top networks in which a given set of genes are involved. We performed functional enrichment analysis with the transcriptome data by taking all the transcripts with a fold-change of 2 or more. Tables  1 and 2 were generated by such analysis. The software can also identify the interactions between miRNAs and target genes. Genes and the miRNAs that had more than 2-fold changes were analyzed together to see whether any of these genes are the targets of the dysregulated miRNAs. These interactions were presented as links in the Circos plot (54). The histogram of histone methylation marks and the heat map of gene expression were also generated by the Circos software.
Author Contributions-X. Y. and M. B. designed the study, performed the experiments, and wrote the paper. P. S. N. and M. N. conceived and coordinated the study, interpreted the results, and wrote the paper. All authors reviewed the results and approved the final version of the manuscript.