|
|
||||||||

* Department of Immunology and
Department of Biological Structure, University of Washington, Seattle, WA 98195
| Abstract |
|---|
|
|
|---|
| Introduction |
|---|
|
|
|---|
Understanding the basis for this epithelial heterogeneity has been hampered by a lack of knowledge regarding TE differentiation. Although a putative progenitor population has been identified in the fetal thymus (8, 9, 10), their persistence in the postnatal thymus has not been established. Furthermore, the developmental potential of these putative progenitor cells is not known. At one extreme, they could be specified to a thymic fate and have a very restricted developmental potential. Alternatively, they could possess multilineage potential that becomes restricted in response to proximal environmental cues. Understanding how TE differentiation is controlled will be of great assistance in determining the basis for mTEC heterogeneity and TRA expression by mTECs.
The purpose of this study was to define the properties of TE that would help refine testable models of mTEC differentiation. One property that has not been closely examined concerns the mitotic activity of TE. Compared with hematogenous cells in the thymus, TE is relatively radiation resistant; this property allows for the generation of radiation bone marrow chimeras. That TE is relatively radiation resistant has led to the notion that TE is composed predominantly of differentiated postmitotic cells, although the deleterious effect of conditioning regimens used in bone marrow transplantation suggests that epithelial turnover and maintenance may also be impaired by these treatments (11, 12). The extent of turnover among TECs, and mTEC in particular, may indicate the extent to which TEC maturation is a postmitotic event.
To account for mTEC heterogeneity, and the expression of some TRAs, we have previously suggested that some epithelial progenitors in the postnatal thymus are multipotential and retain the ability to follow differentiation programs of other epithelial lineages (13, 14). According to this view, subsets of mTECs would express transcription factors typically associated with multipotential progenitor cells and/or those known to be involved in the specification of other epithelial lineages. The absence of these transcription factors in mTECs would require a reassessment of this model.
Finally, as a potential indicator of mTEC heterogeneity, the expression of TRAs by these cells must be assessed with greater precision. Morphological and immunohistochemical data have shown that expression of individual TRAs is restricted to scattered single cells (4, 5, 15) and small multicellular foci (16) of mTECs. These data support the notion that mTECs are a heterogeneous mixture of cells. In contrast, analyses of TRA expression by mTECs, using large pools of cells, have assumed that TRA-expressing mTECs represent a homogeneous population. Use of large populations imposes an important limitation on interpreting these results, because they could reflect either a transcriptional profile common to a homogeneous population of cells or the summed expression of a heterogeneous mTEC population. Although the former interpretation is currently favored, our recent demonstration that TRAs associated with respiratory epithelium are preferentially expressed by a small subset of mTECs favors the latter interpretation (16). A more detailed analysis of the pattern of TRA expression by mTECs would facilitate efforts to determine the mechanisms controlling this process.
In this study, we report that mTECs in adult thymus are a mitotically active population, raising the possibility that morphological heterogeneity and TRA expression by mTECs could be associated with ongoing epithelial differentiation and turnover. The mTEC population expresses transcription factors that are associated with multipotency, consistent with the possibility that epithelial progenitor cells in thymus have a degree of developmental plasticity. The observation that mTECs also express transcription factors involved in the differentiation of other epithelial cell lineages, lineages that are represented by the spectrum of TRAs expressed within the mTEC population, raises the possibility that the expression of some of these structural molecules by mTECs is regulated in the same manner as in the parent tissue. Finally, RT-PCR analysis of TRA expression by small pools and individual mTECs clearly demonstrated a high degree of heterogeneity within the mTEC population. When analyzed on this scale, the pattern of TRA expression by mTECs appears to be more compatible with developmental regulation rather than the chromosomal proximity and accessibility of TRA-encoding genes.
| Materials and Methods |
|---|
|
|
|---|
Adult B6 mice (510 wk) were obtained either directly from the National Cancer Institute, from Charles River Laboratories, or from our colony maintained at the University of Washington SPF facility. Age-matched 4- to 6-wk-old BALB/c mice for BrdU labeling experiments were obtained from the National Cancer Institute. All mice were handled and used in accordance with protocols approved by the University of Washingtons Institutional Animal Care and Use Committee.
Sorting and flow cytometry
Enzymatic dissociation of thymi was performed as described (16). Briefly, thymi were removed and diced in HBSS (Ca2+- and Mg2+-free; CMF) + 2% FBS + HEPES using forceps and scalpel, washed, and sequentially digested using collagenase D (Roche) in CMF HBSS followed by a mixture of collagenase and neutral dispase (Roche) in CMF HBSS to obtain a single-cell suspension. Single-cell suspensions were treated with rat anti-Fc
R II/III mAb (clone 24G2) before staining. After staining, cells were washed twice in HBSS CMF + 2% FBS and resuspended in the same containing 7-aminoactinomycin D (Molecular Probes) at a final concentration of 2 µg/ml. Cells were sorted on either a BD Vantage or Aria cell sorter directly into Trizol (Invitrogen; for analyzing large populations) or into lysis/reverse transcriptase (RT) buffer in 96-well plates for small pools and single-cell analyses. Analysis of mTEC cell cycle status was performed on an Amnis Imagestream (Amnis) using the same dissociation method and primary Abs. Draq5 (Biostatus) was used to stain DNA.
Primary Abs were: monoclonal anti-murine epithelial cell adhesion molecule (Epcam; clone G8.8) either directly conjugated to Alexa 647 (Molecular Probes) or conjugated to digoxigenin; anti-CD45 PE (clone 30-F11; eBiosciences); anti-MHC class II (clone Y3P, conjugated to digoxigenin; a generous gift of A.Y. Rudensky); a mixture of purified anti-BP-1 Abs (clones 6C3 and CDR1, directly conjugated to Alexa 647), CD4 FITC (BD Pharmingen); and CD8 APC (eBiosciences). Secondary Ab was anti-digoxigenin-FITC (Boehringer Mannheim).
Isolation of RNA
RNA samples in Trizol were vortexed and passed five times each through three (18-, 22-, and 25-gauge) needles. RNA was then isolated using standard Trizol isolation protocol and resuspended in 100 µl of 1x DNase Turbo buffer (Ambion). We added 100 U of SUPERasin (Ambion) to each sample, followed by 5 U of DNase I (Ambion). Each sample was incubated at 37°C for 30 min. The DNase-treated RNA was then reisolated using RNeasy columns (Qiagen), concentrated, and linearly amplified following a modified Eberwine procedure as described in Ref. 17 . From 3 to 5 µg of each amplified cRNA was used to generate dscDNA for PCR analysis with random primers and Arrayscript RT (Ambion) to generate the first strand. PCR analysis was conducted using primers for multiple molecules. All primers were designed using Primer Express (Applied Biosystems) to have a Tm of 60°C and to work well under the following conditions: 25- or 50-µl reactions with a final concentration of 10 mM Tris, 2.5 mM MgCl2, 0.2 mM dNTPs, 0.4 µM primers, and 0.625 U Taq DNA polymerase in buffer B (Promega). All primers, where possible, were designed to span introns. No polymerase controls were used to ensure cDNA specificity for amplicons that do not span introns. All primer sequences are available upon request. Cycling conditions were 94°C for 5 min; cycles of 94°C for 30 s, 60°C for 30 s, 72°C for 45 s, 72°C for 10 min; 4°C hold. For cDNA from large mTEC populations, HPRT, Epcam, aire, and CD45 reactions were amplified through 28 cycles; all other PCR were amplified through 3842 cycles. PCR products were examined on 1.8% Tris-buffered EDTA (TBE) gels stained with Sybrsafe DNA gel stain (Invitrogen and Molecular Probes).
BrdU labeling studies
Age-matched cohorts of BALB/c mice (NCI) were maintained either on BrdU drinking water (0.8 µg/ml; prepared and changed daily) or normal drinking water in the University of Washington SPF facility for up to 5 wk. After 5 wk, experimental groups were returned to normal drinking water for an additional 5 wk. At indicated intervals, 10 mice from each group were sacrificed, and thymi were dissociated to form a single-cell suspension and enriched for epithelial cells by anti-CD45 bead depletion as described in Ref. 18 . After bead depletion, the remaining cells were washed, suspended in FBS, and spun onto aminoalkylsilane-coated slides with a cytocentrifuge. Slides were air dried overnight, fixed in a 75% acetone-25% ethanol for 3 min at room temperature, rehydrated in PBS (pH 7.4), and then stained with primary Abs in PBS containing O.1% BSA + 10% FBS for 1 h at room temperature. Slides were washed then stained with secondary Abs (if biotin or digoxigenin conjugated; all fluorochrome-conjugated reagents and anti-BrdU mAbs were applied after HCl treatment), washed, and fixed at room temperature for 20 min in PBS containing 4% paraformaldehyde. Slides were washed three times in PBS before immersion in 4 N HCl for 30 min at room temperature. Slides were rinsed in PBS and transferred to 0.2 M sodium borate (pH 8.5) for 3 min to neutralize any residual acid. Slides were then incubated in PBS + 0.1% BSA for 5 min, before addition of Alexa 488-conjugated or biotin-conjugated anti-BrdU mAb (Invitrogen and Molecular Probes) and all fluorochrome-conjugated secondary/tertiary Abs and/or staining reagents.
Primary Abs were: G8.8 (purified or digoxigenin conjugate); 6C3 (anti-BP-1; purified or digoxigenin conjugate); rabbit anti-pancytokeratin (CK; Sigma-Aldrich); and biotin- or Alexa 488-conjugated anti-BrdU mAb (Invitrogen Life Technologies and Molecular Probes). Secondary/tertiary Abs and reagents were: anti-rat Ig Alexa546 conjugate; anti-rat Ig Alexa 488 conjugate; streptavidin-Alexa 546 conjugate (Molecular Probes); anti-rabbit biotin conjugate (Sigma-Aldrich), anti-digoxigenin FITC conjugate (Roche).
Generation of pooled cDNAs from sorted mTECs
Generation of cDNA from pools of mTECs sorted from adult (6- to 12-wk-old B6 mice) thymus for TRA expression profile analysis was conducted according to the procedures detailed in Ref. 19 . Briefly, cells in
0.5 µl were mixed in 4 µl of lysis/RT buffer (0.6 mM dNTP mix, 1x first strand buffer, 0.13 µM first primer (ACTCACTATAGGGAGGCGGT24), 0.5 U/reaction SUPERasin RNase inhibitor (Ambion), 0.5 U/reaction Ribolock RNase inhibitor (Fermentas), 0.5% Nonidet P-40) heated to 70°C for 3 min 45 s, and cooled on ice for 5 min before addition of 0.5 µl of a 1:1 mixture of lysis/RT buffer-SuperscriptII RT (200 U/µl; Invitrogen Life Technologies). Plates were vortexed briefly, spun down, and incubated at 37°C. After 30 min, samples were heated to 70°C for 10 min to inactivate the RT enzyme and chilled on ice. Next, 5 µl of 2x TdT mix (6 mM dATP, 3 mM CoCl2, 2x tailing buffer + 25 U/reaction TdT (Roche) was mixed into each reaction and incubated at 37°C. After 15 min, samples were heated to 65°C for 15 min to inactivate TdT. Three parallel PCRs were run consisting of 3.3-µl aliquots each of the TdT reactions added to 21.7 µl of primary PCR mix (0.9 mM dNTP mix, 10 µM first primer, 1x HotStarTaq buffer, 0.5x Q solution containing 2.5 U HotStarTaq (Qiagen), and 0.125 U Pfu Turbo (Stratagene) per reaction), mixed, and run on a thermocycler using the following program: 94°C for 15 min; 1 cycle of 94°C for 45 s; 50°C 2 min, 72°C for 3 min and 15 s; 36 cycles of 94°C for 45 s; 68°C for 45 s; 72°C for 3 min; 72°C for 10 min; 4°C hold. Primary PCRs were combined and purified using Qiaquick PCR purification columns (Qiagen) and eluted in a total volume of 100 µl. Secondary PCR were prepared using 40 µl of purified primary PCR products mixed with 60 µl of secondary PCR mix (0.2 mM dNTP mix, 5 µM second primer (ACTCACTATAGGGAGGCGGTTT), 1x HotStarTaq buffer, 0.5x Q solution, 2 U Taq (Promega) and 0.25 U Pfu Turbo per reaction) and cycled with the following program: 94°C for 5 min; 32 cycles of 94°C for 30 s; 55°C for 45 s; 72°C for 3 min; 72°C for 10 min; 4°C hold. Secondary reactions were column purified and eluted in 100-µl volumes. PCR analysis for multiple genes were performed using 2 µl of purified secondary PCR products in standard 25-µl PCR mixtures with 60°C annealing temperatures (38 cycles for all genes tested). Gene-specific PCR products were visualized on 1.8% TBE gels stained with Sybrsafe DNA stain.
Nested gene-specific PCR on pooled and individual sorted MHC IIhigh mTECs was performed using the method described in Refs. 20 . Briefly, MHC IIhigh Epcam+CD45 cells from adult (6- to 12-wk-old B6 mice) thymus were sorted directly into 4.5 µl of lysis/RT buffer (1x RT buffer, dNTP, 12 µM each gene-specific RT primer + RNase inhibitors + 0.5% Nonidet P-40) on ice. After collection, samples were heated to 70°C for 3 min for 45 s and then cooled to 4°C for 3 min before addition of RT enzyme (0.5 µl of SSII-buffer, 1:1). Samples were mixed and incubated at 42°C for 90 min in a thermocycler. After RT reactions were complete, plates were transferred to ice for 5 min before addition of 45 µl of the outer nest PCR mixture (1x PCR buffer, 2.5 mM MgCl2, 0.24 µM gene-specific outer nesting primers, +2.5 U Taq per reaction). Thermocycler conditions for the outer nesting PCR: 94°C for 5 min; 36 cycles of 94°C for 30 s; 60°C for 30 s; 72°C 1 min; 72°C extension for 10 min; 4°C hold. Inner nesting PCR were performed using 1 µl of the outer nest reaction added to 49 µl of individual inner nesting PCR mixture (1x PCR buffer, 2.5 mM MgCl2, 0.4 µM gene-specific inner nesting primers, 1.25 U Taq). Thermocycler conditions for the inner nesting PCR: 94°C for 5 min; 36 cycles of 94°C for 30 s; 60°C for 30 s; 72°C for 45 s; 72°C extension for 10 min; 4°C hold. Inner nesting PCR products were visualized on 1.8% TBE gels stained with Sybrsafe DNA stain (Invitrogen Life Technologies and Molecular Probes).
| Results |
|---|
|
|
|---|
One of the elements necessary for a developmental model of TEC heterogeneity and TRA expression is ongoing differentiation and turnover of mTECs in adult thymus. We used two independent means to assess mitotic activity in mTECs, the first being a set of BrdU labeling experiments. Cohorts of age-matched adult BALB/c mice (46 wk old) received BrdU in their drinking water for a 5-wk pulse period and then returned to normal drinking water for a 5-wk chase period. At various intervals during this regimen, enzymatically dissociated thymi from test and control groups were enriched for epithelial cells by anti-CD45 bead depletion (18), and then assessed for CK and BrdU content by immunofluorescence microscopy (Fig. 1A). Triple staining experiments showed that almost all epithelial cells recovered were mTECs (as defined by staining results with Abs to Epcam (a medullary marker), BP-1 (a cortical marker), and CK; data not shown). Due to low recovery of cortical TECs using this method, we were unable to draw any conclusions regarding the mitotic activity of this population. The results of three independent experiments were consistent and showed significant accumulation of BrdU-labeled mTECs through the 5-wk labeling period, followed by a corresponding decrease in the number of cells retaining the BrdU label after cessation of BrdU treatment (Fig. 1B). The kinetics of the loss of the BrdU label following cessation of BrdU administration mirrored the kinetics of the labeling phase. We conclude that the proliferation documented here reflects epithelial turnover, because the adult thymus has already reached maximal size and is not undergoing prepubertal expansion (21).
|
12% of mTECs isolated from unmanipulated adult (610 wk) B6 thymus have synthesized DNA or recently divided. These results are consistent with the BrdU labeling data and clearly indicate that there is ongoing mitotic activity within the mTEC population of adult mice. It also appears that the rate of epithelial turnover is rather low, which would be consistent with the relative radiation resistance of this population that allows the generation of radiation bone marrow chimeras. Sorted mTECs express molecules known to regulate thymic organogenesis and TEC development
Although ongoing epithelial turnover implies ongoing mTEC differentiation, mitosis could be attributed to self-renewal of adult mTECs. We wanted to determine whether the known molecular regulators of thymic organogenesis and mTEC differentiation are expressed within this mitotically active population. We generated cDNA from 1 to 2 x 105 EpcamhighCD45BP-1 mTECs isolated from adult B6 thymus. We first evaluated these cDNAs for the expression of aire, Epcam, CD45, and a number of TRAs previously shown to be expressed by mTECs, including casein
(Csng), Purkinje cell protein 4 (Pcp4), and insulin II (Ins2) (3, 22). These PCR results showed that all of these molecules (except CD45) were expressed within the sorted mTEC population (Fig. 2A) and verified that we were examining a cell population comparable with those used in previous studies (3, 22).
|
receptor, have been shown to be expressed by postnatal thymic epithelium, the roles of Eya1, Hoxa3, and Pax9 in thymic epithelial development have been largely defined by the impact of their targeted inactivation on thymic organogenesis (reviewed in Ref. 23). Here we have determined that all of these known regulators of thymic epithelial development were expressed in bulk, sorted populations of adult mTECs (Fig. 2B). This expression pattern raised the possibility that Eya1, Hoxa3, and Pax9 may continue to play a role in the development and maturation of mTECs in the adult thymus. These results also demonstrate that the molecular regulators of thymic organogenesis and differentiation continue to be expressed within the adult mTEC population and suggest that mTEC differentiation from a resident progenitor population persists in adult thymus. Molecules characteristic of multipotential cells are expressed by mTECs
Based on the previous demonstration that epithelial foci in adult thymus ultrastructurally and phenotypically resemble other endodermally derived epithelium (reviewed in Ref.413), we have proposed the presence of multipotential progenitor epithelial cells in thymus (13, 14). Because the nature of and product-precursor relationships for mTEC subsets remain ill defined, we could not test this hypothesis directly. However, we could determine whether or not some of the transcription factors previously demonstrated to be characteristic of multipotential cell populations were expressed within the mTEC population. The molecules we examined are largely restricted to multipotential cell populations (such as embryonic stem cells and blastocysts) and are either considered to be characteristic markers of these populations or have been directly implicated in the maintenance of pluripotentiality in the cells that express them. These include Stella (developmental pluripotency associated 3; Dppa3) (24), Ehox (embryonic homeobox) (25), Foxd3 (also known as Genesis) (26), and undifferentiated transcription factor 1 (27). We determined that each of these molecules is expressed within the mTEC population (Fig. 3). We also evaluated the expression of Sox2, Nanog, and Oct4 (Pouf51) by mTECs because they have been reported to compose the core transcriptional regulatory circuit active in multipotential cells (28). As shown in Fig. 3, each of these molecules was also expressed within the mTEC population. In the case of Dppa3, FoxD3, and undifferentiated transcription factor 1, these results confirm array data that had identified these molecules as preferentially expressed by mTECs (29). Although it remains to be determined whether or not the expression of this group of transcription factors can be attributed to a restricted subset of mTECs with unique developmental potential, this work identifies new markers that may be useful for identifying mTEC progenitors and establishing product-precursor relationships among mTEC subsets. Given the reported function of these transcription factors in other cell populations referenced above, we hypothesize that their expression will likely be restricted to a discrete subset of immature mTECs.
|
Having observed expression of multiple regulators of thymic development and transcription factors characteristic of multipotential cells within the mitotically active mTEC population, we next wanted to determine whether mTECs also expressed transcription factors that regulate the differentiation of peripheral epithelial tissues where TRAs are characteristically expressed. We determined that in addition to FoxA1 and FoxA2, which are centrally involved in the development of a number of endodermal derivatives (30, 31, 32, 33, 34, 35, 36), mTECs expressed other transcription factors that are critical for the development of endodermal and neural tissues (Cdx1 (37, 38), Crx (39), Pdx1 (40, 41, 42, 43), Hnf-6 (44, 45, 46), and Foxg1 (47, 48); Fig. 4A). Multiple TRAs characteristic of the tissues that these transcription factors regulate, including blue cone opsin (retina), fatty-acid binding protein 1 (liver), gastric inhibitory polypeptide (Gip; intestinal epithelium), Ins2 (pancreas), and major urinary protein 1 (liver), were also expressed by mTECs (Fig. 2, Fig. 4B, and data not shown). These results show that the spectrum of TRAs expressed within the mTEC population extends beyond structural or terminal TRAs to include upstream developmental mediators characteristic of peripheral epithelial tissues.
|
Several studies have suggested that TRAs are broadly and promiscuously expressed by mature mTECs as a unique property of the mTEC lineage (22, 29); however, other studies have suggested that mTECs are highly heterogeneous and that the expression of individual TRAs is highly restricted to relatively rare mTECs (4, 5, 15). To distinguish between these two possibilities, we assessed the TRA expression profiles for small pools of sorted mTECs.
Using the method of Chiang and Melton (19), we generated cDNA from small pools (20 cells/pool) of sorted EpcamhighCD45BP-1 mTECs from the thymi of adult B6 mice. We analyzed the expression patterns of aire and 15 individual TRAs in pools of cells previously verified by PCR to be Epcam+ and MHC II+. As shown in Fig. 5, each individual TRA, with the notable exception of casein
(Csnb), was rarely expressed in the pools. The ubiquitous (and atypical) expression of Csnb suggests that this molecule is regulated differently than the other TRAs we assessed and may indicate that Csnb is expressed as part of a global mTEC expression profile. When all TRAs were considered, including Csnb, the average number of TRAs expressed per pool was closest to 2 of 15, with no pool expressing >4 of 15 TRAs. If Csnb is omitted from the analysis, the average dropped to 1 of 14 TRAs expressed per pool; no single pool expressed >3 of 14 TRAs, and 5 pools expressed 0 of 14 TRAs.
|
Expression of individual TRAs by functionally mature mTECs is also restricted to relatively rare cells
It has been proposed that mTECs progressively derepress TRA-containing loci as they mature, resulting in the highest levels of aire and TRA expression in the most mature mTECs (29), with high levels of MHC II protein expression used as the criterion for mature mTEC. We generated amplified cDNA from small pools (25 cells/pool) of sorted CD45Epcam+MHC IIhigh mTECs (Fig. 6A) from adult B6 mice and evaluated these for the expression of 17 TRAs. Although the frequency of TRA expression was higher in the MHC IIhigh pools than in the pools of unfractionated CD45Epcam+ mTECs analyzed in Fig. 5, the expression profiles for all of these TRAs (again, with the exception of Csnb) indicate that the frequency of expression for each individual TRA is also quite low in MHC IIhigh mTECs (Fig. 6B). A sample of the raw data represented in Fig. 6B is shown in Fig. 6C. All nine pools shown in Fig. 6 expressed aire and Csnb, whereas expression of thyroglobulin was observed only in pool 8. All pools except pool 3 were free of hematogenous cell contamination as determined by CD45 expression. When all TRAs were considered, including Csnb, the average number of TRAs expressed per pool was
3 of 17; no pool expressed >6 of 17 TRAs. If Csnb is omitted from the analysis, the average drops to 2 of 16 TRAs expressed per pool; no single pool expressed >5 of 16 TRAs, and 4 pools expressed 0 of 16 TRAs. The expression frequency for any individual TRA in the pools (excepting Csnb) ranged from 0 of 9 to 5 of 9. As above, we calculated that the percentage of MHC IIhigh mTECs expressing a given TRA ranges from <0.4 to 2% (0 of 225 cells to 5 of 225 cells). This value was also similar to previous estimates for the percentage of mTECs that express an individual TRA. These results show that even in a population of mature mTECs, the frequency of mTECs that expressed any individual TRA was low.
|
,
, and
(Csnb, Csng, and Csnk) that we observed within the pools are particularly noteworthy, given that these genes reside within the same locus on chromosome 5 (49). Expression of all three by mTECs has been put forth as a strong example of epigenetic derepression of TRAs (29), whereby the genes within the casein locus would be promiscuously opened in mature mTECs, resulting in the coordinated expression of these TRAs; the mechanisms proposed to account for this include run through transcription of the open locus (29). Our analysis of TRA expression within small pools of sorted mTEC clearly showed that Csng and Csnk appear to be expressed independently from both Csnb and each other; indeed, we did not find Csng and Csnk coexpressed within any pool (Figs. 5 and 6B). In contrast, Csnb was highly expressed in mTECs and may be expressed by every mTEC or a large, as yet undefined subset of mTECs. Further single-cell analysis will be required to determine whether that is the case.
Although we saw no evidence for coordinated expression of casein genes, we did observe that the pools that expressed multiple TRAs coexpressed Gip, pancreatic polypeptide (Ppy), and somatostatin (Sst). Although chromatin de-repression does not readily account for the coordinated expression of these genes (Gip and Ppy are located over 6 Mb apart on chromosome 11, Sst is located on chromosome 16), developmental mechanisms may explain the possible coordinated expression of these TRAs. The pancreatic transcription factor Pdx1 (pancreatic and duodenal homeobox1, or somatostatin transactivating factor-1) is known to directly regulate the transcription of both Gip (43) and Sst (40); the development of the cells that express Ppy in pancreas is also dependent on Pdx1 (41). It will be important to determine the extent to which other developmentally related genes are coordinately expressed in mTECs and the extent to which expression of developmentally related TRAs is coincident with the transcription factors that regulate them in peripheral tissues.
Nested PCR analysis of cDNA from pools and sorted single cells confirms the observed TRA expression patterns
We wanted to confirm the small pool data described above using a different method (20) and extend our observations to the single-cell level. To accomplish this, nested gene-specific PCR was performed with cDNA from both small pools (25 cells/pool; Fig. 7A) and individual (Fig. 7B) sorted MHC IIhighEpcam+CD45 cells. Although the use of nested PCR with cDNAs generated using gene-specific primers allowed for increased sensitivity, the spectrum of molecules that could be analyzed was significantly smaller. Amplified cDNA from all pools and single cells examined were first validated via secondary PCR for expression of MHC II. We then tested these pools and single cells for the presence of transcripts for aire and the aire-dependent TRA Gip (22). As shown in Fig. 7, Gip was expressed by only 1 of 23 MHC II+ pools (Fig. 7A and data not shown), and 0 of 40 MHC II+ individual mTECs (Fig. 7B and data not shown). These results were consistent with the pool data in Figs. 5 and 6 and confirm that expression of an individual TRA is limited to rare individual mTECs.
|
The expression of many of the TRAs assessed in Figs. 5 and 6 have been shown to be aire dependent in mTECs, including Csng, elastase 3B, Gip, iFABP, Ins2, major urinary protein (mup), Pcp4, and salivary protein 1 (Spt1; Ref. 22). Although almost all of these pools expressed aire (only pool 13 of the Epcamhigh mTECs did not; see Fig. 5), few of these pools expressed aire-dependent TRAs, and none of the pools expressed all, or even a majority of them. No individual aire-dependent TRA was expressed in >5 of the 21 pools that express aire. These results indicate that the expression of aire in an mTEC was not sufficient for the simultaneous coexpression of aire-dependent TRAs.
The data generated from the pools analyzed by gene-specific nesting (Fig. 7) was consistent with this finding. We found that MHC IIhighEpcam+CD45 mTECs express aire at a relatively low frequency (5 of 40 single cells; Fig. 4B and data not shown). As a population, these MHC IIhigh mTECs have been previously shown to express aire at relatively high levels (29); our finding that only 5 of 40 single cells and 15 of 23 pools tested express aire shows that only a minority of mTECs within the MHC IIhigh population actively express aire, a frequency consistent with immunohistochemical staining patterns for aire observed in adult thymus (Ref. 50 and A. G. Farr and G. O. Gillard, unpublished observations). As in our earlier experiments, expression of aire within an individual cell or pool did not correspond well with Gip expression. Gip, an aire-dependent TRA, was expressed within only 1 of the 15 aire+ pools, and by none of the 5 aire+ single cells. Whether the expression of Gip requires that aire be coexpressed within the same cell still remains to be determined, because the frequency of Gip-expressing cells was so low that examination of a significantly large cohort of sorted single cells will be required to address this question.
| Discussion |
|---|
|
|
|---|
First, generation of differentiated mTEC from a progenitor population predicts a level of epithelial turnover. We have shown that proliferation of TEC continues in the postnatal thymus, albeit at a relatively low level. A basal level of proliferation of mTECs indicates that this population turns over and likely has the capacity for self-renewal. On the basis of other differentiating systems, one could envision low rate of turnover of the least mature progenitor cells, followed by a high rate of turnover and expansion of a transit amplifying population, which then becomes postmitotic and terminally differentiates. The capacity of TECs to proliferate in response to exogenous keratinocyte growth factor in vivo and in vitro (51) suggest that this signaling pathway may play a critical role in regulating TE turnover in postnatal thymus.
The continuing differentiation of mTECs in the postnatal thymus presents alternative, but not mutually exclusive mechanisms for the generation of mTEC heterogeneity and the regulation of TRA expression. Differentiating epithelial cells could access multiple developmental pathways before specification to a terminal, perhaps thymic fate choice. The ability of immature cell populations to express transcription factors characteristic of alternate lineages before commitment to a specified terminal lineage has been observed in multiple settings (20, 52, 53, 54, 55) and may also be occurring in thymic epithelium. A multipotential thymic epithelial progenitor population could respond to developmental cues within the thymic microenvironment to adopt alternate epithelial cell fates. Thymic epithelium shares space with a complex mesenchymal compartment, together forming an environment replete with elements of several signaling pathways (bone morphogenetic protein (56, 57), Wnt (58), fibroblast growth factor (51, 59), Notch (60), and Hedgehog (61, 62)) that play critical roles in the specification and differentiation of a range of other epithelial tissues.
Invoking multipotential epithelial progenitor cells as a basis for mTEC heterogeneity and TRA expression in the adult thymus is partially based on reports that described an epithelial progenitor population in the fetal thymus (8, 9) and the assumption that these cells persist postnatally. mTEC progenitors that are multipotential and not lineage restricted would be predicted to share the molecular signature of other multipotential populations. We have shown that a number of transcription factors normally associated with multipotential cells, including Oct4, nanog, and Sox2, are expressed within the mTEC population. Defining the expression pattern of these transcription factors within the mTEC population may prove useful in parsing mTEC heterogeneity and facilitate the identification of the progenitor epithelial cell population in the adult thymus.
That TECs can adopt alternative epithelial fates is supported by several lines of experimental evidence. In the absence of functional foxn1, the thymic rudiment of nude mice displays the morphology and organization of early fetal lung and a respiratory phenotype (63). One interpretation of this finding is that the progenitor population in the nude thymus is multipotential and that foxn1 is required to specify a thymic epithelial fate. There are also indications that mTEC in normal mice adopt alternative developmental fates. Cystic epithelial structures that are lined, in part, by ciliated epithelial cells and display a wide range of respiratory system-related genes, including plunc, surfactants, and CC10, are commonly found in normal thymus (16). Foci of epithelial cells resembling small thyroid follicles within the normal thymus have also been described (13). Although these latter instances may represent extreme examples of the developmental potential of thymic epithelial progenitor cells where terminal differentiation has occurred, they highlight the potential plasticity of epithelial cells in the thymus.
Our results demonstrating the expression of peripheral developmental transcription factors by mTECs is significant in two respects. First, these findings sustantially expand the number of candidate transcription factors that may be involved in mTEC differentiation, thus providing additional tools for analyzing thymic epithelial heterogeneity. Second, accepting the premise that mTEC heterogeneity arises as a consequence of differentiation and that TRA expression is a reflection of this process, it is reasonable to propose that conventional, lineage-restricted developmental transcriptional networks could regulate expression of downstream TRAs in mTECs. We have previously shown that Ttf1 (Nkx2.1), which is centrally involved in lung and thyroid development, is expressed by mTEC (16), as are multiple lung- and thyroid-specific TRAs (3, 16, 22, 29). We have shown here that mTEC also express a number of transcription factors that play key roles in the differentiation of other peripheral tissues; the mTEC population also expresses a number of TRAs characteristic of these peripheral tissues. Because some of these transcription factors are required for efficient expression of TRAs in peripheral tissues, their functional inactivation might result in lacunar, lineage-restricted deficiencies in TRA expression by mTEC. The observation that multiple Pdx1-dependent and lineage-related TRAs are coexpressed within small pools of sorted mTEC supports this notion that tissue-specific transcriptional networks may be functioning in mTEC as in their native peripheral epithelial tissues.
We favor the interpretation that expression of these tissue-specific transcription factors (and their downstream targets) results from conserved developmental mechanisms during mTEC differentiation; however, an alternative interpretation is that the derepression of chromatin in mature mTECs leads to the promiscuous expression of multiple developmental regulators (and their downstream targets) by mature mTECs. Given the gaps in our knowledge regarding the extent of mTEC heterogeneity, and the developmental relationships among mTEC populations, we cannot presently distinguish between these two interpretations.
Our examination of individual and small pools of sorted mTECs revealed that expression of a given TRA is limited to rare mTECs and that individual mTECs do not express a broad range of TRAs. These data strongly support the notion that mTECs represent a diverse population and that TRA expression reflects that diversity. Previous reports have examined TRA expression in pools of 104106 mTECs (3, 22, 29); by using bulk populations, the diversity of the mTEC population is not apparent because the results represent averaged expression from many cells. If the TRA expression results for the pools shown in Fig. 6 are summed, almost all of the TRAs examined were expressed (12 of 17 TRA expressed by 225 MHC IIhigh mTEC). On the basis of such a summation, one might conclude that the genes within the casein locus are coordinately expressed; examination of smaller samples reveals that these genes are being independently regulated within individual mTECs.
These small pool/single-cell analyses also provide insight into the relationship between the expression of aire and aire-dependent TRAs. The poor correlation between aire and aire-dependent TRA expression indicates that aire expression is not sufficient for the simultaneous coexpression of TRAs in mTEC. This may reflect a requirement for other cofactors or indicate a temporal separation between aire and TRA expression. It is possible that the importance of aire in regulating some TRAs reflects a more conventional regulatory transcription factor activity that influences the ability of immature mTEC to initiate or sustain developmental hierarchies. Alternatively, aire might regulate the temporal pattern of mTEC turnover, differentiation, or commitment, which could, in turn, impact the duration of TRA expression or the temporal window during which immature mTECs could follow alternate lineage pathways. In each of these scenarios, lack of aire would result in lacunar defects in mTEC differentiation and a reduced spectrum of TRA expression.
Our data establish that the mTEC population is dynamic and diverse and that it expresses a wide range of transcription factors known to regulate the developmental potential and lineage specification of multiple cell types. These results provide the basis for a novel framework for understanding and investigating the related processes of thymic epithelial cell potential, differentiation, and function.
| Acknowledgments |
|---|
| Disclosures |
|---|
|
|
|---|
| Footnotes |
|---|
1 This work was supported by the National Institutes of Health Grants AI 24137 and AI 059575. G.O.G. was supported in part by training grants from the National Institutes of Health and the Cancer Research Institute. ![]()
2 Address correspondence and reprint requests to Dr. Andrew G. Farr, Department of Biological Structure, Box 357420, University of Washington, Seattle, WA 98195-7420. E-mail address: farr{at}u.washington.edu ![]()
3 Abbreviations used in this paper: TE, thymic epithelium; aire, autoimmune regulator; CK, cytokeratin; Csnb, casein
; Csng, casein
; Csnk, casein
; Dppa3, developmental pluripotency associated 3; Epcam, epithelial cell adhesion molecule; FoxN1, forkhead N1; Gip, gastric inhibitory protein; Ins2, insulin II; mTEC, medullary thymic epithelial cell; Pdx1, pancreatic and duodenal homeobox 1; Ppy, pancreatic polypeptide; Pcp4, Purkinje cell protein 4; Sst, somatostatin; Spt1, salivary protein 1; Spt2, salivary protein 2; TRA, tissue-restricted Ag; CMF, Ca2+ and Mg2+ free; TBE, Tris-buffered EDTA; mup, major urinary protein; RT, reverse transcriptase. ![]()
Received for publication January 20, 2006. Accepted for publication March 1, 2006.
| References |
|---|
|
|
|---|
is required for notochord development in the mouse embryo. Cell 78: 575-588. [Medline]
cells of pancreas, duodenal epithelium and pancreatic exocrine and endocrine progenitors during ontogeny. Development 121: 11-18. [Abstract]
in the developing murine liver and pancreas. Dev. Biol. 192: 228-246. [Medline]