Skip to main content

Thank you for visiting nature.com. You are using a browser version with limited support for CSS. To obtain the best experience, we recommend you use a more up to date browser (or turn off compatibility mode in Internet Explorer). In the meantime, to ensure continued support, we are displaying the site without styles and JavaScript.

  • Article
  • Published:

TCF-1 and HEB cooperate to establish the epigenetic and transcription profiles of CD4+CD8+ thymocytes

Abstract

Thymocyte development requires a complex orchestration of multiple transcription factors. Ablating either TCF-1 or HEB in CD4+CD8+ thymocytes elicits similar developmental outcomes including increased proliferation, decreased survival, and fewer late Tcra rearrangements. Here, we provide a mechanistic explanation for these similarities by showing that TCF-1 and HEB share ~7,000 DNA-binding sites genome wide and promote chromatin accessibility. The binding of both TCF-1 and HEB was required at these shared sites for epigenetic and transcriptional gene regulation. Binding of TCF-1 and HEB to their conserved motifs in the enhancer regions of genes associated with T cell differentiation promoted their expression. Binding to sites lacking conserved motifs in the promoter regions of cell-cycle-associated genes limited proliferation. TCF-1 displaced nucleosomes, allowing for chromatin accessibility. Importantly, TCF-1 inhibited Notch signaling and consequently protected HEB from Notch-mediated proteasomal degradation. Thus, TCF-1 shifts nucleosomes and safeguards HEB, thereby enabling their cooperation in establishing the epigenetic and transcription profiles of CD4+CD8+ thymocytes.

This is a preview of subscription content, access via your institution

Access options

Buy this article

Prices may be subject to local taxes which are calculated during checkout

Fig. 1: TCF-1 binds accessible chromatin of highly expressed genes.
Fig. 2: TCF-1 binding overlaps with LEF-1, HEB, Runx1, and Ikaros binding.
Fig. 3: TCF-1 or HEB deficiency elicits similar developmental impairments.
Fig. 4: TCF-1 and HEB cobind accessible chromatin of highly expressed genes.
Fig. 5: TCF-1 and HEB promote chromatin accessibility at cobound sites.
Fig. 6: TCF-1 and HEB cooperatively influence DP-thymocyte gene expression.
Fig. 7: TCF-1 promotes HEB stability by inhibiting Notch signaling.

Similar content being viewed by others

Data availability

Sequencing data supporting the findings of this study have been deposited in the Sequence Read Archive (SRA) database under SRA accession number SRP142342. All other relevant data are available from the corresponding authors upon reasonable request.

References

  1. Klein, L., Kyewski, B., Allen, P. M. & Hogquist, K. A. Positive and negative selection of the T cell repertoire: what thymocytes see (and don’t see). Nat. Rev. Immunol. 14, 377–391 (2014).

    Article  CAS  Google Scholar 

  2. Maillard, I., Fang, T. & Pear, W. S. Regulation of lymphoid development, differentiation, and function by the Notch pathway. Annu. Rev. Immunol. 23, 945–974 (2005).

    Article  CAS  Google Scholar 

  3. Weber, B. N. et al. A critical role for TCF-1 in T-lineage specification and differentiation. Nature 476, 63–68 (2011).

    Article  CAS  Google Scholar 

  4. Germar, K. et al. T-cell factor 1 is a gatekeeper for T-cell specification in response to Notch signaling. Proc. Natl Acad. Sci. USA 108, 20060–20065 (2011).

    Article  CAS  Google Scholar 

  5. García-Ojeda, M. E. et al. GATA-3 promotes T-cell specification by repressing B-cell potential in pro-T cells in mice. Blood 121, 1749–1759 (2013).

    Article  Google Scholar 

  6. Scripture-Adams, D. D. et al. GATA-3 dose-dependent checkpoints in early T cell commitment. J. Immunol. 193, 3470–3491 (2014).

    Article  CAS  Google Scholar 

  7. Kueh, H. Y. et al. Asynchronous combinatorial action of four regulatory factors activates Bcl11b for T cell commitment. Nat. Immunol. 17, 956–965 (2016).

    Article  CAS  Google Scholar 

  8. Kumata, K. et al. Development of [11C]MFTC for PET imaging of fatty acid amide hydrolase in rat and monkey brains. ACS Chem. Neurosci. 6, 339–346 (2015).

    Article  CAS  Google Scholar 

  9. Hosokawa, H. et al. Transcription factor PU.1 represses and activates gene expression in early T cells by redirecting partner transcription factor binding. Immunity 48, 1119–1134.e1117 (2018).

    Article  CAS  Google Scholar 

  10. Longabaugh, W. J. R. et al. Bcl11b and combinatorial resolution of cell fate in the T-cell gene regulatory network. Proc. Natl Acad. Sci. USA 114, 5800–5807 (2017).

    Article  CAS  Google Scholar 

  11. Wang, D. et al. The basic helix-loop-helix transcription factor HEBAlt is expressed in pro-T cells and enhances the generation of T cell precursors. J. Immunol. 177, 109–119 (2006).

    Article  CAS  Google Scholar 

  12. Yücel, R., Kosan, C., Heyd, F. & Möröy, T. Gfi1:green fluorescent protein knock-in mutant reveals differential expression and autoregulation of the growth factor independence 1 (Gfi1) gene during lymphocyte development. J. Biol. Chem. 279, 40906–40917 (2004).

    Article  Google Scholar 

  13. David-Fung, E. S. et al. Transcription factor expression dynamics of early T-lymphocyte specification and commitment. Dev. Biol. 325, 444–467 (2009).

    Article  CAS  Google Scholar 

  14. Verbeek, S. et al. An HMG-box-containing T-cell factor required for thymocyte differentiation. Nature 374, 70–74 (1995).

    Article  CAS  Google Scholar 

  15. Schilham, M. W. et al. Critical involvement of Tcf-1 in expansion of thymocytes. J. Immunol. 161, 3984–3991 (1998).

    CAS  PubMed  Google Scholar 

  16. Barndt, R., Dai, M. F. & Zhuang, Y. A novel role for HEB downstream or parallel to the pre-TCR signaling pathway during alpha beta thymopoiesis. J. Immunol. 163, 3331–3343 (1999).

    CAS  PubMed  Google Scholar 

  17. Clevers, H. Wnt/beta-catenin signaling in development and disease. Cell 127, 469–480 (2006).

    Article  CAS  Google Scholar 

  18. Mosimann, C., Hausmann, G. & Basler, K. Beta-catenin hits chromatin: regulation of Wnt target gene activation. Nat. Rev. Mol. Cell Biol. 10, 276–286 (2009).

    Article  CAS  Google Scholar 

  19. Clevers, H. & Nusse, R. Wnt/β-catenin signaling and disease. Cell 149, 1192–1205 (2012).

    Article  CAS  Google Scholar 

  20. Xing, S. et al. Tcf1 and Lef1 transcription factors establish CD8+ T cell identity through intrinsic HDAC activity. Nat. Immunol. 17, 695–703 (2016).

    Article  CAS  Google Scholar 

  21. Johnson, J. L. et al. Lineage-determining transcription factor TCF-1 initiates the epigenetic identity of T cells. Immunity 48, 243–257.e210 (2018).

    Article  CAS  Google Scholar 

  22. Kee, B. L. E and ID proteins branch out. Nat. Rev. Immunol. 9, 175–184 (2009).

    Article  CAS  Google Scholar 

  23. Williams, C. J. et al. The chromatin remodeler Mi-2beta is required for CD4 expression and T cell development. Immunity 20, 719–733 (2004).

    Article  CAS  Google Scholar 

  24. Huang, Z. et al. Transcriptional regulation of CD4 gene expression by T cell factor-1/beta-catenin pathway. J. Immunol. 176, 4880–4887 (2006).

    Article  CAS  Google Scholar 

  25. Ioannidis, V., Beermann, F., Clevers, H. & Held, W. The β-catenin–TCF-1 pathway ensures CD4+CD8+ thymocyte survival. Nat. Immunol. 2, 691–697 (2001).

    Article  CAS  Google Scholar 

  26. Sharma, A., Berga-Bolaños, R. & Sen, J. M. T cell factor-1 controls the lifetime of CD4+ CD8+ thymocytes in vivo and distal T cell receptor α-chain rearrangement required for NKT cell development. PLoS One 9, e115803 (2014).

    Article  Google Scholar 

  27. D’Cruz, L. M., Knell, J., Fujimoto, J. K. & Goldrath, A. W. An essential role for the transcription factor HEB in thymocyte survival, Tcra rearrangement and the development of natural killer T cells. Nat. Immunol. 11, 240–249 (2010).

    Article  Google Scholar 

  28. Steinke, F. C. et al. TCF-1 and LEF-1 act upstream of Th-POK to promote the CD4+ T cell fate and interact with Runx3 to silence CD4 in CD8+ T cells. Nat. Immunol. 15, 646–656 (2014).

    Article  CAS  Google Scholar 

  29. Jones-Mason, M. E. et al. E protein transcription factors are required for the development of CD4+ lineage T cells. Immunity 36, 348–361 (2012).

    Article  CAS  Google Scholar 

  30. Dose, M. et al. β-catenin induces T-cell transformation by promoting genomic instability. Proc. Natl Acad. Sci. USA 111, 391–396 (2014).

    Article  CAS  Google Scholar 

  31. Li, L. et al. A far downstream enhancer for murine Bcl11b controls its T-cell specific expression. Blood 122, 902–911 (2013).

    Article  CAS  Google Scholar 

  32. Winandy, S., Wu, P. & Georgopoulos, K. A dominant mutation in the Ikaros gene leads to rapid development of leukemia and lymphoma. Cell 83, 289–299 (1995).

    Article  CAS  Google Scholar 

  33. Tinsley, K. W. et al. Ikaros is required to survive positive selection and to maintain clonal diversity during T-cell development in the thymus. Blood 122, 2358–2368 (2013).

    Article  CAS  Google Scholar 

  34. Egawa, T., Tillman, R. E., Naoe, Y., Taniuchi, I. & Littman, D. R. The role of the Runx transcription factors in thymocyte differentiation and in homeostasis of naive T cells. J. Exp. Med. 204, 1945–1957 (2007).

    Article  CAS  Google Scholar 

  35. Zhang, J. et al. Harnessing of the nucleosome-remodeling-deacetylase complex controls lymphocyte development and prevents leukemogenesis. Nat. Immunol. 13, 86–94 (2011).

    Article  Google Scholar 

  36. Yu, M. et al. Direct recruitment of polycomb repressive complex 1 to chromatin by core binding transcription factors. Mol. Cell 45, 330–343 (2012).

    Article  CAS  Google Scholar 

  37. Okamura, R. M. et al. Redundant regulation of T cell differentiation and TCRalpha gene expression by the transcription factors LEF-1 and TCF-1. Immunity 8, 11–20 (1998).

    Article  CAS  Google Scholar 

  38. Held, W., Clevers, H. & Grosschedl, R. Redundant functions of TCF-1 and LEF-1 during T and NK cell development, but unique role of TCF-1 for Ly49 NK cell receptor acquisition. Eur. J. Immunol. 33, 1393–1398 (2003).

    Article  CAS  Google Scholar 

  39. Yu, S. et al. The TCF-1 and LEF-1 transcription factors have cooperative and opposing roles in T cell development and malignancy. Immunity 37, 813–826 (2012).

    Article  CAS  Google Scholar 

  40. Whyte, W. A. et al. Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell 153, 307–319 (2013).

    Article  CAS  Google Scholar 

  41. Wang, S. et al. Target analysis by integration of transcriptome and ChIP-seq data with BETA. Nat. Protoc. 8, 2502–2515 (2013).

    Article  CAS  Google Scholar 

  42. Nie, L., Xu, M., Vladimirova, A. & Sun, X. H. Notch-induced E2A ubiquitination and degradation are controlled by MAP kinase activities. EMBO J. 22, 5780–5792 (2003).

    Article  CAS  Google Scholar 

  43. Nie, L. et al. Notch-induced Asb2 expression promotes protein ubiquitination by forming non-canonical E3 ligase complexes. Cell Res. 21, 754–769 (2011).

    Article  CAS  Google Scholar 

  44. Li, X., Gounari, F., Protopopov, A., Khazaie, K. & von Boehmer, H. Oncogenesis of T-ALL and nonmalignant consequences of overexpressing intracellular NOTCH1. J. Exp. Med. 205, 2851–2861 (2008).

    Article  CAS  Google Scholar 

  45. Stritesky, G. L., Jameson, S. C. & Hogquist, K. A. Selection of self-reactive T cells in the thymus. Annu. Rev. Immunol. 30, 95–114 (2012).

    Article  CAS  Google Scholar 

  46. Ben Baruch-Morgenstern, N. et al. Paired immunoglobulin-like receptor A is an intrinsic, self-limiting suppressor of IL-5-induced eosinophil development. Nat. Immunol. 15, 36–44 (2014).

    Article  CAS  Google Scholar 

  47. Adey, A. et al. Rapid, low-input, low-bias construction of shotgun fragment libraries by high-density in vitro transposition. Genome Biol. 11, R119 (2010).

    Article  CAS  Google Scholar 

  48. Jones, M. E. & Zhuang, Y. Acquisition of a functional T cell receptor during T lymphocyte development is enforced by HEB and E2A transcription factors. Immunity 27, 860–870 (2007).

    Article  CAS  Google Scholar 

  49. Wojciechowski, J., Lai, A., Kondo, M. & Zhuang, Y. E2A and HEB are required to block thymocyte proliferation prior to pre-TCR expression. J. Immunol. 178, 5717–5726 (2007).

    Article  CAS  Google Scholar 

  50. Winandy, S., Wu, L., Wang, J. H. & Georgopoulos, K. Pre-T cell receptor (TCR) and TCR-controlled checkpoints in T cell differentiation are set by Ikaros. J. Exp. Med. 190, 1039–1048 (1999).

    Article  CAS  Google Scholar 

  51. Huang, H. et al. Differentiation-dependent interactions between RUNX-1 and FLI-1 during megakaryocyte development. Mol. Cell. Biol. 29, 4103–4115 (2009).

    Article  CAS  Google Scholar 

  52. Heinz, S. et al. Simple combinations of lineage-determining transcription factors prime cis-regulatory elements required for macrophage and B cell identities. Mol. Cell 38, 576–589 (2010).

    Article  CAS  Google Scholar 

  53. Trapnell, C. et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat. Protoc. 7, 562–578 (2012).

    Article  CAS  Google Scholar 

  54. Shen, L., Shao, N., Liu, X. & Nestler, E. ngs.plot: quick mining and visualization of next-generation sequencing data by integrating genomic databases. BMC Genomics 15, 284 (2014).

    Article  Google Scholar 

  55. Freese, N. H., Norris, D. C. & Loraine, A. E. Integrated genome browser: visual analytics platform for genomics. Bioinformatics 32, 2089–2095 (2016).

    Article  CAS  Google Scholar 

  56. Tripathi, S. et al. Meta- and orthogonal Integration of Influenza “OMICs” data defines a role for UBR4 in virus budding. Cell Host Microbe 18, 723–735 (2015).

    Article  CAS  Google Scholar 

  57. Quinlan, A. R. & Hall, I. M. BEDTools: a flexible suite of utilities for comparing genomic features. Bioinformatics 26, 841–842 (2010).

    Article  CAS  Google Scholar 

  58. Kent, W. J., Zweig, A. S., Barber, G., Hinrichs, A. S. & Karolchik, D. BigWig and BigBed: enabling browsing of large distributed datasets. Bioinformatics 26, 2204–2207 (2010).

    Article  CAS  Google Scholar 

  59. Zhang, Y. et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 9, R137 (2008).

    Article  Google Scholar 

  60. Chen, K. et al. DANPOS: dynamic analysis of nucleosome position and occupancy by sequencing. Genome Res. 23, 341–351 (2013).

    Article  CAS  Google Scholar 

Download references

Acknowledgements

We thank H. Kawamoto (RIKEN Research Center for Allergy and Immunology) for reagents; M. Mandal for help with ATAC–seq; B. Kee, M. Clark, A. Khan, and A. Bendelac for advice; and M. Krishnan for helpful suggestions. This work was supported by National Institutes of Health grants R21AI076720, R01AI108682, and U54CA193419, and an ASH bridge grant (to F.G.), and R01CA160436 (to K.K). Further support came from the Chicago Biomedical Consortium (to F.G.) and UL1TR002003 (to M.M.-C.). A.O.E. was supported by an NIH minority supplement; P.S.M. was supported by Institutional NRSA T32 HL07605 and is currently supported as an LLS Fellow; J.Q. was supported by an AAI Careers in Immunology Fellowship; and M.K.O. is a T32HD007009 recipient.

Author information

Authors and Affiliations

Authors

Contributions

A.O.E. designed and performed the experiments, interpreted the experiments, and wrote the manuscript; S.A., L.H., P.S.M., S.M., J.Q., and M.K.O. performed the experiments; M.M.-C. conducted and supervised bioinformatics analyses including the ATAC–seq bioinformatics; K.K. provided advice, interpreted the data, and helped write the manuscript; M.D. and F.G. designed and oversaw the study, interpreted the experiments, and wrote the manuscript.

Corresponding authors

Correspondence to Marei Dose or Fotini Gounari.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Integrated supplementary information

Supplementary Figure 1 WT, Cd4-Cre Tcf7fl/+, and Cd4-Cre Tcf12fl/fl DP thymocytes have comparable fractions of CD71+FSChi cells.

a. Representative CD71 expression versus FSC on gated WT, Cd4-Cre Tcf7fl/+, Cd4-Cre Tcf7fl/– and Cd4-Cre Tcf12fl/fl DP thymocytes. Numbers indicate the percent of cells in each quadrant. b. Cumulative data of CD71+FSChi DP thymocytes in indicated mice (number of biologically independent samples; WT = 9, Cd4-Cre Tcf12fl/fl = 6, Cd4-Cre Tcf7fl/–n = 9, Cd4-Cre Tcf7fl/+n = 3). P-values determined by Ordinary one-way Anova test. Error bars, s.d.

Supplementary Figure 2 HEB binds accessible chromatin of highly expressed genes.

a. K-means clustered heatmap centered on HEB binding sites in WT thymocytes ( ± 1.5 kb) showing enrichment of indicated histone modifications, RNA PolII, and accessibility (ATAC-seq) in DP thymocytes. b. Comparative enrichment histograms of permissive histone modifications (H3Ac, H3K4me1, H3K4me2, H3K4me3 and H3K27ac) and chromatin accessibility (ATAC-seq) at HEB binding sites ( ± 1.5 kb) in K-means clusters shown in a. c. Genomic distribution of HEB peaks in wildtype thymocytes. d. Average expression in DP thymocytes of genes in the indicated groups E = enhancer, P = promoter. Numbers are mean log2 FPKM of means and s.d.. (**** = p ≤ 0.0001, Kruskal-Wallis test). All genes n = 23,360; HEB unbound n = 16,689; HEB bound n = 6,320; HEB p n = 3,519 genes; HEB Active e n = 451 genes; HEB Poised e n = 1255 genes; TCF-1 P + E n = 1,095 genes).

Supplementary Figure 3 Representative TCF-1 and HEB cobinding at promoters and enhancers.

a. ChIP-seq tracks of TCF-1 and HEB co-binding to the promoter of E2f5. b. ChIP-seq tracks of TCF-1 and HEB co-binding to the poised enhancer of Ifng. c. ChIP-seq tracks of TCF-1 and HEB co-binding at the active enhancer of Runx3. d. ChIP-seq tracks of TCF-1 and HEB co-binding at the promoter and enhancer of Rag2. Relative enrichments of H3K4Me1 and H3K4Me3, as well as distance from the transcriptional start site were used in all tracks to identify the promoter, poised enhancer, and active enhancer of the associated gene.

Supplementary Figure 4 TCF-1- and HEB-bound super-clusters overlap at genes involved in T cell development.

a. TCF-1 superclusters representing adjacent TCF-1 binding sites in WT thymocytes within 12.5 kb regions (Homer). Regions were ranked (x-axis) based on the number of sequenced tags (y-axis) per region. The top 271 regions and genes annotated to those regions are shown in red. b. HEB superclusters representing adjacent HEB binding sites in WT thymocytes within 12.5 kb regions (Homer). Regions were ranked (x-axis) based on number of sequenced tags (y-axis) per region. The top 213 regions and genes annotated to these regions are shown in red. c. Overlapping TCF-1 and HEB superclusters and pathways enriched within these clusters (pathways and log10 enrichment values determine by Metascape).

Supplementary Figure 5 TCF-1 and HEB promote chromatin accessibility and affect each other’s DNA binding.

a. Heat map of HEB ChIP-seq enrichment (left) and chromatin accessibility (right) in WT and Cd4-Cre Tcf7fl/– thymocytes ( ± 1.5 kb) around indicated WT HEB binding sites. ATAC-seq was performed in sorted DP thymocytes. b. Comparison of total HEB bound sites in WT and Cd4-Cre Tcf7fl/– thymocytes. c. Histogram representation of HEB enrichment (top) and chromatin accessibility (bottom) around WT HEB binding sites in WT and Cd4-Cre Tcf7fl/– thymocytes. d. Heatmap of TCF-1 ChIP-seq enrichment (left) and chromatin accessibility (right) in WT and Cd4-Cre Tcf12fl/fl thymocytes ( ± 1.5 kb) around indicated WT TCF-1 binding sites. ATAC-seq was performed in sorted DP thymocytes. e. Comparison of total TCF-1 bound sites in WT and Cd4-Cre Tcf12fl/– thymocytes. f. Histogram representation of TCF-1 enrichment (top) and chromatin accessibility (bottom) around WT TCF-1 binding sites in WT and Cd4-Cre Tcf12fl/fl thymocytes.

Supplementary Figure 6 TCF-1 and HEB promote accessibility at sites containing their conserved motifs.

a. Motif enrichment analysis (HOMER) indicating the top 5 enriched motifs (based on p-value, HOMER) in chromatin accessible sites (ATAC-seq) in sorted WT DP thymocytes (51, 452 sites, p≤10e-5). b. Motif enrichment analysis indicating the top 5 enriched motifs in sorted Cd4-Cre Tcf12fl/fl DP thymocytes (42,237 sites, p≤10e-5). c. Motif enrichment analysis indicating the top 5 enriched motifs in chromatin accessible sites in sorted Cd4-Cre Tcf7fl/– DP thymocytes (48,479 sites, p≤10e-5). d. Venn diagram indicating commonly accessible (34,866) and uniquely accessible sites in WT, Cd4-Cre Tcf12fl/fl, and Cd4-Cre Tcf7fl/– DP thymocytes. e. Motif enrichment analysis indicating the top 5 enriched motifs in 7,241 sites only accessible in WT DP thymocytes. f. Motif enrichment analysis indicating the top 5 enriched motifs in 3,485 sites only accessible in Cd4-Cre Tcf12fl/fl DP thymocytes. g. Motif enrichment analysis indicating the top 5 enriched motifs in 4,295 sites only accessible in Cd4-Cre Tcf7fl/– DP thymocytes.

Supplementary Figure 7 TCF-1 is highly enriched at genes associated with Notch signaling and protein ubiquitination.

TCF-1 and HEB binding within core Notch Signaling Cascade genes Dtx1 (a), Lfng (b), and Notch1 (c), that are significantly upregulated in Cd4-Cre Tcf7fl/– DP thymocytes in comparison to WT and Cd4-Cre Tcf12fl/fl DP thymocytes. d. ChIP-seq enrichments of TCF-1, H3K27Ac, H3K4Me1, and chromatin accessibility (ATAC-seq) at TCF-1 binding sites in Notch Signaling/Ubiquitination genes compared to TCF-1-HEB co-bound promoters, poised enhancers, and active enhancers. e. Log2 fold change in chromatin accessibility between WT and Cd4-Cre Tcf7fl/– DP thymocytes at genomic regions described in d.

Supplementary information

Supplementary Information

Supplementary Figures 1–7 and Supplementary Tables 1 and 2

Reporting Summary

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Emmanuel, A.O., Arnovitz, S., Haghi, L. et al. TCF-1 and HEB cooperate to establish the epigenetic and transcription profiles of CD4+CD8+ thymocytes. Nat Immunol 19, 1366–1378 (2018). https://doi.org/10.1038/s41590-018-0254-4

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41590-018-0254-4

This article is cited by

Search

Quick links

Nature Briefing

Sign up for the Nature Briefing newsletter — what matters in science, free to your inbox daily.

Get the most important science stories of the day, free in your inbox. Sign up for Nature Briefing