miR-155 exerts posttranscriptional control of autoimmune regulator (Aire) and tissue-restricted antigen genes in medullary thymic epithelial cells

Background The autoimmune regulator (Aire) gene is critical for the appropriate establishment of central immune tolerance. As one of the main controllers of promiscuous gene expression in the thymus, Aire promotes the expression of thousands of downstream tissue-restricted antigen (TRA) genes, cell adhesion genes and transcription factor genes in medullary thymic epithelial cells (mTECs). Despite the increasing knowledge about the role of Aire as an upstream transcriptional controller, little is known about the mechanisms by which this gene could be regulated. Results Here, we assessed the posttranscriptional control of Aire by miRNAs. The in silico miRNA-mRNA interaction analysis predicted thermodynamically stable hybridization between the 3’UTR of Aire mRNA and miR-155, which was confirmed to occur within the cellular milieu through a luciferase reporter assay. This finding enabled us to hypothesize that miR-155 might play a role as an intracellular posttranscriptional regulator of Aire mRNA. To test this hypothesis, we transfected a murine mTEC cell line with a miR-155 mimic in vitro, which reduced the mRNA and protein levels of Aire. Moreover, large-scale transcriptome analysis showed the modulation of 311 downstream mRNAs, which included 58 TRA mRNAs. Moreover, miR-155 mimic-transfected cells exhibited a decrease in their chemotaxis property compared with control thymocytes. Conclusion Overall, the results indicate that miR-155 may posttranscriptionally control Aire mRNA, reducing the respective Aire protein levels; consequently, the levels of mRNAs encode tissue-restricted antigens were affected. In addition, miR-155 regulated a crucial process by which mTECs allow thymocytes’ migration through chemotaxis. Supplementary Information The online version contains supplementary material available at 10.1186/s12864-022-08631-4.

the expression levels in each individual mTEC, and this process ensures self-representation and self-reactive T cell elimination through negative selection [4,5]. These processes guarantee that central tolerance is extended to virtually all tissues and organs in the body.
One of the main controllers of PGE in mTECs is the autoimmune regulator (Aire) gene, which encodes the AIRE protein that acts as a noncanonical transcription factor that promotes the expression of tissue-restricted antigens (TRAs) [6][7][8][9]. According to Takaba et al. [10] 28.9% of the TRAs expressed in mTECs are under the control of Aire and are referred to as Aire-dependent genes. In addition, another 12.3% of TRAs are under the control of Aire and Forebrain Embryonic Zinc Finger-Like Protein 2 (Fezf2), and the other 58.8% are not under the control of Aire and are referred to as Aireindependent genes.
The Aire gene is expressed in only a few tissues, but it is mostly expressed in the thymus and, more specifically, in mTECs [11]. Additionally, only a subpopulation of mature mTECs with the CD80 hi MHC-II hi phenotype express the Aire gene; these cells are referred to as mTEC hi Aire + cells [1,5]. The AIRE protein acts as a transcriptional modulator of the chromatin of mTECs, promoting the expression of several Aire-dependent TRAs [2,3,11]. The importance of Aire in controlling the maturation and proliferation of mTECs, as well as its participation in the process of adhesion between mTECs and thymocytes during negative selection, have been shown [2,3,[11][12][13].
Immature CD80 low MHC-II low mTECs, which do not express Aire, express a smaller set of TRAs, while the mTEC hi Aire + cell population can express up to 62% of the protein-coding genes (mainly encoding TRAs) in mice [2,3,5].
In addition to TRAs, which are extremely important for testing the self-reactivity of thymocytes, studies have shown that several Aire-dependent genes have different functions in mTECs; these genes participate in processes such as the differentiation and maturation of mTECs [14], interactions of mTECs with thymocytes and dendritic cells [15,16], induction of thymocyte migration within the thymus, apoptosis of mTECs and even the processing of mRNAs (alternative splicing), thus contributing to the expression of TRA isoforms in mTECs [1-3, 5, 6, 17].
Our group, and independently, the former group of Kyewski/Ucar, has contributed to showing that Aire is involved in the expression of miRNAs in mTECs [12,18]. Recently, it was shown that Aire transactivates the expression of HLA-G [19] and that it can regulate the length of the 3′ untranslated regions (3'UTRs) of downstream mRNAs [20,21].
The AIRE protein is composed of five main subdomains: the CARD (caspase recruiter domain) domain, which acts in the homo-multimerization of the protein, since AIRE acts as a tetramer in the nucleus of mTECs; the NLS (nuclear localization signal) domain, which functions in the translocation of the protein to the nucleus; the SAND domain, which is involved in the protein-protein interactions that enable AIRE to act in conjunction with its partner proteins; and 2 PHD domains, namely, PHD-1 and PHD-2, separated by a proline-rich region (PRR), which act as a histone code reader and assists in chromatin decondensation [2,6,22]. AIRE acts as a transcriptional modulator because it does not recognize any particular motif in DNA to promote transcription. The AIRE protein interacts with the RNA Pol II molecule that has already started the transcription process. RNA Pol II is "anchored" to chromatin at specific T-box (TTA TTA ) and G-box (ATT GGT TA) sites in the promoter regions of Aire-dependent genes and is then released by AIRE, promoting the extension of the transcription of downstream genes [23].
AIRE cooperates with other proteins (called Aire partners) to promote the expression of Aire-dependent genes in mTECs. Sirtuin 1 (SIRT-1) promotes the deacetylation of the lysine residues of AIRE to promote its activation [24]. Once activated, AIRE recognizes and binds to epigenetic markers involved in transcriptional repression, such as H3K4me0 [25]. Topoisomerase II (TOP2a) associates with AIRE and promotes double-strand DNA breaks. DNA PKs and other partners are activated by double-strand DNA breaks and promote the relaxation of chromatin [17]. The protein complex of bromodomain-containing 4 (Brd4) [26] and positive transcription elongation factor b (pTEFb) [27] are associated with AIRE and play a role in releasing RNA Pol II from being anchored to chromatin. Then, CREB-binding protein (CBP) acetylates the lysine residues of AIRE, leading to their inactivation [2,6].
Aire, as stated above, is an important upstream transcriptional controller in mTECs, but little is known about the elements or mechanisms that control this gene. The structural characteristics of AIRE and its association with its partners control the molecular action of the AIRE protein in a posttranslational manner. Here, we ask which element(s) regulate the Aire gene posttranscriptionally, that is, after the transcription of Aire mRNA and when it is located in the cytoplasm. Considering the predicted significant thermodynamic stability of the interaction of the Aire mRNA 3'UTR and miR-155, we hypothesized that Aire may be posttranscriptionally controlled in mTECs.

Prediction of the miR-155-Aire 3'UTR interaction and validation by luciferase reporter gene assay (LRGA)
The mmu-miR-155-5p (miR-155) has the following sequence: UUA AUG CUA AUU GUG AUA GGGGU. Analyses performed with the RNAhybrid algorithm predicted that the Aire mRNA 3'UTR strongly hybridizes to this miRNA, exhibiting a minimal free energy (mfe) equal to − 22 kcal/mol, which gives the structure thermodynamic stability (Fig. 1).
The miRNA-mRNA interaction was then validated using the LRGA. The results showed that the luciferase expression decreased in cells cotransfected with the pMIR-Glo/Aire 3'UTR vector and the miR-155 mimic. This result shows that there was an interaction between the miR-155 mimic and the 3'UTR of Aire mRNA but not with a mutant 3'UTR or mimic negative control (Fig. 2). The interaction observed between miR-155, and Aire mRNA does not fall into the classic canonic seedmatched site as it occurs from the fourth miRNA nucleotide hybridizing with the mRNA 3'UTR. Our results show that mutations along the 3'UTR of Aire mRNA weaken miR-155 interaction.

Quantification of the intracellular levels of miR-155 and Aire mRNA
To assess the effect of miR-155 on Aire mRNA, mTEC 3.10 cells were transfected with the miR-155-5p mimic. RT-qPCR analysis of miR-155 was then performed to confirm the transfection efficiency and assess the intracellular levels of miR-155 at 12, 24, and 48 h posttransfection.
Compared with the control, transfection of the miR-155-5p mimic increased the intracellular level of miR-155, and peak expression was achieved at 12 h posttransfection. Twenty-four hours posttransfection, the levels of miR-155 were notably decreased and were statistically comparable to the levels in the control cells. The results showed that 12 h posttransfection was the time point that was most appropriate for evaluating the effect of increasing the levels of miR-155 in mTEC 3.10 cells (Fig. 3).
Samples of total RNA (cDNA) from mTEC 3.10 cells transfected with the miR-155 mimic were used for RT-qPCR analysis of Aire mRNA. We again focused on the  Posttranscriptional interaction between miR-155 with Aire 3'UTR was assayed by luciferase reporter gene assay (LRGA). pMIR-Aire-wt-3'utr (wild type 3'UTR) or pMIR-Aire-mut-3'utr (mutant 3'UTR) luciferase vector constructs were transfected into human HEK-293 T cells to demonstrate the possibility of occurrence of the miRNA-mRNA interactions within the cell milieu. The mutant 3'UTR was used to demonstrate the need of sequence specificity for the interaction to occur. Data presented are from three independent experiments, and standard error mean (SEM) time points of 12, 24, and 48 h posttransfection. The results showed a significant reduction in Aire mRNA expression 12 h after the peak level of miR-155 expression was observed, that is, at 24 h posttransfection. Forty-eight hours after miR-155 transfection, the Aire mRNA levels were comparable to those in the control cells (Fig. 4).

AIRE protein expression analysis Western blotting
The results obtained by RT-qPCR showed that the mRNA levels of Aire decreased 24 h after transfection with the miR-155-5p mimic. Thus, we decided to evaluate the protein levels of AIRE by western blot at this time point. AIRE protein expression in mTEC 3.10 cells was reduced (mean reduction of 37.3%) 24 h after transfection with the miR-155-5p mimic ( Fig. 5 and Supplemental Fig. 1).

Immunofluorescence
To evaluate the effect of miR-155-5p mimic transfection on AIRE protein levels in subcellular compartments, we used immunofluorescence. mTEC 3.10 cells were observed at 24 h posttransfection. Results show that the AIRE protein was located in the nucleus and visualized as red dots. Twenty-four hours posttransfection, there was a reduction of approximately 15% in the fluorescence signal intensity. The number of Aire + mTEC 3.10 cells remained comparable to that of the control group i.e., 4-5 Aire + cells in both groups. The results confirmed that transfection with the miR-155 mimic interfered with the AIRE protein expression level but not with the number of Aire + mTECs (Fig. 6).

Differentially expressed mRNAs in mTEC 3.10 cells transfected with the miR-155 mimic
The mTECs transfected with the miR-155-5p mimic exhibited the most significant effects on the mRNA and protein levels of Aire 24 h posttransfection. Therefore, these samples were selected for large-scale gene expression analysis (mRNA transcriptome) using the microarray technique and the limma linear model package for data analysis. mRNAs were considered to be differentially expressed when they exhibited a fold change of ≥1.5. As a cutoff, an adjusted p-value ≤0.05 (Benjamini-Hochberg FDR) was used. The differentially expressed mRNAs were then visualized on a heatmap (Fig. 7).
We observed the modulation of 311 downstream mRNAs, of which, 103 were induced and 208 were repressed. Of the 311 modulated mRNAs, 249 were annotated in the DAVID genome database platform. Of these, 58 (23.3%) encode TRAs, which is consistent with Sansom Sans et al. [1], which showed that approximately 25% of Aire-dependent genes encode TRAs. Comparing our expression data from miR-155 transfected cells with data from Aire KO mice described by Sansom et al. [1], of the 58 differentially expressed TRAs (this work), 35 were differentially expressed of which 29 presented similar expression profile ( Table 1). The definition of TRA mRNAs was as previously described [1], i.e., an mRNA encoding a TRA should have its expression increased in a maximum of five tissues compared to the average expression levels. As miRNAs have several mRNA targets, in this work, we evaluated which of the mRNAs with altered expression could be targeted by miR-155 according to MirWalk predictions ( Table 2).

Functional enrichment of differentially expressed mRNAs
To identify the effect of miR-155 on mTEC 3.10 cells, in silico predictions of the functions of the differentially expressed mRNAs was performed using the DAVID genome database platform (https:// david. ncifc rf. gov/). Initially, the mRNAs were separated into two groups, namely, the upregulated and downregulated groups, and they were grouped based on their functions in biological processes, cellular components, or molecular function. Categories that included at least three mRNAs and a score < 0.05 with Benjamini-Hochberg correction were further considered (Fig. 8).

Thymocyte migration
To evaluate the effect of the miR-155-5p mimic on the chemotactic properties of mTEC 3.10 cells, we performed a thymocyte migration assay since five of the downregulated mRNAs (Erbb4, Efna1, Prkcl, Prss37, and Matk) are related to cell migration. Conditioned culture media from the control or miR-155-5p mimic-transfected mTEC cell cultures (24 h posttransfection) were used. It was possible to observe a reduction of approximately 34.7% in the number of thymocytes that migrated through the inserts when cultured with conditioned culture medium from transfected cells compared to conditioned culture medium from control cells (Fig. 9).

Discussion
In this study, we investigated the posttranscriptional control of the autoimmune regulator (Aire) gene using an in vitro model system with a murine medullary thymic epithelial cell line. The control of the promiscuous gene expression (PGE) of autoantigens in mTECs is mainly performed by the Aire transcriptional regulator, along with the Fezf2 transcription factor, but to date, little is known about the mechanisms by which this gene is controlled [2,4,36].
Two important axes that control the expression of eukaryotic genes, which are controlled by transcription factors at the chromatin level [37] or posttranscriptionally by miRNAs in the cytoplasm [38][39][40], might act in this case. The miRNAs finely control the expression of proteins by interacting with their respective mRNAs.
Here, we explored mechanisms of posttranscriptional control. We hypothesize that Aire is regulated by miR-155 in mTECs. This hypothesis was generated after we Table 1 Differently expressed tissue-restricted (TRA) mRNAs (fold change ≥1.5). In bold are represented the TRA mRNAs whose expression profiling is shared and asterisk (*) not-shared with mTECs from Aire KO mice (Sansom et al. 2010). Underlined are mRNAs to cell migration process, according to DAVID database

Adamtsl5
Dhrs11 Myl3 Up-regulated TRAs performed in silico prediction of miRNA-Aire mRNA interactions and observed a strong hybridization of miR-155 with the 3'UTR of the mRNA that encodes the AIRE protein (Supplemental Fig. 2). This prediction was validated by a luciferase assay, making it possible to observe the miRNA-mRNA interaction in the intracellular milieu of human HEK-293 T cells. However, we were interested in validating the miR-155-Aire mRNA interaction in the mTEC cell milieu.
We used the murine Aire + mTEC 3.10 cell line [41], which was transfected with the miR-155-5p mimic. In the short term, this transfection led to a reduction in the protein expression of AIRE, as detected by western blotting. The AIRE protein levels returned to normal level, comparable to the levels in the control cells, 48 h after the transfection of the miR-155-5p mimic. The analyses showed that the dose and time course of the miR-155-5p mimic were sufficient to affect the mRNA and protein levels of Aire, which suggests possible consequences for the transcriptome of these cells due to the importance of Aire in PGE in mTECs. The functional categories were identified by using DAVID genome database platform according to GO processes: biological processes (BP), cellular components (CC) and molecular function (MF). The rank is based on the enrichment score, which represents mean p-value. Only those mRNA-groups yielding < 0.05 Benjamini corrected p-value and containing at least three mRNAs are considered to be significant Fig. 9 Decrease of thymocyte migration towards mTEC 3.10 cells (chemotaxis) following miR-155-5p mimic transfection. Differences in the numbers of migrated thymocytes through porous membrane in conditioned medium from transfected or not transfected cells were quantified, and data are shown in bar graph form. The migration degree is plotted as migration index, in which values correspond to mean ± s.d. Data are representative of three independent experiments and were significantly different between control cells versus cells transfected with miR-155-5p mimic. (**P < 0.01, Student t-test) Different time points after transfection of the miR-155-5p mimic, as well as different concentrations of the miR-155-5p mimic, were tested. Twelve hours of transfection with 2 nM miR-155-5p mimic proved to be the best conditions that achieved the highest concentration of intracellular miR-155-5p. At 24 h after the transfection of the miR-155-5p mimic, there was a drastic reduction in the level of Aire transcriptional expression, suggesting degradation of Aire mRNA by the RISC complex.
The levels of miR-155 returned to normal, that is, to levels comparable to those observed in control cultures, 24 h after transfection. Regarding the mRNA levels of Aire, we observed that 48 h after transfection, there was no longer a significant difference in its expression between the control and transfected cultures. This defined a "window" for making posttransfection observations: 12 h for miR-155 or 24 h for Aire mRNA. This window may be a result of the model system, which is analogous to small interfering RNA (siRNA), whose target inhibition is partial and transient [39,42]. However, the levels of overexpression of miR-155 and inhibition of Aire mRNA achieved in this model system were sufficient to answer the central question of the work.
Immunofluorescence showed no change in the number of Aire + mTECs after transfection with the miR-155-5p mimic, but a reduction in the intensity of AIRE protein expression in each cell was observed. This result is consistent with the interpretation of an inhibitory effect of miR-155 on the translation of Aire mRNA rather than an elimination of Aire + cells.
After observing the reduction in the mRNA and protein expression of Aire after the transfection of the miR-155-5p mimic, we proceeded to the analysis of the gene expression of mTEC 3.10 cells on a large scale, considering the impact of Aire and miR155-5p itself on the transcriptome of these cells.
Total RNA extracted from mTECs 24 h after transfection or from control cells was used to analyze gene expression by the microarray method.
In total, 311 mRNAs exhibited altered expressions; of these, 103 were upregulated, and 208 were downregulated. Of these 311 mRNAs, 249 encode previously described proteins. However, the other 62 corresponded to predicted mRNAs that have not been recorded in public databases.
Interestingly, of the 58 TRA mRNAs, ten were overexpressed in mTEC 3.10 cells transfected with the miR-155-5-mimic. These mRNAs must represent a portion of genes that, instead of being induced, are suppressed by Aire at the transcriptional level [1,5].
Considering that miRNAs, because they are small single-stranded RNAs, can have several targets, we evaluated which of the differentially expressed mRNAs might be targets of miR-155 are as predicted by the MirWalk algorithm.
The results showed that the reduction in Aire expression caused by transfection with the miR-155-5p mimic was sufficient to cause the modulation of TRA mRNAs, which are essential for self-representation and negative selection in the thymus [2,43].
Many studies have shown that the importance of Aire in mTECs extends beyond the control of TRA expression during PGE [1,2,5,6,16,17]. The functional enrichment of the differentially expressed genes allowed us to identify mRNAs related to cell migration.
Pezzi et al. [16] showed that Aire controls the expression of mRNAs that encode proteins associated with cell adhesion (mTEC-thymocyte interaction), a process that is essential for negative selection. Among the downregulated mRNAs, we identified Cldn6 and Hapln3, which are both related to the cell adhesion process.
St-Pierre et al. [5] demonstrated that the deletion of Aire in cTECs or in mTECs low , that is, cells that do not express Aire, does not alter the complexity of splicing in these cells. Thus, AIRE can also participate in mRNA processing and contribute to isoform diversity of both Aire-dependent and Aire-independent TRA genes [17,44]. In this study, we identified mRNAs that encode proteins that participate in spliceosome assembly, such as SNRPN, PEPF39, SF3A2, and SNRPA.
Another possible role of Aire is in inducing the apoptotic process in mTECs. Although this function has not yet been fully confirmed, the induction of apoptosis by Aire would explain the rapid turnover of Aire + mTECs and suggests that after the apoptotic process, these cells would be phagocytosed and that TRAs would be presented by antigen-presenting cells present in the thymic medulla, such as macrophages or dendritic cells, which assist in the negative selection process [2,6,45,46].
Among the downregulated mRNAs, six were identified as being related to MAPK (mitogen-activated protein kinase) activation, such as Erbb4, Rreb1, Dmtf1, Trp73, Tbl1x, and Brca1. MAPKs regulate a diversity of biological processes, such as cell differentiation and apoptosis.
The mRNAs encoding ADAMTSL5 and ARSI proteins were also inhibited, and these mRNAs participate in protein metabolism and macromolecule degradation, respectively, which are crucial processes in mTECs since these cells express a vast array of ectopic proteins.
We also observed the repression of mRNAs related to the cell migration process, including Erbb4, Efna1, Prkcl, Prss37, and Matk. Laan et al. [47] showed that Aire increases the expression of chemokine ligands involved in the migration of thymocytes, such as CCR7 and CCR4, strongly suggesting a role for Aire in the chemotactic properties of mTECs and possibly influencing, even indirectly, the migration of thymocytes to the thymus and within the thymus [13].
This prompted us to analyze the biological (functional) effect caused by transfection with the miR-155-5p mimic. The rationale was that if the transfection of miR-155-5p mimic causes a reduction in the expression of mRNAs that encode migration-related proteins, the chemotactic properties of the treated mTECs would be altered.
To test this hypothesis, we performed a thymocyte migration assay. The assay was carried out using inserts from the transwell culture plates that allow thymocytes to pass towards conditioned culture media (media from cultures of control mTECs or miR-155-5p mimic transfected mTECs). The transwell migration process was stopped after three hours to fix and stain the thymocytes in the inserts during the migration process and thus evaluate the migration rate.
Transfection of mTECs with the miR-155-5p mimic effectively reduced the mRNA and protein expression of Aire. The consequence of this was significant for the transfected cells, as we found hundreds of downstream mRNAs that were also reduced due to the reduction in Aire. Of the 311 reduced downstream mRNAs, only 15 are predicted targets of miR-155.
The results showed a significant reduction in thymocyte migration in the conditioned medium compared with cells transfected with the miR-155-5p mimic.
Recent findings by [48] demonstrate that miR-155 safeguards mTEC maturation by interfering with the TGFβ signaling pathway and that mature mTECs in the thymic medulla participate in the T regulatory (Treg) differentiation. Morimoto et al. [49] argued that impaired Tregs production is associated with the development of Airedependent autoimmunity because these cells develop under T-cell receptor (TCR)-TRA recognition. In the lack of Aire in mTECs, the TRA expression is impaired and does not stimulate Tregs accurately. These authors demonstrate that under Aire deficiency, mTECs abnormally express CTLA4 co-stimulatory molecule, which removes the CD80/CD86 ligands from dendritic cells, rendering them incapable of providing the differentiation of Tregs. Considering that miR-155 can: 1) posttranscriptionally control Aire in mTECs (this work), 2) control the TGFβ signaling pathway [48], and 3) that Aire deficiency is associated with impaired production of Tregs [49], we hypothesized that miR-155 would occupy a hierarchical level of control even more upstream in mTECs, either controlling TGFβ signaling pathway [48] or controlling Aire. These findings open the perspective for future studies to understand better the possible dual role of miR-155 in mTECs that can either favor or decrease the production of Tregs, depending on the target at a given moment of mTEC development.

Conclusions
This work showed that Aire might be regulated at the posttranscriptional level by miR-155. Data confirms, at the biological level and using an in vitro model system, that the decrease in the expression of Aire and consequently in the expression of mRNAs that encode chemotactic proteins might have functional consequences regarding the chemotactic property of mTECs.

In silico predictions of miRNA-mRNA interactions
For the analysis of possible miRNAs that interact with Aire mRNA, a search was performed in public databases to predict miRNA-mRNA interactions using the Mir-Walk platform, which is available online at (http:// zmf. umm. uni-heide lberg. de/ apps/ zmf/ mirwa lk2/ miRre tsysself. html). The interaction of miR-155 with the 3'UTR of Aire mRNA was predicted using the PITA algorithm, which is available at https:// genie. weizm ann. ac. il/ pubs/ mir07/ mir07_ predi ction. html.
The thermodynamic stability of the molecular interaction was then validated using the RNAhybrid algorithm (https:// bibis erv2. cebit ec. uni-biele feld. de/ rnahy brid), which analyzes the secondary structure of the interaction and determines the most favorable hybridization site between a determined miRNA and the exact sequence of the 3'UTR of its target mRNA. We only selected interactions with Gibbs free energy ≤ − 20 kcal/mol, which are the most stable.
The algorithms of the MirWalk platform were used to predict other mRNAs, in addition to Aire mRNA, which could be targeted by miR-155.

Validation of miRNA-mRNA interactions through luciferase reporter gene assay (LRGA)
In this study, we used an oligonucleotide encompassing The oligonucleotide was cloned into the polycloning site of the pmirGLO vector (Promega Corporation) between the XhoI/XbaI restriction sites, resulting in the miRNA target region assuming the correct 5′ to 3′ orientation immediately downstream of the luciferase gene.
For the selected target, we introduced mutations along the 3'UTR sequence. The constructs, named "pMIR-Aire-wt-3'utr" to indicate the inclusion of the wild-type sequence and "pMIR-Aire-mut-3'utr" to indicate the inclusion of the mutant sequence, were selected by colony polymerase chain reaction (PCR) using a pair of primers flanking the vector polycloning site. We used Escherichia coli DH5α for cloning.
For the LRGA, 0.2 μg of each pmirGLO construct was transfected into human HEK-293 T cells [50] (6 × 10 4 cells/well) with 1.6 pmol of miR-155 or scrambled miRNA mimic negative control (both from Thermo Scientific Dharmacon, Waltham, MA) in a 96-well plate. Transfections were performed using the Attractene Transfection Reagent (Qiagen, Hilden, Germany) according to the manufacturer's instructions. The transfected cells were incubated at 37 °C in a 5% CO 2 incubator; 24 h after transfection, the cells were lysed in Passive Lysis Buffer. The firefly and Renilla luciferase activities were measured in a Synergy 2 luminometer (BioTek Instruments, Inc., Winooski, VT) using the Dual-Luciferase Reporter System (Promega Corporation) according to the manufacturer's instructions.
The LRGA results are presented as the standard error of the mean (SEM). The differences were evaluated by one-way ANOVA followed by Student's t-test (two groups). P < 0.05 was considered to be statistically significant.
The miR-155-5p mimic-lipofectamine mixture was prepared by mixing 1 μL of 10 μM miR-155 mimic with 500 μL of Opti-MEM I Reduced Serum Medium Gibco (Paisley, Scotland, UK). Next, 10 μl of Lipofectin Reagent was added. This mixture was vigorously stirred for 10 seconds and then incubated for 10 minutes at room temperature.
The miR-155-5p mimic-lipofectamine complex was added to the culture flasks in a total volume of 5 ml per bottle, resulting in a concentration of 2 nM miR155 mimic. For the control, mTEC 3.10 cells were cultured under similar conditions, except that the miR-155 mimic was added alone. The cells were cultured for different periods (12,24, and 48 h posttransfection). Independent experiments and controls were repeated three times.
Total RNA was extracted from control or transfected cells using the mirVana miRNA Isolation Kit (Ambion, Grand Island, NY, USA) according to the manufacturer's instructions. The integrity of the RNA samples was evaluated through microfluidic gel electrophoresis in RNA 6000 nanochips and in an Agilent 2100 Bioanalyzer device (Agilent Technologies, Santa Clara, CA, USA). We used only RNA preparations with an RNA integrity number (RIN) ≥ 7.5 that were free of phenol (A 260 /A 220 approx. 0.7) and free of protein (A 260 /A 280 approx. 2.0).

Reverse transcription for cDNA synthesis
mRNA reverse transcription reactions were performed with components of the Applied Biosystems (Foster City, CA, USA) High Capacity cDNA Reverse Transcription Kit following the manufacturer's protocol. For these reactions, 500 ng of total RNA was used. For the reverse transcription of miR-155, the TaqMan MicroRNA Reverse Transcription Kit (Foster City, CA) was used following the manufacturer's protocol. For these reactions, 50 ng of total RNA was used.

Real-time quantitative PCR (RT-qPCR)
The expression levels of Aire mRNA and miR-155 were assessed by reverse transcription quantitative real-time polymerase chain reaction (RT-qPCR). For these reactions, TaqMan Gene Expression Assay probes specific for Aire, Gapdh (endogenous control for mRNA), miRNA-155, and snoRNA-202 (endogenous control for miRNA) were used. The TaqMan Fast Advanced Master Mix kit (Applied Biosystems) was used following the manufacturer's protocol.
The expression level of miR-155 was normalized to that of snoRNA202, which is commonly used as an endogenous control for miRNA expression analysis.
The expression levels of the Aire gene were normalized to those of the constitutive gene Gapdh (ENS-MUSG00000057666), which is commonly used as a reference.
The relative quantification of gene expression was performed using the 2-ΔΔCT method [53]. For the statistical analysis of the data, Student's t-test and one-way ANOVA were performed using GraphPad Prism statistical software (www. scola ry. com/ tools/ graph pad-prism) according to the manufacturer's instructions.

Extraction and quantification of total proteins from mTEC 3.10 cells
Approximately 5 × 10 5 cells were seeded in T75 culture polystyrene flasks containing 10 mL of RPMI 1640 medium and 10% fetal bovine serum (FBS). The cells were transfected with the miR-155-5p mimic and cultured for 24, 36, or 48 h. The cells were detached from the culture flasks by the conventional trypsin treatment. The cells were washed three times with 1x PBS and centrifuged at 1000 x g for 5 minutes. Then, 300 μL of modified RIPA buffer [50 mM TrisCl pH 7.4, 150 mM NaCl, 1 mM EDTA, 0.1% SDS, 0.1% protease inhibitor cocktail (Sigma-Aldrich, Darmstadt, Hesse, Germany)], and the cells were lysed using a Polytron homogenizer. The homogenate was then centrifuged at 10,000 x g for 20 minutes at 4 °C. The total protein content of the supernatant was determined by using conventional Bradford reagent (Bio Rad, Hercules, CA, USA).

Western blot analysis of the AIRE protein
The electrophoresis of total proteins in denaturing polyacrylamide gels (SDS-PAGE) was performed using the discontinuous buffer system described by Laemmli and Favre [54]. Polyacrylamide mini gels (10 × 8 cm; 0.075 cm thick) were used. The concentration of the stacking gel was 4%, and the concentration of the separation gel was 10%. The protein samples were prepared for electrophoresis at a ratio of 3:1 (vol/vol) of sample to 4x concentrated sample buffer (240 mM Tris-HCl pH 6.8; 8% SDS; 40% glycerol; 20% 2-mercaptoethanol and bromophenol blue). The samples were heated for 5 minutes at 95 °C in a heat block. The electrode buffer contained 25 mM Tris-HCl, pH 8.3; 192 mM glycine; and 0.1% SDS. The protein molecular weight standard Precision Plus Protein ™ Dual Color Standards (BioRad, Hercules, CA, USA) was used.
The protein bands, once separated, were transferred from the gel to a PVDF membrane (BioRad, Hercules, CA, USA) using 25 mM Tris-HCl, 190 mM glycine, 20% methanol, and 1% SDS as transfer buffer at 80 V for 2 hours. For this step, we used the BioRad Mini Trans-Blot Cell device (Hercules, CA, USA). After the transfer, the membrane was stained with 0.2% Ponceau S in 3% trichloroacetic acid for 5 minutes to evaluate the transfer efficiency. Then, the membrane was destained with deionized water. The membrane was blocked with skim milk (MOLICO) diluted to 5% in TBS-Tween (50 mM Tris-HCl, pH 8.0; 150 mM NaCl; 0.1% Tween-20) for 1 hour at room temperature. After blocking, the membrane was incubated with the polyclonal anti-Aire D17 (goat antibody) primary antibody from Santa Cruz Biotechnology, Inc. (Dallas, Texas, USA), which was diluted in TBS-T (1:500), for 18 h, followed by three washes in TBS-T for 5 minutes each.
The membrane was incubated with the anti-goat IgG secondary antibody conjugated to peroxidase (horseradish-peroxidase -HRP), which was diluted 1:10,000 in TBS-T, for 1 hour. After incubation, the membrane was washed in TBS-T, as before. The development was performed in an ImageQuant LAS 500 Ge Life Sciences apparatus (Piscataway, NJ, USA) using the appropriate substrate of peroxidase (Luminata Forte, Merck Millipore). The reaction was conducted in the dark, and we exposed the membrane for different time intervals. The conversion of the intensity of the protein bands in the images into numerical values was performed with the aid of the ImageJ 1.49 program (http:// imagej. nih. gov/ ij/ index. html). The GAPDH (glyceraldehyde 3-phosphate dehydrogenase) protein was used as an internal reference in the analyses. Thus, the membrane was also incubated with the polyclonal (rabbit) anti-GAPDH primary antibody from Cell Signaling Technology (Beverly, MA, USA). In that case, a secondary anti-rabbit IgG antibody conjugated to peroxidase was used. The immunostaining procedure was the same as that used for Aire.

Immunolocalization of the AIRE protein
To assess the intracellular localization of the AIRE protein, we used the immunofluorescence of mTEC 3.10 cells (control or transfected with the miR-155 mimic).
For this experiment, we used the anti-AIRE D17 primary antibody from Santa Cruz Biotechnology (Dallas, Texas, USA). For this assay, 3 × 10 4 mTECs (3.10) were plated on 13-mm 2 glass coverslips placed in a 24-well plate. After transfection with the miR-155-5p mimic, the culture medium was removed from the wells, and the coverslips were washed five times with PBS.
The cells were fixed on the coverslips with 300 μL of 4% paraformaldehyde in PBS for 15 minutes. The cells were then subjected to three washes with PBS for 5 minutes each and then permeabilized with PBS containing 0.5% Triton X-100 for 5 minutes.
After three more washes with PBS (5 minutes each), the cells were blocked for 1 hour with 2% BSA in PBS. A diluted (1:50 in PBS containing 1% BSA) primary antibody (anti-AIRE D17) was added, followed by incubation at room temperature for 1 hour. The cells were washed three times with PBS (10 minutes each). Then, the Novex TM mouse anti-goat IgG (H + L) rhodamine red-conjugated secondary antibody (Life Technologies Corporation, Carlsbad, CA, USA), which was diluted 1:500 in PBS containing 1% BSA, was added and incubated for 1 hour. The coverslips were washed with PBS (3 × 10 minutes).
To visualize the cytoplasmic region, the actin filaments were labelled with phalloidin conjugated to AlexaFluor 488 (Life Technologies; Eugene, Oregon, USA) according to the manufacturer's instructions. The coverslips were mounted with ProlongGold Antifade Mountant (Life TechnologiesTM -Eugene, Oregon, USA) containing DAPI (to visualize the nuclei). The coverslips were visualized and photographed using an ApoTome Zeiss fluorescence microscope (Oberkochen, Germany).

Microarray hybridizations and data analysis
This study followed a protocol that was previously established in our laboratory for microarray hybridizations [30]. Briefly, total RNA preparations were used to synthesize dscDNAs from mRNAs, and then, cyanine 3 (Cy3)-CTP-labeled complementary amplified RNAs (cRNAs) were synthesized. For these reactions, we used the Agilent Linear Amplification kit (Agilent Technologies) according to the instructions of the manufacturer. The (Cy3)-cRNA samples were hybridized to Agilent mouse 4 × 44 K-format oligonucleotide microarrays (Agilent Technologies). The slides were washed according to the manufacturer's instructions and scanned with an Agilent DNA microarray scanner.
The hybridization signals from the scanned microarray slides were extracted using Agilent Feature Extraction software, version 10.7.1.1. The expression profiles of independent samples of the control or miR-155-5p mimictransfected mTEC 3.10 cells were analyzed by comparing the microarray hybridization signals of the respective samples. The numerical and quantitative microarray data were normalized to the 75th percentile and analyzed in the R statistical environment (version 3.3.1) (https:// www.r-proje ct. org). In the data preprocessing phase, we used the arrayQualityMetrics [55] and Agi4x44PreProcess [56] tools, which have algorithms that allow the qualitative evaluation of arrays in addition to the correction of the background and normalization of the data.
For the analysis of differentially expressed mRNAs, we used the functions of the Limma package [57], which applies a linear model to the gene-wise statistical analysis. For an analysis of multiple tests, empirical Bayes and the Benjamini-Hochberg correction were applied.
In this work, we considered a p-value ≥0.05 with correction by FDR (Benjamini-Hochberg) and fold-change ≥1.5 to indicate differentially expressed mRNAs.
The differentially expressed transcripts were then subjected to hierarchical clustering, and heat maps were constructed to evaluate the expression pattern of the mRNAs. The Euclidean distance and the complete linkage method were used to group the samples and the RNAs. The functional enrichment of the differentially expressed (DE) mRNAs was done through the Database for Annotation, Visualization, and Integrated Discovery (DAVID) annotation tool (https:// david. ncifc rf. gov/). This tool was used for the identification of the main biological processes and pathways represented by DE mRNAs. A functional category was considered significant if it had at least three mRNAs and a score of p < 0.005 with Benjamini-Hochberg correction. Tissuerestricted antigen (TRA) mRNAs were defined as Samson et al. [1]. Briefly, an mRNA encoding a TRA was identified according to its higher expression in up to five different tissues concerning the expression mean of the whole set of the differentially expressed mRNAs.

Thymocyte migration assay
For the separation of thymocytes, female C57BL/6 mice aged 5-6 weeks old were killed in a CO 2 chamber, and the thymuses were quickly removed by thoracic surgery. The organs were fragmented with the aid of two surgical forceps in RPMI 1640 culture medium in a Petri dish. The thymocytes were separated through a 10 μm pore size nylon membrane (Sefar Inc. Depew, NY, USA). The thymocyte suspension was centrifuged for 5 minutes at 1000 x g, and the cells were washed twice with PBS and then resuspended in RPMI 1640 medium.
For the thymocyte migration assay, we used a protocol previously published [58] with modifications as follow. Transwell chambers and 6.5 mm inserts with a 0.5 μm in diameter polycarbonate membrane (Corning Inc. Corning, NY, USA) were used. The inserts were placed in 24-well plates and incubated in PBS for 1 hour at 37 °C and then blocked in PBS containing 10 μg BSA/ mL for 1 minute.
The inserts were then transferred to wells containing conditioned culture medium, i.e., medium collected from mTEC 3.10 control cell culture or mTEC 3.10 transfected with miR-155 mimic (cell cultures were grown for 24 hours). The conditioned medium was filtered through 0.45 μ filters before use. 100 μL of conditioned medium from the control cultures or the transfected cultures were placed in each well. A suspension containing 2.5 × 10 6 thymocytes was then deposited on each of the inserts. The inserts were then placed inside the wells and incubated at 37 °C in an atmosphere with 5% CO 2 for 3 hours to allow the migration of thymocytes.
The inserts were then placed in wells containing 70% ethanol for 10 minutes to fix the thymocytes in passing through the pores. The outer surface of each insert was stained with Giemsa for 10 minutes and finally washed in deionized water to remove excess dye and then left to dry in room air. The inserts were visualized and photographed in a stereomicroscope model Stemi 508 Zeiss (Oberkochen, Germany).
The images were analyzed considering the average pixel of the migrated thymocytes (pixels/area). Thymocyte migration ranged from zero (no migration) to 255 pixels, which was the maximum pixel value for migrated and stained thymocytes. Migration values were compared between thymocytes that migrated in conditioned culture medium versus conditioned culture medium of cells transfected with miR-155 mimic.