Draft genome of the aardaker (Lathyrus tuberosus L.), a tuberous legume
BMC Genomic Data volume 23, Article number: 70 (2022)
Lathyrus tuberosus is a nitrogen-fixing member of the Fabaceae which forms protein-rich tubers. To aid future domestication programs for this legume plant and facilitate evolutionary studies of tuber formation, we have generated a draft genome assembly based on Pacific Biosciences sequence reads.
Genomic DNA from L. tuberosus was sequenced with PacBio’s HiFi sequencing chemistry generating 12.8 million sequence reads with an average read length of 14 kb (approximately 180 Gb of sequence data). The reads were assembled to give a draft genome of 6.8 Gb in 1353 contigs with an N50 contig length of 11.1 Mb. The GC content of the genome assembly was 38.3%. BUSCO analysis of the genome assembly indicated a genome completeness of at least 96%. The genome sequence will be a valuable resource, for example, in assessing genomic consequences of domestication efforts and developing marker sets for breeding programs. The L. tuberosus genome will also aid in the analysis of the evolutionary history of plants within the nitrogen-fixing Fabaceae family and in understanding the molecular basis of tuber evolution.
Our current modus operandi for adapting our food production to changing climates is to continuously improve our existing crops to projected future environments. This is a sensible strategy which is of great importance for future food security. A complimentary approach which receives little attention is to select species which have properties we deem useful for future food production and convert these into crops. The lack of research into this approach hampers our ability to make full use of the functional and biological diversity which surrounds us. In addition to diversifying our crop portfolio to improve food supply, one key requirement for a sustainable future is to move away from animal-derived protein to plant-derived protein by growing more protein-rich crops. Peas, beans, and nuts are obvious candidates. However, for wider adoption of protein-rich plant-based diets we also need alternatives to beans and nuts. Protein-rich tubers are good candidates, and one of the plants that produce such protein-rich tubers is Lathyrus tuberosus, a nitrogen fixing member of the Fabaceae which produces tubers with up to 20% protein . L. tuberosus is native to Eurasia and North Africa with a wide geographical distribution extending from Mediterranean to boreal environments. For centuries, it was cultivated or harvested from the wild on a small to medium scale throughout its range for food (leaves, seeds, and tubers) [2,3,4,5]; however, large scale adoption of L. tuberosus as a crop was hampered by poor yields. L. tuberosus is diploid with seven chromosomes and an estimated genome size of 6 Gb [6, 7]. To aid future domestication programs for this legume plant, we have generated a draft genome assembly based on Pacific Biosciences (PacBio) HiFi reads.
L. tuberosus seeds were obtained commercially from Vreeken’s Zaden (Dordrecht, Netherlands) and were grown in Wageningen (Netherlands). Formal identification of the plant material was performed by one of the authors (PJF). A tuber was sent to the botanical garden of the Ruhr-University Bochum (Germany) where it is maintained as a living collection (sample ID: Lathyrus tuberosus NL20). For DNA extraction, shoot tips were collected from the plants grown in Wageningen in September 2020 and immediately frozen in liquid nitrogen. High molecular weight genomic DNA was extracted with the NucleoBond HMW kit (Macherey-Nagel, Germany). PacBio HiFi sequencing libraries were prepared with the SMRTbell Express Template Prep Kit 2.0 (Pacific Biosciences, USA), size-fractionated with the SageELF system (Sage Science, USA) and sequenced on a PacBio Sequel II in six SMRT cells resulting in approximately 12.8 million HiFi reads with an average read length of 14 kb (Table 1, Data set 1, ), providing approximately 30-fold coverage of the genome of L. tuberosus that was previously estimated at 6 Gb using flow cytometry . The sequence reads were assembled with hifiasm , and subsequently the purge_haplotigs tool was used to remove duplicated contigs . The resulting 1668 contigs were searched for putative mitochondrial or chloroplast sequences using BLASTN  against the chloroplast and mitochondrial sequences from Lathyrus sativus and Pisum sativum, respectively [14, 15]. The PacBio sequence reads were mapped against the resulting putative mitochondrial or chloroplast contigs with graphmap  and mapped reads together with the L. sativus or P. sativum organelle genomes were used for similarity-assisted reassembly of the putative mitochondrial or chloroplast contigs with AlignGraph2  resulting in one putative chloroplast contig and five putative mitochondrial contigs. The final L. tuberosus genome assembly (including the reassembled chloroplast and mitochondrial contigs) consists of 1353 contigs with a total length of 6.8 Gb, a contig N50 of 11.1 Mb and a GC content of 38.3% (Table 1, Data file 1, Data set 2, [8, 10]), including five mitochondrial contigs (total length of 476 kb, GC content 45.2%) and a single chloroplast contig (total length of 124 kb, GC content 35.2%) (Table 1, Data file 2, ).
BUSCO (Benchmarking Universal Single-Copy Orthologs) analysis of the assembly showed 96–100% completeness depending on the BUSCO library used for the analysis (Table 1, Data file 3, ). Between 18 and 30% of BUSCO groups were duplicated, most likely due to unphased heterozygous regions (see section Limitations).
The L. tuberosus draft genome sequence will be a valuable resource in future domestication programs, e.g. for developing marker sets for breeding programs. In addition, the L. tuberosus genome will aid in the analysis of the evolutionary history of plants within the nitrogen-fixing Fabaceae family.
The assembly still contains a relatively high degree of duplicated BUSCO groups (up to 30%), most likely due to unphased heterozygous regions. This might also explain the larger assembly size (6.8 Gb) compared to previous estimates by flow cytometry (6 Gb) . L. tuberosus is thought to be an obligate outcrosser, thus obtaining homozygous material is likely to be challenging. Therefore, the duplicated regions might be addressed in future studies, e.g. by using single cell sequencing of pollen grains (gametes) to generate a set of recombinant haploid genotypes, which could be used to phase heterozygous loci (gamete binning, ).
Availability of data and materials
Benchmarking Universal Single-Copy Orthologs
National Center for Biotechnology Information
Sequence Read Archive
Hossaert-Palauqui M, Delbos M. Lathyrus tuberosus L.. Biologie et perspectives d'amélioration. Journal d'agriculture traditionnelle et de botanique appliquée. 1983;30:49–58.
Hanelt P, editor. Mansfeld's encyclopedia of agricultural and horticultural crops: (except ornamentals). Berlin: Springer; 2001.
Dénes A, Papp N, Babai D, Czúcz B, Molnár Z. Wild plants used for food by Hungarian ethnic groups living in the Carpathian Basin. Acta Soc Bot Pol. 2012;81:381–96.
Yildirim E, Dursun A, Turan M. Determination of the nutrition contents of the wild plants used as vegetables in Upper Coruh Valley. Turk J Bot. 2001;25:367–71.
Smýkal P, Erdős L. European tuberous Lathyrus species. Legume Perspect. 2020;19:36–8.
Fisk EL. The chromosomes of Lathyrus tuberosus. Proc Natl Acad Sci U S A. 1931;17:511–3.
Veselý P, Bures P, Smarda P, Pavlícek T. Genome size and DNA base composition of geophytes: the mirror of phenology and ecology? Ann Bot. 2012;109:65–75.
Flood PJ, Nowrousian M, Huettel B, Woehle C, Becker K, Wollenweber TE, Begerow D, Grefen C. Data files for draft genome of the plant Lathyrus tuberosus. figshare. (2022). https://doi.org/10.6084/m9.figshare.19535053.v2.
NCBI Sequence Read Archive. (2022). https://identifiers.org/ncbi/insdc.sra:SRR18139057.
NCBI GenBank. (2022). https://identifiers.org/ncbi/bioproject:PRJNA810344.
Cheng H, Concepcion GT, Feng X, Zhang H, Li H. Haplotype-resolved de novo assembly using phased assembly graphs with hifiasm. Nat Methods. 2021;18:170–5.
Roach MJ, Schmidt SA, Borneman AR. Purge Haplotigs: allelic contig reassignment for third-gen diploid genome assemblies. BMC Bioinform. 2018;19:460.
Altschul SF, Madden TL, Schaffer AA, Zhang J, Zhang Z, Miller W, et al. Gapped BLAST and PSI-BLAST: a new generation of protein database search programs. Nucleic Acids Res. 1997;25:3389–402.
Bogdanova VS, Shatskaya NV, Mglinets AV, Kosterin OE, Vasiliev GV. Discordant evolution of organellar genomes in peas (Pisum L.). Mol Phylogenet Evol. 2021;160:107136. https://doi.org/10.1016/j.ympev.2021.107136.
Magee AM, Aspinall S, Rice DW, Cusack BP, Semon M, Perry AS, et al. Localized hypermutation and associated gene losses in legume chloroplast genomes. Genome Res. 2010;20:1700–10.
Sović I, Šikić M, Wilm A, Fenlon SN, Chen S, Nagarajan N. Fast and sensitive mapping of nanopore sequencing reads with GraphMap. Nat Commun. 2016;7:11307.
Huang S, He X, Wang G, Bao E. AlignGraph2: similar genome-assisted reassembly pipeline for PacBio long reads. Brief Bioinform. 2021;22:bbab022. https://doi.org/10.1093/bib/bbab022.
Campoy JA, Sun H, Goel M, Jiao WB, Folz-Donahue K, Wang N, et al. Gamete binning: chromosome-level and haplotype-resolved genome assembly enabled by high-throughput single-cell sequencing of gamete genomes. Genome Biol. 2020;21:306.
The authors would like to thank Prof. Dr. Karl Köhrer (Heinrich-Heine-Universität Düsseldorf) for support from the Genomics & Transcriptomics Labor as well as the West German Genome Center, the DFG for the funding of the sequencing device (project ID: 388941457), and to acknowledge support by the DFG Open Access Publication Funds of the Ruhr-Universität Bochum. Computational support of the Zentrum für Informations- und Medientechnologie, especially the HPC team (High Performance Computing) at the Heinrich-Heine University is acknowledged.
Open Access funding enabled and organized by Projekt DEAL. MN received funding from the German Research Foundation (DFG, NO 407/7–2). The funders had no role in the design of this study, during its execution, analyses, interpretation of the data, or writing the manuscript.
Ethics approval and consent to participate
Consent for publication
The authors declare no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
Flood, P.J., Nowrousian, M., Huettel, B. et al. Draft genome of the aardaker (Lathyrus tuberosus L.), a tuberous legume. BMC Genom Data 23, 70 (2022). https://doi.org/10.1186/s12863-022-01083-5