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:

A molecular framework underlying the compound leaf pattern of Medicago truncatula

Abstract

Compound leaves show more complex patterns than simple leaves, and this is mainly because of a specific morphogenetic process (leaflet initiation and arrangement) that occurs during their development. How the relevant morphogenetic activity is established and modulated to form a proper pattern of leaflets is a central question. Here we show that the trifoliate leaf pattern of the model leguminous plant Medicago truncatula is controlled by the BEL1-like homeodomain protein PINNATE-LIKE PENTAFOLIATA1 (PINNA1). We identify PINNA1 as a determinacy factor during leaf morphogenesis that directly represses transcription of the LEAFY (LFY) orthologue SINGLE LEAFLET1 (SGL1), which encodes an indeterminacy factor key to the morphogenetic activity maintenance. PINNA1 functions alone in the terminal leaflet region and synergizes with another determinacy factor, the C2H2 zinc finger protein PALMATE-LIKE PENTAFOLIATA1 (PALM1), in the lateral leaflet regions to define the spatiotemporal expression of SGL1, leading to an elaborate control of morphogenetic activity. This study reveals a framework for trifoliate leaf-pattern formation and sheds light on mechanisms generating diverse leaf forms.

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: Phenotypic characterization of pinna1 mutants.
Fig. 2: Positional cloning and characterization of PINNA1.
Fig. 3: Expression pattern of PINNA1 and subcellular localization of the encoded protein.
Fig. 4: PINNA1 negatively regulates SGL1 expression.
Fig. 5: PINNA1 synergizes with PALM1 to regulate lateral leaflet development.
Fig. 6: Model for PINNA1, PALM1 and SGL1 action in the pattern formation of compound leaf in M. truncatula.

Similar content being viewed by others

Data availability

M. truncatula sequence data from this article can be found in the M. truncatula genome version Mt4.0 under the following accession numbers: Medtr3g112290 (PINNA1), Medtr5g014400 (PALM1), Medtr3g098560 (SGL1). The relevant accession number for PINNA1 in GenBank is MN265866. Other data for the current study are available from the paper and its Supplementary Information, the indicated public databases, or from the corresponding authors upon request.

References

  1. Warman, L., Moles, A. T. & Edwards, W. Not so simple after all: searching for ecological advantages of compound leaves. Oikos 120, 813–821 (2011).

    Google Scholar 

  2. Efroni, I., Eshed, Y. & Lifschitz, E. Morphogenesis of simple and compound leaves: a critical review. Plant Cell 22, 1019–1032 (2010).

    CAS  PubMed  PubMed Central  Google Scholar 

  3. Bar, M. & Ori, N. Compound leaf development in model plant species. Curr. Opin. Plant Biol. 23, 61–69 (2015).

    PubMed  Google Scholar 

  4. Hareven, D., Gutfinger, T., Parnis, A., Eshed, Y. & Lifschitz, E. The making of a compound leaf: genetic manipulation of leaf architecture in tomato. Cell 84, 735–744 (1996).

    CAS  PubMed  Google Scholar 

  5. Hay, A. & Tsiantis, M. The genetic basis for differences in leaf form between Arabidopsis thaliana and its wild relative Cardamine hirsuta. Nat. Genet. 38, 942–947 (2006).

    CAS  PubMed  Google Scholar 

  6. Bharathan, G. et al. Homologies in leaf form inferred from KNOXI gene expression during development. Science 296, 1858–1860 (2002).

    CAS  PubMed  Google Scholar 

  7. Vollbrecht, E., Veit, B., Sinha, N. & Hake, S. The developmental gene Knotted-1 Is a member of a maize homeobox gene family. Nature 350, 241–243 (1991).

    CAS  PubMed  Google Scholar 

  8. Lincoln, C., Long, J., Yamaguchi, J., Serikawa, K. & Hake, S. A knotted1-like homeobox gene in Arabidopsis is expressed in the vegetative meristem and dramatically alters leaf morphology when overexpressed in transgenic plants. Plant Cell 6, 1859–1876 (1994).

    CAS  PubMed  PubMed Central  Google Scholar 

  9. Hofer, J. et al. UNIFOLIATA regulates leaf and flower morphogenesis in pea. Curr. Biol. 7, 581–587 (1997).

    CAS  PubMed  Google Scholar 

  10. Gourlay, C. W., Hofer, J. M. & Ellis, T. H. Pea compound leaf architecture is regulated by interactions among the genes UNIFOLIATA, COCHLEATA, AFILA, and TENDRIL-LESS. Plant Cell 12, 1279–1294 (2000).

    CAS  PubMed  PubMed Central  Google Scholar 

  11. Champagne, C. E. et al. Compound leaf development and evolution in the legumes. Plant Cell 19, 3369–3378 (2007).

    CAS  PubMed  PubMed Central  Google Scholar 

  12. Wang, H. et al. Control of compound leaf development by FLORICAULA/LEAFY ortholog SINGLE LEAFLET1 in Medicago truncatula. Plant Physiol. 146, 1759–1772 (2008).

    CAS  PubMed  PubMed Central  Google Scholar 

  13. Hofer, J. et al. Tendril-less regulates tendril formation in pea leaves. Plant Cell 21, 420–428 (2009).

    CAS  PubMed  PubMed Central  Google Scholar 

  14. Wang, Z. et al. Multiple components are integrated to determine leaf complexity in Lotus japonicus. J. Integr. Plant Biol. 55, 419–433 (2013).

    CAS  PubMed  Google Scholar 

  15. Hofer, J. M. & Noel Ellis, T. H. Developmental specialisations in the legume family. Curr. Opin. Plant Biol. 17, 153–158 (2014).

    PubMed  Google Scholar 

  16. Jeong, S. C., Kim, J. H. & Bae, D. N. Genetic analysis of the Lf1 gene that controls leaflet number in soybean. Theor. Appl. Genet. 130, 1685–1692 (2017).

    CAS  PubMed  Google Scholar 

  17. Zhou, C. et al. Transforming compound leaf patterning by manipulating REVOLUTA in Medicago truncatula. Plant J. 100, 562–571 (2019).

    CAS  PubMed  Google Scholar 

  18. Jiao, K. et al. Genetic control of compound leaf development in the mungbean (Vigna radiata L.). Hortic. Res. 6, 23 (2019).

    PubMed  PubMed Central  Google Scholar 

  19. Israeli, A. et al. Multiple auxin-response regulators enable stability and variability in leaf development. Curr. Biol. 29, 1746–1759 (2019).

    CAS  PubMed  Google Scholar 

  20. Naz, A. A. et al. Trifoliate encodes an MYB transcription factor that modulates leaf and shoot architecture in tomato. Proc. Natl Acad. Sci. USA 110, 2401–2406 (2013).

    CAS  PubMed  Google Scholar 

  21. Vlad, D. et al. Leaf shape evolution through duplication, regulatory diversification, and loss of a homeobox gene. Science 343, 780–783 (2014).

    CAS  PubMed  Google Scholar 

  22. Bar, M. et al. CLAUSA is a MYB transcription factor that promotes leaf differentiation by attenuating cytokinin signaling. Plant Cell 28, 1602–1615 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

  23. Blein, T. et al. A conserved molecular framework for compound leaf development. Science 322, 1835–1839 (2008).

    CAS  PubMed  Google Scholar 

  24. Kierzkowski, D. et al. A growth-based framework for leaf shape development and diversity. Cell 177, 1405–1418 (2019).

    CAS  PubMed  PubMed Central  Google Scholar 

  25. Chen, J. et al. Control of dissected leaf morphology by a Cys(2)His(2) zinc finger transcription factor in the model legume Medicago truncatula. Proc. Natl Acad. Sci. USA 107, 10754–10759 (2010).

    CAS  PubMed  Google Scholar 

  26. Tadege, M. et al. Large-scale insertional mutagenesis using the Tnt1 retrotransposon in the model legume Medicago truncatula. Plant J. 54, 335–347 (2008).

    CAS  PubMed  Google Scholar 

  27. Reiser, L. et al. The bell1 gene encodes a homeodomain protein involved in pattern-formation in the Arabidopsis ovule primordium. Cell 83, 735–742 (1995).

    CAS  PubMed  Google Scholar 

  28. Smith, H. M. S. & Hake, S. The interaction of two homeobox genes, BREVIPEDICELLUS and PENNYWISE, regulates internode patterning in the Arabidopsis inflorescence. Plant Cell 15, 1717–1727 (2003).

    CAS  PubMed  PubMed Central  Google Scholar 

  29. Kumar, R. et al. The Arabidopsis BEL1-LIKE HOMEODOMAIN proteins SAW1 and SAW2 act redundantly to regulate KNOX expression spatially in leaf margins. Plant Cell 19, 2719–2735 (2007).

    CAS  PubMed  PubMed Central  Google Scholar 

  30. Kimura, S., Koenig, D., Kang, J., Yoong, F. Y. & Sinha, N. Natural variation in leaf morphology results from mutation of a novel KNOX gene. Curr. Biol. 18, 672–677 (2008).

    CAS  PubMed  Google Scholar 

  31. Tsuda, K. et al. KNOTTED1 cofactors, BLH12 and BLH14, regulate internode patterning and vein anastomosis in maize. Plant Cell 29, 1105–1118 (2017).

    CAS  PubMed  PubMed Central  Google Scholar 

  32. Chou, K. C. & Shen, H. B. Cell-PLoc: a package of Web servers for predicting subcellular localization of proteins in various organisms. Nat. Protoc. 3, 153–162 (2008).

    CAS  PubMed  Google Scholar 

  33. Tsuda, K. & Hake, S. Diverse functions of KNOX transcription factors in the diploid body plan of plants. Curr. Opin. Plant Biol. 27, 91–96 (2015).

    CAS  PubMed  Google Scholar 

  34. Magnani, E. & Hake, S. KNOX lost the OX: the Arabidopsis KNATM gene defines a novel class of KNOX transcriptional regulators missing the homeodomain. Plant Cell 20, 875–887 (2008).

    CAS  PubMed  PubMed Central  Google Scholar 

  35. Peng, J. et al. Regulation of compound leaf development in Medicago truncatula by Fused Compound Leaf1, a class M KNOX gene. Plant Cell 23, 3929–3943 (2011).

    CAS  PubMed  PubMed Central  Google Scholar 

  36. Liu, Y. et al. BEL1-LIKE HOMEODOMAIN6 and KNOTTED ARABIDOPSIS THALIANA7 interact and regulate secondary cell wall formation via repression of REVOLUTA. Plant Cell 26, 4843–4861 (2014).

    CAS  PubMed  PubMed Central  Google Scholar 

  37. Lee, J. H., Lin, H., Joo, S. & Goodenough, U. Early sexual origins of homeoprotein heterodimerization and evolution of the plant KNOX/BELL family. Cell 133, 829–840 (2008).

    CAS  PubMed  Google Scholar 

  38. Furumizu, C., Alvarez, J. P., Sakakibara, K. & Bowman, J. L. Antagonistic roles for KNOX1 and KNOX2 genes in patterning the land plant body plan following an ancient gene duplication. PLoS Genet. 11, e1004980 (2015).

    PubMed  PubMed Central  Google Scholar 

  39. Chen, J. et al. Conserved genetic determinant of motor organ identity in Medicago truncatula and related legumes. Proc. Natl Acad. Sci. USA 109, 11723–11728 (2012).

    CAS  PubMed  Google Scholar 

  40. Jefferson, R. A., Kavanagh, T. A. & Bevan, M. W. Gus fusions; beta-glucuronidase as a sensitive and versatile gene fusion marker in higher plants. EMBO J. 6, 3901–3907 (1987).

    CAS  PubMed  PubMed Central  Google Scholar 

  41. Shao, J. L., Long, Y. S., Chen, G., Xie, J. & Xu, Z. F. The reversed terminator of octopine synthase gene on the Agrobacterium Ti plasmid has a weak promoter activity in prokaryotes. Mol. Biol. Rep. 37, 2157–2162 (2010).

    CAS  PubMed  Google Scholar 

  42. Kapila, J., DeRycke, R., VanMontagu, M. & Angenon, G. An Agrobacterium-mediated transient gene expression system for intact leaves. Plant Sci. 124, 227–227 (1997).

    CAS  Google Scholar 

  43. Coen, E. S. et al. Floricaula—a homeotic gene required for flower development in Antirrhinum majus. Cell 63, 1311–1322 (1990).

    CAS  PubMed  Google Scholar 

  44. Cosson, V., Durand, P., d’Erfurth, I., Kondorosi, A. & Ratet, P. Medicago truncatula transformation using leaf explants. Methods Mol. Biol. 343, 115–127 (2006).

    CAS  PubMed  Google Scholar 

  45. Zhou, Y. et al. Control of plant stem cell function by conserved interacting transcriptional regulators. Nature 517, 377–380 (2015).

    CAS  PubMed  Google Scholar 

  46. Qi, X. T., Chai, X. Q., Liu, J. & Chai, T. Y. A new method for EMSA by modifying DIG high prime DNA labeling and detection starter kit II. Hereditas 28, 721–725 (2006).

    CAS  PubMed  Google Scholar 

  47. Zheng, S. Y. et al. OsAGO2 controls ROS production and the initiation of tapetal PCD by epigenetically regulating OsHXK1 expression in rice anthers. Proc. Natl Acad. Sci. USA 116, 7549–7558 (2019).

    CAS  PubMed  Google Scholar 

  48. Gendrel, A. V., Lippman, Z., Martienssen, R. & Colot, V. Profiling histone modification patterns in plants using genomic tiling microarrays. Nat. Methods 2, 213–218 (2005).

    CAS  PubMed  Google Scholar 

  49. Ikeda, T. et al. BELL1-like homeobox genes regulate inflorescence architecture and meristem maintenance in rice. Plant J. 98, 465–478 (2019).

    CAS  PubMed  Google Scholar 

  50. Thompson, J. D., Gibson, T. J. & Higgins, D. G. Multiple sequence alignment using ClustalW and ClustalX. Curr. Protoc. Bioinformatics 00, 2.3.1–2.3.22 (2003).

    Google Scholar 

  51. Trifinopoulos, J., Nguyen, L. T., von Haeseler, A. & Minh, B. Q. W-IQ-TREE: a fast online phylogenetic tool for maximum likelihood analysis. Nucleic Acids Res. 44, W232–W235 (2016).

    CAS  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Z. Gu for assistance with SEM; Y. Hu for providing the plasmid pGREEN800IILUC; members of the J.C. and Z. Xu laboratory for their valuable input; and S. Khederzadeh for reading this manuscript. These studies were founded by National Natural Science Foundation of China grants U1702234, 31471171 and XDB27030106 (to J.C.), and the Yunnan Applied Basic Research Project 2016FB041. J.C. is also supported by the One Hundred Talent Project of the Chinese Academy of Sciences, High-end Scientific and Technological Talents in Yunnan Province (2015HA031, 2015HA032), the CAS 135 program and Core Botanical Gardens, Chinese Academy of Sciences.

Author information

Authors and Affiliations

Authors

Contributions

J.C. designed research; L.H. and Yu L. performed most of the experiments; J.Q., Y.X. and Ye L. participated in RNA in situ hybridization; X. Zhang, Y. Li and X. Zheng participated in genetic transformation; H.H., Y.M., S.Z. and Q.B. participated in EMSA, yeast two-hybrid, BiFC, pull-down and ChIP assays; J.W., K.S.M., M.T. and Y.X. contributed new reagent and analytic tools; L.H. and J.C. analysed data; L.H. and J.C. wrote the manuscript with input from M.T., Ye L., J.Q., B.Z. and D.W.

Corresponding author

Correspondence to Jianghua Chen.

Ethics declarations

Competing interests

The authors declare no competing interests.

Additional information

Peer review information Nature Plants thanks Yuling Jiao, Miltos Tsiantis, Hirokazu Tsukaya and the other, anonymous, reviewers for their contribution to the peer review of this work.

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

Extended data

Extended Data Fig. 1 Phenotypic analysis of WT and the pinna1-1 mutant.

a,b, Images of 3-week-old (a) and 7-week-old (b) WT and pinna1-1 plants. These experiments were repeated independently three times with similar results. Scale bars, 4 cm. c, Phenotypic variation of pinna1-1 leaves with a pie chart showing the percentage distribution. This experiment was repeated independently three times with similar results. Scale bar, 1 cm. d, SEM images of shoot apices with P1 to P3 leaf primordia. Yellow dotted circles indicate no difference in proximal regions of the P3 TL primordia between WT and pinna1-1. This experiment was repeated independently three times with similar results. Scare bares, 20 μm.

Extended Data Fig. 2 Molecular cloning of PINNA1 gene and cosegregation analysis.

a, Simplified schema of the Tnt1 insertion site-based mapping strategy. Red and blue short-lines represent the recovered Tnt1 FSTs of two different mutant lines. b, Images of representative WT, pinna1-1, -2, -3, -4 and -5 leaves. This experiment was repeated independently three times with similar results. Scale bar, 1 cm. c, The detailed positions (red numbers indicate position from the first ATG codon) of Tnt1 insertion sites in pinna1 alleles. Yellow shaded letters indicate the exon sequences. The consensus GT/AG splice sequences are boxed. d, PINNA1 gene structure, the Tnt1 insertion site of pinna1-1 and the primers (arrowheads) used for genotyping. e, Cosegregation analysis of a F2 population derived from a cross of pinna1-1 x WT (R108). 41 out of 172 F2 individuals showing mutant phenotype were homozygous for Tnt1 insertion. The primer pair Fw100-RvTAG100 was used for detecting intact PINNA1 genomic fragments (upper); the primer pair Fw394-TntR2 was used for verifying the existence of Tnt1 insertion within the PINNA1 genomic region (lower). This experiment was repeated independently three times with similar results.

Extended Data Fig. 3 Genetic complementation of pinna1.

a, 3-week-old plants of pinna1-2 and two 35S::PINNA1-GFP/pinna1-2 transgenic lines with partially rescued (#2-3) and completely rescued (#2-1) phenotypes. This experiment was repeated independently three times with similar results. Scare bar, 2 cm. b,c, Normal RT-PCR (b) and RT-qPCR (c) analysis of PINNA1 expression level in young leaves of seven T0 independent lines (35S::PINNA1-GFP) compared to pinna1-1 and pinna1-2. The normal RT-PCR experiment was repeated independently two times with similar results. The RT-qPCR experiment was repeated three times using independent biological replicates with similar results, and one representative result is shown (mean ± SD, n = 3). P, 35S::PINNA1-GFP plasmid. M, molecular weight markers. PINNA1 gene structure, the positions of Tnt1 insertions in pinna1-1 and pinna1-2, and the RT-qPCR primer-pair are shown in the top panel. d, A 7-week-old 35S::PINNA1-GFP/pinna1-2 transgenic plant (#2-12) with a completely rescued phenotype. This experiment was repeated independently three times with similar results. Scare bar, 5 cm. e, Western blot analysis of T1 progeny from the 35S:PINNA1-GFP/pinna1-2 transgenic plant (#2-12) using an anti-GFP antibody. A DR5::GFP transgenic plant used as positive control. 3 to 10, independent T1 plants. This experiment was repeated independently three times with similar results. f, A 7-week-old pinna1-2/pPINNA1:GFP-PINNA1 transgenic line (#G1) with a completely rescued phenotype. This experiment was repeated independently three times with similar results. Scare bar, 5 cm. g, Normal RT-PCR analysis of PINNA1 expression level in shoot buds of two independent pinna1-2/pPINNA1::GFP-PINNA1 transgenic lines. P, pPINNA1:GFP-PINNA1 plasmid. This experiment was repeated independently three times with similar results.

Source data

Extended Data Fig. 4 Maximum likelihood tree of TALE-superfamily homeodomain proteins from several angiosperms.

Legume species (M. truncatula, pea, soybean and L. japonicas), related core eudicots (S. lycopersicum, C. hirsute and A. thaliana), basal eudicots (Aquilegia coerulea), monocots (Oryza sativa and Zea mays), and basal angiosperms (Amborella trichopoad) are chosen for analysis. This is a detailed representation of the tree shown in Fig. 2e. Bootstrap values are out of 2,000 replicate.

Extended Data Fig. 5 Expression pattern of PINNA1 during leaf development.

Longitudinal serial sections of WT apex showing PINNA1 expression in leaf primordium but not in the shoot apical meristem. This experiment was repeated independently three times with similar results. TL, terminal leaflet; LL, lateral leaflet. Scale bars, 50 μm.

Extended Data Fig. 6 Genetic interaction between pinna1-1 and sgl1-2 in the elp1-3 genetic background.

elp1-3 developed normal trifoliate form leaves a, while pinna1-1 elp1-3 developed typical pinnate form leaves (indicated as “Pinnate”) b. sgl1-2 elp1-3 developed simple-like leaves having a single leaflet c, and pinna1-1 sgl1-2 elp1-3 strongly resembled the sgl1-2 elp1-3 by mostly developing simple-like leaves of single leaflet (indicated as “SL”) (except a few dissected leaves having two or three clustered leaflets that indicated as DL) d. SL, single leaflet; DL, dissected leaf. Proportions of leaf types in each genotype are shown in pie charts. Blue, trifoliate; red, pinnate; green, single; orange, dissected leaves with clustered leaflets. These experiments were repeated independently three times with similar results. Scale bars, 1 cm.

Extended Data Fig. 7 RNA in situ hybridization of SGL1 in palm1 elp1 and pinna1 palm1 elp1.

a, Serial longitudinal sections of the palm1-3 elp1-3 P4 leaf primordium. b,c, Serial longitudinal sections of the pinna1-1 palm1-3 elp1-3 leaf primordia at P4 (b) and at P5 (c) stages. LP, leaf primordia. These experiments were repeated independently three times with similar results. Scale bars, 50 μm.

Extended Data Fig. 8 BiFC analysis of the in vivo interaction between PINNA1 and PALM1.

Epidermal cells of the N. benthamiana leaf were co-transformed with GFPn-PINNA1/GFPc-PALM1, GFPn-PINNA1/GFPc, or GFPn/GFPc-PALM1. GFPn, N-terminal portion of GFP; GFPc, C-terminal portion of GFP. When GFPn- and GFPc-tagged proteins are in close proximity, the GFPn and GFPc portions interact and produce GFP fluorescence. This experiment was repeated independently three times with similar results. Scale bars, 10 μm.

Extended Data Fig. 9 Quantitative analysis of PALM1 and PINNA1 expression level in pinna1 and palm1 respectively.

a, RT-qPCR of PALM1 expression in the lateral leaflet parts of the first emerged young leaves (P7) from WT and pinna1-1 shoots. Data shown as mean ± SD (n = 3 biologically independent replicates). b, RT-qPCR of PINNA1 expression in different parts of the first emerged young leaves (P7) from WT and palm1-3 shoots. Data shown as mean ± SD (n = 3 biologically independent replicates).

Source data

Extended Data Fig. 10 RNA in situ hybridization of SGL1 in WT, palm1 and pinna1.

Shown are serial longitudinal sections of the P3 leaf primordium. An upregulated SGL1 transcription was found in palm1-3 P3 LL primordia (cyan arrows). This experiment was repeated independently three times with similar results.

Supplementary information

Supplementary Information

Supplementary Figs. 1–5 and Table 1.

Reporting Summary

Source data

Source Data Fig. 1

Statistical source data.

Source Data Fig. 3

Statistical source data.

Source Data Fig. 4

Statistical source data.

Source Data Fig. 4

Unprocessed gels.

Source Data Fig. 5

Statistical source data.

Source Data Fig. 5

Unprocessed western blots.

Source Data Extended Data Fig. 2

Unprocessed gels.

Source Data Extended Data Fig. 3

Statistical source data.

Source Data Extended Data Fig. 3

Unprocessed western blots and gels.

Source Data Extended Data Fig. 9

Statistical source data.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

He, L., Liu, Y., He, H. et al. A molecular framework underlying the compound leaf pattern of Medicago truncatula. Nat. Plants 6, 511–521 (2020). https://doi.org/10.1038/s41477-020-0642-2

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1038/s41477-020-0642-2

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