The absence of the autoimmune regulator gene (AIRE) impairs the three-dimensional structure of medullary thymic epithelial cell spheroids
BMC Molecular and Cell Biology volume 23, Article number: 15 (2022)
Besides controlling the expression of peripheral tissue antigens, the autoimmune regulator (AIRE) gene also regulates the expression of adhesion genes in medullary thymic epithelial cells (mTECs), an essential process for mTEC-thymocyte interaction for triggering the negative selection in the thymus. For these processes to occur, it is necessary that the medulla compartment forms an adequate three-dimensional (3D) architecture, preserving the thymic medulla. Previous studies have shown that AIRE knockout (KO) mice have a small and disorganized thymic medulla; however, whether AIRE influences the mTEC-mTEC interaction in the maintenance of the 3D structure has been little explored. Considering that AIRE controls cell adhesion genes, we hypothesized that this gene affects 3D mTEC-mTEC interaction. To test this, we constructed an in vitro model system for mTEC spheroid formation, in which cells adhere to each other, establishing a 3D structure.
The comparisons between AIRE wild type (AIREWT) and AIRE KO (AIRE−/−) 3D mTEC spheroid formation showed that the absence of AIRE: i) disorganizes the 3D structure of mTEC spheroids, ii) increases the proportion of cells at the G0/G1 phase of the cell cycle, iii) increases the rate of mTEC apoptosis, iv) decreases the strength of mTEC-mTEC adhesion, v) promotes a differential regulation of mTEC classical surface markers, and vi) modulates genes encoding adhesion and other molecules.
Overall, the results show that AIRE influences the 3D structuring of mTECs when these cells begin the spheroid formation through controlling cell adhesion genes.
The three-dimensional (3D) structure of the thymus is composed by two histologically distinct but closely related compartments. The major component of the thymic stromal cells is epithelial cells, which generate the 3D organized meshwork structure. The cortex is mainly formed by cortical thymic epithelial cells (cTECs), exhibiting the CD45−, EpCam+, Ly51+, UEA1− phenotype, and the medulla is primarily formed by medullary thymic epithelial cells (mTECs), which in mice present the CD45−, EpCam+, Ly51−, UEA1+ phenotype [1,2,3]. Other cell types include thymic seeding progenitors, thymocytes at different stages of maturation, including double-negative (CD4−CD8−), double-positive (CD4+CD8+), and single-positive (CD4+CD8− and CD4−CD8+) thymocytes, thymic B cells, macrophages, dendritic cells, Hassall’s corpuscles, and mesenchymal thymic cells, which are required for epithelial cell proliferation, composing the unique 3D structure of the organ .
The thymic 3D architecture offers a particular microenvironment for the positive and negative selection of developing thymocytes [4,5,6]. The thymic cortex is implicated in the positive selection of thymocytes that express functional T-cell receptors (α/β TCRs). The medulla is involved in the negative selection of autoreactive thymocytes, which avidly recognize self-peripheral tissue antigens (PTAs) presented on the surface of mTECs. This unique pattern of gene expression allows mTECs to express more than 19,000 protein-encoding genes, including the “ectopic” genes that encode PTAs, a phenomenon known as promiscuous gene expression (PGE) [4, 7,8,9,10]. At present, no other cell type is known that expresses such a large set of genes [11,12,13]. The positive and negative selection processes are crucial for the induction of central immune tolerance that prevents aggressive autoimmunity [9, 14].
The autoimmune regulator gene (AIRE) is the primary transcriptional regulator of PTAs in mTECs, whose encoded Aire protein unleashes stalled RNA Pol II in the chromatin to proceed with the elongation phase of the transcription [15, 9, 10]. A second controller of PGE in mTECs is the forebrain embryonic zinc finger-like protein 2 (Fezf2) that plays a role as a classical transcription factor and thus binds directly to DNA in specific promoter regions [4, 16].
Given the non-specificity of RNA Pol II in transcribing genes and given that the Aire protein in mTECs unleashes this enzyme on chromatin, AIRE ends up controlling the expression of a diverse range of mRNAs. In addition to PTAs, AIRE also regulates the expression of genes involved in cell adhesion, such as the extracellular matrix (ECM) constituent Lama1, the CAM family adhesion molecules Vcam1 and Icam4, which control mTEC-thymocyte adhesion [17, 18]. Cell adhesion corresponds to an essential biological process in the structure and function of the thymus. The adequate adhesion of TECs increases the efficiency for T-cell development [19, 20]. Besides adhesion molecules, the formation of a dense cellular network is necessary for the maturation of thymocytes, which is composed of ECM proteins, such as laminins, integrins, collagens, and fibronectins, as well as soluble molecules such as hormones, cytokines, chemokines, and growth factors that are also mediated by AIRE .
AIRE −/− mice have small and disorganized medulla of the thymus , and little is known about whether Aire influences the mTEC-mTEC adhesion in a 3D thymus structure. Two experimental strategies enable the study of the in vitro formation of the thymus’s 3D structure that reproduces its microenvironment with the extracellular matrix and thymic epithelial cells. One of these strategies is the thymus re-aggregation [19, 20], and the other is the 3D organotypic culture . In the re-aggregation model, the thymus tissue is devoid of cells (decellularized), but retains most of the microenvironment of the extracellular matrix that could support the re-aggregation between TECs. When transplanted into athymic nude mice, the re-aggregate thymus organoids can receive lymphocyte progenitors derived from bone marrow and develop a diverse and functional T-cell repertoire [19, 20].
In the search for an experimental model that could mimic the thymic microenvironment closest to that found in vivo, Pinto et al.  adapted a 3D organotypic co-culture, preserving the main characteristics of mTECs, such as proliferation and differentiation. This strategy helped to identify molecular components and pathways involved in the mTEC differentiation and promiscuous gene expression. Other studies used stem cells to form “thymospheres,” permitting the study of the factors necessary for thymocyte development [24, 25]. Recently, a 3D culture substrate has been reported that allowed TECs to survive and proliferate, using electrospun fibrous meshes (eFMs) functionalized with fibronectin. The mTECs presented increased proliferation, viability, and protein synthesis when cultured on fibronectin-functionalized eFMs (FN-eFMs) .
In the present study, we asked whether AIRE regulates the adhesion between mTECs during the in vitro formation of a 3D structure. For this, we developed a model system in which mTECs are grown in non-adherent agarose micro-wells and adhere to form spheroids. The comparison between AIRE wild type (AIREWT) versus AIRE−/− mTECs allowed us to evaluate the influence of AIRE in the adhesion between these cells during the spheroid formation.
The absence of AIRE disorganizes the initial phases of spheroid formation
Figure 1A shows the beginning of spheroid growth (0 h), comparing AIREWT mTEC 3.10 vs. AIRE−/− mTEC 3.10E6. The cultures started from the seeding of mTEC cells (0 h) until the complete 3D spheroid formation (24 h). A compact real-time video of the dynamics of the spheroid formation, recorded from 0 to 24 h, is shown in the supplemental material video 1 (available at www.rge.fmrp.usp.br/passos/esferoides). The results show that the mTEC cell lines (AIREWT or AIRE−/−) provide morphologic differences during the development of spheroids. The AIREWT spheroids consolidate and are firmly compact at 12 h of culture, whereas the AIRE−/− spheroids showed a delay in cell-cell adhesion and consolidate only after 20 h of culture. This finding indicates that in the absence of AIRE, mTEC-mTEC adhesion is impaired and the spheroid formation is delayed.
To better characterize the 3D spheroid model system, we constructed a growth curve counting cells at each 12 h of culture. At the first 12 h culture, AIRE−/− spheroids had a cell count higher than AIREWT spheroids. Differently, the development of the AIREWT spheroids drew a classical growth curve, permitting the identification of the exponential, stationary, and decay phases. The slope of the AIRE−/− spheroid growth curve indicates that in the absence of AIRE, the growth is increased at the beginning, and then the growth abruptly decreases. This suggests that in addition to controlling adhesion between cells, AIRE might influence the cell cycle of mTECs (Fig. 1B).
Besides the growth curve, we evaluated the cell viability during spheroid growth. The AIRE−/− spheroids increased viability during the first 12 h and decreased thereafter, whereas the AIREWT progressively reached viability starting from 12 h and maintaining thereafter (Fig. 1C).
AIRE −/− spheroids influence the cell cycle kinetics and the cell death
After 12 and 24 h of adhesion, AIREWT and AIRE−/− spheroids were dissociated and analyzed for the distribution of cells at the different cell cycle phases.
The cell counting performed at 12 h of spheroid formation indicated a significant increase in the percentage of cells at the G0/G1 phase, 46% in AIREWT and 57% in AIRE−/ - cells, showing an accumulation of AIRE−/− cells at G0/G1 (Fig. 2A). Interestingly, a significant reduction was observed in the sub-G1 fraction (7%) for AIRE−/− cells compared to AIREWT (22%). Within 24 h of adhesion, a similar pattern was seen, with a significant increase in the percentage of cells at G0/G1 phase (65%) observed for AIRE −/− cells, as well as a lower percentage (5%) of the subG1 fraction compared to AIREWT (46% G0/G1; 28% SubG1). In addition, a significant increase (18%) in AIRE−/− cells was observed at the G2/M transition compared to AIREWT (13%). Taken together, these results show that the absence of AIRE causes a significant alteration in the cell cycle progression, as demonstrated by the high number of cells undergoing the G0/G1 and G2 phases.
Since AIRE −/− spheroids showed lower percentage of sub-G1 fraction than AIREWT, and since these fractions indicate a certain percentage of cell death, we analyzed the percentage of apoptosis and necrosis. After 12 h of spheroid formation, AIRE−/ - cells showed a significant increase in the population of apoptotic cells (15%) compared to the AIREWT(1.96%) (Fig. 2B). On the other hand, a significant reduction of necrotic cells was observed in AIRE−/− (8.3%) spheroids analyzed after 12 h compared to AIREWT(15.26%). At 0 h (prior to spheroid formation) and 24 h, we did not find significant differences between AIRE−/ - and AIREWT spheroids regarding the percentage and the profile of cell death.
Spheroid dead-cell center and spheroid morphology are impaired in the absence of AIRE
In the center of the spheroids, dead cells were observed for both AIREWT or AIRE−/− cells along the growth; however, AIRE−/− spheroids exhibited an increased rate of dead cells along the growth, maximizing at 48 h, as shown in Fig. 3A. Quantitative analysis allowed us to determine the relative fluorescence intensity for both live (green-colored) and dead (red-colored) cells (Fig. 3B).
The high-resolution optical microscopy of the 1 μm histologic cuts allowed the observation of the internal structure of spheroids. The AIREWT or AIRE−/− spheroids showed significant differences in their areas, especially at the first 24 h of growth. Although AIRE−/− spheroids present a larger area than AIREWT spheroids, AIRE−/− cells are dispersed and compact similarly to AIREWT spheroids only after 24 h of adhesion (Fig. 4A, B).
Scanning electron microscopy (SEM) accessed the detailed external surface of the spheroids. Twenty-four-hour AIREWT spheroids were more compact, exhibiting a well-defined contour, while AIRE−/− spheroids showed an irregular surface (Fig. 4C).
The loss of AIRE impairs the strength of the cell-cell adhesion
To observe the strength of cell-cell adhesion, we induced an enzymatic spheroid disaggregation at 12 h of culture. Compared to AIREWT, AIRE−/− spheroids dissociated faster (Fig. 5A) and exhibited an increased number of disaggregated cells along the 25 min of culture (Fig. 5B), indicating that the lack of AIRE gene decreased the strength and spheroid adhesion.
Phenotypic characterization of spheroid medullary cell surface markers
AIRE WT spheroids showed: i) the characteristic CD45− Ly51− medullary phenotype, observed in 67.1% of mTECs, ii) 86.1% of the cell population expressed EpCAM and UEA-1, iii) 0.38% of cells were double positive for CD80+ MHC-II (Fig. 6A). In contrast, the AIRE−/− spheroids featured: i) the characteristic CD45− Ly51− medullary phenotype in 69.4% of mTECs, ii) 19.4% of the cell population expressed EpCAM+ and UEA-1+, iii) 9.77% of cells were double positive for CD80+ MHC-II+ (Fig. 6A). Comparing the expression of each marker separately, we observed that the AIREWT-, AIRE−/− spheroids exhibited: i) decreased expression of EpCAM, UEA-1, CD45, CD80 molecules ii) a closely similar expression of MHCII and Ly51 molecules (Fig. 6B).
Spheroid RNA-Seq analysis
Comparative transcriptome analyses were performed for AIREWT and AIRE−/− spheroids at 12 h. We identified a set of 1210 differentially expressed (DE) mRNAs, of which approximately 83% (1104) corresponded to protein-coding mRNAs, and about 17% (106) corresponded to non-protein-coding RNAs (Fig. 7A). Among the 1104 DE mRNAs, 606 were upregulated and 498 downregulated (Fig. 7B). Figure 7C shows the hierarchical clustering of the DE mRNAs.
Figure 8 shows the differences observed for the top-ten biological functions of downregulated (Fig. 8A) and upregulated mRNAs (Fig. 8B). Figure 9A shows the differentially expressed genes involved in the positive regulation of epithelial cell proliferation. Figure 9B confirms by RT-PCR the differential expression of some top-ten DE mRNAs, including the downregulated PARVB that encodes a protein associated with the cell adhesion pathway and the PLCB2 and P2RX7 that are involved in the calcium signaling pathway in AIRE−/− when compared to AIREWT spheroids. In contrast, the VEGFC and ID1 transcripts, associated with the regulation of epithelial cell proliferation pathway, were upregulated in AIRE−/− spheroids. We took advantage of the use of the STRING algorithm, which establishes interactions among the encoded proteins based on the validated data retrieved from the literature, to understand the interaction of these genes. Figure 9C shows that the Aire protein may interact with other proteins associated with the immunological response (Foxn1 and CD80), autoimmune disorders (Ptpn22), adhesion receptor (Itgb7), collagen (Col6a1 and Col6a2), and apoptosis (Tfeb).
Considering that: i) the AIRE gene, besides controlling the expression of genes that encode tissue-specific antigens (TSAs) [7, 9, 14], also controls the expression of adhesion molecule genes [17, 18], ii) studies regarding the role of AIRE in the negative selection have primarily focused on the interaction mTEC-thymocyte, and iii) little attention has been devoted to the role of AIRE on the mTEC-mTEC adhesion, in the present study we used 3D spheroids to evaluate the role of the AIRE gene on the mTEC-mTEC adhesion. We showed that the absence of the AIRE gene: i) disorganizes the 3D structure of mTEC spheroids, ii) increases the proportion of cells at the G0/G1 phase of the cell cycle, iii) increases the rate of mTEC apoptosis, iv) decreases the strength of mTEC-mTEC adhesion, v) promotes a differential regulation of mTEC classical surface markers, and vi) modulates genes encoding adhesion and other molecules.
Several strategies have been developed to study the functional thymic properties, including ex vivo reaggregation thymic organ culture (RTOC) [19, 20], 2D mTEC-thymocyte adhesion [17, 18], transwell thymocyte migration , 3D fetal thymus organ culture (FTOC) , and 3D organotypic culture . Considering that the 3D intrathymic environment is necessary for the development of a functional thymus scaffolding, development of thymic epithelial cells, chemoattraction of pro-T-cells, commitment to the T lineage, and generation of a repertoire of T-cell receptors (TCR) and natural T regulatory cells [29,30,31], an adequate thymus structure is mandatory to establish organ functionality.
Since the consequences of AIRE deletion in the intercellular adhesion remain to be elucidated, we studied the interaction between mTECs when these cells were seeded in a non-adherent substrate such as agarose, where 3D spheroids may be observed at the various stages of development. Once in an agarose substrate, mTECs aggregate and adhere to each other to form spheroids, which may be helpful for the study of mTEC-mTEC adhesion. Then, we evaluated the consequences arising from the absence of AIRE, evaluating the: i) mTEC-mTEC physical interaction at regular periods, ii) growth curve of mTEC spheroids, iii) dynamics of 3D mTEC spheroid formation and spheroid dissociation, iv) cell-cycle phases and the rate of apoptosis/necrosis, v) morphology of mTEC spheroids, vi) expression of surface mTEC molecules and vii) gene expression profiles, comparing AIRE WT and AIRE−/− cells.
We reported that AIRE interferes with the early process of spheroid formation, since its absence yielded a delay on the spheroid growth and on cell viability (Fig. 1A-C). Although AIRE is a pro-apoptotic factor in mature mTECs , in this study we observed that during the 3D spheroid formation AIRE−/− spheroids presented increased apoptosis rate. Indeed, the absence of AIRE produced a significant impact on the progression of cells through the cell cycle (Fig. 2A), mainly at 12 h. The higher percentage of cells undergoing the G0/G1 phase in AIRE−/− spheroids are possibly related to the higher percentage of apoptotic cells when compared with AIREWT spheroids (Fig. 2B). These results suggest that the lack of AIRE function may affect the cascade of signaling pathways, leading to alterations in the cell cycle progression (checkpoint mechanism) and apoptosis induction in proliferating cells. Corroborating this idea, we observed that the several genes related to apoptotic processes were upregulated in the AIRE−/− spheroids, including the PERP, PDCD4, TNFRSF21, SLC5A8, S100A9, PPID, CYCS, CARD14, ALDOC, CKAP2, STK17B, POLB, CDCA7, GADD45B, APLP1, CASP1, CASP4, BNIPL, TNS4, AIM2, IL24, STEAP3, FAM3B, SRGN, EPB41L3, EPHA7, FGFR2 and S100A8 genes. Notably, the CDCA7 and GADD45B genes are involved in cell division and mediate Fas-induced apoptosis, respectively.
Regarding the physical structure of the 3D spheroids, dead cells were observed in the center areas of both spheroid types; however, the number of these cells progressively increased along time, predominating at 48 h in AIRE−/− spheroids (Fig. 3A-B). Besides, in the model system of this study, it was possible to observe a distinct morphology of the mTEC interactions; i. e., the AIREWT forms a typical spheroid at 12 h, whereas the AIRE−/− initially forms an ellipsoid form instead of a spheroid shape (Fig. 4A/B). Noteworthy, the scanning electron microscopy showed that AIREWT spheroids are well-compacted with a well-defined contour. In contrast, AIRE−/− spheroids exhibited an irregular shape and ill-defined surfaces (Fig. 4C), indicating that AIRE influences the internal and external spheroid structure.
Since the 3D AIRE−/− spheroids exhibited a delayed formation at 12 h, as shown in Fig. 1A, we also evaluate the strength of spheroid dissociation, using a dispase assay. Indeed, AIRE−/− spheroids dissociated faster than AIREWT spheroids and exhibited increased number of disaggregated mTECs (Fig. 5 A-B). Taken in concert, these in vitro results corroborate the in vivo finding of a disorganized thymic medulla in AIRE−/− mouse  and add information regarding the dynamics of mTEC growth and mTEC-mTEC interactions.
In addition to the spheroid structure, we also evaluate the role of AIRE on the mTEC surface markers (Fig. 6A/B). Both AIREWT or AIRE−/− spheroids maintain their characteristic CD45− Ly51− medullary profile, i.e., confirming that they are not hematopoietic cells (CD45+) nor cTECs (Ly51+). The absence of AIRE decreased the expression of the EpCAM and UEA-1 cells, which are classical markers for epithelial cells and medullary thymic cells, respectively. Noteworthy, the expression of double positive (CD80+ MHCII+) was increased in AIRE−/− spheroids (0.38% in AIREWT to 9.77% in AIRE−/−). Three major stages can be observed along mTEC maturation: i) the immature mTEC presents MHC-IIlow, CD80low and Airelow, ii) during TRAs, mTECs exhibit their mature stage, presenting MHC-IIhigh, CD80high and Airehigh, and iii) at the post-Aire stage, mTECs are MHCintermediary, CD80intermediary, accompanied by an increased rate of apoptotic cells . As observed in Fig. 6A/B, the AIREWT spheroids present features of an immature mTEC constitutively exhibiting low expression of AIRE, as previously reported by our group . In contrast, AIRE−/− spheroids present an intermediate expression of CD80 together with increased rate of dead cells, features that are characteristic of the post-Aire stage. These findings indicate that in the absence of Aire, mTECs present a terminally differentiated mTEC during the first 12 h of 3D spheroid culture, which may affect the interaction with thymocytes.
Considering that AIRE controls more than 3300 genes in mTECs , we initially evaluated the mRNA transcript profiles of AIREWT or AIRE−/− spheroids at 12 h of mTEC-mTEC adhesion, showing a distinct pattern of gene expression (Fig. 7A-C). Among the DE mRNAs, we identified genes associated with cell adhesion, positive regulation of cell proliferation, apoptotic process, and response to hypoxia (Fig. 8A/B). Among the DE mRNAs related to cell adhesion, we observed a set of transcripts that encode proteins belonging to the cadherins (PCDHGB7), collagen (COL15A1), integrin (CIB3), fibronectin (FLRT3), or extracellular matrix protein families (MMP19), which were downregulated in AIRE−/− spheroids and are putatively involved on the dysregulation of the spheroid formation, as observed in this study. In contrast, VEGFC and ID1 genes, involved in the positive regulation of epithelial cell proliferation, were upregulated (Fig. 9A/B), a finding that may be associated with the increased number of cells at 12 h of growth (Fig. 1B). The network generated by the STRING tool (Fig. 9C) revealed that Aire interacts with important proteins involved in the immunological response, autoimmune disorders, adhesion receptor, formation of collagen, and apoptosis, indicating that the lack of Aire impairs the 3D spheroid conformation.
Noteworthy, we observed a set of repressed mRNAs associated with the calcium signaling pathway, and among these, the ADCY3, ADCY1, NOS1, P2RX7, ADRB2, PLCB4, ATP2B4, ATP2A3, P2RX3, CACNA1G, ADRA1B, PLCB1, PLCB2, and F2R transcripts are associated with the intercellular calcium waves. Since intercellular communication used by mTECs is mediated by intercellular calcium waves, requiring functional gap junctions and P2 receptors [34, 35], these genes are potential candidates to be further studied.
Even though operating through different ways, AIRE and mTEC-mTEC adhesion culminate into related activities that regulate a cascade of mRNAs that encodes cell adhesion molecules and other relevant molecules that maintain the 3D spheroid formation, opening new perspectives for studies of molecular mechanisms that control the 3D thymic medulla organization and mTEC adhesion and communication.
While most studies have focused on the mTEC-thymocyte interaction to evaluate the negative selection, this study focused on the role of AIRE on the mTEC-mTEC interaction as the primary step to support the adequate thymic medullary structure. Considering the 3D spheroid model, this study reported that the absence of AIRE disorganizes the 3D structure of mTEC spheroids, promotes a differentially regulation of mTEC classical surface markers, and modulates genes encoding adhesion and other molecules.
Materials and methods
Medullary thymic epithelial cell lines
We use the AIREWT murine (Mus musculus) mTEC 3.10 line (EpCAM+, CD45−, Ly51−, UEA-1+), as previously described [17, 36, 37], and the AIRE−/− mTEC 3.10E6 cell clone. The mTEC 3.10E6 clone obtained by the CRISPR-Cas9 system  was characterized as a carrier of indel mutations affecting both AIRE alleles (compound heterozygous). In the AIRE allele 1, there were two types of mutations: a T > G substitution (mRNA nucleotide, nt, position 351), followed by a nine-bp deletion (GCTGGTCCC, mRNA nt positions 352–360) that transcribed a 1647 nt AIRE mRNA. In allele 2, a single G deletion at mRNA nt position 352 transcribed a 1655 nt AIRE mRNA. As both alleles of the mTEC 3.10E6 clone produced a nonfunctional AIRE protein, they were considered knockout; i. e., AIRE−/−.
We use a precast agarose mold with non-adherent 600 μm diameter microwells, making the mTEC cells adhere once seeded in these compartments. The AIREWT mTEC 3.10 and the AIRE−/− mTEC 3.10E6 cell lines were initially cultured as monolayers in RPMI 1640 medium (Gibco, Darmstadt, Germany) supplemented with 10% inactivated fetal bovine serum in 75 cm2 polystyrene plastic bottles (Corning, New York, NY) in an incubator at 37 °C with 5% CO2 atmosphere. After acquiring confluence, mTECS were trypsinized and seeded in agarose molds, using 2% low electroendosmosis agarose (Sigma-Aldrich, Saint Louis, MO), sterilized in 70% ethanol, washed twice in sterile PBS, and irradiated under germicidal UV light for 15 min. Spheroid growth was observed through a Cytosmart® (Lonza Group AG, Basel, Switzerland) inverted microscope to produce a real-time movie of the culture.
To better characterize the 3D culture model, we quantified the cells that form the spheroids at different time points (determinations at every 12 h). At each time point, the spheroids were removed from the agarose microwells and dissociated using trypsin. Isolated cells were counted using a Cellometer® Auto T4 Bright Field Cell Counter (Nexcelom Bioscience, Lawrence, MA). Triplicates were performed for each time point. We have thus drawn a spheroid growth curve to identify the exponential, stationary growth, and decline phases.
Analysis of cell cycle
Analysis of cell cycle kinetics was performed by flow cytometry using the Guava cell cycle analysis. AIREWT and AIRE −/− spheroids were stained with propidium iodide (Sigma-Aldrich), a fluorescent nuclear marker. After 12 and 24 h of spheroid formation, the samples were collected, dissociated with ACCUMAX (Sigma-Aldrich) and stored at − 20 °C in 70% cold ethanol kept for 18 h. The samples were then resuspended in 200 uL of the propidium iodide (5 μg/mL) solution and maintained for 30 min at 37 °C. The cells were analyzed with a Guava EasyCyte Mini System flow cytometer (Merck Millipore, Burlington, MA) using the Guava CytoSoft software, version 4.2.1 (Merck Millipore). At least three independent experiments were analyzed.
Cell death assay
AIRE WT and AIRE −/− cells (2 × 105 cells/well) prior to spheroid formation (0 h) and spheroids with 12 and 24 h of adhesion were stained with annexin V-PE and 7AAD. After disassociation of the spheroids with ACCUMAX, cells were labeled with the eBioscience™ Annexin V Apoptosis Detection reagents (Invitrogen™, Waltham, MA), according to the manufacturer’s protocol. Analysis of apoptosis and necrosis was carried out on a Guava PCA Instrument using CytoSoft Software. At least three independent experiments were analyzed.
We used the LIVE/DEAD® Viability/Cytotoxicity Kit (Thermo-Fisher, Waltham, MA) to assess the proportion of live and dead cells in the spheroids, following the manufacturer’s instructions.
Histological analysis of spheroids
The spheroids were fixed in 10% formaldehyde buffered in PBS, then dehydrated, passing through a battery of aqueous ethanol solution (75 to 100% ethanol) and included in historesin (HistoResin standard kit, Biosystems Switzerland AG, Muttenz, Switzerland). The resin blocks were cut to 1 μm thickness and deposited on microscope slides. After deparaffinization and rehydration, the slides were stained with hematoxylin-eosin (H&E) for further microscopic examination of the spheroid morphology.
Scanning electron microscopy (SEM)
Spheroids were fixed with 1.6% glutaraldehyde in 0.2 M buffered (pH 7.4) sodium cacodylate overnight at 4°C, then washed three times with 0.2 M sodium cacodylate buffer (pH 7.4). Spheroids were immersed into 1% osmium tetroxide (Sigma-Aldrich) for secondary fixation in 0.2 M buffered (pH 7.4) sodium cacodylate. Fixed spheroids were subsequently dehydrated in 25-100% ethanol and immersed in hexamethyldisilazane (HMDS, Sigma-Aldrich) for 10 min at room temperature. The spheroids were freeze-dried until the HMDS had evaporated, mounted, coated with palladium gold, and examined under a scanning electron microscope (Jeol JSM − 6610 LV, Hitachi High-Tech America, Tokyo, Japan) at an acceleration voltage of 5 kV. The SEM procedures were done at the Electron Microscopy Facility, Ribeirão Preto Medical School, University of São Paulo, Ribeirão Preto, SP, Brazil.
Spheroid disaggregation assay
AIRE WT and AIRE −/− spheroids were cultured and after 12 h of adhesion the spheroids were incubated with dispase (Roche Diagnostics, Basel, Switzerland) for a maximum of 30 min to release cells as monolayers. Every 5 min released monolayers were carefully washed with PBS twice and subject to mechanical stress by pipetting with 1 mL pipettes, and the cells were counted using the Cellometer® Auto T4 Bright Field Cell Counter (Nexcelom Bioscience, Lawrence, MA). Triplicates were performed for each time point.
Phenotyping by flow cytometry
We used a BD FACSymphony™ A1 Cell Analyzer (Becton Dickinson, Franklin Lakes, NJ) flow cytometer to evaluate the medullary phenotype of AIREWTor AIRE−/− spheroids. A sample of 1 × 106 mTECs separated from spheroids were labeled with 1:250 dilution of anti-mouse CD45-PerCP, anti-mouse CD326 (EpCam), anti-mouse Ly51-PE (BD Biosciences, San Jose, CA), anti-mouse CD80-APC, and anti-mouse MHCII-PE antibodies in a final volume of 200 μL of cell suspension. Labeling was also performed with a 1:250 dilution of Lectin agglutinin I (Ulex europaeus UEA-I)-FITC (Vector Labs, Burlingame, CA) and the viability marker DAPI- BD Pharmingen™ .
Total RNA preparation and cDNA synthesis
Total RNA of AIREWT or AIRE−/− spheroids was prepared using the mirVana kit® (Ambion, Austin, TX) according to the manufacturer’s instructions. Microfluidic electrophoresis evaluated RNA integrity using Agilent RNA 6000 nano chips and an Agilent 2100 Bioanalyzer (Agilent Technologies Santa Clara, CA). Only RNA samples that were free of proteins and phenol and had an RNA Integrity Number (RIN) ≥7.0 were selected for cDNA synthesis or RNA sequencing. For cDNA synthesis, the SuperScript® reverse transcriptase enzyme (Invitrogen Corporation, Carlsbad, CA) was used according to the manufacturer’s instructions.
Transcriptome analysis through RNA-Seq
We followed a protocol previously described by St-Pierre et al. . Briefly, paired-end (2 × 150 bp) sequencing was performed by an Illumina HiSeq 2500 sequencer (Illumina, San Diego, CA), using a TruSeq Stranded Total RNA Library Prep Kit (Illumina). The quality of raw FASTQ sequences was first analyzed through a FASTQC program (https://www.bioinformatics.babraham.ac.uk/projects/fastqc/). Then, FASTQ sequences were mapped to the Mus musculus reference genome (mm21) using the STAR 2.5 Spliced Aligner program (https://github.com/alexdobin/STAR), which outputs BAM file containing the sequences and their genomic references and a GTF file with gene annotations used for further determinations of the number of reads per transcript through the HTSeq Count program (http://htseq.readthedocs.io). For each RNA sample analyzed, we recovered a list of genes and their respective number of transcripts that served as input for determinations of the differentially expressed (DE) mRNAs through the EdgeR package (https://bioconductor.org/packages/release/bioc/html/edgeR.html) within the R platform (https://www.r-project.org). EdgeR calculates the fold change (FC) for each mRNA, considering a contrast matrix for a given experimental condition. This study defined WT as a contrast variable and as DE the mRNAs exhibiting a p-value < 0.05 and a false discovery rate (FDR Benjamini–Hochberg correction) FC ≥ 2.0. The DE mRNAs were hierarchically clustered, and a heat map was constructed to evaluate the expression profiling. The raw RNA-seq data from this work is available on Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) under accession number PRJNA763914.
Functional enrichment of differentially expressed mRNAs
The list of the DE mRNAs was analyzed in terms of functional enrichment through the Database for Annotation, Visualization, and Integrated Discovery (DAVID) annotation tool (https://david.ncifcrf.gov/) to the identification of the main biological processes and pathways represented by DE mRNAs. A functional category was considered significant if it comprised at least five mRNAs and a score of p < 0.005 with Benjamini–Hochberg correction. The STRING algorithm (Search Tool for the Retrieval of Interacting Proteins database: https://string-b.org/cgi/input?sessionId=brAz10J8qufA&input_page_show_search=on), which establishes interactions among the gene encoded proteins based on the validated data retrieved from the literature, was used to construct networks.
Reverse transcription quantitative real-time PCR (RT-qPCR)
The validation of the transcriptional expression of selected DE mRNAs was assayed by reverse transcription-quantitative real-time PCR (RT-qPCR) of the respective cDNAs. The expression level of each target mRNA was normalized to the housekeeping mRNA Hprt, which is commonly used as a reference. The Primer Blast (https://www.ncbi.nlm.nih.gov/tools/primer-blast/) web tool was used to select pairs of oligonucleotide primers spanning an intron/exon junction with an optimal melting temperature of 60 °C. The respective sequences were retrieved from the NCBI GenBank database (https://www.ncbi.nlm.nih.gov/). The forward (F) and reverse (R) primer sequences (presented in the 5′–3′ orientation) were as follows:
Hprt (NM_013556.2) F = 5' GCCCCAAAATGGTTAAGGTT 3' R = 5' CAAGGGCATATCCAACAACA 3', Plcb2 (NM_177568.2) F = 5'TGGAGTTCCTGGATGTCACG3' R = 5'GCAGGAAGTGGTTGTCTGGA3', Parvb (NM_133167.3) F = 5' TCTTTCTTGGGCAAGTTGGG3' R = 5' CCATTGGAGAGTTGATGGCG3', P2rx7 (NM_011027.4) F = 5'GCACGAATTATGGCACCGTC3' R = 5'TAACAGGCTCTTTCCGCTGG3', Vegfc (NM_009506.2) F = 5'GCTGATGTCTGTCCTGTACCC3' R = 5'ACTGTCCCCTGTCCTGGTAT3', Id1 (NM_001369018.1) F = 5'CCTGAACGGCGAGATCAGTG3' R = 5'AAGTAAGGAAGGGGGACACC3'.
Gene expression was quantified using a StepOne Real-Time PCR System apparatus (Applied Biosystems). The analyses were performed using the cycle threshold (Ct) method, which allows for quantitative analysis of the expression of a factor using the formula 2 − ΔΔCt, in which ΔCt = Ct target gene − Ct of the housekeeping gene Hprt, and ΔΔCt = ΔCt sample − ΔCt. Experiments were performed in three independent replicates.
Depending on the distribution of the variables in each comparison, statistical analyses were performed using the unpaired T or the Mann-Whitney U-tests, considering significant p-values < 0.05.
Availability of data and materials
The raw RNA-seq data from this work is available at Gene Expression Omnibus (GEO) (https://www.ncbi.nlm.nih.gov/geo/) under accession number PRJNA763914 available at (https://www.ncbi.nlm.nih.gov/bioproject/PRJNA763914).
Muñoz JJ, Zapata AG. Thymus Ontogeny and Development. In: Passos G, editor. Thymus Transcriptome and Cell Biology. Cham: Springer; 2019. https://doi.org/10.1007/978-3-030-12040-5_2.
Matsumoto M, Rodrigues PM, Sousa L, Tsuneyama K, Matsumoto M, Alves NL. The Ins and Outs of Thymic Epithelial Cell Differentiation and Function. In: Passos G, editor. Thymus Transcriptome and Cell Biology. Cham: Springer; 2019. https://doi.org/10.1007/978-3-030-12040-5_3.
Abramson J, Anderson G. Thymic epithelial cells. Annu Rev Immunol. 2017;35(2016):85–118. https://doi.org/10.1146/annurev-immunol-051116-052320.
Takaba H, Takayanagi H. The Mechanisms of T Cell Selection in the Thymus. Trends Immunol. 2017;38(11):805–16. https://doi.org/10.1016/j.it.2017.07.010.
Cosway EJ, Lucas B, James KD, Parnell SM, Carvalho-Gaspar M, White AJ, et al. Redefining thymus medulla specialization for central tolerance. J Exp Med. 2017;214(11):3183–95. https://doi.org/10.1084/jem.20171000.
Yoganathan K, Chen ELY, Singh J, Zúñiga-Pflücker JC. T-Cell Development: From T-Lineage Specification to Intrathymic Maturation. In: Passos G, editor. Thymus Transcriptome and Cell Biology. Cham: Springer; 2019. https://doi.org/10.1007/978-3-030-12040-5_4.
Derbinski J, Pinto S, Rösch S, Hexel K, Kyewski B. Promiscuous gene expression patterns in single medullary thymic epithelial cells argue for a stochastic mechanism. Proc Natl Acad Sci U S A. 2008;105(2):657–62. https://doi.org/10.1073/pnas.0707486105.
Abramson J, Giraud M, Benoist C, Mathis D. Aire’s Partners in the Molecular Control of Immunological Tolerance. Cell. 2010;140(1):123–35. https://doi.org/10.1016/j.cell.2009.12.030.
Passos GA, Speck-Hernandez CA, Assis AF, Mendes-da-Cruz DA. Update on Aire and thymic negative selection. Immunology. 2018;153(1):10–20. https://doi.org/10.1111/imm.12831.
Passos GA, Genari AB, Assis AF, Monteleone- AC, Donadi EA, Oliveira EH, et al. In: Passos G, editor. The Thymus as a Mirror of the Body’s Gene Expression. Cham: Thymus Transcriptome and Cell Biology. Springer; 2019. https://doi.org/10.1007/978-3-030-12040-5_9.
St-Pierre C, Brochu S, Vanegas JR, Dumont-Lagacé M, Lemieux S, Perreault C. Transcriptome sequencing of neonatal thymic epithelial cells. Sci Rep. 2013;3. https://doi.org/10.1038/srep01860.
Sansom SN, Shikama-Dorn N, Zhanybekova S, Nusspaumer G, Macaulay IC, Deadman ME, et al. Population and single-cell genomics reveal the Aire dependency, relief from Polycomb silencing, and distribution of self-antigen expression in thymic epithelia. Genome Res. 2014;24(12):1918–31. https://doi.org/10.1101/gr.171645.113.
Abramson J, Goldfarb Y. AIRE: From promiscuous molecular partnerships to promiscuous gene expression. Eur J Immunol. 2016;46(1):22–33. https://doi.org/10.1002/eji.201545792.
Irla M. Thymic Crosstalk: An Overview of the Complex Cellular Interactions That Control the Establishment of T-Cell Tolerance. In: Passos G, editor. Thymus Transcriptome and Cell Biology. Cham: Springer; 2019. https://doi.org/10.1007/978-3-030-12040-5_6.
Giraud M, Yoshid H, Abramson J, Rahl PB, Young RA, Mathis D, et al. Aire unleashes stalled RNA polymerase to induce ectopic gene expression in thymic epithelial cells. Proc Natl Acad Sci U S A. 2012;109(2):535–40. https://doi.org/10.1073/pnas.1119351109.
Takaba H, Morishita Y, Tomofuji Y, Danks L, Nitta T, Komatsu N, et al. Fezf2 Orchestrates a Thymic Program of Self-Antigen Expression for Immune Tolerance. Cell. 2015;163(4):975–87. https://doi.org/10.1016/j.cell.2015.10.013.
Pezzi N, Assis AF, Cotrim-Sousa LC, Lopes GS, Mosella MS, Lima DS, et al. Aire knockdown in medullary thymic epithelial cells affects Aire protein, deregulates cell adhesion genes and decreases thymocyte interaction. Mol Immunol. 2016;77:157–73. https://doi.org/10.1016/j.molimm.2016.08.003.
Speck-Hernandez CA, Assis AF, Felicio RF, Cotrim-Sousa L, Pezzi N, Lopes GS, et al. Aire disruption influences the medullary thymic epithelial cell transcriptome and interaction with thymocytes. Front Immunol. 2018;9. https://doi.org/10.3389/fimmu.2018.00964.
White A, Jenkinson E, Anderson G. Reaggregate thymus cultures. J Vis Exp. 2008;18. https://doi.org/10.3791/905.
Tajima A, et al. Promoting 3-D Aggregation of FACS Purified Thymic Epithelial Cells with EAK 16-II/EAKIIH6 Self-assembling Hydrogel. J Vis Exp. 2016:112.
Mendes-Da-Cruz DA, Stimamiglio MA, Muñoz JJ, Alfaro D, Terra-Granado E, Garcia-Ceca J, et al. Developing T-cell migration: Role of semaphorins and ephrins. FASEB J. 2012;26(11):4390–9. https://doi.org/10.1096/fj.11-202952.
Hubert FX, Kinkel SA, Crewther PE, Cannon PZ, Webster KE, Link M, et al. Aire-deficient C57BL/6 mice mimicking the common human 13-base pair deletion mutation present with only a mild autoimmune phenotype. J Immunol. 2009;182(6):3902–18. https://doi.org/10.4049/jimmunol.0802124.
Pinto S, Stark HJ, Martin I, Boukamp P, Kyewski B. 3D Organotypic Co-Culture Model Supporting Medullary Thymic Epithelial Cell Proliferation, Differentiation and Promiscuous Gene Expression. J Vis Exp. 2015;2015(101):1–7. https://doi.org/10.3791/52614.
Ucar A, Ucar O, Klug P, Matt S, Brunk F, Hofmann TG, et al. Adult thymus contains foxN1- epithelial stem cells that are bipotent for medullary and cortical thymic epithelial lineages. Immunity. 2014;41(2):257–69. https://doi.org/10.1016/j.immuni.2014.07.005.
Sheridan JM, Keown A, Policheni A, Roesley SNA, Rivlin N, Kadouri N, et al. Thymospheres Are Formed by Mesenchymal Cells with the Potential to Generate Adipocytes, but Not Epithelial Cells. Cell Rep. 2017;21(4):934–42. https://doi.org/10.1016/j.celrep.2017.09.090.
Silva CS, Pinto RD, Amorim S, Pires RA, Correia-Neves M, Reis RL, et al. Fibronectin-Functionalized Fibrous Meshes as a Substrate to Support Cultures of Thymic Epithelial Cells. Biomacromolecules. 2020;21(12):4771–80. https://doi.org/10.1021/acs.biomac.0c00933.
Pérez AR, Berbert LR, Lepletier A, Revelli S, Bottasso O, Silva-Barbosa SD, et al. TNF-α is involved in the abnormal thymocyte migration during experimental trypanosoma cruzi infection and favors the export of immature cells. PLoS One. 2012;7(3). https://doi.org/10.1371/journal.pone.0034360.
Anderson G, Jenkinson EJ. Fetal Thymus Organ Culture. Cold Spring Harb Protoc. 2007;2007(8):pdb.prot4808. https://doi.org/10.1101/pdb.prot4808.
Kyewski B, Klein L. A central role for central tolerance. Annu Rev Immunol. 2006;24:571–606. https://doi.org/10.1146/annurev.immunol.23.021704.115601.
Hogquist KA, Baldwin TA, Jameson SC. Central tolerance: Learning self-control in the thymus. Nat Rev Immunol. 2005;5(10):772–82. https://doi.org/10.1038/nri1707.
Brunk F, Michel C, Holland-Letz T, Slynko A, Kopp-Schneider A, Kyewski B, et al. Dissecting and modeling the emergent murine TEC compartment during ontogeny. Eur J Immunol. 2017;47(7):1153–9. https://doi.org/10.1002/eji.201747006.
Sun L, Li H, Luo H, Zhao Y. Thymic epithelial cell development and its dysfunction in human diseases. Biomed Res Int. 2014;2014. https://doi.org/10.1155/2014/206929.
St-Pierre C, Trofimov A, Brochu S, Lemieux S, Perreault C. Differential Features of AIRE-Induced and AIRE-Independent Promiscuous Gene Expression in Thymic Epithelial Cells. J Immunol. 2015;195(2):498–506. https://doi.org/10.4049/jimmunol.1500558.
Da R, Bisaggio C, Nihei OK, Persechini PM, Fundação WS, Cruz O. Characterization of P2 receptors in thymic epithelial cells. Cell Mol Biol. 2017;47(1):19–31 https://www.researchgate.net/publication/12042294.
Nihei OK, Campos de Carvalho AC, Spray DC, Savino W, Alves LA. A novel form of cellular communication among thymic epithelial cells: Intercellular calcium wave propagation. Am J Phys Cell Phys. 2003;285(5 54-5):1304–13. https://doi.org/10.1152/ajpcell.00568.2002.
Hirokawa K, Utsuyama M, Moriizumi E, Handa S. Analysis of the thymic microenvironment by monoclonal antibodies with special reference to thymic nurse cells. Thymus. 1986;8(6):349–60 PMID: 3101237.
Mizuochi T, Kasai M, Kokuho T, Kakiuchi T, Hirokawa K. Medullary but not cortical thymic epithelial cells present soluble antigens to helper T cells. J Exp Med. 1992;175:1601–5. https://doi.org/10.1084/jem.175.6.1601.
We thank Dra. Roberta R. Costa Rosales from the Laboratório Multiusuário de Microscopia Confocal, Ribeirão Preto Medical School, University of São Paulo, Mrs. Maria Dolores and Mr. José Augusto Maulim from the Laboratório Multiusuário de Microscopia Eletrônica, Ribeirão Preto Medical School, and Mrs. Patrícia Viana Bonini Palma and Mrs. Camila Bonaldo, from the Flow Cytometry Laboratory, Ribeirão Preto Blood Center, Ribeirão Preto, SP, Brazil for their technical assistance.
This work was funded by the following Brazilian Research Support Agencies, Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP, grants # 17/10780-4 to GAP and EAD, and # 2019/23448-3 to ACM-C), Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq, grants # 305787/2017-9 to GAP, and 302060/2019-7 to EAD) and Coordenação de Aperfeiçoamento de Pessoal de Nível Superior (CAPES, Finance Code 001).
Ethics approval and consent to participate
Consent for publication
Authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Monteleone-Cassiano, A.C., Dernowsek, J.A., Mascarenhas, R.S. et al. The absence of the autoimmune regulator gene (AIRE) impairs the three-dimensional structure of medullary thymic epithelial cell spheroids. BMC Mol and Cell Biol 23, 15 (2022). https://doi.org/10.1186/s12860-022-00414-9
- AIRE gene
- Medullary thymic epithelial cells
- Cell adhesion
- 3D cell culture