Abstract
Spinach (Spinacia oleracea L.) is an economically important vegetable crop. Here we describe the complete mitochondrial DNA sequence of spinach, which has a length of 329,613 bp and a GC content of 43.4%. It is separated by a pair of directly repeated elements of 7286 bp, to form a large single copy region of 229,375 bp and a small single copy region of 85,666 bp. The genome contains 29 protein-coding genes, 4 pseudogenes, 24 tRNA genes, and 3 rRNA genes. A phylogenetic analysis revealed that spinach was closely related to Beta vulgaris (sugar beet), both belonging to the Amaranthaceae family.
Spinach (Spinacia oleracea L.) is widely cultivated as an economically important green leafy vegetable crop for fresh consumption and processing (van Treuren et al. Citation2011), and is a rich source of several important health-promoting nutrients including carotenoids, folates and vitamin C. Spinach is a member of the Amaranthaceae family, which includes several other important crops such as beets and quinoa. Spinach is native to central Asia and thought to have originated in Persia (Iran), and the genus Spinacia contains two wild species, S. turkestanica and S. tetrandra (Ryder Citation1979). However, no complete mitochondrial genome sequences are currently available for Spinacia species.
In this study, we determined the complete mitochondrial genome sequence of S. oleracea. Genomic DNA was extracted from young leaves of an inbred line (SP75) that was collected from Development and Collaborative Innovation Center of Plant Germplasm Resources at Shanghai Normal University campus (30°50′29.4”N, 121°31′16.7”E), and used to construct Illumina DNA libraries. High-throughput sequencing generated ∼169 Gb high-quality cleaned sequences, which were used previously to assemble a draft nuclear genome (Xu et al. Citation2017). Meantime, these sequences were also used for de novo assembly using SOAPdenovo2 (Luo et al. Citation2012) and the final assembled contig that was from S. oleracea mitochondrion was identified based on the k-mer depth analysis and the comparison to other plant mitogenome sequences. The mitochondrial genome of S. oleracea (GenBank accession number KY768855) is a typical circular DNA molecule with 329,613 bp in length. It includes a pair of directly repeated elements with a length of 7286 bp. The directly repeated elements separate the genome into a large single copy region of 229,375 bp and a small single copy region of 85,666 bp. The structure of the S. oleracea mitochondrial genome is consistent with that reported before, which was determined based on the physical mapping of mitochondrial DNA (Stern & Palmer Citation1986). The overall GC content of the genome is 43.4%, and the GC content in exons is slightly lower (41.8%). The genome contains 29 protein-coding genes and 4 pseudogenes, of which 6 have multiple exons. Two of the multi-exon genes, nad2 and nad5, were trans-spliced in the mitochondrial genome. The total number of exons in the genome is 44, with a total size of 27,708 bp (about 8.4% of the total mitogenome) and an average length of 630 bp. In addition, the genome contains 24 tRNA genes coding for 16 amino acids, and 3 rRNA genes including a 5S rRNA, an 18S rRNA, and a 26S rRNA.
A neighbour-joining phylogenetic tree of S. oleracea and 12 other representative plant species was constructed with single-copy orthologous genes in the mitogenomes using MEGA6 (Tamura et al. Citation2013). As shown in the phylogenetic tree, nearly all nodes were strongly supported with high bootstrap values (). The result suggested that as expected, S. oleracea was placed closely to Beta vulgaris (sugar beet), both of which belong to the Amaranthaceae family of the order of Caryophyllales, which constitutes the basal clade in core eudicots (Dohm et al. Citation2014).
Figure 1. Neighbour-joining phylogenetic tree of Spinacia oleracea and 12 other plant species based on the single-copy orthologous genes in the mitochondrial genomes. GenBank accession numbers: Arabidopsis thaliana (NC_001284); Beta vulgaris (NC_015099); Brassica napus (NC_008285); Capsicum annuum (NC_024624); Citrullus lanatus (GQ856147); Cucumis sativus (NC_016005); Lotus japonicus (NC_016743); Nicotiana tabacum (NC_006581); Vigna radiata (NC_015121); Vitis vinifera (NC_012119); Oryza sativa (DQ167399); Zea mays (NC_007982).
Disclosure statement
The authors report no conflicts of interest. The authors alone are responsible for the content and writing of the article.
Additional information
Funding
References
- Dohm JC, Minoche AE, Holtgräwe D, Capella-Gutiérrez S, Zakrzewski F, Tafer H, Rupp O, Sörensen TR, Stracke R, Reinhardt R, et al. 2014. The genome of the recently domesticated crop plant sugar beet (Beta vulgaris). Nature. 505:546–549.
- Luo R, Liu B, Xie Y, Li Z, Huang W, Yuan J, He G, Chen Y, Pan Q, Liu Y, et al. 2012. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler. GigaScience. 1:18.
- Ryder EJ. 1979. Spinach. In: Leafy salad vegetables. Netherlands: Springer; p. 195–227.
- Stern DB, Palmer JD. 1986. Tripartite mitochondrial genome of spinach: physical structure, mitochondrial gene mapping, and locations of transposed chloroplast DNA sequences. Nucleic Acids Res. 14:5651–5666.
- Tamura K, Stecher G, Peterson D, Filipski A, Kumar S. 2013. MEGA6: molecular evolutionary genetics analysis version 6.0. Mol Biol Evol. 30:2725–2272.
- van Treuren R, Coquin P, Lohwasser U. 2011. Genetic resources collections of leafy vegetables (lettuce, spinach, chicory, artichoke, asparagus, lamb’s lettuce, rhubarb and rocket salad): composition and gaps. Genet Resour Crop Ev. 59:981–997.
- Xu C, Jiao C, Sun H, Cai X, Wang X, Ge C, Zheng Y, Liu W, Sun X, Xu Y, et al. 2017. Draft genome of spinach and transcriptome diversity of 120 Spinacia accessions. Nat Commun. 8:15275.