De Novo Transcriptome Sequencing of the Orange-Fleshed Sweet Potato and Analysis of Differentially Expressed Genes Related to Carotenoid Biosynthesis

First systematic review of the last 30 years of research on sweetpotato: elucidating the frontiers and hotspots

Sweetpotato is an economically important crop, and it has various advantages over other crops in addressing global food security and climate change. Although substantial articles have been published on the research of various aspects of sweetpotato biology, there are no specific reports to systematically crystallize the research achievements. The current review takes the lead in conducting a keyword-centric spatiotemporal dimensional bibliometric analysis of articles on sweetpotato research using CiteSpace software to comprehensively clarify the development status, research hotspot, and development trend in the past 30 years (1993-2022). Quantitative analysis was carried out on the publishing countries, institutions, disciplines, and scholars to understand the basic status of sweetpotato research; then, visual analysis was conducted on high-frequency keywords, burst keywords, and keyword clustering; the evolution of major research hotspots and the development trend in different periods were summarized. Finally, the three main development stages-preliminary stage (1993-2005), rapid stage (2006-2013), and diversified mature stage (2014-2022)-were reviewed and analyzed in detail. Particularly, the development needs of sweetpotato production in improving breeding efficiency, enhancing stress tolerance, coordinating high yield with high quality and high resistance, and promoting demand were discussed, which will help to comprehensively understand the development dynamics of sweetpotato research from different aspects of biological exploration.

The transcription factor IbNAC29 positively regulates the carotenoid accumulation in sweet potato

Carotenoid is a tetraterpene pigment beneficial for human health. Although the carotenoid biosynthesis pathway has been extensively studied in plants, relatively little is known about their regulation in sweet potato. Previously, we conducted the transcriptome database of differentially expressed genes between the sweet potato (Ipomoea batatas) cultivar 'Weiduoli' and its high-carotenoid mutant 'HVB-3'. In this study, we selected one of these candidate genes, IbNAC29, for subsequent analyses. IbNAC29 belongs to the plant-specific NAC (NAM, ATAF1/2, and CUC2) transcription factor family. Relative IbNAC29 mRNA level in the HVB-3 storage roots was ~1.71-fold higher than Weiduoli. Additional experiments showed that the contents of α-carotene, lutein, β-carotene, zeaxanthin, and capsanthin are obviously increased in the storage roots of transgenic sweet potato plants overexpressing IbNAC29. Moreover, the levels of carotenoid biosynthesis genes in transgenic plants were also up-regulated. Nevertheless, yeast one-hybrid assays indicated that IbNAC29 could not directly bind to the promoters of these carotenoid biosynthesis genes. Furthermore, the level of IbSGR1 was down-regulated, whose homologous genes in tomato can negatively regulate carotene accumulation. Yeast three-hybrid analysis revealed that the IbNAC29-IbMYB1R1-IbAITR5 could form a regulatory module. Yeast one-hybrid, electrophoretic mobility shift assay, quantitative PCR analysis of chromatin immunoprecipitation and dual-luciferase reporter assay showed that IbAITR5 directly binds to and inhibits the promoter activity of IbSGR1, up-regulating carotenoid biosynthesis gene IbPSY. Taken together, IbNAC29 is a potential candidate gene for the genetic improvement of nutritive value in sweet potato.

Comparative transcriptome analysis implied a ZEP paralog was a key gene involved in carotenoid accumulation in yellow-fleshed sweetpotato

The mechanisms of carotenoid accumulation in yellow-fleshed sweetpotato cultivars are unclear. In this study, we compared the transcriptome profiles of a yellow-fleshed cultivar, Beniharuka (BH) and two of its spontaneous white-fleshed mutants (WH2 and WH3) to reveal the genes involved in yellow flesh. As a result of RNA sequencing, a total of 185 differentially expressed genes (DEGs) were commonly detected in WH2 and WH3 compared to BH. Of these genes, 85 DEGs and 100 DEGs were commonly upregulated and downregulated in WH2 and WH3 compared to BH, respectively. g1103.t1, a paralog of zeaxanthin epoxidase (ZEP), was only DEG common to WH2 and WH3 among 38 genes considered to be involved in carotenoid biosynthesis in storage roots. The expression level of g1103.t1 was also considerably lower in five white-fleshed cultivars than in five yellow-fleshed cultivars. Analysis of carotenoid composition in the storage roots showed that the epoxidised carotenoids were drastically reduced in both WH2 and WH3. Therefore, we propose that the ZEP paralog, g1103.t1, may be involved in carotenoid accumulation through the epoxidation of β-carotene and β-cryptoxanthin in sweetpotato.

Proteomic and metabolic profile analysis of low-temperature storage responses in Ipomoea batata Lam. tuberous roots

BACKGROUND

Sweetpotato (Ipomoea batatas L.) is one of the seven major food crops grown worldwide. Cold stress often can cause protein expression pattern and substance contents variations for tuberous roots of sweetpotato during low-temperature storage. Recently, we developed proteometabolic profiles of the fresh sweetpotatoes (cv. Xinxiang) in an attempt to discern the cold stress-responsive mechanism of tuberous root crops during post-harvest storage.

RESULTS

For roots stored under 4 °C condition, the CI index, REC and MDA content in roots were significantly higher than them at control temperature (13 °C). The activities of SOD, CAT, APX, O₂^.- producing rate, proline and especially soluble sugar contents were also significantly increased. Most of the differentially expressed proteins (DEPs) were implicated in pathways related to metabolic pathway, especially phenylpropanoids and followed by starch and sucrose metabolism. L-ascorbate peroxidase 3 and catalase were down-regulated during low temperature storage. α-amylase, sucrose synthase and fructokinase were significantly up-regulated in starch and sucrose metabolism, while β-glucosidase, glucose-1-phosphate adenylyl-transferase and starch synthase were opposite. Furthermore, metabolome profiling revealed that glucosinolate biosynthesis, tropane, piperidine and pyridine alkaloid biosynthesis as well as protein digestion and absorption played a leading role in metabolic pathways of roots. Leucine, tryptophan, tyrosine, isoleucine and valine were all significantly up-regulated in glucosinolate biosynthesis.

CONCLUSIONS

Our proteomic and metabolic profile analysis of sweetpotatoes stored at low temperature reveal that the antioxidant enzymes activities, proline and especially soluble sugar content were significantly increased. Most of the DEPs were implicated in phenylpropanoids and followed by starch and sucrose metabolism. The discrepancy between proteomic (L-ascorbate peroxidase 3 and catalase) and biochemical (CAT/APX activity) data may be explained by higher H₂O₂ levels and increased ascorbate redox states, which enhanced the CAT/APX activity indirectly. Glucosinolate biosynthesis played a leading role in metabolic pathways. Leucine, tryptophan, tyrosine, isoleucine and valine were all significantly up-regulated in glucosinolate biosynthesis.

Transcriptome sequencing and whole genome expression profiling of hexaploid sweetpotato under salt stress

BACKGROUND

Purple-fleshed sweetpotato (PFSP) is one of the most important crops in the word which helps to bridge the food gap and contribute to solve the malnutrition problem especially in developing countries. Salt stress is seriously limiting its production and distribution. Due to lacking of reference genome, transcriptome sequencing is offering a rapid approach for crop improvement with promising agronomic traits and stress adaptability.

RESULTS

Five cDNA libraries were prepared from the third true leaf of hexaploid sweetpotato at seedlings stage (Xuzi-8 cultivar) treated with 200 mM NaCl for 0, 1, 6, 12, 48 h. Using second and third generation technology, Illumina sequencing generated 170,344,392 clean high-quality long reads that were assembled into 15,998 unigenes with an average length 2178 base pair and 96.55% of these unigenes were functionally annotated in the NR protein database. A number of 537 unigenes failed to hit any homologs which may be considered as novel genes. The current results indicated that sweetpotato plants behavior during the first hour of salt stress was different than the other three time points. Furthermore, expression profiling analysis identified 4, 479, 281, 508 significantly expressed unigenes in salt stress treated samples at the different time points including 1, 6, 12, 48 h, respectively as compared to control. In addition, there were 4, 1202, 764 and 2195 transcription factors differentially regulated DEGs by salt stress at different time points including 1, 6, 12, 48 h of salt stress. Validation experiment was done using 6 randomly selected unigenes and the results was in agree with the DEG results. Protein kinases include many genes which were found to play a vital role in phosphorylation process and act as a signal transductor/ receptor proteins in membranes. These findings suggest that salt stress tolerance in hexaploid sweetpotato plants may be mainly affected by TFs, PKs, Protein Detox and hormones related genes which contribute to enhance salt tolerance.

CONCLUSION

These transcriptome sequencing data of hexaploid sweetpotato under salt stress conditions can provide a valuable resource for sweetpotato breeding research and focus on novel insights into hexaploid sweetpotato responses to salt stress. In addition, it offers new candidate genes or markers that can be used as a guide to the future studies attempting to breed salt tolerance sweetpotato cultivars.

Illumina® Sequencing Reveals Candidate Genes of Carotenoid Metabolism in Three Pummelo Cultivars (Citrus Maxima) with Different Pulp Color

Pummelo (Citrus maxima) is one of important fruit trees, which belongs to Citrus species. The fruits of different pummelo cultivars have different colors and differ in the contents of carotenoid. Our results clearly showed that ‘Huangjinmiyou’ (HJMY) has the highest content of β-carotene, followed by ‘Hongroumiyou’ (HRMY) and ‘Guanximiyou’ (GXMY). Lycopene is dominantly accumulated in HRMY. However, the molecular mechanism underlying the carotenoid accumulation in pummelo flesh is not fully understood. In this study, we used the RNA-Seq technique to investigate the candidate genes of carotenoid metabolism in the flesh of pummelo cv. GXMY and its mutants HRMY and HJMY in three development periods of fruit. After data assembly and bioinformatic analysis, a total of 357 genes involved in biosynthesis of secondary metabolites were isolated, of which 12 differentially expressed genes (DEGs) are involved in carotenoid biosynthesis. Among these 12 DEGs, phytoene synthase (PSY2), lycopene β-cyclase (LYCB2), lycopene Ɛ-cyclase (LYCE), carotenoid cleavage dioxygenases (CCD4), 9-cis-epoxycarotenoid dioxygenase (NCED2), aldehyde oxidase 3 (AAO3), and ABA 8′-hydroxylases (CYP707A1) are the most distinct DEGs in three pummelo cultivars. The co-expression analysis revealed that the expression patterns of several transcription factors such as bHLH, MYB, ERF, NAC and WRKY are highly correlated with DEGs, which are involved in carotenoid biosynthesis. In addition, the expression patterns of 22 DEGs were validated by real-time quantitative PCR (RT-qPCR) and the results are highly concordant with the RNA-Seq results. Our results provide a global vision of transcriptomic profile among three pummelo cultivars with different pulp colors. These results would be beneficial to further study the molecular mechanism of carotenoid accumulation in pummelo flesh and help the breeding of citrus with high carotenoid content.

Multi-species transcriptome analyses for the regulation of crocins biosynthesis in Crocus

BACKGROUND

Crocins are soluble apocarotenoids that mainly accumulate in the stigma tissue of Crocus sativus and provide the characteristic red color to saffron spice, in addition to being responsible for many of the medicinal properties of saffron. Crocin biosynthesis and accumulation in saffron is developmentally controlled, and the concentration of crocins increases as the stigma develops. Until now, little has been known about the molecular mechanisms governing crocin biosynthesis and accumulation. This study aimed to identify the first set of gene regulatory processes implicated in apocarotenoid biosynthesis and accumulation.

RESULTS

A large-scale crocin-mediated RNA-seq analysis was performed on saffron and two other Crocus species at two early developmental stages coincident with the initiation of crocin biosynthesis and accumulation. Pairwise comparison of unigene abundance among the samples identified potential regulatory transcription factors (TFs) involved in crocin biosynthesis and accumulation. We found a total of 131 (up- and downregulated) TFs representing a broad range of TF families in the analyzed transcriptomes; by comparison with the transcriptomes from the same developmental stages from other Crocus species, a total of 11 TF were selected as candidate regulators controlling crocin biosynthesis and accumulation.

CONCLUSIONS

Our study generated gene expression profiles of stigmas at two key developmental stages for apocarotenoid accumulation in three different Crocus species. Differential gene expression analyses allowed the identification of transcription factors that provide evidence of environmental and developmental control of the apocarotenoid biosynthetic pathway at the molecular level.

Transcriptome analysis reveals the effects of grafting on sweetpotato scions during the full blooming stages

BACKGROUND

Sweetpotato (Ipomoea batatas) is a hexaploid plant and generally most genotypes do not flower at all in sub-tropics. Heterografting was carried out between sweetpotato cultivar 'Xushu 18' and Japanese morning glory (Ipomoea nil). With sweetpotato as 'scion' and I. nil as 'rootstock', sweetpotato was induced flowering in the autumn. However, little is known about the molecular mechanisms underlying sweetpotato responses to grafting, especially during the full blooming stages.

OBJECTIVES

To investigate the poorly understood molecular responses underlying the grafting-induced phenotypic processes in sweetpotato at full anthesis.

METHODS

In this study, to explore the transcriptome diversity and complexity of sweetpotato, PacBio Iso-Seq and Illumina RNA-seq analysis were combined to obtain full-length transcripts and to profile the changes in gene expression of five tissues: scion flowers (SF), scion leaves (SL), scion stems (SS), own-rooted leaves (OL) and own-rooted stems (OS).

RESULTS

A total of 138,151 transcripts were generated with an average length of 2255 bp, and more than 72% (100,396) of the transcripts were full-length. During full blooming, to examine the difference in gene expression of sweetpotato under grafting and natural growth conditions, 7905, 7795 and 15,707 differentially expressed genes were detected in pairwise comparisons of OS versus SS, OL versus SL and SL versus SF, respectively. Moreover, differential transcription of genes associated with anthocyanin biosynthesis, light pathway and photosynthesis, ethylene signal transduction pathway was observed in scion responses to grafting.

CONCLUSION

Our study is useful in understanding the molecular basis of grafting-induced flowering in grafted sweetpotatoes, and will lay a foundation for further research on sweetpotato breeding in the future.

Genome-wide transcriptome profiling of transgenic hop (Humulus lupulus L.) constitutively overexpressing HlWRKY1 and HlWDR1 transcription factors

Background

The hop plant (Humulus lupulus L.) is a valuable source of several secondary metabolites, such as flavonoids, bitter acids, and essential oils. These compounds are widely implicated in the beer brewing industry and are having potential biomedical applications. Several independent breeding programs around the world have been initiated to develop new cultivars with enriched lupulin and secondary metabolite contents but met with limited success due to several constraints. In the present work, a pioneering attempt has been made to overexpress master regulator binary transcription factor complex formed by HlWRKY1 and HlWDR1 using a plant expression vector to enhance the level of prenylflavonoid and bitter acid content in the hop. Subsequently, we performed transcriptional profiling using high-throughput RNA-Seq technology in leaves of resultant transformants and wild-type hop to gain in-depth information about the genome-wide functional changes induced by HlWRKY1 and HlWDR1 overexpression.

Results

The transgenic WW-lines exhibited an elevated expression of structural and regulatory genes involved in prenylflavonoid and bitter acid biosynthesis pathways. In addition, the comparative transcriptome analysis revealed a total of 522 transcripts involved in 30 pathways, including lipids and amino acids biosynthesis, primary carbon metabolism, phytohormone signaling and stress responses were differentially expressed in WW-transformants. It was apparent from the whole transcriptome sequencing that modulation of primary carbon metabolism and other pathways by HlWRKY1 and HlWDR1 overexpression resulted in enhanced substrate flux towards secondary metabolites pathway. The detailed analyses suggested that none of the pathways or genes, which have a detrimental effect on physiology, growth and development processes, were induced on a genome-wide scale in WW-transgenic lines.

Conclusions

Taken together, our results suggest that HlWRKY1 and HlWDR1 simultaneous overexpression positively regulates the prenylflavonoid and bitter acid biosynthesis pathways in the hop and thus these transgenes are presented as prospective candidates for achieving enhanced secondary metabolite content in the hop.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-5125-8) contains supplementary material, which is available to authorized users.

A Sweetpotato Auxin Response Factor Gene (IbARF5) Is Involved in Carotenoid Biosynthesis and Salt and Drought Tolerance in Transgenic Arabidopsis

Auxin response factors (ARFs) compose a family of transcription factors and have been found to play major roles in the process of plant growth and development. However, their roles in plant carotenoid biosynthesis and responses to abiotic stresses are rarely known to date. In the present study, we found that the IbARF5 gene from sweetpotato (Ipomoea batatas (L.) Lam.) line HVB-3 increased the contents of carotenoids and enhanced the tolerance to salt and drought in transgenic Arabidopsis. The transgenic Arabidopsis plants exhibited the increased abscisic acid (ABA) and proline contents and superoxide dismutase (SOD) activity and the decreased H₂O₂ content. Furthermore, it was found that IbARF5 positively regulated the genes associated with carotenoid and ABA biosynthesis and abiotic stress responses. These results suggest that IbARF5 is involved in carotenoid biosynthesis and salt and drought tolerance in transgenic Arabidopsis. This study provides a novel ARF gene for improving carotenoid contents and salt and drought tolerance of sweetpotato and other plants.

Involvement of an ABI-like protein and a Ca2+-ATPase in drought tolerance as revealed by transcript profiling of a sweetpotato somatic hybrid and its parents Ipomoea batatas (L.) Lam. and I. triloba L

Previously, we obtained the sweetpotato somatic hybrid KT1 from a cross between sweetpotato (Ipomoea batatas (L.) Lam.) cv. Kokei No. 14 and its drought-tolerant wild relative I. triloba L. KT1 not only inherited the thick storage root characteristic of Kokei No. 14 but also the drought-tolerance trait of I. triloba L. The aim of this study was to explore the molecular mechanism of the drought tolerance of KT1. Four-week-old in vitro-grown plants of KT1, Kokei No. 14, and I. triloba L. were subjected to a simulated drought stress treatment (30% PEG6000) for 0, 6, 12 and 24 h. Total RNA was extracted from samples at each time point, and then used for transcriptome sequencing. The gene transcript profiles of KT1 and its parents were compared to identify differentially expressed genes, and drought-related modules were screened by a weighted gene co-expression network analysis. The functions of ABI-like protein and Ca2+-ATPase, two proteins screened from the cyan and light yellow modules, were analyzed in terms of their potential roles in drought tolerance in KT1 and its parents. These analyses of the drought responses of KT1 and its somatic donors at the transcriptional level provide new annotations for the molecular mechanism of drought tolerance in the somatic hybrid KT1 and its parents.

Transcriptome profiling and digital gene expression analysis of sweet potato for the identification of putative genes involved in the defense response against Fusarium oxysporum f. sp. batatas

Sweet potato production is constrained by Fusarium wilt, which is caused by Fusarium oxysporum f. sp. batatas (Fob). The identification of genes related to disease resistance and the underlying mechanisms will contribute to improving disease resistance via sweet potato breeding programs. In the present study, we performed de novo transcriptome assembly and digital gene expression (DGE) profiling of sweet potato challenged with Fob using Illumina HiSeq technology. In total, 89,944,188 clean reads were generated from 12 samples and assembled into 101,988 unigenes with an average length of 666 bp; of these unigenes, 62,605 (61.38%) were functionally annotated in the NCBI non-redundant protein database by BLASTX with a cutoff E-value of 10-5. Clusters of Orthologous Groups (COG), Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) annotations were examined to explore the unigenes' functions. We constructed four DGE libraries for the sweet potato cultivars JinShan57 (JS57, highly resistant) and XinZhongHua (XZH, highly susceptible), which were challenged with pathogenic Fob. Genes that were differentially expressed in the four libraries were identified by comparing the transcriptomes. Various genes that were differentially expressed during defense, including chitin elicitor receptor kinase 1 (CERK), mitogen-activated protein kinase (MAPK), WRKY, NAC, MYB, and ethylene-responsive transcription factor (ERF), as well as resistance genes, pathogenesis-related genes, and genes involved in salicylic acid (SA) and jasmonic acid (JA) signaling pathways, were identified. These data represent a sequence resource for genetic and genomic studies of sweet potato that will enhance the understanding of the mechanism of disease resistance.

Transcript profile analysis reveals important roles of jasmonic acid signalling pathway in the response of sweet potato to salt stress

Sweet potato is an important food and bio-energy crop, and investigating the mechanisms underlying salt tolerance will provide information for salt-tolerant breeding of this crop. Here, the root transcriptomes of the salt-sensitive variety Lizixiang and the salt-tolerant line ND98 were compared to identify the genes and pathways involved in salt stress responses. In total, 8,744 and 10,413 differentially expressed genes (DEGs) in Lizixiang and ND98, respectively, were involved in salt responses. A lower DNA methylation level was detected in ND98 than in Lizixiang. In both genotypes, the DEGs, which function in phytohormone synthesis and signalling and ion homeostasis, may underlie the different degrees of salt tolerance. Significant up-regulations of the genes involved in the jasmonic acid (JA) biosynthesis and signalling pathways and ion transport, more accumulation of JA, a higher degree of stomatal closure and a lower level of Na⁺ were found in ND98 compared to Lizixiang. This is the first report on transcriptome responses to salt tolerance in sweet potato. These results reveal that the JA signalling pathway plays important roles in the response of sweet potato to salt stress. This study provides insights into the mechanisms and genes involved in the salt tolerance of sweet potato.

Challenges to genome sequence dissection in sweetpotato

The development of next generation sequencing (NGS) technologies has enabled the determination of whole genome sequences in many non-model plant species. However, genome sequencing in sweetpotato (Ipomoea batatas (L.) Lam) is still difficult because of the hexaploid genome structure. Previous studies suggested that a diploid wild relative, I. trifida (H.B.K.) Don., is the most possible ancestor of sweetpotato. Therefore, the genetic and genomic features of I. trifida have been studied as a potential reference for sweetpotato. Meanwhile, several research groups have begun the challenging task of directly sequencing the sweetpotato genome. In this manuscript, we review the recent results and activities of large-scale genome and transcriptome analysis related to genome sequence dissection in sweetpotato under the sections as follows: I. trifida genome and transcript sequencing, genome sequences of I. nil (Japanese morning glory), transcript sequences in sweetpotato, chloroplast sequences, transposable elements and transfer DNA. The recent international activities of de novo whole genome sequencing in sweetpotato are also described. The large-scale publically available genome and transcript sequence resources and the international genome sequencing streams are expected to promote the genome sequence dissection in sweetpotato.

Improvement for agronomically important traits by gene engineering in sweetpotato

Sweetpotato is the seventh most important food crop in the world. It is mainly used for human food, animal feed, and for manufacturing starch and alcohol. This crop, a highly heterozygous, generally self-incompatible, outcrossing polyploidy, poses numerous challenges for the conventional breeding. Its productivity and quality are often limited by abiotic and biotic stresses. Gene engineering has been shown to have the great potential for improving the resistance to these stresses as well as the nutritional quality of sweetpotato. To date, an Agrobacterium tumefaciens-mediated transformation system has been developed for a wide range of sweetpotato genotypes. Several genes associated with salinity and drought tolerance, diseases and pests resistance, and starch, carotenoids and anthocyanins biosynthesis have been isolated and characterized from sweetpotato. Gene engineering has been used to improve abiotic and biotic stresses resistance and quality of this crop. This review summarizes major research advances made so far in improving agronomically important traits by gene engineering in sweetpotato and suggests future prospects for research in this field.

A genome-wide BAC-end sequence survey provides first insights into sweetpotato (Ipomoea batatas (L.) Lam.) genome composition

Background

Sweetpotato, Ipomoea batatas (L.) Lam., is an important food crop widely grown in the world. However, little is known about the genome of this species because it is a highly heterozygous hexaploid. Gaining a more in-depth knowledge of sweetpotato genome is therefore necessary and imperative. In this study, the first bacterial artificial chromosome (BAC) library of sweetpotato was constructed. Clones from the BAC library were end-sequenced and analyzed to provide genome-wide information about this species.

Results

The BAC library contained 240,384 clones with an average insert size of 101 kb and had a 7.93–10.82 × coverage of the genome, and the probability of isolating any single-copy DNA sequence from the library was more than 99%. Both ends of 8310 BAC clones randomly selected from the library were sequenced to generate 11,542 high-quality BAC-end sequences (BESs), with an accumulative length of 7,595,261 bp and an average length of 658 bp. Analysis of the BESs revealed that 12.17% of the sweetpotato genome were known repetitive DNA, including 7.37% long terminal repeat (LTR) retrotransposons, 1.15% Non-LTR retrotransposons and 1.42% Class II DNA transposons etc., 18.31% of the genome were identified as sweetpotato-unique repetitive DNA and 10.00% of the genome were predicted to be coding regions. In total, 3,846 simple sequences repeats (SSRs) were identified, with a density of one SSR per 1.93 kb, from which 288 SSRs primers were designed and tested for length polymorphism using 20 sweetpotato accessions, 173 (60.07%) of them produced polymorphic bands. Sweetpotato BESs had significant hits to the genome sequences of I. trifida and more matches to the whole-genome sequences of Solanum lycopersicum than those of Vitis vinifera, Theobroma cacao and Arabidopsis thaliana.

Conclusions

The first BAC library for sweetpotato has been successfully constructed. The high quality BESs provide first insights into sweetpotato genome composition, and have significant hits to the genome sequences of I. trifida and more matches to the whole-genome sequences of Solanum lycopersicum. These resources as a robust platform will be used in high-resolution mapping, gene cloning, assembly of genome sequences, comparative genomics and evolution for sweetpotato.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3302-1) contains supplementary material, which is available to authorized users.