A chromosome-scale genome assembly of a diploid alfalfa, the progenitor of autotetraploid alfalfa

Abstract

Alfalfa (Medicago sativa L.) is one of the most important and widely cultivated forage crops. It is commonly used as a vegetable and medicinal herb because of its excellent nutritional quality and significant economic value. Based on Illumina, Nanopore and Hi-C data, we assembled a chromosome-scale assembly of Medicago sativa spp. caerulea (voucher PI464715), the direct diploid progenitor of autotetraploid alfalfa. The assembled genome comprises 793.2 Mb of genomic sequence and 47,202 annotated protein-coding genes. The contig N50 length is 3.86 Mb. This genome is almost twofold larger and contains more annotated protein-coding genes than that of its close relative, Medicago truncatula (420 Mb and 44,623 genes). The more expanded gene families compared with those in M. truncatula and the expansion of repetitive elements rather than whole-genome duplication (i.e., the two species share the ancestral Papilionoideae whole-genome duplication event) may have contributed to the large genome size of M. sativa spp. caerulea. Comparative and evolutionary analyses revealed that M. sativa spp. caerulea diverged from M. truncatula ~5.2 million years ago, and the chromosomal fissions and fusions detected between the two genomes occurred during the divergence of the two species. In addition, we identified 489 resistance (R) genes and 82 and 85 candidate genes involved in the lignin and cellulose biosynthesis pathways, respectively. The near-complete and accurate diploid alfalfa reference genome obtained herein serves as an important complement to the recently assembled autotetraploid alfalfa genome and will provide valuable genomic resources for investigating the genomic architecture of autotetraploid alfalfa as well as for improving breeding strategies in alfalfa.

References

[1] Zhou, Q. et al. MYB transcription factors in alfalfa (Medicago sativa): genome-wide identification and expression analysis under abiotic stresses. PeerJ 7, e7714 (2019).

[2] Liu, Z. et al. Global transcriptome sequencing using the Illumina platform and the development of EST-SSR markers in autotetraploid alfalfa. PLoS ONE 8, e83549 (2013).

[3] Li, X. & Brummer, E. C. Applied genetics and genomics in alfalfa breeding. Agronomy 2, 40-61 (2012).

[4] United States Department of Agriculture-National Agriculture Statistics Service . Crop Production Historical Track Records, April 2018. https://downloads.usda.library.cornell.edu/usda-esmis/files/c534fn92g/6q182n624/v405sd06x/htrcp-04-12-2018.pdf. (2019).

[5] Russelle, M. P. & Birr, A. S. Large-Scale assessment of symbiotic dinitrogen fixation by crops. Agron. J. 96, 1754-1760 (2004).

[6] Bora, K. S. & Sharma, A. Phytochemical and pharmacological potential of Medicago sativa: a review. Pharm. Biol. 49, 211-220 (2011).

[7] Brinker, F. Herb Contraindications and Drug Interactions. Eclectic Medical Publications (Eclectic Medical Publications, 2010).

[8] Malinow, M. R., McLaughlin, P., Naito, H. K., Lewis, L. A. & McNulty, W. P. Effect of alfalfa meal on shrinkage (regression) of atherosclerotic plaques during cholesterol feeding in monkeys. Atherosclerosis 30, 27-43 (1978).

[9] Malinow, M. R., McLaughlin, P. & Stafford, C. Alfalfa seeds: effects on cholesterol metabolism. Experientia 36, 562-564 (1980).

[10] Seida, A., El-Hefnawy, H., Abou-Hussein, D., Mokhtar, F. A. & Abdel-Naim, A. Evaluation of Medicago sativa L. sprouts as antihyperlipidemic and antihyperglycemic agent. Pak. J. Pharm. Sci. 28, 2061-2074 (2015).

[11] Sadeghi, L., Tanwir, F. & Yousefi, B. V. Antioxidant effects of alfalfa can improve iron oxide nanoparticle damage: Invivo and invitro studies. Regul. Toxicol. Pharmacol. 81, 39-46 (2016).

[12] Hong, Y. H., Chao, W. W., Chen, M. L. & Lin, B. F. Ethyl acetate extracts of alfalfa (Medicago sativa L.) sprouts inhibit lipopolysaccharide-induced inflammation in vitro and in vivo. J. Biomed. Sci. 16, 64 (2009).

[13] Zhang, C. & Shi, S. Physiological and Proteomic Responses of Contrasting Alfalfa (Medicago sativa L.) Varieties to PEG-Induced Osmotic Stress. Front. Plant Sci. 9, 242 (2018).

[14] Pan, X. et al. Current Situation and Prospect of Alfalfa Industry. J. Green Sci. Technol. 4, 104-107 (2017). (in Chinese)

[15] Rusk & Nicole Cheap third-generation sequencing. Nat. Methods 6, 244- 244 (2009).

[16] Choi & Chul, S. On the study of microbial transcriptomes using second- and third-generation sequencing technologies. J. Microbiol. 54, 527-536 (2016).

[17] Matheson, N. K., Small, D. M. & Copeland, L. β-d-mannanases in germinating lucerne (alfalfa) seeds. Carbohyd. Res. 82, 325-331 (1980).

[18] Yu, C. Y., Dong, J. G., Hu, S. W. & Xu, A. X. Exposure to trace amounts of sulfonylurea herbicide tribenuron-methyl causes male sterility in 17 species or subspecies of cruciferous plants. BMC Plant Biol. 17, 95 (2017).

[19] 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).

[20] Huang, G. et al. Genome sequence of Gossypium herbaceum and genome updates of Gossypium arboreum and Gossypium hirsutum provide insights into cotton A-genome evolution. Nat. Genet. 52, 516-524 (2020).

[21] Ling, H. Q. et al. Draft genome of the wheat A-genome progenitor Triticum urartu. Nature 496, 87-90 (2013).

[22] Edger, P. P. et al. Author correction: origin and evolution of the octoploid strawberry genome. Nat. Genet. 51, 765 (2019).

[23] Small, E. & Jomphe, M. A synopsis of the genus Medicago (Leguminosae). Can. J. Bot. 67, 3260-3294 (2011).

[24] Hu, J., Fan, J., Sun, Z. & Liu, S. NextPolish: a fast and efficient genome polishing tool for long-read assembly. Bioinformatics 36, 2253-2255 (2020).

[25] Roach, M. J., Schmidt, S. A. & Borneman, A. R. Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinforma. 19, 460 (2018).

[26] Simao, F. A. et al. BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs. Bioinformatics 31, 3210-3212 (2015).

[27] Burton, J. N. et al. Chromosome-scale scaffolding of de novo genome assemblies based on chromatin interactions. Nat. Biotechnol. 31, 1119-1125 (2013).

[28] Kreplav, J. et al. A reference genome for pea provides insight into legume genome evolution. Nat. Genet. 51, 1411-1422 (2019).

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

[30] Karen et al. Strategies for optimizing BioNano and Dovetail explored through a second reference quality assembly for the legume model, Medicago truncatula. BMC Genomics 18, 578 (2017).

[31] Wang, J. et al. Hierarchically aligning 10 legume genomes establishes a family-level genomics platform. Plant Physiol. 174, 284-300 (2017).

[32] Young, N. D., Debellé, F., Oldroyd, G. E. D., Geurts, R. & Roe, B. A. The Medicago genome provides insight into the evolution of rhizobial symbioses. Nature 480, 520-524 (2011).

[33] Cannon, S. B. et al. Multiple polyploidy events in the early radiation of nodulating and nonnodulating legumes. Mol. Biol. Evol. 32, 193-210 (2014).

[34] Lozano, R., Hamblin, M. T., Prochnik, S. & Jannink, J. L. Identification and distribution of the NBS-LRR gene family in the Cassava genome. BMC Genomics 16, 360 (2015).

[35] Xiang, L. et al. Genome-wide comparative analysis of NBS-encoding genes in four Gossypium species. BMC Genomics 18, 292 (2017).

[36] Reddy, M. 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).

[37] Barros, J., Temple, S. & Dixon, R. A. Development and commercialization of reduced lignin alfalfa. Curr. Opin. Biotech. 56, 48-54 (2019).

[38] Shadle, G. et al. Down-regulation of hydroxycinnamoyl CoA: Shikimate hydroxycinnamoyl transferase in transgenic alfalfa affects lignification, development and forage quality. Phytochemistry 68, 1521-1529 (2007).

[39] Bhattarai, K. et al. Agronomic performance and lignin content of HCT down-regulated alfalfa (Medicago sativa L.). Bioenerg. Res. 11, 505-515 (2018).

[40] Kang, M. et al. A chromosome-scale genome assembly of Isatis indigotica, an important medicinal plant used in traditional Chinese medicine. Hortic. Res. 7, 18 (2020).

[41] Tiley, G. P. & Burleigh, J. G. The relationship of recombination rate, genome structure, and patterns of molecular evolution across angiosperms. BMC Evol. Biol. 15, 194 (2015).

[42] Vitte, C. & Bennetzen, J. L. Analysis of retrotransposon structural diversity uncovers properties and propensities in angiosperm genome evolution. Proc. Natl Acad. Sci. USA 103, 17638-17643 (2006).

[43] Wu, S., Han, B. & Jiao, Y. Genetic contribution of paleopolyploidy to adaptive evolution in angiosperms. Mol. Plant 13, 59-71 (2020).

[44] Tang, H. et al. Unraveling ancient hexaploidy through multiply-aligned angiosperm gene maps. Genome Res. 18, 1944-1954 (2008).

[45] Li, S. F. et al. Chromosome-level genome assembly, annotation and evolutionary analysis of the ornamental plant Asparagus setaceus. Hortic. Res. 7, 48 (2020).

[46] Wu, H. et al. A high-quality Actinidia chinensis (kiwifruit) genome. Hortic. Res. 6, 117 (2019).

[47] Barker, D. G. et al. Medicago truncatula, a model plant for studying the molecular genetics of the Rhizobium-legume symbiosis. Plant Mol. Biol. Rep. 8, 40-49 (1990).

[48] Cook, D. R. Medicago truncatula-A model in the making! Curr. Opin. Plant Biol. 2, 301-304 (1999).

[49] Kreplav, J. et al. A reference genome for pea provides insight into legume genome evolution. Nat. Genet. 51, 1411-1422 (2019).

[50] Luo, R. et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. Gigascience 1, 18 (2012).

[51] Lin, H. H., Liao, Y. C. & Dutilh, B. E. Evaluation and validation of assembling corrected Pacbio long reads for microbial genome completion via hybrid approaches. PLoS ONE 10, e0144305 (2015).

[52] Miele, A. & Dekker, J. Mapping cis- and trans-chromatin interaction networks using chromosome conformation capture (3C). Methods Mol. Biol. 464, 105-121 (2009).

[53] Li, H. & Durbin, R. Fast and accurate long-read alignment with Burrows-Wheeler transform. Bioinformatics 26, 589-595 (2010).

[54] Benson, G. Tandem repeats finder: a program to analyze DNA sequences. Nucleic Acids Res. 27, 573-580 (1999).

[55] Tempel, S. Using and understanding RepeatMasker. Methods Mol. Biol. 859, 29-51 (2012).

[56] Jurka, J. et al. Repbase update, a database of eukaryotic repetitive elements. Cytogenet. Genome Res. 110, 462-467 (2005).

[57] Stanke, M. et al. AUGUSTUS: ab initio prediction of alternative transcripts. Nucleic Acids Res. 34, 435-439 (2006).

[58] Majoros, W. H., Pertea, M. & Salzberg, S. L. TigrScan and GlimmerHMM: two open source ab initio eukaryotic gene-finders. Bioinformatics 20, 2878-2879 (2004).

[59] Blanco, E., Parra, G. & Guigó, R. Using geneid to Identify Genes. Curr. Protoc. Bioinforma. 65, 56 (2018).

[60] Birney, E., Clamp, M. & Durbin, R. GeneWise and genomewise. Genome Res. 14, 988-995 (2004).

[61] Keilwagen, J. et al. Using intron position conservation for homology-based gene prediction. Nucleic Acids Res. 44, e89 (2016).

[62] Grabherr, M. G. et al. Trinity: reconstructing a full-length transcriptome without a genome from RNA-Seq data. Nat. Biotechnol. 29, 644-652 (2011).

[63] Haas, B. J. et al. Automated eukaryotic gene structure annotation using EVidenceModeler and the program to assemble spliced alignments. Genome Biol. 9, R7 (2008).

[64] Lowe, T. M. & Eddy, S. R. tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 25, 955-964 (1997).

[65] Camacho, C. et al. BLAST+: architecture and applications. BMC Bioinforma. 10, 421 (2009).

[66] Nawrocki, E. P. & Eddy, S. R. Infernal 1.1: 100-fold faster RNA homology searches. Bioinformatics 29, 2933-2935 (2013).

[67] Kalvari, I. et al. Rfam 13.0: shifting to a genome-centric resource for non-coding RNA families. Nucleic Acids Res. 46, 335-342 (2018).

[68] Jones, P. et al. InterProScan 5: genome-scale protein function classification. Bioinformatics 30, 1236-1240 (2014).

[69] Emms, D. M. & Kelly, S. OrthoFinder: solving fundamental biases in whole genome comparisons dramatically improves orthogroup inference accuracy. Genome Biol. 16, 157 (2015).

[70] Han, M. V., Thomas, G. W., Lugo-Martinez, J. & Hahn, M. W. Estimating gene gain and loss rates in the presence of error in genome assembly and annotation using CAFE 3. Mol. Biol. Evol. 30, 1987-1997 (2013).

[71] Talavera, G. & Castresana, J. Improvement of phylogenies after removing divergent and ambiguously aligned blocks from protein sequence alignments. Syst. Biol. 56, 564-577 (2007).

[72] Alexandros, S. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics 30, 9 (2014).

[73] Sanderson, M. J. r8s: inferring absolute rates of molecular evolution and divergence times in the absence of a molecular clock. Bioinformatics 19, 301-302 (2003).

[74] Wang, Y. et al. MCScanX: a toolkit for detection and evolutionary analysis of gene synteny and collinearity. Nucleic Acids Res. 40, 49 (2012).

[75] Yang, Z. PAML 4: phylogenetic analysis by maximum likelihood. Mol. Biol. Evol. 24, 1586-1591 (2007).

[76] Lupas, A., Van Dyke, M. & Stock, J. Predicting coiled coils from protein sequences. Science 252, 1162-1164 (1991).

[77] McDonnell, A. V., Jiang, T., Keating, A. E. & Berger, B. Paircoil2: improved prediction of coiled coils from sequence. Bioinformatics 22, 356-358 (2006).

[78] Vanholme, R., De Meester, B., Ralph, J. & Boerjan, W. Lignin biosynthesis and its integration into metabolism. Curr. Opin. Biotechnol. 56, 230-239 (2019).

[79] Lampugnani, E. R. et al. Cellulose synthesis-central components and their evolutionary relationships. Trends Plant Sci. 24, 402-412 (2019).

[80] Voorrips, R. E. MapChart: software for the graphical presentation of linkage maps and QTLs. J. Hered. 93, 77-78 (2002).

[81] Kim, D., Langmead, B. & Salzberg, S. L. HISAT: a fast spliced aligner with low memory requirements. Nat. Methods 12, 357-360 (2015).

[82] Pertea, M. et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat. Biotechnol. 33, 290-295 (2015).