The lncRNA Malat1 inhibits miR-15/16 to enhance cytotoxic T cell activation and memory cell formation
Abstract
Proper activation of cytotoxic T cells via the T cell receptor and the costimulatory receptor CD28 is essential for adaptive immunity against viruses, intracellular bacteria, and cancers. Through biochemical analysis of RNA:protein interactions, we uncovered a non-coding RNA circuit regulating activation and differentiation of cytotoxic T cells composed of the long non-coding RNA Malat1 (Metastasis Associated Lung Adenocarcinoma Transcript 1) and the microRNA family miR-15/16. miR-15/16 is a widely and highly expressed tumor suppressor miRNA family important for cell proliferation and survival. miR-15/16 play important roles in T cell responses to viral infection, including the regulation of antigen-specific T cell expansion and memory. Comparative Argonaute-2 high-throughput sequencing of crosslinking immunoprecipitation (AHC) combined with gene expression profiling in normal and miR-15/16-deficient mouse T cells revealed a large network of hundreds of direct miR-15/16 target mRNAs, many with functional relevance for T cell activation, survival and memory formation. Among these targets, Malat1 contained the largest absolute magnitude miR-15/16-dependent AHC peak. This binding site was among the strongest lncRNA:miRNA interactions detected in the T cell transcriptome. We used CRISPR targeting with homology directed repair to generate mice with a 5-nucleotide mutation in the miR-15/16-binding site in Malat1. This mutation interrupted Malat1:miR-15/16 interaction, and enhanced the repression of other miR-15/16 target genes, including CD28. Interrupting Malat1 interaction with miR-15/16 decreased cytotoxic T cell activation, including the expression of interleukin 2 (IL-2) and a broader CD28-responsive gene program. Accordingly, Malat1 mutation diminished memory cell persistence in mice following LCMV Armstrong and Listeria monocytogenes infection. This study marks a significant advance in the study of long non-coding RNAs in the immune system by ascribing cell-intrinsic, sequence-specific in vivo function to Malat1. These findings have implications for T cell-mediated autoimmune diseases, antiviral and anti-tumor immunity, as well as lung adenocarcinoma and other malignancies where Malat1 is overexpressed.
eLife assessment
This is an important study that revealed a new noncoding RNA regulatory circuit involved in T cell function. The authors provide compelling evidence, that is more rigorous than the state-of-the-art, using genetically engineered mice and cell-based experiments. The interpretation of the results should be tempered due to the small effect size observed.
https://doi.org/10.7554/eLife.87900.3.sa0Introduction
Cytotoxic T cells are indispensable for mounting an adaptive immune response against intracellular pathogens and the clearance of mutated cells such as cancer. Over the course of a viral infection, cytotoxic T cells are primed by antigen presenting cells and undergo extensive rounds of intense proliferation (Murali-Krishna et al., 1998). As they clonally expand, these cells differentiate into effector and memory cells and acquire effector functions, including the production of cytotoxins and critical cytokines such as IL-2, tumor necrosis factor alpha (TNFɑ), and interferon gamma (IFNγ) (Bachmann et al., 1999). As the infection is cleared, the cytotoxic T cell response contracts as many cells die via apoptosis, yielding a long lived pool of memory cells poised for secondary expansion and protection against reinfection with the same pathogen (Badovinac et al., 2002; Wherry and Ahmed, 2004). The factors that control this expansion, differentiation, and contraction have been intensely researched in the past decades, focusing in large part on proteins such as transcription factors, signaling enzymes, and cytokines (Chen et al., 2018). However, it is of recent interest to understand how non-protein-coding regions of the genome contribute to the regulation of these cells as well. These non-protein-coding elements can be regulatory in nature such as enhancers (Roychoudhuri et al., 2016; Shapiro et al., 1997), or RNA species that are transcribed but not translated. Two such species of interest to our present study are microRNAs (miRNAs) and long non-coding RNAs (lncRNAs).
MicroRNAs are short (21 nucleotide) RNAs which, when loaded into Argonaute (Ago) proteins, can target mRNAs that contain complementary seed sequences in their 3′ untranslated regions (UTRs) for translation inhibition and degradation (Bartel, 2018; Djuranovic et al., 2012; Eichhorn et al., 2014). In particular, the miR-15/16 family are potent regulators of cell cycle and survival (Liu et al., 2008). Previous work from our group utilized conditional deletion of the Mirc10 and Mirc30 loci that each contain two of the four major miR-15/16 family members (miR-15a, miR-15b, miR-16-1, and miR-16-2) driven by CD4-cre transgene (hereafter referred to as miR-15/16Δ/Δ mice or T cells) to demonstrate that this miRNA family has important effects on cytotoxic T cells (Gagnon and Ansel, 2019). In response to viral infection, miR-15/16Δ/Δ mice generate more viral-antigen-specific T cells, and these cells preferentially differentiate into memory cells that express CD127 and CD27 (Gagnon and Ansel, 2019).
In contrast to the well-defined roles of miRNAs, lncRNAs as a class do not have a single defined function. They are broadly defined as RNAs transcribed by polymerase II, over 200 nucleotides (nt) in length, that lack a translated open reading frame (Su et al., 2021; Wilusz et al., 2009). Some lncRNAs bind to chromatin and regulate nearby or distant genes, and others scaffold transcription factors and other protein complexes (Kopp and Mendell, 2018). Most relevant to the present study, some lncRNAs act as competing endogenous RNAs (ceRNAs) that bind to miRNAs, preventing their binding to and subsequent repression of mRNA targets (Poliseno et al., 2010; Su et al., 2021). However, rigorous connections between lncRNA physiological functions and molecular mechanisms of action have been hampered by the lack of precise tools to facilitate their study (Ponting and Haerty, 2022). As lncRNAs lack an open reading frame, insertion or deletion (indel) mutations do not reliably create null alleles. Instead, investigators have frequently deleted large genomic regions around the promoter or excised the genomic locus of the transcript altogether. These approaches risk disrupting cis-regulatory elements or topologically associated domains that control other nearby genes. Other studies that identified physiological roles for lncRNAs rely on RNAi-mediated lncRNA degradation. Both this approach and genetic manipulations that block lncRNA transcription leave a mechanistic gap between sequence features and downstream function. Studies demonstrating sequence-dependent function of lncRNAs are comparatively rare (Carrieri et al., 2012; Elguindy and Mendell, 2021; Faghihi et al., 2008; Gong and Maquat, 2011; Kleaveland et al., 2018; Lee et al., 1999; Yoon et al., 2012).
One lncRNA that has garnered much attention is the Metastasis Associated Lung Adenocarcinoma Transcript 1 (Malat1). Malat1 was first identified as highly expressed in both malignant tumors and healthy lung and pancreas with high interspecies conservation (Ji et al., 2003). Malat1 is mostly localized to the nucleus, where it is found within nuclear speckles, although it is not necessary for their formation (Nakagawa et al., 2012). These characteristics have generated numerous hypotheses about the function of Malat1 including that it scaffolds epigenetic and splicing factors and acts as a ceRNA to inhibit a variety of miRNAs (Su et al., 2021; Tripathi et al., 2010). In the immune system, Malat1 has been studied in dendritic cells, macrophages, T helper cells, T regulatory cells, and cytotoxic T cells (Hewitson et al., 2020; Kanbar et al., 2022; Masoumi et al., 2019; Wu et al., 2018). Malat1 is actively regulated during cytotoxic T cell differentiation with higher expression in memory precursors and lower expression in short-lived effector cells during LCMV infection (Kakaradov et al., 2017). Similarly, T cells stimulated in vitro express decreased levels of Malat1 over time (Masoumi et al., 2019). Early reports using large genetic deletions in the mouse Malat1 locus detected widespread gene expression changes but failed to identify a role for Malat1 in mouse development or cytotoxic T cell responses (Nakagawa et al., 2012; Yao et al., 2018; Zhang et al., 2012), though a recent report using RNAi detected altered T cell responses (Kanbar et al., 2022). Clearly, there remains much to learn about the sequence-specific mechanisms by which Malat1 tunes pathways and functions essential to cytotoxic T cells.
In the present study, we present a non-coding RNA circuit that regulates the activating signals from CD28 and IL-2, leading to low Bcl-2 expression and loss of memory cells following LCMV infection. We do so by identifying candidate interactions between miRNAs and lncRNAs via a targeted biochemical approach, and by creating CRISPR-targeted transgenic mice with precise mutation of Malat1 to interrogate the physiological function of the interaction between miR-15/16 and Malat1. We use this novel mouse to identify CD28-responsive gene programs affected by this circuit, and provide a new sequence-specific function of Malat1 in vivo.
Results
Malat1 is highly bound by miR-15/16
To identify candidate non-coding ceRNAs in cytotoxic T cells, we performed Argonaute-2 high-throughput sequencing of crosslinking immunoprecipitation (Ago2 HITS-CLIP, AHC). Integrating sequence reads across different classes of transcribed genomic annotations revealed that lncRNAs are bound extensively by Ago2, but the median lncRNA had 7.8 times fewer aligned AHC sequence reads compared to 3′ UTRs where miRNAs canonically bind to mRNAs (Figure 1A). This difference likely reflects the comparatively low expression of many lncRNAs, and it highlights the relatively low occupancy of Ago2 on most of these transcripts. Nevertheless, there were individual lncRNAs that stood out as highly bound across the transcript. To prioritize lncRNAs for further investigation, we manually curated the transcripts with the largest number of aligned AHC sequence reads. Among the top 10 most highly bound transcripts, 3 overlapped with protein-coding genes, 2 were on the mitochondrial chromosome, and 3 were repetitive annotations. Strikingly, all of the top 10 transcripts were either annotated as rRNA or contained an rRNA repeat element, except the eighth most bound transcript, Malat1 (Table 1). When the same analysis was repeated with rRNA repeats masked, Malat1 was the second most highly bound transcript (Table 2).
Malat1 has been proposed to inhibit miRNAs as a ceRNA (Chen et al., 2017; Luan et al., 2016; Qiao et al., 2018; Wang et al., 2019; Wu et al., 2018; Xiao et al., 2015; Xie et al., 2017). We used the Piranha peak calling algorithm to identify sites with the highest degree of miRNA binding, as indicated by AHC sequence read number and density (Figure 1B). Within cytotoxic T cells, the algorithm identified 55 AHC peaks in Malat1, the largest of which was extremely pronounced (Figure 1C). Compared with other lncRNA-binding peaks, this peak had the 36th most total aligned reads and the 142nd highest read density (15th and 16th, respectively, when rRNA reads are masked). Even when compared to peaks in 3′ UTRs of mRNAs, the largest peak in Malat1 was the 100th most (98th percentile) bound peak in terms of read density (Figure 1B).
miRTarget, the custom miRNA-binding prediction algorithm (Liu and Wang, 2019), identified an 8-mer seed-binding sequence of the miR-15/16 family centered within the most densely bound region of the called peak (Figure 1C). miR-15/16 has not been previously shown to interact with Malat1 in mice, so to determine whether this peak is miR-15/16 dependent, we performed AHC with cultured CD8+ miR-15/16Δ/Δ T cells. In the absence of miR-15/16, Ago2 binding to Malat1 was preserved throughout the whole transcript except for the peak containing the predicted miR-15/16 8-mer seed-binding sequence (Figure 1C). We conclude that miRNAs of the miR-15/16 family bind abundantly to Malat1. Furthermore, miR-15/16 binding occurred in a region of high evolutionary conservation (Figure 1D). The corresponding region containing the miR-15/16-binding site in human MALAT1 was highly enriched in two publicly available AHC datasets obtained using the 293 human embryonic kidney (HEK) cell line, consistent with Malat1 regulating miR-15/16 in a colorectal cancer cell line (Karginov and Hannon, 2013; Ji et al., 2019; Li et al., 2018; Figure 1E). We conclude that miR-15/16 bind abundantly to Malat1 in mouse and human cells.
We hypothesized that Malat1 may inhibit the function of miR-15/16 in cytotoxic T cells. To directly address this question in mice, we used CRISPR-Cas9 with homology directed repair to generate mice in which five nucleotides of the miR-15/16 seed-binding sequence within Malat1 were scrambled (Figure 1F). Mice homozygous for this mutation are subsequently referred to as Malat1scr/scr. To confirm the targeted functional outcome of this mutation, AHC was performed on cultured CD8+ T cells isolated from Malat1scr/scr mice. Ago2 binding was preserved across the Malat1 transcript except at the mutated miR-15/16-binding site, where AHC sequence read density was greatly reduced (Figure 1C). This finding confirmed that miR-15/16 binds to Malat1 at this site in a sequence-dependent manner, and provided us with a tool for highly specific investigation of the functional consequences of Malat1:miR-15/16 interaction.
Malat1 inhibits miR-15/16 availability and activity
Malat1:miR-15/16 interaction could lead to regulation and/or degradation of the lncRNA, the miRNAs, or both. To assess whether miR-15/16 degrades Malat1 we compared the expression of Malat1 by mRNA sequencing in primary mouse CD8+ T cells. Malat1 expression was unchanged in miR-15/16Δ/Δ cells as well as in the Malat1scr/scr cells (Figure 2A). Previous studies investigating Malat1 and other lncRNAs as ceRNA inhibitors of miRNAs have suggested that inhibition occurs by either target RNA-directed miRNA degradation (TDMD) or stoichiometric sequestration of the miRNA from protein-coding mRNA targets (Han et al., 2020; Su et al., 2021). Therefore, we also tested the possibility that Malat1 lowers miR-15/16 abundance by TDMD or a related mechanism. However, in freshly isolated mouse CD8+ T cells, miR-15b and miR-16 were unchanged and only miR-15a, a family member with lower expression, was modestly decreased in Malat1scr/scr cells (Figure 2B). We conclude that Malat1 and the miR-15/16 family do not influence each other’s absolute abundance in this setting. To test whether Malat1 affects miR-15/16 function, we examined miR-15/16 target binding in our AHC data. We first defined an experimentally supported list of TargetScan predicted miR-15/16-binding sites with at least one AHC read in both WT and Malat1scr/scr cells. Using this list, we then examined the read depth at these sites in both Malat1scr/scr and WT cells (603 sites contained in 479 genes). In WT cells, these sites constituted on average 3.2% of the binding in each 3′ UTR, and this figure increased to 3.4% in Malat1scr/scr cells, indicating that Ago2 occupancy preferentially increased at these sites when miR-15/16 binding to Malat1 was eliminated (Figure 2C). To confirm that binding to these target sites was miR-15/16 dependent, we assessed their AHC read depths in miR-15/16Δ/Δ cells as well. As predicted, binding was greatly reduced in miR-15/16Δ/Δ cells, representing, on average, 1.0% of binding to a given 3′ UTR, a 69% reduction in Ago2 binding compared to WT (Figure 2D). In contrast, no significant differences were observed for binding at predicted sites for the highly expressed miRNA families of miR-101, Let-7, miR-21, and miR-142, (Figure 2—figure supplement 1). Thus, we conclude that the Malat1scr allele specifically negatively regulates the first requirement of miR-15/16 function – binding to mRNA targets.
We next sought to investigate whether this increased binding resulted in decreased target mRNA expression. To do so, we analyzed mRNA sequencing data from primary CD8+ T cells 24 hr after stimulation. We generated empirical CDF plots from these data comparing target gene expression in miR-15/16fl/fl to miR-15/16Δ/Δ cells as well as WT to Malat1scr/scr cells. For each of these comparisons, we then compared the distribution for miR-15/16 target genes, as defined above, to the distribution for all other expressed genes. In the case of WT to Malat1scr/scr cells, the distribution was shifted in favor of decreased target gene expression in the Malat1scr/scr cells (Figure 2E). In the case of miR-15/16fl/fl to miR-15/16Δ/Δ, the distribution was shifted in favor of increased target gene expression in the miR-15/16Δ/Δ cells, indicating that this gene set is relieved of miRNA-induced repression (Figure 2F). These data indicate that the increased availability of miR-15/16 leads to increased repression of mRNA target genes when Malat1:miR-15/16 interaction is ablated.
We further investigated genes that displayed reciprocal expression changes in miR-15/16Δ/Δ and Malat1scr/scr cells. Of 479 genes in the bound target list, 432 (90%) were either downregulated in Malat1scr/scr cells or upregulated in miR-15/16Δ/Δ cells, compared to controls. Among these genes, the expression of 298 (62%) were decreased in Malat1scr/scr cells, and 294 (61%) were increased in miR-15/16Δ/Δ cells, with 160 (33%) both increased in miR-15/16Δ/Δ and decreased in Malat1scr/scr cells (Figure 2G). Gene ontology enrichment analysis revealed multiple pathways associated with growth factor and antigen receptor signaling affected by the Malat1:miR-15/16 circuit in cytotoxic T cells (Figure 2H). Many of these modules, including the T cell activation module, were identified because they share key signaling proteins. For example, gene expression and AHC data support direct miR-15/16 targeting of the costimulatory receptor Cd28, the alpha subunit of phosphatidylinositol 3-kinase (Pik3r1) and c-Jun N-terminal kinase (JNK1), encoded by Mapk8 (Figure 2—figure supplement 2 and Figure 3). Targeting of these genes by miR-15/16 is consistent with their known tumor suppressor role (Cimmino et al., 2005; Gagnon and Ansel, 2019) and with Malat1’s association with cancer cell proliferation and metastasis (Ji et al., 2003), providing new mechanistic insight into those observations.
Malat1 enhances CD28 expression and downstream CD8 T cell activation
Given that the Malat1:miR-15/16 circuit regulated genes essential to T cell activation, and in particular CD28, we next investigated the functional consequences of perturbing this circuit in cytotoxic T cells. First, we looked directly at CD28. AHC in WT cells detected a prominent peak in the CD28 3′ UTR at the TargetScan predicted binding site for miR-15/16. This peak was absent in miR-15/16Δ/Δ cells, whereas other binding peaks were preserved, empirically verifying that this binding event is miR-15/16 dependent. AHC in Malat1scr/scr cells indicated a modest increase in Ago2 binding at this site compared to WT cells relative to Ago2 binding in the whole 3′ UTR (Figure 3A). These data, along with our previous demonstration that the CD28 3′ UTR is miR-15/16 responsive, indicate that CD28 is part of a module of miR-15/16 target genes that are highly likely to be affected by the Malat1scr allele (Gagnon and Ansel, 2019).
Indeed, CD28 expression and T cell activation were regulated by the Malat1:miR-15/16 circuit. Flow cytometric measurement in primary splenic CD8+ T cells revealed that CD28 protein expression was decreased in Malat1scr/scr mice and enhanced in miR-15/16Δ/Δ mice (Figure 3B, C). To investigate activation-induced gene expression, we performed RNA-seq on splenic CD8+ T cells stimulated for 24 hr with plate-bound ɑCD3 with or without ɑCD28 crosslinking antibodies diagrammed in (Figure 3D). CD28 stimulation enhances distinct activation-induced gene expression changes in T cells (Martínez-Llordella et al., 2013), and these CD28-responsive genes were altered in miR-15/16Δ/Δ and Malat1scr/scr cells (Figure 3E–H). The previously defined set of 164 genes that are upregulated in WT cells stimulated with ɑCD3 + ɑCD28 compared to ɑCD3 alone was also significantly upregulated in the ɑCD3 + ɑCD28 condition compared to ɑCD3 alone for each genotype tested in our experiments (Figure 3—figure supplement 1). Importantly, the Malat1:miR-15/16 circuit affected this gene set in the ɑCD3 alone condition, with Malat1scr/scr cells exhibiting decreased expression compared to WT cells, and miR-15/16Δ/Δ cells exhibiting increased expression compared to miR-15/16fl/fl cells (Figure 3E, G). This trend was preserved, but to a lesser degree, in the ɑCD3 and ɑCD28 conditions (Figure 3F, H). Thus, Malat1:miR-15/16 interaction enhanced expression of a costimulation-responsive gene expression program in activated T cells, and it further enhanced expression of that module even when CD28 costimulation was directly engaged.
Unsupervised hierarchical clustering of these samples based on the expression of the 164 costimulation-responsive gene set further underscored the costimulatory-like effect of the Malat1:miR-15/16 circuit. In comparing gene expression in Malat1scr/scr and WT samples, three major groups emerged. The group with the lowest average costimulation-responsive gene expression contained only samples stimulated with ɑCD3 alone and primarily Malat1scr/scr samples. The group with intermediate expression was the largest group with an even representation of Malat1scr/scr samples and WT samples. While samples in this intermediate group were from both stimulation conditions, the Malat1sc/scr samples tended to be from the ɑCD3 + ɑCD28 condition and the WT samples had an even representation from both stimulation conditions. The group with the highest average expression contained predominantly ɑCD3 + ɑCD28 stimulated samples with an even representation of Malat1scr/scr and WT samples (Figure 3I). Thus, Malat1:miR-15/16 interaction and engagement of CD28 signaling additively induced costimulation-responsive genes. This observation is further supported when clustering the miR-15/16Δ/Δ and miR-15/16fl/fl samples. The two largest clusters were divided nearly exclusively by genotype. The group with lower average expression of costimulation-responsive genes contained predominantly mir-15/16fl/fl samples stimulated with ɑCD3 alone. The group with the higher average expression contained a majority of miR-15/16Δ/Δ samples, and the only miR-15/16fl/fl samples within the group received ɑCD3 + ɑCD28 stimulation. Within the higher expression group, the miR-15/16fl/fl samples stimulated with ɑCD3 + ɑCD28 subclustered with miR-15/16Δ/Δ samples stimulated with ɑCD3 alone (Figure 3J). A third, much smaller cluster was composed of both samples from a single outlier biological replicate. Overall, these data show that miR-15/16 restrict costimulation-responsive gene expression, and that Malat1:miR-15/16 interaction limits this effect.
In addition to the proximal changes in gene expression, downstream functional outcomes of CD28 costimulation were also affected. Early activation genes have been well described in T cells, with CD69 responding to many cues including TCR and CD28 ligation, and Nur77 responding very specifically to TCR signals (Ashouri and Weiss, 2017; Vandenberghe et al., 1993). In similar fashion as above, we assessed the expression of these proteins 2 and 4 hr after stimulation with plate-bound crosslinking antibodies on primary CD8+ T cells. CD69 exhibited decreased expression in Malat1scr/scr cells at both 2 and 4 hr post stimulation (2 hr: p = 0.0014 and 4 hr: p = 0.0598) and increased expression in miR-15/16Δ/Δ cells at 4 hr post stimulation only (p < 0.0001) (Figure 4A, B), in accordance with costimulation-responsive gene module expression (Figure 3). There was no difference, however, in expression of Nur77 across Malat1scr/scr, WT, and miR-15/16Δ/Δ cells, indicating that TCR signals were equivalent (Figure 4A, C). Another key consequence of CD28 ligation is the production of IL-2 (Fraser et al., 1991; Maggirwar et al., 1997; Vandenberghe et al., 1993). Therefore, we assessed cytokine production in the supernatants of these cultures 16 hr after stimulation. In line with previous findings in LCMV-infected mice (Gagnon and Ansel, 2019), miR-15/16Δ/Δ cells made significantly more IL-2 and TNFɑ after ɑCD3 stimulation both with and without ɑCD28, with a trend toward increased IFNγ as well (Figure 4D–F). Malat1scr/scr cells stimulated with ɑCD3 + ɑCD28 exhibited a reciprocal trend specifically for IL-2 (23% decrease compared to WT, p = 0.26) (Figure 4D). Together, these data show that the Malat1:miR-15/16 circuit regulates key functional outcomes of CD28-mediated costimulation, from proximal gene expression changes to early activation protein expression and cytokine secretion.
Malat1 enhances cytotoxic memory T cell differentiation
miR-15/16 restrict memory T cell differentiation, cell cycle, and cell survival during the response to LCMV Armstrong infection (Gagnon and Ansel, 2019). CD28 costimulation is essential for IL-2 production, memory cell formation, and memory recall responses in vivo (Borowski et al., 2007; Eberlein et al., 2012; Fuse et al., 2008; Grujic et al., 2010; Suresh et al., 2001). IL-2 is also essential for these same processes and in particular is required in a CD8+ T cell-intrinsic nature for the formation of CD127+ KLRG1− memory cells (Blattman et al., 2003; Kahan et al., 2022; Pipkin et al., 2010; Toumi et al., 2022; Whyte et al., 2022). Since Malat1 inhibits miR-15/16 activity and this circuit impacts proper T cell activation and IL-2 production after CD28 costimulation, we hypothesized that Malat1scr/scr cells would exhibit poor memory formation and survival.
We first examined steady-state memory populations of polyclonal T cells in unchallenged, young mice. Malat1scr/scr mice had normal naive and central memory T cell populations, but a reduced percentage and number of effector memory (CD44+CD62L–) cells in the spleen (Figure 5—figure supplement 1). Therefore, we investigated the cell-intrinsic nature of the Malat1scr/scr memory cell impairment during a viral challenge known to induce a large memory response. We bred the Malat1scr/scr mice with P14 TCR transgenic mice that express an antigen receptor specific for the immunodominant LCMV GP33 peptide. We then transferred Malat1scr/scr and WT P14 T cells into congenic CD45.1 hosts followed by LCMV Armstrong infection and tracked the acute and memory responses in the spleen and liver, the sites of primary LCMV infection (Figure 5A; Matloubian et al., 1993).
The total numbers of P14 cells were similar in recipients of transferred WT or Malat1scr/scr P14 cells in both the spleen and liver at day 7, and in the spleen at day 31. However, there was a reduction in Malat1scr/scr P14 cells in the liver at day 31 (Figure 5B, C). Despite having no effect on cell expansion at the peak of infection and mixed cell number results in the memory phase, the Malat1scr allele had a distinct effect on the cellular phenotype across organs and time points. Malat1scr/scr P14 cells preferentially displayed a phenotype associated with terminally differentiated effector memory cells (t-TEM), defined by KLRG1 expression and lack of CD127 expression as shown by Milner et al., 2020, and a corresponding reduction in the percentage of KLRG1− CD127+ memory cells (Figure 5D–F). Previous work has shown that substantial heterogeneity exists in the memory pool and that CD27 and CD43 can be useful in delineating functional differences such as recall potential between different memory cells. For instance, CD27+ memory cells tend to produce more IL-2 than CD27− cells (Milner et al., 2020), and CD43+ CD27+ cells are more effective at clearing Listeria upon re-challenge (Hikono et al., 2007; Olson et al., 2013). We therefore assessed these markers on the transferred Malat1scr/scr and WT P14 cells. Consistent with increased KLRG1+ CD127− t-TEM cells there were proportionally more CD27− CD43− cells in both the spleen and liver at day 31 among the Malat1scr/scr transferred cells compared to WT (Figure 5G, I). CD27+ CD43− cells were unchanged. The proportional reduction in CD127+ KLRG1− cells corresponded with a decreased proportion of CD43+ CD27+ cells. Notably, the total numbers of all CD27− subsets were unchanged. The proportional differences we observed were primarily driven by a reduction in the number of CD43+ CD27+ Malat1scr/scr P14 cells compared with wild-type P14 cells (45% reduction in liver (p = 0.0058), 30% reduction in spleen (p = 0.113)) (Figure 5G–I). This phenotype is exactly reciprocal to the increase in CD27+ memory cells previously documented in miR-15/16Δ/Δ mice (Gagnon and Ansel, 2019). To test whether the effect of the Malat1scr allele was epistatic to miR-15/16, we bred Malat1scr/scr to miR-15/16Δ/Δ mice to generate triple mutant Malat1scr/scr miR-15/16Δ/Δ mice. The poly-clonal LCMV response in these mice phenocopied that of miR-15/16Δ/Δ with a WT allele of Malat1, indicating that the observed effects of the Malat1scr allele are epistatic to miR-15/16 (Figure 5—figure supplement 2). We conclude that Malat1 inhibits miR-15/16 in responding CD8+ T cells during LCMV infection, leading to fewer CD43+ CD27+ memory cells. Compared to LCMV, Listeria monocytogenes (LM) infection induces lower expression of multiple costimulatory ligands so the antigen-specific response is more sensitive to the ablation of CD28 costimulation (Welten et al., 2015). Therefore, we hypothesized that cytotoxic T cell response may be more impaired by the Malat1scr allele during LM infection. We first sought to understand the effect of miR-15/16 on cytotoxic T cells during primary LM infection. To do so, we directly infected polyclonal mir-15/16fl/fl and miR-15/16Δ/Δ mice with LM expressing the LCMV GP33 peptide (LM-GP33) (Figure 6A). In contrast to the increased antigen-specific CD8+ T cell numbers in LCMV-infected miR-15/16Δ/Δ mice (Gagnon and Ansel, 2019), we observed similar cell numbers in LM-GP33-infected miR-15/16Δ/Δ and miR-15/16fl/fl mice (Figure 6B). However, the effect on memory cell differentiation was still present with decreased percentages of CD127− KLRG1+ t-TEMs and increased percentages of CD127+ KLRG1− memory cells (Figure 6C, D). Again, miR-15/16Δ/Δ antigen-specific cells had a significantly lower proportion of CD43− CD27− cells and a trend toward proportionally more CD43+ CD27+ cells (p = 0.127) (Figure 6E, F).
Given that antigen-specific miR-15/16Δ/Δ T cells exhibited enhanced memory cell differentiation during LM-GP33 infection, we further tested the role of Malat1:miR-15/16 interaction in this model using the P14 adoptive transfer system (Figure 6F). Transferred Malat1scr/scr and WT P14 cells accumulated in the spleen and the liver to similar numbers at both days 7 and 31 following LM-GP33 infection (Figure 6H). In contrast to LCMV, there was no change in the proportion or number of CD127− KLRG1+ Malat1scr/scr P14 cells at day 7 in the spleen or liver (Figure 6I). However, in LM-GP33 infection at day 7, CD43 and CD27 expression did not explicitly mirror the phenotype observed by CD127 and KLRG1 expression. In the spleen, there was a proportional and numerical increase in CD43− CD27− P14 cells with the Malat1scr/scr genotype (Figure 6J). There was a decreased proportion of CD43+ CD27+ cells (p = 0.0531), but this was despite a trend toward increased numbers of these cells (p = 0.1056) in the Malat1scr/scr P14 cells (Figure 6K). In the liver, there were no significant trends in any of these populations at day 7 (Figure 6J, K).
Although some of the phenotypes observed with LCMV were attenuated or absent during acute infection with LM-GP33 at day 7, Malat1:miR-15/16 interaction had more pronounced effects on memory cell populations at later times post infection (Figure 6L, M). On day 31, KLRG1+ cells were increased in proportion and numbers in the Malat1scr/scr P14 cells (Figure 6N). In the liver, a similar trend in KLRG1+ cell proportion existed (p = 0.233), but the numeric effect was entirely absent (Figure 6N). Surprisingly, this effect was restricted to the KLRG1+ populations, as CD127+ KLRG1− memory cells were unaffected in proportion and number in both the spleen and liver (Figure 6O). In both organs CD43− CD27− t-TEMs were proportionally increased and CD43+ CD27+ memory cells were proportionally decreased in Malat1scr/scr P14 cells (Figure 6P, Q). The numeric underpinnings of these proportional changes were different in each organ, with a 94% increase in CD43− CD27− cells in the spleen, and a 40% loss of the CD43+ CD27+ population in the liver (Figure 6P, Q). Taken together, these data demonstrate that miR-15/16 restrict memory cell differentiation, and reveal the ability of Malat1:miR-15/16 interaction to enhance memory cell differentiation across infection contexts.
Malat1 enhances cytotoxic T cell IL-2 production and survival
Finally, we investigated how the Malat1:miR-15/16 RNA circuit regulates memory cell differentiation and accumulation, giving consideration to the many direct targets of miR-15/16 and the emergent indirect effects on IL-2 and other costimulation-responsive gene expression in CD8 T cells. The pro-survival protein Bcl2 is the first characterized target of miR-15/16 (Cimmino et al., 2005). Higher abundance of Bcl2 in memory cells (as compared with short-lived effector cells) aids their survival and persistence, counteracting their increased expression of pro-apoptotic factors such as Bim (Kurtulus et al., 2011). As such, changes in the balance of pro- and anti-apoptotic factors can have selective effects on the accumulation of memory T cells. Bcl2 expression in Malat1scr/scr P14 cells was reduced at day 31 in the CD43− CD27− population in both LCMV and LM-GP33 infection, and also in CD43+ CD27+ cells in LCMV infection (Figure 7A–C). This reduction is consistent with increased direct miR-15/16 action on the Bcl2 3′ UTR in the Malat1scr/scr cells, but it may also be affected by IL-2 and other costimulation responsive factors. Also consistent with reduced Bcl2 expression, increased proportions of dead Malat1scr/scr cells were observed in both the KLRG1+ and KLRG1− cell fractions at day 7 during LM-GP33 infection (Figure 7D, E). Very few dead cells were detected among P14 cells of either genotype at day 31. IL-2 is critical for the induction of Bcl2 in memory T cells to promote their survival (Toumi et al., 2022). Furthermore, IL-2 production by cytotoxic T cells predisposes them to differentiate into a memory phenotype, likely by an autocrine/paracrine mechanism (Kahan et al., 2022; Kurtulus et al., 2011; Wojciechowski et al., 2007). Given the enhanced ability of miR-15/16Δ/Δ cells to produce IL-2 and the trend toward impaired IL-2 production in Malat1scr/scr cells in vitro (Figure 4D), we assessed IL-2 production by Malat1scr/scr T cells in LCMV and LM-GP33-infected mice. Corroborating our in vitro observations, a lower proportion of Malat1scr/scr P14 cells produced IL-2 when stimulated ex vivo on day 7 in both infection models (Figure 7F, G). This defect was also shared by KLRG1+ and KLRG1− populations and limited to this critical early time point. Higher proportions of memory cells were able to produce IL-2 ex vivo at day 31 after infection with either LCMV or LM-GP33 infection, and Malat1scr/scr P14 cells produced equivalent amounts of IL-2 (Figure 7H).
As CD28 costimulation is key to IL-2 production, the observed early defect is consistent with less robust activation while antigen is still present (Shapiro et al., 1997). Taken together, these results indicate that Malat1scr/scr cells in the context of a viral or bacterial infection receive relatively poor activating cues and subsequently produce less IL-2 early during infection, contributing to a less robust pro-survival and pro-memory state. These findings illustrate how Malat1 regulation of miR-15/16 and its large target gene network can act through multiple connected nodes to coordinate gene expression programs essential to cytotoxic T cell responses (Figure 8).
Discussion
LncRNAs are a large, diverse class of gene products that perform important physiological functions through a variety of molecular mechanisms. However, functional requirements are typically tested using RNAi to degrade the lncRNA or by disrupting their transcription entirely with no paradigm to dissect sequence-specific functions. ‘Sponging’ miRNAs (acting as a ceRNA) is among the most frequently proposed mechanisms of lncRNA function. Networks of non-coding ceRNAs, miRNAs, and target genes likely do shape gene expression programs in many biological contexts. Yet definitive evidence that a lncRNA:miRNA interaction has a physiological effect in a living organism has remained elusive. Guided by a biochemical approach, we investigated the requirements for Malat1 interaction with miR-15/16 by surgically altering just five nucleotides within the endogenous 8 kb Malat1 transcript in mice. Using T cells from these animals, we rigorously attribute changes in costimulation-responsive gene expression and in vivo defects in memory T cell formation to this Malat:miR-15/16 circuit.
Malat1 is a pleiotropic lncRNA implicated in a multitude of processes, including scaffolding splicing and epigenetic regulators, binding to chromatin, and interacting with several miRNAs in different cell types (Arun et al., 2020). Its extremely high expression and interspecies conservation nominated Malat1 as a strong candidate to have an impact on the activity of even highly abundant miRNAs like miR-15/16 that have large effects on gene expression and cell behavior. Less abundant miRNAs can also serve essential functions (Wigton et al., 2021), and they are likely to be more susceptible to inhibition by Malat1 and other ceRNAs. The same biochemically driven sequence-specific approach used here could be applied to probe the physiological effects of other lncRNA:miRNA interactions, and it could also be extended to probe the requirement for interaction with Ezh2 and other proteins.
Another important future direction is to better understand how exactly Malat1 inhibits miR-15/16, given that TDMD does not appear to result from this interaction. Most Malat1 resides in the nucleus, whereas miRNAs reside and function in the cytoplasm. Advances in miRNA fluorescent in situ hybridization (miR-FISH) or the implementation of proximity-CLIP (Benhalevy et al., 2018) may help clarify whether Malat1 nuclear sequestration plays a role in its ceRNA function. Detailed mechanistic understanding of this experimentally tractable circuit may illuminate how ceRNAs may be leveraged by the cell or by cellular engineers to manipulate the miRNA pool and target gene expression.
There is a growing interest in Malat1 function within the immune system, where it is highly expressed in many cell types and often regulated in interesting ways, including differential expression in short lived effector and memory precursor cells produced by asymmetric division of activated T cells during LCMV infection (Kakaradov et al., 2017). Malat1 overexpression in dendritic cells promotes IL-10 production by T regulatory cells (Wu et al., 2018), and Malat1−/− mice exhibited reduced IL-10 and Maf expression by T helper cells in vivo and in vitro (Hewitson et al., 2020). Two prior reports addressed Malat1 function in cytotoxic T cell responses. In one report, no significant differences were observed in Malat1−/− mice (Yao et al., 2018). It is difficult to interpret negative data generated from whole-body Malat1 deficiency given the multitude of cell types in which Malat1 may act. A second recent report, using RNAi to suppress Malat1 expression specifically in antigen-specific CD8+ T cells, detected enhanced t-Tem differentiation (Kanbar et al., 2022). This observation, which contrasts with our findings using Malat1scr/scr mice, was suggested to occur via an epigenetic mechanism involving Malat1 interaction with Ezh2. This mechanism and the ceRNA function that we describe here are not mutually exclusive, and it is possible to observe divergent phenotypes from the different interventions used to probe Malat1 function in T cells. Our Malat1scr/scr mouse combined with the cell transfer system clearly demonstrates a cytotoxic T cell-intrinsic requirement for Malat1 to inhibit miR-15/16 to enhance activation and protect memory cell persistence. This is consistent with our previous work on the miR-15/16 family, and fitting with the increased expression of Malat1 in memory precursors (Gagnon and Ansel, 2019; Kakaradov et al., 2017). The present study builds upon these findings and provides detailed insight into the gene regulatory programs which may drive these phenotypes. Cytotoxic T cells lacking CD28 generate similar numbers of viral specific cells at the height of infection, but fail to maintain equivalent numbers during the memory phase (Borowski et al., 2007). Furthermore, a major outcome of CD28 ligation is the production of IL-2, an essential T cell growth factor, and this IL-2 production is required in a cell-intrinsic manner for generation of CD127+ memory cells (Kahan et al., 2022). Here, we demonstrate that not only does the Malat1-miR-15/16 circuit regulate CD28 costimulation, IL-2 production, and Bcl2 expression, but also the concordant outcome of maintaining CD43+ CD27+ memory cells. miR-15/16 targets a large number of genes, so it is unlikely that the effect on Bcl2 or CD28 alone is entirely responsible for the observed changes in CD43− CD27− t-TEM and CD43+ CD27+ memory cells. Our results are consistent with memory cells requiring a heightened level of pro-survival cues to survive the contraction bottleneck after the peak of infection. In this manner elevated Malat1 levels can act as an enhancer of pro-survival and activating cues to prevent excessive cell death in this sensitive memory population. The magnitude of this effect was modest in acute LCMV and Listeria infection, two models that feature robust pathogen clearance, allowing assessment of memory T cells in the absence of chronic antigen persistence. Further work is needed to probe other settings in which Malat1:miR-15/16 interaction may have a bigger impact on the outcome of immune responses.
The importance of cytotoxic T cells is evident in the context of the global COVID-19 pandemic and the advent of CAR-T cell therapies reaching the clinic. Knowledge of the networks that regulate these cells and the critical nodes within these networks have the potential to augment technologies and therapies from vaccines to anti-tumor immunotherapy. Transcription factors and miRNAs have been extensively studied in this fashion as critical nodes in the regulation of gene transcription and translation. In the study presented here, we show Malat1 acts upstream of one such node, the miR-15/16 family. From this one interaction, Malat1 has the potential to combinatorially and synergistically regulate gene networks essential to cytotoxic T cells. This concept is easily extended when Malat1’s ability to regulate multiple other miRNA families is considered. If Malat1’s sequence-specific function is further defined, then editing or expression of specific Malat1 sequences could be used to tune multiple miRNA families in concert while leaving other functions of this enigmatic transcript un-touched.
Materials and methods
Mice
WT C57BL/6 mice were bred in our facility. miR-15/16Δ/Δ and miR-15/16fl/fl mice were derived as described in Gagnon and Ansel, 2019. Malat1scr/scr mice were generated from WT C57BL/6 mouse zygotes electroporated with CRISPR-Cas9 RNPs and HDR template (guide and template sequences below) as described previously (Chen et al., 2016). B6.SJL-Ptprca Pepcb/BoyJ (CD45.1) strain #002014 were purchased from the Jackson Laboratory. WT P14 mice carrying the TCR transgene specific to the LCMV GP33 peptide were obtained from the Waterfield lab and were backcrossed to C57BL/6 to maintain the line. Malat1scr/scr mice were bred to the P14 line to generate Malat1scr/scr P14 mice. Male and female age- and sex-matched mice were used between 5 and 12 weeks of age. All mice were housed and bred in specific pathogen-free conditions in the Animal Barrier Facility at the University of California, San Francisco. Animal experiments were approved by the Institutional Animal Care and Use Committee (IACUC) of the University of California, San Francisco (Protocol AN200003).
AGO2 HITS-CLIP
Request a detailed protocolCD8+ T cells were isolated as below and stimulated for 3 days with αCD3 and αCD28 antibodies and grown with kool aid complete media. Subsequently, the cells were rested and expanded in kool aid complete media for 2 days supplemented with 100 U/ml recombinant Human IL-2 (R&D Systems Cat# 202-IL-010/CF). Subsequently, 1 × 106 cells were used to prepare NGS libraries as described in Gagnon and Ansel, 2019 and Loeb et al., 2012. Samples were sequenced on HI-Seq 2500 (Illumina). Eleven nucleotide adaptors were trimmed from each read and resultant sequences were aligned to the mm10 genome using bowtie2 (Langmead and Salzberg, 2012). To assure lack of miR-15/16 binding in Malat1scr/scr mice was not do to errors in alignment, reads from Malat1scr/scr cells were aligned to a mm10 genome that contained a single modification changing Malat1 from the WT allele to the Malat1scr allele. To determine maximum binding depth across the genome and to manipulate aligned files the samtools package was used (Li et al., 2009). To assess Ago2 binding, aligned HITS-CLIP reads were integrated across the follow genomic annotations: lncRNA genes from mouse Noncode v6 (http://www.noncode.org/download.php); miRNA target-binding sites from TargetScan V7.2 (https://www.targetscan.org/cgi-bin/targetscan/data_download.mmu80.cgi); mm10 introns, exons, 3′ UTRs, and coding genes from UCSC genome table browser (http://genome.ucsc.edu/cgi-bin/hgGateway?db=mm10). To assess differences in target binding at specific 3′ UTR sites between WT and Malat1scr/scr cells, reads within each 3′ UTR were normalized to the total reads contained in the given 3′ UTR for each genotype. To identify regions of significant Ago2 binding above background, so-called ‘peaks’, we used the Piranha algorithm (Uren et al., 2012) (https://github.com/smithlabcode/piranha) on our HITS-CLIP. To remove sites from analysis where confounding features were present, such as miRNA expression sites or rRNA repeats, bedtools intersect was used to remove these features from annotation files (Quinlan, 2014). Human AHC data from publicly available datasets were trimmed and aligned to the Hg38 genome as previously described (Karginov and Hannon, 2013; Li et al., 2018).
Cell transfer and infections
Request a detailed protocolSpleens were harvested from WT or Malat1scr/scr P14 CD45.2 mice into phosphate-buffered saline (PBS) and passed through 70 μm strainers to generate a single-cell suspension. Samples were then stained for live/dead (apc-cy7), Thy1.2 (BV605), CD8ɑ (apc), and TCRVɑ2 (pe). Whole splenocytes were then transferred retro-orbitally (r.o.) into WT BoyJ (CD45.1) 7-week-old male recipients such that each recipient received 20,000 Thy1.2+ CD8ɑ+ TCRVɑ2+ cells in 200 μl PBS. To initiate LCMV infections, mice were injected intraperitoneally (i.p.) with 2 × 105 plaque forming units (p.f.u.) LCMV armstrong in 200 μl plain RPMI. LCMV virus was produced in house as described in Shehata et al., 2018. To initiate L. monocytogenes infection, mice were injected r.o. with 2 × 105 colony forming units (c.f.u.) Listeria monocytogenes-GP33 in 100 μl PBS. Listeria monocytogenes-GP33 was prepared in house as described in Allen et al., 2020. Blood was collected via submandibular bleeds with goldenrod 4 mm lancets collected into sample tubes coated with K2 EDTA (BD Ref# 365974) and 60 μl of blood was lysed with 500 μl of ACK lysis buffer. To assay spleens and livers, mice were sacrificed and organs harvested. To assay intracellular cytokine production, splenocytes were plated into 96-well u-bottom plates in complete kool-AID media and stimulated for 6 hr with a final concentration of 0.2 mg/ml GP33-41 (KAVYNFATM) in the presence of Brefeldin A.
Flow cytometry
Request a detailed protocolSpleens and livers were harvested into 2% fetal bovine serum (FBS) PBS with and passed through 70 μm strainers to generate a single-cell suspension. Samples were then spun at 450 r.c.f. for 5 min and livers were resuspended in 20% Percoll. Percoll suspensions were spun at 741 r.c.f. for 20 min at 25°C and the supernatants discarded. Cell pellets for both spleens and livers were resuspended in 1 ml ACK lysis buffer and incubated at 4°C for 5 min. Lysis was stopped with 5 ml 2% FBS PBS and samples again spun at 450 r.c.f. for 5 min. Samples were resuspended in 2% FBS and aliquoted into v-bottom 96 well plates. Live dead staining was then performed using the fixable viability dye at 1:2000 in PBS. Subsequently, cells were blocked with mouse Fc block 1:100 in 2% FBS PBS and stained for surface proteins with directly conjugated antibodies diluted 1:100 in 2% FBS PBS. Stains were incubated for 20 min at 4°C protected from light. For surface stains alone, LCMV samples were fixed with 4% paraformaldehyde (PFA) for 5 min at 4°C. For intracellular stains, samples were fixed and permeabilized according the the Transcription Factor FIxation Kit (Invitrogen Cat#00-5521-00). Intracellular antibodies were diluted in 1:100 in permeabilization buffer and incubated at 4°C for 30 min. All samples were spun at 821 r.c.f. for 5 min prior to being resuspended in 2% FBS PBS with 1:10 AccuCount beads (spherotech Cat #ACBP-100–10) for analysis on either the BD LSR II or the BD LSRFortessa flow cytometry analyzer.
Cell isolation and in vitro functional assays
Request a detailed protocolCD8+ T cells were isolated from spleens using negative selection from the EasySep Mouse CD8+ T Cell Isolation Kit (cat# 19853). Cells were counted using trypan blue staining and the nexelcom cellometer spectrum. Cultures were started via stimulation with plate-bound αCD3 (1 μg/ml, clone 2C11) and αCD28 (1 μg/ml, clone 37.51); plates coated overnight in PBS with Ca2+ and Mg2+ at 4°C. Cells were plated in Kool AID complete media (Dulbecco's Modified Eagle Medium (DMEM) high glucose media supplemented with 10% FBS, pyruvate, nonessential amino acids, minimum essential medium (MEM) vitamins, l-arginine, l-asparagine, l-glutamine, folic acid, beta mercaptoethanol, penicillin, and streptomycin) and spun at 450 r.c.f. for 5 min at 25°C to begin stimulation. For functional assays cells were harvested 2 or 4 hr after stimulation for flow cytometry, or 24 hr after stimulation for mRNA-seq. For supernatant cytokine expression, plates were spun at 450 r.c.f. for 5 min at 4°C 16 hr after stimulation and cell-free culture medium was collected and frozen at −80°C. TNFα (Cat #), IFNγ (Cat #), and IL-2 (Cat #) were analyzed by ELISA.
miRNA qPCR
Request a detailed protocolSpleens were harvested and single-cell suspensions generated by passing through 70 μm strainers. CD8+ T cells were isolated as above, and 2 × 106 cells were pelleted at 450 r.c.f. for 5 min at 4°C. Cell pellet was resuspended in 700 μl Trizol reagent (Ambion cat #15596018) and kept at –80 °C. RNA was was isolated using using the Direct-zol-96 RNA Kit (zymogen cat #R2054). This RNA was used as input to the Mir-X miRNA qRT-PCR TB Green Kit (Takara cat# 638316) to generate miRNA cDNA. Specific primers to miR-15a, miR-15b, and miR-16 were used to quantify those miRNA species on the QPCR MACHINE (eppendorf realplex2). ribosomalRNA 5.8 s was used as a housekeeping control for each sample. Each sample was run in duplicate. To quantify miRNA expression technical duplicates were averaged and then normalized to rRNA 5.8 s by subtracting the rRNA 5.8 s Ct value from miRNA Ct value (ΔCt). Expression values reported are generated by 2−ΔCt.
mRNA sequencing
Request a detailed protocol1 × 106 CD8+ T cells were harvested 24 hr after αCD3 ± αCD28 stimulation as described above. Cells were pelleted at 450 r.c.f. for 5 min at 4°C. Cell pellet was resuspended in 700 μl Trizol reagent (ambion ref #15596018) and kept at −80°C. RNA was was isolated using using the Direct-zol-96 RNA Kit (zymogen cat #R2054). The integrity of total RNA was checked on Fragment Analyzer (Agilent, Cat. No. DNF-472), only RNA with RQN number of above 7 was used for library construction. The starting quantity of 100 ng of total RNA was used according to vendor instructions with Universal plus mRNA with Nu Quant (TECAN, Cat. No. 0520), final library PCR amplification was 17 cycles. After library completion, individual libraries were pooled equally by volume, and quantified on Fragment Analyzer (Agilent, Cat. No. DNF-474). Quantified library pool was diluted to 1 nM and further diluted as per protocol and sequenced on Illumina MiniSeq (Illumina, Cat. No. FC-420-1001) to check for quality of reads. Finally, individual libraries were normalized according to MiniSeq output reads, specifically by % protein-coding genes and were sequenced on one lane of NovaSeq6000 S4 PE100 (Illumina, Cat. No. 20028313). Reads were aligned to the mouse genome GRCm38 and quantified using the STAR aligner software version 2.7.2b. Read normalization and differential expression analysis were performed in the R computing environment version 3.6.1 using the software DESeq2 version 1.26. For RNA sequencing analysis, FDR-corrected p-values were used to evaluate significant differences between experimental groups using a significance threshold of 0.05. Lowly expressed genes that had fewer than 2 reads per million on average across all samples were removed from the analysis. Empirical cumulative density plots were made using the ggplot package for R. Heatmaps and hierarchical clustering were performed via the gplots package for R. Gene ontology analysis was performed using the statistical overrepresentation test for Panther pathways (Mi et al., 2019).
Materials availability
Request a detailed protocolMalat1scr/scr mice will be provided upon request.
Statistical and analytical software
Request a detailed protocolAll flow cytometry data were analyzed using FlowJo (version 10). Statistical analyses and plotting were performed using GraphPad Prism (Version 9.2.0) and R (version 4.2.1).
Data availability and software
Request a detailed protocolAgo2 HITS-CLIP data uploaded to NCBI GEO accession #GSE216565.
mRNA sequencing data uploaded to NCBI GEO accession #GSE216113.
Code used to analyze HITS-CLIP data summed over annotations and peaks in WT cells can be found at https://github.com/AnselLab/WT_HITS_CLIP_Analyses, (copy archived at Wheeler, 2023a).
Code used to analyze relative binding densities of HITS-CLIP reads can be found at https://github.com/AnselLab/Malat1_miR1516_AGO2_HITS_CLIP, (copy archived at Wheeler, 2023b).
Code used to analyze mRNAseq data for both miR-15/16 target expression and CD28 responsive genes can be found at https://github.com/AnselLab/Malat1_miR1516_CD3CD28_RNAseq, (copy archived at Wheeler, 2023c).
Appendix 1
Data availability
Ago2 HITS-CLIP data deposited in GEO under accession code #GSE216565. mRNA sequencing data deposited in GEO under accession code #GSE216113.
-
NCBI Gene Expression OmnibusID GSE216565. The lncRNA Malat1 Inhibits miR-15/16 to Enhance Cytotoxic T Cell Activation and Memory Cell Formation [Ago2 HITS-CLIP].
-
NCBI Gene Expression OmnibusID GSE216113. The lncRNA Malat1 Inhibits miR-15/16 to Enhance Cytotoxic T Cell Activation and Memory Cell Formation.
-
NCBI Gene Expression OmnibusID GSE44404. Remodeling of Ago2-mRNA interactions upon cellular stress reflects miRNA complementarity and correlates with altered translation rates.
-
NCBI Gene Expression OmnibusID GSE102321. Transcriptome-wide identification and validation of interactions between the miRNA machinery and HuR on mRNA targets.
References
-
Programmed contraction of CD8(+) T cells after infectionNature Immunology 3:619–626.https://doi.org/10.1038/ni804
-
Therapeutic use of IL-2 to enhance antiviral T-cell responses in vivoNature Medicine 9:540–547.https://doi.org/10.1038/nm866
-
Memory CD8+ T cells require CD28 costimulationThe Journal of Immunology 179:6494–6503.https://doi.org/10.4049/jimmunol.179.10.6494
-
Highly efficient mouse genome editing by CRISPR ribonucleoprotein electroporation of zygotesJournal of Biological Chemistry 291:14457–14467.https://doi.org/10.1074/jbc.M116.733154
-
Long non-coding RNA MALAT1 regulates ZEB1 expression by sponging miR-143-3p and promotes hepatocellular carcinoma progressionJournal of Cellular Biochemistry 118:4836–4843.https://doi.org/10.1002/jcb.26158
-
MicroRNA regulation of CD8+ T cell responsesNon-Coding RNA Investigation 3:24.https://doi.org/10.21037/ncri.2019.07.02
-
Malat1 suppresses immunity to infection through promoting expression of maf and IL-10 in Th cellsJournal of Immunology 204:2949–2960.https://doi.org/10.4049/jimmunol.1900940
-
Activation phenotype, rather than central- or effector-memory phenotype, predicts the recall efficacy of memory CD8+ T cellsThe Journal of Experimental Medicine 204:1625–1636.https://doi.org/10.1084/jem.20070322
-
The long noncoding RNA Malat1 regulates CD8+ T cell differentiation by mediating epigenetic repressionThe Journal of Experimental Medicine 219:e20211756.https://doi.org/10.1084/jem.20211756
-
Bcl-2 allows effector and memory CD8+ T cells to tolerate higher expression of BimJournal of Immunology 186:5729–5737.https://doi.org/10.4049/jimmunol.1100102
-
Fast gapped-read alignment with Bowtie 2Nature Methods 9:357–359.https://doi.org/10.1038/nmeth.1923
-
Tsix, a gene antisense to Xist at the X-inactivation centreNature Genetics 21:400–404.https://doi.org/10.1038/7734
-
The sequence alignment/Map format and SAMtoolsBioinformatics 25:2078–2079.https://doi.org/10.1093/bioinformatics/btp352
-
Transcriptome-wide Identification and Validation of Interactions between the miRNA Machinery and HuR on mRNA TargetsJournal of Molecular Biology 430:285–296.https://doi.org/10.1016/j.jmb.2017.12.006
-
miR-16 family induces cell cycle arrest by regulating multiple cell cycle genesNucleic Acids Research 36:5391–5404.https://doi.org/10.1093/nar/gkn522
-
Regulation of the interleukin-2 CD28-responsive element by NF-ATp and various NF-kappaB/Rel transcription factorsMolecular and Cellular Biology 17:2605–2614.https://doi.org/10.1128/MCB.17.5.2605
-
CD28-inducible transcription factor DEC1 is required for efficient autoreactive CD4+ T cell responseThe Journal of Experimental Medicine 210:1603–1619.https://doi.org/10.1084/jem.20122387
-
Genome-wide analysis of human long noncoding RNAs: a provocative reviewAnnual Review of Genomics and Human Genetics 23:153–172.https://doi.org/10.1146/annurev-genom-112921-123710
-
BEDTools: the swiss-army tool for genome feature analysisCurrent Protocols in Bioinformatics 47:11.https://doi.org/10.1002/0471250953.bi1112s47
-
The role of a ceRNA regulatory network based on lncRNA MALAT1 site in cancer progressionBiomedicine & Pharmacotherapy = Biomedecine & Pharmacotherapie 137:111389.https://doi.org/10.1016/j.biopha.2021.111389
-
Role of CD28-B7 interactions in generation and maintenance of CD8 T cell memoryJournal of Immunology 167:5565–5573.https://doi.org/10.4049/jimmunol.167.10.5565
-
Ligation of the CD5 or CD28 molecules on resting human T cells induces expression of the early activation antigen CD69 by a calcium- and tyrosine kinase-dependent mechanismImmunology 78:210–217.
-
The identification and characterization of novel transcripts from RNA-seq dataBriefings in Bioinformatics 17:678–685.https://doi.org/10.1093/bib/bbv067
-
SoftwareWt_Hits_Clip_Analyses, version swh:1:rev:3ca9c23bf63268d7635cce0b6bb28d4f422f7de1Software Heritage.
-
SoftwareMalat1_Mir1516_Ago2_Hits_Clip, version swh:1:rev:6c64ec2466d2551745f798a321ed701069548468Software Heritage.
-
SoftwareMalat1_Mir1516_Cd3Cd28_Rnaseq, version swh:1:rev:cb1770dc40c609737a07aec5d78522836b60ebc5Software Heritage.
-
Memory CD8 T-cell differentiation during viral infectionJournal of Virology 78:5535–5545.https://doi.org/10.1128/JVI.78.11.5535-5545.2004
-
Context-dependent effects of IL-2 rewire immunity into distinct cellular circuitsThe Journal of Experimental Medicine 219:e20212391.https://doi.org/10.1084/jem.20212391
-
MicroRNA-directed pathway discovery elucidates an miR-221/222-mediated regulatory circuit in class switch recombinationThe Journal of Experimental Medicine 218:e20201422.https://doi.org/10.1084/jem.20201422
-
Long noncoding RNAs: functional surprises from the RNA worldGenes & Development 23:1494–1504.https://doi.org/10.1101/gad.1800909
-
Bim/Bcl-2 balance is critical for maintaining naive and memory T cell homeostasisThe Journal of Experimental Medicine 204:1665–1675.https://doi.org/10.1084/jem.20070618
-
LincRNA-p21 suppresses target mRNA translationMolecular Cell 47:648–655.https://doi.org/10.1016/j.molcel.2012.06.027
Article and author information
Author details
Funding
National Heart, Lung, and Blood Institute (R01 HL109102)
- K Mark Ansel
National Heart, Lung, and Blood Institute (P01 HL107202)
- K Mark Ansel
National Institute of Diabetes and Digestive and Kidney Diseases (P30 DK063720)
- Matthew H Spitzer
- Alexander Marson
- K Mark Ansel
Sandler Foundation (Sandler Asthma Basic Research Center investigator support)
- K Mark Ansel
Howard Hughes Medical Institute (Gilliam Fellowship)
- Priscila Muñoz-Sandoval
- K Mark Ansel
Hooper Foundation (Graduate student fellowships)
- Benjamin D Wheeler
- Wandi S Zhu
The funders had no role in study design, data collection, and interpretation, or the decision to submit the work for publication.
Acknowledgements
We thank Shomyseh Sanjabi, Nadia Roan, and Jason Neidleman for generous donation of the LCMV armstrong virus; Michael Waterfield for generous donation of P14 mice; Marlys Fasset, Simon Zhou, and Eric Wigton for maintaining mutant mouse colonies through the COVID-19 pandemic; Walter Eckelbar, Eunice Wan, Charina Julian, Lenka Maliskova, and Andrew Schroeder for assistance with mRNA-seq and the sequencing of our HITS-CLIP libraries. This work was supported by the NIH (HL109102 and HL107202), the Sandler Asthma Basic Research Center, and the Hooper Foundation. Priscila Muñoz-Sandoval is a Howard Hughes Medical Institute Gilliam Fellow. We acknowledge the PFCC (RRID:SCR_018206) for assistance generating flow cytometry data. Research reported here was supported in part by the DRC Center Grant NIH P30 DK063720.
Ethics
All mice were housed and bred in specific pathogen-free conditions in the Animal Barrier Facility at the University of California, San Francisco. Animal experiments were approved by the Institutional Animal Care and Use Committee (IACUC) of the University of California, San Francisco, protocol number AN2000003.
Version history
- Sent for peer review:
- Preprint posted:
- Reviewed Preprint version 1:
- Reviewed Preprint version 2:
- Version of Record published:
Cite all versions
You can cite all versions using the DOI https://doi.org/10.7554/eLife.87900. This DOI represents all versions, and will always resolve to the latest one.
Copyright
© 2023, Wheeler et al.
This article is distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use and redistribution provided that the original author and source are credited.
Metrics
-
- 1,018
- views
-
- 96
- downloads
-
- 6
- citations
Views, downloads and citations are aggregated across all versions of this paper published by eLife.
Download links
Downloads (link to download the article as PDF)
Open citations (links to open the citations from this article in various online reference manager services)
Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)
Further reading
-
- Evolutionary Biology
- Immunology and Inflammation
The incessant arms race between viruses and hosts has led to numerous evolutionary innovations that shape life’s evolution. During this process, the interactions between viral receptors and viruses have garnered significant interest since viral receptors are cell surface proteins exploited by viruses to initiate infection. Our study sheds light on the arms race between the MDA5 receptor and 5’ppp-RNA virus in a lower vertebrate fish, Miichthys miiuy. Firstly, the frequent and independent loss events of RIG-I in vertebrates prompted us to search for alternative immune substitutes, with homology-dependent genetic compensation response (HDGCR) being the main pathway. Our further analysis suggested that MDA5 of M. miiuy and Gallus gallus, the homolog of RIG-I, can replace RIG-I in recognizing 5’ppp-RNA virus, which may lead to redundancy of RIG-I and loss from the species genome during evolution. Secondly, as an adversarial strategy, 5’ppp-RNA SCRV can utilize the m6A methylation mechanism to degrade MDA5 and weaken its antiviral immune ability, thus promoting its own replication and immune evasion. In summary, our study provides a snapshot into the interaction and coevolution between vertebrate and virus, offering valuable perspectives on the ecological and evolutionary factors that contribute to the diversity of the immune system.
-
- Chromosomes and Gene Expression
- Immunology and Inflammation
Genome-wide association studies (GWAS) have identified hundreds of genetic signals associated with autoimmune disease. The majority of these signals are located in non-coding regions and likely impact cis-regulatory elements (cRE). Because cRE function is dynamic across cell types and states, profiling the epigenetic status of cRE across physiological processes is necessary to characterize the molecular mechanisms by which autoimmune variants contribute to disease risk. We localized risk variants from 15 autoimmune GWAS to cRE active during TCR-CD28 co-stimulation of naïve human CD4+ T cells. To characterize how dynamic changes in gene expression correlate with cRE activity, we measured transcript levels, chromatin accessibility, and promoter–cRE contacts across three phases of naive CD4+ T cell activation using RNA-seq, ATAC-seq, and HiC. We identified ~1200 protein-coding genes physically connected to accessible disease-associated variants at 423 GWAS signals, at least one-third of which are dynamically regulated by activation. From these maps, we functionally validated a novel stretch of evolutionarily conserved intergenic enhancers whose activity is required for activation-induced IL2 gene expression in human and mouse, and is influenced by autoimmune-associated genetic variation. The set of genes implicated by this approach are enriched for genes controlling CD4+ T cell function and genes involved in human inborn errors of immunity, and we pharmacologically validated eight implicated genes as novel regulators of T cell activation. These studies directly show how autoimmune variants and the genes they regulate influence processes involved in CD4+ T cell proliferation and activation.