The genome of Medicago polymorpha provides insights into its edibility and nutritional value as a vegetable and forage legume

Abstract

Medicago polymorpha is a nutritious and palatable forage and vegetable plant that also fixes nitrogen. Here, we reveal the chromosome-scale genome sequence of M. polymorpha using an integrated approach including Illumina, PacBio and Hi-C technologies. We combined PacBio full-length RNA-seq, metabolomic analysis, structural anatomy analysis and related physiological indexes to elucidate the important agronomic traits of M. polymorpha for forage and vegetable usage. The assembled M. polymorpha genome consisted of 457.53 Mb with a long scaffold N50 of 57.72 Mb, and 92.92% (441.83 Mb) of the assembly was assigned to seven pseudochromosomes. Comparative genomic analysis revealed that expansion and contraction of the photosynthesis and lignin biosynthetic gene families, respectively, led to enhancement of nutritious compounds and reduced lignin biosynthesis in M. polymorpha. In addition, we found that several positively selected nitrogen metabolism-related genes were responsible for crude protein biosynthesis. Notably, the metabolomic results revealed that a large number of flavonoids, vitamins, alkaloids, and terpenoids were enriched in M. polymorpha. These results imply that the decreased lignin content but relatively high nutrient content of M. polymorpha enhance its edibility and nutritional value as a forage and vegetable. Our genomic data provide a genetic basis that will accelerate functional genomic and breeding research on M. polymorpha as well as other Medicago and legume plants.

References

[1] Sato, S. et al. Genome structure of the legume, Lotus japonicus. DNA Res. 15, 227-239 (2008).

[2] Xie, M. et al. A reference-grade wild soybean genome. Nat. Commun. 10, 1216 (2019).

[3] Liu, Q. et al. Assembly and annotation of a draft genome sequence for Glycine latifolia, a perennial wild relative of soybean. Plant J. 95, 71-85 (2018).

[4] Young, N. et al. The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature 480, 520-524 (2011).

[5] Pecrix, Y. et al. Whole-genome landscape of Medicago truncatula symbiotic genes. Nat. Plants 4, 1017-1025 (2018).

[6] Chen, H. et al. Allele-aware chromosome-level genome assembly and efficient transgene-free genome editing for the autotetraploid cultivated alfalfa. Nat. Commun. 11, 2494 (2020).

[7] Shen, C. et al. The Chromosome-level genome sequence of the autotetraploid alfalfa and resequencing of core germplasms provide genomic resources for alfalfa research. Mol. Plant 13, 1250-1261 (2020).

[8] Jain, M. et al. A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.). Plant J. 74, 715-729 (2013).

[9] Kang, Y. J. et al. Genome sequence of mungbean and insights into evolution within Vigna species. Nat. Commun. 5, 5443 (2014).

[10] Bertioli, D. J. et al. The genome sequences of Arachis duranensis and Arachis ipaensis, the diploid ancestors of cultivated peanut. Nat. Genet. 48, 438-446 (2016).

[11] Biazzi, E. et al. Genome-wide association mapping and genomic selection for alfalfa (Medicago sativa) forage quality traits. PLoS ONE 12, e0169234 (2017).

[12] Oba, M. & Allen, M. S. Evaluation of the importance of the digestibility of neutral detergent fiber from forage: effects on dry matter intake and milk yield of dairy cows. J. Dairy. Sci. 82, 589-596 (1999).

[13] Boerjan, W., Ralph, J. & Baucher, M. Lignin biosynthesis. Annu. Rev. Plant Biol. 54, 519-546 (2003).

[14] Reddy, M. S. S. et al. Targeted down-regulation of cytochrome P450 enzymes for forage quality improvement in alfalfa (Medicago sativa L.). Proc. Natl Acad. Sci. USA 102, 16573-16578 (2005).

[15] Bounejmate, M., Robson, A. D. & Beale, P. E. Annual Medicago species in Morocco. II. distribution in relation to soil and climate. Aust. J. Agr. Res. 43, 751-763 (1992).

[16] Ewing, M. A. Annual pastures legumes. A vital component stabilising and rehabilitating low rainfall mediterranean ecosystems. Arid. Soil. Res. Rehab. 13, 327-343 (1999).

[17] Loi, A. et al. The adaptation of Medicago polymorpha to a range of edaphic and environmental conditions: effect of temperature on growth, and acidity stress on nodulation and nod gene induction. Aus. J. Exp. Agric. 33, 25-30 (1993).

[18] Wei, Z. W. et al. A New Medicago polymorpha Cultivar‘ Huaiyang Jinhuacai’. Acta Hortic. Sin. 42, 2335-2336 (2015).

[19] Del Pozo, A. et al. Ecotypic differentiation in Medicago polymorpha L. along an environmental gradient in central Chile. I. Phenology, biomass production and reproductive patterns. Plant Ecol. 159, 119-130 (2002).

[20] Hojilla‐Evangelista, M. P. et al. Extraction, composition, and functional properties of dried alfalfa (Medicago sativa L.) leaf protein. J. Sci. Food Agr. 97, 882-888 (2017).

[21] Guo, D. Improvement of in-rumen digestibility of alfalfa forage by genetic manipulation of lignin O-methyltransferases. Transgenic. Res. 10, 457-464 (2001).

[22] Vanholme, R. et al. Lignin engineering. Curr. Opin. Plant Biol. 11, 278-285 (2008).

[23] Zhao, Q. et al. Laccase is necessary and nonredundant with peroxidase for lignin polymerization during vascular development in Arabidopsis. Plant Cell 25, 3976-3987 (2013).

[24] Sung, J. et al. Photosynthesis, metabolite composition and anatomical structure of Oryza sativa and two wild relatives, O. grandiglumis and O. alta. Rice Sci. 24, 218-227 (2017).

[25] Yu, L., Tan, Y. & Zhu, L. Dietary vitamin B2 intake and breast cancer risk: a systematic review and meta-analysis. Arch. Gynecol. Obstet. 295, 721-729 (2017).

[26] Kurosu, M. Biologically Active Molecules from Soybeans (Intech. Press, 2011).

[27] Panche, A. N., Diwan, A. D. & Chandra, S. R. Flavonoids: an overview. J. Nutr. Sci. 5, e47 (2016).

[28] Ahmad, J. et al. Forage yield potential and quality attributes of alfalfa (Medicago sativa L.) under various agro-management techniques. J. Anim. Plant Sci. 26, 465-474 (2016).

[29] Chan, L. et al. Draft genome analysis provides insights into the fiber yield, crude protein biosynthesis, and vegetative growth of domesticated ramie (Boehmeria nivea L. Gaud). DNA Res. 25, 173-181 (2018).

[30] Belton, J. M. et al. Hi-C: a comprehensive technique to capture the conformation of genomes. Methods 58, 268-276 (2012).

[31] Tang, Q. et al. Comparative transcriptomics of 5 high-altitude vertebrates and their low-altitude relatives. Gigascience 6, 1-9 (2017).

[32] Li, H. & Durbin, R. Fast and accurate short read alignment with Burrows-Wheeler transform. Bioinformatics 25, 1754-1760 (2009).

[33] Akdemir, K. C. & Chin, L. HiCPlotter integrates genomic data with interaction matrices. Genome Biol. 16, 198 (2015).

[34] Pertea, M. et al. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat. Protoc. 11, 1650 (2016).

[35] Wang, D. et al. Identification of nutritional components in black sesame determined by widely targeted metabolomics and traditional Chinese medicines. Molecules 23, 1180 (2018).