Combined transient ablation and single cell RNA sequencing reveals the development of medullary thymic epithelial cells

  1. Kristen L Wells
  2. Corey N Miller
  3. Andreas R Gschwind
  4. Wu Wei
  5. Jonah D Phipps
  6. Mark S Anderson  Is a corresponding author
  7. Lars M Steinmetz  Is a corresponding author
  1. Stanford University, United States
  2. University of California, San Francisco, United States
  3. European Molecular Biology Laboratory, Germany

Abstract

Medullary thymic epithelial cells (mTECs) play a critical role in central immune tolerance by mediating negative selection of autoreactive T cells through the collective expression of the peripheral self-antigen compartment, including tissue-specific antigens (TSAs). Recent work has shown that gene expression patterns within the mTEC compartment are remarkably heterogenous and include multiple differentiated cell states. To further define mTEC development and medullary epithelial lineage relationships, we combined lineage tracing and recovery from transient in vivo mTEC ablation with single cell RNA-sequencing in Mus musculus. The combination of bioinformatic and experimental approaches revealed a non-stem transit-amplifying population of cycling mTECs that preceded Aire expression. Based on our findings, we propose a branching model of mTEC development wherein a heterogeneous pool of transit-amplifying cells gives rise to Aire- and Ccl21a-expressing mTEC subsets. We further use experimental techniques to show that within the Aire-expressing developmental branch, TSA expression peaked as Aire expression decreased, implying Aire expression must be established before TSA expression can occur. Collectively, these data provide a higher order roadmap of mTEC development and demonstrate the power of combinatorial approaches leveraging both in vivo models and high-dimensional datasets.

Data availability

RNA-seq data that support the findings of this study have been deposited in the GEO database under accession numbers GSE137699

The following data sets were generated

Article and author information

Author details

  1. Kristen L Wells

    Genetics, Stanford University, Stanford, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-7466-8164
  2. Corey N Miller

    Diabetes Center, University of California, San Francisco, San Francisco, United States
    Competing interests
    The authors declare that no competing interests exist.
  3. Andreas R Gschwind

    Genetics, Stanford University, Stanford, United States
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-0769-6907
  4. Wu Wei

    Genetics, Stanford University, Stanford, United States
    Competing interests
    The authors declare that no competing interests exist.
  5. Jonah D Phipps

    Diabetes Center, University of California, San Francisco, San Francisco, United States
    Competing interests
    The authors declare that no competing interests exist.
  6. Mark S Anderson

    Diabetes Center, University of California, San Francisco, San Francisco, United States
    For correspondence
    Mark.Anderson@ucsf.edu
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3093-4758
  7. Lars M Steinmetz

    Genome Biology Unit, European Molecular Biology Laboratory, Heidelberg, Germany
    For correspondence
    lars.steinmetz@embl.de
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-3962-2865

Funding

National Science Foundation (DGE- 1656518)

  • Kristen L Wells

National Institutes of Health (P01 HG000205)

  • Lars M Steinmetz

National Institutes of Health (R01 AI097457)

  • Mark S Anderson

National Institutes of Health (R01 AI097457)

  • Mark S Anderson

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Ethics

Animal experimentation: Mice were maintained in the University of California San Francisco (UCSF) specific pathogen- free animal facility in accordance with the guidelines established by the Institutional Animal Care and Use Committee (IACUC) and Laboratory Animal Resource Center and all experimental procedures were approved by the Laboratory Animal Resource Center at UCSF. The animal protocol number associated with the study is AN180637-02B.

Copyright

© 2020, Wells et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 4,073
    views
  • 612
    downloads
  • 70
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Kristen L Wells
  2. Corey N Miller
  3. Andreas R Gschwind
  4. Wu Wei
  5. Jonah D Phipps
  6. Mark S Anderson
  7. Lars M Steinmetz
(2020)
Combined transient ablation and single cell RNA sequencing reveals the development of medullary thymic epithelial cells
eLife 9:e60188.
https://doi.org/10.7554/eLife.60188

Share this article

https://doi.org/10.7554/eLife.60188

Further reading

    1. Genetics and Genomics
    2. Microbiology and Infectious Disease
    Chinmaya Jena, Saillesh Chinnaraj ... Nishad Matange
    Research Advance

    Evolution of gene expression frequently drives antibiotic resistance in bacteria. We had previously (Patel and Matange, eLife, 2021) shown that, in Escherichia coli, mutations at the mgrB locus were beneficial under trimethoprim exposure and led to overexpression of dihydrofolate reductase (DHFR), encoded by the folA gene. Here, we show that DHFR levels are further enhanced by spontaneous duplication of a genomic segment encompassing folA and spanning hundreds of kilobases. This duplication was rare in wild-type E. coli. However, its frequency was elevated in a lon-knockout strain, altering the mutational landscape early during trimethoprim adaptation. We then exploit this system to investigate the relationship between trimethoprim pressure and folA copy number. During long-term evolution, folA duplications were frequently reversed. Reversal was slower under antibiotic pressure, first requiring the acquisition of point mutations in DHFR or its promoter. Unexpectedly, despite resistance-conferring point mutations, some populations under high trimethoprim pressure maintained folA duplication to compensate for low abundance DHFR mutants. We find that evolution of gene dosage depends on expression demand, which is generated by antibiotic and exacerbated by proteolysis of drug-resistant mutants of DHFR. We propose a novel role for proteostasis as a determinant of copy number evolution in antibiotic-resistant bacteria.

    1. Developmental Biology
    2. Genetics and Genomics
    Svanhild Nornes, Susann Bruche ... Sarah De Val
    Research Article Updated

    The establishment and growth of the arterial endothelium require the coordinated expression of numerous genes. However, regulation of this process is not yet fully understood. Here, we combined in silico analysis with transgenic mice and zebrafish models to characterize arterial-specific enhancers associated with eight key arterial identity genes (Acvrl1/Alk1, Cxcr4, Cxcl12, Efnb2, Gja4/Cx37, Gja5/Cx40, Nrp1, and Unc5b). Next, to elucidate the regulatory pathways upstream of arterial gene transcription, we investigated the transcription factors binding each arterial enhancer compared to a similar assessment of non-arterial endothelial enhancers. These results found that binding of SOXF and ETS factors was a common occurrence at both arterial and pan-endothelial enhancers, suggesting neither are sufficient to direct arterial specificity. Conversely, FOX motifs independent of ETS motifs were over-represented at arterial enhancers. Further, MEF2 and RBPJ binding was enriched but not ubiquitous at arterial enhancers, potentially linked to specific patterns of behaviour within the arterial endothelium. Lastly, there was no shared or arterial-specific signature for WNT-associated TCF/LEF, TGFβ/BMP-associated SMAD1/5 and SMAD2/3, shear stress-associated KLF4, or venous-enriched NR2F2. This cohort of well-characterized and in vivo-verified enhancers can now provide a platform for future studies into the interaction of different transcriptional and signaling pathways with arterial gene expression.