Introduction
Practicing traditional food habits and using traditional ingredients are of major importance for maintaining a diet with good nutritional value. South Korea is well known for its fermented foods, particularly banchan (fermented side dishes such as kimchi), which are deeply rooted in Korean food culture. Moreover, Korean cuisine has unique characteristics that are widely accepted to provide various health benefits (), with Korean food culture involving high consumption of vegetables due to the characteristics of its long agricultural history. According to the Korean Ministry of Agriculture definition of the standards and fundamentals of Korean food, one major constraint is that only food prepared with ingredients produced or cultivated in Korea can be considered Korean food. For example, kimchi prepared from imported Chinese cabbage cannot be considered Korean food; the same applies to other banchan. As part of the process of preserving Korean food culture, we have initiated the development of genetic resources for Korean varieties of cucumber (Cucumis sativus var. sativus L.), which is widely cultivated in Korea for both fresh and processed consumption. Cucumber originated in India and spread to other parts of the world through adaptation to various environmental factors and indigenous food habits (). This process has led it to become the sixth most widely cultivated vegetable crop in the world, with 2.1 million hectares under cultivation (FAOSTAT, 2020). South Korea is the 16th largest producer of cucumber in the world, with three major cultivar groups being grown: the Baekdadagi-type, Nakhap-type, and Gasi-type cultivars (). In this study, we aimed to obtain detailed insights into the genetics of cucumber varieties by constructing chromosome-scale genome assemblies for two Korean cucumber inbred lines: JEF (semi-white Baekdadagi-type, mainly used for kimchi and other fermented foods) and KWS (Korean solid green, Nakhap-type, a slicer used fresh for salads and gimbap or Korean cold noodles).
As shown by previous studies of model plants and crops, a single reference genome is inadequate to capture the variation among different genetic lineages. For example, significant structural variation among maize inbred lines has been identified through analysis of multiple genomes (). Furthermore, the cost of assembling multiple genomes has been significantly reduced by third-generation sequencing technologies and computational methods, leading to the construction of chromosome-scale genome assemblies for various crops with the aim of obtaining detailed insights into gene–trait associations (). The first version of the cucumber draft genome was released in 2009 for inbred line 9,930, a lineage of the ‘Chinese Long’ cultivar (); the genome has since been updated to version 3 () and the chromosomal level Northern American cucumber genome published in 2012 (). Further insight into variations among and within varieties has recently been provided by the publication of information on the genome of the pickling cucumber “Borszczagowski” (line B10) (). As the chromosome-scale haploid genome assembly of “Chinese Long” line 9,930 (2n = 2x = 14, haploid number 7) is readily available to the public, we used it as our reference for the construction of chromosome-scale assemblies for the two Korean highly inbred lines.
Value of the Data
These new genomes will serve as an additional genetic resource that can be used as a basis and reference for more detailed study into genetic variation and domestication history among Korean cucumber varieties. In addition, they may be valuable for conducting comparative analysis among and within the species in the genus Cucumis, which could improve the genome selection process in molecular-assisted breeding.
Materials and Methods
Sample Collection and Genomic DNA Extraction
The inbreed lines (i.e., JEF and KWS) are obtained from the leading varieties “Joeun Baekdadagi” and “Gyeoulsal-i Cheongjang” from Fomer Heungnong Seeds Co. After selecting the individual that best characteristics represent of each group in the F2 populations, two inbreeds were raised through self-fertilization. The resulted breed line i.e., JEF is gynoecious, which is semi-white fruit skin color with white spine and KWS is monoecious which is uniform dark green skin color with black spine (Figure 1A). The Cucumis sativus breeding line plants were directly harvested in June 2018 in a field in Wanju, Jeollabuk-do, South Korea (35°90′ N, 127°15′ E), near the National Institute of Agricultural Sciences. Sampled fruits are shown in Figure 1A, and the complete work flow followed in this study is given in Supplementary Figure S1.
FIGURE 1
DNA Sequencing and de novo Genome Assembly
Total DNA was isolated from the samples individually according to sequencing protocols. The isolated DNA was sequenced using two different sequencing systems, PacBio Sequel (Pacific Biosciences, Menlo Park, CA, United States) and Illumina HiSeq 2,500 (Illumina, San Diego, CA, United States), which are widely used in long- and short-read sequencing. For Illumina sequencing, DNA was prepared using the TruSeq Nano DNA Library Prep Kit (Illumina). For PacBio sequencing, DNA was prepared using the SMRTbell Express Template Prep Kit (Pacific Biosciences; catalog no. 101–357–000). The experimental procedures were fully conducted by DNA Link (Seoul, Korea), an authorized service provider in South Korea. The Illumina paired-end sequences were initially subjected to filtering of technical artifacts (i.e., base-calling errors [Phred quality score ≤ Q20]) and adapters using Trimmomatic v. 0.32 (). These Illumina reads were used for error correction of PacBio reads in CLC Assembly Cell v. 5.1.1 (Qiagen, Hilden, Germany). The corrected PacBio reads were used to prepare the initial draft version of the cucumber genomes in FALCON-Unzip v. 0.30, a haplotype assembler program (). Finally, using the RaGOO method (), the genome contigs were clustered and reordered according to their alignment with chromosomal units in the reference genome (‘Chinese Long’ 9,930). The assembled genomes were assessed for completeness using BUSCO v. 4.1.4 with the Viridiplantae_odb10 reference dataset ().
Reference Mapping of Bacterial and Organelle Genes
To prepare a clean reference genome, it was necessary to remove bacterial contamination and organelle genomes from the database. The complete GenBank database, which contains draft and reference genomes of bacteria and organelles (mitochondria and plastids), was used as the reference to determine which reads should be removed from the raw sequences. All reference mapping of preprocessed reads was conducted using Bowtie 2 v. 2.2.8 (). Details regarding reference paths and sizes are given in Supplementary Table S1, and mapping statistics are given in Supplementary Table S2.
Genome Size Estimation
All the Illumina-preprocessed sequences from the paired-end library were subjected to genome size estimation based on k-mers. The k-mer frequencies (k-mer size = 17) were obtained using Jellyfish v. 2.0 (), and the genome size was calculated from the following formulas: genome coverage depth = (k-mer coverage depth × average read length)/(average read length–k-mer size +1); genome size = total base number/genome coverage depth. Here, the k-mer coverage depth is the major peak of the k-mer distribution.
Prediction and Classification of Repeat Regions
Repeat regions in the cucumber genomes were predicted using RepeatModeler (www.repeatmasker.org/RepeatModeler/) and classified into subclasses using the repbase v. 20.08 reference database (www.girinst.org/repbase/) (). Finally, the repeats were masked in the genome using RepeatMasker v. 4.0.5 (www.repeatmasker.org) with RMBlastn v. 2.2.27+. The results are shown in Supplementary Figure S3.
RNA Sequencing
The mRNA library from the collected samples was prepared according to the TruSeq Stranded mRNA Prep Kit protocol (Illumina). The isolated mRNA was sequenced using the Illumina sequencer (Supplementary Tables S4 and Supplementary Tables S5).
Gene Prediction and Annotation
The genes from the cucumber draft genomes were predicted using an in-house gene prediction tool that includes three modules: an evidence-based gene modeler (EVM), an ab-initio gene modeler, and a consensus gene modeler. The Illumina-sequenced transcriptomes were mapped to the respective repeat-masked draft genomes using TopHat, and Trinity v2.5.1 method was used to assemble the transcripts and mark gene structural boundaries (). The ab-initio gene modeler and EVM, which included Exonerate (), Geneid and AUGUSTUS (), were trained with several genomes. The final gene and transcript models were optimized using a consensus gene modeler and annotated using Trinotate v. 3.0.1 ().
Comparative Genome Analysis
Total proteins from the two cucumber genomes were subjected to ortholog analysis to provide insight into the differences between cucumber proteins and those of other plants. In total, 14 genomes from Cucurbitaceae (including the two assembled in this study) were used in the ortholog analysis, with Brassicaceae as outliers (Figure 1D and Supplementary Table S3). The complete proteins of the selected genomes were also subjected to ortholog analysis using OrthoMCL (). The single-copy genes from the given genomes were subjected to phylogenetic tree reconstruction using BEAST (Bayesian Evolutionary Analysis Sampling Trees) to assess the evolutionary time and the degree of similarity among the given genomes (). Furthermore, to assess the gain and loss of genes in the given genomes, the proteins were analyzed using CAFE v. 3.1 ().
Preliminary Analysis Report
Initially, the sizes of the cucumber genomes were estimated to be 267.7 (JEF) and 276.4 MB (KWS) (Figure 1B) based on ∼50 GB of short-read sequences (Table 1A and Supplementary Table S4), but 230.8 MB (JEF) and 231.1 MB (KWS) based on the representative scaffolds assembled from ∼30 GB of error-corrected long-read sequences (Table 1A,B). The N50s of the assembled genomes were 30.5 MB (JEF) and 31.3 MB (KWS), and 40% of the assembled contigs were covered by repeats, in which the long terminal repeat (LTR) elements dominated, accounting for 36% of contigs (Supplementary Figure S3). In total, 25,968 genes were predicted from the JEF genome and 26,011 from KWS, with average sizes of 4,111 and 4,114 bases respectively, and BUSCO scores of 97.88 and 98.35% completeness respectively. (Table 1C). Finally, 66.54% of JEF genes and 65.96% of KWS genes had homologous sequences in GenBank, while 60.25% of JEF genes and 59.82% of KWS genes had gene ontology descriptions (Table 1D). The two genomes were scaffolded onto the reference “Chinese Long” 9,930 genome using the RaGOO method. Overall, these genome assemblies have ∼5 MB of additional bases compared with the reference and similar BUSCO completeness scores, indicating that they are of good quality. Additionally, an average of 99% of both DNA and RNA sequences were mapped to the reference assembly as an additional measure to ensure the quality of the new assemblies (Supplementary Figure S2). The ortholog analysis revealed genome-specific genes, as well as gain and loss of genes, in the selected cucumber genomes (Figure 1C and Supplementary Figure S4). In addition, the RNA samples were collected from five different developmental stages, revealing that both genomes contain genes expressed differentially in different organs or at different stages (Supplementary Figures S5 and Supplementary Figures S6).
TABLE 1
| JEF | KWS | Cucumber_9,930_v3 (GCF_000,004,075.3) | |
|---|---|---|---|
| (A) Sequencing | |||
| Short Read | 72.7 GB (315.10X) | 75.9 GB (328.69X) | |
| Long Reads | 31.7 GB (137.48X) | 37.0 GB (160.38X) | |
| (B) Assembly | |||
| Genome size estimation | 267,736,921 | 276,372,239 | — |
| Total length, bp | 230,754,408 | 231,006,969 | 226,211,662 |
| Total length/Estimation | 86.19% | 83.59% | — |
| No. of contigs | 7 Chr +54 unplaced | 7 Chr +57 unplaced | 7 Chr +77 unplaced |
| Scaffold N50 (Contig N50) | 30,569,742 (7,362,017) | 31,270,087 (8,654,608) | 31,125,843 |
| N (%) | 0.01% | 0.01% | 0.02% |
| GC (%) | 33.23% | 33.25% | 32.82% |
| Repeats (MB) | 93.47 (40.50%) | 94.41 (40.87%) | |
| Repeats against references (9,930) 91.00 (39.44%) | 91.00 (39.44%) | 91.69 (39.69%) | 84.02 (37.14%)–in our method |
| BUSCO | 99.06% | 98.82% | 98.82% |
| (C) Structural annotations | |||
| No. of genes | 25,968 | 26,011 | 24,317 |
| Average gene length (bp) | 4,111.58 | 4,114.07 | 4,068.49 |
| BUSCO (Viridiplantae) | 97.88% | 98.35% | 100.00% |
| (D) Functional annotations | |||
| No. hits | 8,690 (33.46%) | 8,853 (34.04%) | |
| Blast hits | 17,278 (66.54%) | 17,158 (65.96%) | |
| Gene Ontology | 15,645 (60.25%) | 15,561 (59.82%) | |
| KEGG | 3,725 (66.54%) | 3,718 (14.29%) | |
Summary of the sequencing to annotation of the cucumber draft genomes along with the reference.
Statements
Data availability statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Author contributions
YS, MJ, and SS: genome assembly and annotation. KS and SS: manuscript preparation. KL, E-AO, JJ, and J-GK: sampling and sequencing. KS and J-GK: funding and modeling the study.
Funding
This work was supported by the Cooperative Research for Agriculture Science and Technology Development (PJ01343202) of the Rural Development Administration, Republic of Korea.
Conflict of interest
YS, MJ and SS were employed by Insilicogen Inc.
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fgene.2021.733188/full#supplementary-material
References
1
AlongeM.SoykS.RamakrishnanS.WangX.GoodwinS.SedlazeckF. J.et al (2019). RaGOO: Fast and Accurate Reference-Guided Scaffolding of Draft Genomes. Genome Biol.20, 224. 10.1186/s13059-019-1829-6
2
BaoW.KojimaK. K.KohanyO. (2015). Repbase Update, a Database of Repetitive Elements in Eukaryotic Genomes. Mobile DNA6, 11. 10.1186/s13100-015-0041-9
3
BolgerA. M.LohseM.UsadelB. (2014). Trimmomatic: a Flexible Trimmer for Illumina Sequence Data. Bioinformatics30, 2114–2120. 10.1093/bioinformatics/btu170
4
BryantD. M.JohnsonK.DiTommasoT.TickleT.CougerM. B.Payzin-DogruD.et al (2017). A Tissue-Mapped Axolotl De Novo Transcriptome Enables Identification of Limb Regeneration Factors. Cel Rep.18, 762–776. 10.1016/j.celrep.2016.12.063
5
ChinC.-S.PelusoP.SedlazeckF. J.NattestadM.ConcepcionG. T.ClumA.et al (2016). Phased Diploid Genome Assembly with Single-Molecule Real-Time Sequencing. Nat. Methods13, 1050–1054. 10.1038/nmeth.4035
6
HanM. V.ThomasG. W. C.Lugo-MartinezJ.HahnM. W. (2013). 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. 10.1093/molbev/mst100
7
HuangS.et al (2009). The Genome of the Cucumber, Cucumis Sativus L. Nat. Genet.41, 1275–1281. 10.1038/ng.475
8
KimS. H.KimM. S.LeeM. S.ParkY. S.LeeH. J.KangS. A.et al (2016). Korean Diet: Characteristics and Historical Background. J. Ethnic Foods3, 26–31. 10.1016/j.jef.2016.03.002
9
LangmeadB.SalzbergS. L. (2012). Fast Gapped-Read Alignment with Bowtie 2. Nat. Methods9, 357–359. 10.1038/nmeth.1923
10
LiL.StoeckertC. J.Jr.RoosD. S. (2003). OrthoMCL: Identification of Ortholog Groups for Eukaryotic Genomes. Genome Res.13, 2178–2189. 10.1101/gr.1224503
11
LiQ.LiH.HuangW.XuY.ZhouQ.WangS.et al (2019). A Chromosome-Scale Genome Assembly of Cucumber (Cucumis Sativus L.). GigaScience8. 10.1093/gigascience/giz072
12
MarçaisG.KingsfordC. (2011). A Fast, Lock-free Approach for Efficient Parallel Counting of Occurrences of K-Mers. Bioinformatics27, 764–770. 10.1093/bioinformatics/btr011
13
OsipowskiP.PawełkowiczM.WojcieszekM.SkarzyńskaA.PrzybeckiZ.PląderW. (2020). A High-Quality Cucumber Genome Assembly Enhances Computational Comparative Genomics. Mol. Genet. Genomics295, 177–193. 10.1007/s00438-019-01614-3
14
ParkG.ChoiY.JungJ. K.ShimE. J.KangM. Y.SimS. C.et al (2021). Genetic Diversity Assessment and Cultivar Identification of Cucumber (Cucumis Sativus L.) Using the Fluidigm Single Nucleotide Polymorphism Assay. Plants (Basel)10. 10.3390/plants10020395
15
SebastianP.SchaeferH.TelfordI. R. H.RennerS. S. (2010). Cucumber (Cucumis Sativus) and Melon (C. Melo) Have Numerous Wild Relatives in Asia and Australia, and the Sister Species of Melon Is from Australia. Proc. Natl. Acad. Sci.107, 14269–14273. 10.1073/pnas.1005338107
16
SeppeyM.ManniM.ZdobnovE. M. (2019). “BUSCO: Assessing Genome Assembly and Annotation Completeness,” in Gene Prediction: Methods and Protocols. Editor KollmarM. (New York, NY: Springer New York), 227–245. 10.1007/978-1-4939-9173-0_14
17
SlaterG. S. C.BirneyE. (2005). Automated Generation of Heuristics for Biological Sequence Comparison. BMC Bioinformatics6, 31. 10.1186/1471-2105-6-31
18
StankeM.SchöffmannO.MorgensternB.WaackS. (2006). Gene Prediction in Eukaryotes with a Generalized Hidden Markov Model that Uses Hints from External Sources. BMC Bioinformatics7, 62. 10.1186/1471-2105-7-62
19
SuchardM. A.LemeyP.BaeleG.AyresD. L.DrummondA. J.RambautA. (2018). Bayesian Phylogenetic and Phylodynamic Data Integration Using BEAST 1.10. Virus. Evol.4. 10.1093/ve/vey016
20
TaoY.JordanD. R.MaceE. S. (2019). Crop Genomics Goes beyond a Single Reference Genome. Trends Plant Sci.24, 1072–1074. 10.1016/j.tplants.2019.10.001
21
TrapnellC.RobertsA.GoffL.PerteaG.KimD.KelleyD. R.et al (2012). Differential Gene and Transcript Expression Analysis of RNA-Seq Experiments with TopHat and Cufflinks. Nat. Protoc.7, 562. 10.1038/nprot.2012.016
22
YangL.KooD. H.LiY.ZhangX.LuanF.HaveyM. J.et al (2012). Chromosome Rearrangements during Domestication of Cucumber as Revealed by High-Density Genetic Mapping and Draft Genome Assembly. Plant J.71, 895–906. 10.1111/j.1365-313x.2012.05017.x
23
YangN.LiuJ.GaoQ.GuiS.ChenL.YangL.et al (2019). Genome Assembly of a Tropical maize Inbred Line Provides Insights into Structural Variation and Crop Improvement. Nat. Genet.51, 1052–1059. 10.1038/s41588-019-0427-6
Summary
Keywords
Korean cucumber, genome, kimchi, slicer, Cucumis sativus
Citation
Song K, Shin Y, Jung M, Subramaniyam S, Lee KP, Oh E-A, Jeong JH and Kim J-G (2021) Chromosome-Scale Genome Assemblies of Two Korean Cucumber Inbred Lines. Front. Genet. 12:733188. doi: 10.3389/fgene.2021.733188
Received
30 June 2021
Accepted
27 October 2021
Published
19 November 2021
Volume
12 - 2021
Edited by
Yiqun Weng, University of Wisconsin-Madison, United States
Reviewed by
Jarkko Salojärvi, Nanyang Technological University, Singapore
Yuhui Wang, Nanjing Agricultural University, China
Shan Wu, Boyce Thompson Institute, United States
Updates
Copyright
© 2021 Song, Shin, Jung, Subramaniyam, Lee, Oh, Jeong and Kim.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Jeong-Gu Kim, jkim5aug@korea.kr
† These authors have contributed equally to this work
This article was submitted to Plant Genomics, a section of the journal Frontiers in Genetics
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.