Abstract
The IL-2 receptor γ common (IL-2Rγc) chain is the shared subunit of the receptors for the IL-2 family of cytokines, which mediate signaling through JAK3 and various downstream pathways to regulate lymphopoiesis. Inactivating mutations in human IL-2Rγc result in SCID, a primary immunodeficiency characterized by greatly reduced numbers of lymphocytes. This study used bioinformatics, expression analysis, gene ablation, and specific pharmacologic inhibitors to investigate the function of two putative zebrafish IL-2Rγc paralogs, il-2rγc.a and il-2rγc.b, and downstream signaling components during early lymphopoiesis. Expression of il-2rγc.a commenced at 16 h post fertilization (hpf) and rose steadily from 4–6 d postfertilization (dpf) in the developing thymus, with il-2rγc.a expression also confirmed in adult T and B lymphocytes. Transcripts of il-2rγc.b were first observed from 8 hpf, but waned from 16 hpf before reaching maximal expression at 6 dpf, but this was not evident in the thymus. Knockdown of il-2rγc.a, but not il-2rγc.b, substantially reduced embryonic lymphopoiesis without affecting other aspects of hematopoiesis. Specific targeting of zebrafish Jak3 exerted a similar effect on lymphopoiesis, whereas ablation of zebrafish Stat5.1 and pharmacologic inhibition of PI3K and MEK also produced significant but smaller effects. Ablation of il-2rγc.a was further demonstrated to lead to an absence of mature T cells, but not B cells in juvenile fish. These results indicate that conserved IL-2Rγc signaling via JAK3 plays a key role during early zebrafish lymphopoiesis, which can be potentially targeted to generate a zebrafish model of human SCID.
Introduction
Interleukin-2 receptor γ common (IL-2Rγc) chain represents the shared component of the receptors for the IL-2 family of cytokines IL-2, IL-4, IL-7, IL-9, IL-15, and IL-21, each of which regulates different aspects of immune development and function (1–4). Specifically, IL-2 stimulates the growth, differentiation, and activation of various T and NK cell populations (5); IL-4 is involved in B cell proliferation, Ig class switching, and Th2 cell development (6); and IL-7 contributes to the development, survival, and homeostatic proliferation of T cells, especially memory T cells (7). IL-9 exerts a wider range of effects, such as mediating the growth and functional activation of T cells and mast cells, and supporting the differentiation of hematopoietic stem cells (8), whereas IL-15 regulates the proliferation, differentiation, and survival of many cell types, including B cells, NK cells, mast cells, neutrophils, eosinophils, monocytes/macrophages, and dendritic cells (9, 10). Finally, IL-21 stimulates lymphoid cell proliferation and the differentiation of B cells to plasma cells, and it regulates apoptosis in B and NK cells (4).
The IL-2Rγc chain acts as the major signal transduction component of the IL-2 receptor family. Ligand binding stimulates the activation of the tyrosine kinase JAK3 associated with the intracellular region of IL-2Rγc, which activates a number of downstream intracellular pathways, including those involving STAT5, PI3K, and ERK, the latter lying downstream of MEK (11, 12). Inactivating mutations in human IL-2Rγc cause T−B+ SCID, characterized by decreased numbers of T cells and a diminished immune response (13). Mice lacking IL-2rγc display a T−B− form of SCID (14), with reduced B and T cells (15), and have proved to be an invaluable model for a range of studies, especially relating to immunity and cancer (16–21).
The zebrafish is now established as an important alternate model for the study of vertebrate development and disease, with particular relevance to hematopoiesis and immunity (22, 23). Like mammals, zebrafish undergo distinct phases of hematopoiesis, with lymphoid progenitors seeding the thymus early in development following the establishment of definitive hematopoiesis (24). Zebrafish also shows broad conservation of cytokine receptors and downstream signaling pathway components (25, 26), particularly those known to facilitate the development of blood and immune cells (27–31).
We have previously identified two putative paralogs of the IL-2Rγc gene in zebrafish, il-2rγc.a and il-2rγc.b (25), but their roles have not been elucidated. In this study, we show that il-2rγc.a, but not il-2rγc.b, is involved in early zebrafish lymphopoiesis. The il-2rγc.a gene is expressed in the developing thymus and adult lymphocytes, with ablation of il-2rγc.a leading to a significant and specific reduction in T lymphopoiesis. This phenotype could be mimicked by targeting zebrafish Jak3 by either morpholino-mediated knockdown or pharmacological inhibition, indicating a conserved role for IL-2Rγc/JAK3 signaling across vertebrate lymphopoiesis, with contributions from the STAT5, PI3K, and MEK/ERK pathways also being identified.
Materials and Methods
Nomenclature conventions
Nomenclature rules for zebrafish, fugu, chicken, mouse, and human genes and proteins differ. Gene and protein names are presented according to the respective nomenclature conventions (zebrafish and fugu: il-2rγc.a, Il-2rγc.a; human and chicken: IL-2Rγc, IL-2Rγc; and mouse: Il-2rγc, IL-2Rγc).
Analysis of il-2rγc paralogs from zebrafish
Expressed sequence tags corresponding to the two putative Danio rerio (zebrafish) IL-2Rγc paralogs, il-2rγc.a and il-2rγc.b (25), were identified using BLASTX (32), with additional sequences obtained from RT-PCR and 5′RACE products. These sequences were assembled using Sequencher (Gene Codes Corporation) and corresponded to full-length reference sequences deposited at GenBank (http://www.ncbi.nlm.nih.gov/genbank/; accession number NP_001121743.1 il-2rγc.a, NP_001116522.1 il-2rγc.b). Multiple-sequence alignment of the encoded protein sequences was performed using the CLUSTALX program (33), along with zebrafish thymic stromal lymphopoietin receptor (Tslpr) (CAM88660.1), Il-13rα1 (CAI94933.1) and Il-13rα2 (NP_001107203.1), Homo sapiens (human) IL-2Rγc (AAA59145.1), thymic stromal lymphopoietin receptor (TSLPR) (NP_071431.2), IL-13Rα1 (EAW89893.1) and IL-13Rα2 (AAH20739.1), Mus musculus (mouse) IL-2Rγc (AAH14720.1), TSLPR (NP_001158207.1), IL-13Rα1 (AAH59939.1) and IL-13Rα2 (EDL14709.1), Gallus gallus (chicken) IL-2Rγc (NP_989858.1), TSLPR (XP_416864.3), IL-13Rα1 (XP_420218.3) and IL-13Rα2 (NP_001041543.1), and Takifugu rubripes (fugu) Il-2rγc (NP_001129354.1) and Il-13rα2 (XP_003971306.1). A phylogenic tree was derived from this alignment using the Neighbor-Joining algorithm and visualized with NJplot (34). Genomic sequences corresponding to the zebrafish il-2rγc.a and il-2rγc.b genes were identified using BLASTN, and the positions of intron/exon boundaries were determined by alignment with the corresponding zebrafish il-2rγc cDNA sequence, applying the ‘GT-AG’ rule (35).
Zebrafish husbandry and manipulation
Wild-type as well as Tg(lck:lck-EGFP)cz1 (Lck:eGFP) (36) and Tg(Cau.ighv-ighm:EGFP)sd19 (IgM1:eGFP) (37) transgenic zebrafish were maintained using standard husbandry practices (38). Wild-type embryos at the 1–8-cell stage were injected with anti-sense morpholinos diluted in 1× Danieau buffer (58 mM NaCl, 0.7 mM KCl, 0.4 mM MgSO4, 0.6 mM CaCl2, 5.0 mM HEPES; pH 9). Morpholinos used were: il-2rγc.aSS1 (5′-CTTGTTTCATTAATACATACCCGCC, 3 mM), il-2rγc.aSS2 (5′-CTTCACTGTTTAACATCACAGTAAC, 7 mM), il-2rγc.bSS (5′- GCAGCACTAAACAAAATATGATGCA, 5 mM), jak3SS (5′-TTAAATGTGTTAGTGTCTCACCACT, 2 mM), stat5.1SS (5′-GTGAACTTGTGACTTACCAGAGTTG, 1 mM), stat5.2SS (5′ GTTGTCATCTGGTGCTCATACCTTC, 1 mM) and standard control morpholino (5′-CCTCTTACCTCAGTTACAATTTATA, 1-7 mM as appropriate). In some experiments, embryos were coinjected with 100 pg/nl mRNA encoding eGFP or constitutively active version of zebrafish Jak3 A572V (this study) or Stat5.1 (H298R/N714F) (39) and standard control and il-2rγc.aSS1 morpholinos. Alternatively, 1-cell embryos were injected with 100 pg/nl mRNA encoding transcription activator-like effector nucleases (TALENs) (40) targeting exon 3 of il-2rγc.a, raised to adulthood intercrossed and progeny screened by RFLP with NdeI to identify mutants. Embryos were also treated with the JAK3 inhibitor tofacitinib (41) at 30 and 60 μM, MEK inhibitor PD98059 (42, 43) at 25 μM, and PI3K inhibitor LY294002 (44) at 15 μM from 56 h postfertilization (hpf) and then fixed at the appropriate time points. Single-cell suspensions were prepared from the thymus of Lck:eGFP and kidney of IgM1:eGFP transgenic zebrafish, as described (45). Sytox Red Dead Cell Stain was used for live cell discrimination (Molecular Probes), and cell sorting was performed with a FACSAria II (BD Biosciences). National and institutional guidelines for the care and use of laboratory animals were followed in all studies.
RT-PCR, quantitative RT-PCR, and 5′RACE
Total RNA was extracted from 30 zebrafish embryos using Trizol reagent (Life Technologies) following the manufacturer’s recommendations and then resuspended in nuclease-free water. This RNA was subjected to semiquantitative RT-PCR with the following primers: il-2rγc.a 5′-CAGGCGTCAGGACCACATACAG and 5′- CTCTCACTATCACTGCTGGACTGG (time course/sequencing), 5′-AGAAGTGCGTTATGTGACCCTG and 5′- TCTGGTCAGTCCTGTAACGAAC (SS1 morpholino titration), 5′- CGAAGACTGTCCTGAATATGAGAC and 5′- AGACTCACTCCACTCGCTCCAG (SS2 morpholino titration and TALEN sequencing), 5′- AGAAGTGCGTTATGTGACCCTG and 5′- TCTGGTCAGTCCTGTAACGAAC (sequencing), 5′- TATGCTGAAAGAATATGTGAAG and 5′- AGACTCACTCCACTCGCTCCAG (sequencing), 5′- CGTCACTGGTCTTGTATGCTG and 5′- GTCGTTTTTCCTCATCAATCTGC (sequencing), il-2rγc.b 5′-TGGAACGAGCACAGCGACAC and 5′- GAAGAACCGCAGGAATCAGC (time course/sequencing), 5′- CAGTCATTTGTCACTCAGACGCTC and 5′- GATGCAGGTTTTACGGAGAGGT (morpholino titration), 5′- CAGTCATTTGTCACTCAGACGCTC and 5′- GAAGAACCGCAGGAATCAGC (sequencing), 5′- CGTCATACAGTGTGTCTCCAGTCTC and 5′- ACAGTATGGATGAGATGAGGATGG (sequencing), jak3 5′- AACTCAGAGACCACCTTCAGCA and 5′- GTGTGACCACCCTTCCTTCC (morpholino titration), stat5.1 5′- CAGGGAGATGCTCTACACCAG and 5′- CTCCGACTTGATGCTCTGC (morpholino titration), stat5.2 5′- CAGCACTTCCCCATTGAGG and 5′- CTCGTGTCAGCCAGGTCTC (morpholino titration), and β-actin 5′-TGGCATCACACCTTCTAC and 5′-AGACCATCACCAGAGTCC. RNA was also subjected to 5′RACE to generate cDNA using a 5′ RACE kit (Life Technologies), according to the manufacturer’s protocol. Gene-specific primers were il-2rγc.a 5′-AGGCTTTTTCAGTTCC (gsp1), 5′- AGACTCACTCCACTCGCTCCAG (gsp2), il-2rγc.a 5′- AGACTCACTCCACTCGCTCCAG (gsp3), 5′GTTGTCGTTCTTCGTAACATTC (gsp4), 5′- AACCTTTCGCTGTGG (new gsp1), il-2rγc.b 5′- AAAGACTGGCTTGGGT (gsp1), il-2rγc.b 5′- GATGCAGGTTTTACGGAGAGGT (gsp2) ,and il-2rγc.b 5′- GAAAGTGTGTCGCTGTGCTC (gsp3).
Total RNA was also extracted from zebrafish embryos, larvae, and pooled adult zebrafish tissues with RNeasy Mini Kit (Qiagen) according to the manufacturer’s protocol for animal tissues. This was subjected to semiquantitative RT-PCR with previously published primers for TCR-β chains (Vb1.5/17.5, Vb12) and Ig heavy chains (igVH1, igVH4) (46, 47). Total RNA was also subjected to quantitative RT-PCR (QRT-PCR) with the following primers: il-2rγc.a 5′-GTCACTGGTCTTGTATGCTGT and 5′-GCTCTCACTATCACTGCTGG, il-2rγc.b 5′-AGAAAGACCCAAGCCAGT and 5′-ATCTTTTTCCTCTCACAGTACC, jak3 5′-AACAGAGCGAGCAGCAGAGAG and 5′-GTGTGACCACCCTTCCTTCC, pik3cγ 5′-AGGGGCACTTGTGATTGAG and 5′-CTTCACTATTTCAATCATTCCAA, β-actin 5′-TGGCATCACACCTTCTAC and 5′-AGACCATCACCAGAGTCC and map2k1 and map2k2 (48). Data were normalized to β-actin, and fold change was calculated using the ΔΔCt method.
Total RNA was extracted from sorted cells using Trizol reagent and DirectZol RNA MiniPrep spin columns (ZymoResearch). cDNA was prepared using QuantiTect Reverse Transcription Kit (Qiagen) and subjected to QRT-PCR with il-2rγc.a, il-2rγc.b, lck, ighm, and ef1a primers (37). Each primer set flanked splice sites such that amplification of contaminating genomic DNA would produce considerably larger fragments in each case. Data were normalized to ef1a, and fold change was calculated using the ΔΔCt method. Each primer set flanked splice sites such that amplification of contaminating genomic DNA would produce considerably larger fragments in each case.
Genomic DNA analysis
Genomic DNA was obtained from pooled F1 embryos with QuickExtract following the manufacturer’s instructions. This DNA was subjected to PCR with specific primers for il-2rγc.a and analyzed by RFLP with NdeI and by Sanger sequencing at the Australian Genome Research Facility.
Whole-mount in situ hybridization (WISH) and histochemistry
Embryos were dechorionated and fixed for 1-2 d in 4% (w/v) paraformaldehyde at 4°C prior to whole-mount in situ hybridization (WISH) with DIG-labeled anti-sense probes, as described (49, 50). Freshly anesthetized embryos were subjected to O-dianisidine staining of hemoglobin, as described (51). Quantitation was achieved by measuring the area of staining relative to eye diameter as determined using CellSens Dimension 1.6 software (Olympus) or counting of individual cells on ∼30 embryos. Data were analyzed for significance with a Student t test, with Welch’s correction used where necessary to account for variation between the data groups.
Transient expression and analysis in human 293T cells
Human 293T cells were grown to 50–80% confluency before transfection with pBKCMV or pBKCMV expressing zebrafish Jak3 and/or Stat5.1 using lipofectamine reagent (Life Technologies). After incubation at 37°C in 10% (v/v) CO2 for 2 d, a total cell lysate was prepared and subjected to Western blot analysis with anti-phospho Stat5.1 (Millipore 05-495) and anti-GAPDH (Millipore CB1001).
Results
Two IL-2Rγc paralogs in zebrafish
We have previously identified two putative zebrafish paralogs of the mammalian IL-2Rγc gene, called il-2rγc.a and il-2rγc.b (25). The full-length sequences of these genes were assembled from expressed sequence tags and predicted mRNA sequences present in gene databases, supplemented by sequencing of RT-PCR and 5′RACE products, and the encoded zebrafish protein sequences deduced. A phylogenetic tree was generated using the human, mouse, chicken, fugu, and zebrafish IL-2Rγc sequences along with those of the closely related thymic stromal lymphopoietin receptor (TSLPR), with the more divergent IL-13Rα sequences used as an out-group (Fig. 1A). This analysis grouped the zebrafish Il-2rγc.a and Il-2rγc.b sequences in a clade with the other IL-2Rγc chains supported by strong bootstrapping values, which was distinct from the related TSLPR clade. Of the two zebrafish sequences, Il-2rγc.b was the more divergent.
Characterization of zebrafish il-2rγc genes and their targeting with morpholinos. (A) Phylogenetic analysis of IL-2Rγc proteins. A phylogenetic tree was constructed from a ClustalX multiple alignment of zebrafish (Dr) Il-2rγc.a, Il-2rγc.b, Il-13rα1, Il-13rα2, and Tslpr along with the human (Hs), mouse (Mm), chicken (Gg), and fugu (Tr) IL-2Rγc, IL-13Rα1, IL-13Rα2, and TSLPR using the Neighbor-Joining algorithm with 1000 replicates. The IL-2Rγc clade is highlighted with a red box and bootstrap values >90% are bolded. (B) Conserved elements within IL-2Rγc intracellular domains. The intracellular regions of human, mouse, chicken, fugu, and zebrafish IL-2Rγc sequences were aligned using ClustalX. The sequences around specific motifs are shown with similar residues colored the same, and conservation across the six sequences indicated (identical *, highly similar:, similar .). (C and D) Expression of il-2rγc genes during zebrafish embryogenesis. Embryos were harvested during embryogenesis at the times indicated and analyzed by RT-PCR (left) and QRT-PCR (right) with primers specific for il-2rγc.a (C) or il-2rγc.b (D) along with β-actin controls with the sizes of amplified products indicated. Representative RT-negative controls are shown (7 dpf-) to confirm the absence of genomic contamination in RT-PCR, with QRT-PCR displayed as a percentage of maximal expression. (E–H) WISH expression of il-2rγc.a and il-2rγc.b. Embryos at 5 dpf were subjected to WISH with il-2rγc.a and il-2rγc.b sense (S) and antisense (AS) RNA probes. Staining in the thymus with il-2rγc.a is indicated with an arrow. (I and J) Adult expression of il-2rγc.a and il-2rγc.b. Adult zebrafish tissues were analyzed with QRT-PCR with primers specific for il2rγc.a (I) or il2rγc.b (J), with levels displayed as a percentage of maximal expression. (K and L) Expression of il-2rγc.a and il-2rγc.b in adult lymphoid cells. GFP+ lymphocytes were sorted from the thymus of Lck:eGFP (K) and kidney of IgM1:eGFP (L) transgenic fish. Shown are the normalized relative transcript levels of il-2rγc.a (black bars), il-2rγc.b (white bars), and lck or ighm (checked bars), relative to unsorted cells, as determined with QRT-PCR (mean ± SEM of three technical replicates). *p < 0.05, **p < 0.001.
Characterization of zebrafish il-2rγc genes and their targeting with morpholinos. (A) Phylogenetic analysis of IL-2Rγc proteins. A phylogenetic tree was constructed from a ClustalX multiple alignment of zebrafish (Dr) Il-2rγc.a, Il-2rγc.b, Il-13rα1, Il-13rα2, and Tslpr along with the human (Hs), mouse (Mm), chicken (Gg), and fugu (Tr) IL-2Rγc, IL-13Rα1, IL-13Rα2, and TSLPR using the Neighbor-Joining algorithm with 1000 replicates. The IL-2Rγc clade is highlighted with a red box and bootstrap values >90% are bolded. (B) Conserved elements within IL-2Rγc intracellular domains. The intracellular regions of human, mouse, chicken, fugu, and zebrafish IL-2Rγc sequences were aligned using ClustalX. The sequences around specific motifs are shown with similar residues colored the same, and conservation across the six sequences indicated (identical *, highly similar:, similar .). (C and D) Expression of il-2rγc genes during zebrafish embryogenesis. Embryos were harvested during embryogenesis at the times indicated and analyzed by RT-PCR (left) and QRT-PCR (right) with primers specific for il-2rγc.a (C) or il-2rγc.b (D) along with β-actin controls with the sizes of amplified products indicated. Representative RT-negative controls are shown (7 dpf-) to confirm the absence of genomic contamination in RT-PCR, with QRT-PCR displayed as a percentage of maximal expression. (E–H) WISH expression of il-2rγc.a and il-2rγc.b. Embryos at 5 dpf were subjected to WISH with il-2rγc.a and il-2rγc.b sense (S) and antisense (AS) RNA probes. Staining in the thymus with il-2rγc.a is indicated with an arrow. (I and J) Adult expression of il-2rγc.a and il-2rγc.b. Adult zebrafish tissues were analyzed with QRT-PCR with primers specific for il2rγc.a (I) or il2rγc.b (J), with levels displayed as a percentage of maximal expression. (K and L) Expression of il-2rγc.a and il-2rγc.b in adult lymphoid cells. GFP+ lymphocytes were sorted from the thymus of Lck:eGFP (K) and kidney of IgM1:eGFP (L) transgenic fish. Shown are the normalized relative transcript levels of il-2rγc.a (black bars), il-2rγc.b (white bars), and lck or ighm (checked bars), relative to unsorted cells, as determined with QRT-PCR (mean ± SEM of three technical replicates). *p < 0.05, **p < 0.001.
Alignment of the human, mouse, chicken, fugu, and zebrafish IL-2Rγc chains confirmed the conservation of key domains in both zebrafish sequences including an extracellular cytokine homology domain, a transmembrane domain and an intracellular domain (data not shown), including a Box 1 motif (Fig. 1B). Importantly, Il-2rγc.a showed conservation of three of the four intracellular tyrosines present in mammalian counterparts, with one additional nonconserved tyrosine, whereas Il-2rγc.b possessed only one conserved tyrosine and showed considerable divergence in its extended cytoplasmic region.
Expression of zebrafish IL-2Rγc paralogs
RT-PCR and QRT-PCR were used to characterize the expression of il-2rγc.a and il-2rγc.b during zebrafish embryogenesis. Expression of il-2rγc.a was first apparent at 16 hpf and then increased particularly from 4 days postfertilization (dpf) during the establishment of lymphopoiesis (Fig. 1C) (24). Conversely, il-2rγc.b showed biphasic expression with strong expression at 8–24 hpf, which waned before increasing from 5 dpf to maximal expression at 6 dpf (Fig. 1D). WISH was performed using probes for il-2rγc.a and il-2rγc.b. Distinct expression within the thymus was only evident for il-2rγc.a (Fig. 1F), which was not observed for il-2rγc.b (Fig. 1H) or sense controls (Fig. 1E, 1G). In adult zebrafish, il-2rγc.a was broadly expressed, but with highest expression in the thymus, kidney, and spleen (Fig. 1I), which represent key lymphoid organs (45, 52), whereas il-2rγc.b expression was highest in the spleen (Fig. 1J). Further analysis in transgenic zebrafish lines revealed that expression of both il-2rγc.a and il-2rγc.b was higher in the lck+ T cell population within the thymus (Fig. 1K) and the igm+ B cell population within the kidney (Fig. 1L).
Knockdown of il-2rγc.a but not il-2rγc.b affects lymphopoiesis
To investigate the functions of zebrafish il-2rγc.a and il-2rγc.b, a morpholino-mediated knockdown strategy was used (53), with morpholinos designed to interfere with splicing of the premRNA in the region encoding the cytokine homology domain (Supplemental Fig. 1A–C). The level of gene knockdown was determined by RT-PCR, which confirmed robust knockdown in embryos injected with 3 mM il-2rγc.aSS1, 7 mM il-2rγc.aSS2, and 5 mM il-2rγc.bSS morpholino, compared with those injected with standard control morpholino at the equivalent concentration (Supplemental Fig. 1D, 1E). Sequence analysis of alternate transcripts seen with il-2rγc.aSS1 and il-2rγc.aSS2 confirmed that these would encode severely truncated proteins in each case (Supplemental Fig. 1F).
Morpholino-injected embryos showed no overt phenotypes; however, analysis by WISH with an early lymphocytic marker, ikaros (54), revealed significantly decreased expression from 3 dpf in the thymus of embryos injected with il-2rγc.aSS1 (Fig. 2B, 2E) and il-2rγc.aSS2 (Fig. 2C, 2F) compared with those injected with the standard control morpholino (Fig. 2A), which continued at 5 dpf (Fig. 2I, 2K, and data not shown). In contrast, ikaros expression in the thymus was not significantly altered in il-2rγc.bSS-injected embryos at any time point tested (Fig. 2D, 2G, 2J, 2L).
Effect of il-2rγc ablation on lymphoid cells. (A–F’) WISH analysis of morphants embryos. Embryos injected with standard control (Co Mo), il-2rγc.aSS1 (γc.a Mo SS1), il-2rγc.aSS2 (γc.a Mo SS2), or il-2rγc.bSS (γc.b Mo SS) morpholinos were subjected to WISH with ikaros (A–L), rag1 (M–X), tcrα (Y–B’), and lck (C’–F’) at the times indicated. The extent of ikaros and rag1 staining in the thymi (shown with the arrows) was determined as a ratio to eye size, averaged for individual embryos, with mean and SEM shown in red for il-2rγc.aSS1 (E, K, Q, and W), il-2rγc.aSS2 (F and R), and il-2rγc.bSS (G, L, S, and X). (G’ and H’) RT-PCR analysis of mutant larvae. Expression of TCR β-chain rearrangements V(D)J-Cβ Vb1.5 and Vb12, Ig H chain gene rearrangements VH1 and VH4 and β-actin as a control were examined in wild-type (++) and il-2rγc.a mutant (----) zebrafish by RT-PCR for expression of TCR and Ig (G’). RT-negative controls yielded no products (data not shown). Genomic DNA was extracted from these larvae and was subjected to PCR with primers surrounding the mutation site and digested with NdeI for genotyping (H’). The arrowheads indicate the cleaved products produced from wild-type alleles, and the arrow shows uncleaved products in embryos carrying mutant alleles. Statistically significant differences are indicated. *p < 0.05, ** p < 0.001. n/s, not significant.
Effect of il-2rγc ablation on lymphoid cells. (A–F’) WISH analysis of morphants embryos. Embryos injected with standard control (Co Mo), il-2rγc.aSS1 (γc.a Mo SS1), il-2rγc.aSS2 (γc.a Mo SS2), or il-2rγc.bSS (γc.b Mo SS) morpholinos were subjected to WISH with ikaros (A–L), rag1 (M–X), tcrα (Y–B’), and lck (C’–F’) at the times indicated. The extent of ikaros and rag1 staining in the thymi (shown with the arrows) was determined as a ratio to eye size, averaged for individual embryos, with mean and SEM shown in red for il-2rγc.aSS1 (E, K, Q, and W), il-2rγc.aSS2 (F and R), and il-2rγc.bSS (G, L, S, and X). (G’ and H’) RT-PCR analysis of mutant larvae. Expression of TCR β-chain rearrangements V(D)J-Cβ Vb1.5 and Vb12, Ig H chain gene rearrangements VH1 and VH4 and β-actin as a control were examined in wild-type (++) and il-2rγc.a mutant (----) zebrafish by RT-PCR for expression of TCR and Ig (G’). RT-negative controls yielded no products (data not shown). Genomic DNA was extracted from these larvae and was subjected to PCR with primers surrounding the mutation site and digested with NdeI for genotyping (H’). The arrowheads indicate the cleaved products produced from wild-type alleles, and the arrow shows uncleaved products in embryos carrying mutant alleles. Statistically significant differences are indicated. *p < 0.05, ** p < 0.001. n/s, not significant.
To confirm these effects, embryos were analyzed with a range of other lymphocyte markers. Significantly decreased expression of rag1, a marker of more mature lymphoid cells (55), was also observed 4–6 dpf in embryos injected with il-2rγc.aSS1 (Fig. 2N, 2Q, 2U, 2W, and data not shown) and il-2rγc.aSS2 (Fig. 2O, 2R, and data not shown), but not those injected with il-2rγc.bSS (Fig. 2P, 2S, 2V, 2X). This effect on lymphopoiesis was confirmed with the decreased expression of other lymphocyte-specific markers, tcrα and lck (36, 56), at 5 dpf in embryos injected with il-2rγc.aSS1 (Fig. 2Z, 2D’) or il-2rγc.aSS2 (Fig. 2A’, 2E’) relative to controls (Fig. 2Y, 2C’), which was not observed with il-2rγc.bSS (Fig. 2B’, 2F’).
To analyze the effects of il-2rγc.a ablation on mature T and B cells, TALEN-mediated genome targeting was used because the effects of morpholinos are only transient (57, 58). Specific TALENs directed to exon 3 (Supplemental Fig. 2A) were injected into one-cell embryos, which were raised to adulthood and intercrossed. One pair produced ∼25% progeny that recapitulated the loss of rag1 at 5 dpf (Supplemental Fig. 2E–H compared with Supplemental Fig. 2B–D), which were shown to have both il-2rγc.a alleles mutated by RFLP (Supplemental Fig. 2I) and sequence analysis (Supplemental Fig. 2J). Siblings from the cross were raised to 28 dpf and analyzed for the expression of rearranged TCR-β and Ig genes, as markers of mature T and B cells, respectively (46) (Fig. 2G’), and also genotyped (Fig. 2H’). This demonstrated the presence of mature T and B cells in wild-type larvae, as expected. In contrast, larvae carrying two il-2rγc.a mutant alleles lacked mature T cells, although mature B cells were present, albeit less consistently than in wild-type larvae.
It was important to determine whether the role of il-2rγc.a was confined to lymphopoiesis or whether it also exerted an effect on hematopoiesis more generally. Therefore, embryos were stained with O-dianisidine to visualize hemoglobin as a measure of erythrocytes (59) and subjected to WISH with lysozyme (lyz) to mark leukocytes (60). No significant differences in either the extent of O-dianisidine staining (Supplemental Fig. 3A, 3B and data not shown) or the number of lyz+ cells (Supplemental Fig. 3C, 3D) were observed in embryos injected with either il-2rγc.aSS1 or il-2rγc.aSS2 compared with standard control morpholino. Similarly, expression of c-myb, a marker of hematopoietic stem–progenitor cells (61), at 3 dpf, was equivalent between il-2rγc.a morphants and controls (Supplemental Fig. 3E, 3F, and data not shown). Collectively, this result indicates a lymphocyte-specific function for Il-2rγc.a.
IL-2Rγc signaling components are conserved
Having identified a conserved functional role for the Il-2rγc.a paralog in zebrafish lymphopoiesis, it was of interest to determine whether the downstream signaling pathways were also conserved. IL-2Rγc signals via JAK3, with zebrafish previously shown to possess a single, highly conserved jak3 gene (26). This signaling component, which was also strongly expressed in the developing thymus (Supplemental Fig. 4), was targeted by two approaches. In the first approach, embryos were treated with the JAK3-specific inhibitor tofacitinib (41) compared with a DMSO control from 56 hpf, prior to the onset of lymphopoiesis (62) (Fig. 3A–3D). Alternatively, embryos were injected with a morpholino targeting a splice site within jak3 (Fig. 3E–3H). Both approaches substantially reduced expression of rag1 (Fig. 3A–3G), along with other lymphocyte markers (data not shown).
Effect of disruption of downstream pathways on lymphoid cells. (A–G, I–U) WISH analysis of manipulated embryos. Embryos at 56 hpf were bathed in DMSO vehicle control (A and M) or inhibitors for JAK3 (B and C), PI3K (N), or MEK (O) or injected at the 1–8-cell stage with standard control (E and I), jak3 (F), stat5.1 (J), or stat5.2 (K) morpholinos or standard control morpholino plus EGFP mRNA (Q) or il-2rγc.aSS1 morpholino plus EGFP mRNA (R) or mRNA encoding constitutively active Jak3 (S) or Stat5.1 (T) and subjected to WISH with rag1 at the times indicated. The rag1 expression in the thymi is indicated with arrows and quantified relative to eye size for embryos in which Jak3 (D and G), Stat5.1 and Stat5.2 (L), PI3K and MEK (P), or Il-2rγc.a (U) were targeted. Each symbol represents the average ratio for individual embryos, with mean and SEM shown in red. (H) Analysis of jak3, stat5.1, and stat5.2 morpholino targeting. Total RNA was prepared from embryos injected with standard control (Co) or jak3 (Jak3) morpholinos at 4 dpf, or stat5.1 (5.1) or stat5.2 (5.2) morpholinos at 5 dpf, and subjected to RT-PCR with gene specific primers (upper panels) and β-actin (middle panels), with the levels of residual wild-type products indicated for each morphant (lower panel). (V) Phosphorylation of zebrafish Stat5.1 by Jak3. Total cell lysates prepared from HEK293T cells transfected with empty vector or vector expressing Stat5.1 or Jak3 plus Stat5.1 and subjected to Western blot analysis with α-pStat5 (upper panel) or α-GAPDH (lower panel). Statistically significant differences are indicated. *p < 0.05, **p < 0.001. n/s, not significant.
Effect of disruption of downstream pathways on lymphoid cells. (A–G, I–U) WISH analysis of manipulated embryos. Embryos at 56 hpf were bathed in DMSO vehicle control (A and M) or inhibitors for JAK3 (B and C), PI3K (N), or MEK (O) or injected at the 1–8-cell stage with standard control (E and I), jak3 (F), stat5.1 (J), or stat5.2 (K) morpholinos or standard control morpholino plus EGFP mRNA (Q) or il-2rγc.aSS1 morpholino plus EGFP mRNA (R) or mRNA encoding constitutively active Jak3 (S) or Stat5.1 (T) and subjected to WISH with rag1 at the times indicated. The rag1 expression in the thymi is indicated with arrows and quantified relative to eye size for embryos in which Jak3 (D and G), Stat5.1 and Stat5.2 (L), PI3K and MEK (P), or Il-2rγc.a (U) were targeted. Each symbol represents the average ratio for individual embryos, with mean and SEM shown in red. (H) Analysis of jak3, stat5.1, and stat5.2 morpholino targeting. Total RNA was prepared from embryos injected with standard control (Co) or jak3 (Jak3) morpholinos at 4 dpf, or stat5.1 (5.1) or stat5.2 (5.2) morpholinos at 5 dpf, and subjected to RT-PCR with gene specific primers (upper panels) and β-actin (middle panels), with the levels of residual wild-type products indicated for each morphant (lower panel). (V) Phosphorylation of zebrafish Stat5.1 by Jak3. Total cell lysates prepared from HEK293T cells transfected with empty vector or vector expressing Stat5.1 or Jak3 plus Stat5.1 and subjected to Western blot analysis with α-pStat5 (upper panel) or α-GAPDH (lower panel). Statistically significant differences are indicated. *p < 0.05, **p < 0.001. n/s, not significant.
IL-2Rγc/JAK3 stimulates various downstream pathways, including those involving STAT5, PI3K, and ERK (12, 63, 64). Zebrafish possess two STAT5 paralogs, stat5.1 and stat5.2 (65), which were ablated separately with morpholinos (Fig. 3I–3L). Only knockdown of Stat5.1 resulted in a significant, albeit small, decrease in rag1 expression compared with the controls. Finally, zebrafish also possess both PI3K and ERK (66, 67), which were investigated with LY294002, a PI3K inhibitor (68), and PD98059, an inhibitor of MEK, which lies upstream of ERK (69). The extent of rag1 expression was significantly decreased for both inhibitors compared with a DMSO control (Fig. 3M–3P). Notably the JAK3 inhibitor reduced the expression of jak3 (1.5 ± 0.1-fold; p = 0.029), as would be expected because this is a marker of developing lymphocytes (Supplemental Fig. 4). The PI3K inhibitor also reduced expression of pi3kcg (3.0 ± 0.3-fold; p < 0.0001), which is known to be highly expressed in lymphocytes (70). In contrast, the MEK inhibitor had no effect on expression of map2k1 (1.00 ± 0.04-fold; p = 0.7662) and increased the expression of map2k2 (1.2-fold ± 0.02; p = 0.002).
Discussion
The IL-2R family has been studied extensively in mammals, where it has an important role in immune cell development (7), with mutations of the shared IL-2Rγc signaling chain leading to a T−B+ SCID in humans (14) and a T−B− SCID in mice (15). A range of bioinformatic approaches confirmed the presence of two IL-2Rγc paralogs in zebrafish, il-2rγc.a, and il-2rγc.b. The Il-2rγc.a protein showed higher sequence conservation, particularly within the intracellular region. Peak expression of il-2rγc.a also occurred coincidentally with embryonic lymphopoiesis, when it was evident in the developing thymus, whereas il-2rγc.b was not expressed at this site. The il-2rγc.a gene was broadly expressed in adult tissue, but showed highest expression in the spleen, thymus, and kidney with its expression demonstrated in the lymphoid cell populations of these thymus and kidney, whereas il-2rγc.b showed peak expression in the spleen. Knockdown of Il-2rγc.a—but not Il-2rγc.b—substantially inhibited embryonic T lymphopoiesis, but did not affect other aspects of hematopoiesis. Mature T cells were also absent in il-2rγc.a ablated larvae, although mature B cells were present, recapitulating human SCID. This finding suggests that Il-2rγc.a has a conserved role in lymphopoiesis, with the function of Il-2rγc.b currently unknown, consistent with the divergent roles observed for many teleost genes (71–73).
Mammalian IL-2Rγc signals specifically through JAK3 (74–76), which then activates a number of downstream pathways. Zebrafish Il-2rγc.a possessed a conserved JAK3 docking site and was coexpressed with Jak3 in the developing thymus. In addition, ablation of zebrafish Jak3 recapitulated the effects of il-2rγc.a knockdown, consistent with JAK3 loss also causing SCID in both humans and mice (77–79), whereas a constitutively active Jak3 could rescue the effects of il-2rγc.a knockdown. Lying downstream of IL-2Rγc/JAK3 are a number of signaling molecules, including STAT5, PI3K, and ERK, each of which have been shown to contribute to lymphopoiesis. Thus, Stat5a−/Stab5b− mice showed reduced B lymphocytes because of decreased proliferative responses to IL-2 (80–83), Erk1−/Erk2− mice exhibited defective thymic T cell maturation (84, 85), whereas PI3K knockout mice displayed a range of immune defects, including impaired proliferation of T cells in p85α or p110δ knockouts and a reduction in thymocyte numbers in p110γ knockouts (86). Ablation or inhibition of these molecules in zebrafish also produced a reduction in lymphocytes suggesting multiple pathways contribute to lymphocyte development. Of the two zebrafish STAT5 proteins, Stat5.1 appeared to play the key role in this regard, whereas a constitutively active Stat5.1 was able to rescue il-2rγc.a morphants partially, with zebrafish Jak3 shown to phosphorylate Stat5.1.
Having determined the importance of Il-2rγc.a/Jak3 signaling in zebrafish lymphopoiesis, it remains of interest to determine which receptor complexes are involved. In zebrafish, ligand-specific receptor subunits have been identified for all IL-2R family components, except IL-2Rα and IL-9Rα (25). Motifs important for downstream pathways are also conserved on several receptor subunits, including the conserved STAT5 docking sites (YXXV/L/M) on Il-2rβ, Il-4rα, and Il-7rα (87–89), and a docking site for SHC, which contributes to activation of the PI3K and ERK pathways, on Il-2rβ (64) (data not shown). Ligands of the IL-2R family identified in zebrafish include Il-4, Il-15, Il-15–like, and Il-21 (29, 90), with potential Il-2 (91–93) and Il-7 ligands (94) characterized in other teleosts, but not Il-9 (95). Therefore, the phenotypes observed following knockdown of il-2rγc.a and jak3 could potentially be due to defective signaling via several cytokines. Zebrafish Il-4 has been shown to stimulate B cell proliferation and Ab production in adult fish (29), but its role in T cell development has not been determined. However, zebrafish carrying an il-7rα mutation showed decreased rag1 expression at 5 dpf (28), indicating that Il-7–mediated signaling represents an important contributor for early T cell lymphopoiesis. This article also showed contributions from both zebrafish Jak1 and Jak3 in mediating these effects, consistent with the data presented in this study.
Finally, disruption of zebrafish Il-2rγc.a signaling resulted in ablation of lymphopoiesis with similar characteristics to human SCID (74, 96). This result provides important proof-of-principle evidence that the generation of a zebrafish SCID model by specific gene targeting is feasible. Although a number of husbandry challenges need to be overcome, such a SCID zebrafish would represent an invaluable new resource to investigate aspects of immunity and provide a platform for transplantation and cancer studies.
Acknowledgements
We thank the Deakin University Animal House staff for superb aquarium management.
Footnotes
This work was supported by an Australian Postgraduate Award (to R.S.), an International Research Scholarship (to F.B.), an Alfred Deakin Postdoctoral Research Fellowship (to C.L.) from Deakin University, and by the resources of the Australian Research Council Linkage Infrastructure Equipment Fund initiative Fish Works: Collaborative Infrastructure for Zebrafish Research (to A.C.W.).
The online version of this article contains supplemental material.
References
Disclosures
The authors have no financial conflicts of interest.